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 (DFT) 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₂, 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.

We present a many-body diffusion quantum Monte Carlo (DMC) study of the bulk and defect properties of NiO. We find excellent agreement with experimental values, within 0.3%, 0.6%, and 3.5% for the lattice constant, cohesive energy, and bulk modulus, respectively. The quasiparticle bandgap was also computed, and the DMC result of 4.72 (0.17) eV compares well with the experimental value of 4.3 eV. Furthermore, DMC calculations of excited states at the L, Z, and the gamma point of the Brillouin zone reveal a flat upper valence band for NiO, in good agreement with Angle Resolved Photoemission Spectroscopy results. To study defect properties, we evaluated the formation energies of the neutral and charged vacancies of oxygen and nickel in NiO. A formation energy of 7.2 (0.15) eV was found for the oxygen vacancy under oxygen rich conditions. For the Ni vacancy, we obtained a formation energy of 3.2 (0.15) eV under Ni rich conditions. These results confirm that NiO occurs as a p-type material with the dominant intrinsic vacancy defect being Ni vacancy.

We have studied the adsorption and diffusion of sulfur at the low-coverage regime of 0.25 ML on the (111), (100), (110), and (211) surfaces of Ni, Cu, Rh, Pd, Ag, Ir, Pt, and Au using density functional theory calculations. Sulfur adsorbed preferentially on three-fold or four-fold high-coordination sites over most of the studied surfaces. On the Ir(110), Pt(110), and Au(110) surfaces, sulfur is more stable on the two-fold sites. Calculations of the minimum energy diffusion pathway show that the energy barrier for the surface diffusion of sulfur depends on the orientation and nature of the metal surfaces. On the (100), sulfur shows the highest diffusion energy, ranging from 0.47 eV in Au(100) to 1.22 eV in Pd(100). In the (110) surface, the diffusion of sulfur is along the channel for Ni, Cu, Rh, Pd, and Ag, and across the channel for Ir, Pt, and Au. In the case of the (211) surfaces, the diffusion is preferentially along the terrace or step-edge sites. Our work provides data for the adsorption of sulfur on many surfaces not previously reported. The present work is a reference point for future computational studies of sulfur and sulfur-containing molecules absorbed on face center cubic metal surfaces.

We have applied the diffusion quantum Monte Carlo (DMC) method to calculate the cohesive energy and the structural parameters of the binary oxides CaO, SrO, BaO, Sc₂O₃, Y₂O₃, and La₂O₃. The aim of our calculations is to systematically quantify the accuracy of the DMC method to study this type of metal oxides. The DMC results were compared with local, semi-local, and hybrid Density Functional Theory (DFT) approximations as well as with experimental measurements. The DMC method yields cohesive energies for these oxides with a mean absolute deviation from experimental measurements of 0.18(2) eV, while with local, semi-local, and hybrid DFT approximations, the deviation is 3.06, 0.94, and 1.23 eV, respectively. For lattice constants, the mean absolute deviations in DMC, local, semi-local, and hybrid DFT approximations are 0.017(1), 0.07, 0.05, and 0.04 Å, respectively. DMC is a highly accurate method, outperforming the DFT approximations in describing the cohesive energies and structural parameters of these binary oxides.