Dipole polarizability computation

Polarization is usually accounted for by computing the interaction between induced dipoles. The induced dipole is computed by multiplying the atomic polarizability by the electric field present at that nucleus. The electric field used is often only that due to the charges of the other region of the system. In a few calculations, the MM charges have been included in the orbital-based calculation itself as an interaction with point charges. 

TOO Schwerdtfeger, P. (2006) Atomic Static Dipole Polarizabilities, in Computational Aspects of Electric Polarizability Calculations Atoms, Molecules and Clusters (ed. G. Maroulis), Imperial College Press, London, pp. 1-32. 

When Jens Oddershede was elected a Fellow of the American Physical Society in 1993, the citation read For contribution to the theory, computation, and understanding of molecular response properties, especially through the elucidation implementation of the Polarization Propagator formalism. Although written more than a decade ago, it is still true today. The common thread that has run through Jens work for the past score of years is development of theoretical methods for studying the response properties of molecules. His primary interest has been in the development and applications of polarization propagator methods for direct calculation of electronic spectra, radiative lifetime and linear and non-linear response properties such as dynamical dipole polarizabilities and... 

Ahlrichs43 has also described in more detail the method used in this work. Werner and Meyer44 have continued work with their version of the PNO-CI and CEPA methods,23 computing the static dipole polarizabilities of CH4. Various basis sets were used in this study, mostly with better results than in earlier work. Meyer also studied the energy surface of CH4+ in earlier work with this method.23... 

Keywords computer modelling, density functional theory, dipole moment, dipole polarizability,... 

Table XXXIII compares the estimated a with those computed. The accuracy of this simple estimator suggests that the isotropic polarizability is significantly determined by the heavy atoms, not the hydrogens, and that multiple bonds are, as expected, less polarizable. The anisotropy of the dipole polarizability, of course, has very much to do with the hydrogens.

The dipole polarizability can be used in place of the dipole moment function, and this will lead to Raman intensities. Likewise, one can compute electrical quadrupole and higher multipole transition moments if these are of interest. 

Atom-Atom Interactions. - The methods applied, usually to interactions in the inert gases, are a natural extension of diatomic molecule calculations. From the interaction potentials observable quantities, especially the virial coefficients can be calculated. Maroulis et al.31 have applied the ab initio finite field method to calculate the interaction polarizability of two xenon atoms. A sequence of new basis sets for Xe, especially designed for interaction studies have been employed. It has been verified that values obtained from a standard DFT method are qualitatively correct in describing the interaction polarizability curves. Haskopoulos et al.32 have applied similar methods to calculate the interaction polarizability of the Kr-Xe pair. The second virial coefficients of neon gas have been computed by Hattig et al.,33 using an accurate CCSD(T) potential for the Ne-Ne van der Waals potential and interaction-induced electric dipole polarizabilities and hyperpolarizabilities also obtained by CCSD calculations. The refractivity, electric-field induced SHG coefficients and the virial coefficients were evaluated. The authors claim that the results are expected to be more reliable than current experimental data. 

Recently, a new theoretical method of calculating potential energy and dipole/polarizability surfaces for van der Waals molecules based on symmetry-adapted perturbation theory (sapt) of intermolecular forces (12)— (15) has been developed (16)-(24). In this method, referred to as many-body symmetry-adapted perturbation theory, all physically important contributions to the potential and the interaction-induced properties, such as electrostatics, exchange, induction, and dispersion are identified and computed separately. By making a perturbation expansion in the intermolecular interaction as well as in the intramolecular electronic correlation, it is possible to sum the correlation contributions to the different physical... 

Burrows and Cohen used their solutions found with Maple to compute the dipole polarizabilities and shielding factors of a general ns state, obtaining results that improve upon recent accurate calculations reported by Montgomery [44] and by Laughlin [43]. 

Note that a distinction is made between electrostatic and polarization energies. Thus the electrostatic term, Ue e, here refers to an interaction between monomer charge distributions as if they were infinitely separated (i.e., t/°le). A perturbative method is used to obtain polarization as a separate entity. The electrostatic and polarization contributions are expressed in terms of multipole expansions of the classical coulomb and induction energies. Electrostatic interactions are computed using a distributed multipole expansion up to and including octupoles at atom centers and bond midpoints. The polarization term is calculated from analytic dipole polarizability tensors for each localized molecular orbital (LMO) in the valence shell centered at the LMO charge centroid. These terms are derived from quantum calculations on the... 

Abstract The modified equation-of-motion coupled cluster approach of Sekino and Bartlett is extended to computations of the mixed electric-dipole/magnetic-dipole polarizability tensor associated with optical rotation in chiral systems. The approach - referred to here as a linearized equation-of-motion coupled cluster (EOM-CCl) method - is a compromise between the standard EOM method and its linear response counterpart, which avoids the evaluation of computationally expensive terms that are quadratic in the field-perturbed wave functions, but still yields properties that are size-extensive/intensive. Benchmark computations on five representative chiral molecules, including (P)-hydrogen peroxide, (5)-methyloxirane, (5 )-2-chloropropioniuile, (/ )-epichlorohydrin, and (75,45)-norbornenone, demonstrate typically small deviations between the EOM-CCl results and those from coupled cluster linear response theory, and no variation in the signs of the predicted rotations. In addition, the EOM-CCl approach is found to reduce the overall computing time for multi-wavelength-specific rotation computations by up to 34%. 

The earliest calculations of the dipole and quadrupole contributions to the crystal field used free-ion values for the dipole polarizability, a,-, and for the quadrupole polarizability (Hutchings and Ray, 1963). It is generally recognized today, however, that the polarizability of an ion in a solid is generally much less than the corresponding value for the free ion (Chakrabarti et al., 1976 Bogomolova et al., 1977). How much less, however, is not clear. For this reason, therefore, computations of dipole and quadrupole contributions to the crystal field should be regarded as imprecise, even if they are performed with reduced polarizability values. 

---