Structure and ultrafast dynamics of liquid water: a quantum mechanics/molecular mechanics molecular dynamics simulations study

A quantum mechanics/molecular mechanics molecular dynamics simulation was performed for liquid water to investigate structural and dynamical properties of this peculiar liquid. The most important region containing a central reference molecule and all nearest surrounding molecules (first coordination shell) was treated by Hartree-Fock (HF), post-Hartree-Fock [second-order Moller-Plesset perturbation theory (MP2)], and hybrid density functional B3LYP [Becke's three parameter functional (B3) with the correlation functional of Lee, Yang, and Parr (LYP)] methods. In addition, another HF-level simulation (2HF) included the full second coordination shell. Site to site interactions between oxygen-oxygen, oxygen-hydrogen, and hydrogen-hydrogen atoms of all ab initio methods were compared to experimental data. The absence of a second peak and the appearance of a shoulder instead in the gO-O graph obtained from the 2HF simulation is notable, as this feature has been observed so far only for pressurized or heated water. Dynamical data show that the 2HF procedure compensates some of the deficiency of the HF one-shell simulation, reducing the difference between correlated (MP2) and HF results. B3LYP apparently leads to too rigid structures and thus to an artificial slow down of the dynamics.     

On the application of accelerated molecular dynamics to liquid water simulations

Our group recently proposed a robust bias potential function that can be used in an efficient all-atom accelerated molecular dynamics (MD) approach to simulate the transition of high energy barriers without any advance knowledge of the potential-energy landscape. The main idea is to modify the potential-energy surface by adding a bias, or boost, potential in regions close to the local minima, such that all transitions rates are increased. By applying the accelerated MD simulation method to liquid water, we observed that this new simulation technique accelerates the molecular motion without losing its microscopic structure and equilibrium properties. Our results showed that the application of a small boost energy on the potential-energy surface significantly reduces the statistical inefficiency of the simulation while keeping all the other calculated properties unchanged. On the other hand, although aggressive acceleration of the dynamics simulation increases the self-diffusion coefficient of water molecules greatly and dramatically reduces the correlation time of the simulation, configurations representative of the true structure of liquid water are poorly sampled. Our results also showed the strength and robustness of this simulation technique, which confirm this approach as a very useful and promising tool to extend the time scale of the all-atom simulations of biological system with explicit solvent models. However, we should keep in mind that there is a compromise between the strength of the boost applied in the simulation and the reproduction of the ensemble average properties.     

Lamellar liquid crystal polymerization of sodium oleate/oleic acid/octadiene/water system

In the lamellar liquid crystallization (LLC) phase of NaOL/ OLA/H2O system, the small angle X-ray diffraction measurements show that the oleic acid is solubilized in the oil layer at first and then into the amphiphile layer. The octadiene added is also located partly in the oil layer and partly in the amphiphile layer in the LLC. With the addition of octadiene as cross-Unking agent, the LLC phase of NaOL/OIA/H2O system was polymerized under the initiation of AIBN with the protection of pure nitrogen at 60℃. Most of the double bond absorption of the monomers in IR spectra disappeared after polymerization. The polymerization takes place not only in the middle of the amphiphile layer between the double bonds of NaOL or OLA and those of octadiene, but also in the oil layer of LLC between the double bonds of OLA and those of octadiene. Interlayer spacing measurements on the copolymer proved d values decreased by about 1-2 nm compared with those of the corresponding system before the polymerization, indicati     

Phase behavior of amphiphiles at liquid crystals/water interface: A coarse-grained molecular dynamics study

We present a coarse-grained molecular dynamics simulation study of phase behavior of amphiphilic monolayers at the liquid crystal （LC）/water interface. The results revealed that LCs at interface can influence the lateral ordering of amphiphiles. Particularly, the amphiphile tails along with perpendicularly penetrated LCs between tails undergo a two-dimension phase transition from liquid-expanded into a liquid-condensed phase as their area density at interface reaches 0.93. While, the liquid-condensed phase of the monolayer never appears at oil/water interface with isotropic shape oil particles. These findings reveal the penetration of anisotropic LC can promote ordered lateral organization of amphiphiles. Moreover, we find the phase transition point is shifted to lower surface coverage of amphiphiles when the LCs have larger affinity to the amphiphile tails.     

Liquid–liquid phase equilibria of the ternary system of water/1,4-dioxane/dihydromyrcene

Liquid–liquid phase equilibria (LLE) for the system of water/1,4-dioxane/dihydromyrene (DHM) were investigated experimentally at 343.15, 348.15, 353.15, 358.15 K and atmospheric pressure. The reliability of the experimental tie-lines has been confirmed by using Othmer–Tobias correlation. The NRTL and improved UNIQUAC models were used to correlate the phase equilibria in the system using the interaction parameters determined from experimental data. The root mean square deviation (RMSD) between the observed and calculated mole percents was 0.58%. It is found that NRTL and improved UNIQUAC used for LLE could provide a good correlation.     

The use of two-phase molecular dynamics simulations to determine the phase behavior and critical point of propane molecular models

Molecular dynamics simulations were performed to determine two-phase configurations of model propane molecules below the critical point and in the near-critical, two-phase region. A postprocessor that uses a Monte Carlo method for determination of volumes attributable to each molecule was used to obtain density histograms of the particles from which the bulk coexisting equilibrium vapor and liquid densities were determined. This method of analyzing coexisting densities in a two-phase simulation is straightforward and can be easily implemented for complex, multisite models. Various degrees of internal flexibility in the propane models have little effect on the coexisting densities at temperatures 40 K or more below the critical point, but internal flexibility (angle bending and bond vibrations) does affect the saturated liquid densities in the near-critical region, changing the critical temperature by approximately 20 K. Shorter cutoffs were also found to affect the phase dome and the location of the critical point.     

Microstructure of concentrated solutions in the water-lithium nitrate-calcium nitrate ternary system from molecular dynamics simulations

By classical molecular dynamics, a series of solutions in the ternary water-lithium nitrate-calcium nitrate system is simulated. The radial distribution functions for different atom pairs are calculated, the coordination numbers are estimated, and the pattern of change of their first solvation shell composition with a change in the component ratio in the model system is examined. The first sphere of water molecule in the solutions with different concentration is considered in detail. Eight main structural types of the nearest-neighbor water molecule surrounding and their major features are identified, the probabilities of the occurrence of these structures in the ternary system are estimated.     

Ion condensation behavior and dynamics of water molecules surrounding the sodium poly(methacrylic acid) chain in water: a molecular dynamics study

All-atom molecular dynamics simulations are used to study the condensation behavior of monovalent (Na(+)) and multivalent (Ca(2+)) salt counterions associated with the co-ions (Cl(-)) surrounding the charged poly(methacrylic acid) (PMAA) chain in water. The study is extended to the influences on chain conformation, local arrangement, and dynamics of water in the highly diluted aqueous solutions. We find that even when the salt ions are monovalent, they attract more than one charged monomer and act as a bridging agent within the chain, as the multivalent salt ions. In principle, the salt ions bridge between not only the "non-adjacent" but also the "adjacent" charged monomers, leading to a more coil-like and a locally stretched conformation, respectively. With an increase in the salt concentration, the amount of coiled-type condensed ions increase and reach a maximum when the chain conformation becomes the most collapsed; whereas, the stretched-type shows an opposite trend. Our results show that the attractive interactions through the condensed salt ions between the non-adjacent monomers are responsible for the conformational collapse. When the salt concentration increases high enough, a significant increase for the stretched-type condensed ions makes an expansion effect on the chain. These stretched-type salt ions, followed by the adsorption of the co-ions and water molecules, tend to form a multilayer organization outside surrounding the PMAA chain. Thus, the expansion degree of the chain conformation is greatly limited. When only the monovalent Na(+) ions are present in the solutions, water molecules are primarily adsorbed into either the condensed Na(+) ions or the COO(-) groups. These adsorbed water molecules form hydrogen bonds with each other and enhance the local bridging behavior associated with the Na(+) condensation on the resultant chain conformation. With an increase in the amount of multivalent Ca(2+) salt ions, more water molecules are bonded directly with the condensed Ca(2+) ions. In this case, only the condensed Ca(2+) ions provide a strong bridging effect within the polymer chain. We observe a significant shift towards a higher frequency of the oxygen vibration spectrum and only a slight shift towards a higher frequency of the hydrogen spectrum for the water molecules associated with the ion condensation.     

Dissipative particle dynamics simulation of phase behavior of aerosol OT/water system

The phase behaviors of the binary mixture of an anionic surfactant aerosol OT (AOT) and water are investigated on a mesoscopic level using dissipative particle dynamics (DPD) computer simulations. With a simple surfactant model, various aggregation structures of AOT in water including the lamellar, viscous isotropic, and reverse hexagonal phases are obtained, which agree well with the experimental phase diagram. Special attention is given on the unusual lamellar regions. Water diffusivity shows much useful information to understand how the phase behaviors varied with concentration and temperature. It is proposed that the anomalous lamellar phenomena at intermediate AOT concentration (about 40%) are due to the formation of a defective structure, pseudoreversed hexagonal phase, which evidently decreases the water diffusivity. After increasing temperature above 328 K, the pseudoreversed hexagonal structure will be partly transformed to a normal lamellar phase structure and the system lamellar ordering is therefore enhanced.     

Elastic bag model for molecular dynamics simulations of solvated systems: application to liquid water and solvated peptides

The fluctuating elastic boundary (FEB) model for molecular dynamics has recently been developed and validated through simulations of liquid argon. In the FEB model, a flexible boundary which consists of particles connected by springs is used to confine the solvated system, thereby eliminating the need for periodic boundary conditions. In this study, we extend this model to the simulation of bulk water and solvated alanine dipeptide. Both the confining potential and boundary particle interaction functions are modified to preserve the structural integrity of the boundary and prevent the leakage of the solute-solvent system through the boundary. A broad spectrum of structural and dynamic properties of liquid water are computed and compared with those obtained from conventional periodic boundary condition simulations. The applicability of the model to biomolecular simulations is investigated through the analysis of conformational population distribution of solvated alanine dipeptide. In most cases we find remarkable agreement between the two simulation approaches.     

The [BMI][Tf2N] ionic liquid/water binary system: a molecular dynamics study of phase separation and of the liquid-liquid interface

We report molecular dynamics (MD) simulations of the aqueous interface of the hydrophobic [BMI][Tf2N] ionic liquid (IL), composed of 1-butyl-3-methylimidazolium cations (BMI+) and bis(trifluoromethylsulfonyl)imide anions (Tf2N-). The questions of water/IL phase separation and properties of the neat interface are addressed, comparing different liquid models (TIP3P vs TIP5P water and +1.0/-1.0 vs +0.9/-0.9 charged IL ions), the Ewald vs the reaction field treatments of the long range electrostatics, and different starting conditions. With the different models, the "randomly" mixed liquids separate much more slowly (in 20 to 40 ns) than classical water-oil mixtures do (typically, in less than 1 ns), finally leading to distinct nanoscopic phases separated by an interface, as in simulations which started with a preformed interface, but the IL phase is more humid. The final state of water in the IL thus depends on the protocol and relates to IL heterogeneities and viscosity. Water mainly fluctuates in hydrophilic basins (rich in O(Tf2N) and aromatic CH(BMI) groups), separated by more hydrophobic domains (rich in CF3(Tf2N) and alkyl(BMI) groups), in the form of monomers and dimers in the weakly humid IL phase, and as higher aggregates when the IL phase is more humid. There is more water in the IL than IL in water, to different extents, depending on the model. The interface is sharper and narrower (approximately 10 A) than with the less hydrophobic [BMI][PF6] IL and is overall neutral, with isotropically oriented molecules, as in the bulk phases. The results allow us to better understand the analogies and differences of aqueous interfaces with hydrophobic (but hygroscopic) ILs, compared to classical organic liquids.     

Thermal Gaussian molecular dynamics for quantum dynamics simulations of many-body systems: application to liquid para-hydrogen

A new method, here called thermal Gaussian molecular dynamics (TGMD), for simulating the dynamics of quantum many-body systems has recently been introduced [I. Georgescu and V. A. Mandelshtam, Phys. Rev. B 82, 094305 (2010)]. As in the centroid molecular dynamics (CMD), in TGMD the N-body quantum system is mapped to an N-body classical system. The associated both effective Hamiltonian and effective force are computed within the variational Gaussian wave-packet approximation. The TGMD is exact for the high-temperature limit, accurate for short times, and preserves the quantum canonical distribution. For a harmonic potential and any form of operator A?, it provides exact time correlation functions C(AB)(t) at least for the case of B, a linear combination of the position, x, and momentum, p, operators. While conceptually similar to CMD and other quantum molecular dynamics approaches, the great advantage of TGMD is its computational efficiency. We introduce the many-body implementation and demonstrate it on the benchmark problem of calculating the velocity time auto-correlation function for liquid para-hydrogen, using a system of up to N = 2592 particles.     

Molecular dynamics simulations of vapor/liquid coexistence using the nonpolarizable water models

The surface tension, vapor-liquid equilibrium densities, and equilibrium pressure for common water models were calculated using molecular dynamics simulations over temperatures ranging from the melting to the critical points. The TIP4P/2005 and TIP4P-i models produced better values for the surface tension than the other water models. We also examined the correlation of the data to scaling temperatures based on the critical and melting temperatures. The reduced temperature (T/T(c)) gives consistent equilibrium densities and pressure, and the shifted temperature T + (T(c, exp) - T(c, sim)) gives consistent surface tension among all models considered in this study. The modified fixed charge model which has the same Lennard-Jones parameters as the TIP4P-FQ model but uses an adjustable molecular dipole moment is also simulated to find the differences in the vapor-liquid coexistence properties between fixed and fluctuating charge models. The TIP4P-FQ model (2.72 Debye) gives the best estimate of the experimental surface tension. The equilibrium vapor density and pressure are unaffected by changes in the dipole moment as well as the surface tension and liquid density.     

Probing the neutral graphene-ionic liquid interface: insights from molecular dynamics simulations

We study basic mechanisms of the interfacial layer formation at the neutral graphite monolayer (graphene)-ionic liquid (1,3-dimethylimidazolium chloride, [dmim][Cl]) interface by fully atomistic molecular dynamics simulations. We probe the interface area by a spherical probe varying the charge (-1e, 0, +1e) as well as the size of the probe (diameter 0.50 nm and 0.38 nm). The molecular modelling results suggest that: there is a significant enrichment of ionic liquid cations at the surface. This cationic layer attracts Cl(-) anions that leads to the formation of several distinct ionic liquid layers at the surface. There is strong asymmetry in cationic/anionic probe interactions with the graphene wall due to the preferential adsorption of the ionic liquid cations at the graphene surface. The high density of ionic liquid cations at the interface adds an additional high energy barrier for the cationic probe to come to the wall compared to the anionic probe. Qualitatively the results from probes with diameter 0.50 nm and 0.38 nm are similar although the smaller probe can approach closer to the wall. We discuss the simulation results in light of available experimental data on the interfacial structure in ionic liquids.     

Density-dependent liquid nitromethane decomposition: molecular dynamics simulations based on ReaxFF

The decomposition mechanism of hot liquid nitromethane at various compressions was studied using reactive force field (ReaxFF) molecular dynamics simulations. A competition between two different initial thermal decomposition schemes is observed, depending on compression. At low densities, unimolecular C-N bond cleavage is the dominant route, producing CH(3) and NO(2) fragments. As density and pressure rise approaching the Chapman-Jouget detonation conditions (～30% compression, >2500 K) the dominant mechanism switches to the formation of the CH(3)NO fragment via H-transfer and/or N-O bond rupture. The change in the decomposition mechanism of hot liquid NM leads to a different kinetic and energetic behavior, as well as products distribution. The calculated density dependence of the enthalpy change correlates with the change in initial decomposition reaction mechanism. It can be used as a convenient and useful global parameter for the detection of reaction dynamics. Atomic averaged local diffusion coefficients are shown to be sensitive to the reactions dynamics, and can be used to distinguish between time periods where chemical reactions occur and diffusion-dominated, nonreactive time periods.     

Analyzing the molecular basis of enzyme stability in ethanol/water mixtures using molecular dynamics simulations

One of the drawbacks of nonaqueous enzymology is the fact that enzymes tend to be less stable in organic solvents than in water. There are, however, some enzymes that display very high stabilities in nonaqueous media. In order to take full advantage of the use of nonaqueous solvents in enzyme catalysis, it is essential to elucidate the molecular basis of enzyme stability in these media. Toward this end, we performed μs-long molecular dynamics simulations using two homologous proteases, pseudolysin, and thermolysin, which are known to have considerably different stabilities in solutions containing ethanol. The analysis of the simulations indicates that pseudolysin is more stable than thermolysin in ethanol/water mixtures and that the disulfide bridge between C30 and C58 is important for the stability of the former enzyme, which is consistent with previous experimental observations. Our results indicate that thermolysin has a higher tendency to interact with ethanol molecules (especially through van der Waals contacts) than pseudolysin, which can lead to the disruption of intraprotein hydrophobic interactions and ultimately result in protein unfolding. In the absence of the C30-C58 disulfide bridge, pseudolysin undergoes larger conformational changes, becoming more open and more permeable to ethanol molecules which accumulate in its interior and form hydrophobic interactions with the enzyme, destroying its structure. Our observations are not only in good agreement with several previous experimental findings on the stability of the enzymes studied in ethanol/water mixtures but also give an insight on the molecular determinants of this stability. Our findings may, therefore, be useful in the rational development of enzymes with increased stability in these media.     

Thermodynamic properties of methane/water interface predicted by molecular dynamics simulations

Molecular dynamics simulations have been performed to examine the thermodynamic properties of methane/water interface using two different water models, the TIP4P/2005 and SPC/E, and two sets of combining rules. The density profiles, interfacial tensions, surface excesses, surface pressures, and coexisting densities are calculated over a wide range of pressure conditions. The TIP4P/2005 water model was used, with an optimized combining rule between water and methane fit to the solubility, to provide good predictions of interfacial properties. The use of the infinite dilution approximation to calculate the surface excesses from the interfacial tensions is examined comparing the surface pressures obtained by different approaches. It is shown that both the change of methane solubilities in pressure and position of maximum methane density profile at the interface are independent of pressure up to about 2 MPa. We have also calculated the adsorption enthalpies and entropies to describe the temperature dependency of the adsorption.     

Mixing coarse-grained and fine-grained water in molecular dynamics simulations of a single system

The use of a supra-molecular coarse-grained (CG) model for liquid water as solvent in molecular dynamics simulations of biomolecules represented at the fine-grained (FG) atomic level of modelling may reduce the computational effort by one or two orders of magnitude. However, even if the pure FG model and the pure CG model represent the properties of the particular substance of interest rather well, their application in a hybrid FG/CG system containing varying ratios of FG versus CG particles is highly non-trivial, because it requires an appropriate balance between FG-FG, FG-CG, and CG-CG energies, and FG and CG entropies. Here, the properties of liquid water are used to calibrate the FG-CG interactions for the simple-point-charge water model at the FG level and a recently proposed supra-molecular water model at the CG level that represents five water molecules by one CG bead containing two interaction sites. Only two parameters are needed to reproduce different thermodynamic and dielectric properties of liquid water at physiological temperature and pressure for various mole fractions of CG water in FG water. The parametrisation strategy for the FG-CG interactions is simple and can be easily transferred to interactions between atomistic biomolecules and CG water.