Theory for the solvation of nonpolar solutes in water

We recently developed an angle-dependent Wertheim integral equation theory (IET) of the Mercedes-Benz (MB) model of pure water [Silverstein et al., J. Am. Chem. Soc. 120, 3166 (1998)]. Our approach treats explicitly the coupled orientational constraints within water molecules. The analytical theory offers the advantage of being less computationally expensive than Monte Carlo simulations by two orders of magnitude. Here we apply the angle-dependent IET to studying the hydrophobic effect, the transfer of a nonpolar solute into MB water. We find that the theory reproduces the Monte Carlo results qualitatively for cold water and quantitatively for hot water.     

Polarizable and flexible model for ethanol

A polarizable, flexible model for ethanol is obtained based on an extensive series of B3LYP/6-311++G(d,p) calculations and molecular dynamics simulations. The ethanol model includes electric-field dependence in both the atomic charges and the intramolecular degrees of freedom. Field-dependent intramolecular potentials have been attempted only once previously, for OH and HH stretches in water [P. Cicu et al., J. Chem. Phys. 112, 8267 (2000)]. The torsional potential involving the hydrogen-bonding hydrogen in ethanol is found to be particularly field sensitive. The methodology for developing field-dependent potentials can be readily generalized to other molecules and is discussed in detail. Molecular dynamics simulations of bulk ethanol are performed and the results are assessed based on comparisons with the self-diffusion coefficient [N. Karger et al., J. Chem. Phys. 93, 3437 (1990)], dielectric constant [J. T. Kindt and C. A. Schmuttenmaer, J. Phys. Chem. 100, 10373 (1996)], enthalpy of vaporization [R. C. Wilhoit and B. J. Zwolinski, J. Phys. Chem. Ref. Data, Suppl. 2, 2 (1973)], and experimental interatomic distributions [C. J. Benmore and Y. L. Loh, J. Chem. Phys. 112, 5877 (2000)]. The simultaneous variation of the atomic charges and the intramolecular potentials requires modified equations of motion and a multiple time step algorithm has been implemented to solve these equations. The article concludes with a discussion of the bulk structure and properties with an emphasis on the hydrogen bonding network.     

Molecular dynamics simulations of surface-initiated melting of nitromethane

The melting of nitromethane initiated at solid-vacuum interfaces has been investigated using molecular dynamics nvt simulations with a realistic force field [D. C. Sorescu et al., J. Phys. Chem. B 104, 8406 (2000)]. The calculated melting point (251+/-5 K) is in good agreement with experiment (244.73 K) and values obtained previously (approximately 255.5 and 266.5+/-8 K) using other simulation methods [P. M. Agrawal et al., J. Chem. Phys. 119, 9617 (2003)]. Analyses of the molecular orientations and diffusion during the simulations as functions of the distance from the exposed surfaces show that the melting is a direct crystal-to-liquid transition, in which the molecules first gain rotational freedom, then mobility. There is a slight dependence of the melting temperature on the exposed crystallographic face.     

Modeling collective behavior of molecules in nanoscale direct deposition processes

We present a theoretical model describing the collective behavior of molecules in nanoscale direct deposition processes such as dip-pen nanolithography. We show that strong intermolecular interactions combined with nonuniform substrate-molecule interactions can produce various shapes of molecular patterns including fractal-like structures. Computer simulations reveal circular and starlike patterns at low and intermediate densities of preferentially attractive surface sites, respectively. At large density of such surface sites, the molecules form a two-dimensional invasion percolation cluster. Previous experimental results showing anisotropic patterns of various chemical and biological molecules correspond to the starlike regime [P. Manandhar et al., Phys. Rev. Lett. 90, 115505 (2003); J.-H. Lim and C. A. Mirkin, Adv. Mater. (Weinheim, Ger.) 14, 1474 (2002); D. L. Wilson et al., Proc. Natl. Acad. Sci. U.S.A. 98, 13660 (2001); M. Su et al., Appl. Phys. Lett. 84, 4200 (2004); R. McKendry et al., Nano Lett. 2, 713 (2002); H. Zhou et al., Appl. Surf. Sci. 236, 18 (2004); G. Agarwal et al., J. Am. Chem. Soc. 125, 580 (2003)].     

Surface relaxation in water clusters: evidence from theoretical analysis of the oxygen 1s photoelectron spectrum

We present a theoretical interpretation of the oxygen 1s photoelectron spectrum published by Ohrwall et al. [J. Chem. Phys. 123, 054310 (2005)]. A water cluster that contains 200 molecules was simulated at 215 K using the polarizable AMOEBA force field. The force field predicts longer O...O distances at the cluster surface than in the bulk. Comparisons to ab initio molecular dynamics (MD) simulations indicate that the force field underestimates the degree of surface relaxation. By comparing cluster lineshape models, computed from MD simulations, to the experimental spectrum we find further evidence of surface relaxation.     

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.     

Computational study of structural and dynamical properties of formamide-water mixtures

A molecular dynamics simulation study of structural and dynamical properties in liquid mixtures of formamide and water is presented. Site-site radial pair distribution functions, local mole fractions, pair energy distributions, and tetrahedral orientational order are the quantities analyzed to investigate the local structure in the simulated mixtures, along with a review of the intermolecular structure in terms of the distribution of hydrogen bonds. Our results indicate that there is a substitution of formamide molecules by water in the hydrogen bonds and a formation of a common hydrogen bond network. By analyzing the extent of tetrahedral order in the liquid as a function of composition, it is observed that whereas the tetrahedral network of liquid water is progressively lost by increasing the formamide concentration, the water structure within the first coordination shell is preserved and somewhat enhanced. The hydrogen-bond mean lifetimes were estimated by performing a time integration of the autocorrelation functions of bond occupation numbers. The lifetimes associated with hydrogen bonds between water, formamide, and interspecies pairs are found to increase with increasing formamide concentration. The lifetimes of the water hydrogen bonds show the largest variations, supporting the picture of an enhancement of the water structure among the nearest neighbors within the first coordination shell. We have used two different force field models for water, SPC/E [J. C. Berendsen et al., J. Phys. Chem. 91, 6269 (1987)] and TIP4P/2005 [J. L. F. Abascal and C. Vega, J. Chem. Phys. 123, 234505 (2005)]. Our results for structural and dynamical properties yield very small differences between those models, the TIP4P/2005 predicting a slightly more structured liquid and, consequently, exhibiting a slightly slower translational and librational dynamics.     

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.     

Density, structure, and dynamics of water: the effect of van der Waals interactions

It is known that ab initio molecular dynamics (AIMD) simulations of liquid water at ambient conditions, based on the generalized gradient approximation (GGA) to density functional theory (DFT), with commonly used functionals fail to produce structural and diffusive properties in reasonable agreement with experiment. This is true for canonical, constant temperature simulations where the density of the liquid is fixed to the experimental density. The equilibrium density, at ambient conditions, of DFT water has recently been shown by Schmidt et al. [J. Phys. Chem. B, 113, 11959 (2009)] to be underestimated by different GGA functionals for exchange and correlation, and corrected by the addition of interatomic pair potentials to describe van der Waals (vdW) interactions. In this contribution we present a DFT-AIMD study of liquid water using several GGA functionals as well as the van der Waals density functional (vdW-DF) of Dion et al. [Phys. Rev. Lett. 92, 246401 (2004)]. As expected, we find that the density of water is grossly underestimated by GGA functionals. When a vdW-DF is used, the density improves drastically and the experimental diffusivity is reproduced without the need of thermal corrections. We analyze the origin of the density differences between all the functionals. We show that the vdW-DF increases the population of non-H-bonded interstitial sites, at distances between the first and second coordination shells. However, it excessively weakens the H-bond network, collapsing the second coordination shell. This structural problem is partially associated to the choice of GGA exchange in the vdW-DF. We show that a different choice for the exchange functional is enough to achieve an overall improvement both in structure and diffusivity.     

Slow dynamics in a primitive tetrahedral network model

We report extensive Monte Carlo and event-driven molecular dynamics simulations of the fluid and liquid phase of a primitive model for silica recently introduced by Ford et al. [J. Chem. Phys. 121, 8415 (2004)]. We evaluate the isodiffusivity lines in the temperature-density plane to provide an indication of the shape of the glass transition line. Except for large densities, arrest is driven by the onset of the tetrahedral bonding pattern and the resulting dynamics is strong in Angell's classification scheme [J. Non-Cryst. Solids 131-133, 13 (1991)]. We compare structural and dynamic properties with corresponding results of two recently studied primitive models of network forming liquids-a primitive model for water and an angular-constraint-free model of four-coordinated particles-to pin down the role of the geometric constraints associated with bonding. Eventually we discuss the similarities between "glass" formation in network forming liquids and "gel" formation in colloidal dispersions of patchy particles.     

Four phases of amorphous water: Simulations versus experiment

Multiplicity of the liquid-liquid phase transitions in supercooled water, first obtained in computer simulations [Brovchenko et al., J. Chem. Phys. 118, 9473 (2003)], has got strong support from the recent experimental observation of the two phase transitions between amorphous ices [Loerting et al., Phys. Rev. Lett. 96, 025702 (2006)]. These experimental results allow assignment of the four amorphous water phases (I-IV) obtained in simulations to the three kinds of amorphous ices. Water phase I (rho approximately 0.90 gcm(3)) corresponds to the low-density amorphous ice, phase III (rho approximately 1.10 gcm(3)) to the high-density amorphous ice, and phase IV (rho approximately 1.20 gcm(3)) to the very-high-density amorphous ice. Phase II of model water with density rho approximately 1.00 gcm(3) corresponds to the normal-density water. Such assignment is confirmed by the comparison of the structural functions of the amorphous phases of model water and real water. In phases I and II the first and second coordination shells are clearly divided. Phase I consists mainly of the four coordinated tetrahedrally ordered water molecules. Phase II is enriched with molecules, which have tetrahedrally ordered four nearest neighbors and up six molecules in the first coordination shell. Majority of the molecules in phase III still have tetrahedrally ordered four nearest neighbors. Transition from phase III to phase IV is characterized by a noticeable drop of tetrahedral order, and phase IV consists mainly of molecules with highly isotropic angular distribution of the nearest neighbors. Relation between the structures of amorphous water phases, crystalline ices, and liquid water is discussed.     

Brownian dynamics simulations of simplified cytochrome c molecules in the presence of a charged surface

Simulations were performed for up to 150 simplified spherical horse heart cytochrome c molecules in the presence of a charged surface, which serves as an approximate model for a lipid membrane. Screened electrostatic and short-ranged attractive as well as repulsive van der Waals forces for interparticle and particle-membrane interactions are utilized in the simulations. At a distance from the membrane, where particle-membrane interactions are negligible, the simulation is coupled to a noninteraction continuum analogous to a heat bath [Geyer et al., J. Chem. Phys. 120, 4573 (2004)]. From the particles' density profiles perpendicular to the planar surface binding isotherms are derived and compared to experimental results [Heimburg et al. (1999)]. Using a negatively charged structureless membrane surface a saturation effect was found for relatively large particle concentrations. Since biological membranes often contain membrane proteins, we also studied the influence of additional charges on our model membrane mimicking bacterial reaction centers. We find that the onset of the saturation occurs for much lower concentrations and is sensitive to the detailed implementation. Therefore we suggest that local distortion of membrane planarity (undulation), or lipid demixing, or the presence of charged integral membrane proteins create preferential binding sites on the membrane. Only then do we observe saturation at physiological concentrations.     

Generalized gradient approximation model exchange holes for range-separated hybrids

We propose a general model for the spherically averaged exchange hole corresponding to a generalized gradient approximation (GGA) exchange functional. Parameters are reported for several common GGAs. Our model is based upon that of Ernzerhof and Perdew [J. Chem. Phys. 109, 3313 (1998)]. It improves upon the former by precisely reproducing the energy of the parent GGA, and by enabling fully analytic evaluation of range-separated hybrid density functionals. Analytic results and preliminary thermochemical tests indicate that our model also improves upon the simple, local-density-based exchange hole model of Iikura et al. [J. Chem. Phys. 115, 3540 (2001)].     

Dynamical properties of the soft sticky dipole-quadrupole-octupole water model: a molecular dynamics study

The dynamical properties of the soft sticky dipole-quadrupole-octupole (SSDQO) water model using SPC/E moments are calculated utilizing molecular dynamics simulations. This new potential for liquid water describes the water-water interactions by a Lennard-Jones term and a sticky potential, which is an approximate moment expansion with point dipole, quadrupole, and octupole moments, and reproduces radial distribution functions of pure liquid water using the moments of SPC/E [Ichiye and Tan, J. Chem. Phys. 124, 134504 (2006)]. The forces and torques of SSDQO water for the dipole-quadrupole, quadrupole-quadrupole, and dipole-octupole interactions are derived here. The simulations are carried out at 298 K in the microcanonical ensemble employing the Ewald method for the long-range dipole-dipole interactions. Here, various dynamical properties associated with translational and rotational motions of SSDQO water using the moments of SPC/E (SSDQO:SPC/E) water are compared with the results from SPC/E and also experiment. The self-diffusion coefficient of SSDQO:SPC/E water is found to be in excellent agreement with both SPC/E and experiment whereas the single particle orientational relaxation time for dipole vector is better than SPC/E water but it is somewhat smaller than experiment. The dielectric constant of SSDQO:SPC/E is essentially identical to SPC/E, and both are slightly lower than experiment. Also, molecular dynamics simulations of the SSDQO water model are found to be about twice as fast as three-site models such as SPC/E.     

Effective electrostatic interactions in solutions of polyelectrolyte stars with rigid rodlike arms

In solutions of star-branched polyelectrolytes, electrostatic interactions between charged arms on neighboring stars can compete with intrastar interactions and rotational entropy to induce anisotropy in the orientational distribution of arms. We explore the influence of arm orientational anisotropy on effective star-star interactions for model stars comprising rigid rodlike arms with evenly spaced charged monomers interacting via an effective screened-Coulomb (Yukawa) potential. Monte Carlo simulation and density-functional theory are used to compute the arm orientational distributions and effective pair potentials between weakly charged stars. For comparison, a torque balance analysis is performed to obtain the configuration and energy of the ground state, in which the torque vanishes on each arm of the two-star system. The degree of anisotropy is found to increase with the strength of electrostatic interactions and proximity of the stars. As two stars begin to overlap, the forward arms are pushed back by interstar arm-arm repulsion, but partially interdigitate due to rotational entropy. At center-center separations approaching complete overlap, the arms relax to an isotropic distribution. For nonoverlapping stars, anisotropy-induced changes in the intra- and interstar arm-arm interactions largely cancel and the effective pair interactions are then well approximated by a simple Yukawa potential, as predicted by linear-response theory for a continuum model of isotropic stars [A. R. Denton, Phys. Rev. E 67, 11804 (2003)]. For overlapping stars, the effective pair interactions in the simple rigid-arm-Yukawa model agree closely with simulations of a molecular model that includes flexible arms and explicit counterions [A. Jusufi et al., Phys. Rev. Lett. 88, 018301 (2002); J. Chem. Phys. 116, 11011 (2002)].