Search papers, labs, and topics across Lattice.
This study develops and benchmarks MACE neural network potentials (NNPs) trained on density functional theory data to model aqueous magnesium solutions, addressing the limitations of classical force fields in capturing key solution properties. The NNPs successfully reproduce the hydration structure, diffusion coefficients, and water-exchange mechanisms, revealing a dissociative exchange process that aligns closely with experimental observations. However, the NNPs fall short in accurately predicting solvation free energies, indicating a need for improved architectures to account for long-range electrostatic effects in ion solvation thermodynamics.
DFT-trained neural network potentials can accurately model hydration and exchange kinetics of magnesium ions, but struggle with solvation free energy predictions.
Magnesium ions play an essential role in many biological processes but remain challenging to model in biomolecular simulations. Despite considerable scientific effort, classical force fields fail to simultaneously reproduce key structural, thermodynamic and kinetic solution properties, likely due to their inability to explicitly account for quantum many-body effects. Here, we develop and systematically benchmark MACE neural network potentials (NNPs) for aqueous MgCl$_2$ solutions trained on revPBE-D3/zd and revPBE0-D3/zd density functional theory reference data and assess their ability to reproduce a broad range of experimental solution properties including the structure of the first hydration shell, diffusion coefficient, activity derivative, water-exchange rate and mechanism as well as solvation free energy. Both NNPs accurately reproduce the octahedral structure of the first hydration shell, ion pairing properties and diffusion coefficients. Combining the NNPs with transition path sampling and other enhanced sampling techniques allows us to capture the rare event of water exchange in the first hydration shell of Mg$^{2+}$ revealing a dissociative exchange mechanism. Transition interface sampling yields exchange rates within one order of magnitude of experiment, representing a substantial improvement over classical dissociative force fields. In contrast, the NNP-derived solvation free energy significantly underestimates the experimental value, revealing a limitation of the present local NNP architectures for describing ion solvation thermodynamics. Our results demonstrate that DFT-trained NNPs can accurately describe Mg$^{2+}$ hydration structure, diffusion, ion pairing, and exchange kinetics, while highlighting the need for explicit long-range electrostatic treatments to achieve quantitative agreement with experimental ion solvation free energies.
