Vacuum ultraviolet laser spectroscopy. V. Rovibronic spectra of Ar2and constants of the ground and excited states

The borderless world of chemical bonding across the van der Waals crust and the valence region

The definition of the van der Waals crust as the spherical section between the atomic radius and the van der Waals radius of an element is discussed and a survey of the application of the penetration index between two interacting atoms in a wide variety of covalent, polar, coordinative or noncovalent bonding situations is presented. It is shown that this newly defined parameter permits the comparison of bonding between pairs of atoms in structural and computational studies independently of the atom sizes.

The ab initio potential energy curves of atom pairs and transport properties of high-temperature vapors of Cu and Si and their mixtures with He, Ar, and Xe gases

The potential energy curves (PECs) for the homonuclear He-He, Ar-Ar, Cu-Cu, and Si-Si dimers, as well as heteronuclear Cu-He, Cu-Ar, Cu-Xe, Si-He, Si-Ar, and Si-Xe dimers, are obtained in quantum Monte Carlo (QMC) calculations. It is shown that the QMC method provides the PECs with an accuracy comparable with that of the state-of-the-art coupled cluster singles and doubles with perturbative triples corrections [CCSD(T)] calculations. The QMC data are approximated by the Morse long range (MLR) and (12-6) Lennard-Jones (LJ) potentials. The MLR and LJ potentials are used to calculate the deflection angles in binary collisions of corresponding atom pairs and transport coefficients of Cu and Si vapors and their mixtures with He, Ar, and Xe gases in the range of temperature from 100 K to 10 000 K. It is shown that the use of the LJ potentials introduces significant errors in the transport coefficients of high-temperature vapors and gas mixtures. The mixtures with heavy noble gases demonstrate anomalous behavior when the viscosity and thermal conductivity can be larger than that of the corresponding pure substances. In the mixtures with helium, the thermal diffusion factor is found to be unusually large. The calculated viscosity and diffusivity are used to determine parameters of the variable hard sphere and variable soft sphere molecular models as well as parameters of the power-law approximations for the transport coefficients. The results obtained in the present work include all information required for kinetic or continuum simulations of dilute Cu and Si vapors and their mixtures with He, Ar, and Xe gases.

Accurate Property Prediction by Second Order Perturbation Theory: The REMP and OO-REMP Hybrids

The prediction of molecular properties such as equilibrium structures or vibrationalwavenumbers is a routine task in computational chemistry. If very high accuracy is required, however, the use of computationally demanding ab initio wavefunction methods is mandatory. We present property calculations utilizing the REMP and OO-REMP hybrid perturbation theories showing that with the latter approach, very accurate results are obtained at second order in perturbation theory. Specifically, equilibrium structures and harmonic vibrational wavenumbers as well as dipole moments of closed and open shell molecules were calculated and compared to the best available experimental results or very accurate calculations.OO-REMP is capable of predicting bond lengths of small closed and open shell molecules with an accuracy of 0.2 pm and 0.5 pm, respectively, often within the range of experimental uncertainty. Equilibrium harmonic vibrational wavenumbers are predicted with an accuracy better than 20 cm−1 . Dipole moments of small closed and open shell molecules are reproduced with a relative error of less than 3 %. Across all investigated properties it turns out that a 20 %:80 % MP:RE mixing ratio consistently provides the best results. This is in line with our previous findings featuring closed and open shell reaction energies.

Rotational spectroscopy of the argon dimer by time-resolved Coulomb explosion imaging of rotational wave packets

We report time-domain rotational spectroscopy of the argon dimer, Ar₂, by implementing time-resolved Coulomb explosion imaging of rotational wave packets. The rotational wave packets are created in Ar₂ with a linearly polarized, nonresonant, ultrashort laser pulse, and their spatiotemporal evolution is fully characterized by measuring angular distribution of the fragmented Ar⁺ promptly ejected from Ar₂²⁺ generated by the more intense probe pulse. The pump-probe measurements have been carried out up to a delay time of 16 ns. The alignment parameters, derived from the observed images, exhibit periodic oscillation lasting for more than 15 ns. The pure rotational spectrum of Ar₂ is obtained by Fourier transformation of the time traces of the alignment parameters. The frequency resolution in the spectrum is about 90 MHz, the highest ever achieved for Ar₂. The rotational constant and the centrifugal distortion constant are determined with much improved precision than the previous experimental results: B₀ = 1.72713 ± 0.00009 GHz and D₀ = 0.0310 ± 0.0005 MHz. The present B₀ value does not match within the quoted experimental uncertainty with that from the VUV spectroscopy, so far accepted as an experimental reference to assess theories. The present improved constants would stand as new references to calibrate state-of-the-art theoretical investigations and an indispensable experimental source for the construction of an accurate empirical intermolecular potential.

Understanding noncovalent bonds and their controlling forces

The fundamental underpinnings of noncovalent bonds are presented, focusing on the σ-hole interactions that are closely related to the H-bond. Different means of assessing their strength and the factors that control it are discussed. The establishment of a noncovalent bond is monitored as the two subunits are brought together, allowing the electrostatic, charge redistribution, and other effects to slowly take hold. Methods are discussed that permit prediction as to which site an approaching nucleophile will be drawn, and the maximum number of bonds around a central atom in its normal or hypervalent states is assessed. The manner in which a pair of anions can be held together despite an overall Coulombic repulsion is explained. The possibility that first-row atoms can participate in such bonds is discussed, along with the introduction of a tetrel analog of the dihydrogen bond.

HFLD: A Nonempirical London Dispersion-Corrected Hartree-Fock Method for the Quantification and Analysis of Noncovalent Interaction Energies of Large Molecular Systems †

A nonempirical quantum mechanical method for the efficient and accurate quantification and analysis of intermolecular interactions is presented and tested on existing benchmark sets. The leading idea here is to focus on the intermolecular part of the correlation energy that contains the all-important London dispersion (LD) interaction. To keep the cost of the method low, essentially at the level of a Hartree-Fock (HF) calculation, the intramolecular part of the correlation energy is neglected. We also neglect the nondispersive parts of the intermolecular correlation energy. This scheme that we denote as Hartree-Fock plus London dispersion (HFLD) can be readily realized on the basis of the recently reported multilevel implementation of the domain-based local pair natural orbital coupled-cluster (DLPNO-CC) theory in conjunction with the well-established local energy decomposition (LED) analysis. The accuracy and efficiency of the HFLD method are evaluated on rare gas dimers, on the S66 and L7 benchmark sets of noncovalent interactions, and on an additional set (LP14) consisting of bulky Lewis pairs held together by intermolecular interactions of various strengths, with interaction energies ranging from -8 to -107 kcal/mol. It is first shown that the LD energy calculated with this approach is essentially identical to that obtained from the full DLPNO-CCSD(T)/LED calculation, with a mean absolute error of 0.2 kcal/mol on the S66 benchmark set. Moreover, in terms of the overall interaction energies, the HFLD method shows an efficiency that is comparable to that of the HF method, while retaining an accuracy between that of the DLPNO-CCSD and DLPNO-CCSD(T) schemes. Since the underlying DLPNO-CCSD method is linear scaling with respect to the system size, the HFLD approach also does not lead to new bottlenecks for large systems. As an illustrative example of its efficiency, the HFLD scheme was applied to the interaction between the substrate and the residues in the active site of the cyclohexanone monooxygenase enzyme. The excellent cost/performance ratio indicates that the HFLD method opens new avenues for the accurate calculation and analysis of noncovalent interaction energies in large molecular systems.

Tracing charge transfer in argon dimers by XUV-pump IR-probe experiments at FLASH

Charge transfer (CT) at avoided crossings of excited ionized states of argon dimers is observed using a two-color pump-probe experiment at the free-electron laser in Hamburg (FLASH). The process is initiated by the absorption of three 27-eV-photons from the pump pulse, which leads to the population of Ar^2+*-Ar states. Due to nonadiabatic coupling between these one-site doubly ionized states and two-site doubly ionized states of the type Ar^+*-Ar⁺, CT can take place leading to the population of the latter states. The onset of this process is probed by a delayed infrared (800 nm) laser pulse. The latter ionizes the dimers populating repulsive Ar²⁺ -Ar⁺ states, which then undergo a Coulomb explosion. From the delay-dependent yields of the obtained Ar²⁺ and Ar⁺ ions, the lifetime of the charge-transfer process is extracted. The obtained experimental value of (531 ± 136) fs agrees well with the theoretical value computed from Landau-Zener probabilities.

Effect of Electron Correlation on Intermolecular Interactions: A Pair Natural Orbitals Coupled Cluster Based Local Energy Decomposition Study

The development of post-Hartree-Fock (post-HF) energy decomposition schemes that are able to decompose the HF and correlation components of the interaction energy into chemically meaningful contributions is a very active field of research. One of the challenges is to provide a clear-cut quantification to the elusive London dispersion component of the intermolecular interaction. London dispersion is well-known to be a pure correlation effect, and as such it is not properly described by mean field theories. In this context, we have recently developed the local energy decomposition (LED) analysis, which provides a chemically meaningful decomposition of the interaction energy between two or more fragments computed at the domain-based local pair natural orbitals coupled cluster (DLPNO-CCSD(T)) level of theory. In this work, this scheme is used in conjunction with other interpretation tools to study a series of molecular adducts held together by intermolecular interactions of different natures. The HF and correlation components of the interaction energy are thus decomposed into a series of chemically meaningful contributions. Emphasis is placed on discussing the physical effects associated with the inclusion of electron correlation. It is found that four distinct physical effects can contribute to the magnitude of the correlation part of intermolecular binding energies (Δ E_int^C): (i) London dispersion, (ii) the correlation correction to the reference induction energy, (iii) the correlation correction to the electron sharing process, and (iv) the correlation correction to the permanent electrostatics. As expected, the largest contribution to the correlation binding energy of neutral, apolar molecules is London dispersion, as in the argon dimer case. In contrast, the correction for the HF induction energy dominates Δ E_int^C in systems in which an apolar molecule interacts with charged or strongly polar species, as in Ar-Li⁺. This effect has its origin in the systematic underestimation of polarizabilities at the HF level of theory. For similar reasons, electron sharing largely contributes to the correlation binding energy of covalently bound molecules, as in the beryllium dimer case. Finally, the correction for HF permanent electrostatics significantly contributes to Δ E_int^C in molecules with strong dipoles, such as water and hydrogen fluoride dimers. This effect originates from the characteristic overestimation of dipole moments at the HF level of theory, leading in some cases to positive Δ E_int^C values. Our results are apparently in contrast to the widely accepted view that Δ E_int^C is typically dominated by London dispersion, at least, in the strongly interacting region. Clearly, post-HF energy decomposition schemes are very powerful tools to analyze, categorize, and understand the various contributions to the intermolecular interaction energy. Hopefully, this will eventually lead to insights that are helpful in designing systems with tailored properties. All analysis tools presented in this work will be available free of charge in the next release of the ORCA program package.

Determination of Interatomic Potentials of He_{2}, Ne_{2}, Ar_{2}, and H_{2} by Wave Function Imaging

We report on a direct method to measure the interatomic potential energy curve of diatomic systems. A cold target recoil ion momentum spectroscopy reaction microscope was used to measure the squares of the vibrational wave functions of H_{2}, He_{2}, Ne_{2}, and Ar_{2}. The Schrödinger equation relates the curvature of the wave function to the potential V(R) and therefore offers a simple but elegant way to extract the shape of the potential.

Predicting vapor liquid equilibria using density functional theory: A case study of argon

Predicting vapor liquid equilibria (VLE) of molecules governed by weak van der Waals (vdW) interactions using the first principles approach is a significant challenge. Due to the poor scaling of the post Hartree-Fock wave function theory with system size/basis functions, the Kohn-Sham density functional theory (DFT) is preferred for systems with a large number of molecules. However, traditional DFT cannot adequately account for medium to long range correlations which are necessary for modeling vdW interactions. Recent developments in DFT such as dispersion corrected models and nonlocal van der Waals functionals have attempted to address this weakness with a varying degree of success. In this work, we predict the VLE of argon and assess the performance of several density functionals and the second order Møller-Plesset perturbation theory (MP2) by determining critical and structural properties via first principles Monte Carlo simulations. PBE-D3, BLYP-D3, and rVV10 functionals were used to compute vapor liquid coexistence curves, while PBE0-D3, M06-2X-D3, and MP2 were used for computing liquid density at a single state point. The performance of the PBE-D3 functional for VLE is superior to other functionals (BLYP-D3 and rVV10). At T = 85 K and P = 1 bar, MP2 performs well for the density and structural features of the first solvation shell in the liquid phase.

Interatomic decay of inner-valence ionized states in ArXe clusters: relativistic approach

In this work we investigate interatomic electronic decay processes taking place in mixed argon-xenon clusters upon the inner-valence ionization of an argon center. We demonstrate that both interatomic Coulombic decay and electron-transfer mediated decay (ETMD) are important in larger rare gas clusters as opposed to dimers. Calculated secondary electron spectra are shown to depend strongly on the spin-orbit coupling in the final states of the decay as well as the presence of polarizable environment. It follows from our calculations that ETMD is a pure interface process taking place between the argon-xenon layers. The interplay of all these effects is investigated in order to arrive at a suitable physical model for the decay of inner-valence vacancies taking place in mixed ArXe clusters.

The dissociation of NO-Ar(A) from around threshold to 200 cm(-1) above threshold

We report an investigation of the dissociation of A state NO-Ar at energies from 23 cm(-1) below the dissociation energy to 200 cm(-1) above. The NO product rotational distributions show population in states that are not accessible with the energy available for excitation from the NO ground state. This effect is observed at photon energies from below the dissociation energy up to approximately 100 cm(-1) above it. Translational energy distributions, extracted from velocity map images of individual rotational levels of the NO product, reveal contributions from excitation of high energy NO-Ar X states at all the excess energies probed, although this diminishes with increasing photon energy and is quite small at 200 cm(-1), the highest energy studied. These translational energy distributions show that there are contributions arising from population in vibrational levels up to the X state dissociation energy. We propose that the reason such sparsely populated levels contribute to the observed dissociation is a considerable increase in the transition moment, via the Franck-Condon factor associated with these highly excited states, which arises because of the quite different geometries in the NO-Ar X and A states. This effect is likely to arise in other systems with similarly large geometry changes.

Vibrational relaxation of NO-(v=1) in icosahedral (Ar)12NO- clusters

Relaxation dynamics of NO(-)(v=1) in icosahedral (Ar)(12)NO(-) clusters are studied using classical dynamics and semiclassical procedures over the temperature range of 100-300 K. The minimum energy of the equilibrium configuration (-9875 cm(-1)) needed in the study is determined by varying the cluster size z in (Ar)(z)NO(-). NO(-)(v=1) is embedded in the cluster, which is filled with low frequency motions: 39 cm(-1) for the argon modes, 77 cm(-1) for the Arc...NO(-) substructure vibration, 109 cm(-1) for the librational frequency of restricted rotation, and 128 cm(-1) for oscillatory local translation. Dynamics calculations show that in the early time period (<20 ps), part of the vibrational energy rapidly transfers to rotation, but most energy transfers to Ar atoms on a long time scale (approximately 1 ns). The long time scale leads to the relaxation rates of 0.403 ns(-1) at 100 K and 0.453 ns(-1) at 300 K. The rates calculated using analytical formulations vary nearly linearly from 0.288 ns(-1) at 100 K to 0.832 ns(-1) at 300 K. Although the temperature dependence is stronger in the latter, both approaches give the rates on a nanosecond time scale. The principal energy transfer pathway is from NO(-) vibration to Ar vibrations via oscillatory local translation, while the NO(-) rotation is in a librational state. The energy transfer probabilities are two orders of magnitude larger than the vibration-to-translation probabilities in the gas phase collision Ar-NO(-)(v=1).

Eigen energies and the statistical distributions of the rovibrational levels of the bosonic van der Waals argon trimer

The eigen energies and the statistical distributions of the rovibrational levels (J = 0-2) of the bosonic van der Waals argon trimer were calculated using a full angular momentum three-dimensional finite element method. The influence of interatomic potentials on the vibrational levels and statistical properties of the trimer was discussed.

Beyond the Lennard-Jones model: a simple and accurate potential function probed by high resolution scattering data useful for molecular dynamics simulations

Scattering data, measured for rare gas-rare gas systems under high angular and energy resolution conditions, have been used to probe the reliability of a recently proposed interaction potential function, which involves only one additional parameter with respect to the venerable Lennard-Jones (LJ) model and is hence called Improved Lennard-Jones (ILJ). The ILJ potential eliminates most of the inadequacies at short- and long-range of the LJ model. Further reliability tests have been performed by comparing calculated vibrational spacings with experimental values and calculated interaction energies at short-range with those obtained from the inversion of gaseous transport properties. The analysis, extended also to systems involving ions, suggests that the ILJ potential model can be used to estimate the behavior of unknown systems and can help to assess the different role of the leading interaction components. Moreover, due to its simple formulation, the physically reliable ILJ model appears to be particularly useful for molecular dynamics simulations of both neutral and ionic systems.

Pulsed-field-ionization zero-kinetic-energy photoelectron spectroscopy of metastable He2: Ionization potential and rovibrational structure of He2+

A supersonic beam of metastable He(*) atoms and He(2) (*) a (3)Sigma(u) (+) molecules has been generated using a pulsed discharge at the exit of a pulsed valve prior to the gas expansion into vacuum. Pulsed-field-ionization zero-kinetic-energy photoelectron spectra of the He(2) (+) X(+) (2)Sigma(u) (+) (v(+)=0-2)<--He(2) (*) a (3)Sigma(u) (+) (v(")=0-2) transitions and photoionization spectra of He(2) (*) in the vicinity of the lowest ionization thresholds have been recorded. The energy level structures of (4)He(2) (+) X(+) (2)Sigma(u) (+) (v(+)< or =2,N(+)< or =23) and (3)He(2) (+) X(+) (2)Sigma(u) (+) (v(+)=0,N(+)< or =11) have been determined, and an accurate set of molecular constants for all isotopomers of He(2) (+) has been derived in a global analysis of all spectroscopic data reported to date on the low vibrational levels of He(2) (+). The analysis of the photoionization spectrum by multichannel quantum defect theory has provided a set of parameters describing the threshold photoionization dynamics.

An optimized algebraic basis for molecular potentials

The computation of vibrational spectra of diatomic molecules through the exact diagonalization of algebraically determined matrices based on powers of Morse coordinates is made substantially more efficient by choosing a properly adapted quantum mechanical basis, specifically tuned to the molecular potential. A substantial improvement is achieved while still retaining the full advantage of the simplicity and numerical light-weightedness of an algebraic approach. In the scheme we propose, the basis is parametrized by two quantities which can be adjusted to best suit the molecular potential through a simple minimization procedure.

Calculation of argon trimer rovibrational spectrum

Rovibrational spectra of Ar3 are computed for total angular momenta up to J=6 using row-orthonormal hyperspherical coordinates and an expansion of the wave function on hyperspherical harmonics. The sensitivity of the spectra to the two-body potential and to the three-body corrections is analyzed. First, the best available semiempirical pair potential (HFDID1) is compared with our recent ab initio two-body potential. The ab initio vibrational energies are typically 1-2 cm-1 higher than the semiempirical ones, which is related to the slightly larger dissociation energy of the semiempirical potential. Then, the Axilrod-Teller asymptotic expansion of the three-body correction is compared with our newly developed ab initio three-body potential. The difference is found smaller than 0.3 cm-1. In addition, we define approximate quantum numbers to describe the vibration and rotation of the system. The vibration is represented by a hyper-radial mode and a two-degree-of-freedom hyperangular mode, including a vibrational angular momentum defined in an Eckart frame. The rotation is described by the total angular momentum quantum number, its projection on the axis perpendicular to the molecular plane, and a hyperangular internal momentum quantum number, related to the vibrational angular momentum by a transformation between Eckart and principal-axes-of-inertia frames. These quantum numbers provide a qualitative understanding of the spectra and, in particular, of the impact of the nuclear permutational symmetry of the system (bosonic with zero nuclear spin). Rotational constants are extracted from the spectra and are shown to be accurate only for the ground hyperangular mode.

The Nature of Variations of Ammonia Proton Affinity in an Argon Environment

The presence of solvent molecules, even inert, may significantly influence processes taking place in the gas phase. The reason for the solvent activity may be found in studies of complexes originating from microsolvation. The structure, thermodynamics, and vibrational properties of NH(4)(+)Ar(n) (n = 0-5) and NH(3)Ar(n) (n = 0-4) complexes are presented in this work with the aim to elucidate the effect of microsolvation on protonation/deprotonation processes. The relation between the nature of interactions in cationic and neutral clusters and the proton affinity is studied as a function of the number of ligands in complexes.

On the influence of microsolvation by argon atoms on the electron affinity properties of water dimer

This work provides a comparison of neutral (H2O)2Ar(n) and negatively charged (H2O)(2-)Ar(n) complexes. The excess electron stabilizes the complexes and leads to the trans to cis rearrangement within the water dimer core. In the case of small complexes (n < or = 4) the microsolvation of the dimer by argon atoms arises on the trans side with respect to the donor water molecule. The stabilization of an excess electron is enhanced by the delocalization of the electronic charge density due to microsolvation. The process of cis to trans rotation is induced by the electric field of the approaching negative charge. The interaction energy decomposition suggests a more ionic character of binding in the negatively charged complexes. The attachment of an electron is controlled by the correlation energy.