No Arabic abstract
We present an accurate study of the static-nucleus electronic energy band gap of solid molecular hydrogen at high pressure. The excitonic and quasiparticle gaps of the $C2/c$, $Pc$, $Pbcn$, and $P6_3/m$ structures at pressures of 250, 300, and 350~GPa are calculated using the fixed-node diffusion quantum Monte Carlo (DMC) method. The difference between the mean-field and many-body band gaps at the same density is found to be almost independent of system size and can therefore be applied as a scissor correction to the mean-field gap of an infinite system to obtain an estimate of the many-body gap in the thermodynamic limit. By comparing our static-nucleus DMC energy gaps with available experimental results, we demonstrate the important role played by nuclear quantum effects in the electronic structure of solid hydrogen. Our DMC results suggest that the metallization of high-pressure solid hydrogen occurs via a structural phase transition rather than band gap closure.
Being the simplest element with just one electron and proton the electronic structure of the Hydrogen atom is known exactly. However, this does not hold for the complex interplay between them in a solid and in particular not at high pressure that is known to alter the crystal as well as the electronic structure. Back in 1935 Wigner and Huntington predicted that at very high pressure solid molecular hydrogen would dissociate and form an atomic solid that is metallic. In spite of intense research efforts the experimental realization, as well as the theoretical determination of the crystal structure has remained elusive. Here we present a computational study showing that the distorted hexagonal P6$_3$/m structure is the most likely candidate for Phase III of solid hydrogen. We find that the pairing structure is very persistent and insulating over the whole pressure range, which suggests that metallization due to dissociation may precede eventual bandgap closure. Due to the fact that this not only resolve one of major disagreement between theory and experiment, but also excludes the conjectured existence of phonon-driven superconductivity in solid molecular hydrogen, our results involve a complete revision of the zero-temperature phase diagram of Phase III.
We have studied solid hydrogen under pressure at low temperatures. With increasing pressure we observe changes in the sample, going from transparent, to black, to a reflective metal, the latter studied at a pressure of 495 GPa. We have measured the reflectance as a function of wavelength in the visible spectrum finding values as high as 0.90 from the metallic hydrogen. We have fit the reflectance using a Drude free electron model to determine the plasma frequency of 30.1 eV at T= 5.5 K, with a corresponding electron carrier density of 6.7x1023 particles/cm3, consistent with theoretical estimates. The properties are those of a metal. Solid metallic hydrogen has been produced in the laboratory.
We present an accurate computational study of the electronic structure and lattice dynamics of solid molecular hydrogen at high pressure. The band-gap energies of the $C2/c$, $Pc$, and $P6_3/m$ structures at pressures of 250, 300, and 350 GPa are calculated using the diffusion quantum Monte Carlo (DMC) method. The atomic configurations are obtained from ab-initio path-integral molecular dynamics (PIMD) simulations at 300 K and 300 GPa to investigate the impact of zero-point energy and temperature-induced motion of the protons including anharmonic effects. We find that finite temperature and nuclear quantum effects reduce the band-gaps substantially, leading to metallization of the $C2/c$ and $Pc$ phases via band overlap; the effect on the band-gap of the $P6_3/m$ structure is less pronounced. Our combined DMC-PIMD simulations predict that there are no excitonic or quasiparticle energy gaps for the $C2/c$ and $Pc$ phases at 300 GPa and 300 K. Our results also indicate a strong correlation between the band-gap energy and vibron modes. This strong coupling induces a band-gap reduction of more than 2.46 eV in high-pressure solid molecular hydrogen. Comparing our DMC-PIMD with experimental results available, we conclude that none of the structures proposed is a good candidate for phases III and IV of solid hydrogen.
We report a numerical study of the equation of state of crystalline body-centered-cubic (BCC) hydrogen, tackled with a variety of complementary many-body wave function methods. These include continuum stochastic techniques of fixed-node diffusion and variational quantum Monte Carlo, and the Hilbert space stochastic method of full configuration-interaction quantum Monte Carlo. In addition, periodic coupled-cluster methods were also employed. Each of these methods is underpinned with different strengths and approximations, but their combination in order to perform reliable extrapolation to complete basis set and supercell size limits gives confidence in the final results. The methods were found to be in good agreement for equilibrium cell volumes for the system in the BCC phase, with a lattice parameter of 3.307 Bohr.
We study the effect of temperature up to 1000K on the structure of dense molecular para-hydrogen and ortho-deuterium, using the path-integral Monte Carlo method. We find a structural phase transition from orientationally disordered hexagonal close packed (hcp) to an orthorhombic structure of Cmca symmetry before melting. The transition is basically induced by thermal fluctuations, but quantum fluctuations of protons (deuterons) are important in determining the transition temperature through effectively hardening the intermolecular interaction. We estimate the phase line between hcp and Cmca phases as well as the melting line of the Cmca solid.