# Non-local X-ray intermolecular radiative decay probes solvation shell of ions in water

Johan Söderström, Lucas M Cornetta, Victor Ekholm, Vincenzo Carravetta, Arnaldo Naves De Brito, Ricardo Marinho, Marcus Agåker, Takashi Tokushima, Conny Såthe, Anirudha Ghosh, et al.

► To cite this version:

Johan Söderström, Lucas M Cornetta, Victor Ekholm, Vincenzo Carravetta, Arnaldo Naves De Brito, et al. .. Non-local X-ray intermolecular radiative decay probes solvation shell of ions in water. Nature Communications, 2025, 16 (1), pp.10046. 10.1038/s41467-025-65581-7. hal-05396918

HAL Id: hal-05396918

https://hal.science/hal-05396918v1

Submitted on 3 Dec 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Non-local X-ray intermolecular radiative decay probes solvation shell of ions in water

Received: 14 February 2025

Accepted: 17 October 2025

Published online: 17 November 2025

<figure><img src="image_3.png" /></figure>

Check for updates

Johan Söderström Ⓔ<fnref n="1" />, Lucas M. Cornetta Ⓔ², Victor Ekholm Ⓔ³, Vincenzo Carravetta⁴, Arnaldo Naves de Brito Ⓔ⁵, Ricardo Marinho Ⓔ⁶,⁷, Marcus Agåker Ⓔ<fnref n="1" />,³, Takashi Tokushima Ⓔ³, Conny Såthe³, Anirudha Ghosh³, Dana Bloß Ⓔ⁸, Andreas Hans Ⓔ⁸, Florian Trinter Ⓔ⁹, Iyas Ismail Ⓔ¹⁰, Debora Vasconcelos<fnref n="1" />, Joel Pinheiro<fnref n="1" />, Yi-Ping Chang Ⓔ¹¹, Manuel Harder¹¹, Zhong Yin Ⓔ¹², Joseph Nordgren<fnref n="1" />, Gunnar Öhrwall Ⓔ³, Hans Ågren<fnref n="1" />,¹³, Jan-Erik Rubensson Ⓔ<fnref n="1" /> & Olle Björneholm<fnref n="1" />✉

Aqueous solutions are crucial in chemistry, biology, environmental science, and technology. The chemistry of solutes is influenced by the surrounding solvation shell of water molecules, which have different chemical properties than bulk water due to their different electronic and geometric structure. It is experimentally challenging to selectively investigate this property-determining electronic and geometric structure. Here, we report experimental results on the non-local X-ray emission process Intermolecular Radiative Decay, for the prototypical ions Na⁺ and Mg²⁺ in water. We show that, in Intermolecular Radiative Decay, an electron from the solvation shell fills a core hole in the solute, and the released energy is emitted as an X-ray photon. We interpret the underlying mechanism using theoretical calculations, and show how Intermolecular Radiative Decay will allow us to meet the challenge of selectively probing the solvation shell from within.

The fundamental importance of water in chemical reactions is intri-
cately linked to the behavior of water molecules in the immediate
vicinity of the reactants<fnref n="1" />⁻³. These solvation shells are interfaces
between the solutes and the bulk solvent. Together with a metal ion,
they form weakly bonded metal-aqua complexes of fundamental
importance in environmental, biological, and applied chemistry. The
structure and dynamics in solvation shells differ substantially from
bulk water; for example, the water molecules around metal cations are
stronger proton donors than the bulk water molecules³,⁴. In the case of
ionic solutes, the water molecules surrounding an ion are oriented
with their positive (negative) end towards the anion (cation), and with

increasing ionic charge, the ion-water distance decreases, leading to
the ion-water interaction changing from ion-dipole towards dative
covalent bonding. The importance of the solvation shell extends
beyond small inorganic ions; for example, the functionality of large
biomolecules such as proteins is affected by their surrounding
hydration shell⁵,⁶.

Several methods are used to determine the molecular geometry
of a solvation shell. X-ray and neutron scattering are general methods
to measure interatomic distances. They are straightforward to apply to
liquids, and they can also provide information about solvation shells⁵,⁷.
The interatomic distance and coordination around specific solutes can

<footnote n="1">¹Department of Physics and Astronomy, Uppsala University, Uppsala, Sweden. ²Instituto de Física, Universidade de São Paulo, São Paulo, Brazil. ³MAX IV Laboratory, Lund University, Lund, Sweden. ⁴Institute of Chemical and Physical Processes, CNR-IPCF, Pisa, Italy. ⁵Institute of Physics Gleb Wataghin, State University of Campinas, Campinas, Brazil. ⁶Department of Physics, University of Brasília, Brasília, Brazil. ⁷Institute of Physics, Federal University of Bahia, Salvador, Brazil. ⁸Institute of Physics, University of Kassel, Kassel, Germany. ⁹Molecular Physics, Fritz-Haber-Institut der Max-Planck-Gesellschaft, Berlin, Germany. ¹⁰Sorbonne Université, CNRS, Laboratoire de Chimie Physique - Matière et Rayonnement, LCPMR, Paris, France. ¹¹European XFEL, Schenefeld, Germany. ¹²International Center for Synchrotron Radiation Innovation Smart, Tohoku University, Sendai, Japan. ¹³Faculty of Chemistry, Wroclaw University of Science and Technology, Wroclaw, Poland. ✉e-mail: Johan.Soderstrom@physics.uu.se; Olle.Bjorneholm@physics.uu.se</footnote>
be inferred from measurements of extended X-ray absorption fine structure (EXAFS)⁷,⁸. Nuclear magnetic resonance (NMR) and vibrational spectroscopy selectively probe the structure and dynamics around solutes²,⁵,⁹. While these techniques successfully provide information about the geometry, the investigation of the electronic structure of solvation shells remains an experimental challenge. An experimental probe of the electronic interaction between the solute with the neighboring solvent could thus lead to an increased understanding of this interaction.

X-ray-based spectroscopies are widely used to probe the electronic structure of matter¹⁰. Still, progress has been hampered by limited selectivity for solvation shells. The element specificity associated with quasi-atomic core holes is insufficient to separate water molecules in the solvation shell from bulk water. In the soft X-ray range, core-hole decay is typically dominated by non-radiative local Auger-Meitner decay, where one valence electron fills the core hole while another is emitted. Non-local decay processes, such as intermolecular Coulombic decay (ICD), are sensitive to local neighbors: the decay of a core hole on a solute involves electrons from neighboring water molecules¹¹. In fact, such core ICD processes have been observed for ions, such as Mg²⁺, Al³⁺, K⁺, and Ca²⁺ in aqueous solutions¹²⁻¹⁵. Very recently, resonant ICD was used to probe the local surroundings of Ca²⁺ in water, including ion pair formation¹⁶. However, due to the short inelastic mean free path of electrons in water, these ICD electrons can only provide information about the surface region of the solution. In radiative core-hole decay, X-rays are emitted as outer electrons fill atom-specific local core holes, which allows for more bulk-sensitive studies than using ICD electrons. The radiative decay of solvated ions is largely unexplored. Probing non-local radiative processes of solvated ions to gain information about nearest neighbors, in analogy to non-radiative ICD, has not been reported before for any type of liquid system.

The electronic configuration of the free metal ions Na⁺ and Mg²⁺ is neon-like, 1s²2s²2p⁶, with the highest occupied states being the 2p core levels; for simplicity, we label these metal ions M of charge q as Mq. Photoionization of the M 1s level in such an Mq ion results in an Mq+1 ion with a 1s core hole. The most energetic radiative decay of the 1s core hole in these Mq+1 ions is the Kα decay, Mq+1 1s⁻¹ → Mq+1 2p⁻¹ + hv. In aqueous solution, Na⁺ and Mg²⁺ ions are surrounded by a solvation shell, which on average consists of six water molecules in an octahedral arrangement at an average ion-oxygen distance of 2.43 and 2.07 Å, respectively¹⁷. For the solvated Na⁺ and Mg²⁺ ions, any radiative decay at higher photon energy than Kα would thus indicate a process involving electrons from the surrounding water molecules. Here, we demonstrate, by measurements and theoretical calculations, how an electron from a neighboring water molecule fills the M 1s core hole, and the released energy is emitted as an X-ray photon, as schematically

shown in Fig. 1. The non-local character of the intermolecular radiative decay (IRD) enables unique investigations of the solvation shell.

Radiative transitions involving a core hole on one atom and a valence orbital mainly located on a neighboring atom have been well known for a long time. In transition metal-oxygen ligand systems, such as MnO₄⁻, CrO₄²⁻, and VO₄³⁻, interatomic radiative decay has been observed from orbitals primarily located on the ligand oxygen atoms into core holes centered on the metal ion, see, e.g., Refs. 18–23. Compared to the ions studied here, which form much weaker bonds with the surrounding water, the systems studied in Refs. 18–23 are held together more strongly either by covalent bonds in molecules, or ionic/covalent bonds in the crystalline systems. These interatomic transitions observed for the strongly bonded molecular systems are still within the same molecule and can therefore be characterized as intramolecular. IRD, on the other hand, is intermolecular, opening up the possibility of probing chemically and biologically relevant systems in an aqueous environment.

In this work, we use the prototypical and well-known Na⁺ and Mg²⁺ ions in water as showcases to demonstrate the existence of IRD. Moreover, by theoretical modeling of the process, we show how the IRD depends on factors such as the electronic structure, ion-water distances, orientational (dis-)order, and chemical composition of the immediate surroundings. This demonstrates the potential of IRD to probe the solvation shell from within, providing uniquely selective information on how it differs from bulk water. We foresee that IRD can find applications in studies of the interaction between weakly interacting neighbors, much in the same way as ICD in the non-radiative channel.

## Results

### IRD: observation and interpretation

Ionization of the 1s level of either of the Mq metal ions results in an Mq+1 1s⁻¹ state, and in the middle of Fig. 2, we present the X-ray emission spectra from the decay of this core-hole state for solvated Na⁺ and Mg²⁺ ions. The main spectral peak for both ions is located toward the lower-photon-energy side of the spectra, ~1041 eV for Na and ~1254 eV for Mg. These peaks are due to the Kα decay, Mq+1 1s⁻¹ → Mq+1 2p⁻¹ + hv, i.e., a local process, in which the M 1s hole is filled by a M 2p electron.

At somewhat higher photon energy than the Kα, around ~1060 eV for Na and ~1295 eV for Mg, there are weaker spectral features, which moreover exhibit substructures. As M 2p is the highest occupied orbital of free Na⁺ and Mg²⁺ ions, these weaker spectral features must involve the neighboring water molecules. This can schematically be written as Mq+1 1s⁻¹ + H₂O → Mq + H₂O⁺val⁻¹ + hv, where val⁻¹ denotes a valence hole on the water molecules in the first solvation shell. Just as the Kα photon energy is determined by the energy difference between the M 1s⁻¹ and M 2p⁻¹ states, the photon energy of such non-local

<figure><img src="image_4.png" /><figcaption>**Fig. 1** | Schematic illustration of the Intermolecular Radiative Decay (IRD) process for a solvated Mq ion. From the left, we show the ground state with a Mq ion surrounded by six water molecules, followed by core photoionization, which results in a core-ionized Mq+1 ion. In the IRD process, an electron from a solvation-</figcaption></figure>

shell water molecule fills the M core hole, and the released energy is emitted as an X-ray photon. In the final state, the Mq ion is back to its ground state, and the hole is in the valence levels of the solvation-shell water molecules.
<figure><img src="image_5.png" /><figcaption>**Fig. 2** | Theoretical and experimental overview of the IRD signal compared to Kα. Top: Schematic illustration of the electronic transitions in the local Kα decay (green box) and non-local IRD decay (red box). Middle: Experimentally measured X-ray emission spectra for Na⁺ and Mg²⁺ after Ionization of the metal ion M in the energy range of Kα and IRD. For both ions, the spectra exhibit two main features, of which the most intense one is due to the local Kα emission, M⁹⁺¹Is⁻¹ → M⁹⁺¹2p⁻¹ + hv. The weaker features at higher photon energies agree well with the energies estimated for the non-local intermolecular radiative decay IRD, M⁹⁺¹Is⁻¹ + H₂O → M⁹⁺ + H₂O⁺val⁻¹ + hv. Bottom: Theoretical X-ray emission spectra for Na⁺ and Mg²⁺ calculated for M⁹⁺[H₂O]₆ clusters. The absolute energy scales of the theoretical spectra are shifted so that the Kα emission aligns with the experimental value. The feature around 1277 eV in the Mg spectrum is due to IRD involving the inner-valence orbital 2a₁ of water, which mainly consists of O 2s. This will not be further discussed, as the corresponding spectral region was not measured. The theoretical curves were obtained using the Hartree-Fock approximation.</figcaption></figure>

radiative decay processes would similarly be given by the energy difference between the M⁹⁺¹Is⁻¹ + H₂O and M⁹⁺ + H₂O⁺val⁻¹ states.

In a first approximation, we can use experimentally determined values for M⁹⁺¹Is⁻¹ from Ref. 14 and H₂O⁺val⁻¹ to roughly estimate the photon energy. For bulk liquid water, the H₂O⁺val⁻¹ states span from -10 to -20 eV binding energy. Assuming that the H₂O⁺val⁻¹ states of the solvation shell are not very different, we use 15 eV as an approximate average value for the H₂O⁺val⁻¹ energy. This yields about 1062 and 1295 eV for the photon energy for such non-local radiative decays for solvated Na⁺ and Mg²⁺ ions, respectively. These values agree quite well with the photon energies of the spectral features observed at higher photon energies, thereby supporting the assignment.

For an accurate interpretation of the observed features, it is relevant to estimate their intensity as predicted by a theoretical model. We have then considered a number of molecular structures and computational models to predict the X-ray emission spectra of solvated Na⁺ and Mg²⁺ ions; details will be provided in the Theory section and in the Supplementary Information. In the lower panel of Fig. 2, we show the spectra computed for a model consisting of Na⁺ and Mg²⁺ ions within a cluster of six water molecules and using the Hartree-Fock (HF) approximation. In addition to the local Kα decay at lower photon energies, spectral features at higher photon energies are obtained for both ions; they coincide very well with the experimentally observed bands. Moreover, we measured the relative intensity of the IRD line vs. the Kα line of ~1% for Na and ~1.5% for Mg. The calculated IRD/Kα

intensity ratios are 1.1 and 1.5%, respectively, i.e., in excellent agreement with the experimental result.

Additionally, an analysis of the involved electronic states, to be further discussed below, reveals that the contributions to the IRD bands are primarily due to decay from hybrid orbitals composed mainly of water and, to a lesser extent, metal ion orbitals.

The IRD process can be seen as a radiative relative of the growing family of non-local decay processes, such as ICD and electron-transfer-mediated decay (ETMD)¹¹,²⁴. IRD bears some resemblance to the non-local valence-hole decay process, radiative charge transfer (RCT). This process, in which a photon is emitted when a valence electron is transferred from the neutral environment to a charged species with one or more valence holes, is well-known from ion collision studies, see, e.g., Refs. 25,26. Closer to the present context, RCT has been indirectly inferred to occur for Ne and Ar dimers²⁷⁻³⁰, and recently directly observed for homogeneous and heterogeneous rare-gas clusters³¹⁻³³. While RCT has been observed in cases where other decay channels are energetically forbidden, the IRD observed here occurs in competition with non-radiative Auger-Meitner decay and local radiative Kα decay. Unlike the valence-hole decay process, RCT, and IRD involve a core hole, which opens the possibility to use them for chemically selective studies. In order to observe IRD a non-zero overlap between the solute and solvent orbitals is required. Furthermore, in the present systems IRD and Kα are well separated; in other solutes, different energetic overlap may complicate interpretation.

We conclude that the spectral features observed at higher photon energies are due to a non-local radiative decay process, IRD, involving the transfer of an electron from the water molecules in the first solvation shell to the core-ionized metal ion. The calculated spectra agree quite well with the experiments in terms of energy positions of the IRD features and their intensity relative to Kα. This gives us confidence to use the calculations in combination with the experimental results to further explore the mechanisms of the IRD decay, how it depends on the surroundings of the ion, and to discuss how IRD can provide an insider's view of the first solvation shell.

## The IRD mechanism

To investigate the mechanisms behind IRD, specifically what influences its intensity and where the electron filling the core hole comes from, we will start by examining a simplified type of model system M⁹⁺[H₂O]₆. This consists of a metal cation M⁹⁺ surrounded by six water molecules in a symmetric D₂h structure, see Supplementary Fig. 3. In the one-particle picture, the electronic structures of these systems can be described by orbitals, here denoted hybrid orbitals ψᵢμ, each of which is built by linear combinations of atomic orbitals (LCAO) χᵢ. In the usual manner, the atomic orbitals are a set of Cartesian-Gaussian functions centered on the metal ion and on all water-related atoms. For simplicity, the χ has been indexed simply by v, in such a way that one can write ψᵢμ = ∑<sub>ν</sub> Cᵢμ<sup>v</sup>χᵣ, where Cᵢμ<sup>v</sup> is the hybrid orbital coefficient. The most relevant hybrid orbitals ψᵢμ of the symmetric model are shown in Supplementary Fig. 4.

As a starting point, we will follow the development of the charge distribution from the initial state, over the core-ionized state, to the different final states reached by radiative decay. Figure 3 shows the Mulliken charges of the M⁹⁺[H₂O]₆ ions in the three states, see also Supplementary Table 1. For a free, isolated Na/Mg ion, the initial state has a charge of +1/+2, and both the core-ionized and the final states have a charge of +2/+3, as indicated in black. For the solvated ions, the total charge develops in the same way, but the presence of the water molecules allows for charge redistribution. In the ground state of the Na⁺[H₂O]₆ ion cluster, the Na ion charge is +0.95, i.e., very close to +1. A small charge of +0.05 is delocalized onto the six water molecules, which means that the Na⁺-water interaction is quite close to an ideal ion-dipole interaction. For the Mg²⁺[H₂O]₆ ion cluster, the charge on the Mg ion is +1.35, with +0.65 on the six water molecules. The Mg²⁺.
water interaction does thus contain charge transfer via hybridization leading to some dative covalent bonding in addition to the ion-dipole interaction.

Core ionization increases the total charge to +2/+3 for Na/Mg. In the core-ionized Na²⁺[H₂O]₆ ion cluster, the Na ion has charge +1.48 instead of +2, implying +0.52 on the waters. In contrast to the ground state, there is now a charge transfer from the water to the ion, amounting to 0.52 of an electron. For the core-ionized Mg³⁺[H₂O]₆ ion, the +3 total charge is even more distributed, with +1.75 on Mg and +1.25 on the water.

The radiative core-hole decay does not change the total charge, but the charge distribution is entirely different after local Kα decay and IRD, as shown in Fig. 3. For Na and Mg, Kα decay increases the charge of the metal ion and decreases the charge on water, compared to the intermediate core-hole state with equal total charge. In contrast, IRD instead decreases the charge of the metal ion and increases the charge on water. The final state after Kα decay has the majority of the charge on the metal ion, but after IRD, the majority of the charge is instead on the water. This clearly shows that the net effect of the IRD process is

that the charge from the water fills the core hole on the metal ion. In the present context, the charge transfer upon core ionization and core-hole decay can be regarded as instantaneous, but the ultrafast dynamics of this may be possible to probe in the not-too-distant future, using, e.g., attosecond X-ray pulses from free-electron lasers.

We can take the discussion of the IRD process further by analyzing what mechanism provides the intensity of the IRD transitions. The spectral intensity of transitions from the different hybrid orbitals ψμ, for a given μ, to the M 1s hole depends on the set of orbital coefficients Cμμ and the dipole transition moment operator t projected on the atomic orbital basis set |χ⟩. The intensity is then proportional to |Tμ|² = |∑ν tM1S,ν Cμμ|².

The IRD bands are due to decay from three different hybrid orbitals ψμ, which consist mainly of 1b₂, 3a₁, and 1b₁ water outer-valence molecular orbitals, respectively. Of the three IRD lines, the one from 3a₁ has >10 times higher intensity than the ones from 1b₂ and 1b₁ in the idealized geometries studied here, see Supplementary Table 2, and we will consequently focus our discussion on this dominant 3a₁ line.

In Table 1, we present the main atomic orbital Cμμ contributions to the hybrid orbitals ψμ, the symmetry-allowed Cartesian component of tM1S,ν—simply denoted by tν hereafter—as well as their product tν Cμμ. Note that in the present structural model, only ungerade hybrid orbitals contribute to the spectra for symmetry reasons. Similarly, the M ns contributions to the hybrid orbitals also do not contribute as tM1S,Mns = 0 due to dipole selection rules. Table 1 forms a basis to discuss the mechanisms behind the different radiative transitions filling the M 1s core hole. Different hybrid orbitals ψμ have different compositions of O 2p and M np contributions, described by the Cμμ orbital coefficients.

The local Kα emission intensity is due to decay from a hybrid orbital for which the by far largest orbital coefficient is for M 2p, making this an essentially atomic orbital centered on the metal ion. For both Na and Mg, the corresponding transition moment tM2p is 5–10 times higher than any of the other transition moments, resulting in a very high intensity of the Kα line.

In contrast to the Kα case, the 3a₁ IRD line is due to decay from a hybrid orbital which is strongly dominated by O 2p, with much smaller M np contributions. The final 3a₁⁻¹ state after the IRD thus primarily consists of a hole delocalized on the water molecules around the metal ion.

As discussed above, the spectral intensity of decay from ψμ is given by |Tμ|² = |∑ν tM1S,ν Cμμ|². As seen in Table 1, for 3a₁ of both Na and Mg, the by far largest tM1S,ν Cμμ term in the sum is due to M np contributions, which are roughly two orders of magnitude larger than

<figure><img src="image_6.png" /><figcaption>**Fig. 3 |** Mulliken charges in the ground, intermediate, and final states for Na⁺[H₂O]₆ (lower panel) and Mg²⁺[H₂O]₆ (upper panel) for the radiative decay of the core hole. The total charge (black) is decomposed into charge on the metal ion (red), and on the water molecules (blue). For both metal ion and water, the final-state charge after both local Kα decay and non-local IRD is shown. The connecting lines are only guides for the eye.</figcaption></figure>

<table><caption>**Table 1 |** The main contributions to the spectral intensity for Kα and 3a₁ IRD for the symmetric Mq[H₂O]₆ model, with Mq = Na⁺ and Mg²⁺, in atomic units</caption><thead><tr><th rowspan="3">Metal atom<br>Transition</th><th colspan="4">Na</th><th colspan="4">Mg</th></tr><tr><th rowspan="2">t<sub>v</sub></th><th colspan="2">K<sub>α</sub></th><th colspan="2">IRD 3a<sub>1</sub></th><th rowspan="2">t<sub>v</sub></th><th colspan="2">K<sub>α</sub></th><th colspan="2">IRD 3a<sub>1</sub></th></tr><tr><th>C<sup>μμ</sup><sub>M2p</sub></th><th>t<sub>v</sub>C<sup>μμ</sup><sub>M2p</sub></th><th>C<sup>μμ</sup><sub>3a<sub>1</sub></sub></th><th>t<sub>v</sub>C<sup>μμ</sup><sub>3a<sub>1</sub></sub></th><th>C<sup>μμ</sup><sub>M2p</sub></th><th>t<sub>v</sub>C<sup>μμ</sup><sub>M2p</sub></th><th>C<sup>μμ</sup><sub>3a<sub>1</sub></sub></th><th>t<sub>v</sub>C<sup>μμ</sup><sub>3a<sub>1</sub></sub></th></tr></thead><tbody><tr><td>Atomic orbital v</td><td>× 10<sup>-3</sup></td><td></td><td></td><td></td><td></td><td>× 10<sup>-3</sup></td><td></td><td></td><td></td><td></td></tr><tr><td>M 2p</td><td>49.11</td><td>0.97</td><td>47.6</td><td>0.08</td><td>3.93</td><td>47.05</td><td>0.99</td><td>46.6</td><td>0.03</td><td>1.41</td></tr><tr><td>M 3p</td><td>-4.53</td><td>0.01</td><td>-0.045</td><td>-0.08</td><td>0.362</td><td>-6.78</td><td>-0.01</td><td>0.068</td><td>-0.35</td><td>2.37</td></tr><tr><td>M 4p</td><td>0.45</td><td>-0.02</td><td>-0.009</td><td>0.11</td><td>0.049</td><td>0.93</td><td>0.01</td><td>0.009</td><td>0.18</td><td>0.167</td></tr><tr><td>O 2p</td><td>0.04</td><td>0.00</td><td>0.002</td><td>0.81</td><td>0.032</td><td>0.06</td><td>0.00</td><td>0.047</td><td>0.79</td><td>0.047</td></tr><tr><td>|T<sub>μ</sub>|<sup>2</sup> × 10<sup>-6</sup></td><td colspan="4">2405.99</td><td colspan="4">2216.23</td><td colspan="2">14.74</td></tr></tbody></table>

The considered value of the M-O distance R was 2.3 Å for Na and 2.1 Å for Mg. In the columns, tν denotes the respective symmetry-allowed Cartesian component of the atomic dipole transition moments tM1S,ν (a. u.) and Cμμν denotes the orbital coefficient of the main atomic orbitals v in the LCAO description of the hybrid orbitals ψμ. In the bottom row, |Tμ|² = |∑ν tμ Cμμν|² is proportional to the total spectral intensity of the Kα and 3a₁ IRD transitions. Note that |Tμ|² also includes minor contributions from M np orbitals with n larger than 4. Data are from molecular orbitals optimized for the ground state by HF.
the ones related to O nl. The different M np contributions differ in the sense that M 2p is fully occupied already in the free ions, while M 3p and M 4p are partially occupied in the solvated ion due to hybridiza- tion with the water molecular orbitals. For Na, the t<sub>M1s,ν</sub> C<sub>μ</sub>'' term associated with Na 2p is ten times higher than any other, whereas for Mg the t<sub>M1s,ν</sub> C<sub>μ</sub>'' term associated with Mg 3p is somewhat higher than for Mg 2p. This is consistent with the above-discussed larger ground- state charge transfer from water to the metal ion for Mg<sup>2+</sup> than for Na<sup>+</sup>, resulting from the increased ion-water hybridization with decreasing ion-water distance. For Na, the intensity thus mainly comes from Na 2p-O 2p hybridization, while for Mg also ground-state charge transfer into Mg 3p is important. The intensity in IRD is therefore mainly due to transitions from the metal atomic orbital components of the hybrid orbitals. The small C<sub>μ</sub>'' metal orbital coefficients of the ψ<sub>μ</sub> hybrid orbitals are compensated by the atomic transition moments t<sub>ν</sub> being much larger for the metal orbitals than for the water orbitals.

Based on the theory presented above, our interpretation is that the electron filling the core hole originates from the water molecules in the first solvation shell, and the IRD signal strength is attributed to hybridization with the metal-ion orbitals. The hybrid orbitals have relatively small metal character, implying that the electron is to a great extent localized on the water molecules. However, the associated high dipole transition moment of the metal component of the hybrid orbital results in a much higher contribution to the intensity than from the water component, which instead provides a very low dipole transition moment. At the same time, the final-state hole is mainly delocalized on the nearest water molecules, meaning that the electron filling the M<sup>q+1</sup> 1s core hole effectively mainly comes from the solvation shell water, and its energy reflects such delocalization. In spite of this non-local character of IRD, the strong contribution to the intensity from the small metal character of the hybrid orbitals makes the one-center approximation a good approximation, as elaborated further in Supplementary Note 2.3. Moreover, by varying the M–O distance, as discussed in Supplementary Note 2.5, the IRD efficiency is seen to decrease rapidly with increasing distance, due to a decrease of the metal-water orbital mixing. This means that IRD, to a good approximation, only involves the first solvation shell.

The above discussion is based on calculations using an idealized,
highly symmetric structure, which is useful to capture the essentials of
the IRD mechanism. For a similar analysis using a more realistic
structure model suggested by molecular dynamics (MD), see Supple-
mentary Note 2.4.

IRD and the electronic structure

As we will show, IRD opens new possibilities for chemically selective
probing the electronic structure of the first solvation shell around a
solvated ion. Starting with bulk liquid water, the electronic structure
has been studied by valence-band photoelectron spectroscopy (PES),
providing the H₂O⁺ val⁻¹ energies, see, e.g., Ref. 34. The valence-band
photoelectron spectrum of liquid water, shown in Fig. 4, consists of
three features. These are associated to ionization of the three outer-
most valence orbitals of the water molecule, Ib₁, 3a₁, and Ib₂, modified
by overlap and screening in the liquid environment³⁵. As a simplistic
first approximation, the solvation-shell H₂O⁺ solv val⁻¹ binding energies
would be expected to increase around a positively charged ion, but we
should also consider that the geometrical arrangement of the close to
six water molecules surrounding a cation is quite different from that of
the tetrahedral-like structure of the close to four water molecules
surrounding a typical water molecule in bulk liquid water. The water-
water orbital interactions can thus be expected to be different in the
solvation shell, and in addition there may be ion-water orbital inter-
actions, including polarization of the orbitals by the electric field of the
ion. We can obtain the corresponding one-hole energies for the
solvation-shell water, E(H₂O⁺ solv val⁻¹), from the IRD spectra as the
difference between the M 1s binding energy of the metal ion, E(M 1s⁻¹),

<figure><img src="image_7.png" /><figcaption>**Fig. 4 | Comparison of one-hole spectra of H₂O, Na⁺, and Mg²⁺.** From the top, liquid-water PES spectrum recorded using a photon energy of 532 eV, i.e., under relatively bulk-sensitive conditions³¹ (blue), Na IRD spectrum (red dots, experiment, and red dashed, theory), Mg IRD spectrum (green dots, experiment, and green dashed, theory). Experimental and theoretical IRD spectra are enlargements of the data shown in Fig. 2. Solid black curves indicate a fitting of the experimental data using a minimal number of Gaussian curves. Original theoretical data for Na and Mg are represented by black bars, and the dashed curves for their convolution with a Gaussian function (FWHM = 0.8 eV) to facilitate the comparison with the experimental spectra. Vertical solid lines indicate calculated binding energies from Ref. 14.</figcaption></figure>

and the photon energy of the IRD transitions, hv(IRD); E(H₂O⁺solv val⁻¹) = E(M 1s⁻¹) – hv(IRD).

Turning to the Na IRD spectrum, shown in Fig. 4, it is seen to
exhibit two peaks with energies very close to the two most intense
peaks of bulk water. As for bulk water, we provisionally interpret these
as corresponding to holes in the 3a₁ and 1b₁ orbitals in the first solva-
tion shell. For Mg, there are three peaks, of which the two stronger
ones exhibit clear shifts to higher one-hole energies relative to
bulk water.

Figure 4 also shows calculated IRD spectra, obtained by HF calculations for Na⁺[H₂O]₆ and Mg²⁺[H₂O]₆ in a representative cluster geometry derived from MD simulations, see structure in the middle of Supplementary Fig. 3. For both ions, the calculations and the experiments agree quite well in terms of energies. For the spectral shape, there is also a good agreement for Mg with 3a₁ dominating over 1b₁ and 1b₂. For Na, the calculated spectrum also peaks in the 3a₁ range, whereas for the experimental spectrum 1b₁ is even stronger than 3a₁. We tentatively ascribe this discrepancy as due to orientational disorder of the water molecules in the solvation shell, as discussed in Supplementary Note 2.6.

It should be noted, however, that the calculations show that the experimentally measured peaks cannot uniquely be assigned to holes in the 1b₁, 3a₁, and 1b₂ molecular orbitals of individual solvation-shell molecules. Instead, the final one-hole states are delocalized over two or more solvation-shell molecules. The orbital composition and energy
of the final one-hole states also vary with the different relative orientations of the water molecules. This means that spectral contributions from holes in the three valence molecular orbitals of water overlap. That said, the relative contribution of 1b₁ holes is highest for the low-energy states, whereas the relative contribution of 3a₁ holes is highest for states of higher energy. We will, therefore, with the above reservations in mind, continue to denote the IRD spectral features as 1b₁, 3a₁, and 1b₂.

We can also compare the observed $H_2O_{solv}^+ val^{-1}$ energies for Na and Mg with calculated energies from Ref. 14, shown as vertical lines in Fig. 4. Contrary to the simplistic expectation, the calculations predict a small shift towards lower one-hole energies relative to pure water for Na. This was explained as due to the structure-breaking properties of Na, causing a disruption of the water hydrogen-bond network.

For Mg, our results and the calculations from Ref. 14 agree upon a shift towards higher energy relative to pure water. For both ions, the calculated energies from Ref. 14 are somewhat lower than our observed ones. This may at least partly be due to the experimentally obtained $H_2O_{solv}^+ val^{-1}$ states being vibrationally excited, whereas the calculated energies are for the vibrational ground states.

The present results based on IRD can also be discussed relative to the O 1s X-ray Emission Spectroscopy (XES) spectra of solvation-shell water obtained by subtracting a suitably scaled spectrum for pure water from a spectrum of a $MgCl_2$ solution³⁶. The two approaches probe the same electronic final states, $H_2O_{solv}^+ val^{-1}$, but in different ways. The spectra obtained by the subtractive procedure in Ref. 36 contain contributions from solvation shells around both $Mg^{2+}$ and $Cl^-$ ions, whereas IRD probes the solvation shells around the different ions selectively. The spectral shapes are also quite different, as the processes are very different, $M^{q+1}1s^{-1} + H_2O_{solv}^+ \rightarrow M^q + H_2O_{solv}^+ val^{-1} + hv$ for the present IRD spectra, versus $H_2O_{solv}^+O 1s^{-1} \rightarrow H_2O_{solv}^+ val^{-1} + hv$ for the differential O 1s XES spectra of Ref. 36. This results in different relative intensities of $H_2O_{solv}^+ val^{-1}$ electronic states, as well as different populations of vibrational states. The spectra obtained by the subtractive procedure are also influenced by the nuclear dynamics during the O 1s core-hole lifetime, a complication not affecting the IRD spectra.

IRD appears to be a general phenomena for ions in liquid surrounding, in addition to these ions with neon-like electronic configuration, we have also observed IRD for transition metal ions ($Cu^{2+}$) and anions ($F^-$) in aqueous surrounding, see Supplementary Note 1.2.

We conclude that IRD can be used to selectively probe the $H_2O_{solv}^+ val^{-1}$ energies of the first solvation shell around ions, revealing how the electronic structure differs from bulk water. This would also include the presence of other species, e.g., contact ion pairs, see Supplementary Note 2.7.

## Discussion

Combining experiments and calculations, we have demonstrated the existence of IRD for $Na^+$ and $Mg^{2+}$ ions in aqueous solution. In IRD, a 1s core hole on Na or Mg is filled by an electron from a neighboring species, and an X-ray photon is emitted. IRD is thus a non-local decay process, which in some sense can be considered a radiative analog of ICD.

The transitions observed in the IRD spectra are interpreted as getting intensity from hybridization between valence orbitals on water molecules and the occupied orbitals on the metal cation. The associated high atomic transition moments result in a much higher contribution to the intensity than from the dominant oxygen character. At the same time, the final-state hole is delocalized on the nearest water molecules, meaning that the electron filling the metal 1s core hole effectively comes from water. For the observed IRD lines, the intensity comes from the relatively small electron density on the metal ion, while the energy position derives from the electron density that is largely delocalized on the water molecules.

These results for hydrated $Na^+$ and $Mg^{2+}$ ions demonstrate the potential of IRD as a more generally applicable method where it is possible to use the solvated ions to selectively probe the solvation shell around each species, thereby opening new possibilities to understand how the solvation shell affects the chemical properties of solutes.

## Methods

### Experimental

The measurements for $Na^+$ and $Mg^{2+}$ in water were performed using a cylindrical liquid jet at the VERITAS beamline at the MAX IV synchrotron radiation facility in Lund, Sweden³⁷. The VERITAS beamline comprises an elliptically polarizing undulator and a collimated plane-grating monochromator with an ellipsoidal refocusing mirror. The beamline is equipped with a large constant-line-spacing grating Rowland spectrometer with a collimating mirror in the non-dispersive direction and a cylindrical grating with 1400 grooves/mm and 67 m radius³⁸. Measurements were done with linear horizontal polarization of the incident radiation.

The samples, 2 M $NaCl$ and 2 M $MgCl_2$ solution, were prepared by dissolving commercially purchased $NaCl$ and $MgCl_2$ (Sigma-Aldrich) with purity of >98% in MilliQ water with a resistivity of 18.2 $MΩ × cm$. The sample was introduced via a vertically mounted liquid jet shooting into a cold trap cooled by liquid nitrogen. The liquid jet was surrounded by a cylindrical differential pumping stage with holes for incident and outgoing photons, allowing the sample to be intersected by the X-rays from the beamline in front of the soft X-ray spectrometer.

Ionization of the 1s level of the $Na^+$ and $Mg^{2+}$ ions was done using photon energies of 1083 and 1318 eV, respectively, i.e., 6 and 8 eV above the respective M 1s thresholds. In both cases, these photon energies are well below the thresholds for the formation of higher-energy states, like M $1s^{-1}2p^{-1}3p$ or M $1s^{-1}2p^{-1}$, ensuring that M $1s^{-1}$ is the only state from which the decay occurs. Emitted photons were detected in the plane of the polarization of the incident radiation at an angle of 90° with respect to the propagation axis. Calibration of the emitted photon energies was done using the energy of the $K_α$ transition, M $1s^{-1} → M 2p^{-1} + hv$, to set the absolute energy at one point of the detector, and a set of elastic peaks to determine the dispersion over the detector. The energy of the $K_α$ transition was determined from the binding energy difference between the M $1s^{-1}$ and M $2p^{-1}$ states of solvated $Na^+$ and $Mg^{2+}$ ions, respectively, obtained by electron spectroscopy¹⁴. For $Na^+$ and $Mg^{2+}$, we thus obtain $hv(K_α)$ as 1041.3 eV (1076.7–35.4) and 1254.3 eV (1309.9–55.6), respectively¹⁴. This procedure relies on the energy difference between the elastic peaks and not their absolute energy, thereby becoming independent of the absolute monochromator energy scale. An upper limit of the spectrometer resolution is the total observed width of the $K_α$ peaks. For both the $Na^+$ and $Mg^{2+}$ edges, this is less than 0.6 eV, which is significantly smaller than the width of the fitted components of the IRD signals (~2.2 eV). The X-ray emission spectra of Na and Mg were measured in the third and second orders of diffraction, respectively.

A pilot study for $Mg^{2+}$ in water was also performed at the P04 beamline of the synchrotron radiation facility PETRA III (DESY, Hamburg, Germany), see Supplementary Note 1.1. The result from P04 is consistent with the results from VERITAS, on which the analysis in this paper is based.

The valence-band photoelectron spectrum of water was recorded at the FlexPES beamline of MAX IV, Lund, Sweden³⁹, using a photon energy of 532 eV. The photon energy resolution was ~0.23 eV, and the spectrometer resolution ~0.25 eV, resulting in a total resolution of ~0.35 eV. To reduce spectral contributions from gas-phase water, a bias voltage of ~50 V was applied to the jet when recording the spectrum.

Preliminary experiments for $Mg^{2+}$ in water were carried out at the P04 beamline of the synchrotron radiation facility PETRA III, DESY, Hamburg⁴⁰, using a setup combining a cylindrical liquid jet and a soft X-ray spectrometer⁴¹. The sample, 2 M $MgCl_2$ solution, was prepared by
dissolving commercially purchased MgCl₂ (Sigma-Aldrich with purity of > 98%) in MilliQ (18.2 MΩ × cm) water. The sample was introduced via a vertically mounted liquid jet shooting into a liquid-nitrogen-cooled cold trap. The liquid jet was intersected by the circularly polarized X-rays from the beamline in front of the soft X-ray spectrometer. The spectrometer was a modified grazing-incidence Rowland-circle spherical grating spectrometer (GRACE) equipped with a 1200 l/mm grating⁴¹. The X-ray emission spectra were measured in fourth order. The beamline and the X-ray detector were protected from the relatively high pressure in the measurement chamber during liquid-jet operation by a differential pumping stage and a filter between the grating and the detector, respectively. Due to the presence of the filter, the internal motions in the spectrometer were impeded. Due to the compromised focusing due to the impeded motions, the spectral resolution was reduced to ~6 eV, which limits the ability to resolve fine substructure in the IRD region.

## Computational

Theoretical investigations were conducted to simulate the experimental spectra and to understand the mechanisms behind the IRD effect; a number of model structures and computational methods have been employed to this end. Two kinds of cluster models were considered: small idealized systems and ion/water clusters inferred by MD simulations of aqueous solutions of NaCl and MgCl₂.

By a reasonable assumption about the interaction between the positive charge of the ions and the dipole of the water molecules, we firstly considered idealized structures in which a few water molecules are oriented with the dipole vector along the M⁹–O axis (R), in a planar geometry for M⁹[H₂O]₄, and in an octahedral geometry for M⁹[H₂O]₆. Using such small systems, as M⁹[H₂O]ₙ, with n = 1, 4, 6, the role of ion-water distance, symmetry, and number of molecules surrounding the cation was studied. Transition energies and radiative transition intensities were computed within the RASSCF level of theory, as implemented in the OpenMolcas package⁴², together with the aug-cc-pVTZ basis set and using the HF wavefunction as reference (HF/RAS). For the different M⁹[H₂O]ₙ cases, the RAS subspaces were built with the M⁹Is orbital in RASI, the three M⁹ 2p orbitals in RAS2, and the 4n water molecular orbitals corresponding to 2a₁, 1b₂, 3a₁, and 1b₁ of the n water molecules in RAS3, in such a way that the ionized system had a total of 2(4 + 4n)–1 active electrons in 4 + 4n orbitals. Relativistic effects based on the Douglas-Kroll decomposition of the electronic Hamiltonian and spin–orbit couplings were included in both calculations.

To go beyond the idealized symmetrical structures and therefore determine a more realistic structural model of Na⁺ or Mg²⁺ ions solvated in water, MD calculations were carried out with the GROMACS code version 5.1.4⁴³ on two samples, both contained in a cubic box with a side length of 50 Å and periodicity conditions in the three dimensions. The sample, which was intended to simulate a situation of extreme dilution, contained a single M⁹ ion and the appropriate number of Cl⁻ counter ions, and 4166 water molecules, while in the second case, the number of M⁹ ions was brought to 150 and the counter ions increased correspondingly, to replicate the concentration of the solution (2 M) used experimentally.

For the M⁹ and Cl⁻ ions, the force field OPLS/AA⁴⁴ has been used, while the water molecule has been simulated with the TIP4P model⁴⁵. Starting from a random distribution of the molecules in the solution, a standard minimization-heating-equilibration protocol was adopted, followed by an extended production run of 10 ns with a time step of 1 fs. Apart from the heating-equilibration steps, the canonical ensemble (NVT) dynamics were conducted at T = 300 K fixed with a velocity-rescaling thermostat. The analysis of the MD simulation carried out on the first sample confirms the result of a previous accurate quantum MD calculation⁴⁶, i.e., that, at high dilution, the M⁹ ion is surrounded by six water molecules with high coordination. In contrast, the second hydration shell is less ordered and more sensitive to rearrangements at

room temperature. From the analysis of the MD trajectories of the larger samples, using the cluster tool of GROMACS, the most typical structures of the solvation shells of M⁹ at the 2 M concentration were identified.

This allowed us to select some specific clusters containing a single M⁹ ion and a variable number of water molecules for which the XES spectrum generated by the ionization of the M⁹ core shell was calculated. For these larger model systems, we adopted the independent particle approximation, where each state is described by a single determinant, namely the less computationally expensive HF and DFT methods, to investigate the dependence of the IRD bands on the number of solvation shells, the orientation of the water molecules closer to the ions, and the presence of ion paring. This was accomplished in the Kohn-Sham DFT ground-state approximation using the StoBe-deMon code⁴⁷, the Perdew-86 density functional, and a TZVP basis set⁴⁸. Such a method can be easily applied to large systems, but is affected by the well-known self-interaction error for the core electron.

In the independent-particle approximation, i.e., neglecting electron correlation, the XES spectrum is obtained by HF calculations for the ground state and the core-hole state. For the latter, the important electronic relaxation is described by optimizing the electronic structure under the constraint of the single occupancy of the core orbital. The final states of the emission process are instead approximated as Koopmans' states, i.e., by removing an electron from one of the frozen valence orbitals of the ground state, thus ignoring the negligible electronic relaxation. The intensity of the radiative process is then calculated, in the long-wavelength approximation, from the transition dipole moment between the initial state and the final states projected onto sets of non-orthogonal orbitals. The HF calculations have been performed by using the DALTON code⁴⁹ with the Ahlrichs-VTZ basis set, while the transition moment has been computed by an homemade code. The restricted open-shell HF approximation, fully including the electron relaxation around the initial core hole, can be considered, despite missing any electron correlation effect, to provide a more reliable simulation of the XES spectra. However, the present calculations found that both DFT methods, with an appropriate energy shift, and HF methods predict quite similar XES spectra.

Comparison of simulated and observed experimental spectra rationalized the mechanisms behind the IRD process. The validity of the one-center approximation for IRD was also considered, see Supplementary Note 2.3, in comparison with full calculations.

## Data availability

The theoretical and processed experimental data needed to re-generate the figures in this study have been deposited at Zenodo Ref. 50. The raw experimental data is available upon request from the authors. The input files and molecular geometries needed to run the theoretical modeling is also available in Ref. 50. Source data are provided with this paper.

## References

1. Atkins, P. W. & Paula, J. Atkin's Physical Chemistry (Oxford University Press, 2002).

2. Fawcett, W. R. *Liquids, Solutions, and Interfaces: From Classical Macroscopic Descriptions to Modern Microscopic Details* (Oxford University Press, 2004).

3. Bunker, B. C. & Casey, W. H. Solvated Ions in Water. In *The Aqueous Chemistry of Oxides* https://doi.org/10.1093/oso/9780199384259.003.0009 (Oxford University Press, 2016).

4. Kropman, M. F. & Bakker, H. J. Dynamics of Water Molecules in Aqueous Solvation Shells. *Science* **291**, 2118–2120 (2001).

5. Laage, D., Elsaesser, T. & Hynes, J. T. Water Dynamics in the Hydration Shells of Biomolecules. *Chem. Rev.* **117**, 10694–10725 (2017).
6. Konstantinovsky, D. et al. Detecting the First Hydration Shell Structure around Biomolecules at Interfaces. ACS Cent. Sci. **8**, 1404–1414 (2022).

7. Bowron, D. T. & Diaz Moreno, S. Using synchrotron X-ray and neutron methods to investigate structural aspects of metal ion solvation and solution structure: An approach using empirical potential structure refinement. Coord. Chem. Rev. **277–278**, 2–14 (2014).

8. Filipponi, A. EXAFS for liquids. J. Phys.: Condens. Matter **13**, R23 (2001).

9. Bakker, H. J. Structural Dynamics of Aqueous Salt Solutions. Chem. Rev. **108**, 1456–1473 (2008).

10. Stöhr, J. *The Nature of X-Rays and Their Interactions with Matter* (Springer, 2023).

11. Jahnke, T. et al. Interatomic and Intermolecular Coulombic Decay. Chem. Rev. **120**, 11295–11369 (2020).

12. Öhrwall, G. et al. Charge Dependence of Solvent-Mediated Inter-molecular Coster-Kronig Decay Dynamics of Aqueous Ions. J. Phys. Chem. B **114**, 17057–17061 (2010).

13. Pokapanich, W. et al. Ionic-Charge Dependence of the Inter-molecular Coulombic Decay Time Scale for Aqueous Ions Probed by the Core-Hole Clock. J. Am. Chem. Soc. **133**, 13430–13436 (2011).

14. Gopakumar, G. et al. Probing aqueous ions with non-local Auger relaxation. Phys. Chem. Chem. Phys. **24**, 8661–8671 (2022).

15. Mudryk, K. et al. How Does $Mg^{2+}_{aq}$ Interact with ATP(aq)? Biomolecular Structure through the Lens of Liquid-Jet Photoemission Spectroscopy. J. Am. Chem. Soc. **146**, 16062–16075 (2024).

16. Dupuy, R. et al. The solvation shell probed by resonant inter-molecular Coulombic decay. Nat. Commun. **15**, 6926 (2024).

17. Persson, I. Structure and size of complete hydration shells of metal ions and inorganic anions in aqueous solution. Dalton Trans. **5**, 15517–15538 (2024).

18. Best, P. E. Electronic Structure of the $MnO_4^-$, $CrO_4^{2-}$, and $VO_4^{3-}$ Ions from the Metal K X-Ray Spectra. J. Chem. Phys. **44**, 3248–3253 (1966).

19. Jones, J. B. & Urch, D. S. Metal-ligand bonding in some vanadium compounds: a study based on X-ray emission data. J. Chem. Soc., Dalton Trans. 1885–1889 (1975).

20. Bergmann, U., Horne, C. R., Collins, T. J., Workman, J. M. & Cramer, S. P. Chemical dependence of interatomic X-ray transition energies and intensities - a study of Mn Kβ" and Kβ$_{2.5}$ spectra. Chem. Phys. Lett. **302**, 119–124 (1999).

21. Pollock, C. J. & DeBeer, S. Valence-to-Core X-ray Emission Spectroscopy: A Sensitive Probe of the Nature of a Bound Ligand. J. Am. Chem. Soc. **133**, 5594–5601 (2011).

22. Lancaster, K. M. et al. X-ray Emission Spectroscopy Evidences a Central Carbon in the Nitrogenase Iron-Molybdenum Cofactor. Science **334**, 974–977 (2011).

23. Pushkar, Y. et al. Direct Detection of Oxygen Ligation to the Mn$_4$Ca Cluster of Photosystem II by X-ray Emission Spectroscopy. Angew. Chem. Int. Ed. **49**, 800–803 (2010).

24. Bloß, D. et al. X-ray radiation damage cycle of solvated inorganic ions. Nat. Commun. **15**, 4594 (2024).

25. Johnsen, R. & Biondi, M. A. Measurements of radiative charge-transfer reactions of doubly and singly charged rare-gas ions with rare-gas atoms at thermal energies. Phys. Rev. A **18**, 996–1003 (1978).

26. Cohen, J. S. & Bardsley, J. N. Calculation of radiative single-charge-transfer cross sections for collisions of $He^{2+}$ with He at low energy. Phys. Rev. A **18**, 1004–1008 (1978).

27. Saito, N. et al. Evidence of radiative charge transfer in argon dimers. Chem. Phys. Lett. **441**, 16–19 (2007).

28. Kreidi, K. et al. Localization of inner-shell photoelectron emission and interatomic Coulombic decay in $Ne_2$. J. Phys. B: At. Mol. Opt. Phys. **41**, 101002 (2008).

29. Yan, S. et al. Observation of interatomic Coulombic decay and electron-transfer-mediated decay in high-energy electron-impact ionization of $Ar_2$. Phys. Rev. A **88**, 042712 (2013).

30. Higuchi, I. et al. Radiative charge transfer and interatomic Coulombic decay following direct double photoionization of neon dimers. J. Phys.: Conf. Ser. **235**, 012015 (2010).

31. Hans, A. et al. Direct evidence for radiative charge transfer after inner-shell excitation and ionization of large clusters. New J. Phys. **20**, 012001 (2018).

32. Hans, A. et al. Electronic Decay of Singly Charged Ground-State Ions by Charge Transfer via van der Waals Bonds. Phys. Rev. Lett. **123**, 213001 (2019).

33. Holzapfel, X. et al. Properties of radiative charge transfer in heterogeneous noble-gas clusters. Phys. Rev. A **105**, 022815 (2022).

34. Thürmer, S. et al. Accurate vertical ionization energy and work function determinations of liquid water and aqueous solutions. Chem. Sci. **12**, 10558–10582 (2021).

35. Winter, B. et al. Full Valence Band Photoemission from Liquid Water Using EUV Synchrotron Radiation. J. Phys. Chem. A **108**, 2625–2632 (2004).

36. Yin, Z. et al. Cationic and Anionic Impact on the Electronic Structure of Liquid Water. J. Phys. Chem. Lett. **8**, 3759–3764 (2017).

37. VERITAS. https://www.maxiv.lu.se/beamlines-accelerators/beamlines/veritas/.

38. Agåker, M. et al. The 10 m collimated Rowland spectrometer at the MAX IV Veritas beamline. J. Synchrotron Radiat. **32**, 1328–1345 (2025).

39. Preobrajenski, A. et al. FlexPES: a versatile soft X-ray beamline at MAX IV (space) Laboratory. J. Synchrotron Radiat. **30**, 831–840 (2023).

40. Viefhaus, J. et al. The Variable Polarization XUV Beamline PO4 at PETRA III: Optics, mechanics and their performance. Nucl. Instrum. Methods Phys. Res. **A** **710**, 151–154 (2013).

41. Nordgren, J. et al. Soft x-ray emission spectroscopy using monochromatized synchrotron radiation. Rev. Sci. Instrum. **60**, 1690–1696 (1989).

42. Fdez. Galván, I. et al. OpenMolcas: From Source Code to Insight. J. Chem. Theory Comput. **15**, 5925–5964 (2019).

43. Abraham, M. J. et al. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to super-computers. SoftwareX **1-2**, 19–25 (2015).

44. Jorgensen, W. L., Maxwell, D. S. & Tirado-Rives, J. Development and Testing of the OPLS All-Atom Force Field on Conformational Energistics and Properties of Organic Liquids. J. Am. Chem. Soc. **118**, 11225–11236 (1996).

45. Jorgensen, W. L., Chandrasekhar, J., Madura, J. D., Impey, R. W. & Klein, M. L. Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. **79**, 926–935 (1983).

46. Di Tommaso, D. Structure and dynamics of the hydrated magnesium ion and of the solvated magnesium carbonates: insights from first principles simulations. Phys. Chem. Chem. Phys. **12**, 894–901 (2010).

47. Hermann, K. et al. StoBedeMon Version 3.3 (2014). StoBe Software: Stockholm, Sweden.

48. Godbout, N., Salahub, D. R., Andzelm, J. & Wimmer, E. Optimization of Gaussian-type basis sets for local spin density functional calculations. Part I. Boron through neon, optimization technique and validation. Can. J. Chem. **70**, 560–571 (1992).

49. Aidas, K. et al. The Dalton quantum chemistry program system. WIREs Comput. Mol. Sci. **4**, 269–284 (2014).

50. Söderström, J. Data for figures https://doi.org/10.5281/zenodo.17107378 (2025).

51. Thürmer, S. et al. Photoelectron Angular Distributions from Liquid Water: Effects of Electron Scattering. Phys. Rev. Lett. **111**, 173005 (2013).
## Acknowledgements

We acknowledge MAX IV Laboratory for time on the VERITAS beamline under proposal 2022-1207. Research conducted at MAX IV, a Swedish national user facility, is supported by the Swedish Research Council under contract 2018-07152, the Swedish Governmental Agency for Innovation Systems under contract 2018-04969, and FORMAS under contract 2019-02496. We acknowledge DESY (Hamburg, Germany), a member of the Helmholtz Association HGF, for the provision of experimental facilities. Parts of this research were carried out at PETRA III and would like to thank Moritz Hoesch and his team for assistance in using beamline PO4. Beamtime was allocated for proposal I-20211465 EC. O.B., J.-E.R., and H.Å. acknowledge funding from the Swedish Research Council (VR) through the projects VR 2023-04346, VR 2021-04017, and VR 2022-03405, respectively, as well as the Swedish Foundation for International Cooperation in Research and Higher Education (STINT) through project 202100-2932. The authors also thank the Swedish National Infrastructure for Computing (NAISS 2022/3-34) at the National Supercomputer Center of Linköping University (Sweden) partially funded by the Swedish Research Council through grant agreement no. 2018-05973. M.H. acknowledges funding by the German Federal Ministry of Education and Research under grant number 13K22XXB DYLIXUT. D.B. and A.H. acknowledge support by the German Federal Ministry of Education and Research (BMBF) through project 05K22RK1-TRANSALP. F.T. acknowledges funding by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) - Project No. 509471550, Emmy Noether Program. A.N.B. acknowledges funding from FAPESP-Brazil - Proj. 2017/11986-5, 2024/00998-6 and FAPEX-UNICAMP 3358/24, 2684/25. A.N.B. and L.M.C. acknowledge FAPESP 2021/06527-7 and 2020/04822-9. R.M. acknowledges funding from the Federal District Research Foundation (FAPDF) for the 410/2022 project.

## Author contributions

J.S., J.-E.R., and O.B. conceived the experiment. J.S., M.A., D.B., F.T., I.I., D.V., J.P., J.-E.R., M.H., Z.Y., J.N., and O.B. performed the initial measurements at PO4. J.S., V.E., O.B., A.N.d.B., M.A., R.M., T.T., C.S., Y.-P.C., and A.G. performed the main experiment at VERITAS. G.Ö. provided the photoemission data. J.S. did the data analysis. L.M.C., V.C., and H.Å. performed the theoretical modeling. J.S., O.B., L.M.C, V.C., and H.Å. wrote the manuscript. J.S., M.A., J.-E.R., A.N.d.B., Z.Y., A.H., J.N., L.M.C., V.C., H.Å., and O.B. discussed the interpretation and provided input on the manuscript.

## Funding

Open access funding provided by Uppsala University.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41467-025-65581-7.

**Correspondence** and requests for materials should be addressed to Johan Söderström or Olle Björneholm.

**Peer review information** *Nature Communications* thanks Majed Chergui, and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. A peer review file is available.

**Reprints and permissions information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2025