On the interatomic Coulombic decay in the Ne dimer

Observation of Nuclear Wave-Packet Interference in Ultrafast Interatomic Energy Transfer

We report the experimental observation of quantum interference in the nuclear wave-packet dynamics driving ultrafast excitation-energy transfer in argon dimers below the threshold of interatomic Coulombic decay (ICD). Using time-resolved photoion-photoion coincidence spectroscopy and quantum dynamics simulations, we reveal that the electronic relaxation dynamics of the inner-valence 3s hole on one atom leading to a 4s or 4p excitation on the other one is influenced by nuclear quantum dynamics in the initial state, giving rise to a deep, periodic modulation on the kinetic-energy-release (KER) spectra of the coincident Ar^{+}-Ar^{+} ion pairs. Moreover, the time-resolved KER spectra show characteristic fingerprints of quantum interference effects during the energy-transfer process. Our findings pave the way to elucidating quantum-interference effects in ultrafast charge- and energy-transfer dynamics in more complex systems, such as molecular clusters and solvated molecules.

Effect of Protonation and Deprotonation on Electron Transfer Mediated Decay and Interatomic Coulombic Decay

Electronically excited atoms or molecules in an environment are often subject to interatomic/intermolecular Coulombic decay (ICD) and/or electron transfer mediated decay (ETMD) mechanisms. A few of the numerous variables that can impact these non-radiative decay mechanisms include bond distance, the number of nearby atoms or molecules, and the polarisation effect. In this paper, we have studied the effect of protonation and deprotonation on the ionization potential (IP), double ionization potential (DIP), and lifetime (or decay width) of the temporary bound state in these non-radiative decay processes. We have chosen LiH-NH3 and LiH-H2 O as test systems. The equation of motion coupled cluster singles and doubles method augmented by complex absorbing potential (CAP-EOM-CCSD) has been used in calculating the energetic position of the decaying state and the system's decay rate. Deprotonation of LiH-NH3 /LiH-H2 O either from the metal center (LiH) or from ammonia/water lowers the IP and DIP compared to the neutral systems. In contrast, protonation increases these quantities compared to neutral systems. The protonation closes the inner valence state relaxation channels for ICD/ETMD. For example, the decay of the O-2s/N-2s state stops in protonated systems (LiH2 + -H2 O, LiH2 + -NH3 , and LiH-NH4 + ). Our study also shows that the efficiency, i. e., the rate of ICD/ETMD, can be altered by protonation and deprotonation. It is expected to have implications for chemical and biological systems.

Complex energies and transition dipoles for shape-type resonances of uracil anion from stabilization curves via Padé

Absorption of slow moving electrons by neutral ground state nucleobases has been known to produce resonance metastable states. There are indications that such metastable states may play a key role in DNA/RNA damage. Therefore, herein, we present an ab initio non-Hermitian investigation of the resonance positions and decay rates for the low lying shape-type states of the uracil anion. In addition, we calculate the complex transition dipoles between these resonance states. We employ the resonance via Padé (RVP) method to calculate these complex properties from real stabilization curves by analytical dilation into the complex plane. This method has already been successfully applied to many small molecular systems, and herein, we present the first application of RVP to a medium-sized system. The presented resonance energies are optimized with respect to the size of the basis set and compared with previous theoretical studies and experimental findings. Complex transition dipoles between the shape-type resonances are computed using the optimal basis set. The ability to calculate ab initio energies and lifetimes of biologically relevant systems paves the way for studying reactions of such systems in which autoionization takes place, while the ability to also calculate their complex transition dipoles opens the door for studying photo-induced dynamics of such biological molecules.

Ultrafast temporal evolution of interatomic Coulombic decay in NeKr dimers

We investigate interatomic Coulombic decay in NeKr dimers after neon inner-valence photoionization [Ne⁺(2s^-1)] using a synchrotron light source. We measure with high energy resolution the two singly charged ions of the Coulomb-exploding dimer dication and the photoelectron in coincidence. By carefully tracing the post-collision interaction between the photoelectron and the emitted ICD electron we are able to probe the temporal evolution of the state as it decays. Although the ionizing light pulses are 80 picoseconds long, we determine the lifetime of the intermediate dimer cation state and visualize the contraction of the nuclear structure on the femtosecond time scale.

Decay Processes in Cationic Alkali Metals in Microsolvated Clusters: A Complex Absorbing Potential Based Equation-of-Motion Coupled Cluster Investigation

We have employed the highly accurate complex absorbing potential based ionization potential equation-of-motion coupled cluster singles and doubles (CAP-IP-EOM-CCSD) method to study the various intermolecular decay processes in ionized metals (Li⁺, Na⁺, K⁺) microsolvated by water molecules. For the Li atom, the electron is ionized from the 1s subshell. However, for Na and K atoms, the electron is ionized from 2s and both 2s and 2p subshells, respectively. We have investigated decay processes for the Li⁺-(H₂O)_n (n = 1-3) systems, as well as Na⁺-(H₂O)_n (n = 1, 2), and K⁺-H₂O. The lithium cation in water can decay only via electron transfer mediated decay (ETMD) as there are no valence electrons in lithium. We have investigated how the various decay processes change in the presence of different alkali metal atoms and how the increasing number of water molecules play a significant role in the decay of microsolvated systems. To see the effect of the environment, we have studied Li⁺-NH₃ in comparison to Li⁺-H₂O. In the case of Na⁺-H₂O, we have studied the impact of bond distance on the decay width. The effect of polarization on decay width was checked for the X⁺-H₂O (X = Li, Na) systems. We used the PCM model to study the polarization effect. We have compared our results with existing theoretical and experimental results wherever available in the literature.

Uniform vs Partial Scaling within Resonances via Padé Based on the Similarities to Other Non-Hermitian Methods: Illustration for the Beryllium 1s22p3s State

Resonance via Padé (RVP) is an efficient method for calculating autoionization resonance states. It is based on the stabilization technique in which the basis set is scaled. The scaling can be uniform (i.e., all basis functions are scaled) or partial. Herein, we compare the two RVP scaling schemes for calculating an autoionization eigenvalue; moreover, the effect of freezing the core electrons is intertwined within this comparison. In order to study the different behavior of the RVP schemes, we associate each RVP scaling scheme with a complex contour of integration. Similarities between RVP and other non-Hermitian methods emerge from the generated contours, which suggest that RVP introduces similar outgoing boundary conditions as the complex scaling (CS), complex basis function (CBF), and reflection-free complex absorbing potential (RF-CAP) methods. A uniform-RVP contour, unlike a partial one, immediately penetrates the complex plane and influences the interaction region. Hence, uniform scaling within RVP destroys the description of the core electrons, as well as the description of the reference state, and yields less reliable results than partial scaling. The 1s²2p3s ¹P autoionization state of Be, at the equation-of-motion coupled-cluster level, is used as our case study model.

Interatomic and Intermolecular Coulombic Decay

Interatomic or intermolecular Coulombic decay (ICD) is a nonlocal electronic decay mechanism occurring in weakly bound matter. In an ICD process, energy released by electronic relaxation of an excited atom or molecule leads to ionization of a neighboring one via Coulombic electron interactions. ICD has been predicted theoretically in the mid nineties of the last century, and its existence has been confirmed experimentally approximately ten years later. Since then, a number of fundamental and applied aspects have been studied in this quickly growing field of research. This review provides an introduction to ICD and draws the connection to related energy transfer and ionization processes. The theoretical approaches for the description of ICD as well as the experimental techniques developed and employed for its investigation are described. The existing body of literature on experimental and theoretical studies of ICD processes in different atomic and molecular systems is reviewed.

Quantum Effects Dominating the Interatomic Coulombic Decay of an Extreme System

LiHe is an extreme open-shell system. It is among the weakest bound systems known, and its mean interatomic distance extends dramatically into the classical forbidden region. Upon 1s → 2p excitation of He, interatomic Coulombic decay (ICD) takes place in which the electronically excited helium atom relaxes and transfers its excess energy to ionize the neighboring lithium atom. A substantial part of the decay is found to be to the dissociation continuum producing Li⁺ and He atoms. The distribution of the kinetic energy released by the ICD products is found to be highly oscillatory. Its analysis reveals that quantum phase shifts between the decaying states and the dissociating final states are controlling this ICD reaction. The semiclassical reflection principle, which commonly explains ICD reactions, fails. The process is expected to be amenable to experiment.

Fano-ADC(2,2) method for electronic decay rates

Fano-ADC is a family of ab initio methods for the prediction of electronic decay widths in excited, singly and doubly ionized systems. It has been particularly successful in elucidating the geometry dependence of the inter-atomic decay widths in clusters and facilitated the prediction of new electronic decay phenomena. However, the available Fano-ADC schemes are limited to the second-order treatment of the initial state and the first-order treatment of the final states of the decay. This confines the applicability of the Fano-ADC approach to first-order decay processes, e.g., normal but not double Auger decay (DAD), and compromises the numerical accuracy of the schemes through the unbalanced treatment of electronic correlation. Here, we introduce the ADC(2,2) approximation for singly ionized states, which describes both initial and final states of the decay up to second order. We use the new scheme to construct the Fano-ADC(2,2) approximation for the decay widths and show that it provides superior accuracy for the decay widths of a series of processes. Moreover, the Fano-ADC(2,2) method provides access to second-order decay processes, such as DAD, which are qualitatively beyond the reach of the previously available Fano-ADC implementations.

Breakdown of the Spectator Concept in Low-Electron-Energy Resonant Decay Processes

We suggest that low-energy electrons, released by resonant decay processes, experience substantial scattering on the electron density of excited electrons, which remain a spectator during the decay. As a result, the angular emission distribution is altered significantly. This effect is expected to be a common feature of low-energy secondary electron emission. In this Letter, we exemplify our idea by examining the spectator resonant interatomic Coulombic decay of Ne dimers. Our theoretical predictions are confirmed by a corresponding coincidence experiment.

Interatomic Coulombic Decay Mediated by Ultrafast Superexchange Energy Transfer

Inner-valence ionized states of atoms and molecules live shorter if these species are embedded in an environment due to the possibility for ultrafast deexcitation known as interatomic Coulombic decay (ICD). In this Letter we show that the lifetime of these ICD active states decreases further when a bridge atom is in proximity to the two interacting monomers. This novel mechanism, termed superexchange ICD, is an electronic correlation effect driven by the efficient energy transfer via virtual states of the bridge atom. The superexchange ICD is discussed in detail on the example of the NeHeNe trimer. We demonstrate that the decay width of the Ne^{+}(2s^{-1})  ^{2}Σ_{g}^{+} resonance increases 6 times in the presence of the He atom at a distance of 4 Å between the two Ne atoms. Using a simple model, we provide a qualitative explanation of the superexchange ICD and we derive analytical expressions for the dependence of the decay width on the distance between the neon atoms.

Impact of metal ion's charge on the interatomic Coulombic decay widths in microsolvated clusters

Interatomic Coulombic decay (ICD) is an efficient electronic decay mechanism of electronically excited atoms and molecules embedded in an environment. For the series of isoelectronic Na(+), Mg(2+), and Al(3+) ions in aqueous solution, ultrashort ICD lifetimes of 3.1 fs, 1.5 fs, and 0.9 fs, respectively, were observed experimentally. The magnitude of the ICD lifetimes and their variation within the series were qualitatively explained by shortening metal-oxygen equilibrium distances and the increasing polarization of the water molecules as the metal's charge grows. We carried out an extensive ab initio investigation of the variation of the ICD widths with the metal-oxygen distances and the number of water neighbors in Na(+)-(H2O)m (m = 1-4) and Mg(2+)-(H2O)n (n = 1-6) clusters including and excluding polarization effects in the decaying state. We demonstrated that the effect of the induced polarization of the water ligand and the equilibrium cation-oxygen distance are equally important in determining the ordering and ratios of the ICD lifetimes in the series. Moreover, we showed that the induced polarization of the water molecules leads to a slower than linear growth of ICD width with the number of equivalent water neighbors; the non-linearity is stronger for Mg(2+). The ab initio ICD widths in microsolvated Na(+)-(H2O)4 and Mg(2+)-(H2O)6 clusters are found to be in good agreement with the experimental values.

Strong enhancement of cage effects in water photolysis caused by interatomic Coulombic decay

The impact of the solvent on the photodissociation of embedded molecules has been intensively investigated in the last decades. Collisions of photofragments with the solvating atoms or molecules can change their kinetic energy distribution or even lead to the de-excitation of the dissociating molecule to a bound electronic state quenching the dissociation. In this article we show that this cage effect is strongly enhanced if interatomic Coulombic decay (ICD) of the excited state becomes allowed. Ab initio calculations in H2O-Cl(-) cluster show that the ultra-fast dissociation of water in the à excited state is strongly quenched by ICD. We found that this very efficient quenching is due to two factors. First, the lifetimes of the à state due to ICD are short ranging between 6 and 30 fs. Second, nuclear dynamics is dominated by the chattering motion of the H atom between O and Cl(-) allowing ICD to act for longer times. We hope that this work will be an important first step in clarifying the impact of ICD on photodissociation of embedded molecules.

Time-resolved measurement of interatomic coulombic decay in Ne2

The lifetime of interatomic Coulombic decay (ICD) [L. S. Cederbaum et al., Phys. Rev. Lett. 79, 4778 (1997)] in Ne2 is determined via an extreme ultraviolet pump-probe experiment at the Free-Electron Laser in Hamburg. The pump pulse creates a 2s inner-shell vacancy in one of the two Ne atoms, whereupon the ionized dimer undergoes ICD resulting in a repulsive Ne+(2p(-1))-Ne+(2p(-1)) state, which is probed with a second pulse, removing a further electron. The yield of coincident Ne+-Ne2+ pairs is recorded as a function of the pump-probe delay, allowing us to deduce the ICD lifetime of the Ne2(+)(2s(-1)) state to be (150±50)  fs, in agreement with quantum calculations.

Kinetic energy release in fragmentation processes following electron emission: a time-dependent approach

A time-dependent approach for the kinetic energy release (KER) spectrum is developed for a fragmentation of a diatomic molecule after an electronic decay process, e.g., Auger process. It allows one to simulate the time-resolved spectra and provides more insight into the molecular dynamics than the time-independent approach. Detailed analysis of the time-resolved emitted electron and KER spectra sheds light on the interrelation between wave packet dynamics and spectra.

Interrelation between the distributions of kinetic energy release and emitted electron energy following the decay of electronic states

In an electronic decay process followed by fragmentation the kinetic energy release and electron spectra can be measured. Classically they are the mirror image of each other, a fact which is often used in practice. Quantum expressions are derived for both spectra and analyzed. It is demonstrated that these spectra carry complementary quantum information and are related to the nuclear dynamics in different participating electronic states. Illustrative examples show that the classical picture of a mirror image can break down and shed light on the underlying physics.

Enhanced production of low energy electrons by alpha particle impact

Radiation damage to living tissue stems not only from primary ionizing particles but to a substantial fraction from the dissociative attachment of secondary electrons with energies below the ionization threshold. We show that the emission yield of those low energy electrons increases dramatically in ion-atom collisions depending on whether or not the target atoms are isolated or embedded in an environment. Only when the atom that has been ionized and excited by the primary particle impact is in immediate proximity of another atom is a fragmentation route known as interatomic Coulombic decay (ICD) enabled. This leads to the emission of a low energy electron. Over the past decade ICD was explored in several experiments following photoionization. Most recent results show its observation even in water clusters. Here we show the quantitative role of ICD for the production of low energy electrons by ion impact, thus approaching a scenario closer to that of radiation damage by alpha particles: We choose ion energies on the maximum of the Bragg peak where energy is most efficiently deposited in tissue. We compare the electron production after colliding He(+) ions on isolated Ne atoms and on Ne dimers (Ne(2)). In the latter case the Ne atom impacted is surrounded by a most simple environment already opening ICD as a deexcitation channel. As a consequence, we find a dramatically enhanced low energy electron yield. The results suggest that ICD may have a significant influence on cell survival after exposure to ionizing radiation.

Intermolecular coulomb decay at weakly coupled heterogeneous interfaces

Surface ejection of H(+)(H(2)O)(n=1-8) from low energy electron irradiated water clusters adsorbed on graphite and graphite with overlayers of Ar, Kr or Xe results from intermolecular Coulomb decay (ICD) at the mixed interface. Inner valence holes in water (2a(1)(-1)), Ar (3s(-1)), Kr (4s(-1)), and Xe (5s(-1)) correlate with the cluster appearance thresholds and initiate ICD. Proton transfer occurs during or immediately after ICD and the resultant Coulomb explosion leads to H(+)(H(2)O)(n=1-8) desorption with kinetic energies that vary with initiating state, final state, and interatomic or molecular distances.

Interatomic electronic decay driven by nuclear motion

The interatomic electronic decay after inner-valence ionization of a neon atom by a single photon in a neon-helium dimer is investigated. The excited neon atom relaxes via interatomic Coulombic decay and the excess energy is transferred to the helium atom and ionizes it. We show that the decay process is only possible if the dimer's bond stretches up to 6.2 Å, i.e., to more than twice the equilibrium interatomic distance of the neutral dimer. Thus, it is demonstrated that the electronic decay, taking place at such long distances, is driven by the nuclear motion.

Interatomic Coulombic decay following photoionization of the helium dimer: observation of vibrational structure

Using synchrotron radiation we simultaneously ionize and excite one helium atom of a helium dimer (He2) in a shakeup process. The populated states of the dimer ion [i.e., He(*+)(n = 2, 3) - He] are found to deexcite via interatomic Coulombic decay. This leads to the emission of a second electron from the neutral site and a subsequent Coulomb explosion. In this Letter we present a measurement of the momenta of fragments that are created during this reaction. The electron energy distribution and the kinetic energy release of the two He+ ions show pronounced oscillations which we attribute to the structure of the vibrational wave function of the dimer ion.

Photo- and auger-electron recoil induced dynamics of interatomic Coulombic decay

At photon energies near the Ne K edge it is shown that for 1s ionization the Auger electron, and for 2s ionization the fast photoelectron, launch vibrational wave packets in a Ne dimer. These wave packets then decay by emission of a slow electron via interatomic Coulombic decay (ICD). The measured and computed ICD electron spectra are shown to be significantly modified by the recoil induced nuclear motion.

Experimental separation of virtual photon exchange and electron transfer in interatomic coulombic decay of neon dimers

We investigate the interatomic Coulombic decay (ICD) of neon dimers following photoionization with simultaneous excitation of the ionized atom (shakeup) in a multiparticle coincidence experiment. We find that, depending on the parity of the excited state, which determines whether ICD takes place via virtual dipole photon emission or overlap of the wave functions, the decay happens at different internuclear distances, illustrating that nuclear dynamics heavily influence the electronic decay in the neon dimer.

Ab initio lifetimes in the interatomic Coulombic decay of neon clusters computed with propagators

Interatomic Coulombic decay (ICD) is a radiationless decay mechanism occurring via electron emission in an inner-valence ionized weakly bound cluster. The ICD has been studied for the neon clusters Nen (n=2,...,5). The decay widths of the neon clusters are calculated using ab initio Green's function method. The non-Dyson version of Green's function is employed. This propagator is analytically continued into the complex energy plane with the aid of a complex absorbing potential, and the decaying states are found as resonance states in this plane.

Observation of resonant interatomic coulombic decay in Ne clusters

We have measured the electron spectra of Ne clusters after excitation with photon energies around the 2s inner valence threshold. At two photon energies below threshold, a resonantly enhanced surplus of low kinetic-energy electrons is observed. The kinetic energy of the peak does not vary with the photon energy and is slightly larger than the transition energy of Interatomic Coulombic Decay (ICD) above threshold. This leads us to assume that an ICD-like process is present. In analogy to the Auger and the resonant Auger decay this new phenomenon is termed resonant ICD.

Properties of resonant interatomic Coulombic decay in Ne dimers

Properties of the interatomic Coulombic decay (ICD) process in Ne dimers have been obtained by tracking the formation of energetic Ne+ ions. The double photoionization cross section, deduced from the Ne+/Ne+ coincidence signal, is dominated by the ICD process and presents a threshold 280 meV below the atomic Ne+2s(-1) threshold. Rydberg excitation of a 2s electron in the dimer creates molecular Rydberg states whose Sigma and Pi symmetries have been resolved. These excited states decay by a resonant ICD process releasing an energetic Ne+ ion and a neutral excited Ne* fragment. Subsequent autoionization of the Ne* fragment explains a double photoionization threshold below the dimer 2s ionization threshold.