Framework for Laplacian-Level Noninteracting Free-Energy Density Functionals.

A framework for orbital-free Laplacian-level meta-generalized-gradient approximation (meta-GGA) for the noninteracting free-energy-density functionals based upon analysis of the fourth-order gradient expansion is developed. The framework presented here provides a new tool for developing advanced orbital-free thermal functionals at the meta-GGA level of theory. A nonempirical meta-GGA functional, which in the slowly varying density limit correctly reduces to the finite-temperature fourth-order gradient expansion for the noninteracting free energy, is constructed. Application to warm dense helium demonstrates that the developed meta-GGA functional drastically increases the accuracy of orbital-free-density functional theory simulations at temperatures below 40 eV as compared to the lower Thomas-Fermi and GGA rungs.

Light-enhanced oxygen degradation of MAPbBr3 single crystal.

Organometal halide perovskites are promising materials for optoelectronic applications, whose commercial realization depends critically on their stability under multiple environmental factors. In this study, a methylammonium lead bromide (MAPbBr3) single crystal was cleaved and exposed to simultaneous oxygen and light illumination under ultrahigh vacuum (UHV). The exposure process was monitored using X-ray photoelectron spectroscopy (XPS) with precise control of the exposure time and oxygen pressure. It was found that the combination of oxygen and light accelerated the degradation of MAPbBr3, which could not be viewed as a simple addition of that by oxygen-only and light-only exposures. The XPS spectra showed significant loss of carbon, bromine, and nitrogen at an oxygen exposure of 1010 Langmuir with light illumination, approximately 17 times of the additive effects of oxygen-only and light-only exposures. It was also found that the photoluminescence (PL) emission was much weakened by oxygen and light co-exposure, while previous reports had shown that PL was substantially enhanced by oxygen-only exposure. Measurements using a scanning electron microscope (SEM) and focused ion beam (FIB) demonstrated that the crystal surface was much roughened by the co-exposure. Density functional theory (DFT) calculations revealed the formation of superoxide and oxygen induced gap state, suggesting the creation of oxygen radicals by light illumination as a possible microscopic driving force for enhanced degradation.

First-principles equation of 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.

Probing atomic physics at ultrahigh pressure using laser-driven implosions.

Spectroscopic measurements of dense plasmas at billions of atmospheres provide tests to our fundamental understanding of how matter behaves at extreme conditions. Developing reliable atomic physics models at these conditions, benchmarked by experimental data, is crucial to an improved understanding of radiation transport in both stars and inertial fusion targets. However, detailed spectroscopic measurements at these conditions are rare, and traditional collisional-radiative equilibrium models, based on isolated-atom calculations and ad hoc continuum lowering models, have proved questionable at and beyond solid density. Here we report time-integrated and time-resolved x-ray spectroscopy measurements at several billion atmospheres using laser-driven implosions of Cu-doped targets. We use the imploding shell and its hot core at stagnation to probe the spectral changes of Cu-doped witness layer. These measurements indicate the necessity and viability of modeling dense plasmas with self-consistent methods like density-functional theory, which impact the accuracy of radiation transport simulations used to describe stellar evolution and the design of inertial fusion targets.

First-principles study of L-shell iron and chromium opacity at stellar interior temperatures.

Recently developed free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter has been applied for studying L-shell opacity of iron and chromium at T=182 eV. We use Mermin-Kohn-Sham density functional theory with a ground-state and a fully-temperature-dependent generalized gradient approximation exchange-correlation (XC) functionals. It is demonstrated that the role of XC at such a high-T is negligible due to the total free energy of interacting systems being dominated by the noninteracting free-energy term, in agreement with estimations for the homogeneous electron gas. Our DFT predictions are compared with the radiative emissivity and opacity of the dense plasma model, with the real-space Green's function method, and with experimental measurements. Good agreement is found between all three theoretical methods, and in the bound-continuum region for Cr when compared with the experiment, while the discrepancy between direct DFT calculations and the experiment for Fe remains essentially the same as for plasma-physics models.

Unraveling the intrinsic atomic physics behind x-ray absorption line shifts in warm dense silicon plasmas.

We present a free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter to cover a wide range of thermodynamic conditions and photon energies including the entire x-ray range. It uses Mermin-Kohn-Sham density functional theory with exchange-correlation (XC) thermal effects taken into account via a fully temperature dependent generalized gradient approximation XC functional. The methodology incorporates a combination of the ab initio molecular dynamics (AIMD) snapshotted Kubo-Greenwood optic data with a single atom in simulation cell calculations to close the photon energy gap between the L and K edges and extend the K-edge tail toward many-keV photon energies. This gap arises in the standard scheme due to a prohibitively large number of bands required for the Kubo-Greenwood calculations with AIMD snapshots. Kubo-Greenwood data on snapshots provide an accurate description of optic properties at low photon frequencies slightly beyond the L edge and x-ray absorption near edges structure (XANES) spectra, while data from periodic calculations with single atoms cover the tail regions beyond the edges. To demonstrate its applicability to mid-Z materials where the standard DFT-based approach is not computationally feasible, we have applied it to opacity calculations of warm dense silicon plasmas. These first-principles calculations revealed a very interesting phenomenon of redshift-to-blueshift in K-L (1s→2p) and K-edge absorptions along both isotherm and isochore, which are absent in most continuum-lowering models of traditional plasma physics. This new physics phenomenon can be attributed to the underlying competition between the screening of deeply bound core electrons and the screening of outer-shell electrons caused by warm-dense-plasma conditions. We further demonstrate that the ratio of 1s→2p to the K-edge x-ray absorptions can be used to characterize warm-dense-plasma conditions. Eventually, based on our absorption calculations, we have established a first-principles opacity table (FPOT) for silicon in a wide range of material densities and temperatures.

Nonempirical Semilocal Free-Energy Density Functional for Matter under Extreme Conditions.

Realizing the potential for predictive density functional calculations of matter under extreme conditions depends crucially upon having an exchange-correlation (XC) free-energy functional accurate over a wide range of state conditions. Unlike the ground-state case, no such functional exists. We remedy that with systematic construction of a generalized gradient approximation XC free-energy functional based on rigorous constraints, including the free-energy gradient expansion. The new functional provides the correct temperature dependence in the slowly varying regime and the correct zero-T, high-T, and homogeneous electron gas limits. Its accuracy in the warm dense matter regime is attested by excellent agreement of the calculated deuterium equation of state with reference path integral Monte Carlo results at intermediate and elevated T. Pressure shifts for hot electrons in compressed static fcc Al and for low-density Al demonstrate the combined magnitude of thermal and gradient effects handled well by this functional over a wide T range.

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.

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.

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

We compare the behavior of the finite-temperature Hartree-Fock model with that of thermal density functional theory using both ground-state and temperature-dependent approximate exchange functionals. The test system is bcc Li in the temperature-density regime of warm dense matter (WDM). In this exchange-only case, there are significant qualitative differences in results from the three approaches. Those differences may be important for Born-Oppenheimer molecular dynamics studies of WDM with ground-state approximate density functionals and thermal occupancies. Such calculations require reliable regularized potentials over a demanding range of temperatures and densities. By comparison of pseudopotential and all-electron results at T=0 K for small Li clusters of local bcc symmetry and bond lengths equivalent to high density bulk Li, we determine the density ranges for which standard projector augmented wave (PAW) and norm-conserving pseudopotentials are reliable. Then, we construct and use all-electron PAW data sets with a small cutoff radius that are valid for lithium densities up to at least 80 g/cm{3}.