Spectral diffusion in a fluctuating charge model of water

We apply the combined electronic structure/molecular dynamics approach of Corcelli, Lawrence, and Skinner [J. Chem. Phys. 120, 8107 (2004)] to the fluctuating charge (SPC-FQ) model of liquid water developed by Rick, Stuart, and Berne [J. Chem. Phys. 101, 6141 (1994)]. For HOD in H(2)O the time scale for the long-time decay of the OD stretch frequency time-correlation function, which corresponds to the time scale for hydrogen-bond rearrangement in the liquid, is about 1.5 ps. This result is significantly longer than the 0.9 ps decay previously calculated for the nonpolarizable SPC/E water model. Our results for the SPC-FQ model are in better agreement with recent vibrational echo experiments.     

Monte Carlo simulations of critical cluster sizes and nucleation rates of water

We have calculated the critical cluster sizes and homogeneous nucleation rates of water at temperatures and vapor densities corresponding to experiments by Wolk and Strey [J. Phys. Chem B 105, 11683 (2001)]. The calculations have been done with an expanded version of a Monte Carlo method originally developed by Vehkamaki and Ford [J. Chem. Phys. 112, 4193 (2000)]. Their method calculates the statistical growth and decay probabilities of molecular clusters. We have derived a connection between these probabilities and kinetic condensation and evaporation rates, and introduce a new way for the calculation of the work of formation of clusters. Three different interaction potential models of water have been used in the simulations. These include the unpolarizable SPC/E [J. Phys. Chem. 91, 6269 (1987)] and TIP4P [J. Chem. Phys. 79, 926 (1983)] models and a polarizable model by Guillot and Guissani [J. Chem. Phys. 114, 6720 (2001)]. We show that TIP4P produces critical cluster sizes and a temperature and vapor density dependence for the nucleation rate that agree well with the experimental data, although the magnitude of nucleation rate is constantly overestimated by a factor of 2 x 10(4). Guissani and Guillot's model is somewhat less successful, but both the TIP4P and Guillot and Guissani models are able to reproduce a much better experimental temperature dependency of the nucleation rate than the classical nucleation theory. Using SPC/E results in dramatically too small critical clusters and high nucleation rates. The water models give different average binding energies for clusters. We show that stronger binding between cluster molecules suppresses the decay probability of a cluster, while the growth probability is not affected. This explains the differences in results from different water models.     

A molecular dynamics study of water nucleation using the TIP4P/2005 model

Extensive molecular dynamics simulations were conducted using the TIP4P/2005 water model of Abascal and Vega [J. Chem. Phys. 123, 234505 (2005)] to investigate its condensation from supersaturated vapor to liquid at 330 K. The mean first passage time method [J. Wedekind, R. Strey, and D. Reguera, J. Chem. Phys. 126, 134103 (2007); L. S. Bartell and D. T. Wu, 125, 194503 (2006)] was used to analyze the influence of finite size effects, thermostats, and charged species on the nucleation dynamics. We find that the Nose?-Hoover thermostat and the one proposed by Bussi et al. [J. Chem. Phys. 126, 014101 (2007)] give essentially the same averages. We identify the maximum thermostat coupling time to guarantee proper thermostating for these simulations. The presence of charged species has a dramatic impact on the dynamics, inducing a marked change towards a pure growth regime, which highlights the importance of ions in the formation of liquid droplets in the atmosphere. It was found a small but noticeable sign preference at intermediate cluster sizes (between 5 and 30 water molecules) corresponding mostly to the formation of the second solvation shell around the ion. The TIP4P/2005 water model predicts that anions induce faster formation of water clusters than cations of the same magnitude of charge.     

Quantum Monte Carlo calculations of the dissociation energy of the water dimer

We report diffusion quantum Monte Carlo (DMC) calculations of the equilibrium dissociation energy D(e) of the water dimer. The dissociation energy measured experimentally, D(0), can be estimated from D(e) by adding a correction for vibrational effects. Using the measured dissociation energy and the modern value of the vibrational energy Mas et al., [J. Chem. Phys. 113, 6687 (2000)] leads to D(e)=5.00+/-0.7 kcal mol(-1), although the result Curtiss et al., [J. Chem. Phys. 71, 2703 (1979)] D(e)=5.44+/-0.7 kcal mol(-1), which uses an earlier estimate of the vibrational energy, has been widely quoted. High-level coupled cluster calculations Klopper et al., [Phys. Chem. Chem. Phys. 2, 2227 (2000)] have yielded D(e)=5.02+/-0.05 kcal mol(-1). In an attempt to shed new light on this old problem, we have performed all-electron DMC calculations on the water monomer and dimer using Slater-Jastrow wave functions with both Hartree-Fock approximation (HF) and B3LYP density functional theory single-particle orbitals. We obtain equilibrium dissociation energies for the dimer of 5.02+/-0.18 kcal mol(-1) (HF orbitals) and 5.21+/-0.18 kcal mol(-1) (B3LYP orbitals), in good agreement with the coupled cluster results.     

Semiempirical self-consistent polarization description of bulk water, the liquid-vapor interface, and cubic ice

We have applied an efficient electronic structure approach, the semiempirical self-consistent polarization neglect of diatomic differential overlap (SCP-NDDO) method, previously parametrized to reproduce properties of water clusters by Chang, Schenter, and Garrett [ J. Chem. Phys. 2008 , 128 , 164111 ] and now implemented in the CP2K package, to model ambient liquid water at 300 K (both the bulk and the liquid-vapor interface) and cubic ice at 15 and 250 K. The SCP-NDDO potential retains its transferability and good performance across the full range of conditions encountered in the clusters and the bulk phases of water. In particular, we obtain good results for the density, radial distribution functions, enthalpy of vaporization, self-diffusion coefficient, molecular dipole moment distribution, and hydrogen bond populations, in comparison to experimental measurements.     

Effect of nitrogen adsorption on the mid-infrared spectrum of water clusters

Experimental Fourier-transform infrared spectra and DFT calculated infrared spectra are compared to investigate the effect of adsorbed nitrogen on the OH-stretch band complex of water clusters. Using a collisional cooling experiment, pure as well as partially and completely N(2)-covered water clusters consisting of 20-200 water molecules have been generated in thermal equilibrium in the aerosol phase within the temperature range of 5-80 K. Computational IR-spectra simulations have been performed for discrete pure and N(2)-covered water clusters including 10, 15, 20, and 30 water molecules. The adsorbed N(2) molecules especially affect the three-coordinated water molecules at the cluster surface which could be observed as a blue shift of the companion O-H band at 2900 cm(-1) and a red shift of the dangling O-H band at 3700 cm(-1) by about 20 cm(-1) in both cases. The most striking effect of the N(2) adsorbate is an intensity increase of the dangling O-H band by a factor of 3-5. Furthermore, the onset temperature of nitrogen adsorption at the water cluster surface was experimentally found to be roughly 30 K for cluster sizes of about 100 water molecules. Experimental and computational results are in good agreement. The presented results are based on and support the work of V. Buch, J. P. Devlin, and co-workers (e.g., J. Phys. Chem. B, 1997; J. Phys. Chem. A, 2003; Int. Rev. Phys. Chem., 2004).     

One-dimensional model for water and aqueous solutions. I. Pure liquid water

Two simplified one-dimensional models for waterlike particles are studied. One is referred to as the primitive model which is a simplified version of a model introduced by Ben-Naim in 1992 [Statistical Thermodynamics for Chemists and Biochemists (Plenum, New York, 1992)]. The second, referred to as the primitive cluster model, is a simplified version of the model used by Lovett and Ben-Naim in 1969 [J. Chem. Phys. 51, 3108 (1969)]. The two models are shown to be nearly equivalent and both exhibit some of the most characteristic behavior of liquid water. It is argued that a key feature of the molecular interactions--the correlation between the strong binding energy and low local density--is essential for the manifestation of the anomalous behavior of liquid water. It is also essential for the understanding of the outstanding behavior of liquid water.     

Thermodynamic anomalies in a lattice model of water: solvation properties

We investigate a lattice-fluid model of water, defined on a three-dimensional body-centered-cubic lattice. Model molecules possess a tetrahedral symmetry, with four equivalent bonding arms. The model is similar to the one proposed by Roberts and Debenedetti [J. Chem. Phys. 105, 658 (1996)], simplified by removing distinction between "donors" and "acceptors." We focus on the solvation properties, mainly as far as an ideally inert (hydrophobic) solute is concerned. As in our previous analysis, devoted to neat water [J. Chem. Phys. 121, 11856 (2004)], we make use of a generalized first-order approximation on a tetrahedral cluster. We show that the model exhibits quite a coherent picture of water thermodynamics, reproducing qualitatively several anomalous properties observed both in pure water and in solutions of hydrophobic solutes. As far as supercooled liquid water is concerned, the model is consistent with the second critical-point scenario.     

Computer simulation study of water using a fluctuating charge model

Hydrogen bonding in small water clusters is studied through computer simulation methods using a sophisticated, empirical model of interaction developed by Ricket al (S W Rick, S J Stuart and B J Berne 1994J. Chem. Phys. 101 6141) and others. The model allows for the charges on the interacting sites to fluctuate as a function of time, depending on their local environment. The charge flow is driven by the difference in the electronegativity of the atoms within the water molecule, thus effectively mimicking the effects of polarization of the charge density. The potential model is thus transferable across all phases of water. Using this model, we have obtained the minimum energy structures of water clusters up to a size often. The cluster structures agree well with experimental data. In addition, we are able to distinctly identify the hydrogens that form hydrogen bonds based on their charges alone, a feature that is not possible in simulations using fixed charge models. We have also studied the structure of liquid water at ambient conditions using this fluctuating charge model.     

Ab initio potential and dipole moment surfaces for water. II. Local-monomer calculations of the infrared spectra of water clusters

We employ recent flexible ab initio potential energy and dipole surfaces [Y. Wang, X. Huang, B. C. Shepler, B. J. Braams, and J. M. Bowman, J. Chem. Phys. 134, 094509 (2011)] to the calculation of IR spectra of the intramolecular modes of water clusters. We use a quantum approach that begins with a partitioned normal-mode analysis of perturbed monomers, and then obtains solutions of the corresponding Schro?dinger equations for the fully coupled intramolecular modes of each perturbed monomer. For water clusters, these modes are the two stretches and the bend. This approach is tested against benchmark calculations for the water dimer and trimer and then applied to the water clusters (H(2)O)(n) for n = 6-10 and n = 20. Comparisons of the spectra are made with previous ab initio harmonic and empirical potential calculations and available experiments.     

The photodissociation of the water dimer in the A band: a twelve-dimensional quasiclassical study

The quasiclassical absorption spectrum of the water dimer in the A band was calculated taking into account motion in all degrees of freedom of the system. The ab initio excited state potentials employed were interpolated by the modified Shepard interpolation method using QMRCI energies and state-averaged MCSCF gradients and Hessians. The ground state vibrational wavefunction was variationally calculated using an adiabatic separation between the high and low frequency normal modes of the system. The calculated spectrum of water dimer shows a clear blueshift with respect to the monomer, but also a small red tail, in agreement with the prediction by Harvey et al. [J. Chem. Phys. 109, 8747 (1998)]. Previous three-dimensional model studies of the photodissociation of the water dimer by Valenzano et al. [J. Chem. Phys. 123, 034303 (2005)] did not show this red tail. A thorough analysis of the dependence of the spectrum on the modes coupled explicitly in the calculation of the spectrum shows that the red tail is due to coupling between the intramolecular stretch vibrations on different monomers.     

Vibronic coupling in asymmetric bichromophores: Theory and application to diphenylmethane

The theory for modeling vibronic interactions in bichromophores was introduced in sixties by Witkowski and Moffitt [J. Chem. Phys. 33, 872 (1960)] and extended by Fulton and Gouterman [J. Chem. Phys. 35, 1059 (1961)]. The present work describes extension of this vibronic model to describe bichromophores with broken vibrational symmetry such as partly deuterated molecules. Additionally, the model is extended to include inter-chromophore vibrational modes. The model can treat multiple vibrational modes by employing Lanczos diagonalization procedure of sparse matrices. The developed vibronic model is applied to simulation of vibronic spectra of flexible bichromophore diphenylmethane and compared to high-resolution experimental spectra [J. A. Stearns, N. R. Pillsbury, K. O. Douglass, C. W. Mu?ller, T. S. Zwier, and D. F. Plusquellic, J. Chem. Phys. 129, 224305 (2008)].     

Using heat capacity and compressibility to choose among two-state models of liquid water

We extend to heat capacity Cp the model of Vedamuthu, Singh, and Robinson (J. Phys. Chem. 1994, 98, 2222). This model and that of Bartell (J. Phys. Chem. B 1997, 101, 7573) fit successfully, even in the supercooled region, the temperature dependence of Cp, volume, and isothermal compressibility kappa(T). The Robinson model is superior for kappa(T). Tanaka's model (J. Chem. Phys. 2000, 112, 799) fails for C(p) even after correction of a derivational error. All three models assume that the liquid consists of low-density component 1 and high-density component 2. We conclude that Robinson's tactics, ignoring the intercomponent equilibrium constant and determining compositions solely from volumes, yield the most reliable compositions and individual-component properties. Our fits of the Robinson model to C(p) yield at 0 degrees C H(2) - H(1) of (135 +/- 35) J/g, H(1) - H(ice) of 0.8DeltaH(fus), and C(2) - C(1) of (0.1 +/- 0.7) J/K.g. The enthalpy difference between the components is largely responsible for the rapid change of C(p) at the lowest supercooled temperatures. We propose an adjustment to Speedy and Angell's (J. Chem. Phys. 1976, 65, 851) experimental values of kappa(T) for supercooled water.     

Insights in quantum dynamical effects in the infrared spectroscopy of liquid water from a semiclassical study with an ab initio-based flexible and polarizable force field

The dynamical properties of liquid water play an important role in many processes in nature. In this paper, we focus on the infrared (IR) absorption spectrum of liquid water based on the linearized semiclassical initial value representation (LSC-IVR) with the local Gaussian approximation (LGA) [J. Liu and W. H. Miller, J. Chem. Phys. 131, 074113 (2009)] and an ab initio based, flexible, polarizable Thole-type model (TTM3-F) [G. S. Fanourgakis and S. S. Xantheas, J. Chem. Phys. 128, 074506 (2008)]. Although the LSC-IVR (LGA) gives the exact result for the isolated three-dimensional shifted harmonic stretching model, it yields a blueshifted peak position for the more realistic anharmonic stretching potential. By using the short-time information of the LSC-IVR correlation function; however, it is shown how one can obtain more accurate results for the position of the stretching peak. Due to the physical decay in the condensed phase system, the LSC-IVR (LGA) is a good and practical approximate quantum approach for the IR spectrum of liquid water. The present results offer valuable insight into future attempts to improve the accuracy of the TTM3-F potential or other ab initio-based models in reproducing the IR spectrum of liquid water.     

Quantum effects in liquid water and ice: model dependence

This paper explores the influence of choice of potential model on the quantum effects observed in liquid water and ice. This study utilizes standard rigid models and a more formal context for the rigid-body centroid molecular dynamics methodology used to perform the quantum simulations is provided. Quantum and classical molecular dynamics simulations are carried out for liquid water and ice Ih at 298 and 220 K, respectively, with the simple point charge/extended and TIP4P-Ew water models. The results obtained for equilibrium and dynamical properties are compared with those recently reported on TIP4P [L. Hernandez de la Pena and P. G. Kusalik, J. Chem. Phys. 121, 5992 (2004); L. Hernandez de la Pena et al., J. Chem. Phys 123, 144506 (2005)]. For the liquid, an energy shift of about 8% and an average molecular uncertainty of about 11 degrees were found independently of the water model. The self-diffusion coefficient consistently increases by more than 50% when going from the classical to the quantum system and quantum dynamics are found to reproduce the experimental isotopic shifts with the models examined. The ice results compare remarkably well with those previously reported for the TIP4P water model; they confirm that quantum effects are considerable and that the quantum mechanical uncertainty and the energy shifts due to quantization are smaller in ice than in liquid water. The relevance of these findings in the context of the construction of water models is briefly discussed.     

A nine-dimensional calculation of the vibrational OH stretching and HOH bending spectrum of the water trimer

We have studied the vibrational high-frequency spectrum of the water trimer computationally. We expand an earlier study [J. Chem. Phys. A 2009, 113, 9124-9132] where we approximated the water trimer as three individually vibrating water monomer units. Some intramolecular potential energy coupling terms are now included in the previous model. The six OH bond lengths and the three HOH bending angles are used as the internal coordinates. The kinetic energy operator is a sum of the kinetic energy operators of the monomer units. We use the coupled cluster method with single, double, and perturbative triple excitations method [CCSD(T)] with augmented correlation consistent polarized valence triple-ζ (aug-cc-pVTZ) basis set to calculate the potential energy surface (PES). The counterpoise correction is included in the one-dimensional part of the PES. We calculate the vibrational energy eigenvalues using the variational method. The corresponding eigenfunctions are used to obtain the absorption intensities.     

Ultrafast vibrational spectroscopy of water and aqueous N-methylacetamide: Comparison of different electronic structure/molecular dynamics approaches

Kwac and Cho [J. Chem. Phys. 119, 2247 (2003)] have recently developed a combined electronic structure/molecular dynamics approach to vibrational spectroscopy in liquids. The method involves fitting ab initio vibrational frequencies for a solute in a cluster of solvent molecules to a linear combination of the electrostatic potentials on the solute atoms due to the charges on the solvent molecules. These authors applied their method to the N-methylacetamide-D/D(2)O system. We (S. A. Corcelli, C. P. Lawrence, and J. L. Skinner, [J. Chem. Phys. 120, 8107 (2004)]) have recently explored a closely related method, where instead of the electrostatic potential, the solute vibrational frequencies are fit to the components of the electric fields on the solute atoms due to the solvent molecules. We applied our method to the HOD/D(2)O and HOD/H(2)O systems. In order to make a direct comparison of these two approaches, in this paper we apply their method to the water system, and our method to the N-methylacetamide system. For the water system we find that the electric field method is superior to the potential approach, as judged by comparison with experiments for the absorption line shape. For the N-methylacetamide system the two methods are comparable.     

Soft sticky dipole-quadrupole-octupole potential energy function for liquid water: an approximate moment expansion

A new, efficient potential energy function for liquid water is presented here. The new model, which is referred here as the soft sticky dipole-quadrupole-octupole (SSDQO) model, describes a water molecule as a Lennard-Jones sphere with point dipole, quadrupole, and octupole moments. It is a single-point model and resembles the hard-sphere sticky dipole potential model for water by Bratko et al. [J. Chem. Phys. 83, 6367 (1985)] and the soft sticky dipole model by Ichiye and Liu [J. Phys. Chem. 100, 2723 (1996)] except now the sticky potential consists of an approximate moment expansion for the dimer interaction potential, which is much faster than the true moment expansion. The object here is to demonstrate that the SSDQO potential energy function can accurately mimic the potential energy function of a multipoint model using the moments of that model. First, the SSDQO potential energy function using the dipole, quadruple, and octupole moments from SPC/E, TIP3P, or TIP5P is shown to reproduce the dimer potential energy functions of the respective multipoint model. In addition, in Monte Carlo simulations of the pure liquid at room temperature, SSDQO reproduces radial distribution functions of the respective model. However, the Monte Carlo simulations using the SSDQO model are about three times faster than those using the three-point models and the long-range interactions decay faster for SSDQO (1/r(3) and faster) than for multipoint models (1/r). Moreover, the contribution of each moment to the energetics and other properties can be determined. Overall, the simplicity, efficiency, and accuracy of the SSDQO potential energy function make it potentially very useful for studies of aqueous solvation by computer simulations.