Diffusion quantum Monte Carlo and density functional calculations of the structural stability of bilayer arsenene

We have studied the structural stability of monolayer and bilayer arsenene (As) in the buckled (b) and washboard (w) phases with diffusion quantum Monte Carlo (DMC) and density functional theory (DFT) calculations. DMC yields cohesive energies of 2.826(2) eV/atom for monolayer b-As and 2.792(3) eV/atom for w-As. In the case of bilayer As, DMC and DFT predict that AA-stacking is the more stable form of b-As, while AB is the most stable form of w-As. The DMC layer-layer binding energies for b-As-AA and w-As-AB are 30(1) and 53(1) meV/atom, respectively. The interlayer separations were estimated with DMC at 3.521(1) Å for b-As-AA and 3.145(1) Å for w-As-AB. A comparison of DMC and DFT results shows that the van der Waals density functional method yields energetic properties of arsenene close to DMC, while the DFT + D3 method closely reproduced the geometric properties from DMC. The electronic properties of monolayer and bilayer arsenene were explored with various DFT methods. The bandgap values vary significantly with the DFT method, but the results are generally qualitatively consistent. We expect the present work to be useful for future experiments attempting to prepare multilayer arsenene and for further development of DFT methods for weakly bonded systems.