Diffusion Quantum Monte Carlo Benchmark of Interlayer Binding and Charge Redistribution in Chemically Distinct Two-Dimensional Van Der Waals Bilayers
Interlayer interactions in two-dimensional materials can generate emergent phenomena absent in their constituent monolayers, including unconventional magnetic order, multiferroicity, and topological magnetic phases. Predicting such emergent behavior requires simultaneously resolving long-range dispersion, short-range orbital hybridization, and electronic correlation interactions that are intrinsically nonlocal and many-body and remain challenging even for advanced density-functional approximations. Here we establish a systematically controlled many-body benchmark for diverse bilayer materials using diffusion Monte Carlo (DMC), spanning single-sheet materials, transition-metal dichalcogenides, and magnetic transition-metal halides. We obtain equilibrium separations, binding energetics, interlayer vibrational modes, and charge redistribution, finding excellent agreement with available experiments while revealing systematic and property-dependent failures across widely used semilocal, meta-GGA, and dispersion-corrected density functionals. Beyond energetics, DMC resolves subtle interlayer charge rearrangements, particularly leading to long-ranged dipolar tails in magnetic Cr trihalides that survive well beyond the regime where semilocal DFT predicts appreciable interlayer polarization, providing a many-body basis for understanding interlayer-coupled ferroic and magnetic phenomena. The resulting energies, response properties, and high-accuracy electron densities constitute a transferable benchmark for developing next-generation density functionals. Finally, we provide a scalable high-performance-computing workflow that enables systematic expansion of many-body benchmark datasets across emerging families of layered quantum materials by the scientific community.