# Heavy N⁺ ion transfer in doubly charged N_2Ar van der Waals cluster

Xiaolong Zhu, Xiaoqing Hu, Shuncheng Yan, Yigeng Peng, Wentian Feng, Dalong Guo, Yong Gao, Shaofeng Zhang, Amine Cassimi, Jiawei Xu, et al.

► To cite this version:

Xiaolong Zhu, Xiaoqing Hu, Shuncheng Yan, Yigeng Peng, Wentian Feng, et al. Heavy N⁺ ion transfer in doubly charged N_2Ar van der Waals cluster. Nature Communications, 2020, 11 (1), 10.1038/s41467-020-16749-w. hal-03346102

HAL Id: hal-03346102

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

Submitted on 19 Sep 2024

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>

ARTICLE

https://doi.org/10.1038/s41467-020-16749-w

OPEN

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

Check for updates

# Heavy N⁺ ion transfer in doubly charged N₂Ar van der Waals cluster

XiaoLong Zhu<sup>1,2,6</sup>, XiaoQing Hu<sup>3,6</sup>, ShunCheng Yan<sup>1,2</sup>, YiGeng Peng<sup>3</sup>, WenTian Feng<sup>1</sup>, DaLong Guo<sup>1,2</sup>, Yong Gao<sup>1</sup>, ShaoFeng Zhang<sup>1,2</sup>, Amine Cassimi<sup>4</sup>, JiaWei Xu<sup>1,2</sup>, DongMei Zhao<sup>1</sup>, DaPu Dong<sup>1,2</sup>, Bang Hai<sup>1,2</sup>, Yong Wu<sup>3,5✉</sup>, JianGuo Wang<sup>3</sup> & X. Ma<sup>1,2✉</sup>

Van der Waals clusters are weakly bound atomic/molecular systems and are an important medium for understanding micro-environmental chemical phenomena in bio-systems. The presence of neighboring atoms may open channels otherwise forbidden in isolated atoms/molecules. In hydrogen-bond clusters, proton transfer plays a crucial role, which involves mass and charge migration over large distances within the cluster and results in its fragmentation. Here we report an exotic transfer channel involving a heavy N⁺ ion observed in a doubly charged cluster produced by 1 MeV Ne<sup>8+</sup> ions: (N<sub>2</sub>Ar)<sup>2+</sup>→N⁺+NAr⁺. The neighboring Ar atom decreases the N<sub>2</sub><sup>2+</sup> barrier height and width, resulting in significant shorter lifetimes of the metastable molecular ion state N<sub>2</sub><sup>2+</sup>(X<sup>1</sup>Σ<sup>+</sup><sub>g</sub>). Consequently, the breakup of the covalent N⁺–N⁺ bond, the tunneling out of the N⁺ ion from the N<sub>2</sub><sup>2+</sup> potential well, as well as the formation of an N⁻Ar⁺ bound system take place almost simultaneously, resulting in a Coulomb explosion of N⁺ and NAr⁺ ion pairs.

<footnote n="1">¹Institute of Modern Physics, Chinese Academy of Sciences, 730000 Lanzhou, China. ²University of Chinese Academy of Sciences, 100049 Beijing, China.</footnote>

<footnote n="3">³Institute of Applied Physics and Computational Mathematics, 100088 Beijing, China. ⁴CIMAP, CEA/CNRS/ENSICAEN/UNICAEN, BP5133, 14070 Caen, France. ⁵HEDPS, Center of Applied Physics and Technology, Peking University, 100871 Beijing, China. ⁶These authors contributed equally: XiaoLong Zhu, XiaoQing Hu.✉email: wu_yong@iapcm.ac.cn; x.ma@impcas.ac.cn</footnote>
Van der Waals (vdW) interactions are ubiquitous in nature and important for many physical and chemical phenomena. A good understanding of intermolecular interaction is essential since weakly bond systems widely exist in biology and play a key role as a reactive intermediate in many processes. Weak bonds may be easily broken but they are very important because they enable to determine and stabilize the shape of biological molecules, such as stabilizing the secondary structure (e.g. alpha-helix and beta-pleated sheet) of proteins¹. Dimers are good prototypes to investigate this weakly bound interaction. Indeed, atomic/molecular dimers have a typical weak binding energy in the order of meV and have a large distance between its two components usually above 7 a.u., which exceeds the covalent bond in molecules (~2 a.u.). Although weakly bound rare gas dimers have been studied extensively during the last two decades²⁻¹⁵, studies of dimers consisting of one or two diatomic molecules remain rare, some examples can be found in refs. ¹⁶⁻²². Therefore, it is essential to understand the structure and fragmentation/dissociation dynamics of charged molecular dimers at the level of molecular motion².

A prominent example is the interatomic Coulombic decay (ICD) theoretically predicted in 1997 (ref. ³), where electronic excitation energy is efficiently transferred over a large internuclear distance within a rare gas dimer to allow for fast nonradiative relaxation. In 2003, ICD was first observed experimentally in $N_2$ rare gas cluster following inner-valence electron ionization⁵. Later, numerous experiments, where photons, electrons, and heavy ions were used as projectiles, have been performed to confirm the ICD and new dissociation mechanisms have been found in different rare gas dimers or larger clusters⁶⁻²¹,²³⁻²⁸, such as radiative charge transfer⁸,¹⁵ or electron-transfer-mediated decay¹⁰,²⁷. For this kind of processes, the energy and charge transfer are mediated by virtual photon/Coulomb interaction, or electron transfer. In addition to the mechanisms mentioned above, an intermolecular proton transfer was reported in the molecular dimers, such as acetylene dimer²⁹, water cluster²³,³⁰, and ammonia clusters³¹. In this process, a proton is transferred from a parent molecule to its partner molecule in a dimer which is bound together by hydrogen bond, and then, leading to a

<figure><img src="image_4.png" /><figcaption>**Fig. 1** Two-dimensional TOF correlation map. The coincidence map shows the time-of-flights between the first and the second fragment ions arriving at the detector. The product ion species are labeled in the plot for horizontal and vertical axes according to their TOFs, respectively. The corresponding reactions labeled with a number are explained in the text.</figcaption></figure>

dissociation breakup. Since the discovery of the double helix form of DNA and the hypothesis of DNA mutation induced by proton transfer³² more than 50 years ago, it has been recognized that proton transfer is crucial in many chemical and biological processes³³. For example, a proton transfer reaction may produce a rare tautomer, which is regarded as the source of spontaneous mutations, and might lead to some important biochemistry processes related to other relevant diseases such as cellular aging and cancer³²,³⁴. On the other hand, there are excited-state proton transfer and multiple proton transfer processes which also play vital roles in biological systems. Nowadays, at the molecular level the simulations employing quantum chemistry could help us to have a thorough and detailed understanding of the mechanism of chemical and biological processes relevant to DNA, heredity, mutation, cancer, and green fluorescence protein, etc.³⁵.

In this article, we report a dissociation breakup mechanism of $N_2^{2+}$Ar cluster induced by heavy $N^+$ ion transfer: $N_2^{2+}$Ar → $N^+$ + $NAr^+$. Here, the strong covalent bond N-N breaks up and a new N-Ar bond is formed, due to $N^+$ transfer from $N_2^{2+}$ center to Ar in doubly charged $N_2^{2+}$Ar. To the best of our knowledge, such a heavy ion transfer process in vdW cluster has never been reported before and the consequent formation of $NAr^+$ is a novel scenario. It is intriguing to build bridges between the heavy ion transfer processes and biological processes at molecular level in the micro-environment of biological system as well as, for example, in the understanding of micro-mechanism of cancer therapy by heavy ion irradiations.

## Results

**Identification of heavy ion transfer channel.** In the present work, doubly charged cluster ($N_2$Ar)$^{2+}$ ions are produced by using 1 MeV $Ne^{8+}$ ions impact on neutral vdW cluster $N_2$Ar target, and the fragment ions are detected in coincidence with the scattered ions. The information of time-of-flight (TOF) and positions of each charged fragment arriving at the detector is used to reconstruct their three-dimensional initial momentum vectors (see Methods). Figure 1 shows the two-dimensional correlation map between the TOFs of the first and the second fragment ions produced in $Ne^{8+}$-$N_2$Ar collisions with one or two electrons stabilized on the projectile. At the top of horizontal axis the ion species and charge states corresponding to the TOF1 are labeled. Similarly, the ion species and charge states corresponding to the TOF2 (vertical axis) are also indicated. Usually, the different reaction channels can be identified according to the charge conservation before and after the collision. As a general rule for TOF correlation map, the slash lines are characteristics of two-body dissociation processes, while the extended area around the slash lines originates from the three-body or multi-body fragmentation³⁶. Taken as an example, the fragment products as labeled in Fig. 1 can be identified mainly from the following dissociation channels of the doubly charged dimers:

$$Ar_2^{2+} \rightarrow Ar^+ + Ar^+ \quad (1)$$

$$(N_2Ar)^{2+} \rightarrow N_2^+ + Ar^+ \quad (2)$$

$$(N_2)_2^{2+} \rightarrow N_2^+ + N_2^+ \quad (3)$$

However, there is one two-body dissociation channel labeled (4) in Fig. 1 which cannot be simply identified as a direct fragmentation of the vdW dimer $N_2$Ar according to the above mentioned rules. From calibration we find that one fragment has a charge to mass ratio of 1/14, while the correlated fragment has a charge to mass ratio of 1/54. It is surprising, because none of the dimer components and probable residual gaseous atoms have such a combination of charge to mass ratios. The coincidence
<figure><img src="image_5.png" /><figcaption>**Fig. 2** Measured KER distributions of N⁺ ion transfer channel. The blue solid circles are experimental results. The error bars are standard statistical uncertainties. The red curve is a fit to the experimental data, and its peak locates at 4.9 eV. Source data are provided as a Source Data file.</figcaption></figure>

measurement between (1/14)⁺ and (1/54)⁺, charge to mass ratio fragments, suggests that this channel corresponds to the following reaction:

$$ (N_2Ar)^{2+} \rightarrow N^+ + NAr^+ \quad (4) $$

In this reaction channel, a significant discovery is that a heavy N⁺ ion (mass = 14) is transferred from one molecule center (N₂⁺) to the other (Ar) in the dimer. One has to note that this exotic dissociation channel means that after the N₂Ar dimer loses two electrons, the strong N-N covalent bond breaks up and, due to the presence of the Ar atom, a new N-Ar bond is formed. We refer to this channel as N⁺ ion transfer channel. From the potential energy curves (PECs) of NAr⁺ (ref. ³⁷), the stable molecular ion NAr⁺ can only be formed from N⁺ ion and neutral Ar atom. This suggests that the channel (4) can only originate from the N₂⁺ Ar parent dimer ion.

**Kinetic energy release of N⁺ ion transfer channel.** The kinetic energy release (KER) distribution of N⁺ ion transfer channel shows a single peak at 4.9 eV as illustrated in Fig. 2. In the reaction (4), N⁺ ion transits from an initially bound N₂⁺ state to a final NAr⁺ bound state. On the other hand, the KER results from three contributions: the first one is the polarization energy of N₂⁺ - Ar bond, which is very weak and can be ignored; the second one is the dissociation energy of N₂⁺, which is determined from the initial electronic metastable state of N₂⁺, and the last one is the binding energy of NAr⁺. According to the PEC of N₂⁺ molecular ion ³⁸, the N₂⁺ molecular ion of N₂Ar dimer could be populated in the low-lying electronic states, namely, X¹Σ⁺g, a³Πu, D³Σ⁺u, A¹Πu, and A³Σ⁻g, etc. However, from the measured KER = 4.9 eV, we infer that the N⁺ ion transfer channel originates essentially from the relaxation of low vibrational levels of N₂⁺(X¹Σ⁺g)Ar and N₂⁺(a³Πu)Ar states, because their dissociation energy is less than 4.9 eV. Here, the neutral Ar in the dimer plays a crucial role in the dissociation process of N₂⁺(X¹Σ⁺g, a³Πu) dication.

Discussion

To understand the reaction mechanism of N⁺ ion transfer channel, let us first consider the PECs of the doubly charged cluster (N₂Ar)²⁺. We mainly focus on the relaxation of

<figure><img src="image_6.png" /><figcaption>**Fig. 3** Calculated potential energy curves for doubly charged cluster N₂⁺Ar ions. The PECs for T-shaped and linear-shaped N₂⁺Ar ions along the N₂-Ar bond shown with solid black curve and dashed red curve, respectively. The dimer geometries for T-shaped and linear-shaped structure are schematically displayed.</figcaption></figure>

N₂⁺(X¹Σ⁺g)Ar, which should be more important than the relaxation of N₂⁺(a³Πu)Ar due to the fact that the N₂⁺(X¹Σ⁺g) state is predominantly located at the dissociation energy region allowed by Frank-Condon transition. The PECs for N₂⁺(X¹Σ⁺g)Ar along the N₂-Ar dimer axis are calculated as shown in Fig. 3, see calculation of PECs in Methods. A comparison of the PECs for the linear-shaped and T-shaped N₂⁺(X¹Σ⁺g)Ar dimer is displayed in Fig. 3. It is clearly shown that the energy minimum for linear-shaped N₂⁺Ar locates at R<sub>N₂-Ar</sub> = 7.2 a.u. and is about 0.2 eV lower than the one for T-shaped, (here, R<sub>X-Y</sub> defines the distance between the Y atom and the center of mass of the X molecule in the whole context). It is well known that neutral N₂Ar cluster has a T-shaped structure¹⁶⁻¹⁸. When the Ne⁸⁺ projectile captures two electrons from N₂-site of neutral N₂Ar cluster, N₂⁺Ar ion should undergo firstly an isomerization process going from an initial T-shape to a linear-shape, and then vibrate in the vicinity of the minimum before dissociation.

With the help of theoretical PECs, we present an intuitive physical picture of the heavy N⁺ ion transfer reaction. Different from the PEC of the monomer N₂⁺ (black solid curve in Fig. 4), the PEC of linear-shaped N₂⁺Ar along the N⁺-N⁺ bond is exhibited by a red dashed curve in Fig. 4, where the influence of neutral Ar atom on the potential energy of N₂⁺ is considered. It is clearly shown that a new potential well, due to the presence of Ar atom, can be found as the R<sub>N-N</sub> (upper horizontal axis in Fig. 4) increases, and its minimum locates at R<sub>N-N</sub> = 7.8 a.u., corresponding to the distance R<sub>N-Ar</sub> = 7.2-7.8/2 = 3.3 a.u. This clearly suggests that a stable NAr⁺ molecular ion could be formed.

However, to form the NAr⁺ ion, the N⁺ ion, initially bound to N₂⁺ molecule, has to overcome the barrier E<sub>c</sub> (i.e. the breakup of N₂⁺ covalent bond) and to be trapped in the NAr⁺ well, which is of course extremely difficult. The only possible pathway is tunneling, namely, the N⁺ ion from the left potential well of N₂⁺(X¹Σ⁺g)Ar (in vibration states v = a) tunnels into the right potential well of N⁺+NAr⁺ (the pink dashed arrow in Fig. 4), and the NAr⁺ is formed in high vibrational levels (v = b) (see the green wave line). Then, the Coulomb explosion between N⁺ and NAr⁺ ion pair takes place along the blue solid curve as shown in Fig. 4. Due to the strong couplings between Coulomb repulsion and vibration of NAr⁺ at close distances, the NAr⁺ at high vibration levels (v = b) can be converted to lower vibration states
<figure><img src="image_7.png" /><figcaption>**Fig. 4** Schematic diagram of heavy ion transfer process through tunneling.</figcaption></figure>

The PECs of N<sub>2</sub><sup>2+</sup>(X<sup>1</sup>Σ<sub>g</sub>)Ar with/without the presence of neutral Ar along R<sub>N-N</sub> (upper horizontal axis) and N<sup>+</sup>-NAr<sup>+</sup> along R<sub>N-NAr</sub> (bottom horizontal axis). The black solid curve for N<sub>2</sub><sup>2+</sup>(X<sup>1</sup>Σ<sub>g</sub><sup>+</sup>) without the neutral Ar, the red dashed curve for N<sub>2</sub><sup>2+</sup>(X<sup>1</sup>Σ<sub>g</sub><sup>+</sup>)Ar at fixed distance

R<sub>N<sub>2</sub>-Ar</sub> = 7.2 a.u. and the blue solid curve indicates the PEC of N<sup>+</sup>-NAr<sup>+</sup> with R<sub>N-NAr</sub> = 3.3 a.u. The gray areas indicate the possibly populated initial and final vibrational states of N<sub>2</sub><sup>2+</sup>Ar and NAr<sup>+</sup> ions, respectively, which contribute to the heavy ion transfer channel. The final NAr<sup>+</sup> potential well is plotted on the left. The pink dashed arrow indicates the N<sup>+</sup> ion tunneling. The green wave line indicates the intermediate state of ion NAr<sup>+</sup>. The origin arrow indicates the fragmentation of N<sup>+</sup> and NAr<sup>+</sup> ion pair via Coulomb explosion. The origin of the potential energy curves is taken such that the potential energy of N<sup>+</sup>, N<sup>+</sup>, and neutral Ar at large separations is zero. The light blue ball represents nitrogen atom while the black ball represents argon atom.

(v = c). The correlation between the KER of N⁺ ion transfer channel and the initial states can be estimated. The initial energy of the doubly charged cluster N₂⁺(X¹Σ⁺g)Ar consists of a sum of minimum potential energy Eᵢ and vibration energy Eₐ (v = a). Similarly, the final potential energy of the system is given by adding the minimum potential energy Eᵦ and the vibration energy Eₑ (v = c). On the other hand, as the energy of the system before and after the heavy ion transfer is conserved, as a result, the KER could be obtained as follows: KER = (Eᵢ+Eₐ) - (Eᵦ + Eₑ). Namely, 4.9 = (4.0 + Eₐ) - (-2.05 + Eₑ), that is, Eₑ - Eₐ = 1.15 eV. Because the final maximum vibration energy is 2.05 eV, while the initial minimum vibration energy is 0 eV, the only possible vibration energies for initial and final states are restricted in the range of 0.0-0.9 and 1.15-2.05 eV, respectively, as shown by the gray areas in Fig. 4.

One has to note that, the monomer N<sub>2</sub><sup>2+</sup> is a metastable ion which has very long lifetimes. The lifetime of N<sub>2</sub><sup>2+</sup>(X<sup>1</sup>Σ<sub>g</sub><sup>+</sup>) state for its vibrational level v = 8 is as long as 312.5 days (2.7 × 10<sup>6</sup> s)<sup>39</sup>, and is even longer for v < 8 states. However, for the dimer N<sub>2</sub><sup>2+</sup>(X<sup>1</sup>Σ<sub>g</sub><sup>+</sup>)Ar, the nearby Ar atom modifies the monomer's potential and leads to significant changes of the lifetimes. We calculated the tunneling lifetimes of N<sub>2</sub><sup>2+</sup>(X<sup>1</sup>Σ<sub>g</sub><sup>+</sup>) in the presence of Ar atom with R<sub>N<sub>2</sub>-Ar</sub> = 7.2 a.u. using the complex absorption potential<sup>40</sup>, see lifetime estimation in Methods. The calculated tunneling lifetimes of N<sub>2</sub><sup>2+</sup>(X<sup>1</sup>Σ<sub>g</sub><sup>+</sup>) (v = 4, 5, 6) are 1110, 245, and 69 ns (1 ns = 10<sup>-9</sup> s), respectively, which means its lifetime can be reduced by more than 12 orders of magnitude in the doubly charged dimer N<sub>2</sub><sup>2+</sup>(X<sup>1</sup>Σ<sub>g</sub><sup>+</sup>)Ar, from hundred days to nanoseconds. This strongly indicates that, because of the

neighboring Ar atom, the originally quasi-stable monomer ions N<sub>2</sub><sup>2+</sup> can dissociate rapidly by tunneling through the potential barrier, form the NAr<sup>+</sup> ion and, hence, create the N<sup>+</sup>-NAr<sup>+</sup> ion pair.

The heavy ion transfer scenario in a charged dimer may be general for molecular dimer ions such as AB<sup>2+</sup>X as long as two conditions are met: (i) the metastable covalent bond A-B is toward the neutral particle X, which is the case for most diatomic molecular ions and also for other polar covalent bonds, such as C-H bond; (ii) the neutral particle X and the neighboring ion A<sup>+</sup>/B<sup>+</sup> can form a bound system XA<sup>+</sup>/XB<sup>+</sup>. This kind of ion transfer process may have potential importance in understanding the underlying mechanisms of biochemistry in tissues, and play an important role in the irradiation damage caused by ion or electron radiation, and furthermore, may directly affect the expression of DNA and protein in vivo.

In conclusion, we present a channel of heavy N⁺ transfer and NAr⁺ ion formation in the two-body dissociation of (N₂Ar)₂⁺ dimer ion, which results from the breakup of linear-shaped dimer ion N₂⁺(X¹Σ⁺g, a³Π_u)Ar. This transfer mechanism is explained by the tunneling of the heavy N⁺ ion and, hence, leads to the formation of NAr⁺ ion. The dissociation of N₂⁺(X¹Σ⁺g) states through the tunneling of N⁺ becomes possible due to the presence of Ar atom in the dimer. Consequently, the N⁺-N⁺ bond breaks up quickly, and the tunneling of N⁺ ion and the formation of N⁺-Ar bond take place almost simultaneously. This is viewed as bound-bound states transition of the heavy ion of N⁺, followed by Coulomb explosion between N⁺ and NAr⁺ ion pairs. Such a mechanism of heavy ion transfer may be general for any molecular dimer ions of type AB<sup>2+</sup>X, and may be of potential importance in various biochemical processes.

## Methods

**Experimental approaches.** Our experiment was conducted at 320 kV platform for multidisciplinary research with highly charged ions in the Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou<sup>41</sup>. A standard Cold Target Recoil Ion Momentum Spectroscopy/Reaction Microscope was employed to determine the momentum of each fragments in coincidence<sup>41-43</sup>. The highly charged Ne<sup>8+</sup> ions were produced by the permanent magnet electron cyclotron resonance ion source and were accelerated to an energy of 1 MeV. Then the Ne<sup>8+</sup> ions were delivered to the reaction chamber, where it crossed a supersonic gas target. Various dimers were generated in the molecular jet, such as N₂Ar, (N₂)₂, and Ar₂, by expanding a mixture gas of N₂ and Ar through a 30-µm nozzle at 20 bars stagnation pressure under room temperature. The gas pressure ratio in the mixture was N₂:Ar ~ 1:1. Charged fragments produced in the collisions were guided by 173 V/cm uniform electric field to a position-sensitive multi-hit detector with a position resolution of 0.1 mm. The information of positions and the TOF of each charged fragment arriving at the detector are obtained; consequently, their three-dimensional initial momentum vectors can be reconstructed by off-line analysis. Downstream, the projectile ions were collected by a Faraday-Cup and the scattered Ne<sup>7+</sup> and Ne<sup>6+</sup> ions were detected by a delay-line anode detector. The fragment momentum calibration is carried out in the same conditions by the two-body dissociation of monomer N₂⁺ ions, which has clearly separated and known peaks<sup>38</sup>.

**Calculations of PECs.** The PECs calculations for N<sub>2</sub><sup>2+</sup>(X<sup>1</sup>Σ<sub>g</sub><sup>+</sup>)Ar along the N<sub>2</sub>-Ar dimer axis are performed using the multi-reference double-excited configuration interaction<sup>44</sup> method with molecular orbitals optimized by complete active space self-consistent field method<sup>45-47</sup>. Here the active space consists of all valence orbitals and totally about 1000 reference configuration state functions are considered in configuration interaction calculations. The orbital wave functions are expanded on aug-cc-pVtZ basis set. The same method is used to calculate the PECs of N<sub>2</sub><sup>2+</sup>(X<sup>1</sup>Σ<sub>g</sub><sup>+</sup>) and NAr<sup>+</sup>. Then, the PEC of linear-shaped N<sub>2</sub><sup>2+</sup>(X<sup>1</sup>Σ<sub>g</sub><sup>+</sup>)Ar along the N-N bond is calculated using the two-body potential model<sup>48</sup>, where only the two-body correlation interactions are considered based on the PECs of N<sub>2</sub><sup>2+</sup> and NAr<sup>+</sup>, and the three-body correlation interaction is neglected. Such a model is reasonable for the PECs calculation of linear-shaped N<sub>2</sub><sup>2+</sup>(X<sup>1</sup>Σ<sub>g</sub><sup>+</sup>)Ar since one of the N<sup>+</sup> ions is far away from Ar (>8 a.u.) and the three-body correlation interaction is very weak.
**Lifetime estimation.** Based on the PECs obtained, the lifetimes of N₂⁺(X¹Σ⁺₆)Ar are calculated using the complex absorption potential (CAP) method⁴⁹,⁵⁰. In the CAP calculations, the nuclear Hamiltonian H(η) is extended by a small absorbing potential

$$H(\eta) = -iW(R) + H_0 \qquad (5)$$

and

$$W(R) = \begin{cases} 0, & R<R_0 \\ (R-R_0)^2, & R>R_0 \end{cases} \qquad (6)$$

in which η is a variation parameter and R₀ = 10 a.u. is used in the present calculation. A set of η values are selected to solve the eigenvalue equation

$$H(\eta)\psi(R) = E(\eta)\psi(R) \qquad (7)$$

and the eigenvalue of resonant state can be obtained with the condition of

$$\eta \left| \frac{\partial E(\eta)}{\partial \eta} \right| = \min. \qquad (8)$$

Finally, the eigenvalue we obtain is as follows:

$$E = E_0 - \frac{1}{2i\Gamma}, \qquad (9)$$

where E₀ is the resonant energy, Γ is the resonant width or decay rate, and τ = 1/Γ is the lifetime.

## Data availability

Source data are provided with this paper. The data supporting this study are also available from the corresponding author upon reasonable request.

## Code availability

The code that supports the theoretical plots within this paper is available from the corresponding authors upon reasonable request.

Received: 27 March 2020; Accepted: 20 May 2020;

Published online: 12 June 2020

## References

1. Campbell, N. & Reece, J. Biology (Benjamin Cummings, 2005).

2. Linnartz, H., Verdes, D. & Maier, J. P. Rotationally resolved infrared spectrum of the charge transfer complex [Ar-N₂]⁺. *Science* **297**, 1166 (2002).

3. Cederbaum, L. S., Zobeley, J. & Tarantelli, F. Giant intermolecular decay and fragmentation of clusters. *Phys. Rev. Lett.* **79**, 4778 (1997).

4. Zobeley, J., Santra, R. & Cederbaum, L. S. Electronic decay in weakly bound heteroclusters: energy transfer versus electron transfer. *J. Chem. Phys.* **115**, 5076 (2001).

5. Marburger et al. Experimental evidence for interatomic Coulombic decay in Ne clusters. *Phys. Rev. Lett.* **90**, 203401 (2003).

6. Jahnke, T. et al. Experimental observation of interatomic Coulombic decay in neon dimers. *Phys. Rev. Lett.* **93**, 163401 (2004).

7. Öhrwall, G. et al. Femtosecond interatomic Coulombic decay in free neon clusters: large lifetime differences between surface and bulk. *Phys. Rev. Lett.* **93**, 173401 (2004).

8. Matsumoto, J. et al. Asymmetry in multiple-electron capture revealed by radiative charge transfer in ar dimers. *Phys. Rev. Lett.* **105**, 263202 (2010).

9. Titze, J. et al. Ionization dynamics of helium dimers in fast collisions with He⁺⁺. *Phys. Rev. Lett.* **106**, 033201 (2011).

10. Sakai, K. et al. Electron-transfer-mediated decay and interatomic Coulombic decay from the triply ionized states in argon dimers. *Phys. Rev. Lett.* **106**, 033401 (2011).

11. O'Keeffe, P. et al. The role of the partner atom and resonant excitation energy in interatomic coulombic decay in rare gas dimers. *J. Phys. Chem. Lett.* **4**, 1797 (2013).

12. Iskandar, W. et al. Atomic site-sensitive processes in low energy ion-dimer collisions. *Phys. Rev. Lett.* **113**, 143201 (2014).

13. Fruhling, U. et al. Time-resolved studies of interatomic Coulombic decay. *J. Electron. Spectrosc. Relat. Phenom.* **204**, 237 (2015).

14. Iskandar, W. Interatomic Coulombic decay as a new source of low energy electrons in slow ion-dimer collisions. *Phys. Rev. Lett.* **114**, 033201 (2015).

15. Ren Xueguang et al. Direct evidence of two interatomic relaxation mechanisms in argon dimers ionized by electron impact. *Nat. Commun.* **7**, 11093 (2016).

16. Wu, J. et al. Structures of N₂Ar, O₂Ar, and O₂Xe dimers studied by Coulomb explosion imaging. *J. Chem. Phys.* **137**, 104308 (2012).

17. Wu, J. et al. Strong field multiple ionization as a route to electron dynamics in a van der Waals cluster. *Phys. Rev. Lett.* **111**, 083003 (2013).

18. Wu Chengyin et al. Communication: Determining the structure of the N₂Ar van der Waals complex with laser-based channel-selected Coulomb explosion. *J. Chem. Phys.* **140**, 141101 (2014).

19. Trinter, F. et al. Resonant Auger decay driving intermolecular Coulombic decay in molecular dimers. *Nature* **505**, 664 (2014).

20. Ding Xiaoyan et al. Ultrafast dissociation of Metastable CO²⁺ in a dimer. *Phys. Rev. Lett.* **118**, 153001 (2017).

21. Mery, A. et al. Role of a neighbor ion in the fragmentation dynamics of covalent molecules. *Phys. Rev. Lett.* **118**, 233402 (2017).

22. Chiang, Ying-Chih et al. Molecular-bond breaking induced by interatomic decay processes. *Phys. Rev. A* **100**, 052701 (2019).

23. Jahnke, T. et al. Ultrafast energy transfer between water molecules. *Nat. Phys.* **6**, 139 (2010).

24. Zettergren, H. et al. Formations of dumbbell C₁₁₈ and C₁₁₉ inside clusters of C₆₀ molecules by collision with α particles. *Phys. Rev. Lett.* **110**, 185501 (2013).

25. Delaunay, R. et al. Molecular growth inside of polycyclic aromatic hydrocarbon clusters induced by ion collisions. *J. Phys. Chem. Lett.* **6**, 1536 (2015).

26. Xueguang, Ren et al. Experimental evidence for ultrafast intermolecular relaxation processes in hydrated biomolecules. *Nat. Phys.* **14**, 1062 (2018).

27. Unger, Isaak et al. Observation of electron-transfer-mediated decay in aqueous solution. *Nat. Chem.* **9**, 708 (2017).

28. Yan, S. et al. Interatomic relaxation processes induced in neon dimers by electron-impact ionization. *Phys. Rev. A* **97**, 010701(R) (2018).

29. Xu, S. et al. Damaging intermolecular energy and proton transfer processes in alpha-particle-irradiated hydrogen-bonded systems. *Angew. Chem. Int. Ed.* **57**, 17023 (2018).

30. Richter, Clemens et al. Competition between proton transfer and intermolecular Coulombic decay in water. *Nat. Commun.* **9**, 4988 (2018).

31. Bart, Oostenrijk et al. The role of charge and proton transfer in fragmentation of hydrogen-bonded nanosystems: the breakup of ammonia clusters upon single photon multi-ionization. *Phys. Chem. Chem. Phys.* **20**, 932 (2018).

32. Löwdin, P. O. Proton tunneling in DNA and its biological implications. *Rev. Mod. Phys.* **35**, 724 (1963).

33. Chojnacki, Henryk Special issue on proton transfer processes. *Int. J. mol. Sci.* **4**, 408 (2003).

34. Jacquemin, D. et al. Assessing the importance of proton transfer reactions in DNA. *Acc. Chem. Res.* **47**, 2467 (2014).

35. Panwang Zhou & Han, Keli Unraveling the detailed mechanism of excited-state proton transfer. *Acc. Chem. Res.* **51**, 1681 (2018).

36. Zhu, X. L. et al. Probing the geometry of Ar₂N₂ cluster by 1 MeV Ne⁸⁺ ions. *Nucl. Instrum. Methods Phys. Res. B* **408**, 42 (2017).

37. Broström, L. et al. The visible photoabsorption spectrum and potential curves of ArN⁺. *J. Chem. Phys.* **94**, 2734 (1991).

38. Lundqvist, M. et al. Doppler-free kinetic energy release spectrum of N₂²⁺. *J. Phys. B* **29**, 1489 (1996).

39. Mathur, D. Long-lived, doubly charged diatomic and triatomic molecular ions. *J. Phys. B* **28**, 3415 (1995).

40. Riss, U. V. & Meyer, H.-D. Calculation of resonance energies and widths using the complex absorbing potential method. *J. Phys. B* **26**, 4503-4535 (1993).

41. Ma, X. et al. Electron emission from single-electron capture with simultaneous single-ionization reactions in 30-keV/u He²⁺-on-argon collisions. *Phys. Rev. A* **83**, 052707 (2011).

42. Dörner, R. et al. Cold target recoil ion momentum spectroscopy: a 'momentum microscope' to view atomic collision dynamics. *Phys. Rep.* **330**, 95 (2000).

43. Ullrich, J. et al. Recoil-ion and electron momentum spectroscopy: reaction-microscopes. *Rep. Prog. Phys.* **66**, 1463 (2003).

44. Krebs, S. & Buenker, R. J. A new table-direct configuration interaction method for the evaluation of Hamiltonian matrix elements in a basis of linear combinations of spin-adapted functions. *J. Chem. Phys.* **103**, 5613 (1995).

45. Knowles, P. J. et al. An efficient second-order MC SCF method for long configuration expansions. *Chem. Phys. Lett.* **115**, 259 (1985).

46. Werner, H. J. et al. Molpro: a general-purpose quantum chemistry program package. *WIREs Comput. Mol. Sci.* **2**, 242 (2012).

47. Werner, H. J. et al. MOLPRO, A Package of Ab Initio Programs. (Cardiff University, Cardiff, U.K. 2012).

48. Qiu, B. & Xiulin, R. Molecular dynamics simulations of lattice thermal conductivity of bismuth telluride using two-body interatomic potentials. *Phys. Rev. B* **80**, 165203 (2009).

49. Moiseyev, N. Quantum theory of resonances: calculating energies, widths and cross-sections by complex scaling. *Phys. Rep.* **302**, 212 (1998).

50. Lin, X. H. et al. Theoretical study of resonances formed in low-energy Li⁻⁺H collisions. *Chem. Phys.* **522**, 10 (2019).

## Acknowledgements

This work was supported by the National Key Research and Development Program of China (Grant Nos. 2017YFA0402300, 2017YFA0402400, and 2017YFA0403200), the Science Challenge Project (Grant No. TZ2016005), and the NSFC of China (Grant Nos. 11934004 and 11974358). The authors thank the staff of 320 kV platform for their technical support. The authors acknowledge the helpful discussions with Dr. Bennaceur Najjari and suggestions from Dr. Michael Schulz.
**Author contributions**

X.M. and X.L.Z. designed the experiment. X.L.Z., W.T.F., D.L.G., Y.G., S.F.Z., D.M.Z., J.W.X., D.P.D., and B.H. performed the experiment. X.L.Z. and S.C.Y. analyzed the experimental data. X.Q.H., Y.G.P., Y.W., and J.G.W. carried out the theoretical calculations and analysis. X.M., X.L.Z., X.Q.H., A.C., and Y.W. contributed to the interpretation of the data. The manuscript was written by X.M., X.L.Z., X.Q.H., and Y.W. All authors discussed and approved the final version of the manuscript.

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary information** is available for this paper at https://doi.org/10.1038/s41467-020-16749-w.

**Correspondence** and requests for materials should be addressed to Y.W. or X.M.

**Peer review information** *Nature Communications* thanks the anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.

**Reprints and permission 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.

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

**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 license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020