Intermolecular configurations dominated by quadrupole-quadrupole electrostatic interactions: explicit correlation treatment of the five-dimensional potential energy surface and infrared spectra for the CO-N2 complex

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.

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.

Photochemistry of phosphenic chloride (ClPO2): isomerization with chlorine metaphosphite (ClOPO) and reduction by carbon monoxide

Phosphenic chloride (ClPO₂) is an elusive congener of nitryl chloride (ClNO₂). By high-vacuum flash pyrolysis of 2-chloro-1,3,2-dioxaphospholane in the gas phase, ClPO₂ has been efficiently generated and subsequently isolated in cryogenic N₂, Ar, and CO matrices (10 K) for a first time study on its photochemistry. Upon 193 nm laser irradiation, ClPO₂ isomerizes to the novel chlorine metaphosphite (ClOPO) by initial cleavage of the Cl-P bond (→ ˙Cl + ˙PO₂) with subsequent Cl-O bond formation inside the N₂ and Ar matrix cages. The reverse transformation becomes feasible under further irradiation at 266 nm. This photochemistry is consistent with the observed absorptions of ClPO₂ and ClOPO at 207 and 250 nm, respectively. When the photolysis was performed in solid CO ice, no isomerization occurs due to CO-trapping of the initially generated ˙Cl atoms by forming caged radical pair ClCO˙⋯˙PO₂. Concomitantly, photolytic reduction of ClPO₂ to ClPO by CO has been observed, yielding a weakly bonded molecular complex consisting of ClPO and CO₂ bonded through short intermolecular C⋯O contact (2.910 Å). The characterization of ClPO, ClPO₂, ClOPO, and the molecular complexes of ClPO₂-CO and ClPO-CO₂ using matrix-isolation IR and UV-vis spectroscopy is supported by the theoretical calculations at the B3LYP/6-311 + G(3df) level, and the photochemistry of ClPO₂ is also compared with the revisited photochemistry of ClNO₂ in the N₂-matrix.

Efficient prediction for high precision CO-N2 potential energy surface by stacking ensemble DNN

High-dimensional potential energy surface (PES) for van der Waals systems with spectroscopic accuracy, is of great importance for quantum dynamics and an extremely challenge job. CO-N₂ is a typical van der Waals system and its high-precision PES may help elucidate weak interaction mechanisms. Taking CO-N₂ potential energies calculated by CCSD(T)-F12b/aug-cc-pVQZ as the benchmark, we establish an accurate, robust, and efficient machine learning model by using only four molecular structure descriptors based on 7966 benchmark potential energies. The highest accuracy is obtained by a stacking ensemble DNN (SeDNN). Its evaluation parameters MAE, RMSE, and R² reach 0.096, 0.163, 0.9999 cm^-1 , respectively, and the spectroscopic accuracy for vibration spectrum is achieved with predicted PES, which shows SeDNN superior goodness-of-fit and prediction performance. An elaborated PES with the reported global minimum has been predicted with the model, which perfectly reproduces CCSD(T) potential energies and the analytical MLR PES [PCCP, 2018, 20, 2036]. The critical points (global minimum, TSI, TSII, and their barriers), potential curve, and entire PES profile are remarkably consistent with CCSD(T) calculations. To further improve the usability of constructing PESs in practice, the size of the training set (energy points) for the model is reduced to 50%, 30%, and 20% of the database, respectively. The results show that even training with the smallest training set (1593 points), the PES only differs 2.555 cm^-1 with the analytic MLR PES. Therefore, the proposed SeDNN is promisingly an alternative efficient tool to construct subtle PES for van der Waals systems.

Vibrational Energy Transfer in CO+N2 Collisions: A Database for V-V and V-T/R Quantum-Classical Rate Coefficients

Knowledge of energy exchange rate constants in inelastic collisions is critically required for accurate characterization and simulation of several processes in gaseous environments, including planetary atmospheres, plasma, combustion, etc. Determination of these rate constants requires accurate potential energy surfaces (PESs) that describe in detail the full interaction region space and the use of collision dynamics methods capable of including the most relevant quantum effects. In this work, we produce an extensive collection of vibration-to-vibration (V-V) and vibration-to-translation/rotation (V-T/R) energy transfer rate coefficients for collisions between CO and N2 molecules using a mixed quantum-classical method and a recently introduced (A. Lombardi, F. Pirani, M. Bartolomei, C. Coletti, and A. Laganà, Frontiers in chemistry, 7, 309 (2019)) analytical PES, critically revised to improve its performance against ab initio and experimental data of different sources. The present database gives a good agreement with available experimental values of V-V rate coefficients and covers an unprecedented number of transitions and a wide range of temperatures. Furthermore, this is the first database of V-T/R rate coefficients for the title collisions. These processes are shown to often be the most probable ones at high temperatures and/or for highly excited molecules, such conditions being relevant in the modeling of hypersonic flows, plasma, and aerospace applications.

Ab initio investigation of the CO-N2 quantum scattering: The collisional perturbation of the pure rotational R(0) line in CO

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 N₂-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.

Analytic intermolecular potential energy surface and first-principles prediction of the rotational profiles for a symmetric top ion-atom complex: A case study of H3O+-Ar

We presented the first three-dimensional (3D) ab initio intermolecular potential energy surface (PES) for the H₃O⁺-Ar complex. The electronic structure computations were carried out at the explicitly correlated coupled cluster theory-F12 with an augmented correlation-consistent triple zeta basis set. Analytic 3D PES was obtained by least-squares fitting the multi-dimensional Morse/Long-Range (mdMLR) potential model to interaction energies, where the mdMLR function form was applied to the nonlinear ion-atom case for the first time. The 3D PES fitting to 1708 points has root-mean-square deviations of 0.19 cm^-1 with only 108 parameters for interaction energies less than 500 cm^-1. With the 3D PES of the H₃O⁺-Ar complex, we employed the combined radial discrete variable representation/angular finite basis representation method and Lanczos algorithm to calculate rovibrational energy levels. The rotational profiles of the O-H anti-stretching vibrational bands of v₃ ⁺(S)←0⁺ and v₃ ^-(A)←0^- for the H₃O⁺-Ar complex were predicted and were in good agreement with the experimental results.

Full quantum calculation of the rovibrational states and intensities for a symmetric top-linear molecule dimer: Hamiltonian, basis set, and matrix elements

The rovibrational energy levels and intensities of the CH₃F-H₂ dimer have been obtained using our recent global intermolecular potential energy surface [X.-L. Zhang et al., J. Chem. Phys. 148, 124302 (2018)]. The Hamiltonian, basis set, and matrix elements are derived and given for a symmetric top-linear molecule complex. This approach to the generation of energy levels and wavefunctions can readily be utilized for studying the rovibrational spectra of other van der Waals complexes composed of a symmetric top molecule and a linear molecule, and may readily be extended to other complexes of nonlinear molecules and linear molecules. To confirm our method, the rovibrational levels of the H₂O-H₂ dimer have been computed and shown to be in good agreement with experiment and with previous theoretical results. The rovibrational Schrödinger equation has been solved using a Lanczos algorithm together with an uncoupled product basis set. As expected, dimers containing ortho-H₂ are more strongly bound than dimers containing para-H₂. Energies and wavefunctions of the discrete rovibrational levels of CH₃F-paraH₂ complexes obtained from the direct vibrationally averaged 5-dimensional potentials are in good agreement with the results of the reduced 3-dimensional adiabatic-hindered-rotor (AHR) approximation. Accurate calculations of the transition line strengths for the orthoCH₃F-paraH₂ complex are also carried out, and are consistent with results obtained using the AHR approximation. The microwave spectrum associated with the orthoCH₃F-orthoH₂ dimer has been predicted for the first time.

Infrared bands of CS2 dimer and trimer at 4.5 μm

We report observation of new infrared bands of (CS₂)₂ and (CS₂)₃ in the region of the CS₂ ν₁ + ν₃ combination band (at 4.5 µm) using a quantum cascade laser. The complexes are formed in a pulsed supersonic slit-jet expansion of a gas mixture of carbon disulfide in helium. We have previously shown that the most stable isomer of (CS₂)₂ is a cross-shaped structure with D_2d symmetry and that for (CS₂)₃ is a barrel-shaped structure with D₃ symmetry. The dimer has one doubly degenerate infrared-active band in the ν₁ + ν₃ region of the CS₂ monomer. This band is observed to have a rather small vibrational shift of -0.844 cm^-1. We expect one parallel and one perpendicular infrared-active band for the trimer but observe two parallel bands and one perpendicular band. Much larger vibrational shifts of -8.953 cm^-1 for the perpendicular band and -8.845 cm^-1 and +16.681 cm^-1 for the parallel bands are observed. Vibrational shifts and possible vibrational assignments, in the case of the parallel bands of the trimer, are discussed using group theoretical arguments.

Predicted infrared spectra in the HF stretching band of the H2-HF complex

The infrared spectra with hydrogen fluoride (HF) and deuterium fluoride (DF) (v₂ = 1 ← 0) for eight isotropic species of H₂-HF complex are predicted, based on our newly constructed high-accuracy ab initio potential energy surface [D. Yang et al., J. Chem. Phys. 148, 184301 (2018)]. The radial discrete variable representation/angular finite basis representation method and Lanczos algorithm were used to determine the ro-vibrational energy levels and wave functions for eight species of H₂-HF complex (para-H₂-HF, ortho-H₂-HF, para-D₂-HF, ortho-D₂-HF, para-H₂-DF, ortho-H₂-DF, para-D₂-DF, and ortho-D₂-DF) with separating the inter- and intra-molecular vibrations. Bound states properties including their dissociation energies and rotational constants were presented. The calculated band origins are all red shifted to the isolated HF molecule and in good agreement with available experimental values. The frequencies and line intensities of ro-vibrational transitions in the HF stretching band were further calculated, and the predicted infrared spectra are consistent with available observed spectra. Among them, the spectra for three isotopic species of H₂-HF (para-H₂-DF, para-D₂-DF, and ortho-D₂-DF) were predicted for the first time.

Analytic Morse/long-range potential energy surfaces and "adiabatic-hindered-rotor" treatment for a symmetric top-linear molecule dimer: A case study of CH3F-H2

A first effective six-dimensional ab initio potential energy surface (PES) for CH₃F-H₂ which explicitly includes the intramolecular Q₃ stretching normal mode of the CH₃F monomer is presented. The electronic structure computations have been carried out at the explicitly correlated coupled cluster level of theory [CCSD(T)-F12a] with an augmented correlation-consistent triple zeta basis set. Five-dimensional analytical intermolecular PESs for ν₃(CH₃F) = 0 and 1 are then obtained by fitting the vibrationally averaged potentials to the Morse/Long-Range (MLR) potential function form. The MLR function form is applied to the nonlinear molecule-linear molecule case for the first time. These fits to 25 015 points have root-mean-square deviations of 0.74 cm^-1 and 0.082 cm^-1 for interaction energies less than 0.0 cm^-1. Using the adiabatic hindered-rotor approximation, three-dimensional PESs for CH₃F-paraH₂ are generated from the 5D PESs over all possible orientations of the hydrogen monomer. The infrared and microwave spectra for CH₃F-paraH₂ dimer are predicted for the first time. These analytic PESs can be used for modeling the dynamical behavior in CH₃F-(H₂)_N clusters, including the possible appearance of microscopic superfluidity.

Ab initio study of the CO–N2 complex: a new highly accurate intermolecular potential energy surface and rovibrational spectrum

We present a highly accurate ab initio intermolecular potential-energy surface and rovibrational spectrum for the CO–N₂ complex.