Neon Dimer Research Articles

The Role of Bond Functions in Describing Intermolecular Electron Correlation for Van der Waals Dimers: A Study of (CH4)2 and Ne2.

We present a study of the intermolecular interactions in van der Waals complexes of methane and neon dimers within the framework of the CCSD method. This approach was implemented and applied to calculate and examine the behavior of the contracted two-particle reduced density matrix (2-RDM). It was demonstrated that the region near the minimum of the two-particle density matrix correlation part, corresponding to the primary bulk of the Coulomb hole contribution, exerts a significant influence on the dispersion interaction energetics of the studied systems. As a result, the bond functions approach was applied to improve the convergence performance for the intermolecular correlation energy results with respect to the size of the atomic basis. For this, substantial acceleration was achieved by introducing an auxiliary basis of bond functions centered on the minima of the 2-RDM. For both methane and neon dimers, this general conclusion was confirmed with a series of CCSD calculations for the 2-RDM and the correlation energies.

The challenge of ab initio calculations of small neon clusters.

Weakly bound neon dimer, trimer and tetramers are studied at HF and CCSD(T) levels using Dunning, ANO and SIGMA-s basis sets. Their ground-state binding energies are studied along with some structural properties. SIGMA-s basis sets have been developed explicitly for this issue but in a manner that can be readily applied to other atoms for the study of larger weakly bound systems. The difficulties for attaining accurate results on these systems are assessed by the computation of total, atomization and correlation energies, as well as equilibrium distances, with several basis sets of increasing size, ranging from non-augmented to double-augmented versions. Extrapolations are proposed to predict stabilization energies and the results are compared with previously published data.

Can the Fluxionality in Borospherene Influence the Confinement-Induced Bonding between Two Noble Gas Atoms?

A density functional theory study is performed to determine the stability and bonding in the neon dimer inside the B30N30 fullerene cage, the fluxional B40 cage, and within non-fluxional cages such as B12N12 and C60. The nature of bonding in the Ne2 encapsulated B40 is compared with the that in other cages in an attempt to determine whether any possible alterations are brought about by the dynamical nature of the host cage apart from the associated confinement effects. The bonding analysis includes the natural bond order (NBO), Bader's Atoms-in-Molecules electron density analysis (AIM), and energy decomposition analysis (EDA), revealing the non-covalent nature of the interactions between the Ne atoms and that between the Ne and the cage atoms. The formation of all the Ne2@cage systems is thermochemically unfavourable, the least being that for the B30N30 cage, which can easily be made favourable at lower temperatures. The Ne-Ne distance is lowest in the smallest cage and increases as the cage size increase due to steric relaxation experienced by the dimer. The dynamical picture of the systems is investigated by performing ab initio molecular dynamics simulations using the atom-centred density matrix propagation (ADMP) technique, which shows the nature of the movement of the dimer inside the cages, and by the fact that since it moves as a single entity, a weak bonding force holds them together, apart from their proven kinetic stability.

Indication of strong interatomic Coulombic decay in slow He2+ –Ne2 collisions

Electron removal in collisions of alpha particles with neon dimers is studied using an independent-atom-independent-electron model based on the semiclassical approximation of heavy-particle collision physics. The dimer is assumed to be frozen at its equilibrium bond length and collision events for the two ion–atom subsystems are combined in an impact parameter by impact parameter fashion for three mutually perpendicular orientations. Both frozen atomic target and dynamic response model calculations are carried out using the coupled-channel two-center basis generator method. We pay particular attention to inner-valence Ne(2s) electron removal, which is associated with interatomic Coulombic decay (ICD), resulting in low-energy electron emission and dimer fragmentation. Our calculations confirm a previous experimental result at 150 keV/amu impact energy regarding the relative strength of ICD compared to direct electron emission. They further indicate that ICD is the dominant Ne+ + Ne+ fragmentation process below 10 keV/amu, suggesting that a strong low-energy electron yield will be observed in the ion–dimer system in a regime in which the creation of continuum electrons is a rare event in the ion–atom problem.

Strengthening the case for interatomic Coulomb decay as a subdominant reaction channel in slow O3+−Ne2 collisions with independent-atom-model coupled-channel calculations

We report on electron removal calculations for 2.81 keV/amu Li$^{3+}$ and O$^{3+}$ ion collisions with neon dimers. The target is described as two independent neon atoms fixed at the dimer's equilibrium bond length, whose electrons are subjected to the time-dependent bare and screened Coulomb potentials of the classically moving Li$^{3+}$ and O$^{3+}$ projectile ions, respectively. Three mutually perpendicular orientations of the dimer with respect to the rectilinear projectile trajectories are considered and collision events for the two ion-atom subsystems are combined in an impact parameter by impact parameter fashion and are orientation-averaged to calculate probabilities and cross sections for the ion-dimer system. The coupled-channel two-center basis generator method is used to solve the ion-atom collision problems. We concentrate the ion-dimer analysis on one-electron and two-electron removal processes which can be associated with interatomic Coulomb decay, Coulomb explosion, and radiative charge transfer. We find that the calculated relative yields are in fair agreement with recent experimental data for O$^{3+}$-Ne$_2$ collisions if we represent the projectile by a screened Coulomb potential, but disagree markedly for a bare Coulomb potential, i.e., for Li$^{3+}$ impact. In particular, our calculations suggest that interatomic Coulomb decay is a significant reaction channel in the former case only, since capture of a Ne($2s$) electron to form hydrogenlike Li$^{2+}$ is unlikely.

Explicitly Correlated Dispersion and Exchange Dispersion Energies in Symmetry-Adapted Perturbation Theory.

The individual interaction energy terms in symmetry-adapted perturbation theory (SAPT) not only have different physical interpretations but also converge to their complete basis set (CBS) limit values at quite different rates. Dispersion energy is notoriously the slowest converging interaction energy contribution, and exchange dispersion energy, while smaller in absolute value, converges just as poorly in relative terms. To speed up the basis set convergence of the lowest-order SAPT dispersion and exchange dispersion energies, we borrow the techniques from explicitly correlated (F12) electronic structure theory and develop practical expressions for the closed-shell Edisp(20)-F12 and Eexch-disp(20)-F12 contributions. While the latter term has been derived and implemented for the first time, the former correction was recently proposed by Przybytek [ J. Chem. Theory Comput. 2018 , 14 , 5105 - 5117 ] using an Ansatz with a full optimization of the explicitly correlated amplitudes. In addition to reimplementing the fully optimized variant of Edisp(20)-F12, we propose three approximate Ansätze that substantially improve the scaling of the method and at the same time avoid the numerical instabilities of the unrestricted optimization. The performance of all four resulting flavors of Edisp(20)-F12 and Eexch-disp(20)-F12 is first tested on helium, neon, argon, water, and methane dimers, with orbital and auxiliary basis sets up to aug-cc-pV5Z and aug-cc-pV5Z-RI, respectively. The double- and triple-ζ basis set calculations are then extended to the entire A24 database of noncovalent interaction energies and compared with CBS estimates for Edisp(20) and Eexch-disp(20) computed using conventional SAPT with basis sets up to aug-cc-pV6Z with midbond functions. It is shown that the F12 treatment is highly successful in improving the basis set convergence of the SAPT terms, with the F12 calculations in an X-tuple ζ basis about as accurate as conventional calculations in bases with cardinal numbers (X + 2) for Edisp(20) and either (X + 1) or (X + 2) for Eexch-disp(20). While the full amplitude optimization affords the highest accuracy for both corrections, the much simpler and numerically stable optimized diagonal Ansatz is a very close second. We have also tested the performance of the simple F12 correction based on the second-order Møller-Plesset perturbation theory, SAPT-F12(MP2) [ Frey , J. A. ; Chem. Rev. 2016 , 116 , 5614 - 5641 ] and observed that it is also quite successful in speeding up the basis set convergence of conventional Edisp(20) + Eexch-disp(20), albeit with some outliers.

Numerical attoclock on atomic and molecular hydrogen

Numerical attoclock is a theoretical model of attosecond angular streaking driven by a very short, nearly a single oscillation, circularly polarized laser pulse. The reading of such an attoclock is readily obtained from a numerical solution of the time-dependent Schr\"odinger equation as well as a semi-classical trajectory simulation. By making comparison of the two approaches, we highlight the essential physics behind the attoclock measurements. In addition, we analyze the predictions of the Keldysh-Rutherford model of the attoclock [Phys. Rev. Lett. 121, 123201 (2018)]. In molecular hydrogen, we highlight a strong dependence of the width of the attoclock angular peak on the molecular orientation and attribute it to the two-center electron interference. This effect is further exemplified in the weakly bound neon dimer.

Dynamics of interatomic Coulombic decay in neon dimers by XUV-pump–XUV-probe spectroscopy

We apply the Monte Carlo wave-packet approach to study the interatomic Coulombic decay (ICD) dynamics of neon dimers after removing a $2s$ electron from one of the Ne atoms by one-photon absorption from an XUV pulse. This method reproduces well both the lifetime for the $2s$ inner-valence vacancy in $\mathrm{Ne}{}_{2}{}^{+}(2{s}^{\ensuremath{-}1})$ and the kinetic energy release (KER) spectra for the coincident ${\mathrm{Ne}}^{+}$ and ${\mathrm{Ne}}^{2+}$ fragments following triple ionization, i.e., two photoionizations and one ionization via ICD, of ${\mathrm{Ne}}_{2}$ measured in an XUV-pump--XUV-probe experiment [K. Schnorr et al., Phys. Rev. Lett. 111, 093402 (2013)]. Comparisons between the calculated and measured nuclear KER spectra give physical insights in the considered process. For example, an analysis of the ratios between the low- and high-energy peaks in the nuclear KER spectra for large delays provides an estimate of the photoionization cross sections for removing a $2p$ electron from the excited states in $\mathrm{Ne}{}_{2}{}^{+}(2{s}^{\ensuremath{-}1})$. Such comparisons also allow an estimate for the ICD rates for the $2s$ inner-valence vacancy in the single-site two hole state in $\mathrm{Ne}{{}_{2}}^{2+}(2{s}^{\ensuremath{-}1}2{p}^{\ensuremath{-}1})$. Finally, the influence of photon statistics of the free electron laser pulses on the nuclear KER spectra is considered.

Double-slit photoelectron interference in strong-field ionization of the neon dimer

Wave-particle duality is an inherent peculiarity of the quantum world. The double-slit experiment has been frequently used for understanding different aspects of this fundamental concept. The occurrence of interference rests on the lack of which-way information and on the absence of decoherence mechanisms, which could scramble the wave fronts. Here, we report on the observation of two-center interference in the molecular-frame photoelectron momentum distribution upon ionization of the neon dimer by a strong laser field. Postselection of ions, which are measured in coincidence with electrons, allows choosing the symmetry of the residual ion, leading to observation of both, gerade and ungerade, types of interference.

Interatomic relaxation processes induced in neon dimers by electron-impact ionization

We report an experimental observation of the interatomic Coulombic decay (ICD) and radiative charge-transfer (RCT) processes in a Ne dimer $(e,2e)$ following a 380-eV electron impact. By detecting the $\mathrm{N}{\mathrm{e}}^{+}\ensuremath{-}\mathrm{N}{\mathrm{e}}^{+}$ cation pair and one of the emitted electrons in coincidence, the fingerprint of the ICD process initiated by the inner-valence ionization of Ne is obtained. Furthermore, the experimental results and ab initio calculations together unambiguously confirm the occurrence of the RCT process, and we show that most of the low-energy electrons produced in ionization of the Ne dimers are due to the ICD, which strongly suggests the importance of the ICD in causing radiation damage in a biological medium.

Application of Matched-Filter Concepts to Unbiased Selection of Data in Pump-Probe Experiments with Free Electron Lasers

Pump-probe experiments are commonly used at Free Electron Lasers (FEL) to elucidate the femtosecond dynamics of atoms, molecules, clusters, liquids and solids. Maximizing the signal-to-noise ratio of the measurements is often a primary need of the experiment, and the aggregation of repeated, rapid, scans of the pump-probe delay is preferable to a single long-lasting scan. The limited availability of beamtime makes it impractical to repeat measurements indiscriminately, and the large, rapid flow of single-shot data that need to be processed and aggregated into a dataset, makes it difficult to assess the quality of a measurement in real time. In post-analysis it is then necessary to devise unbiased criteria to select or reject datasets, and to assign the weight with which they enter the analysis. One such case was the measurement of the lifetime of Intermolecular Coulombic Decay in the weakly-bound neon dimer. We report on the method we used to accomplish this goal for the pump-probe delay scans that constitute the core of the measurement; namely we report on the use of simple auto- and cross-correlation techniques based on the general concept of “matched filter”. We are able to unambiguously assess the signal-to-noise ratio (SNR) of each scan, which then becomes the weight with which a scan enters the average of multiple scans. We also observe a clear gap in the values of SNR, and we discard all the scans below a SNR of 0.45. We are able to generate an average delay scan profile, suitable for further analysis: in our previous work we used it for comparison with theory. Here we argue that the method is sufficiently simple and devoid of human action to be applicable not only in post-analysis, but also for the real-time assessment of the quality of a dataset.

Delocalization of a Vacancy across Two Neon Atoms Bound by the van der Waals Force.

We experimentally study 2p photoionization of neon dimers (Ne_{2}) at a photon energy of hν=36.56 eV. By postselection of ionization events which lead to a dissociation into Ne^{+}+Ne we obtain the photoelectron angular emission distribution in the molecular frame. This distribution is symmetric with respect to the direction of the charged vs neutral fragment. It shows an inverted Cohen-Fano double slit interference pattern of two spherical waves emitted coherently but with opposite phases from the two atoms of the dimer.

Transport and dissociation of neon dimer cations in neon gas: a joint dynamical and Monte Carlo modeling

A hybrid dynamical method based on the classical treatment of nuclei and the quantum treatment of electrons was used to calculate momentum transfer and dissociation cross-sections for collisions of neon dimer cations with neon atoms. For the inclusion of nuclear quantum effects, a semi-empirical factor was introduced to correct the hybrid momentum transfer cross-sections at low collision energies. Both uncorrected and quantum corrected hybrid cross-sections were used to calculate the mobility, and longitudinal and transverse characteristic diffusion energies over a wide range of the reduced electric field. Furthermore, the dissociation rate constant was calculated and compared to measured data. In addition, an approximate inverse method based on an effective isotropic interaction potential was also used to calculate the momentum transfer cross-sections and related transport data.

Defining the contributions of permanent electrostatics, Pauli repulsion, and dispersion in density functional theory calculations of intermolecular interaction energies.

In energy decomposition analysis of Kohn-Sham density functional theory calculations, the so-called frozen (or pre-polarization) interaction energy contains contributions from permanent electrostatics, dispersion, and Pauli repulsion. The standard classical approach to separate them suffers from several well-known limitations. We introduce an alternative scheme that employs valid antisymmetric electronic wavefunctions throughout and is based on the identification of individual fragment contributions to the initial supersystem wavefunction as determined by an energetic optimality criterion. The density deformations identified with individual fragments upon formation of the initial supersystem wavefunction are analyzed along with the distance dependence of the new and classical terms for test cases that include the neon dimer, ammonia borane, water-Na(+), water-Cl(-), and the naphthalene dimer.

Spin-orbit ZORA and four-component Dirac-Coulomb estimation of relativistic corrections to isotropic nuclear shieldings and chemical shifts of noble gas dimers.

Hartree-Fock and density functional theory with the hybrid B3LYP and general gradient KT2 exchange-correlation functionals were used for nonrelativistic and relativistic nuclear magnetic shielding calculations of helium, neon, argon, krypton, and xenon dimers and free atoms. Relativistic corrections were calculated with the scalar and spin-orbit zeroth-order regular approximation Hamiltonian in combination with the large Slater-type basis set QZ4P as well as with the four-component Dirac-Coulomb Hamiltonian using Dyall's acv4z basis sets. The relativistic corrections to the nuclear magnetic shieldings and chemical shifts are combined with nonrelativistic coupled cluster singles and doubles with noniterative triple excitations [CCSD(T)] calculations using the very large polarization-consistent basis sets aug-pcSseg-4 for He, Ne and Ar, aug-pcSseg-3 for Kr, and the AQZP basis set for Xe. For the dimers also, zero-point vibrational (ZPV) corrections are obtained at the CCSD(T) level with the same basis sets were added. Best estimates of the dimer chemical shifts are generated from these nuclear magnetic shieldings and the relative importance of electron correlation, ZPV, and relativistic corrections for the shieldings and chemical shifts is analyzed.

Two-photon resonant excitation of interatomic coulombic decay in neon dimers

The recent availability of intense and ultrashort extreme ultraviolet sources opens up the possibility of investigating ultrafast electronic relaxation processes in matter in an unprecedented regime. In this work we report on the observation of two-photon excitation of interatomic Coulombic decay (ICD) in neon dimers using the tunable intense pulses delivered by the free electron laser FERMI. The unique characteristics of FERMI (narrow bandwidth, spectral stability, and tunability) allow one to resonantly excite specific ionization pathways and to observe a clear signature of the ICD mechanism in the ratio of the ion yield created by Coulomb explosion. The present experimental results are explained by ab initio electronic structure and nuclear dynamics calculations.

Time-resolved study of ICD in Ne dimers using FEL radiation

Interatomic Coulombic Decay (ICD) is a relaxation phenomenon, which takes place in weakly bound atomic and molecular systems, typically within a few to hundreds of femtoseconds depending on the system and the particular decay mechanism. The creation of ICD-active states requires the production of highly excited systems, usually populated by innershell ionization or excitation. To this end, XUV and X-ray radiation from synchrotrons was conventionally applied for the majority of experiments due to the desired state-selective ionization of certain sub-shells. The advent of Free-Electron Lasers (FELs) has enabled an entirely new class of experiments, which finally allow to trace ICD directly in the time domain due to the femtosecond pulse duration. Within this paper, the first time-resolved ICD measurement using an XUV-pump–XUV-probe scheme will be discussed in detail. The experiment was performed on neon dimers and ICD was triggered by removing a 2s electron from one of the neon atoms using a 58eV pulse from the FEL in Hamburg (FLASH). The onset of ICD was probed with a delayed copy of the trigger pulse that further ionized one of the two Ne+ ions emerging after ICD. Thus, the delay-dependent yield of coincident Ne++Ne2+ ion pairs contains the lifetime of the 2s-innershell vacancy decaying via ICD. The result of 150fs ± 50fs is in good agreement with theory but only for those calculations that explicitly take nuclear motion into account.

Interatomic Coulombic Decay in Low Energy Ion-Dimer Collisions

SynopsisLow energy collisions between ions and neon dimers are investigated, focusing attention on the dissociation of the ionized molecular target. The results show a clear dependance of the processes with the projectile charge state and an unexpected appearance of inner-shell single electron capture from one site of the neon dimer leading to interatomic coulombic decay process and low energy electrons emission.

The fragmentation of neon dimer induced by low energy O6+ ions

SynopsisWe have investigated the fragmentation channels Ne+-Ne+ and Ne2+-Ne+ from Ne2 dimer induced by O6+ ions at impact energy of 15 keV/u. The Ne+-Ne+ fragmentation channel mainly results from three pathways: the interatomic coulombic decay of Ne2+(2p4)-Ne initial state, the direct columbic explosion and the radiation charge transfer. Nevertheless, the fragmentation channel Ne2+-Ne+ is predominant by the direct coulombic explosion.

Interatomic Coulombic decay as a new source of low energy electrons in slow ion-dimer collisions.

We provide the experimental evidence that the single electron capture process in slow collisions between O^{3+} ions and neon dimer targets leads to an unexpected production of low-energy electrons. This production results from the interatomic Coulombic decay process, subsequent to inner-shell single electron capture from one site of the neon dimer. Although pure one-electron capture from the inner shell is expected to be negligible in the low collision energy regime investigated here, the electron production due to this process overtakes by 1 order of magnitude the emission of Auger electrons by the scattered projectiles after double-electron capture. This feature is specific to low charge states of the projectile: similar studies with Xe^{20+} and Ar^{9+} projectiles show no evidence of inner-shell single-electron capture. The dependence of the process on the projectile charge state is interpreted using simple calculations based on the classical over the barrier model.