Alex Borgoo

Quantum similarity study of atomic density functions: Insights from information theory and the role of relativistic effects

A novel quantum similarity measure (QSM) is constructed based on concepts from information theory. In an application of QSM to atoms, the new QSM and its corresponding quantum similarity index (QSI) are evaluated throughout the periodic table, using the atomic electron densities and shape functions calculated in the Hartree-Fock approximation. The periodicity of Mendeleevs table is regained for the first time through the evaluation of a QSM. Evaluation of the information theory based QSI demonstrates, however, that the patterns of periodicity are lost due to the renormalization of the QSM, yielding chemically less appealing results for the QSI. A comparison of the information content of a given atom on top of a group with the information content of the elements in the subsequent rows reveals another periodicity pattern. Relativistic effects on the electronic density functions of atoms are investigated. Their importance is quantified in a QSI study by comparing for each atom, the density functions evaluated in the Hartree-Fock and Dirac-Fock approximations. The smooth decreasing of the relevant QSI along the periodic table illustrates in a quantitative way the increase of relativistic corrections with the nuclear charge.

Information theoretical study of chirality: enantiomers with one and two asymmetric centra.

In this work, the Kullback-Leibler information deficiency is probed as a chirality measure. It is argued that the information deficiency, calculated using the shape functions of the R and S enantiomers, considering one as reference for the other, gives an information theory based expression useful for quantifying chirality. The measure is evaluated for five chiral halomethanes possessing one asymmetric carbon atom with hydrogen, fluorine, chlorine, bromine, and iodine as substituents. To demonstrate the general applicability, a study of two halogen-substituted ethanes possessing two asymmetric carbon atoms has been included as well. The basic expression of the sum of the local information deficiency over all atoms can be decomposed into separate summations over coinciding and noncoinciding atoms, or into a global and a mixing entropy term, or into a local entropy contribution for each atom individually based on the Hirshfeld partitioning. Avnirs continuous chirality measure (CCM) has been computed and confronted with the information deficiency. Finally, the relationship between chirality and optical rotation is used to study the proposed measure. The results illustrate Mezeys holographic electron density theorem with an intuitively appealing division of the strength of propagation of the atomic chirality from an asymmetric carbon atom throughout the molecule. The local information deficiency of the carbon atom is proposed as a measure of chirality; more precisely, the difference in information between the R and the S enantiomer turns out to be a quantitative measure of the chirality of the system. It may be evaluated as the arithmetic mean of the different alignments, or considering only the alignment resulting in the highest similarity value, or using the QSSA alignment.

Effective homogeneity of the exchange-correlation and non-interacting kinetic energy functionals under density scaling.

Correlated electron densities, experimental ionisation potentials, and experimental electron affinities are used to investigate the homogeneity of the exchange-correlation and non-interacting kinetic energy functionals of Kohn-Sham density functional theory under density scaling. Results are presented for atoms and small molecules, paying attention to the influence of the integer discontinuity and the choice of the electron affinity. For the exchange-correlation functional, effective homogeneities are highly system-dependent on either side of the integer discontinuity. By contrast, the average homogeneity-associated with the potential that averages over the discontinuity-is generally close to 4/3 when the discontinuity is computed using positive affinities for systems that do bind an excess electron and negative affinities for those that do not. The proximity to 4/3 becomes increasingly pronounced with increasing atomic number. Evaluating the discontinuity using a zero affinity in systems that do not bind an excess electron instead leads to effective homogeneities on the electron abundant side that are close to 4/3. For the non-interacting kinetic energy functional, the effective homogeneities are less system-dependent and the effect of the integer discontinuity is less pronounced. Average values are uniformly below 5/3. The study provides information that may aid the development of improved exchange-correlation and non-interacting kinetic energy functionals.

Density Scaling of Noninteracting Kinetic Energy Functionals

The influence of imposing an approximate density scaling condition on a noninteracting kinetic energy functional is investigated. A simple generalized gradient approximation (GGA) is presented, which satisfies both the density scaling condition and the usual coordinate scaling condition; the remaining multiplicative constant is determined from an energy criterion. In post-Kohn-Sham calculations, noninteracting kinetic energies of the closed-shell molecules of the G1 set determined using the GGA are a modest improvement over those determined using the corresponding local functional, which does not satisfy the density scaling condition. Potential energy curves of CO, F2, and P2 exhibit binding with the GGA, compared to purely repulsive curves with the local functional. Adjusting the exponent in the GGA form in order to optimize energy accuracy violates the density scaling condition, and two of the diatomics no longer exhibit binding. Results are compared with those from other local/GGA functionals in the literature.

Negative Electron Affinities from DFT: Influence of Asymptotic Exchange-Correlation Potential and Effective Homogeneity under Density Scaling

The influence of the asymptotic exchange-correlation potential and density-scaling homogeneity on negative electron affinities determined using the approach of Tozer and De Proft [J. Phys. Chem. A2005, 109, 8923] is investigated. Application of an asymptotic correction to the potential improves the accuracy for several of the systems with the most negative affinities, reflecting their diffuse lowest unoccupied orbitals. For systems with modest affinities, it reduces the accuracy marginally. Enforcing a near-exact effective homogeneity through a simple shift in the potential leads to improved correlation with experimental values but significantly overestimated affinities. Optimal effective homogeneities are therefore determined, and a simple scheme is proposed for enforcing an average optimal value. Application of the scheme to a series of organic molecules maintains the excellent correlation with the experimental values while significantly reducing the absolute errors.

Molecular Binding in Post-Kohn-Sham Orbital-Free DFT.

Molecular binding in post-Kohn-Sham orbital-free DFT is investigated, using noninteracting kinetic energy functionals that satisfy the uniform electron gas condition and which are inhomogeneous under density scaling. A parameter is introduced that quantifies binding, and a series of functionals are determined from fits to near-exact effective homogeneities and/or Kohn-Sham noninteracting kinetic energies. These are then used to investigate the relationship between binding and the accuracy of the effective homogeneity and noninteracting kinetic energy at the equilibrium geometry. For a series of 11 molecules, the binding broadly improves as the effective homogeneity improves, although the extent to which it improves is dependent on the accuracy of the noninteracting kinetic energy; optimal binding appears to require both to be accurate simultaneously. The use of a Thomas-Fermi-von Weizsäcker form, augmented with a second gradient correction, goes some way toward achieving this, exhibiting molecular binding on average. The findings are discussed in terms of the noninteracting kinetic potential and the Hellmann-Feynman theorem. The extent to which the functionals can reproduce the system-dependence of the near-exact effective homogeneity is quantified, and potential energy curves are presented for selected molecules. The study provides impetus for including density scaling homogeneity considerations in the design of noninteracting kinetic energy functionals.

Alternative Kullback−Leibler Information Entropy for Enantiomers

In our series of studies on quantifying chirality, a new chirality measure is proposed in this work based on the Kullback-Leibler information entropy. The index computes the extra information that the shape function of one enantiomer carries over a normalized shape function of the racemate, while in our previous studies the shape functions of the R and S enantiomers were used considering one as reference for the other. Besides being mathematically more elegant (symmetric, positive definite, zero in the case of a nonchiral system), this new index bears a more direct relation with chirality oriented experimental measurements such as circular dichroism (CD) and optical rotation measurements, where the racemate is frequently used as a reference. The five chiral halomethanes holding one asymmetric carbon atom and H, F, Cl, Br, and I as substituents have been analyzed. A comparison with our calculated optical rotation and with Avnirs Continuous Chirality Measure (CCM) is computed. The results show that with this index the emphasis lies on the differences between the noncoinciding substituents.

Charge-constrained auxiliary-density-matrix methods for the Hartree–Fock exchange contribution

Three new variants of the auxiliary-density-matrix method (ADMM) of Guidon, Hutter, and VandeVondele [J. Chem. Theory Comput. 6, 2348 (2010)] are presented with the common feature that they have a simplified constraint compared with the full orthonormality requirement of the earlier ADMM1 method. All ADMM variants are tested for accuracy and performance in all-electron B3LYP calculations with several commonly used basis sets. The effect of the choice of the exchange functional for the ADMM exchange-correction term is also investigated.

Magnetic-Field Density-Functional Theory (BDFT): Lessons from the Adiabatic Connection

We study the effects of magnetic fields in the context of magnetic field density-functional theory (BDFT), where the energy is a functional of the electron density ρ and the magnetic field B. We show that this approach is a worthwhile alternative to current-density functional theory (CDFT) and may provide a viable route to the study of many magnetic phenomena using density-functional theory (DFT). The relationship between BDFT and CDFT is developed and clarified within the framework of the four-way correspondence of saddle functions and their convex and concave parents in convex analysis. By decomposing the energy into its Kohn-Sham components, we demonstrate that the magnetizability is mainly determined by those energy components that are related to the density. For existing density functional approximations, this implies that, for the magnetizability, improvements of the density will be more beneficial than introducing a magnetic-field dependence in the correlation functional. However, once a good charge density is achieved, we show that high accuracy is likely only obtainable by including magnetic-field dependence. We demonstrate that adiabatic-connection (AC) curves at different field strengths resemble one another closely provided each curve is calculated at the equilibrium geometry of that field strength. In contrast, if all AC curves are calculated at the equilibrium geometry of the field-free system, then the curves change strongly with increasing field strength due to the increasing importance of static correlation. This holds also for density functional approximations, for which we demonstrate that the main error encountered in the presence of a field is already present at zero field strength, indicating that density-functional approximations may be applied to systems in strong fields, without the need to treat additional static correlation.

Kohn–Sham energy decomposition for molecules in a magnetic field

ABSTRACT We study the total molecular electronic energy and its Kohn–Sham components within the framework of magnetic-field density-functional theory (BDFT), an alternative to current-dependent density-functional theory (CDFT) for molecules in the presence of magnetic fields. For a selection of closed-shell dia- and paramagnetic molecules, we investigate the dependence of the total electronic energy and its Kohn–Sham components on the magnetic field. Results obtained from commonly used density-functional approximations are compared with those obtained from Lieb optimisations based on magnetic-field dependent relaxed coupled-cluster singles-and-doubles (CCSD) and second-order Møller–Plesset (MP2) densities. We show that popular approximate exchange–correlation functionals at the generalised-gradient-approximation (GGA), meta-GGA, and hybrid levels of theory provide a good qualitative description of the electronic energy and its Kohn–Sham components in a magnetic field–in particular, for the diamagnetic molecules. The performance of Hartree–Fock theory, MP2 theory, CCSD theory and BDFT with different exchange–correlation functionals is compared with coupled-cluster singles-doubles-perturbative-triples (CCSD(T)) theory for the perpendicular component of the magnetisability. Generalisations of the TPSS meta-GGA functional to systems in a magnetic field work well–the cTPSS functional, in particular, with a current-corrected kinetic-energy density, performs excellently, providing an accurate and balanced treatment of dia- and paramagnetic systems and outperforming MP2 theory. GRAPHICAL ABSTRACT