Density-Based Partitioning Methods for Ground-State Molecular Calculations

Calculating Molecular Polarizabilities Using Exact Frozen Density Embedding with External Orthogonality

Frozen density embedding (FDE) with freeze-thaw cycles is a formally exact embedding scheme. In practice, this method is limited to systems with small density overlaps when approximate nonadditive kinetic energy functionals are used. It has been shown that the use of approximate nonadditive kinetic energy functionals can be avoided when external orthogonality (EO) is enforced, and FDE can then generate exact results even for strongly overlapping subsystems. In this work, we present an implementation of exact FDEc-EO (coupled FDE TDDFT with EO) for the calculation of polarizabilities in the Amsterdam density functional program package. EO is enforced using the level-shift projection operator method, which ensures that orbitals between fragments are orthogonal. For pure functionals, we show that only the symmetric EO contributions to the induced density matrix are needed. This leads to a simplified implementation for the calculation of polarizability that can exactly reproduce the supermolecular TDDFT results. We further discuss the limitation of exact FDEc-EO in interpreting subsystem polarizabilities due to the nonunique partitioning of the total density. We show that this limitation is due to the fact that subsystem polarizability partitioning is dependent on how the subsystems are initially polarized. As supermolecular virtual orbitals are exactly reproduced, this dependence is attributed to the description of the occupied orbitals. In contrast, for excitations of subsystems that are localized within one subsystem, we show that the excitation energies are stable with respect to the order of polarization. This observation shows that impacts from the nonunique nature of exact FDE on subsystem properties can be minimized by better fragmentation of the supermolecular systems if the property is localized. For global properties like polarizability, this is not the case, and nonuniqueness remains independent of the fragmentation used.

Stretching Bonds without Breaking Symmetries in Density Functional Theory

Kohn-Sham density functional theory (KS-DFT) stands out among electronic structure methods due to its balance of accuracy and computational efficiency. However, to achieve chemically accurate energies, standard density functional approximations in KS-DFT often need to break underlying symmetries, a long-standing "symmetry dilemma". By employing fragment spin densities as the main variables in calculations (rather than total molecular densities, as in KS-DFT), we present an embedding framework in which this symmetry dilemma is understood and partially resolved. The spatial overlap between fragment densities is used as the main ingredient to construct a simple, physically motivated approximation to a universal functional of the fragment densities. This "overlap approximation" is shown to significantly improve semilocal KS-DFT binding energies of molecules without artificially breaking either charge or spin symmetries. The approach is shown to be applicable to covalently bonded molecules and to systems of the "strongly correlated" type.

Orbital-Free Density Functional Theory: An Attractive Electronic Structure Method for Large-Scale First-Principles Simulations

Kohn-Sham Density Functional Theory (KSDFT) is the most widely used electronic structure method in chemistry, physics, and materials science, with thousands of calculations cited annually. This ubiquity is rooted in the favorable accuracy vs cost balance of KSDFT. Nonetheless, the ambitions and expectations of researchers for use of KSDFT in predictive simulations of large, complicated molecular systems are confronted with an intrinsic computational cost-scaling challenge. Particularly evident in the context of first-principles molecular dynamics, the challenge is the high cost-scaling associated with the computation of the Kohn-Sham orbitals. Orbital-free DFT (OFDFT), as the name suggests, circumvents entirely the explicit use of those orbitals. Without them, the structural and algorithmic complexity of KSDFT simplifies dramatically and near-linear scaling with system size irrespective of system state is achievable. Thus, much larger system sizes and longer simulation time scales (compared to conventional KSDFT) become accessible; hence, new chemical phenomena and new materials can be explored. In this review, we introduce the historical contexts of OFDFT, its theoretical basis, and the challenge of realizing its promise via approximate kinetic energy density functionals (KEDFs). We review recent progress on that challenge for an array of KEDFs, such as one-point, two-point, and machine-learnt, as well as some less explored forms. We emphasize use of exact constraints and the inevitability of design choices. Then, we survey the associated numerical techniques and implemented algorithms specific to OFDFT. We conclude with an illustrative sample of applications to showcase the power of OFDFT in materials science, chemistry, and physics.

A Perspective on Sustainable Computational Chemistry Software Development and Integration

The power of quantum chemistry to predict the ground and excited state properties of complex chemical systems has driven the development of computational quantum chemistry software, integrating advances in theory, applied mathematics, and computer science. The emergence of new computational paradigms associated with exascale technologies also poses significant challenges that require a flexible forward strategy to take full advantage of existing and forthcoming computational resources. In this context, the sustainability and interoperability of computational chemistry software development are among the most pressing issues. In this perspective, we discuss software infrastructure needs and investments with an eye to fully utilize exascale resources and provide unique computational tools for next-generation science problems and scientific discoveries.

Split electrons in partition density functional theory

Partition density functional theory is a density embedding method that partitions a molecule into fragments by minimizing the sum of fragment energies subject to a local density constraint and a global electron-number constraint. To perform this minimization, we study a two-stage procedure in which the sum of fragment energies is lowered when electrons flow from fragments of lower electronegativity to fragments of higher electronegativity. The global minimum is reached when all electronegativities are equal. The non-integer fragment populations are dealt with in two different ways: (1) An ensemble approach (ENS) that involves averaging over calculations with different numbers of electrons (always integers) and (2) a simpler approach that involves fractionally occupying orbitals (FOO). We compare and contrast these two approaches and examine their performance in some of the simplest systems where one can transparently apply both, including simple models of heteronuclear diatomic molecules and actual diatomic molecules with two and four electrons. We find that, although both ENS and FOO methods lead to the same total energy and density, the ENS fragment densities are less distorted than those of FOO when compared to their isolated counterparts, and they tend to retain integer numbers of electrons. We establish the conditions under which the ENS populations can become fractional and observe that, even in those cases, the total charge transferred is always lower in ENS than in FOO. Similarly, the FOO fragment dipole moments provide an upper bound to the ENS dipoles. We explain why and discuss the implications.

Numerically stable inversion approach to construct Kohn-Sham potentials for given electron densities within a Gaussian basis set framework

We present a Kohn-Sham (KS) inversion approach to construct KS exchange-correlation potentials corresponding to given electron densities. This method is based on an iterative procedure using linear response to update potentials. All involved quantities, i.e., orbitals, potentials, and response functions, are represented by Gaussian basis functions. In contrast to previous KS inversion methods relying on Gaussian basis sets, the method presented here is numerically stable even for standard basis sets from basis set libraries due to a preprocessing of the auxiliary basis used to represent an exchange-correlation charge density that generates the exchange-correlation potential. The new KS inversion method is applied to reference densities of various atoms and molecules obtained by full configuration interaction or CCSD(T) (coupled cluster singles doubles perturbative triples). The considered examples encompass cases known to be difficult, such as stretched hydrogen or lithium hydride molecules or the beryllium isoelectronic series. For the stretched hydrogen molecule, potentials of benchmark quality are obtained by employing large basis sets. For the carbon monoxide molecule, we show that the correlation potential from the random phase approximation (RPA) is in excellent qualitative and quantitative agreement with the correlation potential from the KS inversion of a CCSD(T) reference density. This indicates that RPA correlation potentials, in contrast to those from semi-local density-functionals, resemble the exact correlation potential. Besides providing exchange-correlation potentials for benchmark purposes, the proposed KS inversion method may be used in density-partition-based quantum embedding and in subsystem density-functional methods because it combines numerical stability with computational efficiency.

Second-Order Active-Space Embedding Theory

Quantum embedding schemes are a promising way to extend multireference computations to large molecules with strong correlation effects localized on a small number of atoms. This work introduces a second-order active-space embedding theory [ASET(2)] which improves upon mean-field frozen embedding by treating fragment-environment interactions via an approximate canonical transformation. The canonical transformation employed in ASET(2) is formulated using the driven similarity renormalization group. The ASET(2) scheme is benchmarked on the N═N bond dissociation in pentyldiazene, the S₀ to S₁ excitation in 1-octene, and the interaction energy of the O₂-benzene complex. The ASET(2) explicit treatment of fragment-environment interactions beyond the mean-field level generally improves the accuracy of embedded computations, and it becomes necessary to achieve an accurate description of excitation energies of 1-octene and the singlet-triplet gap of the O₂-benzene complex.

Surface science using coupled cluster theory via local Wannier functions and in-RPA-embedding: The case of water on graphitic carbon nitride

A first-principles study of the adsorption of a single water molecule on a layer of graphitic carbon nitride is reported employing an embedding approach for many-electron correlation methods. To this end, a plane-wave based implementation to obtain intrinsic atomic orbitals and Wannier functions for arbitrary localization potentials is presented. In our embedding scheme, the localized occupied orbitals allow for a separate treatment of short-range and long-range correlation contributions to the adsorption energy by a fragmentation of the simulation cell. In combination with unoccupied natural orbitals, the coupled cluster ansatz with single, double, and perturbative triple particle-hole excitation operators is used to capture the correlation in local fragments centered around the adsorption process. For the long-range correlation, a seamless embedding into the random phase approximation yields rapidly convergent adsorption energies with respect to the local fragment size. Convergence of computed binding energies with respect to the virtual orbital basis set is achieved employing a number of recently developed techniques. Moreover, we discuss fragment size convergence for a range of approximate many-electron perturbation theories. The obtained benchmark results are compared to a number of density functional calculations.

Direct orbital selection within the domain-based local pair natural orbital coupled-cluster method

Domain-based local pair natural orbital coupled cluster (DLPNO-CC) has become increasingly popular to calculate relative energies (e.g., reaction energies and reaction barriers). It can be applied within a multi-level DLPNO-CC-in-DLPNO-CC ansatz to reduce the computational cost and focus the available computational resources on a specific subset of the occupied orbitals. We demonstrate how this multi-level DLPNO-CC ansatz can be combined with our direct orbital selection (DOS) approach [M. Bensberg and J. Neugebauer, J. Chem. Phys. 150, 214106 (2019)] to automatically select orbital sets for any multi-level calculation. We find that the parameters for the DOS procedure can be chosen conservatively such that they are transferable between reactions. The resulting automatic multi-level DLPNO-CC method requires no user input and is extremely robust and accurate. The computational cost is easily reduced by a factor of 3 without sacrificing accuracy. We demonstrate the accuracy of the method for a total of 61 reactions containing up to 174 atoms and use it to predict the relative stability of conformers of a Ru-based catalyst.

Projector-Based Quantum Embedding for Molecular Systems: An Investigation of Three Partitioning Approaches

Projector-based embedding is a relatively recent addition to the collection of methods that seek to utilize chemical locality to provide improved computational efficiency. This work considers the interactions between the different proposed procedures for this method and their effects on the accuracy of the results. The interplay between the embedded background, projector type, partitioning scheme, and level of atomic orbital (AO) truncation are investigated on a selection of reactions from the literature. The Huzinaga projection approach proves to be more reliable than the level-shift projection when paired with other procedural options. Active subsystem partitioning from the subsystem projected AO decomposition (SPADE) procedure proves slightly better than the combination of Pipek-Mezey localization and Mulliken population screening (PMM). Along with these two options, a new partitioning criteria is proposed based on subsystem von Neumann entropy and the related subsystem orbital occupancy. This new method overlaps with the previous PMM method, but the screening process is computationally simpler. Finally, AO truncation proves to be a robust option for the tested systems when paired with the Huzinaga projection, with satisfactory results being acquired at even the most severe truncation level.

From Kohn-Sham to Many-Electron Energies via Step Structures in the Exchange-Correlation Potential

Accurately describing excited states within Kohn–Sham (KS) density functional theory (DFT), particularly those which induce ionization and charge transfer, remains a great challenge. Common exchange-correlation (xc) approximations are unreliable for excited states owing, in part, to the absence of a derivative discontinuity in the xc energy (Δ), which relates a many-electron energy difference to the corresponding KS energy difference. We demonstrate, analytically and numerically, how the relationship between KS and many-electron energies leads to the step structures observed in the exact xc potential in four scenarios: electron addition, molecular dissociation, excitation of a finite system, and charge transfer. We further show that steps in the potential can be obtained also with common xc approximations, as simple as the LDA, when addressed from the ensemble perspective. The article therefore highlights how capturing the relationship between KS and many-electron energies with advanced xc approximations is crucial for accurately calculating excitations, as well as the ground-state density and energy of systems which consist of distinct subsystems.

Local embedding of coupled cluster theory into the random phase approximation using plane waves

We present an embedding approach to treat local electron correlation effects in periodic environments. In a single consistent framework, our plane wave based scheme embeds a local high-level correlation calculation [here, Coupled Cluster (CC) theory], employing localized orbitals, into a low-level correlation calculation [here, the direct Random Phase Approximation (RPA)]. This choice allows for an accurate and efficient treatment of long-range dispersion effects. Accelerated convergence with respect to the local fragment size can be observed if the low-level and high-level long-range dispersions are quantitatively similar, as is the case for CC in RPA. To demonstrate the capabilities of the introduced embedding approach, we calculate adsorption energies of molecules on a surface and in a chabazite crystal cage, as well as the formation energy of a lattice impurity in a solid at the level of highly accurate many-electron perturbation theories. The absorption energy of a methane molecule in a zeolite chabazite is converged with an error well below 20 meV at the CC level. As our largest periodic benchmark system, we apply our scheme to the adsorption of a water molecule on titania in a supercell containing more than 1000 electrons.

Electronic couplings for photo-induced processes from subsystem time-dependent density-functional theory: The role of the diabatization

Subsystem time-dependent density-functional theory (sTDDFT) making use of approximate non-additive kinetic energy (NAKE) functionals is known to be capable of describing excitation energy transfer processes in a variety of applications. Here, we show that sTDDFT, especially when combined with projection-based embedding (PbE), can be employed for the entire range of photo-induced electronic couplings essential for modeling photophysical properties of complex chemical and biological systems and therefore represents a complete toolbox for this class of problems. This means that it is capable of capturing the interaction/coupling associated with local- and charge-transfer (CT) excitons. However, this requires the choice of a reasonable diabatic basis. We therefore propose different diabatization strategies of the virtual orbital space in PbE-sTDDFT and show how CT excitations can be included in sTDDFT using NAKE functionals via a phenomenological approach. Finally, these electronic couplings are compared to couplings from a multistate fragment excitation difference (FED)-fragment charge difference (FCD) diabatization procedure. We show that both procedures, multistate FED-FCD and sTDDFT (with the right diabatization procedure chosen), lead to an overall good agreement for the electronic couplings, despite differences in their general diabatization strategy. We conclude that the entire range of photo-induced electronic couplings can be obtained using sTDDFT (with the right diabatization procedure chosen) in a black-box manner.

Nonlocal Subsystem Density Functional Theory

By invoking a divide-and-conquer strategy, subsystem DFT dramatically reduces the computational cost of large-scale, ab initio electronic structure simulations of molecules and materials. The central ingredient setting subsystem DFT apart from Kohn-Sham DFT is the nonadditive kinetic energy functional (NAKE). Currently employed NAKEs are at most semilocal (i.e., they only depend on the electron density and its gradient), and as a result of this approximation, so far large-scale simulations only included systems composed of weakly interacting subsystems. In this work, we advance the state-of-the-art by introducing fully nonlocal NAKEs in subsystem DFT simulations for the first time. A benchmark analysis based on the S22-5 test set shows that nonlocal NAKEs considerably improve the computed interaction energies and electron densities compared to commonly employed GGA NAKEs, especially when increasing intersubsystem electron density overlap is considered. Most importantly, we resolve the long-standing problem of too attractive interaction energy curves typically resulting from the use of GGA NAKEs.

Density functional theory based embedding approaches for transition-metal complexes

Transition metal species are commonly discussed by considering the metal atom embedded in a ligand environment. This apparently makes them interesting targets for modern embedding strategies based on Kohn-Sham density functional theory (DFT), which aim at modelling accurate predictions for large systems by combining different quantum chemical methods. In this perspective, we will focus on subsystem density functional theory and projection-based embedding. We review the developments in the field for transition metal species, demonstrate benefits, drawbacks and analyse error sources of the different strategies using the example of chromium hexacarbonyle, before giving a perspective where the field is currently heading.

Self-Consistent Density-Functional Embedding: A Novel Approach for Density-Functional Approximations

In the present work, we introduce a self-consistent density-functional embedding technique, which leaves the realm of standard energy-functional approaches in density functional theory and targets directly the density-to-potential mapping that lies at its heart. Inspired by the density matrix embedding theory, we project the full system onto a set of small interacting fragments that can be solved accurately. Based on the rigorous relation of density and potential in density functional theory, we then invert the fragment densities to local potentials. Combining these results in a continuous manner provides an update for the Kohn-Sham potential of the full system, which is then used to update the projection. We benchmark our approach for molecular bond stretching in one and two dimensions and show that, in these cases, the scheme converges to accurate approximations for densities and Kohn-Sham potentials. We demonstrate that the known steps and peaks of the exact exchange-correlation potential are reproduced by our method with remarkable accuracy.

Kohn-Sham potentials from electron densities using a matrix representation within finite atomic orbital basis sets

We revisit the static response function-based Kohn-Sham (KS) inversion procedure for determining the KS effective potential that corresponds to a given target electron density within finite atomic orbital basis sets. Instead of expanding the potential in an auxiliary basis set, we directly update the potential in its matrix representation. Through numerical examples, we show that the reconstructed density rapidly converges to the target density. Preliminary results are presented to illustrate the possibility of obtaining a local potential in real space from the optimized potential in its matrix representation. We have further applied this matrix-based KS inversion approach to density functional embedding theory. A proof-of-concept study of a solvated proton transfer reaction demonstrates the method's promise.

Non-additive non-interacting kinetic energy of rare gas dimers

Approximations of the non-additive non-interacting kinetic energy (NAKE) as an explicit functional of the density are the basis of several electronic structure methods that provide improved computational efficiency over standard Kohn-Sham calculations. However, within most fragment-based formalisms, there is no unique exact NAKE, making it difficult to develop general, robust approximations for it. When adjustments are made to the embedding formalisms to guarantee uniqueness, approximate functionals may be more meaningfully compared to the exact unique NAKE. We use numerically accurate inversions to study the exact NAKE of several rare-gas dimers within partition density functional theory, a method that provides the uniqueness for the exact NAKE. We find that the NAKE decreases nearly exponentially with atomic separation for the rare-gas dimers. We compute the logarithmic derivative of the NAKE with respect to the bond length for our numerically accurate inversions as well as for several approximate NAKE functionals. We show that standard approximate NAKE functionals do not reproduce the correct behavior for this logarithmic derivative and propose two new NAKE functionals that do. The first of these is based on a re-parametrization of a conjoint Perdew-Burke-Ernzerhof (PBE) functional. The second is a simple, physically motivated non-decomposable NAKE functional that matches the asymptotic decay constant without fitting.

Extending the density functional embedding theory to finite temperature and an efficient iterative method for solving for embedding potentials

A key element in the density functional embedding theory (DFET) is the embedding potential. We discuss two major issues related to the embedding potential: (1) its non-uniqueness and (2) the numerical difficulty for solving for it, especially for the spin-polarized systems. To resolve the first issue, we extend DFET to finite temperature: all quantities, such as the subsystem densities and the total system's density, are calculated at a finite temperature. This is a physical extension since materials work at finite temperatures. We show that the embedding potential is strictly unique at T > 0. To resolve the second issue, we introduce an efficient iterative embedding potential solver. We discuss how to relax the magnetic moments in subsystems and how to equilibrate the chemical potentials across subsystems. The solver is robust and efficient for several non-trivial examples, in all of which good quality spin-polarized embedding potentials were obtained. We also demonstrate the solver on an extended periodic system: iron body-centered cubic (110) surface, which is related to the modeling of the heterogeneous catalysis involving iron, such as the Fischer-Tropsch and the Haber processes. This work would make it efficient and accurate to perform embedding simulations of some challenging material problems, such as the heterogeneous catalysis and the defects of complicated spin configurations in electronic materials.

Fragment-based treatment of delocalization and static correlation errors in density-functional theory

One of the most important open challenges in modern Kohn-Sham (KS) density-functional theory (DFT) is the correct treatment of systems involving fractional electron charges and spins. Approximate exchange-correlation functionals struggle with such systems, leading to pervasive delocalization and static correlation errors. We demonstrate how these errors, which plague density-functional calculations of bond-stretching processes, can be avoided by employing the alternative framework of partition density-functional theory (PDFT) even using the local density approximation for the fragments. Our method is illustrated with explicit calculations on simple systems exhibiting delocalization and static-correlation errors, stretched H2 (+), H2, He2 (+), Li2 (+), and Li2. In all these cases, our method leads to greatly improved dissociation-energy curves. The effective KS potential corresponding to our self-consistent solutions displays key features around the bond midpoint; these are known to be present in the exact KS potential, but are absent from most approximate KS potentials and are essential for the correct description of electron dynamics.