A systematic analysis of the hydration structure of Cs + ions in solution is derived from simulations carried out using a series of molecular models built upon a hierarchy of approximate representations of many-body effects in ion-water interactions. It is found that a pairwise-additive model, commonly used in biomolecular simulations, provides poor agreement with experimental X-ray spectra, indicating an incorrect description of the underlying hydration structure. Although the agreement with experiment improves in simulations with a polarizable model, the predicted hydration structure is found to lack the correct sequence of water shells. Progressive inclusion of explicit many-body effects in the representation of Cs+-water interactions as well as account for nuclear quantum effects is shown to be necessary for quantitatively reproducing the experimental X-ray spectra. Besides emphasizing the importance of many-body effects, these results suggest that molecular models rigorously derived from many-body expansions hold promise for realistic simulations of aqueous solutions.

<div> <div> <div> <p>A systematic analysis of the hydration structure of Cs+ ions in solution is derived from simulations carried out using a series of molecular models built upon a hierarchy of approximate representations of many-body effects in ion-water interactions. It is found that a pairwise-additive model, commonly used in biomolecular simulations, provides poor agreement with experimental X-ray spectra, indicating an incorrect description of the underlying hydration structure. Although the agreement with experiment improves in simulations with a polarizable model, the predicted hydration structure is found to lack the correct sequence of water shells. Progressive inclusion of explicit many- body effects in the representation of Cs⁺-water interactions as well as account for nuclear quantum effects is shown to be necessary for quantitatively reproducing the experimental spectra. Besides emphasizing the importance of many-body effects, these results suggests that molecular models rigorously derived from many-body expansions hold promise for realistic simulations of aqueous solutions. </p> </div> </div> </div>

download file view on ChemRxiv x-wat2_manuscript.pdf (1.73 MiB) download file view on ChemRxiv x-wat2_si.pdf (819.24 KiB) download file view on ChemRxiv toc_figure.pdf (91.20 KiB)

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.

Disorder-driven degradation phenomena, such as structural phase transformations and surface reconstructions, can significantly reduce the lifetime of Li-ion batteries, especially those with nickel-rich layered-oxide cathodes. We develop a general free energy model for layered-oxide ion-intercalation materials as a function of the degree of disorder, which represents the density of defects in the host crystal. The model accounts for defect core energies, long-range dipolar electrostatic forces, and configurational entropy of the solid solution. In the case of nickel-rich oxides, we hypothesize that nickel with a high concentration of defects is driven into the bulk by electrostatic forces as oxidation reactions at the solid-electrolyte interface remove oxygen ions at high potentials (>4.4V vs. Li/Li+). The model is used in battery cycling simulations to describe the extent of cathode degradation when using different voltage cutoffs, in agreement with experimental observations that lower-voltage cycling can substantially reduce cathode degradation. The theory provides a framework to guide the development of cathode compositions, coatings, and electrolytes to enhance rate capability and enhance battery lifetime. The general theory of cation-disorder formation may also find applications in electrochemical water treatment and ion separations, such as lithium extraction from brines, based on competitive ion intercalation in battery materials.