Quantum study of the CH3+ photodissociation in full-dimensional neural network potential energy surfaces.

C H 3 + , a cornerstone intermediate in interstellar chemistry, has recently been detected for the first time by using the James Webb Space Telescope. The photodissociation of this ion is studied here. Accurate explicitly correlated multi-reference configuration interaction ab initio calculations are done, and full-dimensional potential energy surfaces are developed for the three lower electronic states, with a fundamental invariant neural network method. The photodissociation cross section is calculated using a full-dimensional quantum wave packet method in heliocentric Radau coordinates. The wave packet is represented in angular and radial grids, allowing us to reduce the number of points physically accessible, requiring to push up the spurious states appearing when evaluating the angular kinetic terms, through projection technique. The photodissociation spectra, when employed in astrochemical models to simulate the conditions of the Orion bar, result in a lesser destruction of CH3+ compared to that obtained when utilizing the recommended values in the kinetic database for astrochemistry.

The role of dimers in complex forming reactions at low temperature: full dimension potential and dynamics of (H2 CO)2 +OH reaction.

The (H 2 ${{}_{2}}$ CO) 2 ${{}_{2}}$ +OH and H 2 ${{}_{2}}$ CO-OH+H 2 ${{}_{2}}$ CO reaction dynamics are studied theoretically for temperatures below 300 K. For this purpose, a full dimension potential energy surface is built, which reproduces well accurate ab initio calculations. The potential presents a submerged reaction barrier, as an example of the catalytic effect induced by the presence of the third molecule. However, quasi-classical and ring polymer molecular dynamics calculations show that the dominant channel is the dimer-exchange mechanism below 200 K, and that the reactive rate constant tends to stabilize at low temperatures, because the effective dipole of either dimer is reduced with respect to that of formaldehyde alone. The reaction complex formed at low temperatures does not live long enough to produce complete energy relaxation, as assumed in statistical theories. These results show that the reactivity of the dimers cannot explain the large rate constants measured at temperatures below 100 K.

Spin-orbit transitions in the N+(PJA3)+H2→ NH+(X2Π,4Σ-) + H(2S) reaction, using adiabatic and mixed quantum-adiabatic statistical approaches.

The cross section and rate constants for the title reaction are calculated for all the spin-orbit states of N+(PJA3) using two statistical approaches, one purely adiabatic and the other one mixing quantum capture for the entrance channel and adiabatic treatment for the products channel. This is made by using a symmetry adapted basis set combining electronic (spin and orbital) and nuclear angular momenta in the reactants channel. To this aim, accurate ab initio calculations are performed separately for reactants and products. In the reactants channel, the three lowest electronic states (without spin-orbit couplings) have been diabatized, and the spin-orbit couplings have been introduced through a model localizing the spin-orbit interactions in the N+ atom, which yields accurate results as compared to ab initio calculations, including spin-orbit couplings. For the products, 11 purely adiabatic spin-orbit states have been determined with ab initio calculations. The reactive rate constants thus obtained are in very good agreement with the available experimental data for several ortho-H2 fractions, assuming a thermal initial distribution of spin-orbit states. The rate constants for selected spin-orbit JA states are obtained, to provide a proper validation of the spin-orbit effects to obtain the experimental rate constants.

Neural network potential energy surface for the low temperature ring polymer molecular dynamics of the H2CO + OH reaction.

A new potential energy surface (PES) and dynamical study of the reactive process of H2CO + OH toward the formation of HCO + H2O and HCOOH + H are presented. In this work, a source of spurious long range interactions in symmetry adapted neural network (NN) schemes is identified, which prevents their direct application for low temperature dynamical studies. For this reason, a partition of the PES into a diabatic matrix plus a NN many-body term has been used, fitted with a novel artificial neural network scheme that prevents spurious asymptotic interactions. Quasi-classical trajectory (QCT) and ring polymer molecular dynamics (RPMD) studies have been carried on this PES to evaluate the rate constant temperature dependence for the different reactive processes, showing good agreement with the available experimental data. Of special interest is the analysis of the previously identified trapping mechanism in the RPMD study, which can be attributed to spurious resonances associated with excitations of the normal modes of the ring polymer.

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 + H2 + → H2 + 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 H2 +(v = 0), the quasi-degenerate H2(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 + H2 + → H2 + + 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 H2(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.

Three states global fittings with improved long range: singlet and triplet states of H.

Full dimensional analytical fits of the coupled potential energy surfaces for the three lower singlet and triplet adiabatic states of H are developed, providing analytic derivatives and non-adiabatic coupling matrix elements. The fits are highly accurate and include an improved description of the long range interactions, including new terms for the description of the long range in the diatomic fits and the atom-diatom dissociation channels. The fits are based on the DIM formalism including three body terms in Hamiltonian matrix elements, each of them obeying S2 permutational symmetry, where the positive charge is placed in either of the three hydrogen atoms, but the full system obeys S3 permutational symmetry, invariant under all permutations of the nuclei. The ab initio points used in the fitting are obtained from a complete basis set extrapolation, made for all electronic states. Total root mean square errors of the fits are 27 and 12 cm-1, for the singlet and triplet states, respectively. The errors in the channels are lower than 2 cm-1 and 6 cm-1 for the H + H and H+ + H2 channels respectively. The new fits have been used to calculate the rovibrational bound states of the lowest singlet and lowest triplet states showing very good agreement with previous calculations in the literature.

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 H3+, 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 H3+, 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.

Zero- and high-pressure mechanisms in the complex forming reactions of OH with methanol and formaldehyde at low temperatures.

A recent Ring Polymer Molecular Dynamics study of the reactions of OH with methanol and formaldehyde, at zero pressure and below 100 K, has shown the formation of long lived complexes, with long lifetimes, longer than 100 ns for the lower temperatures studied, 20-100 K (del Mazo-Sevillano et al., 2019). These long lifetimes support the existence of multi collision events with the He buffer-gas atoms under experimental conditions, as suggested by several transition state theory studies of these reactions. In this work we study these secondary collisions, as a dynamical approach to study pressure effects on these reactions. For this purpose, the potential energy surfaces of He with H2CO, OH, H2O and HCO are calculated at highly accurate ab initio level. The stability of some of the complexes is studied using Path Integral Molecular dynamics techniques, determining that OH-H2CO complexes can be formed up to 100 K or higher temperatures, while the weaker He-H2CO complexes dissociate at approximately 50 K. The predicted IR intensity spectra shows new features which could help the identification of the OH-H2CO complex. Finally, the He-H2CO + OH and OH-H2CO + He collisions are studied using quassi-classical trajectories, finding that the cross section to produce HCO + H2O products increases with decreasing collision energy, and that it is ten times higher in the He-H2CO + OH case.

Non-adiabatic Quantum Dynamics of the Dissociative Charge Transfer He++H2 → He+H+H.

We present the non-adiabatic, conical-intersection quantum dynamics of the title collision where reactants and products are in the ground electronic states. Initial-state-resolved reaction probabilities, total integral cross sections, and rate constants of two H2 vibrational states, v 0 = 0 and 1, in the ground rotational state (j 0 = 0) are obtained at collision energies E coll ≤ 3 eV. We employ the lowest two excited diabatic electronic states of HeH 2 + and their electronic coupling, a coupled-channel time-dependent real wavepacket method, and a flux analysis. Both probabilities and cross sections present a few groups of resonances at low E coll, whose amplitudes decrease with the energy, due to an ion-induced dipole interaction in the entrance channel. At higher E coll, reaction probabilities and cross sections increase monotonically up to 3 eV, remaining however quite small. When H2 is in the v 0 = 1 state, the reactivity increases by ~2 orders of magnitude at the lowest energies and by ~1 order at the highest ones. Initial-state resolved rate constants at room temperature are equal to 1.74 × 10-14 and to 1.98 × 10-12 cm3s-1 at v 0 = 0 and 1, respectively. Test calculations for H2 at j 0 = 1 show that the probabilities can be enhanced by a factor of ~1/3, that is ortho-H2 seems ~4 times more reactive than para-H2.

Quantum Roaming in the Complex-Forming Mechanism of the Reactions of OH with Formaldehyde and Methanol at Low Temperature and Zero Pressure: A Ring Polymer Molecular Dynamics Approach.

The quantum dynamics of the title reactions are studied using the ring polymer molecular dynamics (RPMD) method from 20 to 1200 K using recently proposed full dimensional potential energy surfaces which include long-range dipole-dipole interactions. A V-shaped dependence of the reaction rate constants is found with a minimum at 200-300 K, in rather good agreement with the current experimental data. For temperatures above 300 K the reaction proceeds following a direct H-abstraction mechanism. However, below 100 K the reaction proceeds via organic-molecule···OH collision complexes, with very long lifetimes, longer than 10-7 s, associated with quantum roaming arising from the inclusion of quantum effects by the use of RPMD. The long lifetimes of these complexes are comparable to the time scale of the tunnelling to form reaction products. These complexes are formed at zero pressure because of quantum effects and not only at high pressure as suggested by transition state theory (TST) calculations for OH + methanol and other OH reactions. The zero-pressure rate constants reproduce quite well measured ones below 200 K, and this agreement opens the question of how important the pressure effects on the reaction rate constants are, as implied in TST-like formalisms. The zero-pressure mechanism is applicable only to very low gas density environments, such as the interstellar medium, which are not repeatable by experiments.

Origin band of the first photoionizing transition of hydrogen isocyanide.

The photoelectron spectrum of the X1Σ+ → X+2Σ+ ionizing transition of hydrogen isocyanide (HNC) is measured for the first time at a fixed photon energy (13 eV). The assignment of the spectrum is supported by wave-packet calculations simulating the photoionization transition spectrum and using ab initio calculations of the potential energy surfaces for the three lowest electronic states of the cation. The photoelectron spectrum allows the retrieval of the fundamental of the CN stretching mode of the cationic ground state ([small nu, Greek, tilde]3 = 2260 ± 80 cm-1) and the adiabatic ionization energy of hydrogen isocyanide: IE(HNC) = 12.011 ± 0.010 eV, which is far below that of HCN (IE(HCN) = 13.607 eV). In light of this latter result, the thermodynamics of the HCN+/HNC+ isomers is discussed and a short summary of the values available in the literature is given.

Low temperature reaction dynamics for CH3OH + OH collisions on a new full dimensional potential energy surface.

Is the rise of the rate constant measured in laval expansion experiments of OH with organic molecules at low temperatures due to the reaction between the reactants or due to the formation of complexes with the buffer gas? This question has importance for understanding the evolution of prebiotic molecules observed in different astrophysical objects. Among these molecules methanol is one of the most widely observed, and its reaction with OH has been studied by several groups showing a fast increase in the rate constant under 100 K. Transition state theory doesn't reproduce this behavior and here dynamical calculations are performed on a new full dimensional potential energy surface developed for this purpose. The calculated classical reactive cross sections show an increase at low collision energies due to a complex forming mechanism. However, the calculated rate constant at temperatures below 100 K remains lower than the observed one. Quantum effects are likely responsible for the measured behavior at low temperatures.

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 H2 + H3+ → H3+ + H2 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 H5+ 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.

Full dimensional potential energy surface and low temperature dynamics of the H2CO + OH → HCO + H2O reaction.

A new method is proposed to analytically represent the potential energy surface of reactions involving polyatomic molecules capable of accurately describing long-range interactions and saddle points, needed to describe low-temperature collisions. It is based on two terms, a reactive force field term and a many-body term. The reactive force field term accurately describes the fragments, long-range interactions among them and the saddle points for reactions. The many-body term increases the desired accuracy everywhere else. This method has been applied to the OH + H2CO → H2O + HCO reaction, giving a barrier of 27.4 meV. The simulated classical rate constants with this potential are in good agreement with recent experimental results [Ocaña et al., Astrophys. J., 2017, submitted], showing an important increase at temperatures below 100 K. The reaction mechanism is analyzed in detail here, and explains the observed behavior at low energy by the formation of long-lived collision complexes, with roaming trajectories, with a capture observed for very long impact parameters, >100 a.u., determined by the long-range dipole-dipole interaction.

The Photodissociation of HCN and HNC: Effects on the HNC/HCN Abundance Ratio in the Interstellar Medium.

The impact of the photodissociation of HCN and HNC isomers is analyzed in different astrophysical environments. For this purpose, the individual photodissociation cross section of HCN and HNC isomers have been calculated in the 7-13.6 eV photon energy range for a temperature of 10 K. These calculations are based on the ab initio calculation of three-dimensional adiabatic potential energy surfaces of the 21 lower electronic states. The cross sections are then obtained using a quantum wave packet calculation of the rotational transitions needed to simulate a rotational temperature of 10 K. The cross section calculated for HCN shows significant differences with respect to the experimental one, and this is attributed to the need of considering non-adiabatic transitions. Ratios between the photodissociation rates of HCN and HNC under different ultraviolet radiation fields have been computed by renormalizing the rates to the experimental one. It is found that HNC is photodissociated faster than HCN by a factor of 2.2, for the local interstellar radiation field, and 9.2, for the solar radiation field at 1 au. We conclude that to properly describe the HNC/HCN abundance ratio in astronomical environments illuminated by an intense ultraviolet radiation field it is necessary to use different photodissociation rates for each of the two isomers, obtained by integrating the product of the photodissociation cross sections and ultraviolet radiation field over the relevant wavelength range.

Photodissociation of HCN and HNC isomers in the 7-10 eV energy range.

The ultraviolet photoabsorption spectra of the HCN and HNC isomers have been simulated in the 7-10 eV photon energy range. For this purpose, the three-dimensional adiabatic potential energy surfaces of the 7 lowest electronic states, and the corresponding transition dipole moments, have been calculated, at multireference configuration interaction level. The spectra are calculated with a quantum wave packet method on these adiabatic potential energy surfaces. The spectra for the 3 lower excited states, the dissociative electronic states, correspond essentially to predissociation peaks, most of them through tunneling on the same adiabatic state. The 3 higher electronic states are bound, hereafter electronic bound states, and their spectra consist of delta lines, in the adiabatic approximation. The radiative lifetime towards the ground electronic states of these bound states has been calculated, being longer than 10 ns in all cases, much longer that the characteristic predissociation lifetimes. The spectra of HCN is compared with the available experimental and previous theoretical simulations, while in the case of HNC there are no previous studies to our knowledge. The spectrum for HNC is considerably more intense than that of HCN in the 7-10 eV photon energy range, which points to a higher photodissociation rate for HNC, compared to HCN, in astrophysical environments illuminated by ultraviolet radiation.

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.

Density-Difference-Driven Optimized Embedding Potential Method To Study the Spectroscopy of Br2 in Water Clusters.

A variant of the density difference driven optimized embedding potential (DDD-OEP) method, proposed by Roncero et al. (J. Chem. Phys. 2009, 131, 234110), has been applied to the calculation of excited states of Br2 within small water clusters. It is found that the strong interaction of Br2 with the lone electronic pair of the water molecules makes necessary to optimize specific embedding potentials for ground and excited electronic states, separately and using the corresponding densities. Diagnosis and convergence studies are presented with the aim of providing methods to be applied for the study of chromophores in solution. Also, some preliminary results obtained for the study of electronic states of Br2 in clathrate cages are presented.

A comparative study of the Au + H2, Au⁺ + H2, and Au⁻ + H2 systems: Potential energy surfaces and dynamics of reactive collisions.

In order to study the Au(-) + H2 collision, a new global potential energy surface (PES) describing the ground electronic state of AuH2(-) system is developed and compared with the PESs of the neutral [Zanchet et al., J. Chem. Phys. 132, 034301 (2010)] and cationic systems [Anaís et al., J. Chem. Phys. 135, 091102 (2011)]. We found that Au(-) - H2 presents a H-Au-H insertion minimum attributed to the stabilization of the LUMO 3b2 orbital, which can be considered as the preamble of the chemisorption well appearing in larger gold clusters. While the LUMO orbital is stabilized, the HOMO 6a1 is destabilized, creating a barrier at the geometry where the energy orbitals' curves are crossing. In the anion, this HOMO is doubly occupied, while in the neutral system is half-filled and completely empty in the cation, explaining the gradual disappearance of the well and the barrier as the number of electrons decreases. The cation presents a well in the entrance channel partially explained by electrostatic interactions. The three systems' reactions are highly endothermic, by 1.66, 2.79, and 3.23 eV for AuH, AuH(+), and AuH(-) products, respectively. The reaction dynamics is studied using quasi-classical trajectory method for the three systems. The one corresponding to the anionic system is new in this work. Collision energies between 1.00 and 8.00 eV, measured for the cation, are in good agreement with the simulated cross section for the AuH(+). It was also found that the total fragmentation, in three atoms, competes becoming dominant at sufficiently high energy. Here, we study the competition between the two different reaction pathways for the anionic, cationic, and neutral species, explaining the differences using a simple model based on the topology of the potential energy surfaces.

Full dimensional potential energy surface for the ground state of H4(+) system based on triatomic-in-molecules formalism.

In this work, we present a global potential energy surface for the ground electronic state of the H4(+) based on ab initio calculations. The final fit is based on triatomics-in-molecules (TRIM) approximation and it includes extra four-body terms for the better description of some discrepancies found on the TRIM model. The TRIM method itself allows a very accurate description of the asymptotic regions. The global fit uses more than 19,000 multireference configuration interaction ab initio points. The global potential energy surface has an overall root mean square error of 0.013 eV for energies up to 2 eV above the global minimum. This work presents an analysis of the stationary points, reactant and product channels, and crossing between the two lowest TRIM adiabatic states. It is as well included a brief description of the two first excited states of the TRIM matrix, concluding that TRIM method is a very good approximation not only for the ground state but also for at least two of the excited states of H4(+) system.