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

Insight into the Binding of Argon to Cyclic Water Clusters from Symmetry-Adapted Perturbation Theory

This work systematically examines the interactions between a single argon atom and the edges and faces of cyclic H2O clusters containing three-five water molecules (Ar(H2O)n=3-5). Full geometry optimizations and subsequent harmonic vibrational frequency computations were performed using MP2 with a triple-ζ correlation consistent basis set augmented with diffuse functions on the heavy atoms (cc-pVTZ for H and aug-cc-pVTZ for O and Ar; denoted as haTZ). Optimized structures and harmonic vibrational frequencies were also obtained with the two-body-many-body (2b:Mb) and three-body-many-body (3b:Mb) techniques; here, high-level CCSD(T) computations capture up through the two-body or three-body contributions from the many-body expansion, respectively, while less demanding MP2 computations recover all higher-order contributions. Five unique stationary points have been identified in which Ar binds to the cyclic water trimer, along with four for (H2O)4 and three for (H2O)5. To the best of our knowledge, eleven of these twelve structures have been characterized here for the first time. Ar consistently binds more strongly to the faces than the edges of the cyclic (H2O)n clusters, by as much as a factor of two. The 3b:Mb electronic energies computed with the haTZ basis set indicate that Ar binds to the faces of the water clusters by at least 3 kJ mol-1 and by nearly 6 kJ mol-1 for one Ar(H2O)5 complex. An analysis of the interaction energies for the different binding motifs based on symmetry-adapted perturbation theory (SAPT) indicates that dispersion interactions are primarily responsible for the observed trends. The binding of a single Ar atom to a face of these cyclic water clusters can induce perturbations to the harmonic vibrational frequencies on the order of 5 cm-1 for some hydrogen-bonded OH stretching frequencies.

Ab initio relativistic potential energy surfaces of benzene-Xe complex with application to intermolecular vibrations

The benzene-Xe (BXe) complex in its electronic ground state is studied using ab initio methods. Since this complex contains the heavy Xe atom, the relativistic effects cannot be neglected. We test two different approaches that describe the scalar relativistic effects in the framework of the coupled-cluster level of theory with single, double, and perturbative triple excitations, used for the interaction energy calculations. The first one is based on the small core pseudopotential (PP), and the second one is based on the explicit treatment of scalar relativistic effects using the Douglas-Kroll-Hess (DKH) Hamiltonian. A few basis sets are tested with the PP and DKH, and for each one, the analytical potential energy surface (PES) is constructed. It is shown that the difference between PESs determined with PP and DKH methods is small, if the orbitals of the 4d subshell in Xe are correlated. We select the most appropriate approach for the calculation of the potential energy surface of BXe, with respect to accuracy and computational cost. The optimal level of theory includes a small Dunning's basis set for the benzene monomer and a larger PP basis set for Xe supplemented by midbond functions. The PES obtained using such an approach provides a reasonable accuracy when compared to the empirical one derived from the microwave spectra of BXe. The empirical and the theoretical values of intermolecular vibrational energies agree within 0.5 cm^-1 up to second overtones. The vibrational energy level pattern of BXe is characterized by a distinct polyad structure.

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.