TY - JOUR
T1 - Transition to metallization in warm dense helium-hydrogen mixtures using stochastic density functional theory within the Kubo-Greenwood formalism
AU - Cytter, Yael
AU - Rabani, Eran
AU - Neuhauser, Daniel
AU - Preising, Martin
AU - Redmer, Ronald
AU - Baer, Roi
N1 - Publisher Copyright:
© 2019 American Physical Society.
PY - 2019/11/1
Y1 - 2019/11/1
N2 - The Kubo-Greenwood (KG) formula is often used in conjunction with Kohn-Sham (KS) density functional theory (DFT) to compute the optical conductivity, particularly for warm dense matter. For applying the KG formula, all KS eigenstates and eigenvalues up to an energy cutoff are required and thus the approach becomes expensive, especially for high temperatures and large systems, scaling cubically with both system size and temperature. Here, we develop an approach to calculate the KS conductivity within the stochastic DFT framework, which requires knowledge only of the KS Hamiltonian but not its eigenstates and values. We show that the computational effort associated with the method scales linearly with system size and reduces in proportion to the temperature, unlike the cubic increase with traditional deterministic approaches. In addition, we find that the method allows an accurate description of the entire spectrum, including the high-frequency range, unlike the deterministic method which is compelled to introduce a high-frequency cutoff due to memory and computational time constraints. We apply the method to helium-hydrogen mixtures in the warm dense matter regime at temperatures of ∼60kK and find that the system displays two conductivity phases, where a transition from nonmetal to metal occurs when hydrogen atoms constitute ∼0.3 of the total atoms in the system.
AB - The Kubo-Greenwood (KG) formula is often used in conjunction with Kohn-Sham (KS) density functional theory (DFT) to compute the optical conductivity, particularly for warm dense matter. For applying the KG formula, all KS eigenstates and eigenvalues up to an energy cutoff are required and thus the approach becomes expensive, especially for high temperatures and large systems, scaling cubically with both system size and temperature. Here, we develop an approach to calculate the KS conductivity within the stochastic DFT framework, which requires knowledge only of the KS Hamiltonian but not its eigenstates and values. We show that the computational effort associated with the method scales linearly with system size and reduces in proportion to the temperature, unlike the cubic increase with traditional deterministic approaches. In addition, we find that the method allows an accurate description of the entire spectrum, including the high-frequency range, unlike the deterministic method which is compelled to introduce a high-frequency cutoff due to memory and computational time constraints. We apply the method to helium-hydrogen mixtures in the warm dense matter regime at temperatures of ∼60kK and find that the system displays two conductivity phases, where a transition from nonmetal to metal occurs when hydrogen atoms constitute ∼0.3 of the total atoms in the system.
UR - http://www.scopus.com/inward/record.url?scp=85075267380&partnerID=8YFLogxK
U2 - 10.1103/PhysRevB.100.195101
DO - 10.1103/PhysRevB.100.195101
M3 - ???researchoutput.researchoutputtypes.contributiontojournal.article???
AN - SCOPUS:85075267380
SN - 2469-9950
VL - 100
JO - Physical Review B
JF - Physical Review B
IS - 19
M1 - 195101
ER -