We have performed density functional theory (DFT) based calculations for aluminum in extreme conditions of both pressure and temperature, up to 5 times compressed ambient density and over 1 000 000 K in temperature. In order to cover such a domain, DFT methods including phonon calculations, quantum molecular dynamics, and orbital-free DFT are employed. The results are then used to construct a SESAME equation of state for the aluminum 1100 alloy, encompassing the fcc, hcp and bcc solid phases as well as the liquid regime. We provide extensive comparison with experiment and based on this we also provide a slightly modified equation of state for the aluminum 6061 alloy.

We propose a means for constructing highly accurate equations of state (EOS) for elemental solids and liquids essentially from first principles, based upon a particular decomposition of the underlying condensed matter Hamiltonian for the nuclei and electrons. We also point out that at low pressures the neglect of anharmonic and electronphonon terms, both contained in this formalism, results in errors of less than 5% in the thermal parts of the thermodynamic functions. Then we explicitly display the forms of the remaining terms in the EOS, commenting on the use of experiment and electronic structure theory to evaluate them. We also construct an EOS for Aluminum and compare the resulting Hugoniot with data up to 5 Mbar, both to illustrate our method and to see whether the approximation of neglecting anharmonicity et al. remains viable to such high pressures. We find a level of agreement with experiment that is consistent with the low-pressure results.

The shock Hugoniot for lithium 6 deuteride ((6)LiD) was calculated via first principles using Kohn-Sham density-functional theory molecular dynamics (KSMD) for temperatures of 0.5-25 eV. The upper limit of 25 eV represents a practical limit where KSMD is no longer computationally feasible due to the number of electronic bands which are required to be populated. To push the Hugoniot calculations to higher temperatures we make use of orbital-free density-functional theory molecular dynamics (OFMD). Thomas-Fermi-Dirac-based OFMD gives a poor description of the electronic structure at low temperatures so the initial state is not well defined. We propose a method of bootstrapping the Hugoniot from OFMD to the Hugoniot from KSMD between 10 and 20 eV, where the two methods are in agreement. The combination of KSMD and OFMD allows construction of a first-principles Hugoniot from the initial state to 1000 eV. Theoretical shock-compression results are in good agreement with available experimental data and exhibit the appropriate high-temperature limits. We show that a unified KSMD-OFMD Hugoniot can be used to assess the quality of the existing equation-of-state (EOS) models and inform better EOS models based on justifiable physics.

Through quantum molecular dynamics (QMD), utilizing both Kohn-Sham (orbital-based) and orbital-free density functional theory, we calculate the equation of state of warm dense iron in the density range 7-30g/cm^{3} and temperatures from 1 to 100 eV. A critical examination of the iron pseudopotential is made, from which we find a significant improvement at high pressure to the previous QMD calculations of Wang et al. [Phys. Rev. E 89, 023101 (2014)10.1103/PhysRevE.89.023101]. Our results also significantly extend the ranges of density and temperature that were attempted in that prior work. We calculate the shock Hugoniot and find very good agreement with experimental results to pressures over 20 TPa. These results are then incorporated with previous studies to generate a five-phase equation of state for iron.

The liquid regime equation of state of silicon dioxide SiO2 is calculated via quantum molecular dynamics in the density range 5 to 15 g/cc and with temperatures from 0.5 to 100 eV, including the α-quartz and stishovite phase Hugoniot curves. Below 8 eV calculations are based on Kohn-Sham density functional theory (DFT), above 8 eV a new orbital-free DFT formulation, presented here, based on matching Kohn-Sham DFT calculations is employed. Recent experimental shock data is found to be in very good agreement with the current results. Finally both experimental and simulation data are used in constructing a new liquid regime equation of state table for SiO2.