Exchange-correlation Free Energy Research Articles

On the Chemical Potential of Many-Body Perturbation Theory in Extended Systems.

Finite-temperature many-body perturbation theory in the grand-canonical ensemble is fundamental to numerous methods for computing electronic properties at nonzero temperature, such as finite-temperature coupled-cluster. In most applications it is the average number of electrons that is known rather than the chemical potential. Expensive correlation calculations must be repeated iteratively in search for the interacting chemical potential that yields the desired average number of electrons. In extended systems with mobile charges the situation is particular, however. Long-ranged electrostatic forces drive the charges such that the average ratio of negative and positive charges is one for any finite chemical potential. All properties per electron are expected to be virtually independent of the chemical potential, as they are in an electric wire at different voltage potentials. This work shows that per electron, the exchange-correlation free energy and the exchange-correlation grand potential indeed agree in the infinite-size limit. Thus, only one expensive correlation calculation suffices for each system size, sparing the search for the interacting chemical potential. This work also demonstrates the importance of regularizing the Coulomb interaction such that each electron on average interacts only with as many electrons as there are electrons in the simulation, avoiding interactions with periodic images. Numerical calculations of the warm uniform electron gas have been conducted with the Spencer-Alavi regularization employing the finite-temperature Hartree approximation for the self-consistent field and linearized finite-temperature direct-ring coupled-cluster doubles for treating correlation.

Meta-GGA exchange-correlation free energy density functional to increase the accuracy of warm dense matter simulations

We discuss strategies for thermalization of the ground-state meta-generalized gradient approximation (meta-GGA) exchange-correlation (XC) functionals. A simple but accurate scheme is implemented via universal additive thermal correction to XC using a perturbativelike self-consistent approach. The additive correction with explicit temperature dependence is applied to the ground-state deorbitalized, strongly constrained, and appropriately normed (SCAN-L) meta-GGA XC leading to a thermal XC functional denoted here as T-SCAN-L. The thermal T-SCAN-L meta-GGA functional shows significant improvement in density functional theory calculation accuracy for warm dense matter by a factor of 3 to 10, achieving an accuracy of total pressure between a few tenths and $\ensuremath{\sim}1$% when compared to traditional XC functionals, as demonstrated by the comparison to path-integral Monte Carlo simulations for helium equation of state. The T-SCAN-L calculations of dc conductivity of warm dense aluminum also give better agreement with experiments over other XC functionals such as PBE and SCAN-L.

Thermodynamics of the uniform electron gas: Fermionic path integral Monte Carlo simulations in the restricted grand canonical ensemble

The uniform electron gas (UEG) is one of the key models for the understanding of warm dense matter—an exotic, highly compressed state of matter between solid and plasma phases. The difficulty in modelling the UEG arises from the need to simultaneously account for Coulomb correlations, quantum effects, and exchange effects, as well as finite temperature. The most accurate results so far were obtained from quantum Monte Carlo (QMC) simulations with a variety of representations. However, QMC for electrons is hampered by the fermion sign problem. Here, we present results from a novel fermionic propagator path integral Monte Carlo in the restricted grand canonical ensemble. The ab initio simulation results for the spin‐resolved pair distribution functions and static structure factor are reported for two isotherms (T in the units of the Fermi temperature). Furthermore, we combine the results from the linear response theory in the Singwi‐Tosi‐Land‐Sjölander scheme with the QMC data to remove finite‐size errors in the interaction energy. We present a new corrected parametrization for the interaction energy and the exchange–correlation free energy in the thermodynamic limit, and benchmark our results against the restricted path integral Monte Carlo by Brown et al. [Phys. Rev. Lett. 110, 146405 (2013)] and configuration path integral Monte Carlo/permutation‐blocking path integral Monte Carlo by Dornheim et al. [Phys. Rev. Lett. 117, 115701 (2016)].

Strongly coupled electron liquid: Ab initio path integral Monte Carlo simulations and dielectric theories

The strongly coupled electron liquid provides a unique opportunity to study the complex interplay of strong coupling with quantum degeneracy effects and thermal excitations. To this end, we carry out extensive \textit{ab initio} path integral Monte Carlo (PIMC) simulations to compute the static structure factor, interaction energy, density response function, and the corresponding static local field correction in the range of $20\leq r_s \leq 100$ and $0.5\leq \theta\leq 4$. We subsequently compare these data to several dielectric approximations, and find that different schemes are capable to reproduce different features of the PIMC results at certain parameters. Moreover, we provide a comprehensive data table of interaction energies and compare those to two recent parametrizations of the exchange-correlation free energy, where they are available. Finally, we briefly touch upon the possibility of a charge-density wave. The present study is complementary to previous investigations of the uniform electron gas in the warm dense matter regime and, thus, further completes our current picture of this fundamental model system at finite temperature. All PIMC data are available online.

Abinitio Exchange-Correlation Free Energy of the Uniform Electron Gas at Warm Dense Matter Conditions.

In a recent Letter [T. Dornheim etal., Phys. Rev. Lett. 117, 156403 (2016)PRLTAO0031-900710.1103/PhysRevLett.117.156403], we presented the first quantum MonteCarlo (QMC) results for the warm dense electron gas in the thermodynamic limit. However, a complete parametrization of the exchange-correlation free energy with respect to density, temperature, and spin polarization remained out of reach due to the absence of (i)accurate QMC results below θ=k_{B}T/E_{F}=0.5 and (ii)QMC results for spin polarizations different from the paramagnetic case. Here we overcome both remaining limitations. By closing the gap to the ground state and by performing extensive QMC simulations for different spin polarizations, we are able to obtain the first completely abinitio exchange-correlation free energy functional; the accuracy achieved is an unprecedented ∼0.3%. This also allows us to quantify the accuracy and systematic errors of various previous approximate functionals.

Spin-resolved correlations in the warm-dense homogeneous electron gas

We have studied spin-resolved correlations in the warm-dense homogeneous electron gas by determining the linear density and spin-density response functions, within the dynamical self-consistent mean-field theory of Singwi et al. The calculated spin-resolved pair-correlation function g σ σ′(r) is compared with the recent restricted path-integral Monte Carlo (RPIMC) simulations due to Brown et al. [Phys. Rev. Lett. 110, 146405 (2013)], while interaction energy E int and exchange-correlation free energy F xc with the RPIMC and very recent ab initio quantum Monte Carlo (QMC) simulations by Dornheim et al. [Phys. Rev. Lett. 117, 156403 (2016)]. g ↑↓(r) is found to be in good agreement with the RPIMC data, while a mismatch is seen in g ↑↑(r) at small r where it becomes somewhat negative. As an interesting result, it is deduced that a non-monotonic T-dependence of g(0) is driven primarily by g ↑↓(0). Our results of E int and F xc exhibit an excellent agreement with the QMC study due to Dornheim et al., which deals with the finite-size correction quite accurately. We observe, however, a visible deviation of E int from the RPIMC data for high densities (~8% at r s = 1). Further, we have extended our study to the fully spin-polarized phase. Again, with the exception of high density region, we find a good agreement of E int with the RPIMC data. This points to the need of settling the problem of finite-size correction in the spin-polarized phase also. Interestingly, we also find that the thermal effects tend to oppose spatial localization as well as spin polarization of electrons.

AbInitio Quantum Monte Carlo Simulation of the Warm Dense Electron Gas in the Thermodynamic Limit.

We perform abinitio quantum MonteCarlo (QMC) simulations of the warm dense uniform electron gas in the thermodynamic limit. By combining QMC data with the linear response theory, we are able to remove finite-size errors from the potential energy over the substantial parts of the warm dense regime, overcoming the deficiencies of the existing finite-size corrections by Brown etal. [Phys. Rev. Lett. 110, 146405 (2013)]. Extensive new QMC results for up to N=1000 electrons enable us to compute the potential energy V and the exchange-correlation free energy F_{xc} of the macroscopic electron gas with an unprecedented accuracy of |ΔV|/|V|,|ΔF_{xc}|/|F|_{xc}∼10^{-3}. A comparison of our new data to the recent parametrization of F_{xc} by Karasiev etal. [Phys. Rev. Lett. 112, 076403 (2014)] reveals significant deviations to the latter.

Exact thermal density functional theory for a model system: Correlation components and accuracy of the zero-temperature exchange-correlation approximation

Thermal density functional theory calculations often use the Mermin-Kohn-Sham scheme, but employ ground-state approximations to the exchange-correlation (XC) free energy. In the simplest solvable nontrivial model, an asymmetric Hubbard dimer, we calculate the exact many-body energies and the exact Mermin-Kohn-Sham functionals for this system and extract the exact XC free energy. For moderate temperatures and weak correlation, we find this approximation to be excellent. We extract various exact free-energy correlation components and the exact adiabatic connection formula.

Gradient corrections to the exchange-correlation free energy

We develop the first order gradient correction to the exchange-correlation free energy of the homogeneous electron gas for use in finite temperature density functional calculations. Based on this we propose and implement a simple temperature dependent extension for functionals beyond the local density approximation. These finite temperature functionals show improvement over zero temperature functionals as compared to path integral Monte Carlo calculations for deuterium and perform without computational cost increase compared to zero temperature functionals and so should be used for finite temperature calculations.

An improved classical mapping method for homogeneous electron gases at finite temperature.

We introduce a modified classical mapping method to predict the exchange-correlation free energy and the structure of homogeneous electron gases (HEG) at finite temperature. With the classical map temperature parameterized on the basis of the quantum Monte Carlo simulation data for the correlation energy and exact results at high and low temperature limits, the new theoretical procedure greatly improves the classical mapping method for correlating the energetic properties HEG over a broad range of thermodynamic conditions. Improvement can also be identified in predicting the long-range components of the spin-averaged pair correlation functions.

Structure and thermodynamic properties of relativistic electron gases.

Relativistic effect is important in many quantum systems but theoretically complicated from both fundamental and practical perspectives. Herein we introduce an efficient computational procedure to predict the structure and energetic properties of relativistic quantum systems by mapping the Pauli principle into an effective pairwise-additive potential such that the properties of relativistic nonquantum systems can be readily predicted from conventional liquid-state methods. We applied our theoretical procedure to relativistic uniform electron gases and compared the pair correlation functions with those for systems of nonrelativistic electrons. A simple analytical expression has been developed to correlate the exchange-correlation free energy of relativistic uniform electron systems.

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 MonteCarlo (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.

Uniform electron gas at finite temperatures

We calculate the free energy of the quantum uniform electron gas for temperatures from near 0 to 100 times the Fermi energy, approaching the classical limit. An extension of the Vashista-Singwi theory to finite temperatures and a self-consistent compressibility sum rule is presented. Comparisons are made to other local-field correction methods, as well as recent quantum Monte Carlo simulation and classical map-based results. Accurate fits to the exchange-correlation free energy from both theory and simulation are given for future practical applications.

Electron effective mass in the strongly correlated two-dimensional uniform electron fluid from finite-temperature calculations

The very-low temperature thermal effective mass m* of paramagnetic and ferromagnetic electrons in a uniform electron fluid in two dimensions is studied. Analytical and numerical evaluations are used to meaningfully define an m*, even in the the Hartree-Fock approximation. The Hartree-Fock m* decreases linearly with the electron-disk radius r_s. Correlation effects lead to strong cancellations between exchange and correlation. Thus the effective mass is enhanced with increasing r_s for the unpolarized fluid, while m* decreases with the r_s of the polarized fluid. The effective mass is calculated from the coefficient of the quadratic temperature dependence of exchange-correlation free energy F_{xc}. This is calculated in a physically transparent manner using a new formula for the effective mass. This uses the T=0 pair-distribution functions of Gori-Giorgi et al., and the temperature derivative of a quantum analogue of the potential of mean-force well known in the statistical mechanics of classical fluids. The results are compared with recent quantum Monte-Carlo simulations at T=0, as well as with other available experimental and theoretical data for the effective mass.

Variational-average-atom-in-quantum-plasmas (VAAQP) code and virial theorem: Equation-of-state and shock-Hugoniot calculations for warm dense Al, Fe, Cu, and Pb

The numerical code VAAQP (variational average atom in quantum plasmas), which is based on a fully variational model of equilibrium dense plasmas, is applied to equation-of-state calculations for aluminum, iron, copper, and lead in the warm-dense-matter regime. VAAQP does not impose the neutrality of the Wigner-Seitz ion sphere; it provides the average-atom structure and the mean ionization self-consistently from the solution of the variational equations. The formula used for the electronic pressure is simple and does not require any numerical differentiation. In this paper, the virial theorem is derived in both nonrelativistic and relativistic versions of the model. This theorem allows one to express the electron pressure as a combination of the electron kinetic and interaction energies. It is shown that the model fulfills automatically the virial theorem in the case of local-density approximations to the exchange-correlation free-energy. Applications of the model to the equation-of-state and Hugoniot shock adiabat of aluminum, iron, copper, and lead in the warm-dense-matter regime are presented. Comparisons with other approaches, including the inferno model, and with available experimental data are given. This work allows one to understand the thermodynamic consistency issues in the existing average-atom models. Starting from the case of aluminum, a comparative study of the thermodynamic consistency of the models is proposed. A preliminary study of the validity domain of the inferno model is also included.

Spin-polarized electron liquid at arbitrary temperatures: Exchange-correlation energies, electron-distribution functions, and the static response functions

We use a recently introduced classical mapping of the Coulomb interactions in a quantum electron liquid [Phys. Rev. Lett. 84, 959 (2000)] to present a unified treatment of the thermodynamic properties and the static response of the finite-temperature electron liquid, valid for arbitrary coupling and spin-polarization. The method is based on using a ``quantum temperature'' ${T}_{q}$ such that the distribution functions of the classical electron liquid at ${T}_{q}$ leads to the same correlation energy as the quantum electron liquid at $T=0.$ The functional form of ${T}_{q}{(r}_{s})$ is presented. The electron-electron pair-distribution functions (PDF's) calculated using ${T}_{q}$ are in good quantitative agreement with available $(T=0)$ quantum Monte Carlo results. The method provides a means of treating strong-coupling regimes of $n,T,$ and $\ensuremath{\zeta}$ currently unexplored by quantum Monte Carlo or Feenberg-functional methods. The exchange-correlation free energies, distribution functions ${g}_{11}(r),$ ${g}_{12}(r),$ ${g}_{22}(r)$ and the local-field corrections to the static response functions as a function of density n, temperature T, and spin polarization $\ensuremath{\zeta}$ are presented and compared with any available finite-$T$ results. The exchange-correlation free energy ${f}_{\mathrm{xc}}(n,T,\ensuremath{\zeta}),$ is given in a parametrized form. It satisfies the expected analytic behavior in various limits of temperature, density, and spin polarization, and can be used for calculating other properties like the equation of state, the exchange-correlation potentials, compressibility, etc. The static local-field correction provides a static response function which is consistent with the PDF's and the relevant sum rules. Finally, we use the finite-$T$ xc-potentials to examine the Kohn-Sham bound- and continuum states at an ${\mathrm{Al}}^{13+}$ nucleus immersed in a hot electron gas to show the significance of the xc-potentials.