Path integral Monte Carlo simulations of warm dense aluminum
Here, we perform first-principles path integral Monte Carlo (PIMC) and density functional theory molecular dynamics (DFT-MD) calculations to explore warm dense matter states of aluminum. Our equation of state (EOS) simulations cover a wide density-temperature range of $0.1-32.4$~g$\,$cm$^{-3}$ and $10^4-10^8$~K. Since PIMC and DFT-MD accurately treat effects of the atomic shell structure, we find two compression maxima along the principal Hugoniot curve attributed to K-shell and L-shell ionization. The results provide a benchmark for widely used EOS tables, such as SESAME, QEOS, and models based on Thomas-Fermi and average-atom techniques. A subsequent multi-shock analysis provides a quantitative assessment for how much heating occurs relative to an isentrope in multi-shock experiments. Finally, we compute heat capacity, pair-correlation functions, the electronic density of states, and Z to reveal the evolution of the plasma structure and ionization behavior.