Search

Inelastic collisions and elementary chemical reactions proceeding through the formation and subsequent decay of an intermediate collision complex, with an associated deep well on the potential energy surface, pose a challenge for accurate fully quantum mechanical approaches, such as the close-coupling method. In this study, we report on the theoretical prediction of temperature-dependent state-to-state rate coefficients for these complex-mode processes, using a statistical quantum method. This statistical adiabatic channel model is benchmarked by a direct comparison using accurate rate coefficients from the literature for a number of systems (H2 + H+, HD + H+, SH+ + H, and CH+ + H) of interest in astrochemistry and astrophysics. For all of the systems considered, an error of less than factor 2 was found, at least for the dominant transitions and at low temperatures, which is sufficiently accurate for applications in the above mentioned disciplines. [ABSTRACT FROM AUTHOR]

We present a full-dimensional ab initio intermolecular potential energy surface (IPES) for the OC–HF van der Waals complex. 3167 ab initio points were computed at the frozen-core (FC) explicitly correlated coupled cluster [FC-CCSD(T)-F12b] level, with the augmented correlation-consistent polarized valence quadruple-zeta basis set plus bond functions. Basis set superposition error correction was also considered by the full counterpoise procedure. Gaussian process regression (GPR) was used to map out the potential energy surface, while a multipole expansion method was employed to smooth the ab initio noise of intermolecular potential in the long range. The global minimum of −1248.364 cm−1 was located at the linear configuration with the C atom pointing toward the H atom of the HF molecule. In addition, a local minimum of −602.026 cm−1 was found at another linear configuration with the O atom pointing toward the H atom of the HF molecule. The eigenstates were calculated on the vibrational averaged four-dimensional IPESs with the mixed radial discrete variable representation/angular finite basis representation method and Lanczos propagation algorithm. The dissociation energy D0 was calculated to be 701.827 cm−1, well reproducing the experimental value of 732 ± 2 cm−1. The dipole moment surfaces were also fitted by GPR from 3132 ab initio points calculated using the coupled cluster method [CCSD(T)] with AVTZ basis set plus bond functions. The frequencies and relative line intensities of rovibrational transitions in the HF (DF) and CO stretching bands were further calculated and compared well with the experimental results. These results indicate the high fidelity of the new IPES. [ABSTRACT FROM AUTHOR]

Linearly polarized synchrotron radiation has been used to record polarization dependent valence shell photoelectron spectra of imidazole in the photon energy range 21–100 eV. These have allowed the photoelectron angular distributions, as characterized by the anisotropy parameter β, and the electronic state intensity branching ratios to be determined. Complementing these experimental data, theoretical photoionization partial cross sections and β-parameters have been calculated for the outer valence shell orbitals. The assignment of the structure appearing in the experimental photoelectron spectra has been guided by vertical ionization energies and spectral intensities calculated by various theoretical methods that incorporate electron correlation and orbital relaxation. Strong orbital relaxation effects have been found for the 15a′, nitrogen lone-pair orbital. The calculations also predict that configuration mixing leads to the formation of several low-lying satellite states. The vibrational structure associated with ionization out of a particular orbital has been simulated within the Franck–Condon model using harmonic vibrational modes. The adiabatic approximation appears to be valid for the X2A″ state, with the β-parameter for this state being independent of the level of vibrational excitation. However, for all the other outer valence ionic states, a disparity occurs between the observed and the simulated vibrational structure, and the measured β-parameters are at variance with the behavior expected at the level of the Franck–Condon approximation. These inconsistencies suggest that the excited electronic states may be interacting vibronically such that the nuclear dynamics occur over coupled potential energy surfaces. [ABSTRACT FROM AUTHOR]

We implement a fully quantum mechanical study of methane trapping on Ir(111), where the phonons, the molecule–surface interaction, and the molecule–phonon coupling are all computed from first-principles. We find that both the surface corrugation and the phonon coupling vary strongly with molecular orientation and that there is a "chemical" aspect to this due to the catalytic nature of the metal. For example, molecules with reactive orientations can approach close to surface sites with low barriers to dissociation. Moreover, lattice motion can lower the barrier to dissociation, leading to unusual behavior for the phonon coupling. We find good agreement with experiment and two recent classical studies if we average our potential energy surface over several orientations of the molecule. We also find reasonable agreement with a recent study of methane diffraction, although we show that diffraction does not play a major role in trapping on the smooth Ir(111) surface and that trapping obeys normal energy scaling, consistent with experiment. We show that the trapping probability can be sensitive to the temperature at both high and low incidence energies. Relaxation and desorption of trapped particles are examined. [ABSTRACT FROM AUTHOR]

An analysis of the structures, some energy related properties, and key aspects of the bonding nature of the microsolvated electron with four water molecules is presented. The study is based on an exhaustive potential energy surface scan of the ground state of (H2O)4− at the UCCSD(T)/6-311(3+,4+)G(d,p)//UMP2/6-311(3+,4+)G(d,p) level. A total of 18 structures, most of them not reported before, spanning in an energy range of 8.8 kcal mol−1 were found. The energetic stability of the clusters is dictated by the effect of the excess electron on their structures, on their partial fragmentation, and on the hydrogen bonds' framework. Vertical detachment energies depend on the number of water molecules holding the excess electron in "direct contact" to their two protons at the same time and, to a lesser extent, also depend on the hydrogen bond sequence in the rest of the structure. In general, hydrogen bonds in (H2O)4− are of closed shell character, and there are other less common interactions assisted by the excess electron. [ABSTRACT FROM AUTHOR]

Modeling protoplanetary disks and other interstellar media that are not in local thermal equilibrium require the knowledge of rovibrational transition rate coefficients of molecules in collision with helium and hydrogen. We present a computational method based on the numerically exact coupled-channel (CC) method for rotational transitions and a multi-channel distorted-wave Born approximation (MC-DWBA) for vibrational transitions to calculate state-to-state rate coefficients. We apply this method to the astrophysically important case of CO2–He collisions, using newly computed ab initio three-dimensional potential energy surfaces for CO2–He with CO2 distorted along the symmetric and asymmetric stretch (ν1 and ν3) coordinates. It is shown that the MC-DWBA method is almost as accurate as full CC calculations, but more efficient. We also made computations with the more approximate vibrational coupled-channel rotational infinite-order sudden method but found that this method strongly underestimates the vibrationally inelastic collision cross sections and rate coefficients for both CO2 modes considered. [ABSTRACT FROM AUTHOR]

Since the seminal work of Tully [J. Chem. Phys. 93, 1061 (1990)], two-level scattering models have been extensively adopted as the standard benchmark systems to assess the performance of different trajectory surface hopping methods for nonadiabatic dynamics simulations. Here, we extend the branching and phase corrections to multilevel systems and combine them with both the traditional fewest switches surface hopping (FSSH) and its variant global flux surface hopping (GFSH) algorithms. To get a comprehensive evaluation of the proposed methods, we construct a series of more challenging and diverse three-level and four-level scattering models and use exact quantum solutions as references. Encouragingly, both FSSH and GFSH with the branching and phase corrections produce excellent and nearly identical results in all investigated systems, indicating that the new surface hopping methods are robust to describe multilevel problems and the reliability is insensitive to the definition of self-consistent hopping probabilities in the adiabatic representation. Furthermore, the branching correction is found to be especially important when dealing with strongly repulsive potential energy surfaces, which are common in realistic systems, thus promising for general applications. [ABSTRACT FROM AUTHOR]

In this work, we propose a new method to calculate molecular nonradiative electronic relaxation rates based on the numerically exact time-dependent density matrix renormalization group theory. This method could go beyond the existing frameworks under the harmonic approximation (HA) of the potential energy surface (PES) so that the anharmonic effect could be considered, which is of vital importance when the electronic energy gap is much larger than the vibrational frequency. We calculate the internal conversion (IC) rates in a two-mode model with Morse potential to investigate the validity of HA. We find that HA is unsatisfactory unless only the lowest several vibrational states of the lower electronic state are involved in the transition process when the adiabatic excitation energy is relatively low. As the excitation energy increases, HA first underestimates and then overestimates the IC rates when the excited state PES shifts toward the dissociative side of the ground state PES. On the contrary, HA slightly overestimates the IC rates when the excited state PES shifts toward the repulsive side. In both cases, a higher temperature enlarges the error of HA. As a real example to demonstrate the effectiveness and scalability of the method, we calculate the IC rates of azulene from S1 to S0 on the ab initio anharmonic PES approximated by the one-mode representation. The calculated IC rates of azulene under HA are consistent with the analytically exact results. The rates on the anharmonic PES are 30%–40% higher than the rates under HA. [ABSTRACT FROM AUTHOR]

Recent studies have indicated that nitramide (NH2NO2) may be formed more plentifully in the atmosphere than previously thought, while also being a missing source of the greenhouse gas nitrous oxide (N2O) via catalyzed isomerization. To validate the importance of NH2NO2 in the Earth's atmosphere, the ground and first electronic excited states of NH2NO2 were characterized and its photochemistry was investigated using multireference and coupled cluster methods. NH2NO2 is non-planar and of singlet multiplicity in the ground state while exhibiting large out-of-plane rotation in the triplet first excited state. One-dimensional cuts of the adiabatic potential energy surface calculated using the MRCI+Q method show low-lying singlet electronic states with minima in their potential along the N–N and N–O bond coordinates. Due to vertical excitation energies in the 225–180 nm region, photochemical processes will not compete in the troposphere, causing N2O production to be the predicted major removal process of NH2NO2. In the upper atmosphere, photodissociation to form NH2NO + O (3P) is suggested to be a major photochemical removal pathway. [ABSTRACT FROM AUTHOR]

Cathodes are critical components of rechargeable batteries. Conventionally, the search for cathode materials relies on experimental trial-and-error and a traversing of existing computational/experimental databases. While these methods have led to the discovery of several commercially viable cathode materials, the chemical space explored so far is limited and many phases will have been overlooked, in particular, those that are metastable. We describe a computational framework for battery cathode exploration based on ab initio random structure searching (AIRSS), an approach that samples local minima on the potential energy surface to identify new crystal structures. We show that by delimiting the search space using a number of constraints, including chemically aware minimum interatomic separations, cell volumes, and space group symmetries, AIRSS can efficiently predict both thermodynamically stable and metastable cathode materials. Specifically, we investigate LiCoO2, LiFePO4, and LixCuyFz to demonstrate the efficiency of the method by rediscovering the known crystal structures of these cathode materials. The effect of parameters, such as minimum separations and symmetries, on the efficiency of the sampling is discussed in detail. The adaptation of the minimum interatomic distances on a species-pair basis, from low-energy optimized structures to efficiently capture the local coordination environment of atoms, is explored. A family of novel cathode materials based on the transition-metal oxalates is proposed. They demonstrate superb energy density, oxygen-redox stability, and lithium diffusion properties. This article serves both as an introduction to the computational framework and as a guide to battery cathode material discovery using AIRSS. [ABSTRACT FROM AUTHOR]

Atomistic simulations based on the first-principles of quantum mechanics are reaching unprecedented length scales. This progress is due to the growth in computational power allied with the development of new methodologies that allow the treatment of electrons and nuclei as quantum particles. In the realm of materials science, where the quest for desirable emergent properties relies increasingly on soft weakly bonded materials, such methods have become indispensable. In this Perspective, an overview of simulation methods that are applicable for large system sizes and that can capture the quantum nature of electrons and nuclei in the adiabatic approximation is given. In addition, the remaining challenges are discussed, especially regarding the inclusion of nuclear quantum effects (NQEs) beyond a harmonic or perturbative treatment, the impact of NQEs on electronic properties of weakly bonded systems, and how different first-principles potential energy surfaces can change the impact of NQEs on the atomic structure and dynamics of weakly bonded systems. [ABSTRACT FROM AUTHOR]

Ab initio molecular dynamics (AIMD) is a valuable technique for studying molecules and materials at finite temperatures where the nuclei evolve on potential energy surfaces obtained from accurate electronic structure calculations. In this work, we present an approach to running AIMD simulations on noisy intermediate-scale quantum (NISQ)-era quantum computers. The electronic energies are calculated on a quantum computer using the variational quantum eigensolver (VQE) method. Algorithms for computation of analytical gradients entirely on a quantum computer require quantum fault-tolerant hardware, which is beyond NISQ-era. Therefore, we compute the energy gradients numerically using finite differences, the Hellmann–Feynman theorem, and a correlated sampling technique. This method only requires additional classical calculations of electron integrals for each degree of freedom without any additional computations on a quantum computer beyond the initial VQE run. As a proof of concept, AIMD simulations are demonstrated for the H2 molecule on IBM quantum devices. In addition, we demonstrate the validity of the method for larger molecules using full configuration interaction wave functions. As quantum hardware and noise mitigation techniques continue to improve, the method can be utilized for studying larger molecular systems. [ABSTRACT FROM AUTHOR]

In this work, a first-principles systematic study of (Pt3Cu)n, n = 1–9, clusters was performed employing the linear combination of Gaussian-type orbital auxiliary density functional theory approach. The growth of the clusters has been achieved by increasing the previous cluster by one Pt3Cu unit at a time. To explore in detail the potential energy surface of these clusters, initial structures were obtained from Born–Oppenheimer molecular dynamics trajectories generated at different temperatures and spin multiplicities. For each cluster size, several dozens of structures were optimized without any constraints. The most stable structures were characterized by frequency analysis calculations. This study demonstrates that the obtained most stable structures prefer low spin multiplicities. To gain insight into the growing pattern of these systems, average bond lengths were calculated for the lowest stable structures. This work reveals that the Cu atoms prefer to be together and to localize inside the cluster structures. Moreover, these systems tend to form octahedra moieties in the size range of n going from 4 to 9 Pt3Cu units. Magnetic moment per atom and spin density plots were obtained for the neutral, cationic, and anionic ground state structures. Dissociation energies, ionization potential, and electron affinity were calculated, too. The dissociation energy and the electron affinity increase as the number of Pt3Cu units grows, whereas the ionization potential decreases. [ABSTRACT FROM AUTHOR]

Mn(I) α-diimine carbonyl complexes have shown promise in the development of luminescent CO release materials (photoCORMs) for diagnostic and medical applications due to their ability to balance the energy of the low-lying metal-to-ligand charge transfer (MLCT) and metal-centered (MC) states. In this work, the excited state dynamics of [Mn(im)(CO)3(phen)]+ (im = imidazole; phen = 1,10-phenanthroline) is investigated by means of wavepacket propagation on the potential energy surfaces associated with the 11 low-lying Sn singlet excited states within a vibronic coupling model in a (quasi)-diabatic representation including 16 nuclear degrees of freedom. The results show that the early time photophysics (<400 fs) is controlled by the interaction between two MC dissociative states, namely, S5 and S11, with the lowest S1–S3 MLCT bound states. In particular, the presence of S1/S5 and S2/S11 crossings within the diabatic picture along the Mn–COaxial dissociative coordinate (qMn–COaxial) favors a two-stepwise population of the dissociative states, at about 60–70 fs (S11) and 160–180 fs (S5), which reaches about 10% within 200 fs. The one-dimensional reduced densities associated with the dissociative states along qMn–COaxial as a function of time clearly point to concurrent primary processes, namely, CO release vs entrapping into the S1 and S2 potential wells of the lowest luminescent MLCT states within 400 fs, characteristics of luminescent photoCORM. [ABSTRACT FROM AUTHOR]

We present a close coupling study of the bending relaxation of H2O by collision with He, taking explicitly into account the bending–rotation coupling within the rigid-bender close-coupling method. A 4D potential energy surface is developed based on a large grid of ab initio points calculated at the coupled-cluster single double triple level of theory. The bound states energies of the He–H2O complex are computed and found to be in excellent agreement with previous theoretical calculations. The dynamics results also compare very well with the rigid-rotor results available in the Basecol database and with experimental data for both rotational transitions and bending relaxation. The bending–rotation coupling is also demonstrated to be very efficient in increasing bending relaxation when the rotational excitation of H2O increases. [ABSTRACT FROM AUTHOR]

It is well known that ultraviolet photoexcitation of iron pentacarbonyl results in rapid loss of carbonyl ligands leading to the formation of coordinatively unsaturated iron carbonyl compounds. We employ ultrafast mid-infrared transient absorption spectroscopy to probe the photodissociation dynamics of gas-phase iron pentacarbonyl following ultraviolet excitation at 265 and 199 nm. After photoexcitation at 265 nm, our results show evidence for sequential dissociation of iron pentacarbonyl to form iron tricarbonyl via a short-lived iron tetracarbonyl intermediate. Photodissociation at 199 nm results in the prompt production of Fe(CO)3 within 0.25 ps via several energetically accessible pathways. An additional 15 ps time constant extracted from the data is tentatively assigned to intersystem crossing to the triplet manifold of iron tricarbonyl or iron dicarbonyl. Mechanisms for formation of iron tetracarbonyl, iron tricarbonyl, and iron dicarbonyl are proposed and theoretically validated with one-dimensional cuts through the potential energy surface as well as bond dissociation energies. Ground state calculations are computed at the CCSD(T) level of theory and excited states are computed with EOM-EE-CCSD(dT). [ABSTRACT FROM AUTHOR]

Investigations into bimolecular reaction kinetics probe the details of the underlying potential energy surface (PES), which can help to validate high-level quantum chemical calculations. We utilize a combined linear Paul ion trap with a time-of-flight mass spectrometer to study isotopologue reactions between acetylene cations ( C 2 H 2 + ) and two isomers of C3H4: propyne (HC3H3) and allene (H2C3H2). In a previous study [Schmid et al., Phys. Chem. Chem. Phys. 22, 20303 (2020)],1 we showed that the two isomers of C3H4 have fundamentally different reaction mechanisms. Here, we further explore the calculated PES by isotope substitution. While isotopic substitution of reactants is a standard experimental tool in the investigation of molecular reaction kinetics, the controlled environment of co-trapped, laser-cooled Ca+ ions allows the different isotopic reaction pathways to be followed in greater detail. We report branching ratios for all of the primary products of the different isotopic species. The results validate the previously proposed mechanism: propyne forms a bound reaction complex with C 2 H 2 + , while allene and C 2 H 2 + perform long-range charge exchange only. [ABSTRACT FROM AUTHOR]

We present a first-principles study of the static and dynamic aspects of the strong Jahn–Teller (JT) and pseudo-JT (PJT) effects in niobium tetrafluoride, NbF4, in the manifold of its electronic ground state, 2E, and its first excited state, 2T2. The complex topography of the full-dimensional multi-sheeted adiabatic JT/PJT surfaces is analyzed computationally at the complete-active-space self-consistent-field (CASSCF) and multireference second-order perturbation levels of electronic structure theory, providing a detailed characterization of minima, saddle points, and minimum-energy conical intersection points. The calculations reveal that the tetrahedral (Td) configuration of NbF4 undergoes strong JT distortions along the bending mode of e symmetry, yielding tetragonal molecular structures of D2d symmetry with Td → D2d stabilization energies of about 2000 cm−1 in the X ̃ 2 E state and about 6400 cm−1 in the à 2 T 2 state. In addition, there exists strong X ̃ 2 E − à 2 T 2 PJT coupling via the bending mode of t2 symmetry, which becomes important near the crossing seam of the X ̃ 2 E and à 2 T 2 potential energy surfaces. A five-state five-mode JT/PJT vibronic-coupling Hamiltonian is constructed in terms of symmetry-invariant polynomial expansions of the X ̃ 2 E and à 2 T 2 diabatic potential energy surfaces in the e and t2 bending coordinates. The parameters of the Hamiltonian are determined by a least-squares fit of its eigenvalues to the CASSCF ab initio data. The vibronic spectra and the time evolution of adiabatic electronic population probabilities are computed with the multi-configuration time-dependent Hartree method. The complexity of the spectra reflects the effects of the exceptionally strong E × e and T2 × e JT couplings and (E + T2) × (e + t2) PJT coupling. The time evolution of the populations of the adiabatic electronic states after the initial preparation of the à 2 T 2 state reveals the femtosecond nonadiabatic dynamics through a multidimensional seam of conical intersection. These results represent the first study of the static and dynamical JT/PJT effects in the X ̃ 2 E and à 2 T 2 electronic states of NbF4. [ABSTRACT FROM AUTHOR]

A first-principles based quantum dynamics study of the Li + LiNa(v = 0, j = 0) → Li2( v ′ , j′) + Na reaction is reported for collision energies spanning the ultracold (1 nK) to cold (1 K) regimes. A full-dimensional ab initio potential energy surface for the ground electronic state of Li2Na is utilized that includes an accurate treatment of the long-range interactions. The Li + LiNa reaction is barrierless and exoergic and exhibits a deep attractive potential well that supports complex formation. Thus, significant reactivity occurs even for collision temperatures approaching absolute zero. The reactive scattering calculations are based on a numerically exact time-independent quantum dynamics methodology in hyperspherical coordinates. Total and rotationally resolved rate coefficients are reported at 56 collision energies and include all contributing partial waves. Several shape resonances are observed in many of the rotationally resolved rate coefficients and a small resonance feature is also reported in the total rate coefficient near 50 mK. Of particular interest, the angular distributions or differential cross sections are reported as a function of both the collision energy and scattering angle. Unique quantum fingerprints (bumps, channels, and ripples) are observed in the angular distributions for each product rotational state due to quantum interference and shape resonance contributions. The Li + LiNa reaction is under active experimental investigation so that these intriguing features could be verified experimentally when sufficient product state resolution becomes feasible for collision energies below 1 K. [ABSTRACT FROM AUTHOR]

We examine rotational transitions of HCl in collisions with H2 by carrying out quantum mechanical close-coupling and quasi-classical trajectory (QCT) calculations on a recently developed globally accurate full-dimensional ab initio potential energy surface for the H3Cl system. Signatures of rainbow scattering in rotationally inelastic collisions are found in the state resolved integral and differential cross sections as functions of the impact parameter (initial orbital angular momentum) and final rotational quantum number. We show the coexistence of distinct dynamical regimes for the HCl rotational transition driven by the short-range repulsive and long-range attractive forces whose relative importance depends on the collision energy and final rotational state, suggesting that the classification of rainbow scattering into rotational and l-type rainbows is effective for H2 + HCl collisions. While the QCT method satisfactorily predicts the overall behavior of the rotationally inelastic cross sections, its capability to accurately describe signatures of rainbow scattering appears to be limited for the present system. [ABSTRACT FROM AUTHOR]

Multi-configurational time-dependent Hartree (MCTDH) calculations using time-dependent grid representations can be used to accurately simulate high-dimensional quantum dynamics on general ab initio potential energy surfaces. Employing the correlation discrete variable representation, sets of direct product type grids are employed in the calculation of the required potential energy matrix elements. This direct product structure can be a problem if the coordinate system includes polar and azimuthal angles that result in singularities in the kinetic energy operator. In the present work, a new direct product-type discrete variable representation (DVR) for arbitrary sets of polar and azimuthal angles is introduced. It employs an extended coordinate space where the range of the polar angles is taken to be [−π, π]. The resulting extended space DVR resolves problems caused by the singularities in the kinetic energy operator without generating a very large spectral width. MCTDH calculations studying the F·CH4 complex are used to investigate important properties of the new scheme. The scheme is found to allow for more efficient integration of the equations of motion compared to the previously employed cot-DVR approach [G. Schiffel and U. Manthe, Chem. Phys. 374, 118 (2010)] and decreases the required central processing unit times by about an order of magnitude. [ABSTRACT FROM AUTHOR]

A new 3D-potential energy surface (3D-PES) for the weakly bound CH3Cl–He complex is mapped in Jacobi coordinates. Electronic structure calculations are performed using the explicitly correlated coupled clusters with single, double, and perturbative triple excitations approach in conjunction with the aug-cc-pVTZ basis set. Then, an analytical expansion of this 3D-PES is derived. This PES shows three minimal structures for collinear C–Cl–He arrangements and for He located in between two H atoms, in the plane parallel to the three H atoms, which is near the center of mass of CH3Cl. The latter form corresponds to the global minimum. Two maxima are also found, which connect the minimal structures. We then evaluated the pressure broadening coefficients of the spectral lines of CH3Cl in a helium bath based on our ab initio potential. Satisfactory agreement with experiments was observed, confirming the good accuracy of our 3D-PES. We also derived the bound rovibronic levels for ortho- and para-CH3Cl–He dimers after quantum treatment of the nuclear motions. For both clusters, computations show that although the ground vibrational state is located well above the intramolecular isomerization barriers, the rovibronic levels may be associated with a specific minimal structure. This can be explained by vibrational localization and vibrational memory effects. [ABSTRACT FROM AUTHOR]

We evaluate the effect of electronic decoherence on intersystem crossing in the photodynamics of thioformaldehyde. First, we show that the state-averaged complete-active-space self-consistent field electronic structure calculations with a properly chosen active space of 12 active electrons in 10 active orbitals can predict the potential energy surfaces and the singlet–triplet spin–orbit couplings quite well for CH2S, and we use this method for direct dynamics by coherent switching with decay of mixing (CSDM). We obtain similar dynamical results with CSDM or by adding energy-based decoherence to trajectory surface hopping, with the population of triplet states tending to a small steady-state value over 500 fs. Without decoherence, the state populations calculated by the conventional trajectory surface hopping method or the semiclassical Ehrenfest method gradually increase. This difference shows that decoherence changes the nature of the results not just quantitatively but qualitatively. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

We employ theoretically "exact" and numerically "accurate" Beyond Born–Oppenheimer (BBO) treatment to construct diabatic potential energy surfaces (PESs) of the benzene radical cation ( C 6 H 6 + ) for the first time and explore the workability of the time-dependent discrete variable representation (TDDVR) method for carrying out dynamical calculations to evaluate the photoelectron (PE) spectra of its neutral analog. Ab initio adiabatic PESs and nonadiabatic coupling terms are computed over a series of pairwise normal modes, which exhibit rich nonadiabatic interactions starting from Jahn–Teller interactions and accidental conical intersections/seams to pseudo Jahn–Teller couplings. Once the electronic structure calculation is completed on the low-lying five doublet electronic states ( X ̃ 2 E 1 g , B ̃ 2 E 2 g , and C ̃ 2 A 2 u ) of the cationic species, diabatization is carried out employing the adiabatic-to-diabatic transformation (ADT) equations for the five-state sub-Hilbert space to compute highly accurate ADT angles, and thereby, single-valued, smooth, symmetric, and continuous diabatic PESs and couplings are constructed. Subsequently, such surface matrices are used to perform multi-state multi-mode nuclear dynamics for simulating PE spectra of benzene. Our theoretical findings clearly depict that the spectra for X ̃ 2 E 1 g and B ̃ 2 E 2 g − C ̃ 2 A 2 u states obtained from BBO treatment and TDDVR dynamics exhibit reasonably good agreement with the experimental results as well as with the findings of other theoretical approaches. [ABSTRACT FROM AUTHOR]

Valence photoelectron spectra and photoelectron angular distributions of trans-dichloroethene have been measured with vibrational resolution at photon energies between 19 eV and 90 eV. Calculations of photoelectron anisotropy parameters, β, and harmonic vibrational modes help provide initial insight into the molecular structure. The photon energy range encompasses the expected position of the atomic Cl 3p Cooper minimum. A corresponding dip observed here in the anisotropy of certain photoelectron bands permits the identification and characterization of those molecular orbitals that retain a localized atomic Cl character. The adiabatic approximation holds for the X2Au state photoelectron band, but vibronic coupling was inferred within the A–B–C and the D–E states by noting various failures of the Franck–Condon model, including vibrationally dependent β-parameters. This is further explored using the linear vibronic coupling model with interaction parameters obtained from ab initio calculations. The A/B photoelectron band is appreciably affected by vibronic coupling, owing to the low-lying conical intersection of the A2Ag and B2Bu states. The C2Bg band is also affected, but to a lesser extent. The adiabatic minima of the D2Au and E2Ag states are almost degenerate, and the vibronic interaction between these states is considerable. The potential energy surface of the D2Au state is predicted to have a double-minimum shape with respect to the au deformations of the molecular structure. The irregular vibrational structure of the resulting single photoelectron band reflects the non-adiabatic nuclear dynamics occurring on the two coupled potential energy surfaces above the energy of their conical intersection. [ABSTRACT FROM AUTHOR]

We address topics related to molecules coupled to quantum radiation. The formalism of light–matter interaction is different for classical and quantum fields, but some analogies remain, such as the formation of light induced crossings. We show that under particular circumstances, the molecular dynamics under quantum or classical fields produce similar results, as long as the radiation is prepared as a Fock state and far from ultra-strong coupling regimes. At this point, the choice of specific initial Fock states is irrelevant since the dynamics scales. However, in realistic multistate molecular systems, radiative scaling may fail due to the presence of simultaneous efficient non-radiative couplings in the dynamics. Polar molecules have permanent dipoles, and within the context of the full quantum Rabi model with a Pauli–Fierz Hamiltonian, they play a crucial role in the polaritonic dynamics since both permanent dipole moments and self-energy terms produce drastic changes on the undressed potential energy surfaces at high coupling strengths. We also gauge the effect of including rotational degrees of freedom in cavity molecular photodynamics. For diatomic molecules, the addition of rotation amounts to transform (both with classical or quantum fields) a light induced crossing into a light induced conical intersection. However, we show that conical intersections due to molecular rotation do not represent the standard properties of well-known efficient intrinsic conical intersections inasmuch they do not enhance the quantum transition rates. [ABSTRACT FROM AUTHOR]

Global coupled three-state two-channel potential energy and property/interaction (dipole and spin–orbit coupling) surfaces for the dissociation of NH3(Ã) into NH + H2 and NH2 + H are reported. The permutational invariant polynomial-neural network approach is used to simultaneously fit and diabatize the electronic Hamiltonian by fitting the energies, energy gradients, and derivative couplings of the two coupled lowest-lying singlet states as well as fitting the energy and energy gradients of the lowest-lying triplet state. The key issue in fitting property matrix elements in the diabatic basis is that the diabatic surfaces must be smooth, that is, the diabatization must remove spikes in the original adiabatic property surfaces attributable to the switch of electronic wavefunctions at the conical intersection seam. Here, we employ the fit potential energy matrix to transform properties in the adiabatic representation to a quasi-diabatic representation and remove the discontinuity near the conical intersection seam. The property matrix elements can then be fit with smooth neural network functions. The coupled potential energy surfaces along with the dipole and spin–orbit coupling surfaces will enable more accurate and complete treatment of optical transitions, as well as nonadiabatic internal conversion and intersystem crossing. [ABSTRACT FROM AUTHOR]

We present a semiclassically approximate quantum treatment of solvation with the purpose of investigating the accuracy of the Caldeira–Leggett model. We do that by simulating the vibrational features of water solvation by means of two different approaches. One is entirely based on the adoption of an accurate ab initio potential to describe water clusters of increasing dimensionality. The other one consists of a model made of a central water molecule coupled to a high-dimensional Caldeira–Leggett harmonic bath. We demonstrate the role of quantum effects in the detection of water solvation and show that the computationally cheap approach based on the Caldeira–Leggett bath is only partially effective. The main conclusion of the study is that quantum methods associated with high-level potential energy surfaces are necessary to correctly study solvation features, while simplified models, even if attractive owing to their reduced computational cost, can provide some useful insights but are not able to come up with a comprehensive description of the solvation phenomenon. [ABSTRACT FROM AUTHOR]

The vibrational quenching cross sections and corresponding low-temperature rate constants for the ν = 1 and ν = 2 states of CN−(1Σ+) colliding with He and Ar atoms have been computed ab initio using new three-dimensional potential energy surfaces. Little work has been carried out so far on low-energy vibrationally inelastic collisions for anions with neutral atoms. The cross sections and rates calculated at energies and temperatures relevant for both ion traps and astrochemical modeling are found by the present calculations to be even smaller than those of the similar C 2 − /He and C 2 − /Ar systems, which are in turn of the order of those existing for the collisions involving neutral diatom–atom systems. The implications of our finding in the present case mainly focus on the possible role of small computed rate constants in the dynamics of molecular cooling and the evolution of astrochemical modeling networks. [ABSTRACT FROM AUTHOR]

Potential energy surfaces for high-energy collisions between an oxygen molecule and a nitrogen atom are useful for modeling chemical dynamics in shock waves. In the present work, we present doublet, quartet, and sextet potential energy surfaces that are suitable for studying collisions of O2( 3 Σ g − ) with N(4S) in the electronically adiabatic approximation. Two sets of surfaces are developed, one using neural networks (NNs) with permutationally invariant polynomials (PIPs) and one with the least-squares many-body (MB) method, where a two-body part is an accurate diatomic potential and the three-body part is expressed with connected PIPs in mixed-exponential-Gaussian bond order variables (MEGs). We find, using the same dataset for both fits, that the fitting performance of the PIP-NN method is significantly better than that of the MB-PIP-MEG method, even though the MB-PIP-MEG fit uses a higher-order PIP than those used in previous MB-PIP-MEG fits of related systems (such as N4 and N2O2). However, the evaluation of the PIP-NN fit in trajectory calculations requires about 5 times more computer time than is required for the MB-PIP-MEG fit. [ABSTRACT FROM AUTHOR]

Minor structural modifications to the DNA and RNA nucleobases have a significant effect on their excited state dynamics and electronic relaxation pathways. In this study, the excited state dynamics of 7-deazaguanosine and guanosine 5′-monophosphate are investigated in aqueous solution and in a mixture of methanol and water using femtosecond broadband transient absorption spectroscopy following excitation at 267 nm. The transient spectra are collected using photon densities that ensure no parasitic multiphoton-induced signal from solvated electrons. The data can be fit satisfactorily using a two- or three-component kinetic model. By analyzing the results from steady-state, time-resolved, computational calculations, and the methanol–water mixture, the following general relaxation mechanism is proposed for both molecules, Lb → La → 1πσ*(ICT) → S0, where the 1πσ*(ICT) stands for an intramolecular charge transfer excited singlet state with significant πσ* character. In general, longer lifetimes for internal conversion are obtained for 7-deazaguanosine compared to guanosine 5′-monophosphate. Internal conversion of the 1πσ*(ICT) state to the ground state occurs on a similar time scale of a few picoseconds in both molecules. Collectively, the results demonstrate that substitution of a single nitrogen atom for a methine (C–H) group at position seven of the guanine moiety stabilizes the 1ππ* Lb and La states and alters the topology of their potential energy surfaces in such a way that the relaxation dynamics in 7-deazaguanosine are slowed down compared to those in guanosine 5′-monophosphate but not for the internal conversion of 1πσ*(ICT) state to the ground state. [ABSTRACT FROM AUTHOR]

The gas-phase reaction of CCl+ with acetonitrile (CH3CN) is studied using a linear Paul ion trap coupled to a time-of-flight mass spectrometer. This work builds on a previous study of the reaction of CCl+ with acetylene [K. J. Catani et al., J. Chem. Phys. 152, 234310 (2020)] and further explores the reactivity of CCl+ with organic neutral molecules. Both of the reactant species are relevant in observations and models of chemistry in the interstellar medium. Nitriles, in particular, are noted for their relevance in prebiotic chemistry and are found in the atmosphere of Titan, one of Saturn's moons. This work represents one of the first studied reactions of a halogenated carbocation with a nitrile and the first exploration of CCl+ with a nitrile. Reactant isotopologues are used to unambiguously assign ionic primary products from this reaction: HNCCl+ and C2H3+. Branching ratios are measured, and both primary products are determined to be equally probable. Quantum chemical and statistical reaction rate theory calculations illuminate pertinent information for interpreting the reaction data, including reaction thermodynamics and a potential energy surface for the reaction, as well as rate constants and branching ratios for the observed products. In particular, the reaction products and potential energy surface stimulate questions regarding the strength and role of the nitrile functional group, which can be further explored with more reactions of this class. [ABSTRACT FROM AUTHOR]

We introduce vibrational heat-bath configuration interaction (VHCI) as an accurate and efficient method for calculating vibrational eigenstates of anharmonic systems. Inspired by its origin in electronic structure theory, VHCI is a selected CI approach that uses a simple criterion to identify important basis states with a pre-sorted list of anharmonic force constants. Screened second-order perturbation theory and simple extrapolation techniques provide significant improvements to variational energy estimates. We benchmark VHCI on four molecules with 12–48 degrees of freedom and use anharmonic potential energy surfaces truncated at fourth and sixth orders. When compared to other methods using the same truncated potentials, VHCI produces vibrational spectra of tens or hundreds of states with sub-wavenumber accuracy at low computational cost. [ABSTRACT FROM AUTHOR]

A new spin-averaged potential energy surface (PES) for non-reactive O2( Σ g − 3 ) + O2( Σ g − 3 ) collisions is presented. The potential is formulated analytically according to the nature of the principal interaction components, with the main van der Waals contribution described through the improved Lennard-Jones model. All the parameters involved in the formulation, having a physical meaning, have been modulated in restricted variation ranges, exploiting a combined analysis of experimental and ab initio reference data. The new PES is shown to be able to reproduce a wealth of different physical properties, ranging from the second virial coefficients to transport properties (shear viscosity and thermal conductivity) and rate coefficients for inelastic scattering collisions. Rate coefficients for the vibrational inelastic processes of O2, including both vibration-to-vibration (V–V) and vibration-to-translation/rotation (V–T/R) energy exchanges, were then calculated on this PES using a mixed quantum–classical method. The effective formulation of the potential and its combination with an efficient, yet accurate, nuclear dynamics treatment allowed for the determination of a large database of V–V and V–T/R energy transfer rate coefficients in a wide temperature range. [ABSTRACT FROM AUTHOR]

Revisiting recently published Raman jet spectra of monomeric formic acid with accurate high order perturbative calculations based on two explicitly correlated coupled-cluster quality potential energy surfaces from the literature, we assign and add 11 new vibrational band centers to the trans-HCOOH database and 53 for its three deuterated isotopologs. Profiting from the synergy between accurate calculations and symmetry information from depolarized Raman spectra, we reassign eight literature IR bands up to 4000 cm−1. Experimental detection of highly excited torsional states (ν9) of trans-HCOOH, such as 4ν9 and ν6 + 2ν9, reveals substantial involvement of the C–O stretch ν6 into the O–H bend/torsion resonance ν5/2ν9, which is part of a larger resonance polyad. Depolarization and isotopic C-D substitution experiments further elucidate the nature of Raman peaks in the vicinity of the O–H stretching fundamental (ν1), which seem to be members of a large set of interacting states that can be identified and described with a polyad quantum number and that gain intensity via resonance mixing with ν1. [ABSTRACT FROM AUTHOR]

We report fully quantum calculations of the collisional perturbation of a molecular line for a system that is relevant for Earth's atmosphere. We consider the N2-perturbed pure rotational R(0) line in CO. The results agree well with the available experimental data. This work constitutes a significant step toward populating the spectroscopic databases with ab initio collisional line-shape parameters for atmosphere-relevant systems. The calculations were performed using three different recently reported potential energy surfaces (PESs). We conclude that all three PESs lead to practically the same values of the pressure broadening coefficients. [ABSTRACT FROM AUTHOR]

Nuclear dynamics in the first six vibronically coupled electronic states of pentafluorobenzene radical cation is studied with the aid of the standard vibronic coupling theory and quantum dynamical methods. A model 6 × 6 vibronic Hamiltonian is constructed in a diabatic electronic basis using symmetry selection rules and a Taylor expansion of the elements of the electronic Hamiltonian in terms of the normal coordinate of vibrational modes. Extensive ab initio quantum chemistry calculations are carried out for the adiabatic electronic energies to establish the diabatic potential energy surfaces and their coupling surfaces. Both time-independent and time-dependent quantum mechanical methods are employed to perform nuclear dynamics calculations. The vibronic spectrum of the electronic states is calculated, assigned, and compared with the available experimental results. Internal conversion dynamics of electronic states is examined to assess the impact of various couplings on the nuclear dynamics. The impact of increasing fluorination of the parent benzene radical cation on its radiative emission is examined and discussed. [ABSTRACT FROM AUTHOR]

Two different ab initio potential energy surfaces are employed to investigate the efficiency of the rotational excitation channels for the polar molecular ion HeH+ interacting with He atoms. We further use them to investigate the quantum dynamics of both the proton-exchange reaction and the purely rotational inelastic collisions over a broad range of temperatures. In current modeling studies, this cation is considered to be one of the possible cooling sources under early universe conditions after the recombination era and has recently been found to exist in the interstellar medium. The results from the present calculations are able to show the large efficiency of the state-changing channels involving rotational states of this cation. In fact, we find them to be similar in size and behavior to the inelastic and reaction rate coefficients obtained in previous studies, where H atoms were employed as projectiles. The same rotational excitation processes, occurring when free electrons are the collision partners of this cation, are also compared with the present findings. The relative importance of the reactive, proton-exchange channel and the purely inelastic channels is also analyzed and discussed. The rotational de-excitation processes are also investigated for the cooling kinetics of the present cation under cold trap conditions with He as the buffer gas. The implications of the present results for setting up more comprehensive numerical models to describe the chemical evolution networks in different environments are briefly discussed. [ABSTRACT FROM AUTHOR]

Phosphorus is of particular interest in astrochemistry because it is a biogenic element together with hydrogen, carbon, nitrogen, oxygen, and sulfur. However, the chemical evolution of such element in the interstellar medium (ISM) is still far from an accurate characterization, with the chemistry of P-bearing molecules being poorly understood. To provide a contribution in this direction, we have carried out an accurate investigation of the potential energy surface for the reaction between the CP radical and methanimine (CH2NH), two species already detected in the ISM. In analogy to similar systems, i.e., CH2NH + X, with X = OH, CN, and CCH, this reaction can occur—from an energetic point of view—under the harsh conditions of the ISM. Furthermore, since the major products of the aforementioned reaction, namely, E- and Z-2-phosphanylidyneethan-1-imine (HN=CHCP) and N-(phosphaneylidynemethyl)methanimine (H2C=NCP), have not been spectroscopically characterized yet, some effort has been made for filling this gap by means of accurate computational approaches. [ABSTRACT FROM AUTHOR]

A new analytical potential energy surface (PES) has been constructed for the Ar2H+ system from a dataset consisting of a large number of ab initio energies computed using the coupled-cluster singles, doubles and perturbative triples method and aug-cc-pVQZ basis set. The long-range interaction is added to the diatomic potentials using a standard long range expansion form to better describe the asymptotic regions. The vibrational states for the most stable structures of the Ar2H+ system have been calculated, and few low lying states are assigned to quantum numbers. Reactive scattering studies have been performed for the Ar + Ar′H+ → Ar′ + ArH+ proton exchange reaction on the newly generated PES. Reaction probability, cross sections, and rate constants are calculated for the Ar + Ar′H+(v = 0, j = 0) collisions within 0.01 eV–0.6 eV of relative translational energy using exact quantum dynamical simulations as well as quasiclassical trajectory (QCT) calculations. The effect of vibrational excitation of the reactants is also explored for the reaction. State averaged rate constants are calculated for the proton exchange reaction at different temperatures using the QCT method. The mechanistic pathways for the reaction are understood by analyzing the quasiclassical trajectories. [ABSTRACT FROM AUTHOR]

The nuclear–electronic orbital (NEO) method is a multicomponent quantum chemistry theory that describes electronic and nuclear quantum effects simultaneously while avoiding the Born–Oppenheimer approximation for certain nuclei. Typically specified hydrogen nuclei are treated quantum mechanically at the same level as the electrons, and the NEO potential energy surface depends on the classical nuclear coordinates. This approach includes nuclear quantum effects such as zero-point energy and nuclear delocalization directly into the potential energy surface. An extended NEO potential energy surface depending on the expectation values of the quantum nuclei incorporates coupling between the quantum and classical nuclei. Herein, theoretical methodology is developed to optimize and characterize stationary points on the standard or extended NEO potential energy surface, to generate the NEO minimum energy path from a transition state down to the corresponding reactant and product, and to compute thermochemical properties. For this purpose, the analytic coordinate Hessian is developed and implemented at the NEO Hartree–Fock level of theory. These NEO Hessians are used to study the SN2 reaction of ClCH3Cl− and the hydride transfer of C4H9+. For each system, analysis of the single imaginary mode at the transition state and the intrinsic reaction coordinate along the minimum energy path identifies the dominant nuclear motions driving the chemical reaction. Visualization of the electronic and protonic orbitals along the minimum energy path illustrates the coupled electronic and protonic motions beyond the Born–Oppenheimer approximation. This work provides the foundation for applying the NEO approach at various correlated levels of theory to a wide range of chemical reactions. [ABSTRACT FROM AUTHOR]

Extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) [A. M. N. Niklasson, Phys. Rev. Lett. 100, 123004 (2008)] is formulated for orbital-free Hohenberg–Kohn density-functional theory and for charge equilibration and polarizable force-field models that can be derived from the same orbital-free framework. The purpose is to introduce the most recent features of orbital-based XL-BOMD to molecular dynamics simulations based on charge equilibration and polarizable force-field models. These features include a metric tensor generalization of the extended harmonic potential, preconditioners, and the ability to use only a single Coulomb summation to determine the fully equilibrated charges and the interatomic forces in each time step for the shadow Born–Oppenheimer potential energy surface. The orbital-free formulation has a charge-dependent, short-range energy term that is separate from long-range Coulomb interactions. This enables local parameterizations of the short-range energy term, while the long-range electrostatic interactions can be treated separately. The theory is illustrated for molecular dynamics simulations of an atomistic system described by a charge equilibration model with periodic boundary conditions. The system of linear equations that determines the equilibrated charges and the forces is diagonal, and only a single Ewald summation is needed in each time step. The simulations exhibit the same features in accuracy, convergence, and stability as are expected from orbital-based XL-BOMD. [ABSTRACT FROM AUTHOR]

"Δ-machine learning" refers to a machine learning approach to bring a property such as a potential energy surface (PES) based on low-level (LL) density functional theory (DFT) energies and gradients close to a coupled cluster (CC) level of accuracy. Here, we present such an approach that uses the permutationally invariant polynomial (PIP) method to fit high-dimensional PESs. The approach is represented by a simple equation, in obvious notation VLL→CC = VLL + ΔVCC–LL, and demonstrated for CH4, H3O+, and trans and cis-N-methyl acetamide (NMA), CH3CONHCH3. For these molecules, the LL PES, VLL, is a PIP fit to DFT/B3LYP/6-31+G(d) energies and gradients and ΔVCC–LL is a precise PIP fit obtained using a low-order PIP basis set and based on a relatively small number of CCSD(T) energies. For CH4, these are new calculations adopting an aug-cc-pVDZ basis, for H3O+, previous CCSD(T)-F12/aug-cc-pVQZ energies are used, while for NMA, new CCSD(T)-F12/aug-cc-pVDZ calculations are performed. With as few as 200 CCSD(T) energies, the new PESs are in excellent agreement with benchmark CCSD(T) results for the small molecules, and for 12-atom NMA, training is done with 4696 CCSD(T) energies. [ABSTRACT FROM AUTHOR]

Accommodation and migration of the ground-state (2s22p4 3P) oxygen atom in the ideal Ar, Kr, and Xe rare gas crystals are investigated using the classical model. The model accounts for anisotropy of interaction between guest and host atoms, spin–orbit coupling, and lattice relaxation. Interstitial and substitutional accommodations are found to be the only thermodynamically stable sites for trapping atomic oxygen. Mixing of electronic states coupled to lattice distortions justifies that its long-range thermal migration follows the adiabatic ground-state potential energy surface. Search for the migration paths reveals a common direct mechanism for interstitial diffusion. Substitutional atoms are activated by the point lattice defects, whereas the direct guest–host exchange meets a higher activation barrier. These three low-energy migration mechanisms provide plausible interpretation for multiple migration activation thresholds observed in Kr and Xe free-standing crystals, confirmed by reasonable agreement between calculated and measured activation energies. An important effect of interaction anisotropy and a minor role of spin–orbit coupling are emphasized. [ABSTRACT FROM AUTHOR]

Rotational excitation of interstellar PN molecules induced by collisions with H2 is investigated. We present the first ab initio four-dimensional potential energy surface (PES) for the PN–H2 van der Waals system. The PES was obtained using an explicitly correlated coupled cluster approach with single, double, and perturbative triple excitations [CCSD(T)-F12b]. The method of interpolating moving least squares was used to construct an analytical PES from these data. The equilibrium structure of the complex was found to be linear, with H2 aligned at the N end of the PN molecule, at an intermolecular separation of 4.2 Å. The corresponding well-depth is 224.3 cm−1. The dissociation energies were found to be 40.19 cm−1 and 75.05 cm−1 for complexes of PN with ortho-H2 and para-H2, respectively. Integral cross sections for rotational excitation in PN–H2 collisions were calculated using the new PES and were found to be strongly dependent on the rotational level of the H2 molecule. These new collisional data will be crucial to improve the estimation of PN abundance in the interstellar medium from observational spectra. [ABSTRACT FROM AUTHOR]

Cross sections and rate coefficients for sulfur monoxide (SO) + H2 collisions are calculated using a full six-dimensional (6D) potential energy surface (PES). The coupled states (CS) approximation is used to compute fine-structure resolved cross sections for rovibrational transitions between states with v = 0–2, where v is the vibrational quantum number of the SO molecule. The CS calculations for Δv = 1 are benchmarked against close-coupling (CC) results for spin-free interactions. For Δv = 0, the present fine-structure resolved CS results are benchmarked against existing CC results obtained with a rigid rotor approximation. In both cases, the agreement is found to be satisfactory, which suggests that the present results may provide reliable estimates for fine-structure resolved rovibrational transitions. These estimates are the first of their kind based on a full 6D PES. Rate coefficients are reported for temperatures between 10 K and 3000 K for both para- and ortho-H2 colliders. A comparison of the para-H2 rates with mass-scaled results for He shows substantial differences that may be important in astrophysical models. [ABSTRACT FROM AUTHOR]

The collisions transferring large portions of energy are often called supercollisions. In the H + C2H2 reactive system, the rovibrationally cold C2H2 molecule can be activated with substantial internal excitations by its collision with a translationally hot H atom. It is interesting to investigate the mechanisms of collisional energy transfer in other important reactions of H with hydrocarbons. Here, an accurate, global, full-dimensional potential energy surface (PES) of H + C2H4 was constructed by the fundamental invariant neural network fitting based on roughly 100 000 UCCSD(T)-F12a/aug-cc-pVTZ data points. Extensive quasi-classical trajectory calculations were carried out on the full-dimensional PES to investigate the energy transfer process in collisions of the translationally hot H atoms with C2H4 in a wide range of collision energies. The computed function of the energy-transfer probability is not a simple exponential decay function but exhibits large magnitudes in the region of a large amount of energy transfer, indicating the signature of supercollisions. The supercollisions among non-complex-forming nonreactive (prompt) trajectories are frustrated complex-forming processes in which the incoming H atom penetrates into C2H4 with a small C–H distance but promptly and directly leaves C2H4. The complex-forming supercollisions, in which either the attacking H atom leaves (complex-forming nonreactive collisions) or one of the original H atoms of C2H4 leaves (complex-forming reactive trajectories), dominate large energy transfer from the translational energy to internal excitation of molecule. The current work sheds valuable light on the energy transfer of this important reaction in the combustion and may motivate related experimental investigations. [ABSTRACT FROM AUTHOR]

The exact time-dependent potential energy surface driving the nuclear dynamics was recently shown to be a useful tool to understand and interpret the coupling of nuclei, electrons, and photons in cavity settings. Here, we provide a detailed analysis of its structure for exactly solvable systems that model two phenomena: cavity-induced suppression of proton-coupled electron-transfer and its dependence on the initial state, and cavity-induced electronic excitation. We demonstrate the inadequacy of simply using a weighted average of polaritonic surfaces to determine the dynamics. Such a weighted average misses a crucial term that redistributes energy between the nuclear and the polaritonic systems, and this term can in fact become a predominant term in determining the nuclear dynamics when several polaritonic surfaces are involved. Evolving an ensemble of classical trajectories on the exact potential energy surface reproduces the nuclear wavepacket quite accurately, while evolving on the weighted polaritonic surface fails after a short period of time. The implications and prospects for application of mixed quantum-classical methods based on this surface are discussed. [ABSTRACT FROM AUTHOR]