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

Ab Initio Potentials for the Ground S0 and the First Electronically Excited Singlet S1 States of Benzene-Helium with Application to Tunneling Intermolecular Vibrational States

We present new ab initio intermolecular potential energy surfaces for the benzene-helium complex in its ground (S0) and first excited (S1) states. The coupled-cluster level of theory with single, double, and perturbative triple excitations, CCSD(T), was used to calculate the ground state potential. The excited state potential was obtained by adding the excitation energies S0 → S1 of the complex, calculated using the equation of motion approach EOM-CCSD, to the ground state potential interaction energies. Analytical potentials are constructed and applied to study the structural and vibrational dynamics of benzene-helium. The binding energies and equilibrium distances of the ground and excited states were found to be 89 cm-1, 3.14 Å and 77 cm-1, 3.20 Å, respectively. The calculated vibrational energy levels exhibit tunneling of He through the benzene plane and are in reasonable agreement with recently reported experimental values for both the ground and excited states [Hayashi, M.; Ohshima, Y. J. Phys. Chem. Lett. 2020, 11, 9745]. Prospects for the theoretical study of complexes with large aromatic molecules and He are also discussed.

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.

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.

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. 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)].