Chemistry, astronomy and physics of H3+

A nearly complete treatment of the effect of non-adiabaticity on rovibrational energies of H3+ (Part III)

In this article, significant contributions of non-adiabaticity for the rovibrational bound states up to 25 000 cm-1 and total angular momentum J = 0-20 of H3+ are investigated. A coupled-perturbed full configuration interaction (CP-FCI) treatment is applied to calculate all couplings between electronic states caused by the nuclear motion. These derivative couplings were evaluated up to the second order by means of a perturbation treatment and include all nuclear Cartesian first and second derivatives of the electronic wavefunctions. In particular, the coupling of special derivatives with respect to r and R in the Jacobi coordinate representation is more significant than thought. The perturbation approach is especially optimal for the treatment of weak non-adiabaticity in case of rovibrational energies in H3+ and had not been available before for H3+ or other triatomics. Using exclusively Gaussian basis functions for CP-FCI appears to be sufficient, because explicit correlated wavefunctions are already used for all other potential energy contributions. Our work is an extension of earlier non-adiabatic investigations based on first derivative couplings of electronic states that led to the concept of geometry-dependent effective nuclear masses and which needs only a single potential energy surface for the dynamics. The implementation allows us to include all non-adiabatic effects up to the order of O(μ-2), μ being the reduced nuclear mass. Our treatment works for any isotopologue and for the whole potential energy curve or surface. By this treatment, a further reduction in deviations to experimental data for most rovibrational levels to less than 0.1 cm-1 is possible. For the related transition frequencies, 1366 of 1720 known rovibrational transitions in H3+ have deviations less than 0.1 cm-1 without using any empirically adjustable parameters or optimizing the nuclear mass for a specific transition. For many questionable assignments (deviations >0.3 cm-1) of observed transitions in H3+, a new labeling is proposed.

Trihydrogen Cation Helium Clusters: A New Potential Energy Surface

We present a new analytical potential energy surface (PES) for the interaction between the trihydrogen cation and a He atom,H 3 + - H e ${{H}_{3}^{+}-He}$ , in its electronic ground state. The proposed PES has been built as a sum of two contributions: a polarization energy term due to the electric field generated by the molecular cation at the position of the polarizable He atom, and an exchange-repulsion and dispersion interactions represented by a sum of "atom-bond" potentials between the three bonds ofH 3 + ${{H}_{3}^{+}}$ and the He atom. All parameters of this new PES have been chosen and fitted from data obtained from high-level ab-initio calculations. Using this new PES plus the Aziz-Slaman potential for the interaction between Helium atoms and assuming pair-wise interactions, we carry out classical Basin-Hopping (BH) global optimization, semiclassical BH with Zero Point Energy corrections, and quantum Diffusion Monte Carlo simulations. We have found the minimum energy configurations of small He clusters doped withH 3 + ${{H}_{3}^{+}}$ ,H 3 + H e N ${{H}_{3}^{+}{\left(He\right)}_{N}}$ , with N=1-16. The study of the energies of these clusters allows us to find a pronounced anomaly for N=12, in perfect agreement with previous experimental findings, which we relate to a greater relative stability of this aggregate.

An efficient algorithm for capturing quantum effects in classical reactive scattering: application to D + H+3 → H2D+ + H

Motivated by a recent semiclassical analysis of chemical reaction thresholds [Bonnet et al., J. Chem. Phys., 2022, 157, 094114], we present an efficient algorithm for including zero-point energy (ZPE) effects in classical reactive scattering. The algorithm is an extension of the quasi-classical trajectory (QCT) Gaussian binning method. We apply it to the astrophysically important D + H+3 reaction, where there are significant quantum effects and where application of other methods is problematic [Braunstein et al., Phys. Chem. Chem. Phys., 2022, 24, 5489]. The rate constants computed with the new, general algorithm closely match recent Ring Polymer Molecular Dynamics (RPMD) [Bulut et al., J. Phys. Chem. A, 2019, 123, 8766] and experimentally derived [Bowen et al., J. Chem. Phys., 2021, 154, 084307] ones spanning ∼4 orders of magnitude from 70 to 1500 K.

Capturing quantum effects with quasi-classical trajectories in the D + H+3 → H2D+ + H reaction

We present quasi-classical trajectory (QCT) cross sections, rate constants, and product state distributions for the D + H+3 → H₂D⁺ + H reaction. Using the same H+4 potential surface, the rate constants obtained from several QCT-based methods correcting for zero-point effects by Gaussian binning the product H₂D⁺ are compared to ring polymer molecular dynamics (RPMD) rate constants [Bulut et al., J. Phys. Chem. A, 2019, 123, 8766] which include quantum effects and to recent experimentally derived rate constants [Bowen et al., J. Chem. Phys., 2021, 154, 084307]. QCT with standard binning predicts rate constants that increase slowly as the temperature decreases from 1500 to 100 K. In contrast, the RPMD rate constants decrease rapidly with decreasing temperature. By 100 K, the QCT standard binning rate constant is more than 3 orders of magnitude larger than the RPMD rate constant. We show that QCT with Gaussian binning and proper normalization captures the zero-point effects and reproduces the RPMD rate constants over a large temperature range. Furthermore, the simple technique of counting only reactive trajectories with vibrational energy above the product zero-point energy matches the RPMD results well down to ∼300 K. The present Gaussian binned rate constants are in fair agreement with new experimentally derived rate constants from 100 to 1500 K. However, because the Gaussian binned rate constants do not include tunneling, important at lower temperatures, and the RPMD and experimentally derived rate constants have significant differences, the roles of the competing effects of zero-point energy, internal excitation of the H+3, and quantum tunneling are not simple and require further study for a consistent picture of the dynamics. Since rate constants for complex forming reactions, such as the title reaction, are difficult to converge with RPMD, alternative QCT-based methods, which include quantum effects and in addition provide product state distributions as described here, are highly desirable.

Near-resonant effects in the quantum dynamics of the H + H2 + → H2 + H+ charge transfer reaction and isotopic variants

The non-adiabatic quantum dynamics of the H + H₂ ⁺ → H₂ + H⁺ charge transfer reactions, and some isotopic variants, is studied with an accurate wave packet method. A recently developed 3 × 3 diabatic potential model is used, which is based on very accurate ab initio calculations and includes the long-range interactions for ground and excited states. It is found that for initial H₂ ⁺(v = 0), the quasi-degenerate H₂(v' = 4) non-reactive charge transfer product is enhanced, producing an increase in the reaction probability and cross section. It becomes the dominant channel from collision energies above 0.2 eV, producing a ratio between v' = 4 and the rest of v's, which that increase up to 1 eV. The H + H₂ ⁺ → H₂ ⁺ + H exchange reaction channel is nearly negligible, while the reactive and non-reactive charge transfer reaction channels are of the same order, except that corresponding to H₂(v' = 4), and the two charge transfer processes compete below 0.2 eV. This enhancement is expected to play an important vibrational and isotopic effect that needs to be evaluated. For the three proton case, the problem of the permutation symmetry is discussed when using reactant Jacobi coordinates.

Atmospheric implications of the lack of H3+ detection at Neptune

H3+ has been detected at all of the solar system giant planets aside from Neptune. Current observational upper limits imply that there is far less H3+ emission at Neptune than rudimentary modelling would suggest. Here, we explore via modelling a range of atmospheric conditions in order to find some that could be consistent with observational constraints. In particular, we consider that the upper atmosphere might be much cooler than it was during the 1989 Voyager 2 encounter, and we examine the impact of an enhanced influx of external material that could act to reduce H3+ density. Resulting ionosphere models that are consistent with existing H3+ observational constraints have an exospheric temperature of 450 K or less, 300 K lower than the Voyager 2 value. Alternatively, if a topside CO influx of 2 × 108 cm-2 s-1 is imposed, the upper atmospheric temperature can be higher, up to 550 K. The potential cooling of Neptune's atmosphere is relevant for poorly understood giant planet thermospheric energetics, and would also impact aerobreaking manoeuvers for any future spacecraft. Such a large CO influx, if present, could imply Triton is a very active moon with prominent atmospheric escape, and/or that Neptune's rings significantly modify its upper atmosphere, and the introduction of so much exogenic material would complicate interpretation of the origin of species observed in Neptune's lower atmosphere. This article is part a discussion meeting issue 'Future exploration of ice giant systems'.

Final-State-Resolved Dynamics of the H3+ + CO → H2 +HCO+/HOC+ Reaction: A Quasi-Classical Trajectory Study

The ion-molecule reaction H₃⁺ + CO → H₂ + HCO⁺/HOC⁺, which initiates the formation of crucial organic molecules, plays a key role in interstellar and circumstellar environments. In this work, the quasi-classical trajectory method is employed to study the reaction dynamics on a recently developed full-dimensional global potential energy surface (PES). The calculated product internal energy distributions and relative internal excited fractions agree reasonably well with the experimental measurements. For the two reaction channels, most of the available energy flows into the vibrational modes of HCO⁺ or HOC⁺ at low collision energies, followed by the translational mode and the rotational modes of HCO⁺ or HOC⁺. As the collision energy increases, the proportion of the product translational energy increases while the proportion of the product vibrational energy decreases. Furthermore, the CH and CO stretching modes and their combination bands are effectively excited for the product HCO⁺ while the bending mode is remarkably excited for the product HOC⁺.

Analysis of QED and non-adiabaticity effects on the rovibrational spectrum of H3 + using geometry-dependent effective nuclear masses

The influence of QED effects (including one- and two-electron Lamb-shift, Araki-Sucher term, one-loop self-energy, and finite nuclear size correction) together with non-adiabatic effects on the rovibrational bound states of H₃ ⁺ has been investigated. Non-adiabaticity is modeled by using geometry-dependent effective nuclear masses together with only one single potential energy surface. In conclusion, for rovibrational states below 20 000 cm^-1, QED and relativistic effects do nearly compensate, and a potential energy surface based on Born-Oppenheimer energies and diagonal adiabatic corrections has nearly the same quality as the one including relativity with QED; the deviations between the two approaches for individual rovibrational states are mostly below 0.02 cm^-1. The inclusion of non-adiabatic effects is important, and it reduces deviations from experiments mostly below 0.1 cm^-1.

Quantum Effects on the D + H3+ → H2D+ + H Deuteration Reaction and Isotopic Variants

The title reaction and its isotopic variants are studied using quasi-classical trajectory (QCT) (without taking into account corrections to account for the possible zero point energy breakdown) and ring polymer molecular dynamics (RPMD) methods with a full dimensional and accurate potential energy surface which presents an exchange barrier of approximately 0.144 eV. The QCT rate constant increases when the temperature decreases from 1500 to 10 K. On the contrary, the RPMD rate constant decreases with decreasing temperature, in semiquantitative agreement with recent experimental results. The present RPMD results are in between the thermal and translational experimental rate constants, extracted from the measured data to eliminate the initial vibrational excitation of H₃⁺, obtained in an arc discharge. The difference between the present RPMD results and experimental values is attributed to the possible existence of non thermal vibrational excitation of H₃⁺, not completely removed by the semiempirical model used for the analysis of the experimental results. Also, it is found that, below 200 K, the RPMD trajectories are trapped, forming long-lived collision complexes, with lifetimes longer than 1 ns. These collision complexes can fragment by either redissociating back to reactants or react to products, in the two cases tunneling through the centrifugal and reaction barriers, respectively. The contribution of the formation of the complex to the total deuteration rate should be calculated with more accurate quantum methods, as has been found recently for reactions of larger systems, and the present four atoms system is a good candidate to benchmark the adequacy of RPMD method at temperatures below 100 K.

The smallest proton-bound dimer H5+: theoretical progress

The protonated hydrogen dimer, H5+, is the smallest system including proton transfer, and has been of long-standing interest since its first laboratory observation in 1962. H5+ and its isotopologues are the intermediate complexes in deuterium fractionation reactions, and are of central importance in molecular astrophysics. The recently recorded infrared spectra of both H5+ and D5+ reveal a rich vibrational dynamics of the cations, which presents a challenge for standard theoretical approaches. Although H5+ is a four-electron ion, which makes highly accurate electronic structure calculations tractable, the construction of ab initio-based potential energy and dipole moment surfaces has proved a hard task. In the same vein, the difficulties in treating the nuclear motion could also become cumbersome due to their high dimensionality, floppiness and/or symmetry. These systems are prototypical examples for studying large-amplitude motions, as they are highly delocalized, interconverting between equivalent minima through internal rotation and proton transfer motions requiring state-of-the-art treatments. Recent advances in the computational vibrational spectroscopy of the H5+ cation and its isotopologues are reported from full quantum spectral simulations, providing important information in a rigorous manner, and open perspectives for further future investigations. This article is part of a discussion meeting issue 'Advances in hydrogen molecular ions: H3+, H5+ and beyond'.

H3+ as an ionospheric sounder of Jupiter and giant planets: an observational perspective

Thirty years of observations of [Formula: see text] on Jupiter have addressed many complex questions about the physics of the ionospheres of the giant planets. Spectroscopy, imaging and imaging spectroscopy in the infrared have allowed investigators to retrieve fundamental parameters of the ionosphere, overcoming the inherent limitations and complexities in radiative transfer, and these results are now introduced as model constraints for upper atmospheric structure and dynamics. This paper will focus on the mid-latitude emissions, which are fainter and less well studied than the auroral regions. A new analysis of VLT/ISAAC spectral imaging observations of Jupiter obtained in 2000 at 3.5 µm is presented and discussed in comparison with previous observations to show the spatial distribution of [Formula: see text] emissions compared with other atmospheric structures. Cylindrical maps of Jupiter in three different selected wavelengths show the spatial variations at different altitudes in the atmosphere, from cloud level up to the ionosphere. Evidence for fluctuations in the [Formula: see text] emissions could be due to the presence of stationary or dynamic processes. If the exact origin of these phenomena remains unidentified, several plausible mechanisms are proposed to explain the observed energy deposition and variability: future observation campaigns should deepen the understanding of these complex phenomena, in order to prepare for the future ESA/JUICE mission. This article is part of a discussion meeting issue 'Advances in hydrogen molecular ions: H₃⁺, H₅⁺ and beyond'.

Modelling H3+ in planetary atmospheres: effects of vertical gradients on observed quantities

Since its detection in the aurorae of Jupiter approximately 30 years ago, the H3+ ion has served as an invaluable probe of giant planet upper atmospheres. However, the vast majority of monitoring of planetary H3+ radiation has followed from observations that rely on deriving parameters from column-integrated paths through the emitting layer. Here, we investigate the effects of density and temperature gradients along such paths on the measured H3+ spectrum and its resulting interpretation. In a non-isothermal atmosphere, H3+ column densities retrieved from such observations are found to represent a lower limit, reduced by 20% or more from the true atmospheric value. Global simulations of Uranus' ionosphere reveal that measured H3+ temperature variations are often attributable to well-understood solar zenith angle effects rather than indications of real atmospheric variability. Finally, based on these insights, a preliminary method of deriving vertical temperature structure is demonstrated at Jupiter using model reproductions of electron density and H3+ measurements. The sheer diversity and uncertainty of conditions in planetary atmospheres prohibits this work from providing blanket quantitative correction factors; nonetheless, we illustrate a few simple ways in which the already formidable utility of H3+ observations in understanding planetary atmospheres can be enhanced. This article is part of a discussion meeting issue 'Advances in hydrogen molecular ions: H3+, H5+ and beyond'.

Hydrogen molecular ions: H3+, H5+ and beyond

Three decades after the spectroscopic detection of H3+ in space, the inspiring developments in physics, chemistry and astronomy of Hn+ (n = 3, 5, 7) systems, which led to this Royal Society Discussion Meeting, are reviewed, the present state of the art as represented by the meeting surveyed and future lines of research considered. This article is part of a discussion meeting issue 'Advances in hydrogen molecular ions: H3+, H5+ and beyond'.

Dissociative recombination of H3+ and D5. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2019;377:20180397. [PMID: 31378170 DOI: 10.1098/rsta.2018.0397] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Compared with earlier years, the dissociative recombination of [Formula: see text] has not been very actively studied in recent years. New results from afterglow experiments are quoted and compared with results from ion storage rings and theory. New results for [Formula: see text] are discussed. This article also contains some historical remarks on hydrogen and its importance for the advancement of physics and chemistry. This article is part of a discussion meeting issue 'Advances in hydrogen molecular ions: H₃⁺, H₅⁺ and beyond'.

A Ring Polymer Molecular Dynamics Approach to Study the Transition between Statistical and Direct Mechanisms in the H2 + H3+ → H3+ + H2 Reaction

Because of its fundamental importance in astrochemistry, the H₂ + H₃⁺ → H₃⁺ + H₂ reaction has been studied experimentally in a wide temperature range. Theoretical studies of the title reaction significantly lag primarily because of the challenges associated with the proper treatment of the zero-point energy (ZPE). As a result, all previous theoretical estimates for the ratio between a direct proton-hop and indirect exchange (via the H₅⁺ complex) channels deviate from the experiment, in particular, at lower temperatures where the quantum effects dominate. In this work, the ring polymer molecular dynamics (RPMD) method is applied to study this reaction, providing very good agreement with the experiment. RPMD is immune to the shortcomings associated with the ZPE leakage and is able to describe the transition from direct to indirect mechanisms below room temperature. We argue that RPMD represents a useful tool for further studies of numerous ZPE-sensitive chemical reactions that are of high interest in astrochemistry.

Preferential Isomer Formation Observed in H3+ + CO by Crossed Beam Imaging

The proton transfer reaction H3(+) + CO is one of the cornerstone chemical processes in the interstellar medium. Here, the dynamics of this reaction have been investigated using crossed beam velocity map imaging. Formyl product cations are found to be predominantly scattered into the forward direction irrespective of the collision energy. In this process, a high amount of energy is transferred to internal product excitation. By fitting a sum of two distribution functions to the measured internal energy distributions, the product isomer ratio is extracted. A small HOC(+) fraction is obtained at a collision energy of 1.8 eV, characterized by an upper limit of 24% with a confidence level of 84%. At lower collision energies, the data indicate purely HCO(+) formation. Such low values are unexpected given the previously predicted efficient formation of both HCO(+) and HOC(+) isomers for thermal conditions. This is discussed in light of the direct reaction dynamics that are observed.

Non-adiabatic couplings and dynamics in proton transfer reactions of Hn (+) systems: Application to H2+H2 (+)→H+H3 (+) collisions

Analytical derivatives and non-adiabatic coupling matrix elements are derived for Hn (+) systems (n = 3-5). The method uses a generalized Hellmann-Feynman theorem applied to a multi-state description based on diatomics-in-molecules (for H3 (+)) or triatomics-in-molecules (for H4 (+) and H5 (+)) formalisms, corrected with a permutationally invariant many-body term to get high accuracy. The analytical non-adiabatic coupling matrix elements are compared with ab initio calculations performed at multi-reference configuration interaction level. These magnitudes are used to calculate H2(v(')=0,j(')=0)+H2 (+)(v,j=0) collisions, to determine the effect of electronic transitions using a molecular dynamics method with electronic transitions. Cross sections for several initial vibrational states of H2 (+) are calculated and compared with the available experimental data, yielding an excellent agreement. The effect of vibrational excitation of H2 (+) reactant and its relation with non-adiabatic processes are discussed. Also, the behavior at low collisional energies, in the 1 meV-0.1 eV interval, of interest in astrophysical environments, is discussed in terms of the long range behaviour of the interaction potential which is properly described within the triatomics-in-molecules formalism.

Theoretical predictions on the role of the internal H3(+) rotation in the IR spectra of the H5(+) and D5(+) cations

The IR spectra of the H5(+) and D5(+) cations in the mid- and far-IR spectral regions have been recently reported by experimentalists. These spectra show very rich vibrational patterns representing a challenge for state-of-the-art theoretical methods to provide definitive interpretations of them. Using a full-dimensional quantum anharmonic treatment, within the MCTDH approach, together with ab initio potential and dipole moment surfaces, the predominant features in the spectra are assigned, completing an important part in previous theoretical and experimental comparisons. The internal rotation of the H3(+) unit by exciting the H3(+)-H2 stretching mode is found to correspond to the new calculated features at 1182, 1876, and 2139 cm(-1) of the H5(+) spectrum, leading to a consistent assignment with the experimental spectra. In the calculated spectra of both H5(+) and D5(+) clusters, the progressions in the H3(+)-H2 stretch of the shared proton and the in- and out-of- plane H3(+) rotation are demonstrated to be the main features. Such states are expected to play a central role in the low temperature hydrogen/deuterium proton hop/exchange H3(+) + H2 reactions.

First-principles simulations of vibrational states and spectra for H5(+) and D5(+) clusters using multiconfiguration time-dependent Hartree approach

Simulations of the infrared (IR) spectra of the H5(+) and D5(+) clusters are carried out in the whole energy range, using a recent, reliable "on the fly" DFT-based potential energy surface, and its corresponding dipole moment surface. For the present study we adopted a recently proposed four-dimensional quantum model to describe the proton transfer motion between the two vibrating H2 or D2 units. Time-dependent and time-independent approaches within the multiconfiguration time-dependent Hartree method are employed for investigating the vibrational dynamics of the complexes. The obtained spectra are compared with recent experimental data available for energies up to 4500 and 3500 cm(-1) for the H5(+) and D5(+), respectively. Even though the present results are based on a reduced dimensional model, the infrared spectra are shown to be in good qualitative accord with those observed experimentally. Also as the reported data are subject to the potential energy surface, comparisons with previous theoretical calculations based on an analytical ab initio parameterized surface are also presented. The differences on the topology of the potentials are discussed in connection with their effect on the spectral features. We found that the main characteristics of the experimentally observed spectra are reproduced by both surfaces, evaluating in this way the sensitivity of such computations on the quality of the underlying potential. This finding serves to connect aspects of the potential surface of these systems to their spectral complexity, and could be indicative to calibrate intrinsic errors in their calculation for future studies.

Full-dimensional quantum calculations of the dissociation energy, zero-point, and 10 K properties of H7+/D7+ clusters using an ab initio potential energy surface

Diffusion Monte Carlo (DMC) and path-integral Monte Carlo computations of the vibrational ground state and 10 K equilibrium state properties of the H7 (+)/D7 (+) cations are presented, using an ab initio full-dimensional potential energy surface. The DMC zero-point energies of dissociated fragments H5 (+)(D5 (+))+H2(D2) are also calculated and from these results and the electronic dissociation energy, dissociation energies, D0, of 752 ± 15 and 980 ± 14 cm(-1) are reported for H7 (+) and D7 (+), respectively. Due to the known error in the electronic dissociation energy of the potential surface, these quantities are underestimated by roughly 65 cm(-1). These values are rigorously determined for first time, and compared with previous theoretical estimates from electronic structure calculations using standard harmonic analysis, and available experimental measurements. Probability density distributions are also computed for the ground vibrational and 10 K state of H7 (+) and D7 (+). These are qualitatively described as a central H3 (+)/D3 (+) core surrounded by "solvent" H2/D2 molecules that nearly freely rotate.

Ground state of the H3(+) molecular ion: physics behind

Five physics mechanisms of interaction leading to the binding of the H3(+) molecular ion are identified. They are realized in a form of variational trial functions, and their respective total energies are calculated. Each of them provides subsequently the most accurate approximation for the Born–Oppenheimer (BO) ground state energy among (two–three–seven)-parametric trial functions being, correspondingly, H2-molecule plus proton (two variational parameters), H2(+)-ion plus H-atom (three variational parameters), and generalized Guillemin–Zener (seven variational parameters). These trial functions are chosen following a criterion of physical adequacy. They include the electronic correlation in the exponential form, exp(γr12), where γ is a variational parameter. Superpositions of two different mechanisms of binding are investigated, and a particular one, which is a generalized Guillemin–Zener plus H2-molecule plus proton (ten variational parameters), provides the total energy at the equilibrium of E = −1.3432 au. The superposition of three mechanisms, generalized Guillemin–Zener plus (H2-molecule plus proton) plus (H2(+)-ion plus H) (14 parameters) leads to the total energy, which deviates from the best known BO energy to 0.0004 au, it reproduces two–three significant digits in exact, non-BO total energy. In general, our variational energy agrees in two–three–four significant digits with the most accurate results available at present as well as major expectation values.