Comparison of density functional approximations and the finite-temperature Hartree-Fock approximation in warm dense lithium

Mixed Resolution-of-the-Identity Compressed Exchange for Hybrid Mixed Deterministic-Stochastic Density Functional Theory from Low to Extreme Temperatures

Exact exchange contributions included in density functional theory calculations have rendered excellent electronic structure results on both cold and extremely hot matter. In this work, we develop a mixed deterministic-stochastic resolution-of-the-identity compressed exchange (mRICE) method for efficient calculation of exact and hybrid electron exchange, suitable for applications alongside mixed stochastic-deterministic density functional theory. mRICE offers accurate calculations of the electronic structure at a largely reduced computation time compared to other compression algorithms, such as Lin's adaptive compressed exchange, for the warm dense matter. mRICE grants flexibility in the number of exchange compression vectors used to resolve the approximated exchange operator kernel, reducing the computation time of the application of the exchange operator to the vectors by up to 40% while maintaining accuracy in electronic structure predictions. We demonstrate mRICE by computing the density of states of warm dense carbon and neon between temperatures of 10 and 50 eV (116,045 and 580,226 K) and comparing timing and accuracy at varying levels of compression. Finally, we carry out mRICE on the difference between the Fock exchange operator and the semilocal exchange potential kernels and show an enhanced convergence of electronic structure calculations at reduced stochastic sampling.

Importance of finite-temperature exchange correlation for warm dense matter calculations

The effects of an explicit temperature dependence in the exchange correlation (XC) free-energy functional upon calculated properties of matter in the warm dense regime are investigated. The comparison is between the Karasiev-Sjostrom-Dufty-Trickey (KSDT) finite-temperature local-density approximation (TLDA) XC functional [Karasiev et al., Phys. Rev. Lett. 112, 076403 (2014)PRLTAO0031-900710.1103/PhysRevLett.112.076403] parametrized from restricted path-integral Monte Carlo data on the homogeneous electron gas (HEG) and the conventional Monte Carlo parametrization ground-state LDA XC [Perdew-Zunger (PZ)] functional evaluated with T-dependent densities. Both Kohn-Sham (KS) and orbital-free density-functional theories are used, depending upon computational resource demands. Compared to the PZ functional, the KSDT functional generally lowers the dc electrical conductivity of low-density Al, yielding improved agreement with experiment. The greatest lowering is about 15% for T=15 kK. Correspondingly, the KS band structure of low-density fcc Al from the KSDT functional exhibits a clear increase in interband separation above the Fermi level compared to the PZ bands. In some density-temperature regimes, the deuterium equations of state obtained from the two XC functionals exhibit pressure differences as large as 4% and a 6% range of differences. However, the hydrogen principal Hugoniot is insensitive to the explicit XC T dependence because of cancellation between the energy and pressure-volume work difference terms in the Rankine-Hugoniot equation. Finally, the temperature at which the HEG becomes unstable is T≥7200 K for the T-dependent XC, a result that the ground-state XC underestimates by about 1000 K.

Orbital-free density functional theory implementation with the projector augmented-wave method

We present a computational scheme for orbital-free density functional theory (OFDFT) that simultaneously provides access to all-electron values and preserves the OFDFT linear scaling as a function of the system size. Using the projector augmented-wave method (PAW) in combination with real-space methods, we overcome some obstacles faced by other available implementation schemes. Specifically, the advantages of using the PAW method are twofold. First, PAW reproduces all-electron values offering freedom in adjusting the convergence parameters and the atomic setups allow tuning the numerical accuracy per element. Second, PAW can provide a solution to some of the convergence problems exhibited in other OFDFT implementations based on Kohn-Sham (KS) codes. Using PAW and real-space methods, our orbital-free results agree with the reference all-electron values with a mean absolute error of 10 meV and the number of iterations required by the self-consistent cycle is comparable to the KS method. The comparison of all-electron and pseudopotential bulk modulus and lattice constant reveal an enormous difference, demonstrating that in order to assess the performance of OFDFT functionals it is necessary to use implementations that obtain all-electron values. The proposed combination of methods is the most promising route currently available. We finally show that a parametrized kinetic energy functional can give lattice constants and bulk moduli comparable in accuracy to those obtained by the KS PBE method, exemplified with the case of diamond.

Accurate homogeneous electron gas exchange-correlation free energy for local spin-density calculations

An accurate analytical parametrization for the exchange-correlation free energy of the homogeneous electron gas, including interpolation for partial spin polarization, is derived via thermodynamic analysis of recent restricted path integral Monte Carlo (RPIMC) data. This parametrization constitutes the local spin density approximation (LSDA) for the exchange-correlation functional in density functional theory. The new finite-temperature LSDA reproduces the RPIMC data well, satisfies the correct high-density and low- and high-T asymptotic limits, and is well behaved beyond the range of the RPIMC data, suggestive of broad utility.