Interatomic Coulombic Decay (ICD)

• Interatomic Coulombic Decay (ICD) is a non-radiative, correlated energy-transfer process where an excited atom or molecule transfers excess energy to a neighbor, resulting in ionization and low-energy electron emission.
• ICD exhibits a strong dependence on interatomic distance, typically following an R⁻⁶ decay law, with its rate further influenced by nuclear motion and the local environment.
• ICD has broad applications in probing ultrafast dynamics in systems such as rare gas clusters, helium nanodroplets, heteronuclear dimers, and engineered semiconductor nanostructures.

Interatomic Coulombic Decay (ICD) is a non-radiative, ultrafast electronic relaxation process in which an excited atom or molecule transfers excess energy to a neighboring center, usually ionizing that neighbor and emitting a low-energy electron. First proposed by Cederbaum et al. in 1997, ICD has been observed in rare gas clusters, water droplets, surfaces, and systems with van-der-Waals or hydrogen bonds, and later analyzed in helium nanodroplets, ion-dimer collisions, quantum-well nanostructures, and gases (Trinter et al., 2013, Falkowski et al., 9 Mar 2026). Across these settings, ICD appears as a family of correlated decay channels rather than a single microscopic mechanism: it can follow inner-valence ionization, resonant Auger decay, multiple resonant excitation, double excitation, electron-impact excitation, or collision-induced electron removal and excitation.

1. Fundamental process and energetic formulation

In its canonical form, ICD occurs when an electronically excited or inner-valence-ionized center cannot efficiently relax locally, but can instead transfer energy to a neighbor. In Ne$_2$, a $2s$ vacancy on one atom is filled by a $2p$ electron, and the excess energy ionizes the neighboring atom, producing Ne$^+(2p^{-1})$–Ne$^+(2p^{-1})$ and an ICD electron (Schnorr et al., 2013). In argon dimers after resonant core excitation and spectator-type resonant Auger decay, the unstable states $3p^{-2}3d$ and $3p^{-2}4d$ decay by transferring energy to the neighboring atom, yielding two Ar$^+$ ions and a low-energy ICD electron (Kimura et al., 2013). In helium nanodroplets, even a doubly excited helium atom can decay non-locally according to

$$\mathrm{He}^{**}\mathrm{He} \to \mathrm{He}^*\mathrm{He}^+ + e^- ,$$

demonstrating that ICD can compete with extremely fast intra-atomic autoionization (Bastian et al., 2023).

A common representation treats the ICD precursor as an electronic resonance,

$E_\text{res} = E_R - i\frac{\Gamma}{2},$

with the lifetime related to the decay width by

$2s$0

For helium-pair calculations in droplets, the same relation is written as

$2s$1

where $2s$2 depends strongly on the interatomic separation $2s$3 (Parravicini et al., 2023, LaForge et al., 2020). This formalism is central because ICD precursors are metastable continuum-embedded states rather than bound eigenstates.

Several experimentally important variants fall under the ICD umbrella. Resonant ICD in helium droplets proceeds after multiple resonant excitations of He atoms (LaForge et al., 2020). Site-specific ICD in Ar$2s$4 can follow resonant Auger decay after a $2s$5 excitation, producing slow electrons at a selected atomic site (Kimura et al., 2013). Two-center resonant photoionization in HeNe uses He as an “antenna” and Ne as a “receiver,” with resonant absorption on He leading to ICD-mediated ionization of Ne (Trinter et al., 2013). These examples establish ICD as a correlation-driven energy-transfer process whose trigger can be inner-valence ionization, core excitation, or resonant excitation of neutral centers.

2. Distance dependence, nuclear motion, and time-domain measurements

A defining property of ICD is its strong dependence on geometry. For large separations, the decay width follows the characteristic law

$2s$6

and in the helium dimer the kinetic energy release is related to the internuclear distance at the moment of decay by

$2s$7

in atomic units (Trinter et al., 2013). Because nuclear motion can proceed on the same timescale as the decay, ICD is often a coupled electron-nuclear process rather than a purely electronic one.

The first direct time-resolved measurement of ICD was performed on He$2s$8 using a post-collision-interaction streaking method. In this system, the helium dimer has an exceptionally weak bond of $2s$9 neV and a mean interatomic separation of $2p$0 Å, yet ICD still occurs. The experiment showed that ICD can happen at any internuclear distance along the contraction of the excited-state wavepacket, and that the probability for ICD increases as the nuclei approach, consistent with the $2p$1 dependence at large distances (Trinter et al., 2013).

A pump-probe experiment on Ne$2p$2 determined the lifetime of the Ne$2p$3 ICD precursor to be $2p$4 fs, in agreement with quantum calculations that include nuclear wavepacket dynamics (Schnorr et al., 2013). A later time-resolved study of two-photon doubly excited Ne$2p$5 identified two classes of ICD precursor states: long-lived states with a lifetime of $2p$6 fs and short-lived states with lifetime less than $2p$7 fs (Takanashi et al., 2019). These measurements established that even within one nominal system, ICD lifetimes can vary strongly with electronic symmetry and the nuclear coordinate sampled before decay.

State resolution sharpened this picture in HeNe. There, vibrationally resolved linewidths yielded ICD lifetimes ranging from about $2p$8 fs for low vibrational states to $2p$9 fs for the $^+(2p^{-1})$0 level, and the lifetime was found to strongly increase with increasing vibrational state (Trinter et al., 2013). The interpretation given in the data is that higher vibrational excitation samples larger mean internuclear distances, thereby reducing the decay width.

Helium nanodroplets showed that distance dependence can be dynamically reshaped by the environment. Standard virtual-photon ICD models predicted time constants in the tens to hundreds of picoseconds, with $^+(2p^{-1})$1 ps for the fastest channel, but direct time-resolved measurements of resonant ICD in multiply excited He droplets found decay times as fast as $^+(2p^{-1})$2 fs and weaker-than-expected dependence on droplet size and laser intensity (LaForge et al., 2020). Time-dependent density functional theory and ab initio calculations traced this acceleration to bubble formation and merging: pairs of excited helium atoms strongly attract each other, their void bubbles overlap and merge, and the atoms are driven to $^+(2p^{-1})$3 Å in under $^+(2p^{-1})$4 fs, where $^+(2p^{-1})$5 becomes very large (LaForge et al., 2020). This result directly contradicted a static-separation picture of the process.

3. Helium nanodroplets as an ICD medium

Helium nanodroplets provide an especially rich ICD environment because the excitations, ions, and neutral background are all mobile. In pure droplets, ICD can be induced by photoexciting the $^+(2p^{-1})$6 state of He$^+(2p^{-1})$7 with XUV or VUV radiation, after which the excited ionic state relaxes non-locally and creates two neighboring He$^+(2p^{-1})$8 ions plus an ICD electron (Shcherbinin et al., 2017, Wiegandt et al., 2018). Coincidence imaging showed that the electron spectra resemble the binary He dimer process, but the ionic products depend strongly on where in the droplet the decay occurs.

Three spatial regimes were distinguished. Surface ICD gives energetic He$^+(2p^{-1})$9 ions with kinetic-energy peaks around $^+(2p^{-1})$0–$^+(2p^{-1})$1 eV, close to free-dimer values. Subsurface ICD produces a broad $^+(2p^{-1})$2 eV He$^+(2p^{-1})$3 feature caused by elastic “billiard ball” collisions with surrounding He atoms that interrupt the Coulomb explosion and restart it at larger separation. Interior ICD stops the ions through multiple collisions, leading to trapped ions and the formation of He$^+(2p^{-1})$4 complexes (Shcherbinin et al., 2017). Direct observation of ICD and subsequent ion-atom scattering in pure $^+(2p^{-1})$5He nanoclusters further showed that a main energy-loss mechanism is a single hard binary collision with one atom of the cluster (Wiegandt et al., 2018).

Resonant ICD in multiply excited droplets added a distinct regime. Across droplet sizes $^+(2p^{-1})$6 atoms and FEL pulse energies from $^+(2p^{-1})$7 to $^+(2p^{-1})$8J, the measured decay times remained ultrafast and only weakly dependent on excitation density, despite earlier predictions of strong sensitivity to interatomic separation (LaForge et al., 2020). Only about $^+(2p^{-1})$9–$3p^{-2}3d$0 of excited helium atoms underwent ICD in the relevant excitation range, with the efficiency rising with excitation density; incomplete depletion at longer delays indicated competing channels such as ejection of excited atoms or excimer formation (LaForge et al., 2020).

The helium-droplet data therefore show two linked features of ICD in a quantum fluid. First, the primary decay remains the same correlated energy-transfer process identified in dimers. Second, the host medium redistributes kinetic energy, changes pair separations dynamically, and opens secondary channels that strongly affect observables without eliminating the underlying ICD signature.

4. Secondary processes, resonant selectivity, and doped droplets

A major development has been the recognition that ICD in droplets is often induced indirectly by secondary electron dynamics. In large He nanodroplets irradiated with XUV photons above the ionization threshold, primary photoionization creates a photoelectron that inelastically scatters and excites a neutral He atom to He$3p^{-2}3d$1, then slows by elastic scattering and can recombine with the parent He$3p^{-2}3d$2 ion to create a second He$3p^{-2}3d$3. These two excited species can meet at the droplet surface and undergo ICD (Ltaief et al., 2023). For droplets with radius $3p^{-2}3d$4 nm, this indirect ICD becomes the main electron-emission channel, and the normalized ICD electron yield reaches $3p^{-2}3d$5 for $3p^{-2}3d$6 nm (Ltaief et al., 2023).

Spectroscopically resolved studies established that this physics extends both below and above threshold. Below the ionization threshold, resonant excitation of multiple centers induces resonant ICD in large droplets even with weak synchrotron radiation. Above threshold, photoelectron scattering and electron-ion recombination populate excited states that subsequently decay by ICD, demonstrating the importance of secondary processes in extended condensed-phase systems (Ltaief et al., 2023). The same study concluded that ICD can serve as a diagnostic tool for monitoring the relaxation dynamics of highly excited and ionized weakly bound nanosystems (Ltaief et al., 2023).

The resonance selectivity of ICD is especially clear in doubly excited helium. High-resolution photoelectron spectra of He nanodroplets excited to the $3p^{-2}3d$7 Fano resonance provided the first experimental observation of ICD induced by double excitation of helium. Even though the doubly excited He atom has an autoionization lifetime of $3p^{-2}3d$8 fs, ICD still proceeds, yielding He$3p^{-2}3d$9He$3p^{-2}4d$0 atom-pair states and an ICD probability relative to total ionization of $3p^{-2}4d$1 (Bastian et al., 2023). The data showed that ICD is absent off resonance, which was interpreted as element- and state-specific spectral selectivity (Bastian et al., 2023).

Doped droplets reveal additional ICD channels and a nontrivial competition between mechanisms. In Li-doped large He nanodroplets exposed to XUV photons above $3p^{-2}4d$2 eV, Li ions are efficiently produced by an ICD process involving metastable He atoms and He$3p^{-2}4d$3 excimers populated by elastic and inelastic scattering of photoelectrons and by electron-ion recombination. This indirect ICD is more efficient than Li ionization by ICD following direct resonant photoexcitation at $3p^{-2}4d$4 eV and by charge-transfer ionization, exceeding them by more than an order of magnitude in the inelastic-scattering regime (Ltaief et al., 2024). In metal-doped He droplets with Li or Rb, coincidence spectroscopy and ab initio Fano-CI-Stieltjes calculations further showed that charge-exchange ICD dominates both short-range and long-range autoionization up to about $3p^{-2}4d$5 Å, while the virtual-photon contribution is negligible in that system (Ltaief et al., 2019). This constitutes an important qualification of the common identification of long-range autoionization with direct virtual-photon ICD.

5. Collisions, heteronuclear dimers, and engineered systems

ICD is also a prominent source of low-energy electrons in collisions. In slow collisions between O$3p^{-2}4d$6 and Ne$3p^{-2}4d$7, single $3p^{-2}4d$8 electron capture from one site of the dimer creates a Ne$3p^{-2}4d$9 state that decays by ICD, producing electrons in the $^+$0–$^+$1 eV range with peaks near $^+$2 eV. The ICD signal was one order of magnitude stronger than the radiative charge-transfer peak and overtook Auger electron emission by the scattered projectiles after double-electron capture (Iskandar et al., 2014). By contrast, analogous measurements with Ar$^+$3 and Xe$^+$4 projectiles showed no evidence of inner-shell single-electron capture and no observable ICD (Iskandar et al., 2014).

A theoretical study of slow He$^+$5-Ne$^+$6 collisions predicted an even stronger ICD regime. Using an independent-atom-independent-electron model with coupled-channel two-center basis generator calculations, the work found that below $^+$7 keV/amu, ICD is the dominant Ne$^+$8+Ne$^+$9 fragmentation process and the ratio of ICD to direct electron emission is predicted to reach $$\mathrm{He}^{**}\mathrm{He} \to \mathrm{He}^*\mathrm{He}^+ + e^- ,$$0–$$\mathrm{He}^{**}\mathrm{He} \to \mathrm{He}^*\mathrm{He}^+ + e^- ,$$1 (Kirchner, 2021). At $$\mathrm{He}^{**}\mathrm{He} \to \mathrm{He}^*\mathrm{He}^+ + e^- ,$$2 keV/amu, ICD accounts for $$\mathrm{He}^{**}\mathrm{He} \to \mathrm{He}^*\mathrm{He}^+ + e^- ,$$3–$$\mathrm{He}^{**}\mathrm{He} \to \mathrm{He}^*\mathrm{He}^+ + e^- ,$$4 of the total Ne$$\mathrm{He}^{**}\mathrm{He} \to \mathrm{He}^*\mathrm{He}^+ + e^- ,$$5+Ne$$\mathrm{He}^{**}\mathrm{He} \to \mathrm{He}^*\mathrm{He}^+ + e^- ,$$6 fragmentation yield depending on the response model (Kirchner, 2021).

Collision-induced ICD in argon dimers involves excitation rather than only vacancy production. In helium-ion–argon-dimer collisions, configurations of the form Ar$$\mathrm{He}^{**}\mathrm{He} \to \mathrm{He}^*\mathrm{He}^+ + e^- ,$$7, especially Ar$$\mathrm{He}^{**}\mathrm{He} \to \mathrm{He}^*\mathrm{He}^+ + e^- ,$$8, were identified as dominant pathways to ICD, with other $$\mathrm{He}^{**}\mathrm{He} \to \mathrm{He}^*\mathrm{He}^+ + e^- ,$$9–$E_\text{res} = E_R - i\frac{\Gamma}{2},$0 excitations also contributing significantly. The calculations also found that a He$E_\text{res} = E_R - i\frac{\Gamma}{2},$1 projectile offers a strong pathway for ICD as the projectile impact energy decreases (Starko et al., 24 Mar 2026).

Heteronuclear dimers provide high spectroscopic selectivity. In HeNe, resonant coupling on the He side enhances HeNe$E_\text{res} = E_R - i\frac{\Gamma}{2},$2 ion production by more than a factor of $E_\text{res} = E_R - i\frac{\Gamma}{2},$3, directly demonstrating two-center resonant photoionization and yielding vibrationally and electronically resolved ICD widths (Trinter et al., 2013). In Ar$E_\text{res} = E_R - i\frac{\Gamma}{2},$4, resonant Auger precursors $E_\text{res} = E_R - i\frac{\Gamma}{2},$5 and $E_\text{res} = E_R - i\frac{\Gamma}{2},$6 both undergo ICD, but with very different timescales: $E_\text{res} = E_R - i\frac{\Gamma}{2},$7–$E_\text{res} = E_R - i\frac{\Gamma}{2},$8 fs for $E_\text{res} = E_R - i\frac{\Gamma}{2},$9 and $2s$00–$2s$01 ps for $2s$02, leading to different KER signatures because the slower state decays after nuclear contraction (Kimura et al., 2013).

Beyond atomic and molecular aggregates, ICD has been proposed in semiconductor nanostructures. In two coupled quantum wells modeled with a one-dimensional effective potential and a single-band effective-mass approximation, the ICD lifetime is on the picosecond timescale and decreases with decreasing inter-well distance. At a tuned distance $2s$03 Å, the emitted ICD electron is trapped in a one-electron resonance state, shortening the lifetime by an order of magnitude even at long distance (Goldzak et al., 2015). This demonstrates that ICD can be engineered through continuum-state structure as well as geometry.

6. Theoretical descriptions, three-body effects, and extended-range ICD

Because ICD precursors are resonances, their description requires non-Hermitian or continuum-adapted electronic-structure methods. Equation-of-motion coupled-cluster theory combined with complex basis functions and Feshbach-Fano projection has been applied to Ne$2s$04, NeAr, NeMg, and $2s$05, covering outgoing-electron energies between $2s$06 eV and $2s$07 eV. Both methods give better results when the outgoing electron is fast, but the characteristic $2s$08 distance dependence of the ICD width is captured much better with complex basis functions (Parravicini et al., 2023). This comparison is methodologically important because it ties the observed distance law directly to how the continuum electron is represented.

Three-body ICD introduces an environmental scattering problem. A virtual-photon description based on Green’s tensors showed that a passive mediator atom can substantially enhance or suppress the ICD rate at large distances. In the non-retarded regime, the mediator modifies the rate through geometry- and polarizability-dependent terms; in the retarded regime, interference between direct and mediated virtual-photon pathways can produce oscillatory enhancement or suppression depending on the mediator position (Bennett et al., 2018). The relevant framework writes the rate in terms of the electromagnetic Green’s tensor,

$2s$09

and thereby makes ICD explicitly environment-dependent (Bennett et al., 2018).

Recent work has broadened the range of ICD even further by identifying a hitherto unrecognized intermolecular Coulombic decay mechanism in gases. In that regime, the short-range Coulombic term becomes negligible, while retardation effects lead to a rate expression containing $2s$10, $2s$11, and $2s$12 terms, with the $2s$13 contribution dominating at micrometer scales (Falkowski et al., 9 Mar 2026). The paper concluded that ICD can be efficiently active in atomic and molecular gases despite very large distances between units, and analyzed population dynamics in which ICD can initially be much faster than radiative decay (Falkowski et al., 9 Mar 2026).

Taken together, these theoretical developments correct several oversimplified views. ICD is not restricted to weakly bound dimers at fixed geometry; not all long-range ICD is virtual-photon dominated; and not all environment-induced autoionization at large distance can be reduced to a single distance law. A plausible implication is that ICD should now be understood as a hierarchy of correlated decay regimes—Coulombic, charge-exchange, mediator-modified virtual-photon, and retardation-dominated—whose relative importance is set by electronic structure, continuum energetics, geometry, and medium response (Ltaief et al., 2019, Bennett et al., 2018, Falkowski et al., 9 Mar 2026).