Spectral characterization of hydrogen bond network dynamics in water

Some computational aspects of the characterization of the complex hydrogen bond network dynamics using power spectral analysis are discussed. In the case of hydrogen-bonded liquids, the tagged molecule potential energy is shown to be a useful quantity for capturing the behavior of the networked liquid on different lengths and time scales. The computation of the tagged potential energy for rigid-body effective pair potentials, such as the TIP5P-E and SPC-E models, is discussed. The more structured nature of the TIP5P-E potential, compared to the SPC/E potential, shows up as differences in the high-frequency librational band of the power spectra of the tagged molecule potential energies. The static distributions of the tagged molecule potential energies are also more structured in the case of TIP5P-E, rather than SPC/E, water. The overall behavior of the key power spectral features remains the same in both the models. The possibility of detailed characterization of the power spectrum, and therefore of the underlying dynamics, using a model-based parametric fitting procedure for the power spectra is also discussed. We show that a parametric fitting can allow one to test alternative models of the dynamics underlying the liquid state dynamics.     

Effect of ionic solutes on the hydrogen bond network dynamics of water: power spectral analysis of aqueous NaCl solutions

To understand the modifications of the hydrogen bond network of water by ionic solutes, power spectra as well as static distributions of the potential energies of tagged solvent molecules and solute ions have been computed from molecular dynamics simulations of aqueous NaCl solutions. The key power spectral features of interest are the presence of high-frequency peaks due to localized vibrational modes, the existence of a multiple time scale or 1/falpha frequency regime characteristic of networked liquids, and the frequency of crossover from 1/falpha type behavior to white noise. Hydrophilic solutes, such as the sodium cation and the chloride anion, are shown to mirror the multiple time scale behavior of the hydrogen bond network fluctuations, unlike hydrophobic solutes which display essentially white noise spectra. While the power spectra associated with tagged H2O molecules are not very sensitive to concentration in the intermediate frequency 1/falpha regime, the crossover to white noise is shifted to lower frequencies on going from pure solvent to aqueous alkali halide solutions. This suggests that new and relatively slow time scales enter the picture, possibly associated with processes such as migration of water molecules from the hydration shell to the bulk or conversion of contact ion pairs into solvent-separated ion pairs which translate into variations in equilibrium transport properties of salt solutions with concentration. For anions, cations, and solvent molecules, the trends in the alpha exponents of the multiple time scale region and the self-diffusivities are found to be strongly correlated.     

Diffusional anomaly and network dynamics in liquid silica

The present study applies the power spectral analysis technique to understand the diffusional anomaly in liquid silica, modeled using the Beest-Kramer-van Santen (BKS) potential. Molecular-dynamics simulations have been carried out to show that power spectrum of tagged particle potential energy of silica shows a regime with 1f(alpha) dependence on frequency f which is the characteristic signature of multiple time scale behaviour in networks. As demonstrated earlier in the case of water [J. Chem. Phys. 122, 104507 (2005)], the variations in the mobility associated with the diffusional anomaly are mirrored in the scaling exponent alpha associated with this multiple time scale behavior. Our results indicate that in the anomalous regime, as the local tetrahedral order decreases with temperature or pressure, the coupling of local modes to network reorganizations increases and so does the diffusivity. This symmetry-dependence of the vibrational couplings is responsible for the connection between the structural and diffusional anomalies.     

Measuring diffusivity in supercooled liquid nanoscale films using inert gas permeation. I. Kinetic model and scaling methods

We describe in detail a diffusion model used to simulate inert gas transport through supercooled liquid overlayers. In recent work, the transport of the inert gas has been shown to be an effective probe of the diffusivity of supercooled liquid methanol in the experimentally challenging regime near the glass transition temperature. The model simulations accurately and quantitatively describe the inert gas permeation desorption spectra. The simulation results are used to validate universal scaling relationships between the diffusivity, overlayer thickness, and the temperature ramp rate for isothermal and temperature programmed desorption. From these scaling relationships we derive simple equations from which the diffusivity can be obtained using the peak desorption time or temperature for an isothermal or set of TPD experiments, respectively, without numerical simulation. The results presented here demonstrate that the permeation of gases through amorphous overlayers has the potential to be a powerful technique to obtain diffusivity data in deeply supercooled liquids.     

Spectral signatures of the diffusional anomaly in water

Power spectra for various tagged particle quantities in bulk extended simple point charge model water [H. J. C. Berendsen, J. R. Grigera, and T. P. Straatsma, J. Phys. Chem. 91, 6269 (1987)] are shown to have a regime with 1f(alpha) dependence on frequency f with alpha lying between 1 and 1.5 if the dynamical changes in the particular observable are sensitive to the multiple time-scale behavior of the hydrogen-bond network. The variations in mobility associated with the diffusional anomaly are mirrored in the scaling exponent alpha associated with this multiple time-scale behavior, suggesting that monitoring of 1f(alpha) behavior is a simple and direct method for linking phenomena on three distinctive length and time scales: the local molecular environment, hydrogen-bond network reorganizations, and the diffusivity. Our results indicate that experimental studies of supercooled water to probe the density dependence of 1f(alpha) spectral features, or equivalent stretched exponential behavior in time-correlation functions, will be of interest.     

Measuring diffusivity in supercooled liquid nanoscale films using inert gas permeation. II. Diffusion of Ar, Kr, Xe, and CH4 through methanol

We present an experimental technique to measure the diffusivity of supercooled liquids at temperatures near their T(g). The approach uses the permeation of inert gases through supercooled liquid overlayers as a measure of the diffusivity of the supercooled liquid itself. The desorption spectra of the probe gas are used to extract the low temperature supercooled liquid diffusivities. In the preceding companion paper, we derived equations using ideal model simulations from which the diffusivity could be extracted using the desorption peak times for isothermal or peak temperatures for temperature programmed desorption experiments. Here, we discuss the experimental conditions for which these equations are valid and demonstrate their utility using amorphous methanol with Ar, Kr, Xe, and CH(4) as probe gases. The approach offers a new method by which the diffusivities of supercooled liquids can be measured in the experimentally challenging temperature regime near the glass transition temperature.     

Ionic melts with waterlike anomalies: thermodynamic properties of liquid BeF2

Thermodynamic properties of liquid beryllium difluoride (BeF(2)) are studied using canonical ensemble molecular dynamics simulations of the transferable rigid ion model potential. The negative slope of the locus of points of maximum density in the temperature-pressure plane is mapped out. The excess entropy, computed within the pair correlation approximation, is found to show an anomalous increase with isothermal compression at low temperatures which will lead to diffusional as well as structural anomalies resembling those in water. The anomalous behavior of the entropy is largely connected with the behavior of the Be-F pair correlation function. The internal energy shows a T(35) temperature dependence. The pair correlation entropy shows a T(-25) temperature dependence only at high densities and temperatures. The correlation plots between internal energy and the pair correlation entropy for isothermal compression show the characteristic features expected of network-forming liquids with waterlike anomalies. The tagged particle potential energy distributions are shown to have a multimodal form at low temperatures and densities similar to those seen in other liquids with three-dimensional tetrahedral networks, such as water and silica.     

Structural correlations and cooperative dynamics in supercooled liquids

The relationships between diffusivity and the excess, pair and residual multiparticle contributions to the entropy are examined for Lennard-Jones liquids and binary glassformers, in the context of approximate inverse power law mappings of simple liquids. In the dense liquid where diffusivities are controlled by collisions and cage relaxations, Rosenfeld-type excess entropy scaling of diffusivities is found to hold for both crystallizing as well as vitrifying liquids. The crucial differences between the two categories of liquids emerge only when local cooperative effects in the dynamics result in significant caging effects in the time-dependent behaviour of the single-particle mean square displacement. In the case of glassformers, onset of such local cooperativity coincides with onset of deviations from Rosenfeld-type excess entropy scaling of diffusivities and increasing spatiotemporal heterogeneity. In contrast, for two- and three-dimensional liquids with a propensity to crystallise, the onset of local cooperative dynamics is sufficient to trigger crystallization provided that the liquid is sufficiently supercooled that the free energy barrier to nucleation of the solid phase is negligible. The state points corresponding to onset of transient caging effects can be associated with typical values, within reasonable bounds, of the excess, pair, and residual multiparticle entropy as a consequence of the isomorph-invariant character of the excess entropy, diffusivity and related static and dynamic correlation functions.     

Solvation of perfluorooctane and octane in water, methanol, acetonitrile, and aqueous mixtures of methanol and acetonitrile

Molecular dynamics simulations are used to examine the local solvation structure of single octane and perfluorooctane molecules in liquid water, methanol, acetonitrile, and aqueous mixtures of methanol and acetonitrile. The motivation is to obtain baseline information about the solvation of perfluorooctane by liquids used as the mobile phase in liquid chromatography and how it differs from the solvation of octane. While octane is uniformly solvated by both water and the second component, perfluorooctane is solvated by methanol and acetonitrile with the exclusion of water from the first solvation layer when the solvent is a mixture.     

Does Marcus theory apply to redox processes in ionic liquids? A simulation study

Simulation studies on a model system of a spherical ion with various charges in two imidazolium ionic liquids and in acetonitrile are compared. The average vertical ionisation potentials as a function of the charge on the ion are similar for all three systems. The Landau free energies of each system as a function of the vertical ionisation potential are computed and are close to being parabolic. Results are shown for the solvent reorganisation energies and for the activation free energies. The similarities of all these quantities are interpreted in terms of continuum models. However, the dynamics are likely to be very different in a polar liquid and in an ionic liquid as in the former case screening occurs by reorientation of molecules and in the latter case it occurs by translation of ions.     

Estimation of infinite dilution activity coefficients in aqueous mixtures via molecular simulation

Molecular dynamics simulations are employed to calculate infinite dilution activity coefficients of water and methanol-like species in binary mixtures using a variant of the Kirkwood coupling parameter method. Differences in residual Gibbs free energies are obtained as integrals over ensemble averages of the derivatives of total potential energy with respect to simple functions of the intramolecular potential energy parameters. The calculated limiting activity coefficients are compared with the experimental values at the same temperature obtained by direct measurement of the water/methanol binary and via extrapolation from vapor–liquid equilibria data.     

Mutual diffusion in binary mixtures of ionic liquids and molecular liquids by dynamic light scattering (DLS)

The present study shows that dynamic light scattering (DLS) is capable of measuring mutual diffusion coefficients for binary mixtures of ionic liquids (ILs) with different molecular liquids over the complete composition range. Evidence is given that the light scattering signals are related to true molecular binary diffusion. The method stands out due to its ability to work non-invasively in macroscopic thermodynamic equilibrium with reasonable accuracy and within convenient measurement periods. Compared with other techniques, mixtures with distinctly higher viscosities can be probed. For exemplary binary mixtures of 1-ethyl-3-methylimidazolium ethylsulfate ([EMIM][EtSO(4)]) with acetone, acetonitrile, dichloromethane, ethanol, or water as well as of 1-ethyl-3-methylimidazolium methanesulfonate ([EMIM][MeSO(3)]) with acetone, water, or methanol, mutual diffusivity data were measured over a wide range of composition at a temperature of 293.15 K. In general, the mutual diffusivity increases with increasing mole fraction of the molecular liquid and similarities to aqueous solutions of classical inorganic salts can be found. The characteristic behavior of the mutual diffusion coefficients is influenced by the nature of the chosen molecular liquid. For IL water mixtures, low light scattering intensities were observed despite the large refractive index difference of the pure components. The reason for this behavior may be the existence of water clusters in the mixtures. Additional measurements for IL acetone mixtures at temperatures ranging from 278.15 K to 323.15 K showed that the temperature dependence of the mutual diffusivity can be represented by Arrhenius functions and is increasing for decreasing mole fractions of acetone.     

First-principles calculation of thermodynamic stability of acids and bases under pH environment: a microscopic pH theory

Despite being one of the most important thermodynamic variables, pH has yet to be incorporated into first-principles thermodynamics to calculate stability of acidic and basic solutes in aqueous solutions. By treating the solutes as defects in homogeneous liquids, we formulate a first-principles approach to calculate their formation energies under proton chemical potential, or pH, based on explicit molecular dynamics. The method draws analogy to first-principle calculations of defect formation energies under electron chemical potential, or Fermi energy, in semiconductors. From this, we propose a simple pictorial representation of the general theory of acid-base chemistry. By performing first-principles molecular dynamics of liquid water models with solutes, we apply the formulation to calculate formation energies of various neutral and charged solutes such as H(+), OH(-), NH(3), NH(4)(+), HCOOH, and HCOO(-) in water. The deduced auto-dissociation constant of water and the difference in the pKa values of NH(3) and HCOOH show good agreement with known experimental values. Our first-principles approach can be further extended and applied to other bio- and electro-chemical molecules such as amino acids and redox reaction couples that could exist in aqueous environments to understand their thermodynamic stability.     

Solvation and solvatochromism in CO2-expanded liquids. 1. Simulations of the solvent systems CO2 + cyclohexane, acetonitrile, and methanol

Molecular dynamics simulations of CO(2)-expanded cyclohexane, acetonitrile, and methanol are reported at various compositions along the experimental bubble-point curve at 298 K. Simulated properties include energies, local compositions, viscosities, diffusion coefficients, and dielectric constants and relaxation times. On the basis of the limited comparisons to experimental data currently available, the results indicate that simple intermolecular potential models previously developed for simulating the pure components provide reasonable representations of the energetics and dynamics of these gas-expanded liquids.     

Simulation studies of ionic liquids: orientational correlations and static dielectric properties

The ionic liquids BMIM+I-, BMIM+BF4-, and BMIM+PF6- were simulated by means of the molecular dynamics method over a time period of more than 100 ns. Besides the common structural analysis, e.g., radial distribution functions and three dimensional occupancy plots, a more sophisticated orientational analysis was performed. The angular correlation functions g(00)110(r) and g(00)101(r) are the first distance dependent coefficients of the pairwise orientational distribution function g(rij,Omega1,Omega2,Omega12). These functions help to interpret the three dimensional plot and reveal interesting insights into the local structure of the analyzed ionic liquids. Furthermore, the collective network of ionic liquids can be characterized by the Kirkwood factor Gkappa(r) [J. Chem. Phys. 7, 911 (1939)]. The short-range behavior (r<10 A) of this factor may be suitable to predict the water miscibility of the ionic liquid. The long-range limit of Gkinfinity is below 1 which demonstrates the strongly coupled nature of the ionic liquid networks. In addition, this factor relates the orientational structure and the dielectric properties of the ionic liquids. The static dielectric constant epsilon(omega=0) for the simulated system is 8.9-9.5. Since in ionic liquids the very same molecule contributes to the total dipole moment as well as carries a net charge, a small, but significant contribution of the cross term between the total dipole moment and the electric current to epsilon(omega=0) is observed.     

Polarizable Atomic Multipole-based Molecular Mechanics for Organic Molecules

An empirical potential based on permanent atomic multipoles and atomic induced dipoles is reported for alkanes, alcohols, amines, sulfides, aldehydes, carboxylic acids, amides, aromatics and other small organic molecules. Permanent atomic multipole moments through quadrupole moments have been derived from gas phase ab initio molecular orbital calculations. The van der Waals parameters are obtained by fitting to gas phase homodimer QM energies and structures, as well as experimental densities and heats of vaporization of neat liquids. As a validation, the hydrogen bonding energies and structures of gas phase heterodimers with water are evaluated using the resulting potential. For 32 homo- and heterodimers, the association energy agrees with ab initio results to within 0.4 kcal/mol. The RMS deviation of hydrogen bond distance from QM optimized geometry is less than 0.06 ?. In addition, liquid self-diffusion and static dielectric constants computed from molecular dynamics simulation are consistent with experimental values. The force field is also used to compute the solvation free energy of 27 compounds not included in the parameterization process, with a RMS error of 0.69 kcal/mol. The results obtained in this study suggest the AMOEBA force field performs well across different environments and phases. The key algorithms involved in the electrostatic model and a protocol for developing parameters are detailed to facilitate extension to additional molecular systems.     

Understanding the interactions between tris(pentafluoroethyl)-trifluorophosphate-based ionic liquid and small molecules from molecular dynamics simulation

Tris(pentafluoroethyl)trifluorophosphate ([FEP])-based ionic liquids have been widely applied in many fields. For better understanding the properties of [FEP]-based ionic liquids, the interactions between 1-hexyl-3-methylimidazolium ([hmim])[FEP] and small molecules were investigated by molecular dynamics simulations in this work. The small molecules are water, methanol and dimethyl ether. The united-atom (UA) force fields were proposed for methanol and dimethyl ether based on AMBER force field. The densities, enthalpies of vaporization, excess molar properties, and diffusion coefficients of the mixtures were calculated, as well as the microscopic structures characterized by radial distribution functions. Both of the results of the excess energies and microscopic properties show that the strongest interaction is between [hmim][FEP] and dimethyl ether, whereas the interaction between [hmim][FEP] and water is the weakest. Moreover, [hmim][FEP] is more hydrophobic than [hmim] hexafluorophosphate ([PF6]), and the three solutes are mainly distributed around [FEP] anion.     

Relating composition, structural order, entropy and transport in multi-component molten salts

Molecular dynamics simulations of the LiF-BeF(2) molten salt mixture are used to establish relationships between composition, structural order, entropy, and transport properties of multi-component ionic liquids. A sharp rise in tetrahedral order associated with formation of the fluoroberyllate network occurs for compositions with BeF(2) concentrations greater than that of the Li(2)BeF(4)-BeF(2) eutectic. The excess entropy of the liquid in this regime, within the pair correlation approximation, is strongly correlated with the local tetrahedral order. The different degree of participation of beryllium, fluorine, and lithium ions in the cooperative dynamics of the fluoroberyllate network can be related to the degree of deviation from Rosenfeld-type excess entropy scaling, with the lithium ions remaining essentially unaffected by the liquid state network. We demonstrate that the deviations from Nernst-Einstein and Stokes-Einstein behaviour emerge only in temperature-composition regimes where tetrahedral order strongly correlates with the pair entropy. Implications for understanding structure-property relationships in other ionic liquids, such as molten salts, oxide melts, and RTILs are considered.     

Analysis of multiple front formation in the wetting of thin-layer plates

Multiple fronts were formed when pure liquids such as methanol or heptane and multi-component liquids such as ethanol-toluene mixtures were flowed through initially dry thin-layer chromatography plates, whether the stationary phase was naked silica or octadecyl-silica. A minor modification to a commercially available thin-layer chromatograph allowed the fronts to be clearly seen, and their dynamics to be measured experimentally. A model is proposed to explain their formation and motion.