Can DFT-Trained Neural Network Potentials Reproduce Structure, Solvation, and Water-Exchange Properties in Aqueous Magnesium Solutions?

Authors

Sebastian Falkner, Pablo Montero de Hijes, Christoph Dellago, Nadine Schwierz

Published in

Journal of chemical theory and computation. Volume 22. Issue 21. Pages 11432-11441. Nov 10, 2026.

Abstract

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 MgCl2 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 Mg2+, revealing a dissociative exchange mechanism. Transition interface sampling yields exchange rates within 1 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 architecture for describing ion solvation thermodynamics. Our results demonstrate that DFT-trained NNPs can accurately describe Mg2+ hydration structure, diffusion, ion pairing, and exchange kinetics while highlighting limitations of local NNP architectures and the importance of long-range electrostatic contributions to achieve quantitative agreement with experimental ion solvation free energies.

PMID:
43006963
Bibliographic data and abstract were imported from PubMed on 10 Nov 2026.

Read full publication at:
Please sign in to see all details.

Stats

1-terrible, 9-excellent. How would you rate this publication? Sign in in to submit your rating.

Post a comment

You need to be signed in to post comments. You can sign in here.

Comments

There are no comments yet.