Non-equilibrium simulations of thermally induced electric fields in water

Author(s)
P. Wirnsberger, D. Fijan, A. Saric, M. Neumann, C. Dellago, D. Frenkel
Abstract

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.

Organisation(s)
Computational and Soft Matter Physics, Department for Teacher Education
External organisation(s)
University of Cambridge, University of Oxford, University College London
Journal
Journal of Chemical Physics
Volume
144
No. of pages
16
ISSN
0021-9606
DOI
https://doi.org/10.1063/1.4953036
Publication date
06-2016
Peer reviewed
Yes
Austrian Fields of Science 2012
103015 Condensed matter, 103029 Statistical physics
Keywords
ASJC Scopus subject areas
Physics and Astronomy(all), Physical and Theoretical Chemistry
Portal url
https://ucris.univie.ac.at/portal/en/publications/nonequilibrium-simulations-of-thermally-induced-electric-fields-in-water(0284f45e-8d54-4930-8919-8577cb122823).html