Explicitly correlated ab initio potential energy surface and predicted rovibrational spectra for H2O-N2 and D2O-N2 complexes

A New Ab Initio Potential Energy Surface and Rovibrational Spectra for the CO-N2O Complex

We constructed a new ab initio potential energy surface (PES) for CO-N2O which includes the intramolecular Q3 normal coordinate for the N2O ν3 antisymmetric stretching vibration. The intermolecular potential was evaluated employing the supermolecular method at the [CCSD(T)]-F12a level, with the aug-cc-pVTZ basis set plus bond functions. By integral over the intramolecular Q3 coordinate, we obtained the vibrationally averaged PESs for the CO-N2O system in the ground and ν3 excited states of N2O. Each PES features one nearly T-shaped global minimum and one skewed T-shaped local minimum. Based on these obtained PESs of CO-N2O, the radial discrete variable representation/angle finite base representation method and the Lanczos algorithm were applied for the calculations of bound states and rovibrational energy levels. The calculated ν3 vibrational band origin shift of the N2O monomer in CO-N2O is 2.7570 cm-1, matching well with the observed value of 2.9048 cm-1. The computed microwave and infrared transition frequencies, as well as the rotational parameters, are consistent with the experimental observations.

Towards the generation of potential energy surfaces of weakly bound medium-sized molecular systems: the case of benzonitrile-He complex

Currently, the explicitly correlated coupled cluster method is used routinely to generate the multi-dimensional potential energy surfaces (mD-PESs) of van der Waals complexes of small molecular systems relevant for atmospheric, astrophysical and industrial applications. Although very accurate, this method is computationally prohibitive for medium and large molecules containing clusters. For instance, the recent detections of complex organic molecules (COMs) in the interstellar medium, such as benzonitrile, revealed the need to establish an accurate enough electronic structure approach to map the mD-PESs of these species interacting with the surrounding gases. As a benchmark, we have treated the case of the polar molecule benzonitrile interacting with helium, where we use post-Hartree-Fock and symmetry-adapted perturbation theory (SAPT) techniques. Accordingly, we show that MP2 and distinguishable-cluster approximation (DCSD) cannot be used for this purpose, whereas accurate enough PESs may be obtained using the corresponding explicitly correlated versions (MP2-F12 or DCSD-F12) with a reduction in computational costs. Alternatively, computations revealed that SAPT(DFT) is as performant as CCSD(T)-F12/aug-cc-pVTZ, making it the method of choice for mapping the mD-PESs of COMs containing clusters. Therefore, we have used this approach to generate the 3D-PES of the benzonitrile-He complex along the intermonomer Jacobi coordinates. As an application, we have incorporated the analytic form of this PES into quantum dynamical computations to determine the cross sections of the rotational (de-)excitation of benzonitrile colliding with helium at a collision energy of 10 cm-1.

Understanding the high-resolution spectral signature of the N2-H2O van der Waals complex in the 2OH stretch region

We present the observation of the N2-H2O van der Waals complex in the 2OH stretch overtone region. The high-resolution jet cooled spectra were measured using a sensitive continuous wave cavity ringdown spectrometer. Several bands were observed and vibrationally assigned in terms of ν1, ν2, and ν3, the vibrational quantum numbers of the isolated H2O molecule, as (ν1'ν2'ν3')←(ν1″ν2″ν3″)=(200)←(000) and (101) ← (000). A combination band involving the excitation of the in-plane bending motion of N2 and the (101) vibration of water is also reported. The spectra were analyzed using a set of four asymmetric top rotors, each associated with a nuclear spin isomer. Several local perturbations of the (101) vibrational state were observed. These perturbations were assigned to the presence of the nearby (200) vibrational state and to the combination of (200) with intermolecular modes.

MLRNet: Combining the Physics-Motivated Potential Models with Neural Networks for Intermolecular Potential Energy Surface Construction

A physics-based machine learning model called MLRNet has been developed to construct the high-accuracy two-body intermolecular potential energy surface (IPES). The outputs of the neural network are integrated into the physically realistic Morse/long-range (MLR) function, which ensures that the MLRNet has meaningful extrapolation at both short and long ranges and solves the asymptotic problem in common neural network potential (NNP) models. The neural network representation of the MLR parameters is more flexible and more efficient than the polynomial expansion in the conventional mdMLR model, especially for systems containing nonrigid monomer(s). The present work illustrates the basic framework of the current MLRNet model, including (i) how to combine the physically meaningful MLR function with different possible NN structures, (ii) the preservation of permutation symmetry, and (iii) the predetermination of the long-range function u_LR. We choose two realistic systems to demonstrate the performance of MLRNet: the three-dimensional IPES of CO₂-He including the CO₂ antisymmetric vibration Q₃ and the six-dimensional IPES of the H₂O-Ar system. In both cases, the fitting errors of the MLRNet are several times smaller than those of the conventional mdMLR model. Both short-range and long-range extrapolation tests were performed to illustrate the extrapolation ability of the MLRNet and its damping function version. Moreover, for the 6-D H₂O-Ar system, the MLRNet only needs 1596 trainable parameters, which is almost equal to the number needed for the 5-D mdMLR model (1509) and half that needed for the PIP-NN model (3501) within similar accuracy, which illustrates the model efficiency in high-dimensional IPES fitting.

Basis sets dependency in constructing spectroscopy-accuracy Ab Initio global electric dipole moment functions

Recently, more attention have been paid on the construction of dipole moment functions (DMF) using theoretical methods. However, the computational methods to construct DMFs are not validated as much as those for potential energy surfaces do. In this letter, using Ar ⋯ He as an example, we tested how spectroscopy-accuracy DMFs can be constructed using ab initio methods. We especially focused on the basis set dependency in this scenario, i.e., the convergence of DMF with the sizes of basis sets, basis set superposition error, and mid-bond functions. We also tested the explicitly correlated method, which converges with smaller basis sets than the conventional methods do. This work can serve as a pictorial sample of all these computational technologies behaving in the context of constructing DMFs.

Spectroscopic study of the tunneling dynamics in N2-water observed in the O-D stretch region

The O-D stretch rovibrational spectra of N₂-D₂O and N₂-DOH were measured and analyzed. A combination band involving the in-plane N₂ bending vibration was also observed. These bands were recorded using a pulsed-slit supersonic jet expansion and a mid-infrared tunable optical parametric oscillator. The spectra were analyzed by considering the feasible tunneling motions, and transitions were fitted to independent asymmetric rotors for each tunneling state. The rotational constants of the four tunneling components of N₂-D₂O were retrieved for the excited vibrational states. A two order of magnitude increase in the tunneling splittings is observed for the asymmetric O-D stretch (ν₃ in D₂O) excitation compared to the symmetric stretch (ν₁ in D₂O) and to the ground vibrational state. This last finding indicates that the ν₃ vibrational state is likely perturbed by a combination state that includes ν₁. Finally, the observation of a local perturbation in the ν₃ vibrational band, affecting the positions of few rovibrational levels, provides an experimental lower limit of the dissociation energy of the complex, D₀ > 120 cm^-1.

Quantum nature of molecular vibrational quenching: Water-molecular hydrogen collisions

Rates of conversions of molecular internal energy to and from kinetic energy by means of molecular collision allow us to compute collisional line shapes and transport properties of gases. Knowledge of ro-vibrational quenching rates is necessary to connect spectral observations to physical properties of warm astrophysical gasses, including exo-atmospheres. For a system of paramount importance in this context, the vibrational bending mode quenching of H₂O by H₂, we show here that the exchange of vibrational to rotational and kinetic energy remains a quantum process, despite the large numbers of quantum levels involved and the large vibrational energy transfer. The excitation of the quantized rotor of the projectile is by far the most effective ro-vibrational quenching path of water. To do so, we use a fully quantum first-principles computation, potential and dynamics, converging it at all stages, in a full coupled channel formalism. We present here rates for the quenching of the first bending mode of ortho-H₂O by ortho-H₂, up to 500 K, in a fully converged coupled channel formalism.