Efficient real space formalism for hybrid density functionals.

We present an efficient real space formalism for hybrid exchange-correlation functionals in generalized Kohn-Sham density functional theory (DFT). In particular, we develop an efficient representation for any function of the real space finite-difference Laplacian matrix by leveraging its Kronecker product structure, thereby enabling the time to solution of associated linear systems to be highly competitive with the fast Fourier transform scheme while not imposing any restrictions on the boundary conditions. We implement this formalism for both the unscreened and range-separated variants of hybrid functionals.

Overcoming the Chemical Complexity Bottleneck in on-the-Fly Machine Learned Molecular Dynamics Simulations.

We develop a framework for on-the-fly machine learned force field molecular dynamics simulations based on the multipole featurization scheme that overcomes the bottleneck with the number of chemical elements. Considering bulk systems with up to 6 elements, we demonstrate that the number of density functional theory calls remains approximately independent of the number of chemical elements, in contrast to the increase in the smooth overlap of the atomic positions scheme.

Self-Consistent Convolutional Density Functional Approximations: Application to Adsorption at Metal Surfaces.

The exchange-correlation (XC) functional in density functional theory is used to approximate multi-electron interactions. A plethora of different functionals are available, but nearly all are based on the hierarchy of inputs commonly referred to as "Jacob's ladder." This paper introduces an approach to construct XC functionals with inputs from convolutions of arbitrary kernels with the electron density, providing a route to move beyond Jacob's ladder.

Strain engineering of Zeeman and Rashba effects in transition metal dichalcogenide nanotubes and their Janus variants: anstudy.

We study the influence of mechanical deformations on the Zeeman and Rashba effects in transition metal dichalcogenide nanotubes and their Janus variants from first principles. In particular, we perform symmetry-adapted density functional theory simulations with spin-orbit coupling to determine the variation in the electronic band structure splittings with axial and torsional deformations. We find significant effects in molybdenum and tungsten nanotubes, for which the Zeeman splitting decreases with increase in strain, going to zero for large enough tensile/shear strains, while the Rashba splitting coefficient increases linearly with shear strain, while being zero for all tensile strains, a consequence of the inversion symmetry remaining unbroken.

Kohn-Sham accuracy from orbital-free density functional theory via Δ-machine learning.

We present a Δ-machine learning model for obtaining Kohn-Sham accuracy from orbital-free density functional theory (DFT) calculations. In particular, we employ a machine-learned force field (MLFF) scheme based on the kernel method to capture the difference between Kohn-Sham and orbital-free DFT energies/forces. We implement this model in the context of on-the-fly molecular dynamics simulations and study its accuracy, performance, and sensitivity to parameters for representative systems.

Calculation of phonons in real-space density functional theory.

We present an accurate and efficient formulation for the calculation of phonons in real-space Kohn-Sham density functional theory. Specifically, employing a local exchange-correlation functional, norm-conserving pseudopotential in the Kleinman-Bylander representation, and local form for the electrostatics, we derive expressions for the dynamical matrix and associated Sternheimer equation that are particularly amenable to the real-space finite-difference method, within the framework of density functional perturbation theory. In particular, the formulation is applicable to insulating as well as metallic systems of any dimensionality, enabling the efficient and accurate treatment of semi-infinite and bulk systems alike, for both orthogonal and nonorthogonal cells.

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.

GPU acceleration of local and semilocal density functional calculations in the SPARC electronic structure code.

We present a Graphics Processing Unit (GPU)-accelerated version of the real-space SPARC electronic structure code for performing Kohn-Sham density functional theory calculations within the local density and generalized gradient approximations. In particular, we develop a modular math-kernel based implementation for NVIDIA architectures wherein the computationally expensive operations are carried out on the GPUs, with the remainder of the workload retained on the central processing units (CPUs). Using representative bulk and slab examples, we show that relative to CPU-only execution, GPUs enable speedups of up to 6× and 60× in node and core hours, respectively, bringing time to solution down to less than 30 s for a metallic system with over 14 000 electrons and enabling significant reductions in computational resources required for a given wall time.

Properties of carbon up to 10 million kelvin from Kohn-Sham density functional theory molecular dynamics.

Accurately modeling dense plasmas over wide-ranging conditions of pressure and temperature is a grand challenge critically important to our understanding of stellar and planetary physics as well as inertial confinement fusion. In this work, we employ Kohn-Sham density functional theory (DFT) molecular dynamics (MD) to compute the properties of carbon at warm and hot dense matter conditions in the vicinity of the principal Hugoniot. In particular, we calculate the equation of state (EOS), Hugoniot, pair distribution functions, and diffusion coefficients for carbon at densities spanning 8 g/cm^{3} to 16 g/cm^{3} and temperatures ranging from 100 kK to 10 MK using the Spectral Quadrature method.

On the bending of rectangular atomic monolayers along different directions: anstudy.

We study the bending of rectangular atomic monolayers along different directions from first principles. Specifically, choosing the phosphorene, GeS, TiS, and AsSmonolayers as representative examples, we perform Kohn-Sham density functional theory calculations to determine the variation in transverse flexoelectric coefficient and bending modulus with the direction of bending. We find that while the flexoelectric coefficient is nearly isotropic, there is significant and complex anisotropy in bending modulus that also differs between the monolayers, with extremal values not necessarily occurring along the principal directions.

Pseudodiagonalization Method for Accelerating Nonlinear Subspace Diagonalization in Density Functional Theory.

In density functional theory, each self-consistent field (SCF) nonlinear step updates the discretized Kohn-Sham orbitals by solving a linear eigenvalue problem. The concept of pseudodiagonalization is to solve this linear eigenvalue problem approximately and specifically utilizing a method involving a small number of Jacobi rotations that takes advantage of the good initial guess to the solution given by the approximation to the orbitals from the previous SCF iteration. The approximate solution to the linear eigenvalue problem can be very rapid, particularly for those steps near SCF convergence.

Real-space density kernel method for Kohn-Sham density functional theory calculations at high temperature.

Kohn-Sham density functional theory calculations using conventional diagonalization based methods become increasingly expensive as temperature increases due to the need to compute increasing numbers of partially occupied states. We present a density matrix based method for Kohn-Sham calculations at high temperatures that eliminates the need for diagonalization entirely, thus reducing the cost of such calculations significantly. Specifically, we develop real-space expressions for the electron density, electronic free energy, Hellmann-Feynman forces, and Hellmann-Feynman stress tensor in terms of an orthonormal auxiliary orbital basis and its density kernel transform, the density kernel being the matrix representation of the density operator in the auxiliary basis.

Accurate parameterization of the kinetic energy functional for calculations using exact-exchange.

Electronic structure calculations based on Kohn-Sham density functional theory (KSDFT) that incorporate exact-exchange or hybrid functionals are associated with a large computational expense, a consequence of the inherent cubic scaling bottleneck and large associated prefactor, which limits the length and time scales that can be accessed. Although orbital-free density functional theory (OFDFT) calculations scale linearly with system size and are associated with a significantly smaller prefactor, they are limited by the absence of accurate density-dependent kinetic energy functionals. Therefore, the development of accurate density-dependent kinetic energy functionals is important for OFDFT calculations of large realistic systems.

Torsional strain engineering of transition metal dichalcogenide nanotubes: anstudy.

We study the effect of torsional deformations on the electronic properties of single-walled transition metal dichalcogenide (TMD) nanotubes. In particular, considering forty-five select armchair and zigzag TMD nanotubes, we perform symmetry-adapted Kohn-Sham density functional theory calculations to determine the variation in bandgap and effective mass of charge carriers with twist. We find that metallic nanotubes remain so even after deformation, whereas semiconducting nanotubes experience a decrease in bandgap with twist-originally direct bandgaps become indirect-resulting in semiconductor to metal transitions.

Torsional moduli of transition metal dichalcogenide nanotubes from first principles.

We calculate the torsional moduli of single-walled transition metal dichalcogenide (TMD) nanotubes usingdensity functional theory (DFT). Specifically, considering forty-five select TMD nanotubes, we perform symmetry-adapted DFT calculations to calculate the torsional moduli for the armchair and zigzag variants of these materials in the low-twist regime and at practically relevant diameters. We find that the torsional moduli follow the trend: MS> MSe> MTe.

Implementation of Perdew-Zunger self-interaction correction in real space using Fermi-Löwdin orbitals.

Most widely used density functional approximations suffer from self-interaction error, which can be corrected using the Perdew-Zunger (PZ) self-interaction correction (SIC). We implement the recently proposed size-extensive formulation of PZ-SIC using Fermi-Löwdin Orbitals (FLOs) in real space, which is amenable to systematic convergence and large-scale parallelization. We verify the new formulation within the generalized Slater scheme by computing atomization energies and ionization potentials of selected molecules and comparing to those obtained by existing FLOSIC implementations in Gaussian based codes.

Development of a Multiphase Beryllium Equation of State and Physics-based Variations.

We construct a family of beryllium (Be) multiphase equation of state (EOS) models that consists of a baseline ("optimal") EOS and variations on the baseline to account for physics-based uncertainties. The Be baseline EOS is constructed to reproduce a set of self-consistent data and theory including known phase boundaries, the principal Hugoniot, isobars, and isotherms from diamond-anvil cell experiments. Three phases are considered, including the known hexagonal closed-packed (hcp) phase, the liquid, and the theoretically predicted high-pressure body-centered cubic (bcc) phase.

Flexoelectricity in atomic monolayers from first principles.

We study the flexoelectric effect in fifty-four select atomic monolayers using ab initio Density Functional Theory (DFT). Specifically, considering representative materials from each of the Group III monochalcogenides, transition metal dichalcogenides (TMDs), Groups IV, III-V, and V monolayers, Group IV dichalcogenides, Group IV monochalcogenides, transition metal trichalcogenides (TMTs), and Group V chalcogenides, we perform symmetry-adapted DFT simulations to calculate transversal flexoelectric coefficients along the principal directions at practically relevant bending curvatures. We find that the materials demonstrate linear behavior and have similar coefficients along both principal directions, with values for TMTs being up to a factor of five larger than those of graphene.

Real-space formulation of the stress tensor for O(N) density functional theory: Application to high temperature calculations.

We present an accurate and efficient real-space formulation of the Hellmann-Feynman stress tensor for O(N) Kohn-Sham density functional theory (DFT). While applicable at any temperature, the formulation is most efficient at high temperature where the Fermi-Dirac distribution becomes smoother and the density matrix becomes correspondingly more localized. We first rewrite the orbital-dependent stress tensor for real-space DFT in terms of the density matrix, thereby making it amenable to O(N) methods.

Bending moduli for forty-four select atomic monolayers from first principles.

We calculate bending moduli along the principal directions for forty-four select atomic monolayers using ab initio density functional theory (DFT). Specifically, considering representative materials from each of Groups IV, III-V, V monolayers, Group IV monochalcogenides, transition metal trichalcogenides, transition metal dichalcogenides and Group III monochalcogenides, we utilize the recently developed Cyclic DFT method to calculate the bending moduli in the practically relevant but previously intractable low-curvature limit. We find that the moduli generally increase with thickness of the monolayer, while spanning three orders of magnitude between the different materials.

On the calculation of the stress tensor in real-space Kohn-Sham density functional theory.

We present an accurate and efficient formulation of the stress tensor for real-space Kohn-Sham density functional theory calculations. Specifically, while employing a local formulation of the electrostatics, we derive a linear-scaling expression for the stress tensor that is applicable to simulations with unit cells of arbitrary symmetry, semilocal exchange-correlation functionals, and Brillouin zone integration. In particular, we rewrite the contributions arising from the self-energy and the nonlocal pseudopotential energy to make them amenable to the real-space finite-difference discretization, achieving up to three orders of magnitude improvement in the accuracy of the computed stresses.

Discrete discontinuous basis projection method for large-scale electronic structure calculations.

We present an approach to accelerate real-space electronic structure methods several fold, without loss of accuracy, by reducing the dimension of the discrete eigenproblem that must be solved. To accomplish this, we construct an efficient, systematically improvable, discontinuous basis spanning the occupied subspace and project the real-space Hamiltonian onto the span. In calculations on a range of systems, we find that accurate energies and forces are obtained with 8-25 basis functions per atom, reducing the dimension of the associated real-space eigenproblems by 1-3 orders of magnitude.

Two-Level Chebyshev Filter Based Complementary Subspace Method: Pushing the Envelope of Large-Scale Electronic Structure Calculations.

We describe a novel iterative strategy for Kohn-Sham density functional theory calculations aimed at large systems (>1,000 electrons), applicable to metals and insulators alike. In lieu of explicit diagonalization of the Kohn-Sham Hamiltonian on every self-consistent field (SCF) iteration, we employ a two-level Chebyshev polynomial filter based complementary subspace strategy to (1) compute a set of vectors that span the occupied subspace of the Hamiltonian; (2) reduce subspace diagonalization to just partially occupied states; and (3) obtain those states in an efficient, scalable manner via an inner Chebyshev filter iteration. By reducing the necessary computation to just partially occupied states and obtaining these through an inner Chebyshev iteration, our approach reduces the cost of large metallic calculations significantly, while eliminating subspace diagonalization for insulating systems altogether.