Effect of Nanoscale Confinement on Freezing of Modified Water at Room Temperature and Ambient Pressure

Understanding the phase behavior of confined water is central to fields as diverse as heterogeneous catalysis, corrosion, nanofluidics, and to emerging energy technologies. Altering the state points (temperature, pressure, etc.) or introduction of a foreign surface can result in the phase transformation of water. At room temperature, ice nucleation is a very rare event and extremely high pressures in the GPa–TPa range are required to freeze water. Here, we perform computer experiments to artificially alter the balance between electrostatic and dispersion interactions between water molecules, and demonstrate nucleation and growth of ice at room temperature in a nanoconfined environment. Local perturbations in dispersive and electrostatic interactions near the surface are shown to provide the seed for nucleation (nucleation sites), which lead to room temperature liquid–solid phase transition of confined water. Crystallization of water occurs over several tens of nanometers and is shown to be independent of the nature of the substrate (hydrophilic oxide vs. hydrophobic graphene and crystalline oxide vs. amorphous diamond‐like carbon). Our results lead us to hypothesize that the freezing transition of confined water can be controlled by tuning the relative dispersive and electrostatic interaction.     

Collective motion artifacts arising in long‐duration molecular dynamics simulations

We tested a variety of molecular dynamics simulation strategies in long‐duration (up to several nanoseconds) constant‐temperature simulations of liquid water under periodic boundary conditions. Such long durations are necessary to achieve adequate conformational sampling in simulations of membrane assemblies and other large biomolecular systems. Under a variety of circumstances, serious artifacts arise in the form of spurious collective behavior that becomes obvious only after the simulation has gone at least several hundred picoseconds. The potential energy of the system drops and the system changes from a liquid to an icy or glassy state. The underlying cause is accumulated center‐of‐mass motion of the system, coupled with velocity rescaling associated with constant‐temperature control. The velocity rescaling in the constant‐temperature algorithm reduces the thermal velocity as the net center‐of‐mass velocity grows, effectively causing the kinetic energy of the system to drain from thermal motions into coordinated motions. We found that the incidence and magnitude of the underlying artifactual motion leading to the spurious transition is mediated by: choice of method for computing electrostatic interactions; choice of ensemble; size of the simulation cell; SHAKE tolerance; frequency of nonbonded pairlist updating; and closeness of coupling to the temperature bath. The appearance of the spurious transition can be avoided by periodically subtracting net center‐of‐mass motion during the dynamics, or by improving the accuracy of the simulation by means of tightening SHAKE tolerance and updating nonbonded pairlists every timestep. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 121–131, 2000     

计算离子液体溶液汽液相平衡的分子热力学模型

用平均球近似理论、微扰理论和UNIFAC基团贡献方法分别考虑离子之间的长程静电作用、离子与溶剂之间的中程静电作用以及所有粒子之间的短程作用,本文提出了一种新的分子热力学模型,可用于离子液体溶液中溶剂活度系数的计算.通过对含烷基咪唑磷酸酯类离子液体与水、甲醇或乙醇组成的9个二元体系的饱和蒸汽压数据进行关联,获得了相关的模型参数,即溶剂的分子直径和基团之间的交互作用能参数.溶剂活度系数及饱和蒸汽压的计算结果与实验值的平均偏差为1.40%,符合良好,因此本模型可望用于含离子液体体系汽液相平衡的预测.     

Characterization of the TIP4P-Ew water model: vapor pressure and boiling point

The liquid-vapor-phase equilibrium properties of the previously developed TIP4P-Ew water model have been studied using thermodynamic integration free-energy simulation techniques in the temperature range of 274-400 K. We stress that free-energy results from simulations need to be corrected in order to be compared to the experiment. This is due to the fact that the thermodynamic end states accessible through simulations correspond to fictitious substances (classical rigid liquids and classical rigid ideal gases) while experiments operate on real substances (liquids and real gases, with quantum effects). After applying analytical corrections the vapor pressure curve obtained from simulated free-energy changes is in excellent agreement with the experimental vapor pressure curve. The boiling point of TIP4P-Ew water under ambient pressure is found to be at 370.3+/-1.9 K, about 7 K higher than the boiling point of TIP4P water (363.7+/-5.1 K; from simulations that employ finite range treatment of electrostatic and Lennard-Jones interactions). This is in contrast to the approximately +15 K by which the temperature of the density maximum and the melting temperature of TIP4P-Ew are shifted relative to TIP4P, indicating that the temperature range over which the liquid phase of TIP4P-Ew is stable is narrower than that of TIP4P and resembles more that of real water. The quality of the vapor pressure results highlights the success of TIP4P-Ew in describing the energetic and entropic aspects of intermolecular interactions in liquid water.     

Thermal hysteresis in dithiadiazolyl and dithiazolyl radicals induced by supercooling of paramagnetic liquids close to room temperature: a study of F3CCNSSN and an interpretation of the behaviour of F3CCSNSCCF3

The trifluoromethyl-substituted dithiadiazolyl and dithiazolyl radicals, F3CCNSSN (1) and F3CCSNSCCF3 (2) associate through pi*-pi* covalent and electrostatic S delta+...N delta- interactions in the solid state, but melt with a dramatic volume increase to generate paramagnetic liquids; these radicals exhibit thermal hysteresis, which arises through a meta-stable super-cooled liquid state, close to room temperature.     

Theory of surface micelles of semifluorinated alkanes

Surface structures of semifluorinated alkanes F(CF(2))(n)(CH(2))(m)H (referred to as FnHm) spread on the air/water interface are investigated theoretically. The study is focused on the disklike surface micelles that were recently identified by AFM and scattering techniques at sufficiently high surface concentrations. We show that (1) the micelles emerge as a result of liquid/liquid (rather than liquid/gas) phase separation in the Langmuir layer; (2) the micelles are islands of the higher-density phase with roughly vertical orientation of FnHm molecules (F-parts extend toward air, H-parts toward water) and the matrix is the lower density-phase where the FnHm diblocks are nearly parallel to the water surface; (3) the micelles and the hexagonal structure they form are stabilized by the electrostatic interactions which are mainly due to the vertical dipole moments of the CF(2)- CH(2) bonds in the vertical phase; and (4) the electrostatic repulsive interactions can serve to suppress the micelle size polydispersity.     

Water under temperature gradients: polarization effects and microscopic mechanisms of heat transfer

We report non-equilibrium molecular dynamics simulations (NEMD) of water under temperature gradients using a modified version of the central force model (MCFM). This model is very accurate in predicting the equation of state of water for a wide range of pressures and temperatures. We investigate the polarization response of water to thermal gradients, an effect that has been recently predicted using Non-Equilibrium Thermodynamics (NET) theory and computer simulations, as a function of the thermal gradient strength. We find that the polarization of the liquid varies linearly with the gradient strength, which indicates that the ratio of phenomenological coefficients regulating the coupling between the polarization response and the heat flux is independent of the gradient strength investigated. This notion supports the NET theoretical predictions. The coupling effect leading to the liquid polarization is fairly strong, leading to polarization fields of ~10(3-6) V m(-1) for gradients of ~10(5-8) K m(-1), hence confirming earlier estimates. Finally we employ our NEMD approach to investigate the microscopic mechanism of heat transfer in water. The image emerging from the computation and analysis of the internal energy fluxes is that the transfer of energy is dominated by intermolecular interactions. For the MCFM model, we find that the contribution from hydrogen and oxygen is different, with the hydrogen contribution being larger than that of oxygen.     

A reoptimization of the five-site water potential (TIP5P) for use with Ewald sums

The five-site transferable interaction potential (TIP5P) for water is most accurate at reproducing experimental data when used with a simple spherical cutoff for the long-ranged electrostatic interactions. When used with other methods for treating long-ranged interactions, the model is considerably less accurate. With small modifications, a new TIP5P-like potential can be made which is very accurate for liquid water when used with Ewald sums, a more physical and increasingly more commonly used method for treating long-ranged electrostatic interactions. The new model demonstrates a density maximum near 4 degrees C, like the TIP5P model, and otherwise is similar to the TIP5P model for thermodynamic, dielectric, and dynamical properties of liquid water over a range of temperatures and densities. An analysis of this and other commonly used water models reveals how the quadrupole moment of a model can influence the dielectric response of liquid water.     

A simple, efficient polarizable coarse-grained water model for molecular dynamics simulations

The development of coarse-grained (CG) models that correctly represent the important features of compounds is essential to overcome the limitations in time scale and system size currently encountered in atomistic molecular dynamics simulations. Most approaches reported in the literature model one or several molecules into a single uncharged CG bead. For water, this implicit treatment of the electrostatic interactions, however, fails to mimic important properties, e.g., the dielectric screening. Therefore, a coarse-grained model for water is proposed which treats the electrostatic interactions between clusters of water molecules explicitly. Five water molecules are embedded in a spherical CG bead consisting of two oppositely charged particles which represent a dipole. The bond connecting the two particles in a bead is unconstrained, which makes the model polarizable. Experimental and all-atom simulated data of liquid water at room temperature are used for parametrization of the model. The experimental density and the relative static dielectric permittivity were chosen as primary target properties. The model properties are compared with those obtained from experiment, from clusters of simple-point-charge water molecules of appropriate size in the liquid phase, and for other CG water models if available. The comparison shows that not all atomistic properties can be reproduced by a CG model, so properties of key importance have to be selected when coarse graining is applied. Yet, the CG model reproduces the key characteristics of liquid water while being computationally 1-2 orders of magnitude more efficient than standard fine-grained atomistic water models.     

Thermal study of accumulation of conformational disorders in the self-assembled monolayers of C8 and C18 alkanethiols on the Au111 surface

The thermal stability of short alkanethiol CH(3)(CH(2))(7)SH (C(8)) and long C(18) self-assembled monolayers (SAMs) is investigated using grazing angle reflection-absorption infrared spectroscopy, cyclic voltammetry, and molecular dynamics simulation. We track the disordering of SAM by untilting and gauche defect accumulation with increasing temperature in the 300-440 K range, a range of interest to tribology. Molecular dynamics simulation with both fully covered and partially covered C(6), C(8), and C(18) monolayers brings out the morphological changes in the SAM, which may be associated with the observed thermal stability characteristics. The molecular dynamics simulations reveal that short-chain C(6) and C(8) alkanethiols are more defective at lower temperature than the long-chain C(18) alkanethiol. With increasing temperature disorder in the SAM, as reflected in both untilting and gauche defect accumulation, tends to saturate at temperatures below 360 K for short-chain SAMs such that any further increase in temperature, until desorption, does not lead to any significant change in conformational order. In contrast the disorder in the long-chain C(18) SAM increases monotonically with temperature beyond 360 K. Thus, in a practical range of temperature, the ability of a SAM to retain order with increasing thermal perturbations is governed by the state of disorder prior to heat treatment. This deduction derived from molecular dynamics simulation helps to rationalize the significant difference we have observed experimentally between the thermal response of short- and long-chain thiol molecules.     

Balancing local order and long-ranged interactions in the molecular theory of liquid water

A molecular theory of liquid water is identified and studied on the basis of computer simulation of the TIP3P model of liquid water. This theory would be exact for models of liquid water in which the intermolecular interactions vanish outside a finite spatial range, and therefore provides a precise analysis tool for investigating the effects of longer-ranged intermolecular interactions. We show how local order can be introduced through quasichemical theory. Long-ranged interactions are characterized generally by a conditional distribution of binding energies, and this formulation is interpreted as a regularization of the primitive statistical thermodynamic problem. These binding-energy distributions for liquid water are observed to be unimodal. The Gaussian approximation proposed is remarkably successful in predicting the Gibbs free energy and the molar entropy of liquid water, as judged by comparison with numerically exact results. The remaining discrepancies are subtle quantitative problems that do have significant consequences for the thermodynamic properties that distinguish water from many other liquids. The basic subtlety of liquid water is found then in the competition of several effects which must be quantitatively balanced for realistic results.     

Multiple time scale behaviors and network dynamics in liquid methanol

Canonical ensemble molecular dynamics simulations of liquid methanol, modeled using a rigid-body, pair-additive potential, are used to compute static distributions and temporal correlations of tagged molecule potential energies as a means of characterizing the liquid state dynamics. The static distribution of tagged molecule potential energies shows a clear multimodal structure with three distinct peaks, similar to those observed previously in water and liquid silica. The multimodality is shown to originate from electrostatic effects, but not from local, hydrogen bond interactions. An interesting outcome of this study is the remarkable similarity in the tagged potential energy power spectra of methanol, water, and silica, despite the differences in the underlying interactions and the dimensionality of the network. All three liquids show a distinct multiple time scale (MTS) regime with a 1/ f (alpha) dependence with a clear positive correlation between the scaling exponent alpha and the diffusivity. The low-frequency limit of the MTS regime is determined by the frequency of crossover to white noise behavior which occurs at approximately 0.1 cm (-1) in the case of methanol under standard temperature and pressure conditions. The power spectral regime above 200 cm (-1) in all three systems is dominated by resonances due to localized vibrations, such as librations. The correlation between alpha and the diffusivity in all three liquids appears to be related to the strength of the coupling between the localized motions and the larger length/time scale network reorganizations. Thus, the time scales associated with network reorganization dynamics appear to be qualitatively similar in these systems, despite the fact that water and silica both display diffusional anomalies but methanol does not.     

Validity of the "sharp-kink approximation" for water and other fluids

The contact angle of a liquid droplet on a solid surface is a direct measure of fundamental atomic-scale forces acting between liquid molecules and the solid surface. In this work, the validity is assessed of a simple equation, which approximately relates the contact angle of a liquid on a surface to its density, its surface tension, and the effective molecule-surface potential. This equation is derived in the sharp-kink approximation, where the density profile of the liquid is assumed to drop precipitously within one molecular diameter of the substrate. It is found that this equation satisfactorily reproduces the temperature-dependence of the contact angle for helium on alkali metal surfaces. The equation also seems be applicable to liquids such as water on solid surfaces such as gold and graphite, on the basis of a comparison of predicted and measured contact angles near room-temperature. Nevertheless, we conclude that, to fully test the equation's applicability to fluids such as water, it remains necessary to measure the contact angle's temperature-dependence. We hypothesize that the effects of electrostatic forces can increase with temperature, potentially driving the wetting temperature much higher and closer to the critical point, or lower, closer to room temperature, than predicted using current theories.     

Imidazoline type stationary phase for hydrophilic interaction chromatography and reversed-phase liquid chromatography

An imidazoline was prepared by solvent-free microwave-assisted organic synthesis and immobilized on porous silica particles by polymerization. The resulting material was composed of both hydrophobic alkyl ester chains and hydrophilic imidazoline rings, which gave it both hydrophilic interaction and reversed-phase characteristics. The titration curve suggests that the new material has buffering capacity and acquires increasing positive charge over the pH range 9-4, and is "zwitterionic" in the upper part of this pH range. Through investigating the effect of column temperature, the water content, pH and ion strength of mobile phase on the retention time of polar compounds in highly organic eluents, it was found that the new material could be used as a hydrophilic interaction liquid chromatography (HILIC) stationary phase which involved a complex retention process consisting of partitioning, surface adsorption and electrostatic interactions. In addition, the retention behavior of aromatic compounds in different mobile phase conditions was also studied, which showed the new material mainly exhibited a partitioning mechanism in the reversed-phase liquid chromatography (RPLC) mode. The separation of six water-soluble vitamins and five aromatic compounds were achieved by using the new material in the HILIC and RPLC modes, respectively.     

Use of noncovalent interactions to form novel liquid crystalline polymeric materials

A novel concept in polymeric liquid crystalline (LC) materials, via electrostatic grafting or complexation of functionalized mesogens and appropriate polymers, is described. It is applied to three specific systems characterized by side-chain LC architecture, two of which involve hydrogen-bonding interactions, the third ionic interactions. In these systems, complexation tends to create greater order in the mesophases present. In the hydrogen-bonded systems, the transition temperature from the isotropic state tends to decrease compared to that of the functionalized mesogen; in the ionically-bonded system, the reverse is true. At low concentrations, the mesogens act as plasticizers.     

Electrostatics in protein binding and function

Protein electrostatic properties stem from the proportion and distribution of polar and charged residues. Polar and charged residues regulate the electrostatic properties by forming short-range interactions, like salt-bridges and hydrogen-bonds, and by defining the over-all electrostatic environment in the protein. Electrostatics play a major role in defining the mechanisms of protein-protein complex formation, molecular recognitions, thermal stabilities, conformational adaptabilities and protein movements. For example:- Functional hinges, or flexible regions of the protein, lack short-range electrostatic interactions; Thermophilic proteins have higher electrostatic interactions than their mesophilic counter parts; Increase in binding specificity and affinity involve optimization of electrostatics; High affinity antibodies have higher, and stronger, electrostatic interactions with their antigens; Rigid parts of proteins have higher and stronger electrostatic interactions. In this review we address the significance of electrostatics in protein folding, binding and function. We discuss that the electrostatic properties are evolutionally selected by a protein to perform an specific function. We also provide bona fide examples to illustrate this. Additionally, using continuum electrostatic and molecular dynamics approaches we show that the "hot-spot" inter-molecular interactions in a very specific antibody-antigen binding are mainly established through charged residues. These "hot-spot" molecular interactions stay intact even during high temperature molecular dynamics simulations, while the other inter-molecular interactions, of lesser functional significance, disappear. This further corroborates the significance of charge-charge interactions in defining binding mechanisms. High affinity binding frequently involves "electrostatic steering". The forces emerge from over-all electrostatic complementarities and by the formation of charged and polar interactions. We demonstrate that although the high affinity binding of barnase-barstar and anti-hen egg white lysozyme (HEL) antibody-HEL complexes involve different molecular mechanisms, it is electrostatically regulated in both the cases. These observations, and several other studies, suggest that a fine tuning of local and global electrostatic properties are essential for protein binding and function.     

Infrared Absorption Intensity Analysis as a New Tool for Investigation of Salt Effect on Proteins

The native protein structures in buffer solution are maintained by the electrostatic force as well as the hydrophobic force, salt ions play an important role in maintaining the protein native structures, and their effect on the protein stability has attracted tremendous interests. Infrared spectroscopy has been generally used in molecular structure analysis due to its fingerprint resolution for different species including macromolecules as proteins. However spectral intensities have received much less attention than the vibrational frequencies. Here we report that the spectral intensities of protein amide I band, the finger prints for the protein secondary structures, are very sensitive to the local electric field known as Onsager reaction field caused by salt ions. IR absorbance thermal titrations have been conducted for a series of samples including simple water soluble amino acids, water soluble monomeric protein cytochrome c and dimeric protein DsbC and its single-site mutant G49R. We found that at lower temperature range （10-20℃）, there exists a thermal activated salting-in process, where the IR intensity increases with a rise in the temperature, corresponding to the ions binding of the hydrophobic surface of protein. This process is absent for the amino acids. When further raising the temperature, the IR intensity decreases, this is interpreted as the thermal activated breaking of the ion-protein surface binding. Applying Van＇t Hoff plot to the thermal titration curves, the thermodynamic parameters such as AH and AS for salting-in and ion unbinding processes can be derived for various protein secondary structural components, revealing quantitatively the extent of hydrophobic interaction as well as the strength of the ion-protein binding.     

Tunable phase‐separation behavior of thermoresponsive polyampholytes through molecular design

Here, we report that carboxylated poly‐l ‐lysine, a polyampholyte, shows lower critical solution temperature (LCST)‐type temperature‐responsive liquid–liquid phase separation and coacervate formation in aqueous solutions. The phase‐separation temperature of polyampholytes is strongly affected by the polymer concentration, balance between the carboxyl and amino groups, hydrophobicity of the side chain, and NaCl concentration in the solution. We concluded that the phase separation was caused by both electrostatic interactions between the carboxyl and amino groups and intermolecular hydrophobic interactions. The addition of NaCl weakened the electrostatic interactions, causing the two phases to remix. The introduction of the hydrophobic moiety decreased the phase‐separation temperature by making the molecular interactions stronger. Finally, temperature‐responsive hydrogels were prepared from the polyampholytes to explore their applicability as biomaterials and in drug delivery systems. The fine‐tuning of the phase‐separation temperature of poly‐l ‐lysine‐based polyampholytes through molecular design should open new avenues for their use in precisely controlled biomedical applications. © 2016 Wiley Periodicals, Inc. J. Polym. Sci., Part A: Polym. Chem. 2017 , 55, 876–884     

A plastic phase of water from computer simulation

We report a member of ices called plastic or rotator phase, in which individual water molecules make facile rotations as in liquid state but are held tightly in an ordered structure. Molecular dynamics simulations of three classical models of water show that a plastic ice phase appears at a temperature when ice VII is heated or liquid water is cooled at high pressures above several gigapascals. A large amount of latent heat is absorbed when ice VII is transformed to the rotator phase at 590 K and 10 GPa, which is a typical characteristic of the plastic transitions for nearly spherical molecules. In addition to the spontaneous formation of plastic phase in the simulations, its existence is supported by robustness of plastic phase for hypothetical water with varying degrees of Coulombic interactions.