Non-equilibrium simulations of thermally induced electric fields in water.

Using non-equilibrium molecular dynamics simulations, it has been recently demonstrated that water molecules align in response to an imposed temperature gradient, resulting in an effective electric field. Here, we investigate how thermally induced fields depend on the underlying treatment of long-ranged interactions. For the short-ranged Wolf method and Ewald summation, we find the peak strength of the field to range between 2 × 10(7) and 5 × 10(7) V/m for a temperature gradient of 5.2 K/Å. Our value for the Wolf method is therefore an order of magnitude lower than the literature value [J. A. Armstrong and F. Bresme, J. Chem. Phys. 139, 014504 (2013); J. Armstrong et al., J. Chem. Phys. 143, 036101 (2015)]. We show that this discrepancy can be traced back to the use of an incorrect kernel in the calculation of the electrostatic field. More seriously, we find that the Wolf method fails to predict correct molecular orientations, resulting in dipole densities with opposite sign to those computed using Ewald summation. By considering two different multipole expansions, we show that, for inhomogeneous polarisations, the quadrupole contribution can be significant and even outweigh the dipole contribution to the field. Finally, we propose a more accurate way of calculating the electrostatic potential and the field. In particular, we show that averaging the microscopic field analytically to obtain the macroscopic Maxwell field reduces the error bars by up to an order of magnitude. As a consequence, the simulation times required to reach a given statistical accuracy decrease by up to two orders of magnitude.

[1]  Chao Zhang,et al.  Computing the dielectric constant of liquid water at constant dielectric displacement , 2015, 1510.00643.

[2]  D. Frenkel,et al.  An enhanced version of the heat exchange algorithm with excellent energy conservation properties. , 2015, The Journal of chemical physics.

[3]  F. Bresme,et al.  Note: How does the treatment of electrostatic interactions influence the magnitude of thermal polarization of water? The SPC/E model. , 2015, The Journal of chemical physics.

[4]  G. Fanourgakis An extension of Wolf's method for the treatment of electrostatic interactions: application to liquid water and aqueous solutions. , 2015, The journal of physical chemistry. B.

[5]  L. G. Macdowell,et al.  Damped reaction field method and the accelerated convergence of the real space Ewald summation. , 2014, The Journal of chemical physics.

[6]  Madan Lamichhane,et al.  Real space electrostatics for multipoles. I. Development of methods. , 2014, The Journal of chemical physics.

[7]  Pengyu Y. Ren,et al.  Classical electrostatics for biomolecular simulations. , 2014, Chemical reviews.

[8]  F. Bresme,et al.  Enhancement of the thermal polarization of water via heat flux and dipole moment dynamic correlations. , 2013, The journal of physical chemistry. B.

[9]  Ikuo Fukuda,et al.  Zero-multipole summation method for efficiently estimating electrostatic interactions in molecular system. , 2013, The Journal of chemical physics.

[10]  F. Bresme,et al.  Water polarization induced by thermal gradients: the extended simple point charge model (SPC/E). , 2013, The Journal of chemical physics.

[11]  D. Bedeaux,et al.  Thermomolecular orientation of nonpolar fluids. , 2012, Physical review letters.

[12]  N. Spaldin,et al.  A beginner's guide to the modern theory of polarization , 2012, 1202.1831.

[13]  F. Bresme,et al.  A comparison of Coulombic interaction methods in non-equilibrium studies of heat transfer in water. , 2011, The Journal of chemical physics.

[14]  F. Bresme,et al.  Water under temperature gradients: polarization effects and microscopic mechanisms of heat transfer. , 2011, Physical chemistry chemical physics : PCCP.

[15]  Tetsu Narumi,et al.  Cutoff radius effect of the isotropic periodic sum and Wolf method in liquid-vapor interfaces of water. , 2011, The Journal of chemical physics.

[16]  Haruki Nakamura,et al.  Molecular dynamics scheme for precise estimation of electrostatic interaction via zero-dipole summation principle. , 2011, The Journal of chemical physics.

[17]  A. Wallqvist,et al.  On the proper calculation of electrostatic interactions in solid-supported bilayer systems. , 2011, The Journal of chemical physics.

[18]  Jocelyn M. Rodgers,et al.  Local molecular field theory for the treatment of electrostatics , 2008, 0809.2845.

[19]  Jocelyn M. Rodgers,et al.  Interplay of local hydrogen-bonding and long-ranged dipolar forces in simulations of confined water , 2008, Proceedings of the National Academy of Sciences.

[20]  Fernando Bresme,et al.  Water polarization under thermal gradients. , 2008, Physical review letters.

[21]  J. Alejandre,et al.  The Wolf method applied to the liquid-vapor interface of water. , 2008, The Journal of chemical physics.

[22]  D. Vanderbilt,et al.  Electric displacement as the fundamental variable in electronic-structure calculations , 2008, 0806.0844.

[23]  Alexander O. Govorov,et al.  Generating heat with metal nanoparticles , 2007 .

[24]  C. Wick,et al.  Simulated surface potentials at the vapor-water interface for the KCl aqueous electrolyte solution. , 2006, The Journal of chemical physics.

[25]  J. D. Gezelter,et al.  Is the Ewald summation still necessary? Pairwise alternatives to the accepted standard for long-range electrostatics. , 2006, The Journal of chemical physics.

[26]  J. Weeks,et al.  Local molecular field theory for effective attractions between like charged objects in systems with strong Coulomb interactions , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[27]  Bernard R Brooks,et al.  Isotropic periodic sum: a method for the calculation of long-range interactions. , 2005, The Journal of chemical physics.

[28]  Stefan M. Kast,et al.  Enhancement of the wolf damped Coulomb potential: Static, dynamic, and dielectric properties of liquid water from molecular simulation , 2002 .

[29]  Berend Smit,et al.  Understanding Molecular Simulation , 2001 .

[30]  G. Hummer,et al.  Treatment of electrostatic interactions in computer simulations and calculation of thermodynamic properties such as free energies and pressures , 1999 .

[31]  R. Fox,et al.  Classical Electrodynamics, 3rd ed. , 1999 .

[32]  M. Berkowitz,et al.  Ewald summation for systems with slab geometry , 1999 .

[33]  S. Phillpot,et al.  Exact method for the simulation of Coulombic systems by spherically truncated, pairwise r−1 summation , 1999 .

[34]  M. Deserno,et al.  How to mesh up Ewald sums. I. A theoretical and numerical comparison of various particle mesh routines , 1998, cond-mat/9807099.

[35]  G. Hummer,et al.  Pressure calculation in polar and charged systems using Ewald summation: Results for the extended simple point charge model of water , 1998, physics/9806038.

[36]  D. Tildesley,et al.  The free surface of water: molecular orientation, surface potential and nonlinear susceptibility , 1997 .

[37]  Eckhard Spohr,et al.  Effect of electrostatic boundary conditions and system size on the interfacial properties of water and aqueous solutions , 1997 .

[38]  B. Brooks,et al.  Effect of Electrostatic Force Truncation on Interfacial and Transport Properties of Water , 1996 .

[39]  M. R. Philpott,et al.  Molecular Dynamics Study of Interfacial Electric Fields , 1996 .

[40]  Bernard R. Brooks,et al.  New spherical‐cutoff methods for long‐range forces in macromolecular simulation , 1994, J. Comput. Chem..

[41]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[42]  Steve Plimpton,et al.  Fast parallel algorithms for short-range molecular dynamics , 1993 .

[43]  O. Steinhauser,et al.  Cutoff size does strongly influence molecular dynamics results on solvated polypeptides. , 1992, Biochemistry.

[44]  K. Suslick,et al.  Interparticle collisions driven by ultrasound. , 1990, Science.

[45]  A. Pohorille,et al.  Comment on "Study on the liquid-vapor interface of water. I. Simulation results of thermodynamic properties and orientational structure". , 1989, The Journal of chemical physics.

[46]  A. Pohorille,et al.  Surface potential of the water liquid-vapor interface. , 1988, The Journal of chemical physics.

[47]  T. Straatsma,et al.  THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS , 1987 .

[48]  Martin Neumann,et al.  Dielectric relaxation in water. Computer simulations with the TIP4P potential , 1986 .

[49]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[50]  S. Nosé A unified formulation of the constant temperature molecular dynamics methods , 1984 .

[51]  Martin Neumann,et al.  Consistent calculation of the static and frequency-dependent dielectric constant in computer simulations , 1984 .

[52]  W. C. Swope,et al.  A computer simulation method for the calculation of equilibrium constants for the formation of physi , 1981 .

[53]  J. Perram,et al.  Simulation of electrostatic systems in periodic boundary conditions. I. Lattice sums and dielectric constants , 1980, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[54]  M. Neumann,et al.  The influence of boundary conditions used in machine simulations on the structure of polar systems , 1980 .

[55]  David van der Spoel,et al.  The Origin of Layer Structure Artifacts in Simulations of Liquid Water. , 2006 .

[56]  Berend Smit,et al.  Accelerating Monte Carlo Sampling , 2002 .

[57]  M. Neumann Computer simulation and the dielectric constant at finite wavelength , 1986 .

[58]  P. P. Ewald Die Berechnung optischer und elektrostatischer Gitterpotentiale , 1921 .