Search

We consider non-relativistic electron correlation energies of heavy noble gas atoms including the superheavy element Og. The corresponding data enables us to quantify fixed-node errors in real space quantum Monte Carlo methods as a function of the atomic number $Z$. We confirm that single-reference trial function nodes lead to an overall trend of mild decrease in recovered correlation energy with the increasing $Z$. This agrees with our previous study that has shown increasing fixed-node biases with the increasing electron density. We also estimate the value of the linear term in the asymptotic expansion of the atomic correlation energy., Comment: Submitted in May 2024

The field of two-dimensional (2D) materials has grown dramatically in the last two decades. 2D materials can be utilized for a variety of next-generation optoelectronic, spintronic, clean energy, and quantum computation applications. These 2D structures, which are often exfoliated from layered van der Waals (vdW) materials, possess highly inhomogeneous electron densities and can possess short- and long-range electron correlations. The complexities of 2D materials make them challenging to study with standard mean-field electronic structure methods such as density functional theory (DFT), which relies on approximations for the unknown exchange-correlation functional. In order to overcome the limitations of DFT, highly accurate many-body electronic structure approaches such as Diffusion Monte Carlo (DMC) can be utilized. In the past decade, DMC has been used to calculate accurate magnetic, electronic, excitonic, and topological properties in addition to accurately capturing interlayer interactions and cohesion and adsorption energetics of 2D materials. This approach has been applied to 2D systems of wide interest including graphene, phosphorene, MoS$_2$, CrI$_3$, VSe$_2$, GaSe, GeSe, borophene, and several others. In this review article, we highlight some successful recent applications of DMC to 2D systems for improved property predictions beyond standard DFT.

We derive many-body single- and multi-reference wave functions for quantum Monte Carlo of periodic systems with an anti-symmetric portion that explicitly integrates over the Brillouin zone of one-particle Bloch states. The wave functions are BCS-like determinants for singlets and pfaffians for polarized states built with appropriate pair orbitals. This ab initio formalism is broadly applicable, eg, to description of quasi-particle band gaps, optical excitations and to systems with complicated Fermi surfaces. It generalizes to spin-dependent interactions using two-component spinor pairs.

We construct correlation-consistent effective core potentials (ccECPs) for a selected set of heavy atoms and f-elements that are of significant current interest in materials and chemical applications, including Y, Zr, Nb, Rh, Ta, Re, Pt, Gd, and Tb. As customary, ccECPs consist of spin-orbit averaged relativistic effective potential (AREP) and effective spin-orbit (SO) terms. For the AREP part, our constructions are carried out within a relativistic coupled-cluster framework while also taking into objective function one-particle characteristics for improved convergence in optimizations. The transferability is adjusted using binding curves of hydride and oxide molecules. We address the difficulties encountered with f-elements, such as the presence of large cores and multiple near-degeneracies of excited levels. For these elements, we construct ccECPs with core-valence partitioning that includes 4f-subshell in the valence space. The developed ccECPs achieve an excellent balance between accuracy, size of the valence space, and transferability and are also suitable to be used in plane wave codes with reasonable energy cutoffs.

A member of the RMn$_6$Sn$_6$ rare-earth family materials, TbMn$_6$Sn$_6$, recently showed experimental signatures of the realization of a quantum-limit Chern magnet. In this work, we use quantum Monte Carlo (QMC) and density functional theory with Hubbard $U$ (DFT$+U$) calculations to examine the electronic structure of TbMn$_6$Sn$_6$. To do so, we optimize accurate, correlation-consistent pseudopotentials for Tb and Sn using coupled-cluster and configuration-interaction (CI) methods. We find that DFT$+U$ and single-reference QMC calculations suffer from the same overestimation of the magnetic moments as meta-GGA and hybrid density functional approximations. Our findings point to the need for improved orbitals/wavefunctions for this class of materials, such as natural orbitals from CI, or for the inclusion of multi-reference effects that capture the static correlations for an accurate prediction of magnetic properties. DFT$+U$ with Mn magnetic moments adjusted to experiment predict the Dirac crossing in bulk to be close to the Fermi level, within $\sim 120$ meV, in agreement with the experiments. Our non-stoichiometric slab calculations show that the Dirac crossing approaches even closer to the Fermi level, suggesting the possible realization of Chern magnetism in this limit., Comment: This is the peer-reviewed, accepted version of the manuscript. 17 pages, 10 figures. https://doi.org/10.1038/s41535-023-00583-6

We construct a new modification of correlation consistent effective potentials (ccECPs) for late $3d$ elements Cr-Zn with Ne-core that are adapted for efficiency and low energy cut-offs in plane wave calculations. The decrease in accuracy is rather minor so that the constructions are in the same overall accuracy class as the original ccECPs. The resulting new constructions work with energy cut-offs at or below $\approx$ 400 Ry and thus make calculations of large systems with transition metals feasible for plane wave codes. We provide also the basic benchmarks for atomic spectra and molecular tests of this modified option that we denote as ccECP-soft.

Accurately predicting the formation energy of a compound, which describes its thermodynamic stability, is a key challenge in materials physics. Here, we employ many-body quantum Monte Carlo (QMC) with single-reference trial functions to compute the formation energy of two electronically disparate compounds, the intermetallic VPt$_2$ and the semiconductor CuI, for which standard density functional theory (DFT) predictions using both the Perdew-Burke Ernzerhof (PBE) and the strongly constrained and appropriately normed (SCAN) density functional approximations deviate markedly from available experimental values. For VPt$_2$, we find an agreement between QMC, SCAN, and PBE0 estimates, which therefore remain in disagreement with the much less exothermic experimental value. For CuI, the QMC result agrees with neither SCAN nor PBE pointing towards DFT exchange-correlation biases, likely related to the localized Cu $3d$ electrons. Compared to the behavior of some density functional approximations within DFT, spin-averaged QMC exhibits a smaller but still appreciable deviation when compared to experiment. The QMC result is slightly improved by incorporating spin-orbit corrections for CuI and solid I$_2$, so that experiment and theory are brought into imperfect but reasonable agreement within about 120~meV/atom.

Layered material $\alpha$-RuCl$_3$ has caught wide attention due to its possible realization of Kitaev's spin liquid and its electronic structure that involves the interplay of electron-electron correlations and spin-orbit effects. Several DFT$+U$ studies have suggested that both electron-electron correlations and spin-orbit effects are crucial for accurately describing the band gap. This work studies the importance of these two effects using fixed-node and fixed-phase diffusion Monte Carlo calculations both in spin-averaged and explicit spin-orbit formalisms. In the latter, the Slater-Jastrow trial function is constructed from two-component spin-orbitals using our recent quantum Monte Carlo (QMC) developments and thoroughly tested effective core potentials. Our results show that the gap in the ideal crystal is already accurately described by the spin-averaged case, with the dominant role being played by the magnetic ground state with significant exchange and electron correlation effects. We find qualitative agreement between hybrid DFT, DFT+$U$, and QMC. In addition, QMC results agree very well with available experiments, and we identify the values of exact Fock exchange mixing that provide comparable gaps. Explicit spin-orbit QMC calculations reveal that the effect of spin-orbit coupling on the gap is minor, of the order of 0.2 eV, which corresponds to the strength of the spin-orbit of the Ru atom.

We introduce new correlation consistent effective core potentials (ccECPs) for the elements I, Te, Bi, Ag, Au, Pd, Ir, Mo, and W with $4d$, $5d$, $6s$ and $6p$ valence spaces. These ccECPs are given as a sum of spin-orbit averaged relativistic effective potential (AREP) and effective spin-orbit (SO) terms. The construction involves several steps with increasing refinements from more simple to fully correlated methods. The optimizations are carried out with objective functions that include weighted many-body atomic spectra, norm-conservation criteria, and spin-orbit splittings. Transferability tests involve molecular binding curves of corresponding hydride and oxide dimers. The constructed ccECPs are systematically better and in a few cases on par with previous effective core potential (ECP) tables on all tested criteria and provide a significant increase in accuracy for valence-only calculations with these elements. Our study confirms the importance of the AREP part in determining the overall quality of the ECP even in the presence of sizable spin-orbit effects. The subsequent quantum Monte Carlo (QMC) calculations point out the importance of accurate trial wave functions which in some cases (mid series transition elements) require treatment well beyond single-reference., Comment: 24 pages, 27 figures

We study beryllium dihydride (BeH$_2$) and acetylene (C$_2$H$_2$) molecules using real-space diffusion Monte Carlo (DMC) method. The molecules serve as perhaps the simplest prototypes that illustrate the difficulties with biases in the fixed-node DMC calculations that might appear with the use of effective core potentials (ECPs) or other nonlocal operators. This is especially relevant for the recently introduced correlation consistent ECPs (ccECPs) for $2s2p$ elements. Corresponding ccECPs exhibit deeper potential functions due to higher fidelity to all-electron counterparts, which could lead to larger local energy fluctuations. We point out that the difficulties stem from issues that are straightforward to address by upgrades of basis sets, use of T-moves for nonlocal terms, inclusion of a few configurations into the trial function and similar. The resulting accuracy corresponds to the ccECP target (chemical accuracy) and it is in consistent agreement with independent correlated calculations. Further possibilities for upgrading the reliability of the DMC algorithm and considerations for better adapted and more robust Jastrow factors are discussed as well., Comment: 13 pages, 5 figures

We study the nodal properties of many-body eigenstates of stationary Schr\"odinger equation that affect the accuracy of real-space quantum Monte Carlo calculations. In particular, we introduce weighted nodal domain averages that provide a new probe of nodal surfaces beyond the usual expectations. Particular choices for the weight function reveal, for example, that the difference between two arbitrary fermionic eigenvalues is given by the nodal hypersurface integrals normalized by overlaps with the bosonic ground state of the given Hamiltonian. Noninteracting and fully interacting Be atom with corresponding almost exact and approximate wave functions are used to illustrate several aspects of these concepts. Variational formulations that employ different weights are proposed for prospective improvement of nodes in variational and fixed-node diffusion Monte Carlo calculations.

We have carried out quantum Monte Carlo (QMC) calculations of silicon crystal focusing on the accuracy and systematic biases that affect the electronic structure characteristics. The results show that 64 and 216 atom supercells provide an excellent consistency for extrapolated energies per atom in the thermodynamic limit for ground, excited, and ionized states. We have calculated the ground state cohesion energy with both $\textit{systematic and statistical errors}$ below $\approx$0.05 eV. The ground state exhibits a fixed-node error of only $1.3(2)\%$ of the correlation energy, suggesting an unusually high accuracy of the corresponding single-reference trial wave function. We obtain a very good agreement between optical and quasiparticle gaps that affirms the marginal impact of excitonic effects. Our most accurate results for band gaps differ from the experiments by about 0.2 eV. This difference is assigned to a combination of residual finite-size and fixed-node errors. We have estimated the crystal Fermi level referenced to vacuum that enabled us to calculate the edges of valence and conduction bands in agreement with experiments., Comment: Peer-reviewed version

We present high-accuracy correlated calculations of small Si$_x$H$_y$ molecular systems both in the ground and excited states. We employ quantum Monte Carlo (QMC) together with a variety of many-body wave function approaches based on basis set expansions. The calculations are carried out in a valence-only framework using recently derived correlation consistent effective core potentials. Our primary goal is to understand the fixed-node diffusion QMC errors in both the ground and excited states with single-reference trial wave functions. Using a combination of methods, we demonstrate the very high accuracy of the QMC atomization energies being within $\approx$ 0.07 eV or better when compared with essentially exact results. By employing proper choices for trial wave functions, we have found that the fixed-node QMC biases for total energies are remarkably uniform ranging between $1-3.5$ % with absolute values at most $\approx$ 0.2 eV across the systems and several types of excitations such as singlets and triplets as well as low-lying and Rydberg-like states. Our results further corroborate that Si systems, and presumably also related main group IV and V elements of the periodic table (Ge, Sn, etc), exhibit some of the lowest fixed-node biases found in valence-only electronic structure QMC calculations., Comment: 10 pages, 1 figure

We review recent advances in the capabilities of the open source ab initio Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for greater efficiency and reproducibility. The auxiliary field QMC (AFQMC) implementation has been greatly expanded to include k-point symmetries, tensor-hypercontraction, and accelerated graphical processing unit (GPU) support. These scaling and memory reductions greatly increase the number of orbitals that can practically be included in AFQMC calculations, increasing accuracy. Advances in real space methods include techniques for accurate computation of band gaps and for systematically improving the nodal surface of ground state wavefunctions. Results of these calculations can be used to validate application of more approximate electronic structure methods including GW and density functional based techniques. To provide an improved foundation for these calculations we utilize a new set of correlation-consistent effective core potentials (pseudopotentials) that are more accurate than previous sets; these can also be applied in quantum-chemical and other many-body applications, not only QMC. These advances increase the efficiency, accuracy, and range of properties that can be studied in both molecules and materials with QMC and QMCPACK.

We present real space quantum Monte Carlo (QMC) calculations of the scandate LaScO$_3$ that proved to be challenging for traditional electronic structure approaches due to strong correlation effects resulting in inaccurate band gaps from DFT and $GW$ methods when compared with existing experimental data. Besides calculating an accurate QMC band gap corrected for supercell size biases and in agreement with numerous experiments, we also predict the cohesive energy of the crystal using the standard fixed-node QMC without any empirical or non-variational parameters. We show that promotion (optical) gap and fundamental gap agree with each other illustrating a clear absence of significant excitonic effects in the ideal crystal. We obtained these results in perfect consistency in two independent tracks that employ different basis sets (plane wave vs. localized gaussians), different codes for generating orbitals (\textsc{Quantum Espresso} vs. \textsc{Crystal}), different QMC codes (\textsc{Qmcpack} vs. \textsc{Qwalk}) and different high-accuracy pseudopotentials (ccECPs vs. Troullier-Martins) presenting the maturity and consistency of QMC methodology and tools for studies of strongly correlated problems., Comment: 12 pages, 8 tables, 5 figures

Very recently, we introduced a set of correlation consistent effective core potentials (ccECPs) constructed within full many-body approaches. By employing significantly more accurate correlated approaches we were able to reach a new level of accuracy for the resulting effective core Hamiltonians. We also strived for simplicity of use and easy transferability into a variety of electronic structure methods in quantum chemistry and condensed matter physics. Here, as a reference for future use, we present exact or nearly-exact total energy calculations for these ccECPs. The calculations cover H-Kr elements and are based on the state-of-the-art configuration interaction (CI), coupled-cluster (CC), and quantum Monte Carlo (QMC) calculations with systematically eliminated/improved errors. In particular, we carry out full CI/CCSD(T)/CCSDT(Q) calculations with cc-pVnZ with up to n=6 basis sets and we estimate the complete basis set limits. Using combinations of these approaches, we achieved an accuracy of $\approx$ 1-10 mHa for K-Zn atoms and $\approx$ 0.1-0.3 mHa for all other elements $-$ within about 1% or better of the ccECP total correlation energies. We also estimate the corresponding kinetic energies within the feasible limit of full CI calculations. In order to provide data for QMC calculations, we include fixed-node diffusion Monte Carlo energies for each element that give quantitative insights into the fixed-node biases for single-reference trial wave functions. The results offer a clear benchmark for future high accuracy calculations in a broad variety of correlated wave function methods such as CI and CC as well is in stochastic approaches such as real space sampling QMC.

Recently, we developed a new method for generating effective core potentials (ECPs) using valence energy isospectrality with explicitly correlated all-electron (AE) excitations and norm-conservation criteria. We apply this methodology to the 3$^{rd}$-row main group elements, creating new correlation consistent effective core potentials (ccECPs) and also derive additional ECPs to complete the ccECP table for H-Kr. For K and Ca, we develop Ne-core ECPs and for the $4p$ main group elements, we construct [Ar]$3d^{10}$-core potentials. Scalar relativistic effects are included in their construction. Our ccECPs reproduce AE spectra with significantly better accuracy than many existing pseudopotentials and show better overall consistency across multiple properties. The transferability of ccECPs is tested on monohydride and monoxide molecules over a range of molecular geometries. For the constructed ccECPs we also provide optimized DZ - 6Z valence Gaussian basis sets.

Quantum Monte Carlo (QMC) is applied to obtain the fundamental (quasiparticle) electronic band gap, $\Delta_f$, of a semiconducting two-dimensional (2D) phosphorene whose optical and electronic properties fill the void between graphene and 2D transition metal dichalcogenides. Similarly to other 2D materials, the electronic structure of phosphorene is strongly influenced by reduced screening, making it challenging to obtain reliable predictions by single-particle density functional methods. Advanced GW techniques, which include many-body effects as perturbative corrections, are hardly consistent with each other, predicting the band gap of phosphorene with a spread of almost 1 eV, from 1.6 to 2.4 eV. Our QMC results, from infinite periodic superlattices as well as from finite clusters, predict $\Delta_f$ to be about 2.4 eV, indicating that available GW results are systematically underestimating the gap. Using the recently uncovered universal scaling between the exciton binding energy and $\Delta_f$, we predict the optical gap of 1.75 eV that can be directly related to measurements even on encapsulated samples due to its robustness against dielectric environment. The QMC gaps are indeed consistent with recent experiments based on optical absorption and photoluminescence excitation spectroscopy. We also predict the cohesion of phosphorene to be only slightly smaller than that of the bulk crystal. Our investigations not only benchmark GW methods and experiments, but also open the field of 2D electronic structure to computationally intensive but highly predictive QMC methods which include many-body effects such as electronic correlations and van der Waals interactions explicitly., Comment: 7 pages, 3 figures main text; 6 pages, 4 figures supplemental text

Recently, we have introduced a new generation of effective core potentials (ECPs) designed for accurate correlated calculations but equally useful for a broad variety of approaches. The guiding principle has been the isospectrality of all-electron and ECP Hamiltonians for a subset of valence many-body states using correlated, nearly-exact calculations. Here we present such ECPs for the 3d transition series Sc to Zn with Ne-core, i.e, with semi-core 3s and 3p electrons in the valence space. Besides genuine many-body accuracy, the operators are simple, being represented by a few gaussians per symmetry channel with resulting potentials that are bounded everywhere. The transferability is checked on selected molecular systems over a range of geometries. The ECPs show a high overall accuracy with valence spectral discrepancies typically $\approx$ 0.01-0.02 eV or better. They also reproduce binding curves of hydride and oxide molecules typically within 0.02-0.03 eV deviations over the full non-dissociation range of interatomic distances., Comment: 21 pages, 13 figures, 13 tables

Very recently, we have introduced correlation consistent effective core potentials (ccECPs) derived from many-body approaches with the main target being its use in explicitly correlated methods but also in mainstream approaches. The ccECPs are based on reproducing excitation energies for a subset of valence states, i.e., achieving a near-isospectrality between the original and pseudo Hamiltonians. In addition, binding curves of dimer molecules have been used for refinement and overall improvement of transferability over a range of bond lengths. Here we apply similar ideas to the second row elements and study several aspects of the constructions in order to find the optimal (or nearly-optimal) solutions within the chosen ECP forms with $3s,3p$ valence space (Ne-core). New constructions exhibit accurate low-lying atomic excitations and equilibrium molecular bonds (on average within $\approx$ $0.03$ eV and $3$ m\AA), however, the errors for Al and Si oxide molecules at short bond lengths are notably larger for both ours and existing ECPs. Assuming this limitation, our ccECPs show a systematic balance between the criteria of atomic spectra accuracy and transferability for molecular bonds. In order to provide another option with much higher uniform accuracy, we also construct He-core ECPs for the whole row with typical discrepancies of $\approx$ 0.01 eV or smaller.

QMCPACK is an open source quantum Monte Carlo package for ab-initio electronic structure calculations. It supports calculations of metallic and insulating solids, molecules, atoms, and some model Hamiltonians. Implemented real space quantum Monte Carlo algorithms include variational, diffusion, and reptation Monte Carlo. QMCPACK uses Slater-Jastrow type trial wave functions in conjunction with a sophisticated optimizer capable of optimizing tens of thousands of parameters. The orbital space auxiliary field quantum Monte Carlo method is also implemented, enabling cross validation between different highly accurate methods. The code is specifically optimized for calculations with large numbers of electrons on the latest high performance computing architectures, including multicore central processing unit (CPU) and graphical processing unit (GPU) systems. We detail the program's capabilities, outline its structure, and give examples of its use in current research calculations. The package is available at http://www.qmcpack.org .

We present a recently developed projector quantum Monte Carlo method for calculations of electronic structure in systems with spin-orbit interactions. The method solves for many-body eigenstates in the presence of spin-orbit using the fixed-phase approximation.The trial wave function is built from two-component spinors and explicit Jastrow correlation factors while the core electrons are eliminated by relativistic effective core potentials with explicit spin-orbit terms. We apply this method to WO and W$_2$ molecules that enable us to build multi-reference wave functions and analyze in detail the impact of both electron correlations and of spin-orbit terms. These developments open new opportunities for calculations of systems with significant spin-orbit effects by many-body quantum Monte Carlo methods., Comment: 19 pages, 2 figures, 2 tables

We outline ideas on desired properties for a new generation of effective core potentials (ECPs) that will allow valence-only calculations to reach the full potential offered by recent advances in many-body wave function methods. The key improvements include consistent use of correlated methods throughout ECP constructions and improved transferability as required for an accurate description of molecular systems over a range of geometries. The guiding principle is the isospectrality of all-electron and ECP Hamiltonians for a subset of valence states. We illustrate these concepts on a few first- and second-row atoms (B, C, N, O, S) and we obtain higher accuracy in transferability than previous constructions while using a semi-local ECPs with a small number of parameters. In addition, the constructed ECPs enable many-body calculations of valence properties with higher (or same) accuracy than their all-electron counterparts with uncorrelated cores. This implies that the ECPs include also some of the impacts of core-core and core-valence correlations on valence properties. The results open further prospects for ECP improvements and refinements., Comment: 14 pages, 14 figures, 11 tables

We study several aspects of the recently introduced fixed-phase spin-orbit diffusion Monte Carlo (FPSODMC) method, in particular, its relation to the fixed-node method and its potential use as a general approach for electronic structure calculations. We illustrate constructions of spinor-based wave functions with the full space-spin symmetry without assigning up or down spin labels to particular electrons, effectively "complexifying" even ordinary real-valued wave functions. Interestingly, with proper choice of the simulation parameters and spin variables, such fixed-phase calculations enable one to reach also the fixed-node limit. The fixed-phase solution provides a straightforward interpretation as the lowest bosonic state in a given effective potential generated by the many-body approximate phase. In addition, the divergences present at real wave function nodes are smoothed out to lower dimensionality, decreasing thus the variation of sampled quantities and making the sampling also more straightforward. We illustrate some of these properties on calculations of selected first-row systems that recover the fixed-node results with quantitatively similar levels of the corresponding biases. At the same time, the fixed-phase approach opens new possibilities for more general trial wave functions with further opportunities for increasing accuracy in practical calculations., Comment: 8 pages, 4 figures, 2 tables

We present a study of mono(benzene)TM and bis(benzene)TM systems, where TM={Mo,W}. We calculate the binding energies by quantum Monte Carlo (QMC) approaches and compare the results with other methods and available experiments. The orbitals for the determinantal part of each trial wave function were generated from several types of DFT in order to optimize for fixed-node errors. We estimate and compare the size of the fixed-node errors for both the Mo and W systems with regard to the electron density and degree of localization in these systems. For the W systems we provide benchmarking results of the binding energies, given that experimental data is not available.

We compare the fixed-phase approximation with the better known, but closely related fixed-node approximation on several testing examples. We found that both approximations behave very similarly with the fixed-phase results being very close to the fixed-node method whenever nodes/phase were of high and comparable accuracy. The fixed-phase exhibited larger biases when the trial wave functions errors in the nodes/phase were intentionally driven to unrealistically large values. We also present a formalism that enables to describe wave functions with the full antisymmetry in spin-spatial degrees of freedom using our recently developed method for systems with spins as fully quantum variables. This opens new possibilities for simulations of fermionic systems in the fixed-phase approximation formalism., Comment: 6 pages, 2 figures

We investigate the inclusion of variable spins in electronic structure quantum Monte Carlo, with a focus on diffusion Monte Carlo with Hamiltonians that include spin-orbit interactions. Following our previous introduction of fixed-phase spin-orbit diffusion Monte Carlo (FPSODMC), we thoroughly discuss the details of the method and elaborate upon its technicalities. We present a proof for an upper-bound property for complex nonlocal operators, which allows for the implementation of T-moves to ensure the variational property. We discuss the time step biases associated with our particular choice of spin representation. Applications of the method are also presented for atomic and molecular systems. We calculate the binding energies and geometry of the PbH and Sn$_2$ molecules, as well as the electron affinities of the 6$p$ row elements in close agreement with experiments., Comment: 14 pages, 5 figures

We develop generalization of the fixed-phase diffusion Monte Carlo method for Hamiltonians which explicitly depend on particle spins such as for spin-orbit interactions. The method is formulated in zero variance manner and is similar to treatment of nonlocal operators in commonly used static- spin calculations. Tests on atomic and molecular systems show that it is very accurate, on par with the fixed-node method. This opens electronic structure quantum Monte Carlo methods to a vast research area of quantum phenomena in which spin-related interactions play an important role., Comment: Version 3: Some text additions. Results and conclusions unchanged. 5 pages, 2 figures

We study lithium systems over a range of number of atoms, e.g., atomic anion, dimer, metallic cluster, and body-centered cubic crystal by the diffusion Monte Carlo method. The calculations include both core and valence electrons in order to avoid any possible impact by pseudo potentials. The focus of the study is the fixed-node errors, and for that purpose we test several orbital sets in order to provide the most accurate nodal hyper surfaces. We compare our results to other high accuracy calculations wherever available and to experimental results so as to quantify the the fixed-node errors. The results for these Li systems show that fixed-node quantum Monte Carlo achieves remarkably high accuracy total energies and recovers 97-99 % of the correlation energy., Comment: 7 pages, 4 figures

Reliable theoretical predictions of noncovalent interaction energies, which are important e.g. in drug-design and hydrogen-storage applications, belong to longstanding challenges of contemporary quantum chemistry. In this respect, the fixed-node diffusion Monte Carlo (FN-DMC) is a promising alternative to the commonly used "gold standard" coupled-cluster CCSD(T)/CBS method for its benchmark accuracy and favourable scaling, in contrast to other correlated wave function approaches. This work is focused on the analysis of protocols and possible tradeoffs for FN-DMC estimations of noncovalent interaction energies and proposes an efficient yet accurate computational protocol using simplified explicit correlation terms with a favorable O(N^3) scaling. It achieves an excellent agreement (mean unsigned error ~0.2 kcal/mol) with respect to the CCSD(T)/CBS data on a number of complexes, including benzene/hydrogen,T-shape benzene dimer, stacked adenine-thymine and a set of small noncovalent complexes A24. The high accuracy and reduced computational costs predestinate the reported protocol for practical interaction energy calculations of large noncovalent complexes, where the CCSD(T)/CBS is prohibitively expensive.

We elucidate the origin of large differences (two-fold or more) in the fixed-node errors between the first- vs second-row systems for single-configuration trial wave functions in quantum Monte Carlo calculations. This significant difference in the fixed-node biases is studied across a set of atoms, molecules, and also Si, C solid crystals. The analysis is done over valence isoelectronic systems that share similar correlation energies, bond patterns, geometries, ground states, and symmetries. We show that the key features which affect the fixed-node errors are the differences in electron density and the degree of node nonlinearity. The findings reveal how the accuracy of the quantum Monte Carlo varies across a variety of systems, provide new perspectives on the origins of the fixed-node biases in electronic structure calculations of molecular and condensed systems, and carry implications for pseudopotential constructions for heavy elements

We outline the basic notions of nodal hypersurface and domain averages for antisymmetric wave functions. We illustrate their properties and analyze the results for a few electron explicitly solvable cases and discuss possible further developments.

Heteronuclear dimers are of significant interest to experiments seeking to exploit ultracold polar molecules in a number of novel ways including precision measurement, quantum computing, and quantum simulation. We calculate highly accurate Born-Oppenheimer total energies and electric dipole moments as a function of internuclear separation for two such dimers, LiSr and KRb. We apply fully-correlated, high-accuracy quantum Monte Carlo methods for evaluating these molecular properties in a many-body framework. We use small-core effective potentials combined with multi-reference Slater-Jastrow trial wave functions to provide accurate nodes for the fixed-node diffusion Monte Carlo method. For reference and comparison, we calculate the same properties with Hartree-Fock and with restricted Configuration Interaction methods, and carefully assess the impact of the recovered many-body correlations on the calculated quantities. For LiSr we find a highly nonlinear dipole moment curve, which may make this molecule's dipole moment tunable through vibrational state control., Comment: 13 pages, 6 figures, 2 tables, 78 references. Submitted to a special issue of Molecular Physics on "Manipulation of Molecules with Electromagnetic Fields." Published June 2013. http://dx.doi.org/10.1080/00268976.2013.788741

We analyze the effect of increasing charge density on the Fixed Node Errors in Diffusion Monte Carlo by comparing FN-DMC calculations of the total ground state energy on a 4 electron system done with a Hartree-Fock based trial wave function to calculations by the same method on the same system using a Configuration Interaction based trial wave function. We do this for several different values of nuclear charge, Z. The Fixed Node Error of a Hartree-Fock trial wave function for a 4 electron system increases linearly with increasing nuclear charge.

We calculate the ground-state properties of unpolarized two-component Fermi gas by the diffusion quantum Monte Carlo (DMC) methods. Using an extrapolation to the zero effective range of the attractive two-particle interaction, we find $E/E_{\rm free}$ to be 0.212(2), 0.407(2), 0.409(3) and 0.398(3) for 4, 14, 38 and 66 atoms, respectively. Our results indicate that the dependence of the total energy on the effective range is sizable and the extrapolation is therefore quite important. In order to test the quality of nodal surfaces and to estimate the impact of the fixed-node approximation we perform released-node DMC calculations for 4 and 14 atoms. Analysis of the released-node and the fixed-node results suggests that the main sources of the fixed-node errors are long-range correlations which are difficult to sample in the released-node approaches due to the fast growth of the bosonic noise. Besides energies, we evaluate the two-body density matrix and the condensate fraction. We find that the condensate fraction for the 66 atom system converges to 0.56(1) after the extrapolation to the zero interaction range., Comment: 7pages, 6 figures

The unitarity limit describes interacting particles where the range of the interaction is zero and the scattering length is infinite. We present precision benchmark calculations for two-component fermions at unitarity using three different ab initio methods: Hamiltonian lattice formalism using iterated eigenvector methods, Euclidean lattice formalism with auxiliary-field projection Monte Carlo, and continuum diffusion Monte Carlo with fixed and released nodes. We have calculated the ground state energy of the unpolarized four-particle system in a periodic cube as a dimensionless fraction of the ground state energy for the non-interacting system. We obtain values 0.211(2) and 0.210(2) using two different Hamiltonian lattice representations, 0.206(9) using Euclidean lattice, and an upper bound of 0.212(2) from fixed-node diffusion Monte Carlo. Released-node calculations starting from the fixed-node result yield a decrease of less than 0.002 over a propagation of 0.4/E_F in Euclidean time, where E_F is the Fermi energy. We find good agreement among all three ab initio methods., Comment: 23 pages, 7 figures, final version to appear in Phys. Rev. A

Recently, a diffusion Monte Carlo algorithm was applied to the study of spin dependent interactions in condensed matter. Following some of the ideas presented therein, and applied to a Hamiltonian containing a Rashba-like interaction, a general variational Monte Carlo approach is here introduced that treats in an efficient and very accurate way the spin degrees of freedom in atoms when spin orbit effects are included in the Hamiltonian describing the electronic structure. We illustrate the algorithm on the evaluation of the spin-orbit splittings of isolated carbon and lead atoms. In the case of the carbon atom, we investigate the differences between the inclusion of spin-orbit in its realistic and effective spherically symmetrized forms. The method exhibits a very good accuracy in describing the small energy splittings, opening the way for a systematic quantum Monte Carlo studies of spin-orbit effects in atomic systems., Comment: 7 pages, 0 figures

The quantum Monte Carlo methods represent a powerful and broadly applicable computational tool for finding very accurate solutions of the stationary Schroedinger equation for atoms, molecules, solids and a variety of model systems. The algorithms are intrinsically parallel and are able to take full advantage of the present-day high-performance computing systems. This review article concentrates on the fixed-node/fixed-phase diffusion Monte Carlo method with emphasis on its applications to electronic structure of solids and other extended many-particle systems., Comment: 34 pages

Quantum Monte Carlo (QMC) is an advanced simulation methodology for studies of manybody quantum systems. In this review, we focus on the electronic structure QMC, i.e., methods relevant for systems described by the electron-ion Hamiltonians. Some of the key QMC achievements include direct treatment of electron correlation, accuracy in predicting energy differences and favorable scaling in the system size. Calculations of atoms, molecules, clusters and solids have demonstrated QMC applicability to real systems with hundreds of electrons while providing 90-95% of the correlation energy and energy differences typically within a few percent of experiments. Advances in accuracy beyond these limits are hampered by the so-called fixed-node approximation which is used to circumvent the notorious fermion sign problem. Many-body nodes of fermion states and their properties have therefore become one of the important topics for further progress in predictive power and efficiency of QMC calculations. Some of our recent results on the wave function nodes and related nodal domain topologies will be briefly reviewed. This includes analysis of few-electron systems and descriptions of exact and approximate nodes using transformations and projections of the highly-dimensional nodal hypersurfaces into the 3D space. Studies of fermion nodes offer new insights into topological properties of eigenstates such as explicit demonstrations that generic fermionic ground states exhibit the minimal number of two nodal domains. Recently proposed trial wave functions based on Pfaffians with pairing orbitals are presented and their nodal properties are tested in calculations of first row atoms and molecules. Finally, backflow "dressed" coordinates are introduced as another possibility for capturing correlation effects and for decreasing the fixed-node bias., Comment: published in: Acta Physica Slovaca 59, No.2, 81-168 (2009) (88 pages) (Received 10 May 2009, accepted 14 May 2009), Keywords: Condensed Matter, Computational Methods, Electronic Structure, Quantum Monte Carlo, Correlated Electrons, Fermion Nodes, Pfaffians; PACS:02.70.Ss, 03.65.Ge, 05.30.Fk, 31.25.-v, 67.10.Db, 71.10.w, 71.15.-m, 71.20.-b

We investigate how the fixed-node diffusion Monte Carlo energy of solids depends on single-particle orbitals used in Slater--Jastrow wave functions. We demonstrate that the dependence can be significant, in particular in the case of 3d transition-metal compounds, which we adopt as examples. We illustrate how exchange-correlation functionals with variable exact-exchange component can be exploited to reduce the fixed-node errors. On the basis of these results we argue that the fixed-node quantum Monte Carlo provides a variational approach for optimization of effective hamiltonians with parameters., Comment: 8 pages, 6 figures, RevTeX4

We calculate the electronic and atomic structure of (MnO)n (n=1-4) using the HF exchange, VWN, PBE and B3LYP exchange-correlation functionals. We also perform diffusion Monte Carlo calculation to evaluate more accurate energies. We ompare these results and discuss the accuracy of the exchange-correlation functionals., Comment: 11 pages

We determine the equation of state of stoichiometric FeO employing the diffusion Monte Carlo method. The fermionic nodes are fixed to those of a wave function having the form of a single Slater determinant. The calculated ambient pressure properties (lattice constant, bulk modulus and cohesive energy) agree very well with available experimental data. At approximately 65 GPa, the lattice structure is found to change from rocksalt type (B1) to NiAs based (inverse B8)., Comment: 5 pages, 3 figures

The persistent current of correlated electrons in a continuous one-dimensional ring with a single scatterer is calculated by solving the many-body Schrodinger equation for several tens of electrons interacting via the electron-electron (e-e) interaction of finite range. The problem is solved by the configuration-interaction (CI) and diffusion Monte Carlo (DMC) methods. The CI and DMC results are in good agreement. In both cases, the persistent current $I$ as a function of the ring length $L$ exhibits the asymptotic dependence $I \propto L^{-1-\alpha}$ typical of the Luttinger liquid, where the power $\alpha$ depends only on the e-e interaction. The numerical values of $\alpha$ agree with the known formula of the renormalisation-group theory., Comment: Conference proceedings. Accepted for publication in Physica E

We describe QWalk, a new computational package capable of performing Quantum Monte Carlo electronic structure calculations for molecules and solids with many electrons. We describe the structure of the program and its implementation of Quantum Monte Carlo methods. It is open-source, licensed under the GPL, and available at the web site http://www.qwalk.org, Comment: 12 pages, 3 figures, describes code available at http://www.qwalk.org

The transition metal-oxygen bond appears very prominently throughout chemistry and solid-state physics. Many materials, from biomolecules to ferroelectrics to the components of supernova remnants contain this bond in some form. Many of these materials' properties depend strongly on fine details of the TM-O bond, which makes accurate calculations of their properties very challenging. Here we report on highly accurate first principles calculations of the properties of TM monoxide molecules within fixed-node Diffusion Monte Carlo and Reptation Monte Carlo., Comment: 6 pages, 2 figures, submitted to Journal of Chemical Physics

We employ the density functional theory to study a structural transition of MnO from B1 (rocksalt) to B8 (NiAs) structures that was observed experimentally at pressures around 100 GPa. We utilize all-electron description as well as norm-conserving pseudopotentials and demonstrate that these two approaches can significantly differ in quantitative predictions. We explicitly show that even small-core pseudopotentials exhibit transferability inaccuracies for quantities sensitive to the energy differences between high- and low-spin polarizations of valence electrons., Comment: 5 pages, 3 figures, REVTeX4

A fermion node is subset of fermionic configurations for which a real wave function vanishes due to the antisymmetry and the node divides the configurations space into compact nodal cells (domains). We analyze the properties of fermion nodes of fermionic ground state wave functions for a number of systems. For several models we demonstrate that noninteracting spin-polarized fermions in dimension two and higher have closed-shell ground state wave functions with the minimal two nodal cells for any system size and we formulate a theorem which sumarizes this result. The models include periodic fermion gas, fermions on the surface of a sphere, fermions in a box. We prove the same property for atomic states with up to 3d half-filled shells. Under rather general assumptions we then derive that the same is true for unpolarized systems with arbitrarily weak interactions using Bardeen-Cooper-Schrieffer (BCS) variational wave function. We further show that pair correlations included in the BCS wave function enable singlet pairs of particles to wind around the periodic box without crossing the node pointing towards the relationship of nodes to transport and many-body phases such as superconductivity. Finally, we point out that the arguments extend also to fermionic temperature dependent/imaginary-time density matrices. The results reveal fundamental properties of fermion nodal structures and provide new insights for accurate constructions of wave functions and density matrices in quantum and path integral Monte Carlo methods., Comment: 17 pages, 10 figures

We study nodes of fermionic ground state wave functions. For 2D and higher we prove that spin-polarized, noninteracting fermions in a harmonic well have two nodal cells for arbitrary system size. The result extends to other noninteracting/mean-field models such as fermions on a sphere, in a periodic box or in Hartree-Fock atomic states. Spin-unpolarized noninteracting states have multiple nodal cells, however, interactions and many-body correlations generally relax the multiple cells to the minimal number of two. With some conditions, this is proved for interacting 2D and higher dimensions harmonic fermion systems of arbitrary size using the Bardeen-Cooper-Schrieffer variational wave function. Implications and extent of these results are briefly discussed., Comment: 4 pages, 2 figures

Electronic structure calculations using quantum Monte Carlo (QMC) methods conclusively show that pure CaB$_6$ is a narrow-gap semiconductor with an $X$-point gap $\sim 1.3$ eV. This should put to rest controversies as to whether the compound is a semimetal or a semiconductor, as also theories suggesting it to be an excitonic insulator. Our extensive Hartree-Fock (HF) and density functional theory (DFT) calculations support the view that the hexaboride ferromagnetism can be induced by contamination with magnetic element atoms like Fe and that presence of La is not essential for this. La impurity, however, gives rise to a metallic state with a small electron-like Fermi surface as seen in ARPES experiments., Comment: 5 pages, 3 figures

Electronic structure of the manganese oxide solid is studied by the quantum Monte Carlo (QMC) methods. The trial wavefunctions are built using orbitals from unrestricted Hartree-Fock and Density Functional Theory, and the electron-electron correlation is recovered by the fixed-node QMC. The correlation effects are significant and QMC estimations of the gap and cohesion show a very good agreement with experiment. Comparison with hybrid functional results points out the importance of the exact exchange for improvement of the Density Functional description of transition metal oxide systems., Comment: 4 pages, 2 figures