Ab initio potential energy surface and infrared spectra of H2-CO and D2-CO van der Waals complexes

Ultrafast photoinduced C-H bond formation from two small inorganic molecules

The formation of carbon-hydrogen (C-H) bonds via the reaction of small inorganic molecules is of great significance for understanding the fundamental transition from inorganic to organic matter, and thus the origin of life. Yet, the detailed mechanism of the C-H bond formation, particularly the time scale and molecular-level control of the dynamics, remain elusive. Here, we investigate the light-induced bimolecular reaction starting from a van der Waals molecular dimer composed of two small inorganic molecules, H2 and CO. Employing reaction microscopy driven by a tailored two-color light field, we identify the pathways leading to C-H photobonding thereby producing HCO+ ions, and achieve coherent control over the reaction dynamics. Using a femtosecond pump-probe scheme, we capture the ultrafast formation time, i.e., 198 ± 16 femtoseconds. The real-time visualization and coherent control of the dynamics contribute to a deeper understanding of the most fundamental bimolecular reactions responsible for C-H bond formation, thus contributing to elucidate the emergence of organic components in the universe.

Theory cracks old data: Rovibrational energy levels of orthoH2-CO derived from experiment

Measurements of rovibrational spectra of clusters provide physical insight only if spectral lines can be assigned to pairs of quantum states, and further insight is obtained if one can deduce the quantitative energy-level pattern. Both steps can be so difficult that some measured spectra remain unassigned, one example is orthoH2-CO. To extend the scope of spectroscopic insights, we propose to use theoretical information in interpretation of spectra. We first performed high accuracy, full-dimensional calculations of the orthoH2-CO spectrum, at the highest practically achievable levels of electronic structure theory and quantum nuclear dynamics. Then, an iterative, theory-guided method developed here allowed us to fully interpret the spectrum of orthoH2-CO, extending the range of van der Waals clusters for which spectroscopy can provide physical insights.

Interaction of methanol with molecular hydrogen: Ab initio potential energy surface and scattering calculations

The potential energy surface (PES) describing the interaction of the methanol molecule with molecular hydrogen has been calculated by the use of the explicitly correlated coupled cluster method, including single, double, and (perturbative) triple excitations [CCSD(T)-F12a] and a correlation-consistent aug-cc-pVTZ basis, with the assumption of fixed molecular geometries. The computed points were fit to a functional form appropriate for time-independent quantum scattering calculations of rotationally inelastic cross sections and rate coefficients. Stationary points on the PES were located, and the global minimum was found to have an energy equal to -254.7 cm-1 relative to the energy of the separated molecules. This PES was used in time-independent close coupling quantum scattering calculations to determine state-to-state cross sections and rate coefficients for rotational transitions within the A- and E-type nuclear spin torsional ground states.

An explicitly correlated six-dimensional potential energy surface for the SiCSi + H2 complex

The first six-dimensional potential energy surface (PES) for the SiCSi + H₂ complex is presented in this work. This surface is developed from a large number of ab initio energies computed at the explicitly correlated coupled-cluster level of theory together with the augmented correlation-consistent polarized valence triple zeta basis set (CCSD(T)-F12/aug-cc-pVTZ). These energies are fitted to an analytical function through a procedure that combines spline, least-squares, and kernel-based methods. Two minimums of similar depths were found at the equilibrium geometry of the SiCSi molecule. The dependence of the PES on the bending angle is analyzed. Furthermore, a reduced four-dimensional PES averaged over the H₂ orientation is presented. Finally, the six-dimensional PES is used for computing the second virial coefficient of the SiCSi + H₂ pair using classical and semi-classical methods.

Stereodynamic control of cold rotationally inelastic CO + HD collisions

Quantum control of molecular collision dynamics is an exciting emerging area of cold collisions. Co-expansion of collision partners in a supersonic molecular beam combined with precise control of their quantum states and alignment/orientation using Stark-induced Adiabatic Raman Passage allows exquisite stereodynamic control of the collision outcome. This approach has recently been demonstrated for rotational quenching of HD in collisions with H₂, D₂, and He and D₂ by He. Here we illustrate this approach for HD(v = 0, j = 2) + CO(v = 0, j = 0) → HD(v' = 0, j') + CO(v' = 0, j') collisions through full-dimensional quantum scattering calculations at collision energies near 1 K. It is shown that the collision dynamics at energies between 0.01-1 K are controlled by an interplay of L = 1 and L = 2 partial wave resonances depending on the final rotational levels of the two molecules. Polarized cross sections resolved into magnetic sub-levels of the initial and final rotational quantum numbers of the two molecules also reveal a significant stereodynamic effect in the cold energy regime. Overall, the stereodynamic effect is controlled by both geometric and dynamical factors, with parity conservation playing an important role in modulating these contributions depending on the particular final state.

Interaction of the HCO radical with molecular hydrogen: Ab initio potential energy surface and scattering calculations

The potential energy surface describing the interaction of the HCO radical with molecular hydrogen has been computed through explicitly correlated coupled cluster calculations including single, double, and (perturbative) triple excitations [RCCSD(T)-F12a], with the assumption of fixed molecular geometries. The computed points were fit to an analytical form suitable for time-independent quantum scattering calculations of rotationally inelastic cross sections and rate coefficients. Since the spin-rotation splittings in HCO are small, cross sections for fine-structure resolved transitions are computed with electron-spin free T matrix elements through the recoupling technique usually employed to determine hyperfine-resolved cross sections. Both spin-free and fine-structure resolved state-to-state cross sections for rotationally inelastic transitions are presented and discussed.

SAPT codes for calculations of intermolecular interaction energies

Symmetry-adapted perturbation theory (SAPT) is a method for calculations of intermolecular (noncovalent) interaction energies. The set of SAPT codes that is described here, the current version named SAPT2020, includes virtually all variants of SAPT developed so far, among them two-body SAPT based on perturbative, coupled cluster, and density functional theory descriptions of monomers, three-body SAPT, and two-body SAPT for some classes of open-shell monomers. The properties of systems governed by noncovalent interactions can be predicted only if potential energy surfaces (force fields) are available. SAPT is the preferred approach for generating such surfaces since it is seamlessly connected to the asymptotic expansion of interaction energy. SAPT2020 includes codes for automatic development of such surfaces, enabling generation of complete dimer surfaces with a rigid monomer approximation for dimers containing about one hundred atoms. These codes can also be used to obtain surfaces including internal degrees of freedom of monomers.

Interaction of the H2S molecule with molecular hydrogen: Ab initio potential energy surface and scattering calculations

The potential energy surface describing the interaction of H₂S with molecular hydrogen has been computed through explicitly correlated coupled cluster calculations including single, double, and (perturbative) triple excitations [CCSD(T)-F12a], with the assumption of fixed molecular geometries. The computed points were fit to an analytical form suitable for time-independent quantum scattering calculations of rotationally inelastic cross sections and rate coefficients. Representative sets of energy dependent state-to-state rotationally inelastic cross sections and rate coefficients are presented and discussed.

First-Principles Many-Body Nonadditive Polarization Energies from Monomer and Dimer Calculations Only: A Case Study on Water

The many-body polarization energy is the major source of nonadditivity in strongly polar systems such as water. This nonadditivity is often considerable and must be included, if only in an average manner, to correctly describe the physical properties of the system. Models for the polarization energy are usually parametrized using experimental data, or theoretical estimates of the many-body effects. Here we show how many-body polarization models can be developed for water complexes using data for the monomer and dimer only using ideas recently developed in the field of intermolecular perturbation theory and state-of-the-art approaches for calculating distributed molecular properties based on the iterated stockholder atoms (ISA) algorithm. We show how these models can be calculated, and we validate their accuracy in describing the many-body nonadditive energies of a range of water clusters. We further investigate their sensitivity to the details of the polarization damping models used. We show how our very best polarization models yield many-body energies that agree with those computed with coupled-cluster methods, but at a fraction of the computational cost.

Complexes of Alkali Metal Cations and Molecular Hydrogen: Potential Energy Surfaces and Bound States

Complexes between metal cations and molecular hydrogen are systems quite amenable for precise spectroscopic and theoretical studies, and at the same time, they are relevant for applications in hydrogen storage and astrochemistry. In this work, we report new intermolecular potential energy surfaces and rovibrational states calculations for complexes involving molecular hydrogen and alkaline metal cations, M⁺-H₂ (M⁺ = Na⁺, K⁺, Rb⁺, Cs⁺). The intermolecular potentials, formulated in an internally consistent way to emphasize differences in the properties of the systems, are represented by simple analytical expressions whose parameters have been optimized from comparison with accurate ab initio calculations. Properties of the low-lying bound states-binding energies, frequencies, and rotational constants-are compared with previous measurements or computations and an overall good agreement is achieved, supporting the reliability of the present formulation. Variations of these properties as a function of the cation size and isotopic substitution, with a proper sequence of ortho and para rotational levels, are also discussed.

Ab initio potential energy surface and microwave spectra for the H2-HCCCN complex

We present a four-dimensional ab initio potential energy surface of the H₂-HCCCN complex at the coupled-cluster singles and doubles with noniterative inclusion of connected triples [CCSD(T)]-F12 level with a large basis set including an additional set of bond functions. The artificial neural networks method was extended to fit the intermolecular potential energy surface. The complex has a planar linear global minimum with the well depth of 199.366 cm^-1 located at R = 5.09 Å, φ = 0°, θ₁ = 0°, and θ₂ = 180°. An additional planar local minimum is also found with a depth of 175.579 cm^-1 that is located at R = 3.37 Å, φ = 0°, θ₁ = 110°, and θ₂ = 104°. The radial discrete variable representation/angular finite basis representation and the Lanczos algorithm were employed to calculate the rovibrational energy levels for four species of H₂-HCCCN (pH₂-HCCCN, oH₂-HCCCN, pD₂-HCCCN, and oD₂-HCCCN). The rotational frequencies and spectroscopic parameters were also determined for four complexes, which agree well with the experimental values.

On the importance of full-dimensionality in low-energy molecular scattering calculations

Scattering of H₂ on CO is of great importance in astrophysics and also is a benchmark system for comparing theory to experiment. We present here a new 6-dimensional potential energy surface for the ground electronic state of H₂-CO with an estimated uncertainty of about 0.6 cm⁻¹ in the global minimum region, several times smaller than achieved earlier. This potential has been used in nearly exact 6-dimensional quantum scattering calculations to compute state-to-state cross-sections measured in low-energy crossed-beam experiments. Excellent agreement between theory and experiment has been achieved in all cases. We also show that the fully 6-dimensional approach is not needed with the current accuracy of experimental data since an equally good agreement with experiment was obtained using only a 4-dimensional treatment, which validates the rigid-rotor approach widely used in scattering calculations. This finding, which disagrees with some literature statements, is important since for larger systems full-dimensional scattering calculations are currently not possible.

Next-Generation Force Fields from Symmetry-Adapted Perturbation Theory

Symmetry-adapted perturbation theory (SAPT) provides a unique set of advantages for parameterizing next-generation force fields from first principles. SAPT provides a direct, basis-set superposition error free estimate of molecular interaction energies, a physically intuitive energy decomposition, and a seamless transition to an asymptotic picture of intermolecular interactions. These properties have been exploited throughout the literature to develop next-generation force fields for a variety of applications, including classical molecular dynamics simulations, crystal structure prediction, and quantum dynamics/spectroscopy. This review provides a brief overview of the formalism and theory of SAPT, along with a practical discussion of the various methodologies utilized to parameterize force fields from SAPT calculations. It also highlights a number of applications of SAPT-based force fields for chemical systems of particular interest. Finally, the review ends with a brief outlook on the future opportunities and challenges that remain for next-generation force fields based on SAPT.

Microwave spectroscopy of the seeded binary and ternary clusters CO-(pH2)2, CO-pH2-He, CO-HD, and CO-(oD2)(N=1,2)

We report the Fourier transform microwave spectra of the a-type J = 1-0 transitions of the binary and ternary CO-(pH2)2, CO-pH2-He, CO-HD, and CO-(oD2)N=1,2 clusters. In addition to the normal isotopologue of CO for all clusters, we observed the transitions of the minor isotopologues, (13)C(16)O, (12)C(18)O, and (13)C(18)O, for CO-(pH2)2 and CO-pH2-He. All transitions lie within 335 MHz of the experimentally or theoretically predicted values. In comparison to previously reported infrared spectra [Moroni et al., J. Chem. Phys. 122, 094314 (2005)], we are able to tentatively determine the vibrational shift for CO-pH2-He, in addition to its b-type J = 1-0 transition frequency. The a-type frequency of CO-pH2-He is similar to that of CO-He2 [Surin et al., Phys. Rev. Lett. 101, 233401 (2008)], suggesting that the pH2 molecule has a strong localizing effect on the He density. Perturbation theory analysis of CO-oD2 reveals that it is approximately T-shaped, with an anisotropy of the intermolecular potential amounting to ∼9 cm(-1).

Spin-orbit quenching of the C+(2P) ion by collisions with para- and ortho-H2

Spin-orbit (de-)excitation of C(+)((2)P) by collisions with H2, a key process for astrochemistry, is investigated. Quantum-mechanical calculations of collisions between C(+) ions and para- and ortho-H2 have been performed in order to determine the cross section for the C(+) (2)P3∕2 → (2)P1∕2 fine-structure transition at low and intermediate energies. The calculation are based on new ab initio potential energy surfaces obtained using the multireference configuration interaction method. Corresponding rate coefficients were obtained for temperatures ranging from 5 to 500 K. These rate coefficients are compared to previous estimations, and their impact is assessed through radiative transfer computation. They are found to increase the flux of the (2)P3∕2 → (2)P1∕2 line at 158 μm by up to 30% for typical diffuse interstellar cloud conditions.

Spectral line shape profile of rovibrational transitions of CO embedded in p-H2 crystals studied by high resolution IR diode laser spectroscopy

Line profiles of rovibrational transitions of CO embedded in p-H(2) crystals were studied by high resolution midinfrared diode laser spectroscopy. The line profile analysis for the R(0)(parallel), R(0)(perpendicular), P(1)(parallel), and P(1)(perpendicular) transitions shows that spectral line shapes are well reproduced by a convolution of Gaussian and Lorentzian functions. The temperature dependence of the Lorentzian Gamma(L)(T) and Gaussian widths Gamma(G)(T) shows that there is a nonzero linewidth contribution to each at the T = 0 K limit. The main part of the Lorentzian width Gamma(L)(T = 0) shows anisotropy in the hcp structure and is explained by spontaneous decay of the rotational excited state energy to phonon modes. A smaller part of Gamma(L)(T = 0) is attributed to inhomogeneous broadening due to the point defects of other CO molecules in the crystal. On the other hand, the Gaussian width Gamma(G)(T = 0) is explained by inhomogeneous broadening due to dislocations. In the T > 0 region, Gamma(L)(T) shows strong temperature dependence but Gamma(G)(T) does not. The center frequencies of the R(0)(perpendicular) and P(1)(parallel) transitions show blueshifts and those of the R(0)(parallel) and P(1)(perpendicular) transitions show redshifts with increasing temperature. This phenomenon is explained by a decrease in the anisotropy in the crystal field, which is caused by the averaging of thermal lattice fluctuations. Furthermore, the contribution of vibration and rotation to the linewidth is discussed.

N(2)O in small para-hydrogen clusters: Structures and energetics

We present the minimum-energy structures and energetics of clusters of the linear N(2)O molecule with small numbers of para-hydrogen molecules with pairwise additive potentials. Interaction energies of (p-H(2))-N(2)O and (p-H(2))-(p-H(2)) complexes were calculated by averaging the corresponding full-dimensional potentials over the H(2) angular coordinates. The averaged (p-H(2))-N(2)O potential has three minima corresponding to the T-shaped and the linear (p-H(2))-ONN and (p-H(2))-NNO structures. Optimization of the minimum-energy structures was performed using a Genetic Algorithm. It was found that p-H(2) molecules fill three solvation rings around the N(2)O axis, each of them containing up to five p-H(2) molecules, followed by accumulation of two p-H(2) molecules at the oxygen and nitrogen ends. The first solvation shell is completed at N = 17. The calculated chemical potential oscillates with cluster size up to the completed first solvation shell. These results are consistent with the available experimental measurements.

One and Two Hydrogen Molecules in the Large Cage of the Structure II Clathrate Hydrate: Quantum Translation−Rotation Dynamics Close to the Cage Wall

We have performed a rigorous theoretical study of the quantum translation-rotation (T-R) dynamics of one and two H2 and D2 molecules confined inside the large hexakaidecahedral (5(12)6(4)) cage of the sII clathrate hydrate. For a single encapsulated H2 and D2 molecule, accurate quantum five-dimensional calculations of the T-R energy levels and wave functions are performed that include explicitly, as fully coupled, all three translational and the two rotational degrees of freedom of the hydrogen molecule, while the cage is taken to be rigid. In addition, the ground-state properties, energetics, and spatial distribution of one and two p-H2 and o-D2 molecules in the large cage are calculated rigorously using the diffusion Monte Carlo method. These calculations reveal that the low-energy T-R dynamics of hydrogen molecules in the large cage are qualitatively different from that inside the small cage, studied by us recently. This is caused by the following: (i) The large cage has a cavity whose diameter is about twice that of the small cage for the hydrogen molecule. (ii) In the small cage, the potential energy surface (PES) for H2 is essentially flat in the central region, while in the large cage the PES has a prominent maximum at the cage center, whose height exceeds the T-R zero-point energy of H2/D2. As a result, the guest molecule is excluded from the central part of the large cage, its wave function localized around the off-center global minimum. Peculiar quantum dynamics of the hydrogen molecule squeezed between the central maximum and the cage wall manifests in the excited T-R states whose energies and wave functions differ greatly from those for the small cage. Moreover, they are sensitive to the variations in the hydrogen-bonding topology, which modulate the corrugation of the cage wall.

Five-dimensional ab initio potential energy surface and predicted infrared spectra of H2–CO2 van der Waals complexes

The authors present a new five-dimensional potential energy surface for H2-CO2 including the Q3 normal mode for the nu3 antisymmetric stretching vibration of the CO2 molecule. The potential energies were calculated using the supermolecular approach with the full counterpoise correction at the CCSD(T) level with an aug-cc-pVTZ basis set supplemented with bond functions. The global minimum is at two equivalent T-shaped coplanar configurations with a well depth of 219.68 cm-1. The rovibrational energy levels for four species of H2-CO2 (paraH2-, orthoH2-, paraD2-, and orthoD2-CO2) were calculated employing the discrete variable representation (DVR) for radial variables and finite basis representation (FBR) for angular variables and the Lanczos algorithm. Our calculations showed that the off-diagonal intra- and intermolecular vibrational coupling could be neglected, and separation of the intramolecular vibration by averaging the total Hamiltonian with the wave function of a specific vibrational state of CO2 should be a good approximation with high accuracy. The calculated band origin shift in the infrared spectra in the nu3 region of CO2 is -0.113 cm-1 for paraH2-CO2 and -0.099 cm-1 for orthoH2-CO2, which agrees well with the observed values of -0.198 and -0.096 cm-1. The calculated rovibrational spectra for H2-CO2 are consistent with the available experimental spectra. For D2-CO2, it is predicted that only a-type transitions occur for paraD2-CO2, while both a-type and b-type transitions are significant for orthoD2-CO2.

Spectroscopy and potential energy surface of the H2–CO2van der Waals complex: experimental and theoretical studies

A 4-D ab initio potential energy surface is calculated for the intermolecular interaction of hydrogen and carbon dioxide, using the CCSD(T) method with a large basis set. The surface has a global minimum with a well depth of 212 cm(-1) and an intermolecular distance of 2.98 A for a planar configuration with both the O-C-O and H-H axes perpendicular to the intermolecular axis. Bound state calculations are performed for the H(2)-CO(2) van der Waals complex with H(2) in both the para and ortho spin states, and the binding energy of paraH(2)-CO(2)(50.4 cm(-1)) is found to be significantly less than that of orthoH(2)-CO(2)(71.7 cm(-1)). The surface supports 7 bound intermolecular vibrational states for paraH(2)-CO(2) and 19 for orthoH(2)-CO(2), and the lower rotational levels with J< or = 4 follow an asymmetric rotor pattern. The calculated infrared spectrum of paraH(2)-CO(2) agrees well with experiment. For orthoH(2)-CO(2), the ground state rotational levels allowed by symmetry are found to have (K(a), K(c))=(even, odd) or (odd, even). This somewhat unexpected fact enables the previously observed experimental spectrum to be assigned for the first time, in good agreement with theory, and indicates that the orientation of hydrogen is perpendicular to the intermolecular axis in the ground state of the orthoH(2)-CO(2) complex.

Hydrogen Molecule in the Small Dodecahedral Cage of a Clathrate Hydrate: Quantum Five-Dimensional Calculations of the Coupled Translation−Rotation Eigenstates

We report quantum five-dimensional (5D) calculations of the energy levels and wave functions of the hydrogen molecule, para-H2 and ortho-H2, confined inside the small dodecahedral (H2O)20 cage of the sII clathrate hydrate. All three translational and the two rotational degrees of freedom of H2 are included explicitly, as fully coupled, while the cage is treated as rigid. The 5D potential energy surface (PES) of the H2-cage system is pairwise additive, based on the high-quality ab initio 5D (rigid monomer) PES for the H2-H2O complex. The bound state calculations involve no dynamical approximations and provide an accurate picture of the quantum 5D translation-rotation dynamics of H2 inside the cage. The energy levels are assigned with translational (Cartesian) and rotational quantum numbers, based on calculated root-mean-square displacements and probability density plots. The translational modes exhibit negative anharmonicity. It is found that j is a good rotational quantum number, while the threefold degeneracy of the j = 1 level is lifted completely. There is considerable translation-rotation coupling, particularly for excited translational states.

A five-dimensional potential energy surface and predicted infrared spectra for the N2O-hydrogen complexes

We present a five-dimensional potential energy surface for the N(2)O-hydrogen complex using supermolecular approach with the full counterpoise correction at the coupled-cluster singles and doubles with noniterative inclusion of connected triple level. The normal mode Q(3) for the nu(3) antisymmetric stretching vibration of the N(2)O molecule was included in the calculations of the potential energies. The radial discrete variable representation/angular finite basis representation method and Lanczos algorithm were employed to calculate the rovibrational energy levels for four species of N(2)O-hydrogen complexes (N(2)O-para-H(2), -ortho-H(2), -ortho-D(2), and -para-D(2)) without separating the inter- and intramolecular vibrations. The calculated band origins are all blueshifted relative to the isolated N(2)O molecule and in good agreement with the experimental values. The calculated rotational spectroscopic constants and molecular structures agree well with the available experimental results. The frequencies and line intensities of the rovibrational transitions in the nu(3) region of N(2)O for the van der Waals ground vibrational state were calculated and compared with the observed spectra. The predicted infrared spectra are consistent with the observed spectra and show that the N(2)O-H(2) complexes are mostly a-type transitions while both a-type and b-type transitions are significant for the N(2)O-D(2) complexes.

HF in clusters of molecular hydrogen: II. Quantum solvation by H2 isotopomers, cluster rigidity, and comparison with CO-doped parahydrogen clusters

We present a comprehensive theoretical study of the quantum solvation of the HF molecule by small clusters of the H2 isotopomers, p-H2, HD, and o-D2, with up to 13 hydrogen solvent molecules. This complements our earlier work on the HF-doped parahydrogen clusters [H. Jiang and Z. Bacic, J. Chem. Phys. 122, 244306 (2005)]. The ground-state properties of the clusters are calculated exactly using the diffusion Monte Carlo method. Detailed information is obtained regarding the size and isotopomer dependences of the energetics, vibrationally averaged structures, and their rigidity. The rigidity of these clusters is investigated further by analyzing the distributions of their principal moments of inertia from the diffusion Monte Carlo simulations. The clusters are found to be rather rigid, especially when compared with the pure parahydrogen clusters of the same size. Extensive comparison is made with the quantum Monte Carlo results for the CO-doped parahydrogen clusters and significant differences are observed in the size evolution of certain properties, notably the chemical potential.