Explicit correlation treatment of the six-dimensional potential energy surface and predicted infrared spectra for OCS–H2

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.

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.

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

An ab initio intermolecular potential energy surface (PES) for the van der Waals complex of H₂O-N₂ that explicitly incorporates the intramolecular Q₂ bending normal mode of the H₂O monomer is presented. The electronic structure computations have been carried out at the explicitly correlated coupled cluster theory [CCSD(T)-F12] with an augmented correlation-consistent triple zeta basis set and an additional bond function. Analytic five-dimensional intermolecular PESs for ν₂(H₂O) = 0 and 1 are obtained by fitting to the multi-dimensional Morse/long-range potential function form. These fits to 40 890 points have the root-mean-square (rms) discrepancy of 0.88 cm^-1 for interaction energies less than 2000.0 cm^-1. The resulting vibrationally averaged PESs provide good representations of the experimental microwave and infrared data: for microwave transitions of H₂O-N₂, the rms discrepancy is only 0.0003 cm^-1, and for infrared transitions of the A₁ symmetry of the H₂O(ν₂ = 1 ← 0)-N₂, the rms discrepancy is 0.001 cm^-1. The calculated infrared band origin shifts associated with the ν₂ bending vibration of water are 2.210 cm^-1 and 1.323 cm^-1 for H₂O-N₂ and D₂O-N₂, respectively, in good agreement with the experimental values of 2.254 cm^-1 and 1.266 cm^-1. The benchmark tests and comparisons of the predicted spectral properties are carried out between CCSD(T)-F12a and CCSD(T)-F12b approaches.

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.

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.

A full-dimensional potential energy surface and quantum dynamics of inelastic collision process for H2-HF

A full-dimensional ab initio potential energy surface for the H₂-HF van der Waals complex was constructed by employing the coupled-cluster singles and doubles with noniterative inclusion of connected triples with augmented correlation-consistent polarised valence quadruple-zeta basis set plus bond functions. Using the improved coupled-states approximation including the nearest neighbor Coriolis couplings, we calculated the state-to-state scattering dynamics for pure rotational and ro-vibrational energy transfer processes. For pure rotational energy transfer, our results showed a different dynamical behavior for para-H₂ and ortho-H₂ in collision with hydrogen fluoride (HF), which is consistent with the previous study. Interestingly, some strong resonant peaks were presented in the cross sections for ro-vibrational energy transfer. In addition, the calculated vibrational-resolved rate constant is in agreement with the experimental results reported by Bott et al. These dynamics data can be further applied to the numerical simulation of HF chemical lasers.

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.

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 thorough understanding of the intermolecular configurations of van der Waals complexes is a great challenge due to their weak interactions, floppiness and anharmonic nature. Although high-resolution microwave or infrared spectroscopy provides one of the most direct and precise pieces of experimental evidence, the origin and key role in determining such intermolecular configurations of a van der Waals system strongly depend on its highly accurate potential energy surface (PES) and a detailed analysis of its ro-vibrational wavefunctions. Here, a new five-dimensional potential energy surface for the van der Waals complex of CO-N₂ which explicitly incorporates the dependence on the stretch coordinate of the CO monomer is generated using the explicitly correlated couple cluster (CCSD(T)-F12) method in conjunction with a large basis set. Analytic four-dimensional PESs are obtained by the least-squares fitting of vibrationally averaged interaction energies for v = 0 and v = 1 to the Morse/Long-Range potential mode (V_MLR). These fits to 7966 points have root-mean-square deviations (RMSD) of 0.131 cm^-1 and 0.129 cm^-1 for v = 0 and v = 1, respectively, with only 315 parameters. Energy decomposition analysis is carried out, and it reveals that the dominant factor in controlling intermolecular configurations is quadrupole-quadrupole electrostatic interactions. Moreover, the rovibrational levels and wave functions are obtained for the first time. The predicted infrared transitions and intensities for the ortho-N₂-CO complex as well as the calculated energy levels for para-N₂-CO are in good agreement with the available experimental data with RMSD discrepancies smaller than 0.068 cm^-1. The calculated infrared band origin shift associated with the fundamental band frequency of CO is -0.721 cm^-1 for ortho-N₂-CO which is in excellent agreement with the experimental value of -0.739 cm^-1. The agreement with experimental values validates the high quality of the PESs and enhances our confidence to explain the observed mystery lines around 2163 cm^-1.