Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userStructural properties and enthalpy of formation of magnesium hydride from quantum Monte Carlo calculations
146
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We have used diffusion Monte Carlo (DMC) calculations to study the structural properties of magnesium hydride (MgH$_2$), including the pressure-volume equation of state, the cohesive energy and the enthalpy of formation from magnesium bulk and hydrogen gas. The calculations employ pseudopotentials and B-spline basis sets to expand the single particle orbitals used to construct the trial wavefunctions. Extensive tests on system size, time step, and other sources of errors, performed on periodically repeated systems of up to 1050 atoms, show that all these errors together can be reduced to below 10 meV per formula unit. We find excellent agreement with the experiments for the equilibrium volume of both the Mg and the MgH$_2$ crystals. The cohesive energy of the Mg crystal is found to be 1.51(1) eV, and agrees perfectly with the experimental value of 1.51 eV. The enthalpy of formation of MgH$_2$ from Mg bulk and H$_2$ gas is found to be $0.85 pm 0.01$ eV/formula unit, or $82 pm 1$ kJ/mole, which is off the experimental one of $76.1 pm 1$ kJ/mole only by 6 kJ/mole. This shows that DMC can almost achieve chemical accuracy (1 kcal/mole) on this system. Density functional theory errors are shown to be much larger, and depend strongly on the functional employed.
rate research
Read More
We show how accurate benchmark values of the surface formation energy of crystalline lithium hydride can be computed by the complementary techniques of quantum Monte Carlo (QMC) and wavefunction-based molecular quantum chemistry. To demonstrate the high accuracy of the QMC techniques, we present a detailed study of the energetics of the bulk LiH crystal, using both pseudopotential and all-electron approaches. We show that the equilibrium lattice parameter agrees with experiment to within 0.03 %, which is around the experimental uncertainty, and the cohesive energy agrees to within around 10 meV per formula unit. QMC in periodic slab geometry is used to compute the formation energy of the LiH (001) surface, and we show that the value can be accurately converged with respect to slab thickness and other technical parameters. The quantum chemistry calculations build on the recently developed hierarchical scheme for computing the correlation energy of a crystal to high precision. We show that the hierarchical scheme allows the accurate calculation of the surface formation energy, and we present results that are well converged with respect to basis set and with respect to the level of correlation treatment. The QMC and hierarchical results for the surface formation energy agree to within about 1 %.
We have applied the many-body ab-initio diffusion quantum Monte Carlo (DMC) method to study Zn and ZnO crystals under pressure, and the energetics of the oxygen vacancy, zinc interstitial and hydrogen impurities in ZnO. We show that DMC is an accurate and practical method that can be used to characterize multiple properties of materials that are challenging for density functional theory approximations. DMC agrees with experimental measurements to within 0.3 eV, including the band-gap of ZnO, the ionization potential of O and Zn, and the atomization energy of O$_2$, ZnO dimer, and wurtzite ZnO. DMC predicts the oxygen vacancy as a deep donor with a formation energy of 5.0(2) eV under O-rich conditions and thermodynamic transition levels located between 1.8 and 2.5 eV from the valence band maximum. Our DMC results indicate that the concentration of zinc interstitial and hydrogen impurities in ZnO should be low under n-type, and Zn- and H-rich conditions because these defects have formation energies above 1.4 eV under these conditions. Comparison of DMC and hybrid functionals shows that these DFT approximations can be parameterized to yield a general correct qualitative description of ZnO. However, the formation energy of defects in ZnO evaluated with DMC and hybrid functionals can differ by more than 0.5 eV.
Zirconia (zirconium dioxide) and hafnia (hafnium dioxide) are binary oxides used in a range of applications. Because zirconium and hafnium are chemically equivalent, they have three similar polymorphs, and it is important to understand the properties and energetics of these polymorphs. However, while density functional theory calculations can get the correct energetic ordering, the energy differences between polymorphs depend very much on the specific density functional theory approach, as do other quantities such as lattice constants and bulk modulus. We have used highly accurate quantum Monte Carlo simulations to model the three zirconia and hafnia polymorphs. We compare our results for structural parameters, bulk modulus, and cohesive energy with results obtained from density functional theory calculations. We also discuss comparisons of our results with existing experimental data, in particular for structural parameters where extrapolation to zero temperature can be attempted. We hope our results of structural parameters as well as for cohesive energy and bulk modulus can serve as benchmarks for density-functional theory based calculations and as a guidance for future experiments.
We have used diffusion Monte Carlo (DMC) simulations to calculate the energy barrier for H$_2$ dissociation on the Mg(0001) surface. The calculations employ pseudopotentials and systematically improvable B-spline basis sets to expand the single particle orbitals used to construct the trial wavefunctions. Extensive tests on system size, time step, and other sources of errors, performed on periodically repeated systems of up to 550 atoms, show that all these errors together can be reduced to $sim 0.03$ eV. The DMC dissociation barrier is calculated to be $1.18 pm 0.03$ eV, and is compared to those obtained with density functional theory using various exchange-correlation functionals, with values ranging between 0.44 and 1.07 eV.
We analyze the problem of eliminating finite-size errors from quantum Monte Carlo (QMC) energy data. We demonstrate that both (i) adding a recently proposed [S. Chiesa et al., Phys. Rev. Lett. 97, 076404 (2006)] finite-size correction to the Ewald energy and (ii) using the model periodic Coulomb (MPC) interaction [L. M. Fraser et al., Phys. Rev. B 53, 1814 (1996); P. R. C. Kent et al., Phys. Rev. B 59, 1917 (1999); A. J. Williamson et al., Phys. Rev. B 55, 4851 (1997)] are good solutions to the problem of removing finite-size effects from the interaction energy in cubic systems, provided the exchange-correlation (XC) hole has converged with respect to system size. However, we find that the MPC interaction distorts the XC hole in finite systems, implying that the Ewald interaction should be used to generate the configuration distribution. The finite-size correction of Chiesa et al. is shown to be incomplete in systems of low symmetry. Beyond-leading-order corrections to the kinetic energy are found to be necessary at intermediate and high densities, and we investigate the effect of adding such corrections to QMC data for the homogeneous electron gas. We analyze finite-size errors in two-dimensional systems and show that the leading-order behavior differs from that which has hitherto been supposed. We compare the efficiency of different twist-averaging methods for reducing single-particle finite-size errors and we examine the performance of various finite-size extrapolation formulas. Finally, we investigate the system-size scaling of biases in diffusion QMC.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
