The hydration structure of Na + and K + ions in solution is systematically investigated using a hierarchy of molecular models that progressively include more accurate representations of many-body interactions. We found that a conventional empirical pairwise additive force field that is commonly used in biomolecular simulations is unable to reproduce the extended X-ray absorption fine structure (EXAFS) spectra for both ions. In contrast, progressive inclusion of many-body effects rigorously derived from the many-body expansion of the energy allows the MB-nrg potential energy functions (PEFs) to achieve nearly quantitative agreement with the experimental EXAFS spectra, thus enabling the development of a molecular-level picture of the hydration structure of both Na + and K + in solution. Since the MB-nrg PEFs have already been shown to accurately describe isomeric equilibria and vibrational spectra of small ion−water clusters in the gas phase, the present study demonstrates that the MB-nrg PEFs effectively represent the long-sought-after models able to correctly predict the properties of ionic aqueous systems from the gas to the liquid phase, which has so far remained elusive.

The development of molecular models with quantum-mechanical accuracy for predictive simulations of biomolecular systems has been a long standing goal in the field of computational biophysics and biochemistry. As a first step towards a transferable force field for biomolecules entirely derived from first principles, we introduce a data-driven many-body energy (MB-nrg) potential energy function (PEF) for N-methyl acetamide (NMA), a peptide bond capped by two methyl groups that is commonly used as a proxy for the protein backbone. The MB-nrg PEF is shown to accurately describe the energetics and structural properties of an isolated NMA molecule, including the normal modes of both cis and trans isomers and the energy variation along the isomerization path, as well as the multidimensional potential energy landscape of the NMA-H2O dimer in the gas phase. Importantly, we show that the MB-nrg PEF is fully transferable, enabling molecular dynamics simulations of NMA in solution with quantum-mechanical accuracy. Comparisons with results obtained with a popular pairwise- additive force field for biomolecules and a classical polarizable PEF demonstrate the ability of the MB-nrg PEF to accurately represent many-body effects in NMAH2O interactions at both short and long distances, which is key to guaranteeing full transferability from the gas to the liquid phase.

The development of molecular models with quantum-mechanical accuracy for predictive simulations of biomolecular systems has been a long-standing goal in the field of computational biophysics and biochemistry. As a first step toward a transferable force field for biomolecules entirely derived from "first-principles", we introduce a data-driven many-body energy (MB-nrg) potential energy function (PEF) for N-methylacetamide (NMA), a peptide bond capped by two methyl groups that is commonly used as a proxy for the protein backbone. The MB-nrg PEF is shown to accurately describe the energetics and structural properties of an isolated NMA molecule, including the normal modes of both cis and trans isomers and the energy variation along the isomerization path, as well as the multidimensional potential energy landscape of the NMA−H 2 O dimer in the gas phase. Importantly, we show that the MB-nrg PEF is fully transferable, enabling molecular dynamics simulations of NMA in solution with quantum-mechanical accuracy. Comparisons with results obtained with a popular pairwise-additive force field for biomolecules and a classical polarizable PEF demonstrate the ability of the MB-nrg PEF to accurately represent many-body effects in NMA−H 2 O interactions at both short and long distances, which is key to guaranteeing full transferability from the gas phase to the liquid phase.

The hydration structure of Na+ and K+ ions in solution is systematically investigated using a hierarchy of molecular models that progressively include more accurate representations of many-body interactions. We found that a conventional empirical pairwise additive force field that is commonly used in biomolecular simulations is unable to reproduce the extended X-ray absorption fine structure (EXAFS) spectra for both ions. In contrast, progressive inclusion of many-body effects rigorously derived from the many-body expansion of the energy allows the MB-nrg potential energy functions (PEFs) to achieve nearly quantitative agreement with the experimental EXAFS spectra, thus enabling the development of a molecular-level picture of the hydration structure of both Na+ and K+ in solution. Since the MB-nrg PEFs have already been shown to accurately describe isomeric equilibria and vibrational spectra of small ion–water clusters in the gas phase, the present study demonstrates that the MB-nrg PEFs effectively represent the long-sought-after models able to correctly predict the properties of ionic aqueous systems from the gas to the liquid phase, which has so far remained elusive.

The hydration structure of Na+ and K+ ions in solution is systematically investigated using a hierarchy of molecular models that progressively include more accurate representations of many-body interactions. We found that a conventional empirical pairwise additive force field that is commonly used in biomolecular simulations is unable to reproduce the extended X-ray absorption fine structure (EXAFS) spectra for both ions. In contrast, progressive inclusion of many-body effects rigorously derived from the many-body expansion of the energy allows the MB-nrg potential energy functions (PEFs) to achieve nearly quantitative agreement with the experimental EXAFS spectra, thus enabling the development of a molecular-level picture of the hydration structure of both Na+ and K+ in solution. Since the MB-nrg PEFs have already been shown to accurately describe isomeric equilibria and vibrational spectra of small ion–water clusters in the gas phase, the present study demonstrates that the MB-nrg PEFs effectively represent the long-sought-after models able to correctly predict the properties of ionic aqueous systems from the gas to the liquid phase, which has so far remained elusive.