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.

Theoretical treatment of IO-X (X = N2, CO, CO2, H2O) complexes

Iodine monoxide (IO) is an important component of the biogeochemical cycle of iodine. For instance, it is present in the troposphere, where it plays a crucial role in the physical chemical processes involving iodine containing compounds. Here, we present a theoretical study on a series of atmospherically relevant complexes of IO with N₂, CO, CO₂ and H₂O, where their structural and spectroscopic properties and their interaction energies are computed. Calculations are carried out by means of ab initio post Hartree-Fock (RCCSD(T) and RMP2) methods and density functional theory DFT (PBE0 and M05-2X) based approaches with and without the inclusion of dispersion correction. After comparison to RCCSD(T), we highlight the good performance of M05-2X(+D3) DFT in describing the bonding between IO and X (X = N₂, CO, CO₂, H₂O). Moreover, we found that the IO-X (X = N₂, CO, CO₂, H₂O) complexes are formed by non-covalent interactions between the two monomers. In sum, we characterized two types of complexes: I-bonded and O-bonded, where the former is more stable. The atmospheric implications of the present findings are also discussed such as in the formation of the iodine oxide particles (IOPs).

Theoretical study of the CO2–O2 van der Waals complex: potential energy surface and applications

A four-dimensional-potential energy surface (4D-PES) of the atmospherically relevant carbon dioxide–oxygen molecule (CO2-O2) van der Waals complex is mapped using the ab initio explicitly correlated coupled cluster method with single,...

Explicitly correlated potential energy surface of the CH3Cl-He van der Waals complex and applications

A new 3D-potential energy surface (3D-PES) for the weakly bound CH₃Cl-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 CH₃Cl. 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 CH₃Cl 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-CH₃Cl-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.

IO(X2Π)-Ar cluster: ab initio potential energy surface and dynamical computations

Iodine oxide (IO) is an important tropospheric molecule. In the present paper, we mapped the potential energy surfaces (PESs) of the doubly degenerate IO(X2Π)-Ar van der Waals system using single- and double-excitation coupled cluster approaches with non-iterative perturbation treatment of triple excitations [RCCSD(T)] extrapolated to the complete basis set (CBS) limit. In addition to bent local minima, we identified a linear Ar-IO complex as a global minimum. Afterwards, we performed scattering calculations on these PESs, considering the non-zero spin-orbit contribution and the Renner-Teller effect. The integral cross-sections exhibit an oscillatory structure vs. the final rotational state, as already observed for the NO(X2Π)-Ar system. Moreover, computations reveal that the Ar-IO complex is stable toward dissociation into IO and Ar. Therefore, it can be found in the atmosphere and participates in iodine compound physical chemical processes occurring there.

Characterization and reactivity of the weakly bound complexes of the [H, N, S](-) anionic system with astrophysical and biological implications

We investigate the lowest electronic states of doublet and quartet spin multiplicity states of HNS(-) and HSN(-) together with their parent neutral triatomic molecules. Computations were performed using highly accurate ab initio methods with a large basis set. One-dimensional cuts of the full-dimensional potential energy surfaces (PESs) along the interatomic distances and bending angle are presented for each isomer. Results show that the ground anionic states are stable with respect to the electron detachment process and that the long range parts of the PESs correlating to the SH(-) + N, SN(-) + H, SN + H(-), NH + S(-), and NH(-) + S are bound. In addition, we predict the existence of long-lived weakly bound anionic complexes that can be formed after cold collisions between SN(-) and H or SH(-) and N. The implications for the reactivity of these species are discussed; specifically, it is shown that the reactions involving SH(-), SN(-), and NH(-) lead either to the formation of HNS(-) or HSN(-) in their electronic ground states or to autodetachment processes. Thus, providing an explanation for why the anions, SH(-), SN(-), and NH(-), have limiting detectability in astrophysical media despite the observation of their corresponding neutral species. In a biological context, we suggest that HSN(-) and HNS(-) should be incorporated into H2S-assisted heme-catalyzed reduction mechanism of nitrites in vivo.

Potential energy surface of the CO2-N2 van der Waals complex

Four-dimensional potential energy surface (4D-PES) of the atmospherically relevant CO2-N2 van der Waals complex is generated using the explicitly correlated coupled cluster with single, double, and perturbative triple excitation (CCSD(T)-F12) method in conjunction with the augmented correlation consistent triple zeta (aug-cc-pVTZ) basis set. This 4D-PES is mapped along the intermonomer coordinates. An analytic fit of this 4D-PES is performed. Our extensive computations confirm that the most stable form corresponds to a T-shape structure where the nitrogen molecule points towards the carbon atom of CO2. In addition, we located a second isomer and two transition states in the ground state PES of CO2-N2. All of them lay below the CO2 + N2 dissociation limit. This 4D-PES is flat and strongly anisotropic along the intermonomer coordinates. This results in the possibility of the occurrence of large amplitude motions within the complex, such as the inversion of N2, as suggested in the recent spectroscopic experiments. Finally, we show that the experimentally established deviations from the C2v structure at equilibrium for the most stable isomer are due to the zero-point out-of-plane vibration correction.