First-principles molecular dynamics of liquid alkali metals based on the quantal hypernetted chain theory

First principles molecular dynamics of dense plasmas

Ab initio molecular dynamics calculations are performed for the equation of state of aluminum, spanning condensed matter and dense plasma regimes. Electronic exchange and correlation are included with either a zero- or finite-temperature local density approximation potential. Standard methods are extended to above the Fermi temperature by using final state pseudopotentials to describe thermally excited ion cores. The predicted Hugoniot equation of state agrees well with earlier plasma theories and with experiment for temperatures from 0 to 3 x 10(6) K.

Nucleus-electron model for states changing from a liquid metal to a plasma and the Saha equation

We extend the quantal hypernetted-chain (QHNC) method, which has been proved to yield accurate results for liquid metals, to treat a partially ionized plasma. In a plasma, the electrons change from a quantum to a classical fluid gradually with increasing temperature; the QHNC method applied to the electron gas is in fact able to provide the electron-electron correlation at an arbitrary temperature. As an illustrating example of this approach, we investigate how liquid rubidium becomes a plasma by increasing the temperature from 0 to 30 eV at a fixed normal ion density 1.03x10(22)/cm(3). The electron-ion radial distribution function (RDF) in liquid Rb has distinct inner-core and outer-core parts. Even at a temperature of 1 eV, this clear distinction remains as a characteristic of a liquid metal. At a temperature of 3 eV, this distinction disappears, and rubidium becomes a plasma with the ionization 1.21. The temperature variations of bound levels in each ion and the average ionization are calculated in Rb plasmas at the same time. Using the density-functional theory, we also derive the Saha equation applicable even to a high-density plasma at low temperatures. The QHNC method provides a procedure to solve this Saha equation with ease by using a recursive formula; the charge population of differently ionized species are obtained in Rb plasmas at several temperatures. In this way, it is shown that, with the atomic number as the only input, the QHNC method produces the average ionization, the electron-ion and ion-ion RDF's, and the charge population that are consistent with the atomic structure of each ion for a partially ionized plasma.