Forces from Stochastic Density Functional Theory under Nonorthogonal Atom-Centered Basis Sets

Time-dependent density-functional-theory calculations of the nonlocal electron stopping range for inertial confinement fusion applications

Nonlocal electron transport is important for understanding laser-target coupling for laser-direct-drive (LDD) inertial confinement fusion (ICF) simulations. Current models for the nonlocal electron mean free path in radiation-hydrodynamic codes are based on plasma-physics models developed decades ago; improvements are needed to accurately predict the electron conduction in LDD simulations of ICF target implosions. We utilized time-dependent density functional theory (TD-DFT) to calculate the electron stopping power (SP) in the so-called conduction-zone plasmas of polystyrene in a wide range of densities and temperatures relevant to LDD. Compared with the modified Lee-More model, the TD-DFT calculations indicated a lower SP and a higher stopping range for nonlocal electrons. We fit these electron SP calculations to obtain a global analytical model for the electron stopping range as a function of plasma conditions and the nonlocal electron kinetic energy. This model was implemented in the one-dimensional radiation-hydrodynamic code lilac to perform simulations of LDD ICF implosions, which are further compared with simulations by the standard modified Lee-More model. Results from these integrated simulations are discussed in terms of the implications of this TD-DFT-based mean-free-path model to ICF simulations.

Stochastic and mixed density functional theory within the projector augmented wave formalism for simulation of warm dense matter

Stochastic density functional theory (DFT) and mixed stochastic-deterministic DFT are burgeoning approaches for the calculation of the equation of state and transport properties in materials under extreme conditions. In the intermediate warm dense matter regime, a state between correlated condensed matter and kinetic plasma, electrons can range from being highly localized around nuclei to delocalized over the whole simulation cell. The plane-wave basis pseudopotential approach is thus the typical tool of choice for modeling such systems at the DFT level. Unfortunately, stochastic DFT methods scale as the square of the maximum plane-wave energy in this basis. To reduce the effect of this scaling and improve the overall description of the electrons within the pseudopotential approximation, we present stochastic and mixed DFT approaches developed and implemented within the projector augmented wave formalism. We compare results between the different DFT approaches for both single-point and molecular dynamics trajectories and present calculations of self-diffusion coefficients of solid density carbon from 1 to 50 eV.

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.

Stochastic Vector Techniques in Ground-State Electronic Structure

We review a suite of stochastic vector computational approaches for studying the electronic structure of extended condensed matter systems. These techniques help reduce algorithmic complexity, facilitate efficient parallelization, simplify computational tasks, accelerate calculations, and diminish memory requirements. While their scope is vast, we limit our study to ground-state and finite temperature density functional theory (DFT) and second-order perturbation theory. More advanced topics, such as quasiparticle (charge) and optical (neutral) excitations and higher-order processes, are covered elsewhere. We start by explaining how to use stochastic vectors in computations, characterizing the associated statistical errors. Next, we show how to estimate the electron density in DFT and discuss highly effective techniques to reduce statistical errors. Finally, we review the use of stochastic vector techniques for calculating correlation energies within the second-order Møller-Plesset perturbation theory and its finite temperature variational form. Example calculation results are presented and used to demonstrate the efficacy of the methods. Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.