Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userGW space-time method: Energy band-gap of solid hydrogen
184
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We implement the GW space-time method at finite temperatures, in which the Greens function G and the screened Coulomb interaction W are represented in the real space on a suitable mesh and in imaginary time in terms of Chebyshev polynomials, paying particular attention to controlling systematic errors of the representation. Having validated the technique by the canonical application to silicon and germanium, we apply it to calculation of band gaps in hexagonal solid hydrogen with the bare Greens function obtained from density functional approximation and the interaction screened within the random phase approximation (RPA). The results, obtained from the asymptotic decay of the full Greens function without resorting to analytic continuation, suggest that the solid hydrogen above 250 GPa can not adopt the hexagonal-closed-pack (hcp) structure. The demonstrated ability of the method to store the full G and W functions in memory with sufficient accuracy is crucial for its subsequent extensions to include higher orders of the diagrammatic series by means of diagrammatic Monte Carlo algorithms.
rate research
Read More
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 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.
Excited-state calculations, notably for quasiparticle band structures, are nowadays routinely performed within the GW approximation for the electronic self-energy. Nevertheless, certain numerical approximations and simplifications are still employed in practice to make the computations feasible. An important aspect for periodic systems is the proper treatment of the singularity of the screened Coulomb interaction in reciprocal space, which results from the slow 1/r decay in real space. This must be done without introducing artificial interactions between the quasiparticles and their periodic images in repeated cells, which occur when integrals of the screened Coulomb interaction are discretised in reciprocal space. An adequate treatment of both aspects is crucial for a numerically stable computation of the self-energy. In this article we build on existing schemes for isotropic screening and present an extension for anisotropic systems. We also show how the contributions to the dielectric function arising from the non-local part of the pseudopotentials can be computed efficiently. These improvements are crucial for obtaining a fast convergence with respect to the number of points used for the Brillouin zone integration and prove to be essential to make GW calculations for strongly anisotropic systems, such as slabs or multilayers, efficient.
Using the GW approximation, we study the electronic structure of the recently synthesized hydrogenated graphene, named graphane. For both conformations, the minimum band gap is found to be direct at the $Gamma$ point, and it has a value of 5.4 eV in the stable chair conformation, where H atoms attach C atoms alternatively on opposite sides of the two dimensional carbon network. In the meta-stable boat conformation the energy gap is 4.9 eV. Then, using a supercell approach, the electronic structure of graphane was modified by introducing either an hydroxyl group or an H vacancy. In this last case, an impurity state appears at about 2 eV above the valence band maximum.
We use the diffusion quantum Monte Carlo to revisit the enthalpy-pressure phase diagram of the various products from the different proposed decompositions of H$_2$S at pressures above 150~GPa. Our results entails a revision of the ground-state enthalpy-pressure phase diagram. Specifically, we find that the C2/c HS$_2$ structure is persistent up to 440~GPa before undergoing a phase transition into the C2/m phase. Contrary to density functional theory, our calculations suggest that the C2/m phase of HS is more stable than the I4$_1$/amd HS structure over the whole pressure range from 150 to 400 GPa. Moreover, we predict that the Im-3m phase is the most likely candidate for H$_3$S, which is consistent with recent experimental x-ray diffraction measurements.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
