Multipole moments, polarizabilities, and hyperpolarizabilities for N2from fourth‐order many‐body perturbation theory calculations

Toward Reliable Dipole Moments without Single Excitations: The Role of Orbital Rotations and Dynamical Correlation

The dipole moment is a crucial molecular property linked to a molecular system's bond polarity and overall electronic structure. To that end, the electronic dipole moment, which results from the electron density of a system, is often used to assess the accuracy and reliability of new electronic structure methods. This work analyses electronic dipole moments computed with the pair coupled cluster doubles (pCCD) ansätze and its linearized coupled cluster (pCCD-LCC) corrections using the canonical Hartree-Fock and pCCD-optimized (localized) orbital bases. The accuracy of pCCD-based dipole moments is assessed against experimental and CCSD(T) reference values using relaxed and unrelaxed density matrices and different basis set sizes. Our test set comprises molecules of various bonding patterns and electronic structures, exposing pCCD-based methods to a wide range of electron correlation effects. Additionally, we investigate the performance of pCCD-in-DFT dipole moments of some model complexes. Finally, our work indicates the importance of orbital relaxation in the pCCD model and shows the limitations of the linearized couple cluster corrections in predicting electronic dipole moments of multiple-bonded systems. Most importantly, pCCD with a linearized CCD correction can reproduce the dipole moment surfaces in singly bonded molecules, which are comparable to the multireference ones.

Recent Advances in Cartesian-Grid DFT in Atoms and Molecules

In the past several decades, density functional theory (DFT) has evolved as a leading player across a dazzling variety of fields, from organic chemistry to condensed matter physics. The simple conceptual framework and computational elegance are the underlying driver for this. This article reviews some of the recent developments that have taken place in our laboratory in the past 5 years. Efforts are made to validate a viable alternative for DFT calculations for small to medium systems through a Cartesian coordinate grid- (CCG-) based pseudopotential Kohn–Sham (KS) DFT framework using LCAO-MO ansatz. In order to legitimize its suitability and efficacy, at first, electric response properties, such as dipole moment (μ), static dipole polarizability (α), and first hyperpolarizability (β), are calculated. Next, we present a purely numerical approach in CCG for proficient computation of exact exchange density contribution in certain types of orbital-dependent density functionals. A Fourier convolution theorem combined with a range-separated Coulomb interaction kernel is invoked. This takes motivation from a semi-numerical algorithm, where the rate-deciding factor is the evaluation of electrostatic potential. Its success further leads to a systematic self-consistent approach from first principles, which is desirable in the development of optimally tuned range-separated hybrid and hyper functionals. Next, we discuss a simple, alternative time-independent DFT procedure, for computation of single-particle excitation energies, by means of “adiabatic connection theorem” and virial theorem. Optical gaps in organic chromophores, dyes, linear/non-linear PAHs, and charge transfer complexes are faithfully reproduced. In short, CCG-DFT is shown to be a successful route for various practical applications in electronic systems.

Molecules under external electric field: On the changes in the electronic structure and validity limits of the theoretical predictions

A detailed analysis of the electronic structure of the ground and first excited spin state of three diatomic molecules ( N2, BH and CO) under static applied electric field is performed at CCSD(T), DFT, MRCI and MRCI(Q) levels of theory. Our findings have revealed that by boosting the applied field one induces changes in the occupation numbers of molecular orbitals, giving rise to changes in the equilibrium geometry and in the HOMO-LUMO energy gap. Specifically, singlet to triplet spin transition can be induced by increasing the applied electric field beyond a critical value. Accordingly, affecting the accuracy of the widely used expression of energy expanded in Taylor series with respect to the applied electric field. © 2018 Wiley Periodicals, Inc.

Finite Field Method for Nonlinear Optical Property Prediction Using Rational Function Approximants

The finite field (FF) method is a quick, easy-to-implement tool for the prediction of nonlinear optical properties. Here, we present and explore a novel variant of the FF method, which uses a rational function to fit a molecule's energy with respect to an electric field. Similarly to previous FF methods, factors crucial for the method's accuracy were tuned. These factors include the number of terms in the function, the distribution of fields used to construct the approximation, and the initial field in the approximation. It was found that the approximant form that best fits the energy has four numerator terms and three denominator terms. To determine a reasonable field distribution, the common ratio of a geometric progression was optimized to √2. Finally, an algorithm for determining a good initial field guess was devised. The optimized FF method was used to compute the polarizability and second hyperpolarizability for a set of 121 molecules and the first hyperpolarizability for a set of 91 molecules. The results from this were compared to a previous polynomial-based FF method. It was found that using a rational function gives higher errors compared to the polynomial model. However, unlike the polynomial model, no subsequent refinement steps were needed to obtain usable results. An overall comparison of the behavior of the two methods also shows that the rational function is less sensitive to the chosen initial field, making it a good choice for new quantum chemistry codes.

A spherical cavity model for quadrupolar dielectrics

The dielectric properties of a fluid composed of molecules possessing both dipole and quadrupole moments are studied based on a model of the Onsager type (molecule in the centre of a spherical cavity). The dielectric permittivity ε and the macroscopic quadrupole polarizability αQ of the fluid are related to the basic molecular characteristics (molecular dipole, polarizability, quadrupole, quadrupolarizability). The effect of αQ is to increase the reaction field, to bring forth reaction field gradient, to decrease the cavity field, and to bring forth cavity field gradient. The effects from the quadrupole terms are significant in the case of small cavity size in a non-polar liquid. The quadrupoles in the medium are shown to have a small but measurable effect on the dielectric permittivity of several liquids (Ar, Kr, Xe, CH4, N2, CO2, CS2, C6H6, H2O, CH3OH). The theory is used to calculate the macroscopic quadrupolarizabilities of these fluids as functions of pressure and temperature. The cavity radii are also determined for these liquids, and it is shown that they are functions of density only. This extension of Onsager's theory will be important for non-polar solutions (fuel, crude oil, liquid CO2), especially at increased pressures.

A finite field method for calculating molecular polarizability tensors for arbitrary multipole rank

A finite field method for calculating spherical tensor molecular polarizability tensors α(lm;l'm') = ∂Δ(lm)/∂ϕ(l'm')* by numerical derivatives of induced molecular multipole Δ(lm) with respect to gradients of electrostatic potential ϕ(l'm')* is described for arbitrary multipole ranks l and l'. Interconversion formulae for transforming multipole moments and polarizability tensors between spherical and traceless Cartesian tensor conventions are derived. As an example, molecular polarizability tensors up to the hexadecapole-hexadecapole level are calculated for water using the following ab initio methods: Hartree-Fock (HF), Becke three-parameter Lee-Yang-Parr exchange-correlation functional (B3LYP), Møller-Plesset perturbation theory up to second order (MP2), and Coupled Cluster theory with single and double excitations (CCSD). In addition, intermolecular electrostatic and polarization energies calculated by molecular multipoles and polarizability tensors are compared with ab initio reference values calculated by the Reduced Variation Space method for several randomly oriented small molecule dimers separated by a large distance. It is discussed how higher order molecular polarizability tensors can be used as a tool for testing and developing new polarization models for future force fields.

Quadrupole moment function and absolute infrared quadrupolar intensities for N2

High level ab initio methods have been used to calculate values of the quadrupole moment of the ground X (1)Sigmag+ state of N2 on a dense radial mesh spanning the interval of 0.8-12.1 a.u. Detailed convergence tests indicate that the resulting equilibrium values of the quadrupole moment Theta(e)=-1.1273 a.u. and its first radial derivative dTheta(R)/dR/e=0.9604 a.u. have absolute uncertainties of 0.3% and 0.8%, respectively, and are more accurate than the best experimental values of these quantities. The calculated quadrupole moment function, together with a recently reported accurate analytic empirical potential energy function [Le Roy et al., J. Chem. Phys. 125, 164310 (2006)], is used to generate values of the radial matrix elements determining the absolute intensities of infrared vibration-rotation transitions of ground-state N2, which take full account of vibration-rotation interactions. These results should improve the reliability of the interpretations of N2 contributions to infrared atmospheric spectra.

The collision-induced polarizability of a pair of hydrogen molecules

Collision-induced light scattering, impulsive stimulated scattering, and subpicosecond-induced birefringence all depend on the transient changes Deltaalpha in molecular polarizabilities that occur when molecules collide. Ab initio results for Deltaalpha are needed to permit comparisons with accurate experimental results for these spectra and for refractive index virial coefficients and dielectric virial coefficients. In this work, we provide results for Deltaalpha for a pair of hydrogen molecules, treated at CCSD(T) level, with an aug-cc-pV5Z (spdf) basis set. Our values replace the best previous ab initio results for the variation of Deltaalpha with intermolecular separation, the self-consistent-field results obtained by Bounds [Mol. Phys. 38, 2099 (1979)] with a relatively small (3s2p) basis set for H2. For the six geometrical configurations studied by Bounds, the inclusion of correlation and improvements in the basis tend to increase both the trace Deltaalpha(0)0 and the anisotropy Deltaalpha2m of the pair polarizability. The change in the anisotropy is relatively small, but our values for the trace differ by factors of 2 or more from Bounds' results. For use in computing experimental line shapes, intensities, and virial coefficients, we have calculated Deltaalpha for 18 different relative orientations of a pair of H2 molecules, with the intermolecular separation R ranging from 2 a.u. (3 a.u. for a linear pair) to 10 a.u. The H2 bond length is fixed at the vibrationally averaged internuclear separation in the ground state r=1.449 a.u. Our results agree well with the CCSD(T) results for Deltaalpha obtained by Maroulis [J. Phys. Chem. A 104, 4772 (2000)] for two pair configurations of H2...H2 (linear and T-shaped) at a fixed internuclear distance of R=6.5 a.u. in a [6s4p1d] basis. As the intermolecular distance increases (for R>or=8 a.u.), the spherical-tensor components of Deltaalpha converge to the results from a long-range model that includes dipole-induced-dipole (DID) interactions, higher-multipole induction, nonuniformity of the local field, hyperpolarization, and van der Waals dispersion. Deviations from the first-order DID model are still evident for R between 8 and 10 a.u. in most orientations of the pair. At shorter range, overlap damping, exchange, and orbital distortion reduce both Deltaalpha0(0) and Deltaalpha(2)0 below their long-range limiting forms.

Multipole moments and polarizabilities of nonpolar diatomic molecules through quantum solvation in solid parahydrogen: The quadrupole moment of N2

It is shown in this paper that from the study of the induced infrared absorption spectra of homonuclear diatomic molecules solvated as impurities in a molecular quantum solid, it is possible to extract information about the rovibrational matrix elements of the multipole moments and polarizability of the embedded molecule. Theoretical expressions are derived for the integrated absorption coefficients of various multipole-field-induced double transitions involving guest-host pairs in a solid para-H(2) matrix. The intensities of some of the quadrupole moment induced transitions involving the N(2)-para-H(2) pair have been measured. From a comparison of the experimental and theoretical intensities, rovibrational matrix elements of the quadrupole moment of N(2) are determined in its ground vibrational state.

Polarizabilities and hyperpolarizabilities for the atoms Al, Si, P, S, Cl, and Ar: Coupled cluster calculations

Accurate static dipole polarizabilities and hyperpolarizabilities are calculated for the ground states of the Al, Si, P, S, Cl, and Ar atoms. The finite-field computations use energies obtained with various ab initio methods including Moller-Plesset perturbation theory and the coupled cluster approach. Excellent agreement with experiment is found for argon. The experimental alpha for Al is likely to be in error. Only limited comparisons are possible for the other atoms because hyperpolarizabilities have not been reported previously for most of these atoms. Our recommended values of the mean dipole polarizability (in the order Al-Ar) are alpha/e(2)a(0) (2)E(h) (-1)=57.74, 37.17, 24.93, 19.37, 14.57, and 11.085 with an error estimate of +/-0.5%. The recommended values of the mean second dipole hyperpolarizability (in the order Al-Ar) are gamma/e(4)a(0) (4)E(h) (-3)=2.02 x 10(5), 4.31 x 10(4), 1.14 x 10(4), 6.51 x 10(3), 2.73 x 10(3), and 1.18 x 10(3) with an error estimate of +/-2%. Our recommended polarizability anisotropy values are Deltaalpha/e(2)a(0) (2)E(h) (-1)=-25.60, 8.41, -3.63, and 1.71 for Al, Si, S, and Cl respectively, with an error estimate of +/-1%. The recommended hyperpolarizability anisotropies are Deltagamma/e(4)a(0) (4)E(h) (-3)=-3.88 x 10(5), 4.16 x 10(4), -7.00 x 10(3), and 1.65 x 10(3) for Al, Si, S, and Cl, respectively, with an error estimate of +/-4%.

Ab initio determination of the electric multipole moments and static (hyper)polarizability of HCCX, X = F, Cl, Br, and I

We report electric multipole moments and (hyper)polarizabilities for the haloethynes HCCX, X = F, Cl, Br, and I. The molecular properties have been obtained from finite-field self-consistent field, Møller-Plesset perturbation theory and coupled cluster calculations with large, carefully optimized basis sets of gaussian-type functions. The mean dipole (hyper)polarizability and the mean quadrupole polarizability near the Hartree-Fock limit are alpha/e(2)a(0) (2)E(h) (-1) = 23.74 (HCCF), 37.26 (HCCCl), 43.97 (HCCBr), 56.44 (HCCI), beta/e(3)a(0) (3)E(h) (-2) = -73.9 (HCCF), -67.0 (HCCCl), -39.5 (HCCBr), 42.7 (HCCI), gamma/e(4)a(0) (4)E(h) (-3) = 4,914 (HCCF), 6,554 (HCCCl), 9,328 (HCCBr), 14,949 (HCCI), and C/e(2)a(0) (4)E(h) (-1) = 160.3 (HCCF), 317.1 (HCCCl), 471.2 (HCCBr), 671.2 (HCCI). Electron correlation has a small effect on the dipole polarizability but affects strongly the hyperpolarizability. Agreement with the available experimental data is more or less fair for HCCF, HCCCl, and HCCBr but less satisfactory for HCCI.