Partitioning of higher multipole polarizabilities: numerical evaluation of transferability

In this work, the partitioning of higher multipole polarizabilities, such as dipole-quadrupole, quadrupole-dipole, and quadrupole-quadrupole polarizabilities, into atomic contributions is studied. Partitioning of higher multipole polarizabilities is necessary in the study of accurate interaction energies where dispersion interactions are of importance. The fractional occupation Hirhsfeld-I (FOHI) method is used to calculate the atomic polarizabilities and is briefly explained together with the methodology for partitioning of the polarizabilities. The atomic multipole polarizabilities are calculated for different sets of molecules, linear alkanes, water clusters, and small organic molecules with different functional groups. It is found that the atomic and group contributions of the dipole and quadrupole polarizabilities are transferable as a function of the functional groups.     

Accuracy of distributed multipoles and polarizabilities: comparison between the LoProp and MpProp models

Localized multipole moments up to the fifth moment as well as localized dipole polarizabilities are calculated with the MpProp and the newly developed LoProp methods for a total of 20 molecules, predominantly derived from amino acids. A comparison of electrostatic potentials calculated from the multipole expansion obtained by the two methods with ab initio results shows that both methods reproduce the electrostatic interaction with an elementary charge with a mean absolute error of approximately 1.5 kJ/mol at contact distance and less than 0.1 kJ/mol at distances 2 A further out when terms up to the octupole moments are included. The polarizabilities are tested with homogenous electric fields and are found to have similar accuracy. The MpProp method gives better multipole moments unless diffuse basis sets are used, whereas LoProp gives better polarizabilities.     

Discrete and continuum contributions to multipole polarizabilities and shielding factors of hydrogen

Discrete and continuum contributions to the multipole polarizabilities and shielding factors of atomic hydrogen are computed. The discrete series show logarithmic convergence which can be accelerated by the u-transform. The continuum contributions increase with increasing multipole order and are already dominant in the quadrupole polarizability and shielding factor. The shielding factors have greater continuum contributions than the polarizabilities.     

Efficient determination of accurate atomic polarizabilities for polarizeable embedding calculations

We evaluate embedding potentials, obtained via various methods, used for polarizable embedding computations of excitation energies of para‐nitroaniline in water and organic solvents as well as of the green fluorescent protein. We found that isotropic polarizabilities derived from DFTD3 dispersion coefficients correlate well with those obtained via the LoProp method. We show that these polarizabilities in conjunction with appropriately derived point charges are in good agreement with calculations employing static multipole moments up to quadrupoles and anisotropic polarizabilities for both computed systems. The (partial) use of these easily‐accessible parameters drastically reduces the computational effort to obtain accurate embedding potentials especially for proteins. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.     

Gaussian induced dipole polarization model

A new induced dipole polarization model based on interacting Gaussian charge densities is presented. In contrast to the original induced point dipole model, the Gaussian polarization model is capable of finite interactions at short distances. Aspects of convergence related to the Gaussian model will be explored. The Gaussian polarization model is compared with the damped Thole-induced dipole model and the point dipole model. It will be shown that the Gaussian polarization model performs slightly better than the Thole model in terms of fitting to molecular polarizability tensors. An advantage of the model based on Gaussian charge distribution is that it can be easily generalized to other multipole moments and provide effective damping for both permanent electrostatic and polarization models. Finally, a method of parameterizing polarizabilities is presented. This method is based on probing a molecule with point charges and fitting polarizabilities to electrostatic potential. In contrast to the generic atom type polarizabilities fit to molecular polarizability tensors, probed polarizabilities are significantly more accurate in terms of reproducing molecular polarizability tensors and electrostatic potential, while retaining conformational transferability.     

Multipole polarizabilities and shielding factors of the hydrogen atom from the hydrodynamic analogy to quantum mechanics. III

Variational iterative solutions for multipole polarizabilities of the H atom have been carried out. The resulting numerical values of the Cauchy coefficients are in harmony with our recent work. The large frequency expansion of the multipole polarizabilities of the H atom is also calculated. The first term in such expansions, the multipole oscillator strength sum rule, is in harmony with earlier work.     

Multipole polarizabilities of Ar

The four lowest multipole polarizabilities of Ar have been calculated by using the complete fourth-order many-body perturbation theory approach and a large GTO/CGTO basis set including a number of diffuse and polarization functions. The present results for the dipole polarizability (α =11.23 au), quadrupole polarizability (C= 26.79 au), dipole-quadrupole polarizability (B = -164.3), and the dipole hyperpolarizability (γ =1329 au) are compared with other theoretical data and experimental values.     

International journal of quantum chemistry—a journal devoted to quantum theory and computations in chemistry,condensed matter physics,and biology. Editorial—program and policies

A procedure based on the polarization propagator technique is used to determine the longitudinal asymptotic electric polarizabilities per unit cell of infinite periodic systems. They are computed ab initio at the random-phase approximation level of accuracy for infinite model hydrogen chains using the STO -3G minimal basis set. This work also focuses on the effect of the number of interacting cells taken into account in the cell index summations as well as the number of k-points needed to perform the integration in the first Brillouin zone. Long-range effects are shown to influence the desired response properties more strongly than do short-range effects. Our direct method is presented as the only way to get the asymptotic longitudinal polarizabilities per unit cell of the systems presenting the largest polarizabilities, the most interesting polymers. © 1993 John Wiley & Sons, Inc.     

Linear-scaling formation of Kohn-Sham Hamiltonian: application to the calculation of excitation energies and polarizabilities of large molecular systems

We present calculations of excitation energies and polarizabilities in large molecular systems at the local-density and generalized-gradient approximation levels of density-functional theory (DFT). Our results are obtained using a linear-scaling DFT implementation in the program system DALTON for the formation of the Kohn-Sham Hamiltonian. For the Coulomb contribution, we introduce a modification of the fast multipole method to calculations over Gaussian charge distributions. It affords a simpler implementation than the original continuous fast multipole method by partitioning the electrostatic Coulomb interactions into "classical" and "nonclassical" terms which are explicitly evaluated by linear-scaling multipole techniques and a modified two-electron integral code, respectively. As an illustration of the code, we have studied the singlet and triplet excitation energies as well as the static and dynamic polarizabilities of polyethylenes, polyenes, polyynes, and graphite sheets with an emphasis on the trends observed with system size.     

Dispersion interaction between helium pair from hydrodynamic analogy to quantum mechanics

Frequency-dependent multipole polarizabilities of the He sequence have been calculated from a hydrodynamic model of quantum mechanics and using an independent-particle model Hamiltonian. Our present scheme is parallel to the uncoupled Hartree–Fock approximation so our values of polarizabilities, multipole transition energies and dispersion force coefficients between the He? He pair are comparable with earlier works using the uncoupled Hartree–Fock approximation.     

First-principles calculation of local atomic polarizabilities

Common methods of determining atomic polarizabilities suffer from the inclusion of nonlocal effects such as charge polarization. A new method is described for determining fully ab initio atomic polarizabilities based on calculating the response of atomic multipoles to the local electrostatic potential. The localized atomic polarizabilities are then used to calculate induction energies that are compared to ab initio induction energies to test their usefulness in practical applications. These polarizabilities are shown to be an improvement over the corresponding molecular polarizabilities, in terms of both absolute accuracy and the convergence of the multipolar induction series. The transferability of localized polarizabilities for the alkane series is also discussed.     

Fundamental relationships in the theory of electric multipole moments and multipole polarizabilities in static fields

New derivations are given of equations relating molecular electric multipole moments and polarizabilities of general order to the electrostatic energy. The unabridged moment convention is shown to yield relatively simple relations between derivatives of the energy with respect to field gradients and the multipole moments and polarizabilities. Care is taken to distinguish various forms of these derivatives, and one form leads to a proof of a general symmetry of polarizability tensors with respect to permutations of rank indices. The condition of internal equilibrium is shown to be fundamental to the existence of this symmetry. The transformation of multipole moment and polarizability tensors under translation of the coordinate origin is expressed in relatively simple general form. The traceless multipole and polarizability tensors of Buckingham and McLean and Yoshimine are obtained as linear combinations of the unabridged tensors and their traces.     

How important are temperature effects for cluster polarizabilities?

State-of-the-art first-principle all-electron density functional theory calculations on small sodium clusters are performed to study the temperature dependency of their polarizabilities. For this purpose Born-Oppenheimer molecular dynamics simulations with more than 100,000 time steps (>200 ps) are recorded employing gradient corrected functionals in combination with a double-zeta valence polarization basis set. For each cluster 18 trajectories between 50 and 900 K are collected. The cluster polarizabilities are then calculated along these trajectories employing a triple-zeta valence polarization basis set augmented with field-induced polarization functions. The analysis of these calculations shows that the temperature dependency of the sodium cluster polarizabilities varies strongly with cluster size. For several clusters characteristic changes in the polarizability per atom as a function of temperature are observed. It is shown that the inclusion of finite temperature effects resolves the long-standing mismatch between calculated and measured sodium cluster polarizabilities.     

Electric field effects on the shielding constants of noble gases: a four-component relativistic Hartree-Fock study

Second derivatives of nuclear shielding constants with respect to an electric field, i.e., shielding polarizabilities, have been calculated for the noble gas atoms from helium to xenon. The calculations have been carried out using the four-component relativistic Hartree-Fock method. In order to assess the importance of the individual relativistic corrections, the shielding polarizabilities have also been calculated at the nonrelativistic Hartree-Fock level, with spin-orbit and scalar (Darwin and mass-velocity) effects having been established by perturbative methods. Electron correlation effects have been estimated using the second-order polarization propagator approach. The relativistic effects on the tensor components of the shielding polarizabilities are found to be larger and changing less regularly with the atomic number than for the shielding constant itself. However, there is a partial cancellation of the contributions to the parallel and perpendicular components of the shielding polarizability and as a consequence the mean shielding polarizability is far less affected than the individual components.     

Multipole polarizabilities of spherically confined hydrogen atom and positronium in screened Coulomb potential

The static multipole polarizabilities of the hydrogen atom interacting with a screened Coulomb potential confined in an impenetrable spherical box are calculated in the sum-over-states formalism. The system eigenenergies and wave functions are produced by employing the generalized pseudospectral method. High-precision results are obtained to resolve the discrepancies between previous predictions, and used to analyze the critical behavior of the system with varying the confinement radius and screening parameter. A scaling law of the multipole polarizabilities with respect to the nuclear charge and electron reduced mass is proposed. Based on the scaling law we extend the investigation to the spatially confined positronium atom in a screened Coulomb potential by utilizing a reduced one-electron model. The results are in good agreement with previous estimates. However, it is conjectured that when the confinement radius is small the confined positronium atom should be treated more reasonably by an effective two-electron model.     

Multipole moments for embedding potentials: Exploring different atomic allocation algorithms

Polarizable quantum mechanical (QM)/molecular mechanics (MM)‐embedding methods are currently among the most promising methods for computationally feasible, yet reliable, production calculations of localized excitations and molecular response properties of large molecular complexes, such as proteins and RNA/DNA, and of molecules in solution. Our aim is to develop a computational methodology for distributed multipole moments and their associated multipole polarizabilities which is accurate, computationally efficient, and with smooth convergence with respect to multipole order. As the first step toward this goal, we herein investigate different ways of obtaining distributed atom‐centered multipole moments that are used in the construction of the electrostatic part of the embedding potential. Our objective is methods that not only are accurate and computationally efficient, but which can be consistently extended with site polarizabilities including internal charge transfer terms. We present a new way of dealing with well‐known problems in relation to the use of basis sets with diffuse functions in conventional atomic allocation algorithms, avoiding numerical integration schemes. Using this approach, we show that the classical embedding potential can be systematically improved, also when using basis sets with diffuse functions, and that very accurate embedding potentials suitable for QM/MM embedding calculations can be acquired. © 2016 Wiley Periodicals, Inc.     

Variational method for calculating molecular multipole polarizabilities using atom centered extended trial functions

An improvement of the variational method for calculating the electronic multipole polarizabilities is proposed. This modification allows the computation of the polarizabilities at any point and the results are compatible with the relations existing for a change of origin. It is applied to H₂, HF, CO and N₂ by using SCF wavefunctions developed on a limited basis. The computed polarizabilities are systematically too large but this discrepancy is attributed to the fact that the ground state is too far from the exact wavefunction.E.R.A. au C.N.R.S. No. 22.