skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Effects of polarizability and charge transfer on water dynamics and the underlying activation energies
A large number of force fields have been proposed for describing the behavior of liquid water within classical atomistic simulations, particularly molecular dynamics. In the past two decades, models that incorporate molecular polarizability and even charge transfer have become more prevalent, in attempts to develop more accurate descriptions. These are frequently parameterized to reproduce the measured thermodynamics, phase behavior, and structure of water. On the other hand, the dynamics of water is rarely considered in the construction of these models, despite its importance in their ultimate applications. In this paper, we explore the structure and dynamics of polarizable and charge-transfer water models, with a focus on timescales that directly or indirectly relate to hydrogen bond (H-bond) making and breaking. Moreover, we use the recently developed fluctuation theory for dynamics to determine the temperature dependence of these properties to shed light on the driving forces. This approach provides key insight into the timescale activation energies through a rigorous decomposition into contributions from the different interactions, including polarization and charge transfer. The results show that charge transfer effects have a negligible effect on the activation energies. Furthermore, the same tension between electrostatic and van der Waals interactions that is found in fixed-charge water models also governs the behavior of polarizable models. The models are found to involve significant energy–entropy compensation, pointing to the importance of developing water models that accurately describe the temperature dependence of water structure and dynamics.  more » « less
Award ID(s):
2102656 2117449
PAR ID:
10414566
Author(s) / Creator(s):
;
Date Published:
Journal Name:
The Journal of Chemical Physics
Volume:
158
Issue:
19
ISSN:
0021-9606
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. (, ArXiv.org)
    Given the piecewise approach to modeling intermolecular interactions for force fields, they can be difficult to parameterize since they are fit to data like total energies that only indirectly connect to their separable functional forms. Furthermore, by neglecting certain types of molecular interactions such as charge penetration and charge transfer, most classical force fields must rely on, but do not always demonstrate, how cancellation of errors occurs among the remaining molecular interactions accounted for such as exchange repulsion, electrostatics, and polarization. In this work we present the first generation of the (many-body) MB-UCB force field that explicitly accounts for the decomposed molecular interactions commensurate with a variational energy decomposition analysis, including charge transfer, with force field design choices that reduce the computational expense of the MB-UCB potential while remaining accurate. We optimize parameters using only single water molecule and water cluster data up through pentamers, with no fitting to condensed phase data, and we demonstrate that high accuracy is maintained when the force field is subsequently validated against conformational energies of larger water cluster data sets, radial distribution functions of the liquid phase, and the temperature dependence of thermodynamic and transport water properties. We conclude that MB-UCB is comparable in performance to MB-Pol, but is less expensive and more transferable by eliminating the need to represent short-ranged interactions through large parameter fits to high order polynomials 
    more » « less
  2. ; (, Nanoscale)
    Nature employs channel proteins to selectively pass water across cell membranes, which inspires the search for bio-mimetic analogues. Carbon nanotube porins (CNTPs) are intriguing mimics of water channels, yet ion transport in CNTPs still poses questions. As an alternative to continuum models, here we present a molecular mean-field model that transparently describes ion coupling, yet unlike continuum models, computes ab initio all required thermodynamic quantities for the KCl salt and H + and OH − ions present in water. Starting from water transfer, the model considers the transfer of free ions, along with ion-pair formation as a proxy of non-mean-field ion–ion interactions. High affinity to hydroxide, suggested by experiments, making it a dominant charge carrier in CNTPs, is revealed as an exceptionally favorable transfer of KOH pairs. Nevertheless, free ions, coexisting with less mobile ion-pairs, apparently control ion transport. The model well explains the observed effects of salt concentration and pH on conductivity, transport numbers, anion permeation and its activation energies, and current rectification. The proposed approach is extendable to other sub-nanochannels and helps design novel osmotic materials and devices. 
    more » « less
  3. ; ; ; (, Computational Materials Science)
    The morphological evolution of nanoporous gold is generally believed to be governed by surface diffusion. This work specifically explores the dependence of mass transport by surface diffusion on the curvature of a gold surface. The surface diffusivity is estimated by molecular dynamics simulations for a variety of surfaces of constant mean curvature, eliminating any chemical potential gradients and allowing the possible dependence of the surface diffusivity on mean curvature to be isolated. The apparent surface diffusivity is found to have an activation energy of ~0.74 eV with a weak dependence on curvature, but is consistent with the values reported in the literature. The apparent concentration of mobile surface atoms is found to be highly variable, having an Arrhenius dependence on temperature with an activation energy that also has a weak curvature dependence. These activation energies depend on curvature in such a way that the rate of mass transport by surface diffusion is nearly independent of curvature, but with a higher activation energy of ~1.01 eV. The curvature dependencies of the apparent surface diffusivity and concentration of mobile surface atoms is believed to be related to the expected lifetime of a mobile surface atom, and has the practical consequence that a simulation study that does not account for this finite lifetime could underestimate the activation energy for mass transport via surface diffusion by ~0.27 eV. 
    more » « less
  4. ; (, Physical Chemistry Chemical Physics)
    As known, small HCl–water nanoclusters display a particular dissociation behaviour, whereby at least four water molecules are required for the ionic dissociation of HCl. In this work, we examine how intermolecular interactions promote the ionic dissociation of such nanoclusters. To this end, a set of 45 HCl–water nanoclusters with up to four water molecules is introduced. Energy decomposition analysis based on absolutely localized molecular orbitals (ALMO-EDA) is employed in order to study the importance of frozen interaction, dispersion, polarization, and charge-transfer for the dissociation. The vertical ALMO-EDA scheme is applied to HCl–water clusters along a proton-transfer coordinate varying the amount of spectator water molecules. The corresponding ALMO-EDA results show a clear preference for the dissociated cluster only in the presence of four water molecules. Our analysis of adiabatic ALMO-EDA results reveals a push–pull mechanism for the destabilization of the HCl bond based on the synergy between forward and backward charge-transfer. 
    more » « less
  5. ; ; (, Frontiers in Materials)
    Ionic transport is a critical property for the glass industry, since emerging applications such as sensors, batteries, and electric melting are based on the phenomenon. Short-range interactions (anion-charge carrier) have not been able to explain the total activation barrier observed experimentally, and, as such, it is critical to understand the larger role of all ions in a glass, not just the carrier and the ‘site’ ions. This research focuses on the role of network formers and their impact on diffusion in glasses, something that current models lack an explicit explanation of. Atomistic simulations with randomly generated parameters for the cation potentials and classical simulations were used to determine the diffusion coefficients and activation energies for synthetic network formers. Using this database, explainable machine learning algorithms were employed to explore network former interactions and determine which parameters are the most influential for ion diffusion. Results suggest that the bond length of the cations changes the geometry of the structure contributing the greatest to cation-modifier interactions. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email