Character of intermolecular interaction in pyridine-argon complex: Ab initio potential energy surface, internal dynamics, and interrelations between SAPT energy components

Long-range interactions of aromatic molecules with alkali-metal and alkaline-earth-metal atoms

The isotropic and anisotropic coefficients C_n ^l,m of the long-range spherical expansion ∼1/Rⁿ (R-the intermolecular distance) of the dispersion and induction intermolecular energies are calculated using the first principles for the complexes containing an aromatic molecule (benzene, pyridine, furan, and pyrrole) and alkali-metal (Li, Na, K, Rb, and Cs) or alkaline-earth-metal (Be, Mg, Ca, Sr, and Ba) atoms in their electronic ground states. The values of the first- and second-order properties of the aromatic molecules are calculated using the response theory with the asymptotically corrected LPBE0 functional. The second-order properties of the closed-shell alkaline-earth-metal atoms are obtained using the expectation-value coupled cluster theory and of the open-shell alkali-metal atoms using analytical wavefunctions. These properties are used for the calculation of the dispersion C_n,disp ^l,m and induction C_n,ind ^l,m coefficients (C_n ^l,m=C_n,disp ^l,m+C_n,ind ^l,m) with n up to 12 using the available implemented analytical formulas. It is shown that the inclusion of the coefficients with n > 6 is important for reproducing the interaction energy in the van der Waals region at R ≈ 6 Å. The reported long-range potentials should be useful for constructing the analytical potentials valid for the whole intermolecular interaction range, which are needed for spectroscopic and scattering studies.

Structures of the (Imidazole)nH+ ... Ar (n=1,2,3) complexes determined from IR spectroscopy and quantum chemical calculations

Here, we present new cryogenic infrared spectra of the (Imidazole)$$_{n}\mathrm{H}^{+}$$ n H + (n=1,2,3) ions. The data was obtained using helium tagging infrared predissociation spectroscopy. The new results were compared with the data obtained by Gerardi et al. (Chem. Phys. Lett. 501:172–178, 2011) using the same technique but with argon as a tag. Comparison of the two experiments, assisted by theoretical calculations, allowed us to evaluate the preferable attachment positions of argon to the (Imidazole)$$_{n}\mathrm{H}^{+}$$ n H + frame. Argon attaches to nitrogen-bonded hydrogen in the case of the (Imidazole)H$$^+$$ + ion, while in (Imidazole)$$_{2}\mathrm{H}^{+}$$ 2 H + and (Imidazole)$$_{3}\mathrm{H}^{+}$$ 3 H + the preferred docking sites for the argon are in the center of the complex. This conclusion is supported by analyzing the spectral features attributed to the N–H stretching vibrations. Symmetry adapted perturbation theory (SAPT) analysis of the non-covalent forces between argon and the (Imidazole)$$_{n}\mathrm{H}^{+}$$ n H + (n=1,2,3) frame revealed that this switch of docking preference with increasing complex size is caused by an interplay between induction and dispersion interactions.

Theoretical study of the complexes of dichlorobenzene isomers with argon. II. SAPT analysis of the intermolecular interaction

The interaction of argon with dichlorobenzene isomers (DCB-Ar) has been analyzed with the help of the symmetry-adapted perturbation theory based on the density functional description of monomer properties (DFT-SAPT). The global potential energy surface (PES) of these complexes determined from the DFT-SAPT interaction energy (E_int) values has been compared to the CCSD(T) (coupled cluster method including single and double excitations with perturbative triple excitations) PES reported in the companion Paper I [J. Makarewicz and L. Shirkov, J. Chem. Phys. 150, 074301 (2019)]. The equilibrium structures and the binding energies found using DFT-SAPT and CCSD(T) methods combined with adequate basis sets are in good agreement. Besides DCB-Ar, we confirmed that DFT-SAPT gives accurate values of these quantities for other complexes containing an aromatic molecule and Ar. However, DFT-SAPT PES of DCB-Ar is flatter than the corresponding CCSD(T) one. As a result, the intermolecular vibrational energies are systematically underestimated. The analytical form of the important interrelations between SAPT components of E_int, established previously by us [J. Makarewicz and L. Shirkov, J. Chem. Phys. 144, 204115 (2016)], has been approved for the DCB-Ar complexes. Simplified SAPT models based on these relations have been employed to explain physical reasons for differences in the structures and the binding energies of DCB-Ar isomers. It is shown that the equilibrium distance of Ar to DCB plane and the binding energy are determined mainly by dispersion energy. The shift of Ar toward Cl is caused by both exchange and dispersion terms.

Theoretical study of the complexes of dichlorobenzene isomers with argon. I. Global potential energy surface for all the isomers with application to intermolecular vibrations

The complexes of para- (p-), meta- (m-), and ortho- (o-)dichlorobenzene (DCB) isomers with argon are studied using an ab initio method. The interaction energy in the ground electronic state of the complexes has been calculated using the CCSD(T) method (coupled cluster method including single and double excitations with perturbative triple excitations) and Dunning's double-ζ (aug-cc-pVDZ) basis set supplemented by midbond functions. Local interaction parameters have been defined and interesting relations fulfilled by them, independent of the DCB isomer, have been revealed. This finding has allowed us to construct the accurate global analytical intermolecular potential energy surface for all the DCB-Ar complexes with the same set of parameters, except for the monomer geometries. Each complex is characterized by two symmetrically equivalent global minima, one located above and the other located below the monomer plane at distances equal to 3.497 Å, 3.494 Å, and 3.485 Å for p-, m-, and o-isomers of DCB bound to Ar, respectively. Additionally, the Ar atom is shifted from the geometrical center of the DCB monomer towards the chlorine atoms by the value x_e of 0.182 Å for m-isomer and 0.458 Å for o-isomer. The calculated binding energy D_e of 460 cm^-1, 465 cm^-1, and 478 cm^-1 for p-, m-, and o-complex, respectively, are related to x_e by simple relations. The intermolecular bending fundamentals calculated from PES depend strongly on the isomer structure. The calculated dissociation energies fit in the intervals estimated by the experiment of Gaber et al. for the S₀ state [Phys. Chem. Chem. Phys. 11, 1628 (2009)].