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.

This paper studies the problem of expected loss minimization given a data distribution that is dependent on the decision-maker's action and evolves dynamically in time according to a geometric decay process. Novel algorithms for both the information setting in which the decision-maker has a first order gradient oracle and the setting in which they have simply a loss function oracle are introduced. The algorithms operate on the same underlying principle: the decision-maker repeatedly deploys a fixed decision over the length of an epoch, thereby allowing the dynamically changing environment to sufficiently mix before updating the decision. The iteration complexity in each of the settings is shown to match existing rates for first and zero order stochastic gradient methods up to logarithmic factors. The algorithms are evaluated on a "semi-synthetic" example using real world data from the SFpark dynamic pricing pilot study; it is shown that the announced prices result in an improvement for the institution's objective (target occupancy), while achieving an overall reduction in parking rates., Comment: Accepted at AAAI 2022

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 study the problem of online resource allocation, where multiple customers arrive sequentially and the seller must irrevocably allocate resources to each incoming customer while also facing a procurement cost for the total allocation. Assuming resource procurement follows an a priori known marginally increasing cost function, the objective is to maximize the reward obtained from fulfilling the customers' requests sans the cumulative procurement cost. We analyze the competitive ratio of a primal-dual algorithm in this setting, and develop an optimization framework for synthesizing a surrogate function for the procurement cost function to be used by the algorithm, in order to improve the competitive ratio of the primal-dual algorithm. Our first design method focuses on polynomial procurement cost functions and uses the optimal surrogate function to provide a more refined bound than the state of the art. Our second design method uses quasiconvex optimization to find optimal design parameters for a general class of procurement cost functions. Numerical examples are used to illustrate the design techniques. We conclude by extending the analysis to devise a posted pricing mechanism in which the algorithm does not require the customers' preferences to be revealed.

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

A thin-film model for a meniscus driven by Rayleigh surface acoustic waves (SAW) is analysed, a problem closely related to the classical Landau-Levich or dragged-film problem where a plate is withdrawn at constant speed from a bath. We consider a mesoscopic hydrodynamic model for a partially wetting liquid, were wettability is incorporated via a Derjaguin (or disjoining) pressure and combine SAW driving with the elements known from the dragged-film problem. For a one-dimensional substrate, i.e., neglecting transverse perturbations, we employ numerical path continuation to investigate in detail how the various occurring steady and time-periodic states depend on relevant control parameters like the Weber number and SAW strength. The bifurcation structure related to qualitative transitions caused by the SAW is analysed with particular attention on the {appearance and interplay of Hopf bifurcations where branches of time-periodic states emerge. The latter correspond to the regular shedding of liquid ridges from the meniscus. The obtained information is relevant to the entire class of dragged-film problems.

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 a brief comparative investigation of the bifurcation structure related to the formation of two-dimensional deposition patterns as described by continuum models of Cahn-Hilliard type. These are, on the one hand a driven Cahn-Hilliard model for Langmuir-Blodgett transfer of a surfactant layer from the surface of a bath onto a moving plate and on the other hand a driven thin-film equation modelling the surface acoustic wave-driven coating of a plate by a simple liquid. In both cases, we present selected two-dimensional steady states corresponding to deposition patterns and discuss the main structure of the corresponding bifurcation diagrams.

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

We introduce algorithms for online, full-information prediction that are competitive with contextual tree experts of unknown complexity, in both probabilistic and adversarial settings. We show that by incorporating a probabilistic framework of structural risk minimization into existing adaptive algorithms, we can robustly learn not only the presence of stochastic structure when it exists (leading to constant as opposed to $\mathcal{O}(\sqrt{T})$ regret), but also the correct model order. We thus obtain regret bounds that are competitive with the regret of an optimal algorithm that possesses strong side information about both the complexity of the optimal contextual tree expert and whether the process generating the data is stochastic or adversarial. These are the first constructive guarantees on simultaneous adaptivity to the model and the presence of stochasticity., Comment: 33 pages, 5 figures

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

We present calculations of electronic and magnetic structures of vanadium-benzene multidecker clusters V$_{n}$Bz$_{n+1}$ ($n$ = 1 - 3) using advanced quantum Monte Carlo methods. These and related systems have been identified as prospective spin filters in spintronic applications, assuming that their ground states are half-metallic ferromagnets. Although we find that magnetic properties of these multideckers are consistent with ferromagnetic coupling, their electronic structures do not appear to be half-metallic as previously assumed. In fact, they are ferromagnetic insulators with large and broadly similar $\uparrow$-/$\downarrow$-spin gaps. This makes the potential of these and related materials as spin filtering devices very limited, unless they are further modified or functionalized.