Time-dependent coupled-cluster calculations of polarizabilities and dispersion energy coefficients

On the geometric dependence of the molecular dipole polarizability in water: A benchmark study of higher-order electron correlation, basis set incompleteness error, core electron effects, and zero-point vibrational contributions

In this work, we investigate how geometric changes influence the static dipole polarizability (α) of a water molecule by explicitly computing the corresponding dipole polarizability surface (DPS) across 3125 total (1625 symmetry-unique) geometries using linear response coupled cluster theory including single, double, and triple excitations (LR-CCSDT) and the doubly augmented triple-ζ basis set (d-aug-cc-pVTZ). Analytical formulae based on power series expansions of this ab initio surface are generated using linear least-squares analysis and provide highly accurate estimates of this quantity as a function of molecular geometry (i.e., bond and angle variations) in a computationally tractable manner. An additional database, which consists of 25 representative molecular geometries and incorporates a more thorough treatment of both basis sets and core electron effects, is provided as a current benchmark for this quantity and the corresponding leading-order C ₆ dispersion coefficient. This database has been utilized to assess the importance of these effects as well as the relative accuracy that can be obtained using several quantum chemical methods and a library of density functional approximations. In addition to high-level electron correlation methods (like CCSD) and our analytical least-squares formulae, we find that the SCAN0, PBE0, MN15, and B97-2 hybrid functionals yield the most accurate descriptions of the molecular polarizability tensor in H₂O. Using first-order perturbation theory, we compute the zero-point vibrational correction to α at the CCSDT/d-aug-cc-pVTZ level and find that this correction contributes approximately 3% to the isotropic (α _iso) and nearly 50% to the anisotropic (α _aniso) polarizability values. In doing so, we find that α _iso = 9.8307 bohr³, which is in excellent agreement with the experimental value of 9.83 ± 0.02 bohr³ provided by Russell and Spackman. The DPS reported herein provides a benchmark-quality quantum mechanical estimate of this fundamental quantity of interest and should find extensive use in the development (and assessment) of next-generation force fields and machine-learning based approaches for modeling water in complex condensed-phase environments.

Vibrationally averaged isotropic dispersion energy coefficients of the parahydrogen dimer

We compare the sum-over-states and coupled cluster linear response formalisms for the determination of imaginary-frequency polarizabilities of H(2). Using both approaches, we compute isotropic dispersion energy coefficients C(n) (n = 6, 8, 10) for H(2)-H(2) molecular pairs over a wide range of H(2) bond lengths. We present vibrationally averaged dispersion energy coefficients for H(2)-H(2), H(2)-D(2), and D(2)-D(2) molecular pairs and examine the coefficients' convergence with respect to basis set.

Intermolecular potential energy surface and second virial coefficients for the nonrigid water-CO dimer

A seven-dimensional potential energy surface is calculated for the interaction of water and carbon monoxide using second-order Moller-Plesset theory, coupled-cluster theory, and extrapolated intermolecular perturbation theory. The effects of stretching the CO molecule and bending the water molecule are included. The minimum energy structure of the water-CO dimer changes from an H-C hydrogen bond to an H-O hydrogen bond when the CO bond length increases by less than 10 pm from its equilibrium value. Second virial coefficients for the water-CO interaction are calculated for a wide range of temperatures and compared with the limited experimental data. Allowing the CO bond length and water bond angle to vary has little effect on the second virial coefficients.

A Comparison of Møller-Plesset and Coupled Cluster Linear Response Theory Methods for the Calculation of Dipole Oscillator Strength Sum Rules and C6 Dispersion Coefficients

Four correlated linear response theory methods - the second order polarization propagator approximation (SOPPA), the second order polarization propagator approximation with coupled cluster singles and doubles amplitudes, SOPPA(CCSD), the CC2 and coupled cluster singles doubles (CCSD) linear response theory - were used to determine the dipole oscillator strength sum rules of the hydrogen halides HX (with X = F, Cl, Br and I) and the C6 dispersion coefficient for all pairs of interacting HX molecules via numerical integration of the Casimir-Polder formula. The dependence of the polarizabilities, their frequency dependence and the C6 coefficients on the level of correlation and the dependence of the C6 coefficients on the two intramolecular bond lengths were studied.