Electrical Conductivities and Low Frequency Opacities in the Warm Dense Matter Regime

In this article, we examine different approaches for calculating low frequency opacities in the warm dense matter regime. The relevance of the average-atom approximation and of different models for calculating opacities, such as the Ziman or Ziman–Evans models is discussed and the results compared to ab initio simulations. We begin by recalling the derivation of the Ziman–Evans resistivity from Kubo's linear response theory, using the local approximation to the solutions of the Lippmann–Schwinger equation. With the help of this approximation, we explicitly introduce an ionic structure factor into the Ziman formula, without resorting to the Born approximation. Both approaches involve the calculation of scattering phase shifts, which we integrate from Calogero equation with an adaptive step numerical scheme based on a Runge–Kutta–Merson solver. We show that if the atomic number $Z$ is not too large, integrating the phase shifts in this way is more time-efficient than using a classical Numerov-type scheme to solve the radial Schrödinger equation. Various approximations are explored for phase shifts to further improve computation time. For the Born approximation, we show that using Born phase shifts directly in the scattering cross-section gives more accurate results than with the integral formula based on the Fourier transform of the electron-ion potential. We also compare an analytical formula based on a Yukawa fit of the electron-ion potential to a numerical integration. The average-atom results are compared with DFT-based molecular dynamics simulations for aluminum in the dilute regime and for copper, aluminum and gold at solid density and different temperatures.