Point defect formation energies in graphene from diffusion quantum Monte Carlo and density functional theory
arXiv:2112.11539 · doi:10.1103/PhysRevB.105.184114
Abstract
Density functional theory (DFT) is widely used to study defects in monolayer graphene with a view to applications ranging from water filtration to electronics to investigation of radiation damage in graphite moderators. To assess the accuracy of DFT in such applications, we report diffusion quantum Monte Carlo (DMC) calculations of the formation energies of some common and important point defects in monolayer graphene: monovacancies, Stone-Wales defects, and silicon substitutions. We find that standard DFT methods underestimate monovacancy formation energies by around 1 eV. The disagreement between DFT and DMC is somewhat smaller for Stone-Wales defects and silicon substitutions. We examine vibrational contributions to the free energies of formation for these defects, finding that vibrational effects are non-negligible. Finally, we compare the DMC atomization energies of monolayer graphene, monolayer silicene, and bulk silicon, finding that bulk silicon is significantly more stable than monolayer silicene by 0.7522(5) eV per atom.
References in corpus (11)
- Electric Field Effect in Atomically Thin Carbon Films
- The electronic properties of graphene
- Chemical functionalization of graphene with defects
- Jastrow correlation factor for atoms, molecules, and solids
- Spin-Valleytronics in Silicene: Quantum-Spin-Quantum-Anomalous Hall Insulators and Single-Valley Semimetals
- Inhomogeneous backflow transformations in quantum Monte Carlo calculations
- The Finite Size Error in Many-body Simulations with long-Ranged Interactions
- Smooth relativistic Hartree-Fock pseudopotentials for H to Ba and Lu to Hg
- Norm-conserving Hartree-Fock pseudopotentials and their asymptotic behavior
- A comparative study of density functional and density functional tight binding calculations of defects in graphene
- Trail-Needs pseudopotentials in quantum Monte Carlo calculations with plane-wave/blip basis sets