Quantum Monte Carlo description of correlated electrons in two-dimensional FeSe

Electronic correlation effects on the structural properties of double-layer FeSe are studied using variation and diffusion quantum Monte Carlo methods. The Slater-Jastrow many-body wavefunction with two different forms for the homogeneous two-body pair-correlation term is used. The ground-state energy of the system is obtained at the thermodynamic limit using two different trial wave functions called JDFT and JSD. Only the Jastrow factor is fully optimized in the JDFT wave function, while the Slater determinant comes from the density functional approximation. In the JSD trial wave function, the Slater determinant and the Jastrow factor are fully optimized simultaneously. We calculated the VMC and DMC energies as a function of interlayer separation for two different in-plane iron-iron bond lengths. Our QMC results indicate that the optimized interlayer separation decreases with increasing iron-iron bond length, driven by stretch. We found that a three- to two-dimensional phase transition increases the electron-electron correlation effects and shifts the system from moderately correlated to strongly correlated. The value of correlation energy implies that the Hubbard value, which is widely used in DFT+U calculations, for two-dimensional FeSe is larger than its value in bulk.