This roadmap presents the state-of-the-art, current challenges and near future developments anticipated in the thriving field of warm dense matter physics. Originating from strongly coupled plasma physics, high pressure physics and high energy density science, the warm dense matter physics community has recently taken a giant leap forward. This is due to spectacular developments in laser technology, diagnostic capabilities, and computer simulation techniques. Only in the last decade has it become possible to perform accurate enough simulations & experiments to truly verify theoretical results as well as to reliably design experiments based on predictions. Consequently, this roadmap discusses recent developments and contemporary challenges that are faced by theoretical methods, and experimental techniques needed to create and diagnose warm dense matter. A large part of this roadmap is dedicated to specific warm dense matter systems and applications in astrophysics, inertial confinement fusion and novel material synthesis.

We present a deorbitalization of the recent simplified, regularized Tao-Mo exchange functional () that is faithful to the parent functional. That is a major gain relative to our earlier deorbitalization which did poorly on molecular heats of formation (). The improvement arises from augmentation of the Mejia-Rodriguez and Trickey deorbitalization strategy () to use a smoothed replacement for the reduced density Laplacian (conventionally denoted q) obtained from that Laplacian itself. The augmentation also rationalizes the improvement obtained from the cutoff of q < 0 that was poorly understood at the time of the previous paper. The new scheme yields deorbitalized chemical region indicators that are much closer to those from the parent, orbital-dependent functional than were obtainable from the previous deorbitalization. It also replicates the good 3d elemental magnetization of the parent functional reasonably well.

Experimental results on two supramolecular complexes in which a CrIII or FeIII d-orbital single-ion magnet center is embedded between a pair of FeII spin-crossover moieties make those two complexes interesting as possible candidates for use in quantum information technologies. We report detailed computational results for their structure and electronic properties and use the resulting data to parametrize a spin Hamiltonian that facilitates comparison with experimental results and their interpretation. Consistent with experimental results on decoherence in [Fe(ox)3]@[Fe2L3]+, we find it to be easy-plane type while the [Cr(ox)3]@[Fe2L3]+ system is easy-axis type.

We propose computationally cheap and accurate approximants to the noninteracting kinetic energy density functional T s [ n ] by leveraging the simplicity and computational efficiency of the Wang-Teter functional [L.-W. Wang and M. Teter, Phys. Rev. B 45 , 13196 (1992)]. It depends on a single parameter, the average electron density rho 0 . We address limitations of the Wang-Teter functional, which include variational instabilities and inability to treat materials with finite band gaps. We introduce three physically motivated methods for determining rho 0 : DEN, minimizing the integrated difference of the self-consistent Wang-Teter electron density from the one from conventional Kohn-Sham density functional theory (DFT); KIN, minimizing the deviation between the Wang-Teter T s [ n ] and the exact value from conventional Kohn-Sham DFT; and finally ENE, minimizing the difference between the Wang-Teter and conventional Kohn-Sham DFT total energies. The crucial result of this work is that our approaches effectively mitigate the drawbacks of the Wang-Teter functional. We provide a thorough analysis of our methods and discuss their potential for large-scale simulations and as templates for density-dependent nonlocal functionals.

The revised, regularized Tao-Mo (rregTM) exchange-correlation density functional approximation (DFA) [A. Patra, S. Jana, and P. Samal, J. Chem. Phys. 153, 184112 (2020) and Jana et al., J. Chem. Phys. 155, 024103 (2021)] resolves the order-of-limits problem in the original TM formulation while preserving its valuable essential behaviors. Those include performance on standard thermochemistry and solid data sets that is competitive with that of the most widely explored meta-generalized-gradient-approximation DFAs (SCAN and r2SCAN) while also providing superior performance on elemental solid magnetization. Puzzlingly however, rregTM proved to be intractable for de-orbitalization via the approach of Mejía-Rodríguez and Trickey [Phys. Rev. A 96, 052512 (2017)]. We report investigation that leads to diagnosis of how the regularization in rregTM of the z indicator functions (z = the ratio of the von-Weizsäcker and Kohn-Sham kinetic energy densities) leads to non-physical behavior. We propose a simpler regularization that eliminates those oddities and that can be calibrated to reproduce the good error patterns of rregTM. We denote this version as simplified, regularized Tao-Mo, sregTM. We also show that it is unnecessary to use rregTM correlation with sregTM exchange: Perdew-Burke-Ernzerhof correlation is sufficient. The subsequent paper shows how sregTM enables some progress on de-orbitalization.

In Paper I [H. Francisco, A. C. Cancio, and S. B. Trickey, J. Chem. Phys. 159, 214102 (2023)], we gave a regularization of the Tao-Mo exchange functional that removes the order-of-limits problem in the original Tao-Mo form and also eliminates the unphysical behavior introduced by an earlier regularization while essentially preserving compliance with the second-order gradient expansion. The resulting simplified, regularized (sregTM) functional delivers performance on standard molecular and solid state test sets equal to that of the earlier revised, regularized Tao-Mo functional. Here, we address de-orbitalization of that new sregTM into a pure density functional. We summarize the failures of the Mejía-Rodríguez and Trickey de-orbitalization strategy [Phys. Rev. A 96, 052512 (2017)] when used with both versions. We discuss how those failures apparently arise in the so-called z' indicator function and in substitutes for the reduced density Laplacian in the parent functionals. Then, we show that the sregTM functional can be de-orbitalized somewhat well with a rather peculiarly parameterized version of the previously used deorbitalizer. We discuss, briefly, a de-orbitalization that works in the sense of reproducing error patterns but that apparently succeeds by cancelation of major qualitative errors associated with the de-orbitalized indicator functions α and z, hence, is not recommended. We suggest that the same issue underlies the earlier finding of comparatively mediocre performance of the de-orbitalized Tao-Perdew-Staroverov-Scuseri functional. Our work demonstrates that the intricacy of such two-indicator functionals magnifies the errors introduced by the Mejía-Rodríguez and Trickey de-orbitalization approach in ways that are extremely difficult to analyze and correct.

Calculation of transition temperatures T1/2 for thermally driven spin-crossover in condensed phases is challenging, even with sophisticated state-of-the-art density functional approximations. The first issue is the accuracy of the adiabatic crossover energy difference ΔEHL between the low- and high-spin states of the bistable metal-organic complexes. The other is the proper inclusion of entropic contributions to the Gibbs free energy from the electronic and vibrational degrees of freedom. We discuss the effects of treatments of both contributions upon the calculation of thermochemical properties for a set of 20 spin-crossover materials using a Hubbard-U correction obtained from a reference ensemble spin-state. The U values obtained from a simplest bimolecular representation may overcorrect, somewhat, the ΔEHL values, hence giving somewhat excessive reduction of the T1/2 results with respect to their U = 0 values in the crystalline phase. We discuss the origins of the discrepancies by analyzing different sources of uncertainties. By use of a first-coordination-sphere approximation and the assumption that vibrational contributions from the outermost atoms in a metal-organic complex are similar in both low- and high-spin states, we achieve T1/2 results with the low-cost, widely used PBE generalized gradient density functional approximation comparable to those from the more costly, more sophisticated r2SCAN meta-generalized gradient approximation. The procedure is promising for use in high-throughput materials screening, because it combines rather low computational effort requirements with freedom from user manipulation of parameters.

Though calculations based on density functional theory (DFT) are used remarkably widely in chemistry, physics, materials science, and biomolecular research and though the modern form of DFT has been studied for almost 60 years, some mathematical problems remain. From a physical science perspective, it is far from clear whether those problems are of major import. For context, we provide an outline of the basic structure of DFT as it is presented and used conventionally in physical sciences, note some unresolved mathematical difficulties with those conventional demonstrations, then pose several questions regarding both the time-independent and time-dependent forms of DFT that could benefit from attention in applied mathematics. Progress on any of these would aid in development of better approximate functionals and in interpretation of DFT.

Effects of finite-temperature quasiparticle self-energy corrections to x-ray absorption spectra are investigated within the finite-temperature quasiparticle local density GW approximation up to temperatures $T$ of order the Fermi temperature. To facilitate the calculations, we parametrize the quasiparticle self-energy using low-order polynomial fits. We show that temperature-driven decrease in the electron lifetime substantially broadens the spectra in the near-edge region with increasing $T$. However, the quasiparticle shift is most strongly modified near the onset of plasmon excitations.

We report APW+lo (augmented plane wave plus local orbital) density functional theory (DFT) calculations of large molecular systems using the domain specific SIRIUS multi-functional DFT package. The APW and FLAPW (full potential linearized APW) task and data parallelism options and the advanced eigen-system solver provided by SIRIUS can be exploited for performance gains in ground state Kohn-Sham calculations on large systems. This approach is distinct from our prior use of SIRIUS as a library backend to another APW+lo or FLAPW code. We benchmark the code and demonstrate performance on several magnetic molecule and metal organic framework systems. We show that the SIRIUS package in itself is capable of handling systems as large as a several hundred atoms in the unit cell without having to make technical choices that result in the loss of accuracy with respect to that needed for the study of magnetic systems.

Kohn-Sham Density Functional Theory (KSDFT) is the most widely used electronic structure method in chemistry, physics, and materials science, with thousands of calculations cited annually. This ubiquity is rooted in the favorable accuracy vs cost balance of KSDFT. Nonetheless, the ambitions and expectations of researchers for use of KSDFT in predictive simulations of large, complicated molecular systems are confronted with an intrinsic computational cost-scaling challenge. Particularly evident in the context of first-principles molecular dynamics, the challenge is the high cost-scaling associated with the computation of the Kohn-Sham orbitals. Orbital-free DFT (OFDFT), as the name suggests, circumvents entirely the explicit use of those orbitals. Without them, the structural and algorithmic complexity of KSDFT simplifies dramatically and near-linear scaling with system size irrespective of system state is achievable. Thus, much larger system sizes and longer simulation time scales (compared to conventional KSDFT) become accessible; hence, new chemical phenomena and new materials can be explored. In this review, we introduce the historical contexts of OFDFT, its theoretical basis, and the challenge of realizing its promise via approximate kinetic energy density functionals (KEDFs). We review recent progress on that challenge for an array of KEDFs, such as one-point, two-point, and machine-learnt, as well as some less explored forms. We emphasize use of exact constraints and the inevitability of design choices. Then, we survey the associated numerical techniques and implemented algorithms specific to OFDFT. We conclude with an illustrative sample of applications to showcase the power of OFDFT in materials science, chemistry, and physics.

We study intramolecular electron transfer in the single-molecule magnetic complex [Mn12O12(O2CR)16 (H2O)4] for R = -H, -CH3, -CHCl2, -C6H5, and -C6H4F ligands as a mechanism for switching of the molecular dipole moment. Energetics is obtained using the density functional theory (DFT) with onsite Coulomb energy correction (DFT + U). Lattice distortions are found to be critical for localizing an extra electron on one of the easy sites on the outer ring in which localized states can be stabilized. We find that the lowest-energy path for charge transfer is for the electron to go through the center via superexchange-mediated tunneling. The energy barrier for such a path ranges from 0.4 to 54 meV depending on the ligands and the isomeric form of the complex. The electric field strength needed to move the charge from one end to the other, thus reversing the dipole moment, is 0.01-0.04 V/Å.

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 300 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 776 entries, the paper represents a broad snapshot of DFT, anno 2022.

Full-potential linearized augmented plane wave (LAPW) and APW plus local orbital (APW+lo) codes differ widely in both their user interfaces and in capabilities for calculations and analysis beyond their common central task of all-electron solution of the Kohn–Sham equations. However, that common central task opens a possible route to performance enhancement, namely to offload the basic LAPW/APW+lo algorithms to a library optimized purely for that purpose. To explore that opportunity, we have interfaced the Exciting-Plus (“EP”) LAPW/APW+lo DFT code with the highly optimized SIRIUS multi-functional DFT package. This simplest realization of the separation of concerns approach yields substantial performance over the base EP code via additional task parallelism without significant change in the EP source code or user interface. We provide benchmarks of the interfaced code against the original EP using small bulk systems, and demonstrate performance on a spin-crossover molecule and magnetic molecule that are of size and complexity at the margins of the capability of the EP code itself.

Current spin-crossover (SCO) energy calculations depend on nearly artisanal skill in picking quantum mechanical approximations and computational methods. That is incompatible with automated (work-flow-driven) screening. An acceptable methodology must be quantum mechanically sound, yield both basic structure and property values, and accurate ΔEHL without steering or tuning. Cost vs. accuracy causes focus on density functional theory (DFT). We show by a near-exhaustive study of schemes for calculating the basic molecular high-low spin energy difference, ΔEHL, that presently there is no combination of a constraint-based, non-empirical density functional approximation (DFA) and a set of well-defined semi-empirical corrections to it adequate for such a protocol. Somewhat successful hybrid DFA calculations of ΔEHL are too costly for high-throughput screening of condensed phases. Lower-cost alternatives combine a generalized gradient approximation (GGA) DFA with a Hubbard-U correction (DFT + U). But we show that neither U=0 nor any currently available unsteered U calculation gives a decent ΔEHL value for [Mn(taa)] without also degrading molecular structure or property predictions. Moving to the SCAN meta-GGA does not solve the problem. The revised-restored SCAN (r2SCAN) meta-GGA together with its deorbitalized version r2SCAN-L give improved but not wholly satisfactory results. We also document and diagnose several non-obvious technical and procedural sensitivities and inter-code differences. In addition to being a formidable challenge to DFA development, the lack we delineate is a major impediment to progress in the development of quantum materials and spintronics.

High-throughput searches for spin-crossover molecules require Hubbard -U corrections to common density functional exchange-correlation (XC) approximations. However, the Ueff values obtained from linear response or based on previous studies overcorrect the spin crossover energies. We demonstrate that employing a linearly mixed ensemble average spin state as the reference configuration for the linear response calculation of Ueff resolves this issue. Validation on a commonly used set of spin-crossover complexes shows that these ensemble Ueff values consistently are smaller than those calculated directly on a pure spin state, irrespective of whether that be low-or high-spin. Adiabatic crossover energies using this methodology for a generalized gradient approximation XC functional are closer to the expected target energy range than with conventional Ueffvalues. Based on the observation that the Ueffcorrection is similar for different complexes that share transition metals with the same oxidation state, we devise a set of recommended averaged Ueff values for high-throughput calculations.

The deviations from linearity of the energy as a function of the number of electrons that arise with current approximations to the exchange-correlation (XC) energy functional have important consequences for the frontier eigenvalues of molecules and the corresponding valence-band maxima for solids. In this work, we present an analysis of the exact theory that allows one to infer the effects of such approximations on the highest occupied and lowest unoccupied molecular orbital eigenvalues. Then, we show the importance of the asymptotic behavior of the XC potential in the generalized gradient approximation (GGA) in the case of the NCAPR functional (nearly correct asymptotic potential revised) for determining the shift of the frontier orbital eigenvalues toward the exact values. Thereby we establish a procedure at the GGA level of refinement that allows one to make a single calculation to determine the ionization potential, the electron affinity, and the hardness of molecules (and its solid counterpart, the bandgap) with an accuracy equivalent to that obtained for those properties through energy differences, a procedure that requires three calculations. For solids, the accuracy achieved for the bandgap lies rather close to that which is obtained through hybrid XC energy functionals, but those also demand much greater computational effort than what is required with the simple NCAPR GGA calculation.

We report the interfacing of the Exciting-Plus (EP) FLAPW DFT code with the SIRIUS multi-functional DFT library. Use of the SIRIUS library enhances EP with additional task parallelism in ground state DFT calculations. Without significant change in the EP source code, the additional eigensystem solver method from the SIRIUS library can be exploited for performance gains in diagonalizing the Kohn-Sham Hamiltonian. We benchmark the interfaced code against the original EP using small bulk systems, and then demonstrate performance on much larger molecular magnet systems that are well beyond the capability of the original EP code.