Orbital-free Density Functional Theory Research Articles

Neural network learned Pauli potential for the advancement of orbital-free density functional theory.

The Pauli kinetic energy functional and its functional derivative, termed Pauli potential, play a crucial role in the successful implementation of orbital-free density functional theory for electronic structure calculations. However, the exact forms of these two quantities are not known. Therefore, perforce, one employs the approximate forms for the Pauli functional or Pauli potential for performing orbital-free density functional calculations. In the present study, we developed a feed-forward neural network-based representation for the Pauli potential using a 1-dimensional (1-D) model system. We expanded density in terms of basis functions, and the coefficients of the expansion were used as input to a feed-forward neural network. Using the neural network-based representation of the Pauli potential, we calculated the ground-state densities of the 1-D model system by solving the Euler equation. We calculated the Pauli kinetic energy using the neural network-based Pauli potential employing the exact relation between the Pauli kinetic energy functional and the potential. The sum of the neural network-based Pauli kinetic energy and the von Weizsäcker kinetic energy resulted in an accurate estimation of the total kinetic energy. The approach presented in this paper can be employed for the calculation of Pauli potential and Pauli kinetic energy, obviating the need for a functional derivative. The present study is an important step in the advancement of application of machine learning-based techniques toward the orbital-free density functional theory-based methods.

An efficient zero-order evolutionary method for solving the orbital-free density functional theory problem by direct minimization.

A differential evolution (DE) global optimization method for all-electron orbital-free density functional theory (OF-DFT) is presented. This optimization method does not need information about function derivatives to find extreme solutions. Results for a series of known orbital-free energy functionals are presented. Ground state energies of atoms (H to Ar) are obtained by direct minimization of the energy functional without using either Lagrange multipliers or damping procedures for reaching convergence. Our results are in agreement with previous OF-DFT calculations obtained using the standard Newton-Raphson and trust region methods. Being a zero-order method, the DE method can be applied to optimization problems dealing with non-differentiable functionals or functionals with non-closed forms.

Orbital-free density functional theory: Differences and similarities between electronic and nuclear systems

Abstract Orbital-free density functional theory (OF-DFT) has been used when studying atoms, molecules, and solids. In nuclear physics, applications of OF-DFT have been quite scarce so far, as DFT has been widely applied to the study of many nuclear properties mostly within the Kohn–Sham (KS) scheme. There are many realizations of nuclear KS-DFT, but computations become very demanding for heavy systems, such as superheavy nuclei and the inner crust of neutron stars, and it is hard to describe exotic nuclear shapes using a finite basis made with a limited number of orbitals. These bottlenecks could, in principle, be overcome by an orbital-free formulation of DFT. This work is a first step towards the realistic application of OF-DFT to nuclei. In particular, we have implemented possible choices for an orbital-free kinetic energy and solved the associated Schrödinger equation either with simple potentials or with simplified nuclear density functionals. While the former choice sheds light on the differences between electronic and nuclear systems, the latter choice allows us to discuss the practical applications to nuclei as well as open questions.

Dragon: A multi-GPU orbital-free density functional theory molecular dynamics simulation package for modeling of warm dense matter

As progress in electronic structure theoretical methods is made, ab initio molecular dynamics (MD) based on orbital-free density functional theory (OF-DFT) is becoming increasingly more successful at substituting the traditional, very accurate but computationally costly Kohn-Sham (KS) approach for simulations of matter at the challenging warm dense matter (WDM) regime. However, despite the significant cost alleviation of eliminating the dependence on the KS orbitals, OF-DFT MD runs require ∼102 to 103 CPU cores running for days, or even weeks, for simulations of systems comprised of 102 to 103 atoms, depending on thermodynamic conditions. We present Dragon, a multi-GPU OF-DFT MD code for fast and efficient simulations of WDM. With a relatively small allocation of resources (4 to 8 GPU devices) it can provide an order of magnitude speedup for simulations containing O(104) atoms and target systems composed of O(105) atoms at conditions within the WDM regime, which is currently outside the capabilities of CPU codes.

Vortex-lattice structures in rotating Bose-Fermi superfluid mixtures

The system of Bose-Fermi superfluid mixture offers a playground to explore rich macroscopic quantum phenomena. In a recent experiment of Yao et al. (2016) [35], 41K-6Li superfluid mixture is implemented. Coupled quantized vortices are generated via rotating the superfluid mixture, and a few unconventional behaviors on the formations of vortex numbers are observed, which can be traced to boson-fermion interactions. Here we provide a theoretical insight into the unconventional behaviors observed in the experiment. To this end, the orbital-free density functional theory is hired, and its utility is validated by making comparison of the numerical results and a microscopic theory for vortex lattices in strongly interacting Fermi superfluids alone. We also predict interesting phenomena which can be readily explored experimentally, including the novel structures of bosonic-fermionic vortices, and attractive interactions between the bosonic and fermionic vortices.

A modified cusp condition for the single density equations of DFT and orbital-free DFT for atoms

A modified cusp condition for the single density equations of DFT and orbital-free DFT for atoms

Assessing the source of error in the Thomas-Fermi-von Weizsäcker density functional.

We investigate the source of error in the Thomas-Fermi-von Weizsäcker (TFW) density functional relative to Kohn-Sham density functional theory (DFT). In particular, through numerical studies on a range of materials, for a variety of crystal structures subject to strain and atomic displacements, we find that while the ground state electron density in TFW orbital-free DFT is close to the Kohn-Sham density, the corresponding energy deviates significantly from the Kohn-Sham value. We show that these differences are a consequence of the poor representation of the linear response within the TFW approximation for the electronic kinetic energy, confirming conjectures in the literature. In so doing, we find that the energy computed from a non-self-consistent Kohn-Sham calculation using the TFW electronic ground state density is in very good agreement with that obtained from the fully self-consistent Kohn-Sham solution.

Efficient single-grid and multi-grid solvers for real-space orbital-free density functional theory

To improve the computational efficiencies of the real-space orbital-free density functional theory, this work develops a new single-grid solver by directly providing the closed-form solution to the inner iteration and using an improved bisection method to accelerate the line search process in the outer iteration, and extended the single-grid solver to a multi-grid solver. Numerical examples show that while maintaining high level of accuracy, the proposed single-grid solver can improve the computational efficiencies by two orders of magnitude comparing with the methods in the literature (For example, the ground state energy of a 45×45×45 FCC Aluminum cluster (376786 atoms) can be computed within 4h27min on a desktop computer.) and the multi-grid solver can improve the computational efficiencies even once for the cases where high-resolution electron densities are needed.

Automatic differentiation for orbital-free density functional theory.

Differentiable programming has facilitated numerous methodological advances in scientific computing. Physics engines supporting automatic differentiation have simpler code, accelerating the development process and reducing the maintenance burden. Furthermore, fully differentiable simulation tools enable direct evaluation of challenging derivatives-including those directly related to properties measurable by experiment-that are conventionally computed with finite difference methods. Here, we investigate automatic differentiation in the context of orbital-free density functional theory (OFDFT) simulations of materials, introducing PROFESS-AD. Its automatic evaluation of properties derived from first derivatives, including functional potentials, forces, and stresses, facilitates the development and testing of new density functionals, while its direct evaluation of properties requiring higher-order derivatives, such as bulk moduli, elastic constants, and force constants, offers more concise implementations than conventional finite difference methods. For these reasons, PROFESS-AD serves as an excellent prototyping tool and provides new opportunities for OFDFT.

On the origins of spontaneous spherical symmetry-breaking in open-shell atoms through polymer self-consistent field theory.

An alternative approach to density functional theory based on self-consistent field theory for ring polymers is applied to neutral atoms hydrogen to neon in their ground-states. The spontaneous emergence of an atomic shell structure and spherical symmetry-breaking of the total electron density are predicted by the model using the ideas of polymer excluded-volume between pairs of electrons to enforce the Pauli-exclusion principle and an exact electron self-interaction correction. The Pauli potential is approximated by neglecting inter-atomic correlations along with other types of correlations, and comparisons to Hartree-Fock theory are made, which also ignores correlations. The model shows excellent agreement with Hartree-Fock theory to within the standards of orbital-free density functional theory for the atomic binding energies and density profiles of the first six elements, providing exact matches for the elements hydrogen and helium. The predicted shell structure starts to deviate significantly past the element neon, and spherical symmetry-breaking is first predicted to occur at carbon instead of boron. The self-consistent field theory energy functional that describes the model is decomposed into thermodynamic components to trace the origin of spherical symmetry-breaking. It is found to arise from the electron density approaching closer to the nucleus in non-spherical distributions, which lowers the energy despite resulting in frustration between the quantum kinetic energy, electron-electron interaction, and the Pauli exclusion interaction. The symmetry-breaking effect is found to have a minimal impact on the binding energies, which suggests that the spherical-averaging approximation used in previous work is physically reasonable when investigating atomic systems. The pair density contour plots display behavior similar to polymer macro-phase separation, where individual electron pairs occupy single lobe structures that together form a dumbbell shape analogous to the 2p orbital shape. It is further shown that the predicted densities satisfy known constraints and produce the same total electronic density profile that is predicted by other formulations of quantum mechanics.

Geometric potential of the exact electron factorization: Meaning, significance, and application

The theoretical and computational description of materials properties is a task of utmost scientific and technological importance. A first-principles description of electron-electron interactions poses an immense challenge that is usually approached by converting the many-electron problem to an effective one-electron problem. There are different ways to obtain an exact one-electron theory for a many-electron system. An emergent method is the exact electron factorization (EEF) -- one of the branches of the exact factorization approach to many-body systems. In the EEF, the Schr\"odinger equation for one electron, in the environment of all other electrons, is formulated. The influence of the environment is reflected in the potential ${v}^{\mathrm{H}}$, which represents the energy of the environment, and in a potential ${v}^{\mathrm{G}}$, which has a geometrical meaning. In this paper, we focus on ${v}^{\mathrm{G}}$ and study its properties in detail. We investigate the geometric origin of ${v}^{\mathrm{G}}$ as a metric measuring the change of the environment, exemplify how translation and scaling of the state of the environment are reflected in ${v}^{\mathrm{G}}$, and explain its shape for homo- and heteronuclear diatomic model systems. Based on the close connection between the EEF and density functional theory, we also use ${v}^{\mathrm{G}}$ to provide an alternative interpretation to the Pauli potential in orbital-free density functional theory.

Development of nonlocal kinetic-energy density functional for the hybrid QM/MM interaction.

Development of the electronic kinetic-energy density functional is a subject of major interest in theoretical physics and chemistry. In this work, the nonlocal kinetic-energy functional is developed in terms of the response function for the molecular system to realize the orbital free density-functional theory (OF-DFT) to be utilized in the hybrid QM/MM (quantum mechanical/molecular mechanical) method. The present approach shows a clear contrast to the previous functionals where the homogeneous electron gas serves as a reference to build the response function. As a benchmark test, we apply the method to a QM water molecule in a dimer system and that embedded in a condensed environment to make comparisons with the results given by the QM/MM calculations employing the Kohn-Sham DFT. It was found that the energetics and the polarization density of the QM solute under the influence of the MM environment can be adequately reproduced with our approach. This work suggests the potential ability of the kinetic-energy functional based on the response functions for the molecular reference systems.

Mass-zero constrained dynamics for simulations based on orbital-free density functional theory.

A new algorithm for efficient and fully time-reversible integration of first-principles molecular dynamics based on orbital-free density functional theory (OFDFT) is presented. The algorithm adapts to this nontrivial case, the recently introduced Mass-Zero (MaZe) constrained dynamics. The formalism ensures that full adiabatic separation is enforced between nuclear and electronic degrees of freedom and, consequently, that the exact Born-Oppenheimer probability for the nuclei is sampled. Numerical integration of the MaZe dynamics combines standard molecular dynamics algorithms, e.g., Verlet or velocity Verlet, with the SHAKE method to impose the minimum conditions on the electronic degrees of freedom as a set of constraints. The developments presented in this work, which include a bespoke adaptation of the standard SHAKE algorithm, ensure that the quasilinear scaling of OFDFT is preserved by the new method for a broad range of kinetic and exchange-correlation functionals, including nonlocal ones. The efficiency and accuracy of the approach are demonstrated via calculations of static and dynamic properties of liquid sodium in the constant energy and constant temperature ensembles.

Development and Applications of the Density-Based Theory of Chemical Reactivity.

Density functional theory, which is well-recognized for its accuracy and efficiency, has become the workhorse for modeling the electronic structure of molecules and extended materials in recent decades. Nevertheless, establishing a density-based conceptual framework to appreciate bonding, stability, function, reactivity, and other physicochemical properties is still an unaccomplished task. In this Perspective, we at first provide an overview of the four pathways currently available in the literature to tackle the matter, including orbital-free density functional theory, conceptual density functional theory, direct use of density-associated quantities, and the information-theoretic approach. Then, we highlight several recent advances of employing these approaches to realize new understandings for chemical concepts such as covalent bonding, noncovalent interactions, cooperation, frustration, homochirality, chirality hierarchy, electrophilicity, nucleophilicity, regioselectivity, and stereoselectivity. Finally, we provide a few possibilities for the future development of this relatively uncharted territory. Opportunities are abundant, and they are all ours for the taking.

Kinetic energy density for open-shell systems: analysis and development of a novel technique

The quest for an approximate yet accurate kinetic energy density (KED) functional is central to the development of orbital-free density functional theory. While a recipe for closed-shell systems has been proposed earlier, we have shown that it cannot be naïvely extended to open-shell atoms. In this present work, we investigated the efficacy of an ad-hoc recipe to compute the kinetic energy densities for open-shell atoms by extending the methodology used for closed-shell systems. We have also analysed the spin-dependent features of Pauli potentials derived from two previously devised enhancement factors. Further, we have proposed an alternate but exact methodology to systematically compute the KED for atoms of arbitrary spin multiplicity.

First-principles equationof state of CHON resin for inertial confinement fusion applications.

A wide-range (0 to 1044.0 g/cm^{3} and 0 to 10^{9}K) equation-of-state (EOS) table for a CH_{1.72}O_{0.37}N_{0.086} quaternary compound has been constructed based on density-functional theory (DFT) molecular-dynamics (MD) calculations using a combination of Kohn-Sham DFT MD, orbital-free DFT MD, and numerical extrapolation. The first-principles EOS data are compared with predictions of simple models, including the fully ionized ideal gas and the Fermi-degenerate electron gas models, to chart their temperature-density conditions of applicability. The shock Hugoniot, thermodynamic properties, and bulk sound velocities are predicted based on the EOS table and compared to those of C-H compounds. The Hugoniot results show the maximum compression ratio of the C-H-O-N resin is larger than that of CH polystyrene due to the existence of oxygen and nitrogen; while the other properties are similar between CHON and CH. Radiation hydrodynamic simulations have been performed using the table for inertial confinement fusion targets with a CHON ablator and compared with a similar design with CH. The simulations show CHON outperforms CH as the ablator for laser-direct-drive target designs.

Accelerating equilibration in first-principles molecular dynamics with orbital-free density functional theory

We introduce a practical hybrid approach that combines orbital-free density functional theory (DFT) with Kohn-Sham DFT for speeding up first-principles molecular dynamics simulations. Equilibrated ionic configurations are generated using orbital-free DFT for subsequent Kohn-Sham DFT molecular dynamics. This leads to a massive reduction of the simulation time without any sacrifice in accuracy. We assess this finding across systems of different sizes and temperature, up to the warm dense matter regime. To that end, we use the cosine distance between the time series of radial distribution functions representing the ionic configurations. Likewise, we show that the equilibrated ionic configurations from this hybrid approach significantly enhance the accuracy of machine-learning models that replace Kohn-Sham DFT. Our hybrid scheme enables systematic first-principles simulations of warm dense matter that are otherwise hampered by the large numbers of atoms and the prevalent high temperatures. Moreover, our finding provides an additional motivation for developing kinetic and noninteracting free-energy functionals for orbital-free DFT.

Orbital-free spherical density functional theory

Recently, it has been proved that a set of spherically symmetric densities uniquely determines the Coulomb external potential. Utilizing this theorem, a previous method is extended to construct an orbital-free density functional theory. Auxilary spherical non-interacting systems and Kohn–Sham-like equations are constructed. A set of spherical generating functions are defined. These functions have two extra variables and are reduced to the set of spherically symmetric densities if the extra variables are equal to zero. Euler equations are rewritten as Schrödinger-like equations for the square root of the spherically symmetric generating functions. Generalized Pauli potentials are defined, and it is shown that these potentials can be calculated in the knowledge of the generating functions. The Euler equations can be solved with these Pauli potentials.

Efficient time-dependent orbital-free density functional theory: Semilocal adiabatic response

Orbital-free density functional theory and its time-dependent extension are efficient ab initio methods for calculating the electronic structure and dynamics of large systems. Through the calculation of the optical spectra of selected clusters, we reach three important conclusions: (1) The quality of the spectra is strongly affected by the quality of the corresponding ground-state electron density; (2) the adiabatic part of the electronic response to external perturbations can be safely evaluated at the semilocal level; and (3) the nonadiabatic, current-dependent part of the time-dependent Pauli potential is key to recover correct spectral envelopes.

Density Embedding Method for Nanoscale Molecule-Metal Interfaces.

In this work, we extend the applicability of standard Kohn-Sham DFT (KS-DFT) to model realistically sized molecule-metal interfaces where the metal slabs venture into the tens of nanometers in size. Employing state-of-the-art noninteracting kinetic energy functionals, we describe metallic subsystems with orbital-free DFT and combine their electronic structure with molecular subsystems computed at the KS-DFT level resulting in a multiscale subsystem DFT method. The method reproduces within a few millielectronvolts the binding energy difference of water and carbon dioxide molecules adsorbed on the top and hollow sites of an Al(111) surface compared to KS-DFT of the combined supersystem. It is also robust for Born-Oppenheimer molecular dynamics simulations. Very large system sizes are approached with standard computing resources thanks to a parallelization scheme that avoids accumulation of memory at the gather-scatter stage. The results as presented are encouraging and open the door to ab initio simulations of realistically sized, mesoscopic molecule-metal interfaces.