Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userFate of density functional theory in high-pressure solid hydrogen
510
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
This paper investigates some of the successes and failures of density functional theory in the study of high-pressure solid hydrogen at low temperature. We calculate the phase diagram, metallization pressure, phonon spectrum, and proton zero-point energy using three popular exchange-correlation functionals: the local density approximation (LDA), the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation, and the semi-local Becke-Lee-Yang-Parr (BLYP) functional. We focus on the solid molecular P$6_3$/m, C2/c, Cmca-12, and Cmca structures in the pressure range from $100<P<500$ GPa over which phases I, II and III are observed experimentally. At the static level of theory, in which proton zero-point energy is ignored, the LDA, PBE and BLYP functionals give very different structural transition and metallization pressures, with the BLYP phase diagram in better agreement with experiment. Nevertheless, all three functionals provide qualitatively the same information about the band gaps of the four structures and the phase transitions between them. Going beyond the static level, we find that the frequencies of the vibron modes observed above 3000 cm$^{-1}$ depend strongly on the choice of exchange-correlation functional, although the low-frequency part of the phonon spectrum is little affected. The largest and smallest values of the proton zero-point energy, obtained using the BLYP and LDA functionals, respectively, differ by more than 10 meV/proton. Including the proton zero-point energy calculated from the phonon spectrum within the harmonic approximation improves the agreement of the BLYP and PBE phase diagrams with experiment. Taken as a whole, our results demonstrate the inadequacy of mean-field-like density functional calculations of solid molecular hydrogen in phases I, II and III and emphasize the need for more sophisticated methods.
rate research
Read More
We use the diffusion quantum Monte Carlo (DMC) method to calculate the ground state phase diagram of solid molecular hydrogen and examine the stability of the most important insulating phases relative to metallic crystalline molecular hydrogen. We develop a new method to account for finite-size errors by combining the use of twist-averaged boundary conditions with corrections obtained using the Kwee-Zhang-Krakauer (KZK) functional in density functional theory. To study band-gap closure and find the metallization pressure, we perform accurate quasi-particle many-body calculations using the $GW$ method. In the static approximation, our DMC simulations indicate a transition from the insulating Cmca-12 structure to the metallic Cmca structure at around 375 GPa. The $GW$ band gap of Cmca-12 closes at roughly the same pressure. In the dynamic DMC phase diagram, which includes the effects of zero-point energy, the Cmca-12 structure remains stable up to 430 GPa, well above the pressure at which the $GW$ band gap closes. Our results predict that the semimetallic state observed experimentally at around 360 GPa [Phys. Rev. Lett. {bf 108}, 146402 (2012)] may correspond to the Cmca-12 structure near the pressure at which the band gap closes. The dynamic DMC phase diagram indicates that the hexagonal close packed $P6_3/m$ structure, which has the largest band gap of the insulating structures considered, is stable up to 220 GPa. This is consistent with recent X-ray data taken at pressures up to 183 GPa [Phys. Rev. B {bf 82}, 060101(R) (2010)], which also reported a hexagonal close packed arrangement of hydrogen molecules.
We investigate the van der Waals interactions in solid molecular hydrogen structures. We calculate enthalpy and the Gibbs free energy to obtain zero and finite temperature phase diagrams, respectively. We employ density functional theory (DFT) to calculate the electronic structure and Density functional perturbation theory (DFPT) with van der Waals (vdW) functionals to obtain phonon spectra. We focus on the solid molecular $C2/c$, $Cmca$-12, $P6_3/m$, $Cmca$, and $Pbcn$ structures within the pressure range of 200 $<$ P $<$ 450 GPa. We propose two structures of the $C2/c$ and $Pbcn$ for phase III which are stabilized within different pressure range above 200 GPa. We find that vdW functionals have a big effect on vibrations and finite-temperature phase stability, however, different vdW functionals have different effects. We conclude that, in addition to the vdW interaction, a correct treatment of the high charge gradient limit is essential. We show that the dependence of molecular bond-lengths on exchange-correlation also has a considerable influence on the calculated metallization pressure, introducing errors of up to 100GPa.
An explicit expression for the quadratic density-response function of a many-electron system is obtained in the framework of the time-dependent density-functional theory, in terms of the linear and quadratic density-response functions of noninteracting Kohn-Sham electrons and functional derivatives of the time-dependent exchange-correlation potential. This is used to evaluate the quadratic stopping power of a homogeneous electron gas for slow ions, which is demonstrated to be equivalent to that obtained up to second order in the ion charge in the framework of a fully nonlinear scattering approach. Numerical calculations are reported, thereby exploring the range of validity of quadratic-response theory.
We reassess the phase diagram of high-pressure solid hydrogen using mean-field and many-body wave function based approaches to determine the nature of phase III of solid hydrogen. To discover the best candidates for phase III, density functional theory calculations within the meta-generalized gradient approximation by means of the strongly constrained and appropriately normed (SCAN) semilocal density functional are employed. We study eleven molecular structures with different symmetries, which are the most competitive phases, within the pressure range of 100 to 500~GPa. The SCAN phase diagram predicts that the $C2/c-24$ and $P6_122-36$ structures are the best candidates for phase III with an energy difference of less than 1~meV/atom. To verify the stability of the competitive insulator structures of $C2/c-24$ and $P6_122-36$, we apply the diffusion Monte Carlo (DMC) method to optimise the percentage $alpha$ of exact-exchange in the trial many-body wave function. We found that the optimised $alpha$ equals to $40 %$, and denote the corresponding exchange and correlation functional as PBE1. The energy gain with respect to the well-known hybrid functional PBE0, where $alpha = 25%$, varies with density and structure. The PBE1-DMC enthalpy-pressure phase diagram predicts that the $P6_122-36$ structure is stable up to 210~GPa, where it transforms to the $C2/c-24$. Hence, we predict that the phase III of high-pressure solid hydrogen is polymorphic.
A theoretical study is reported of the molecular-to-atomic transition in solid hydrogen at high pressure. We use the diffusion quantum Monte Carlo method to calculate the static lattice energies of the competing phases and a density-functional-theory-based vibrational self-consistent field method to calculate anharmonic vibrational properties. We find a small but significant contribution to the vibrational energy from anharmonicity. A transition from the molecular Cmca-12 direct to the atomic I4_1/amd phase is found at 374 GPa. The vibrational contribution lowers the transition pressure by 91 GPa. The dissociation pressure is not very sensitive to the isotopic composition. Our results suggest that quantum melting occurs at finite temperature.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
