Developing machine-learned potentials to simultaneously capture the dynamics of excess protons and hydroxide ions in classical and path integral simulations

Austin O. Atsango,Tobias Morawietz,Ondrej Marsalek,Thomas E. Markland
DOI: https://doi.org/10.1063/5.0162066
2023-08-12
Abstract:The transport of excess protons and hydroxide ions in water underlies numerous important chemical and biological processes. Accurately simulating the associated transport mechanisms ideally requires utilizing ab initio molecular dynamics simulations to model the bond breaking and formation involved in proton transfer and path-integral simulations to model the nuclear quantum effects relevant to light hydrogen atoms. These requirements result in a prohibitive computational cost, especially at the time and length scales needed to converge proton transport properties. Here, we present machine-learned potentials (MLPs) that can model both excess protons and hydroxide ions at the generalized gradient approximation and hybrid density functional theory levels of accuracy and use them to perform multiple nanoseconds of both classical and path-integral proton defect simulations at a fraction of the cost of the corresponding ab initio simulations. We show that the MLPs are able to reproduce ab initio trends and converge properties such as the diffusion coefficients of both excess protons and hydroxide ions. We use our multi-nanosecond simulations, which allow us to monitor large numbers of proton transfer events, to analyze the role of hypercoordination in the transport mechanism of the hydroxide ion and provide further evidence for the asymmetry in diffusion between excess protons and hydroxide ions.
Chemical Physics
What problem does this paper attempt to address?
The paper aims to address the following key issues: ### Research Background and Objectives - **Simulation Challenges**: In aqueous solutions, the transport mechanisms of excess protons (H₃O⁺) and hydroxide ions (OH⁻) are crucial for various chemical and biological processes. However, accurately simulating the transport mechanisms of these particles typically requires ab initio molecular dynamics (AIMD) to describe the breaking and forming of chemical bonds, as well as path integral simulations to account for the nuclear quantum effects of hydrogen atoms. These methods, while precise, are computationally very expensive, especially when long time-scale simulations are needed to converge the properties of proton transport. - **Research Objectives**: Develop machine learning potentials (MLPs) to simultaneously capture the dynamic behavior of excess protons and hydroxide ions, and use these MLPs for classical and path integral simulations to significantly reduce computational costs. ### Specific Goals - **Construct MLPs**: Train MLPs using datasets sampled from AIMD simulations to achieve density functional theory (DFT) level accuracy, capable of simultaneously describing the behavior of excess protons and hydroxide ions. - **Validate MLPs**: Verify whether the developed MLPs can reproduce the results of AIMD simulations, including vibrational density of states (VDOS), free energy distribution along the proton-sharing coordinate (δ), and other key properties. - **Simulate and Analyze**: Use the validated MLPs to perform multi-nanosecond classical and path integral simulations to explore the role of hypercoordination in the transport mechanism of hydroxide ions and further confirm the asymmetry in diffusion rates between excess protons and hydroxide ions. ### Major Contributions - **Improved Simulation Efficiency**: The proposed MLPs can effectively simulate the behavior of excess protons and hydroxide ions at a much lower cost than AIMD. - **Exploration of Physical Mechanisms**: Through long time-scale simulations, the specific mechanisms of excess proton and hydroxide ion transport in aqueous solutions are revealed, particularly the impact of hypercoordination and its contribution to the differences in diffusion rates. In summary, this paper overcomes the computational challenges faced by traditional AIMD simulations by developing efficient MLPs, thereby providing an in-depth exploration of the transport mechanisms of excess protons and hydroxide ions in aqueous solutions.