🚨 New paper in
@NatureComms! 🚨
Do halide ions adsorb at the air/water interface, or do they prefer to remain in bulk water? 🌊
This deceptively simple question has been debated for more than a century, with classical theories, continuum models, molecular simulations, and experiments often producing very different pictures. 🤔
Our new paper is now online:
👉
nature.com/articles/s41467-0…
Vibrational sum-frequency generation (
#vSFG) studies have added to the controversy, leading to opposing molecular interpretations, from extreme enrichment of the heavier halides at the surface to a stratified interface with an ion-depleted topmost layer [J. Phys. Chem. B 108, 2252 (2004); J. Phys. Chem. B 108, 5051 (2004); Nat. Commun. 5, 4083 (2014); Nat. Chem. 16, 644 (2024)]. ⚡️
Because vSFG probes atomic ions only indirectly through the response of interfacial water, translating the spectra into molecular structure depends critically on computational models whose accuracy and predictive power across different environments and properties are not always independently established.
We address this problem using MB-pol for water and a controlled hierarchy of ion–water models built on increasingly complete and physically correct approximations to the many-body expansion. 🖥️
By keeping the description of water unchanged and systematically improving only the ion–water interactions, we reconstruct the major pictures proposed over the past century and identify the physical approximations behind each one:
1️⃣ Purely pairwise-additive (2B)-MB-nrg potentials recover the generic ion exclusion predicted by the image-charge theories of Wagner [Phys. Z. 25, 474 (1924)] and Onsager–Samaras [J. Chem. Phys. 2, 528 (1934)].
2️⃣ The classical polarizable TTM-nrg potentials recover the strong adsorption of the heavier halides predicted by earlier polarizable models, including those of Jungwirth and Tobias [J. Phys. Chem. B 106, 6361 (2002)].
3️⃣ Adding an accurate short-range 2-body ion–water interaction to classical many-body polarization, as in (2B+NB)-MB-nrg, recovers the more moderate qualitative picture of Yan Levin’s extended continuum theory, which combines ionic polarizability with cavitation and interfacial solvation [Phys. Rev. Lett. 102,147803 (2009); Phys. Rev. Lett. 103, 257802 (2009)].
4️⃣ With the full (2B+3B+NB)-MB-nrg
#datadriven #manybody potentials, strong adsorption disappears: F⁻, Cl⁻, and Br⁻ remain bulk-favored, while I⁻ retains only a shallow interfacial minimum of order kT.
The resulting picture supports neither extreme surface enrichment nor a completely ion-depleted topmost layer. Instead, only iodide displays a weak intrinsic surface preference in the dilute single-ion limit.
What we found is that short-range quantum-mechanical many-body effects beyond classical polarization are decisive.
Why is iodide different? Not because its direct interaction with water becomes more favorable at the surface. Instead, a localized reorganization of interfacial water offsets part of the ion–water desolvation penalty. Entropy opposes this energetic stabilization, leaving only a weak surface preference.
At a time when “physics-aware” molecular models are receiving so much attention, this paper also highlights a philosophy that has guided our work for more than a decade. Our MB-pol and MB-nrg
#datadriven #manybody potentials combine physically correct long-range interactions with data-driven representations of short-range quantum mechanics trained on CCSD(T) data.
This paper also completes a major chapter in our work on single-ion hydration. After MB-pol provided a realistic and quantitative description of water from gas-phase clusters to liquid water and ice, MB-nrg extended the same framework to hydrated single ions, connecting ion–water dimers and clusters with structural, thermodynamic, dynamical, and spectroscopic properties in bulk solutions and now at aqueous interfaces. 😎
Huge congratulations to Henry, who led this project, and to Richa (
@RichaRashm), Xuanyu (
@philipzxy), and Saswata (
@quantum_winger) for their contributions! 🎉
The next challenge is to develop a quantitative description of finite-concentration electrolyte solutions, where counterions, ion–ion correlations, ion pairing, and concentration-dependent restructuring take us beyond the Debye–Hückel theory.
We believe our
#datadriven #manybody potentials can get us there. Stay tuned! 🏄♀️
Interested in performing
#datadriven #manybody simulations? Check out
#MBX:
✅
mbxsimulations.com
✅
github.com/paesanilab/MBX
✅
groups.google.com/g/mbx-user…
Big thanks also to
@NSF for supporting this research and to
@ACCESSforCI for providing computational resources! 🙏
@UCSanDiego @UCSDPhySci @UCSDChem @HDSIUCSD @SDSC_UCSD