An efficient hybrid orbital representation for quantum Monte Carlo calculations

Enhanced Twist-Averaging Technique for Magnetic Metals: Applications Using Quantum Monte Carlo

We propose an improved twist-averaging (TA) scheme for quantum Monte Carlo methods that use converged Kohn-Sham or Hartree-Fock orbitals as the reference. This TA technique is tailored to sample the Brillouin zone of magnetic metals, although it naturally extends to nonmagnetic (NM) conducting systems. The proposed scheme aims to reproduce the reference magnetization and achieves charge neutrality by construction, thus avoiding the large energy fluctuations and the postprocessing needed to correct the energies. It shows the most robust convergence of total energy and magnetism to the thermodynamic limit (TDL) when compared to four other TA schemes. Diffusion Monte Carlo applications are shown on NM Al and ferromagnetic α-Fe. The cohesive energy of Al in the TDL shows an excellent agreement with the experimental result. Furthermore, the magnetic moments in α-Fe exhibit rapid convergence with an increasing number of twists.

DFT+U and quantum Monte Carlo study of electronic and optical properties of AgNiO2 and AgNi1-xCoxO2 delafossite

As the only semimetallic d10-based delafossite, AgNiO2 has received a great deal of attention due to both its unique semimetallicity and its antiferromagnetism in the NiO2 layer that is coupled with a lattice distortion. In contrast, other delafossites such as AgCoO2 are insulating. Here we study how the electronic structure of AgNi1-xCoxO2 alloys vary with Ni/Co concentration, in order to investigate the electronic properties and phase stability of the intermetallics. While the electronic and magnetic structure of delafossites have been studied using density functional theory (DFT), earlier studies have not included corrections for strong on-site Coulomb interactions. In order to treat these interactions accurately, in this study we use Quantum Monte Carlo (QMC) simulations to obtain accurate estimates for the electronic and magnetic properties of AgNiO2. By comparison to DFT results we show that these electron correlations are critical to account for. We show that Co doping on the magnetic Ni sites results in a metal-insulator transition near x ∼0.33, and reentrant behavior near x ∼ 0.66.

Finite-Size Error Cancellation in Diffusion Monte Carlo Calculations of Surface Chemistry

The accurate prediction of reaction mechanisms in heterogeneous (surface) catalysis is one of the central challenges in computational chemistry. Quantum Monte Carlo methods─Diffusion Monte Carlo (DMC) in particular─are being recognized as higher-accuracy, albeit more computationally expensive, alternatives to Density Functional Theory (DFT) for energy predictions of catalytic systems. A major computational bottleneck in the broader adoption of DMC for catalysis is the need to perform finite-size extrapolations by simulating increasingly large periodic cells (supercells) to eliminate many-body finite-size effects and obtain energies in the thermodynamic limit. Here, we show that it is possible to significantly reduce this computational cost by leveraging the cancellation of many-body finite-size errors that accompanies the evaluation of energy differences when calculating quantities like adsorption (binding) energies and mapping potential energy surfaces. We analyze the cancellation and convergence of many-body finite-size errors in two well-known adsorbate/slab systems, H₂O/LiH(001) and CO/Pt(111). Based on this analysis, we identify strategies for obtaining binding energies in the thermodynamic limit that optimally utilize error cancellation to balance accuracy and computational efficiency. Using one such strategy, we then predict the correct order of adsorption site preference on CO/Pt(111), a challenging problem for a wide range of density functionals. Our accurate and inexpensive DMC calculations are found to unambiguously recover the top > bridge > hollow site order, in agreement with experimental observations. We proceed to use this DMC method to map the potential energy surface of CO hopping between Pt(111) adsorption sites. This reveals the existence of an L-shaped top-bridge-hollow diffusion trajectory characterized by energy barriers that provide an additional kinetic justification for experimental observations of CO/Pt(111) adsorption. Overall, this work demonstrates that it is routinely possible to achieve order-of-magnitude speedups and memory savings in DMC calculations by taking advantage of error cancellation in the calculation of energy differences that are ubiquitous in heterogeneous catalysis and surface chemistry more broadly.