Energetic Information from Information-Theoretic Approach in Density Functional Theory as Quantitative Measures of Physicochemical Properties.

The Hohenberg-Kohn theorem of density functional theory (DFT) stipulates that energy is a universal functional of electron density in the ground state, so energy can be thought of having encoded essential information for the density. Based on this, we recently proposed to quantify energetic information within the framework of information-theoretic approach (ITA) of DFT (J. Chem. Phys. 2022, 157, 101103). In this study, we systematically apply energetic information to a variety of chemical phenomena to validate the use of energetic information as quantitative measures of physicochemical properties. To that end, we employed six ITA quantities such as Shannon entropy and Fisher information for five energetic densities, yielding twenty-six viable energetic information quantities. Then, they are applied to correlate with physicochemical properties of molecular systems, including chemical bonding, conformational stability, intermolecular interactions, acidity, aromaticity, cooperativity, electrophilicity, nucleophilicity, and reactivity. Our results show that different quantities of energetic information often behave differently for different properties but a few of them, such as Shannon entropy of the total kinetic energy density and information gain of the Pauli energy density, stand out and strongly correlate with several properties across different categories of molecular systems. These results suggest that they can be employed as quantitative measures of physicochemical properties. This work not only enriches the body of our knowledge about the relationship between energy and information, but also provides scores of newly introduced explicit density functionals to quantify physicochemical properties, which can serve as robust features for building machine learning models in future studies.

GBasis: A Python library for evaluating functions, functionals, and integrals expressed with Gaussian basis functions.

GBasis is a free and open-source Python library for molecular property computations based on Gaussian basis functions in quantum chemistry. Specifically, GBasis allows one to evaluate functions expanded in Gaussian basis functions (including molecular orbitals, electron density, and reduced density matrices) and to compute functionals of Gaussian basis functions (overlap integrals, one-electron integrals, and two-electron integrals). Unique features of GBasis include supporting evaluation and analytical integration of arbitrary-order derivatives of the density (matrices), computation of a broad range of (screened) Coulomb interactions, and evaluation of overlap integrals of arbitrary numbers of Gaussians in arbitrarily high dimensions. For circumstances where the flexibility of GBasis is less important than high performance, a seamless Python interface to the Libcint C package is provided. GBasis is designed to be easy to use, maintain, and extend following many standards of sustainable software development, including code-quality assurance through continuous integration protocols, extensive testing, comprehensive documentation, up-to-date package management, and continuous delivery. This article marks the official release of the GBasis library, outlining its features, examples, and development.

Spin-Polarized Conceptual Density Functional Theory from the Convex Hull.

We present a new, nonarbitrary, internally consistent, and unambiguous framework for spin-polarized conceptual density-functional theory (SP-DFT). We explicitly characterize the convex hull of energy, as a function of the number of electrons and their spin, as the only accessible ground states in spin-polarized density functional theory. Then, we construct continuous linear and quadratic models for the energy. The nondifferentiable linear model exactly captures the simplicial geometry of the complex hull about the point of interest and gives exact representations for the conceptual DFT reactivity indicators. The continuous quadratic energy model is the paraboloid of maximum curvature, which most tightly encloses the point of interest and neighboring vertices. The quadratic model is invariant to the choice of coordinate system (i.e., {N, S} vs {Nα, Nß}) and reduces to a sensible formulation of spin-free conceptual DFT in the appropriate limit. Using the quadratic model, we generalize the Parr functions {P+(r), P-(r)} (and their derivatives with respect to number of electrons) to this new spin-polarized framework, integrating the Parr function concept into the context of (spin-polarized) conceptual DFT, and extending it to include higher-order effects.

Grid: A Python library for molecular integration, interpolation, differentiation, and more.

Grid is a free and open-source Python library for constructing numerical grids to integrate, interpolate, and differentiate functions (e.g., molecular properties), with a strong emphasis on facilitating these operations in computational chemistry and conceptual density functional theory. Although designed, maintained, and released as a stand-alone Python library, Grid was originally developed for molecular integration, interpolation, and solving the Poisson equation in the HORTON and ChemTools packages. Grid is designed to be easy to use, extend, and maintain; this is why we use Python and adopt many principles of modern software development, including comprehensive documentation, extensive testing, continuous integration/delivery protocols, and package management. We leverage popular scientific packages, such as NumPy and SciPy, to ensure high efficiency and optimized performance in grid development. This article is the official release note of the Grid library showcasing its unique functionality and scope.

Flexible Ansatz for N-Body Perturbation Theory.

We propose a new perturbation theory framework that can be used to help with the projective solution of the Schrödinger equation for arbitrary wave functions. This Flexible Ansatz for N-body Perturbation Theory (FANPT) is based on our previously proposed Flexible Ansatz for the N-body Configuration Interaction (FANCI). We derive recursive FANPT expressions, including arbitrary orders in the perturbation hierarchy. We show that the FANPT equations are well-behaved across a wide range of conditions, including static correlation-dominated configurations and highly nonlinear wave functions.

Coupled cluster-inspired geminal wavefunctions.

Electron pairs have an illustrious history in chemistry, from powerful concepts to understanding structural stability and reactive changes to the promise of serving as building blocks of quantitative descriptions of the electronic structure of complex molecules and materials. However, traditionally, two-electron wavefunctions (geminals) have not enjoyed the popularity and widespread use of the more standard single-particle methods. This has changed recently, with a renewed interest in the development of geminal wavefunctions as an alternative to describing strongly correlated phenomena. Hence, there is a need to find geminal methods that are accurate, computationally tractable, and do not demand significant input from the user (particularly via cumbersome and often ill-behaved orbital optimization steps). Here, we propose new families of geminal wavefunctions inspired by the pair coupled cluster doubles ansatz. We present a new hierarchy of two-electron wavefunctions that extends the one-reference orbital idea to other geminals. Moreover, we show how to incorporate single-like excitations in this framework without leaving the quasiparticle picture. We explore the role of imposing seniority restrictions on these wavefunctions and benchmark these new methods on model strongly correlated systems.

Chiral Jahn-Teller Distortion in Quasi-Planar Boron Clusters.

In this work, we have observed that some chiral boron clusters (B16-, B20-, B24-, and B28-) can simultaneously have helical molecular orbitals and helical spin densities; these seem to be the first compounds discovered to have this intriguing property. We show that chiral Jahn-Teller distortion of quasi-planar boron clusters drives the formation of the helical molecular spin densities in these clusters and show that elongation/enhancement in helical molecular orbitals can be achieved by simply adding more building blocks via a linker. Aromaticity of these boron clusters is discussed. Chiral boron clusters may find potential applications in spintronics, such as molecular magnets.

The tale of HORTON: Lessons learned in a decade of scientific software development.

HORTON is a free and open-source electronic-structure package written primarily in Python 3 with some underlying C++ components. While HORTON's development has been mainly directed by the research interests of its leading contributing groups, it is designed to be easily modified, extended, and used by other developers of quantum chemistry methods or post-processing techniques. Most importantly, HORTON adheres to modern principles of software development, including modularity, readability, flexibility, comprehensive documentation, automatic testing, version control, and quality-assurance protocols. This article explains how the principles and structure of HORTON have evolved since we started developing it more than a decade ago. We review the features and functionality of the latest HORTON release (version 2.3) and discuss how HORTON is evolving to support electronic structure theory research for the next decade.

Fragment-Based Deep Learning for Simultaneous Prediction of Polarizabilities and NMR Shieldings of Macromolecules and Their Aggregates.

Simultaneous prediction of the molecular response properties, such as polarizability and the NMR shielding constant, at a low computational cost is an unresolved issue. We propose to combine a linear-scaling generalized energy-based fragmentation (GEBF) method and deep learning (DL) with both molecular and atomic information-theoretic approach (ITA) quantities as effective descriptors. In GEBF, the total molecular polarizability can be assembled as a linear combination of the corresponding quantities calculated from a set of small embedded subsystems in GEBF. In the new GEBF-DL(ITA) protocol, one can predict subsystem polarizabilities based on the corresponding molecular wave function (thus electron density and ITA quantities) and DL model rather than calculate them from the computationally intensive coupled-perturbed Hartree-Fock or Kohn-Sham equations and finally obtain the total molecular polarizability via a linear combination equation. As a proof-of-concept application, we predict the molecular polarizabilities of large proteins and protein aggregates. GEBF-DL(ITA) is shown to be as accurate enough as GEBF, with mean absolute percentage error <1%. For the largest protein aggregate (>4000 atoms), GEBF-DL(ITA) gains a speedup ratio of 3 compared with GEBF. It is anticipated that when more advanced electronic structure methods are used, this advantage will be more appealing. Moreover, one can also predict the NMR chemical shieldings of proteins with reasonably good accuracy. Overall, the cost-efficient GEBF-DL(ITA) protocol should be a robust theoretical tool for simultaneously predicting polarizabilities and NMR shieldings of large systems.

Some Recent Advances in Density-Based Reactivity Theory.

Establishing a chemical reactivity theory in density functional theory (DFT) language has been our intense research interest in the past two decades, exemplified by the determination of steric effect and stereoselectivity, evaluation of electrophilicity and nucleophilicity, identification of strong and weak interactions, and formulation of cooperativity, frustration, and principle of chirality hierarchy. In this Featured Article, we first overview the four density-based frameworks in DFT to appreciate chemical understanding, including conceptual DFT, use of density associated quantities, information-theoretic approach, and orbital-free DFT, and then present a few recent advances of these frameworks as well as new applications from our studies. To that end, we will introduce the relationship among these frameworks, determining the entire spectrum of interactions with Pauli energy derivatives, performing topological analyses with information-theoretic quantities, and extending the density-based frameworks to excited states. Applications to examine physiochemical properties in external electric fields and to evaluate polarizability for proteins and crystals are discussed. A few possible directions for future development are followed, with the special emphasis on its merger with machine learning.

Energy is not a convex function of particle number for r-k interparticle potentials with k > log34.

The energy of a many-particle system is not convex with respect to particle number for r-k interparticle repulsion potentials if k > log34 ≈ 1.262. With such potentials, some finite electronic systems have ionization potentials that are less than the electron affinity: they have negative band gap (chemical hardness). Although the energy may be a convex function of the number of electrons (for which k = 1), it suggests that finding an analytic proof of convexity will be very difficult. The bound on k is postulated to be tight. An apparent signature of non-convex behavior is that the Dyson orbital corresponding to the lowest-energy mode of electron attachment has a vanishingly small amplitude.

Accurate and Efficient Prediction of Post-Hartree-Fock Polarizabilities of Condensed-Phase Systems.

To accurately and efficiently predict the molecular response properties (such as polarizability) at post-Hartree-Fock levels for condensed-phase systems under periodic boundary conditions (PBC) is still an unaccomplished and ongoing task. We demonstrate that static isotropic polarizabilities can be cost-effectively predicted at post-Hartree-Fock levels by combining the linear-scaling generalized energy-based fragmentation (GEBF) and information-theoretic approach (ITA) quantities. In PBC-GEBF, the total molecular polarizability of an extended system is obtained as a linear combination of the corresponding quantities of a series of small embedded subsystems of several monomers. Here, we show that in the PBC-GEBF-ITA framework, one can obtain the molecular polarizabilities and establish linear relations to ITA quantities. Once these relations are established for smaller subsystems, one can predict the polarizabilities of larger subsystems directly from the molecular wavefunction (or electron density) via ITA quantities. Alternatively, one can determine the total molecular polarizability via a linear combination equation in PBC-GEBF. We have corroborated that this newly proposed PBC-GEBF-ITA protocol is much more efficient than the original PBC-GEBF approach but is not much less accurate and that this conclusion holds for both many-body perturbation theory and the coupled cluster calculations. Good efficiency and transferability of the PBC-GEBF-ITA protocol are demonstrated for periodic systems with several hundred atoms in a unit cell.

Topological analysis of information-theoretic quantities in density functional theory.

We have witnessed considerable research interest in the recent literature about the development and applications of quantities from the information-theoretic approach (ITA) in density functional theory. These ITA quantities are explicit density functionals, whose local distributions in real space are continuous and well-behaved. In this work, we further develop ITA by systematically analyzing the topological behavior of its four representative quantities, Shannon entropy, two forms of Fisher information, and relative Shannon entropy (also called information gain or Kullback-Leibler divergence). Our results from their topological analyses for 103 molecular systems provide new insights into bonding interactions and physiochemical properties, such as electrophilicity, nucleophilicity, acidity, and aromaticity. We also compare our results with those from the electron density, electron localization function, localized orbital locator, and Laplacian functions. Our results offer a new methodological approach and practical tool for applications that are especially promising for elucidating chemical bonding and reactivity propensity.

An information-theoretic approach to basis-set fitting of electron densities and other non-negative functions.

The numerical ill-conditioning associated with approximating an electron density with a convex sum of Gaussian or Slater-type functions is overcome by using the (extended) Kullback-Leibler divergence to measure the deviation between the target and approximate density. The optimized densities are non-negative and normalized, and they are accurate enough to be used in applications related to molecular similarity, the topology of the electron density, and numerical molecular integration. This robust, efficient, and general approach can be used to fit any non-negative normalized functions (e.g., the kinetic energy density and molecular electron density) to a convex sum of non-negative basis functions. We present a fixed-point iteration method for optimizing the Kullback-Leibler divergence and compare it to conventional gradient-based optimization methods. These algorithms are released through the free and open-source BFit package, which also includes a L2-norm squared optimization routine applicable to any square-integrable scalar function.

Investigating the Linear Response Function under Approximations Following the Coupled-Perturbed Approach for Atoms and Molecules.

The linear response kernel also referred to as linear response function (LRF) in the framework of conceptual density functional theory has gained tremendous success in time-dependent density functional theory. Comparatively less attention has been devoted to the LRF from a chemical reactivity perspective in its time- or frequency-independent context, although it has recently been used to qualitatively describe electron delocalization, (anti-)aromaticity, inductive and mesomeric effects, etc. Despite these successes, which were obtained by approximating the LRF using the independent particle approximation deriving from a coupled-perturbed Kohn-Sham computation, the robustness of this LRF approach needs to be assessed. In this work, we compute the LRF at four levels of approximation (independent particle approximation, random phase approximation, Hartree-Fock approximation, and the (exact) DFT (density functional theory) expression) using functionals from the first four rungs of Jacob's ladder of exchange-correlation energy functionals. To scrutinize the impact of these approximations, new visualization strategies are discussed and systematized. The overall conclusion is that the independent particle approximation yields qualitatively correct results (ergo previous conceptual applications of the LRF are trustworthy), but for quantitative results, LRF expressions including coulomb and exchange(-correlation) terms should be included. With respect to functionals, density-gradient contributions to the exchange-correlation kernel are less than 10% and may be omitted safely where that is preferable computationally.

Excited-State Polarizabilities: A Combined Density Functional Theory and Information-Theoretic Approach Study.

Accurate and efficient determination of excited-state polarizabilities (α) is an open problem both experimentally and computationally. Following our previous work, (Phys. Chem. Chem. Phys. 2023, 25, 2131-2141), in which we employed simple ground-state (S0) density-related functions from the information-theoretic approach (ITA) to accurately and efficiently evaluate the macromolecular polarizabilities, in this work we aimed to predict the lowest excited-state (S1) polarizabilities. The philosophy is to use density-based functions to depict excited-state polarizabilities. As a proof-of-principle application, employing 2-(2'-hydroxyphenyl)benzimidazole (HBI), its substituents, and some other commonly used ESIPT (excited-state intramolecular proton transfer) fluorophores as model systems, we verified that either with S0 or S1 densities as an input, ITA quantities can be strongly correlated with the excited-state polarizabilities. When transition densities are considered, both S0 and S1 polarizabilities are in good relationships with some ITA quantities. The transferability of the linear regression model is further verified for a series of molecules with little or no similarity to those molecules in the training set. Furthermore, the excitation energies can be predicted based on multivariant linear regression equations of ITA quantities. This study also found that the nature of both the ground-state and excited-state polarizabilities of these species are due to the spatial delocalization of the electron density.

Fanpy: A python library for prototyping multideterminant methods in ab initio quantum chemistry.

Fanpy is a free and open-source Python library for developing and testing multideterminant wavefunctions and related ab initio methods in electronic structure theory. The main use of Fanpy is to quickly prototype new methods by making it easier to convert the mathematical formulation of a new wavefunction ansätze to a working implementation. Fanpy is designed based on our recently introduced Flexible Ansatz for N-electron Configuration Interaction (FANCI) framework, where multideterminant wavefunctions are represented by their overlaps with Slater determinants of orthonormal spin-orbitals. In the simplest case, a new wavefunction ansatz can be implemented by simply writing a function for evaluating its overlap with an arbitrary Slater determinant. Fanpy is modular in both implementation and theory: the wavefunction model, the system's Hamiltonian, and the choice of objective function are all independent modules. This modular structure makes it easy for users to mix and match different methods and for developers to quickly explore new ideas. Fanpy is written purely in Python with standard dependencies, making it accessible for various operating systems. In addition, it adheres to principles of modern software development, including comprehensive documentation, extensive testing, quality assurance, and continuous integration and delivery protocols. This article is considered to be the official release notes for the Fanpy library.

Efficient and accurate density-based prediction of macromolecular polarizabilities.

Accurately and efficiently predicting macromolecules' polarizabilities is an open problem. In this work, we employ a few simple density-based quantities from the information-theoretic approach (ITA) to predict polarizability of proteins. We first build quantitative structure/property relationships between molecular polarizabilities and ITA quantities. We then verify the broad applicability of ITA quantities for polarizability prediction for inorganic, organic, and biological systems with both localized and delocalized electronic structure. As a proof-of-concept application, we predict the molecular polarizabilities of complex proteins. Based on the linear regression equations for 20 natural amino acid residues, 400 dipeptides, and 8000 tripeptides, one then predicts the molecular polarizability of a larger peptide or even a protein once the molecular wavefunction is obtained. Because it is extremely costly to determine the wavefunction for a macromolecule like a protein, we propose to combine the ITA with the linear-scaling generalized energy-based fragmentation (GEBF) method to predict the macromolecular polarizability. In GEBF, the total molecular polarizability is obtained as a linear combination of the corresponding quantities from a series of small subsystems. We can predict them based on the subsystem wavefunction and linear regression equations rather than compute them from the nearly-intractable coupled-perturbed Hartree-Fock or Kohn-Sham equations for the whole macromolecule. Computational results showcase that the GEBF-ITA protocol should be an inexpensive yet accurate theoretical tool for predicting macromolecular polarizabilities.

Uranium(iv) alkyl cations: synthesis, structures, comparison with thorium(iv) analogues, and the influence of arene-coordination on thermal stability and ethylene polymerization activity.

Reaction of [(XA2)U(CH2SiMe3)2] (1; XA2 = 4,5-bis(2,6-diisopropylanilido)-2,7-di-tert-butyl-9,9-dimethylxanthene) with 1 equivalent of [Ph3C][B(C6F5)4] in arene solvents afforded the arene-coordinated uranium alkyl cations, [(XA2)U(CH2SiMe3)(η n -arene)][B(C6F5)4] {arene = benzene (2), toluene (3), bromobenzene (4) and fluorobenzene (5)}. Compounds 2, 3, and 5 were crystallographically characterized, and in all cases the arene is π-coordinated. Solution NMR studies of 2-5 suggest that the binding preferences of the [(XA2)U(CH2SiMe3)]+ cation follow the order: toluene ≈ benzene > bromobenzene > fluorobenzene. Compounds 2-4 generated in C6H5R (R = H, Me or Br, respectively) showed no polymerization activity under 1 atm of ethylene. By contrast, 5 and 5-Th (the thorium analogue of 5) in fluorobenzene at 20 and 70 °C achieved ethylene polymerization activities between 16 800 and 139 200 g mol-1 h-1 atm-1, highlighting the extent to which common arene solvents such as toluene can suppress ethylene polymerization activity in sterically open f-element complexes. However, activation of [(XA2)An(CH2SiMe3)2] {M = U (1) or Th (1-Th)} with [Ph3C][B(C6F5)4] in n-alkane solvents did not afford an active polymerization catalyst due to catalyst decomposition, illustrating the critical role of PhX (X = H, Me, Br or F) coordination for alkyl cation stabilization. Gas phase DFT calculations, including fragment interaction calculations with energy decomposition and ETS-NOCV analysis, were carried out on the cationic portion of 2'-Th, 2', 3' and 5' (analogues of 2-Th, 2, 3 and 5 with hydrogen atoms in place of ligand backbone methyl and tert-butyl groups), providing insight into the nature of actinide-arene bonding, which decreases in strength in the order 2'-Th > 2' ≈ 3' > 5'.

DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science.

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.