BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Fernando A. Reboredo

Fernando A. Reboredo appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.mtrl-scicond-mat.str-elphysics.comp-phcond-mat.othercond-mat.supr-conphysics.chem-phcond-mat.mes-hallcond-mat.quant-gascond-mat.soft

Catalog footprint

What is connected

22works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2021arXiv

Candidate structure for the H$_2$-PRE phase of solid hydrogen

Experimental progress finally reached the metallic solid hydrogen phase, which was predicted by Wigner and Huntington over 80 years ago. However, the different structures in the phase diagram are still been debated due to the difficulty of diffraction experiments for high-pressured hydrogen. The determination of crystal structures under extreme condition is both of the basic condensed matter physics, and in planetary science: the behavior of giant gaseous planets (e.g. Jupiter, Saturn...) strongly depends on the properties of inner high-pressured hydrogen. This work describes new possible structures appearing under high pressures of 400$\sim$600 GPa. We applied a structural search using particle swarm optimization with density functional theory (DFT) to propose several candidate structures. For these structures, we performed fixed-node diffusion Monte Carlo simulations combined with DFT zero-point energy corrections to confirm their relative stability. We found $P2_{1}/c$-8 as a promising candidate structure for the H$_2$-PRE phase. $P2_{1}/c$-8 is predicted the most stable at 400 and 500~GPa. $P2_{1}/c$-8 reproduces qualitatively the IR spectrum peaks observed in the H$_2$-PRE phase.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Diffusion Monte Carlo Study on Relative Stabilities of Boron Nitride Polymorphs

Although Boron nitride (BN) is a well-known compound widely used for engineering and scientific purposes, the phase stability of its polymorphs, one of its most fundamental properties, is still under debate. The ab initio determination of the ground state of the BN polymorphs, such as hexagonal and zinc-blende, is difficult because of the elusive Van der Waals interaction, which plays a decisive role in some of the polymorphs, making quantitative prediction highly challenging. Hence, despite multiple theoretical studies, there has been no consensus on the ground state yet, primarily due to contradicting reports. In this study, we apply a state-of-the-art ab initio framework - fixed-node diffusion Monte Carlo (FNDMC), to four well known BN polymorphs, namely hexagonal, rhombohedral, wurtzite, and zinc-blende BNs. Our FNDMC calculations show that hBN is thermodynamically the most stable among the four polymorphs at 0 K as well as at 300K. This result agrees with the experimental data of Corrigan~{\it et al.} and Fukunaga. The conclusions are consistent with those obtained using other high-level methods, such as coupled cluster. We demonstrate that the FNDMC is a powerful method to address polymorphs that exhibit bonds of various forms. It also provides valuable information, like reliable reference energies, when reliable experimental data are missing or difficult to access. Our findings should promote the application of FNDMC for other van der Waals materials.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Range of biquadratic and triquadratic Heisenberg effective couplings deduced from multiorbital Hubbard models

We studied a multi-orbital Hubbard model at half-filling for two and three orbitals per site on a two-site cluster via full exact diagonalization, in a wide range for the onsite repulsion $U$, from weak to strong coupling, and multiple ratios of the Hund coupling $J_H$ to $U$. The hopping matrix elements among the orbitals were also varied extensively. At intermediate and large $U$, we mapped the results into a Heisenberg model. For two orbitals per site, the mapping is into a $S=1$ Heisenberg model where by symmetry both nearest-neighbor $(\mathbf{S}_{i}\cdot\mathbf{S}_{j})$ and $(\mathbf{S}_{i}\cdot\mathbf{S}_{j})^{2}$ are allowed, with respective couplings $J_1$ and $J_2$. For the case of three orbitals per site, the mappping is into a $S=3/2$ Heisenberg model with $(\mathbf{S}_{i}\cdot\mathbf{S}_{j})$, $(\mathbf{S}_{i}\cdot\mathbf{S}_{j})^{2}$, and $(\mathbf{S}_{i}\cdot\mathbf{S}_{j})^{3}$ terms, and respective couplings $J_1$, $J_2$, and $J_3$. The strength of these coupling constants in the Heisenberg models depend on the $U$, $J_H$, and hopping amplitudes of the underlying Hubbard model. Our study allows to establish bounds on how large the ratios $J_2/J_1$ and $J_3/J_1$ can be. We show that those ratios are severely limited and, as a crude guidance, we conclude that $J_2/J_1$ is less than 0.4 and $J_3/J_1$ is less than 0.2, establishing bounds on effective models for strongly correlated Hubbard systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion Quantum Monte Carlo

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cohesive energy and structural parameters of binary oxides of groups IIA and IIIB from diffusion quantum Monte Carlo

We have applied the diffusion quantum Monte Carlo (DMC) method to calculate the cohesive energy and the structural parameters of the binary oxides CaO, SrO, BaO, Sc2O3, Y2O3 and La2O3. The aim of our calculations is to systematically quantify the accuracy of the DMC method to study this type of metal oxides. The DMC results were compared with local, semi-local and hybrid Density Functional Theory (DFT) approximations as well as with experimental measurements. The DMC method yields cohesive energies for these oxides with a mean absolute deviation from experimental measurements of 0.18(2) eV, while with local, semi-local and hybrid DFT approximations the deviation is 3.06, 0.94 and 1.23 eV, respectively. For lattice constants, the mean absolute deviation in DMC, local, semi-local and hybrid DFT approximations, are 0.017(1), 0.07, 0.05 and 0.04 Å, respectively. DMC is highly accurate method, outperforming the DFT approximations in describing the cohesive energies and structural parameters of these binary oxides.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Strain Control of Oxygen Vacancies in Epitaxial Strontium Cobaltite Films

The ability to manipulate oxygen anion defects rather than metal cations in complex oxides can facilitate creating new functionalities critical for emerging energy and device technologies. However, the difficulty in activating oxygen at reduced temperatures hinders the deliberate control of important defects, oxygen vacancies. Here, strontium cobaltite (SrCoOx) is used to demonstrate that epitaxial strain is a powerful tool for manipulating the oxygen vacancy concentration even under highly oxidizing environments and at annealing temperatures as low as 300 C. By applying a small biaxial tensile strain (2%), the oxygen activation energy barrier decreases by ~30%, resulting in a tunable oxygen deficient steady-state under conditions that would normally fully oxidize unstrained cobaltite. These strain-induced changes in oxygen stoichiometry drive the cobaltite from a ferromagnetic metal towards an antiferromagnetic insulator. The ability to decouple the oxygen vacancy concentration from its typical dependence on the operational environment is useful for effectively designing oxides materials with a specific oxygen stoichiometry.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Structural Stability and Defect Energetics of ZnO from Diffusion Quantum Monte Carlo

We have applied the many-body ab-initio diffusion quantum Monte Carlo (DMC) method to study Zn and ZnO crystals under pressure, and the energetics of the oxygen vacancy, zinc interstitial and hydrogen impurities in ZnO. We show that DMC is an accurate and practical method that can be used to characterize multiple properties of materials that are challenging for density functional theory approximations. DMC agrees with experimental measurements to within 0.3 eV, including the band-gap of ZnO, the ionization potential of O and Zn, and the atomization energy of O$_2$, ZnO dimer, and wurtzite ZnO. DMC predicts the oxygen vacancy as a deep donor with a formation energy of 5.0(2) eV under O-rich conditions and thermodynamic transition levels located between 1.8 and 2.5 eV from the valence band maximum. Our DMC results indicate that the concentration of zinc interstitial and hydrogen impurities in ZnO should be low under n-type, and Zn- and H-rich conditions because these defects have formation energies above 1.4 eV under these conditions. Comparison of DMC and hybrid functionals shows that these DFT approximations can be parameterized to yield a general correct qualitative description of ZnO. However, the formation energy of defects in ZnO evaluated with DMC and hybrid functionals can differ by more than 0.5 eV.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Ab initio quantum Monte Carlo calculations of spin superexchange in cuprates: the benchmarking case of Ca$_2$CuO$_3$

In view of the continuous theoretical efforts aimed at an accurate microscopic description of the strongly correlated transition metal oxides and related materials, we show that with continuum quantum Monte Carlo (QMC) calculations it is possible to obtain the value of the spin superexchange coupling constant of a copper oxide in a quantitatively excellent agreement with experiment. The variational nature of the QMC total energy allows us to identify the best trial wave function out of the available pool of wave functions, which makes the approach essentially free from adjustable parameters and thus truly ab initio. The present results on magnetic interactions suggest that QMC is capable of accurately describing ground state properties of strongly correlated materials.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Energy density matrix formalism for interacting quantum systems: a quantum Monte Carlo study

We develop an energy density matrix that parallels the one-body reduced density matrix (1RDM) for many-body quantum systems. Just as the density matrix gives access to the number density and occupation numbers, the energy density matrix yields the energy density and orbital occupation energies. The eigenvectors of the matrix provide a natural orbital partitioning of the energy density while the eigenvalues comprise a single particle energy spectrum obeying a total energy sum rule. For mean-field systems the energy density matrix recovers the exact spectrum. When correlation becomes important, the occupation energies resemble quasiparticle energies in some respects. We explore the occupation energy spectrum for the finite 3D homogeneous electron gas in the metallic regime and an isolated oxygen atom with ground state quantum Monte Carlo techniques implemented in the QMCPACK simulation code. The occupation energy spectrum for the homogeneous electron gas can be described by an effective mass below the Fermi level. Above the Fermi level evanescent behavior in the occupation energies is observed in similar fashion to the occupation numbers of the 1RDM. A direct comparison with total energy differences shows a quantitative connection between the occupation energies and electron addition and removal energies for the electron gas. For the oxygen atom, the association between the ground state occupation energies and particle addition and removal energies becomes only qualitative. The energy density matrix provides a new avenue for describing energetics with quantum Monte Carlo methods which have traditionally been limited to total energies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Oxygen diffusion pathways in brownmillerite SrCoO2.5: Influence of structure and chemical potential

To design and discover new materials for next-generation energy materials such as solid-oxide fuel cells (SOFCs), a fundamental understanding of their ionic properties and behaviors is essential. The potential applicability of a material for SOFCs is critically determined by the activation energy barrier of oxygen along various diffusion pathways. In this work, we investigate interstitial-oxygen (Oi) diffusion in brownmillerite oxide SrCoO2.5, employing a first-principles approach. Our calculations indicate highly anisotropic ionic diffusion pathways, which result from its anisotropic crystal structure. The one-dimensional-ordered oxygen vacancy channels are found to provide the easiest diffusion pathway with an activation energy barrier height of 0.62 eV. The directions perpendicular to the vacancy channels have higher energy barriers for Oint diffusion. In addition, we have studied migration barriers for oxygen vacancies that could be present as point defects within the material. This in turn could also facilitate the transport of oxygen. Interestingly, for oxygen vacancies, the lowest barrier height was found to occur within the octahedral layer with an energy of 0.82 eV. Our results imply that interstitial migration would be highly one-dimensional in nature. Oxygen vacancy transport, on the other hand, could preferentially occur in the two-dimensional octahedral plane.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Successes and failures of Hubbard-corrected density functional theory: The case of Mg doped LiCoO$_2$

We have evaluated the successes and failures of the Hubbard-corrected density functional theory (DFT+U) approach to study Mg doping of LiCoO$_2$. We computed the effect of the U parameter on the energetic, geometric and electronic properties of two possible doping mechanisms: (1) substitution of Mg onto a Co (or Li) site with an associated impurity state and, (2) formation of impurity-state-free complexes of substitutional Mg and point defects in LiCoO$_2$. We find that formation of impurity states results in changes on the valency of Co in LiCoO$_2$. Variation of the Co U shifts the energy of the impurity state, resulting in energetic, geometric and electronic properties that depend significantly on the specific value of U. In contrast, the properties of the impurity-state-free complexes are insensitive to U. These results identify reasons for the strong dependence on the doping properties on the chosen value of U and for the overall difficulty of achieving agreement with the experimentally known energetic and electronic properties of doped transition metal oxides such as LiCoO$_2$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Generalizing the self-healing diffusion Monte Carlo approach to finite temperature: a path for the optimization of low-energy many-body bases

A statistical method is derived for the calculation of thermodynamic properties of many-body systems at low temperatures. This method is based on the self-healing diffusion Monte Carlo method for complex functions [F. A. Reboredo J. Chem. Phys. 136, 204101 (2012)] and some ideas of the correlation function Monte Carlo approach [D. M. Ceperley and B. Bernu, J. Chem. Phys. 89, 6316 (1988)]. In order to allow the evolution in imaginary time to describe the density matrix, we remove the fixed-node restriction using complex antisymmetric guiding wave functions. In the process we obtain a parallel algorithm that optimizes a small subspace of the many-body Hilbert space to have maximum overlap with the subspace spanned by the lowest-energy eigenstates of a many-body Hamiltonian. We show in a model system that the partition function is progressively maximized within this subspace. We show that the subspace spanned by the small basis systematically converges towards the subspace spanned by the lowest energy eigenstates. Possible applications of this method to calculate the thermodynamic properties of many-body systems near the ground state are discussed. The resulting basis can be also used to accelerate the calculation of the ground or excited states with Quantum Monte Carlo.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Ground state and Spin-Wave dynamics in Brownmillerite SrCoO2.5, A combined Hybrid Functional and LSDAU study

We theoretically investigate the ground state magnetic properties of the brownmillerite phase of SrCoO2.5. Strong correlations within Co d electrons are treated within the local spin density approximations of Density Functional theory (DFT) with Hubbard U corrections (LSDAU) and results are compared with the Heyd Scuzeria Ernzerhof (HSE) functional. The parameters computed with a U value of 7.5 eV are found to match closely to those computed within the HSE functional. A G type antiferromagnetic structure is found to be the most stable one, consistent with experimental observation. By mapping the total energies of different magnetic configurations onto a Heisenberg Hamiltonian we compute the magnetic exchange interaction parameters, J, between the nearest neighbor Co atoms. The J s obtained are then used to compute the spin wave frequencies and inelastic neutron scattering intensities. Among four spin wave branches, the lowest energy mode was found to have the largest scattering intensity at the magnetic zone center, while the other modes becomes dominant at different momenta. These predictions can be tested by experimentally.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Superconducting instability of a non-magnetic metallic band in an antiferromagnetic background

It is shown that a non-magnetic metallic band in the presence of an antiferromagnetic background coupled only by the exchange interaction develops a superconducting instability similar to the one described by BCS theory plus additional terms that strongly renormalize the spin excitation spectra. A Bardeen-Pines-like hole-hole interaction Hamiltonian, mediated by magnetic excitations, is deduced from a microscopic model of a fermion band and a spin band that interact with each other only via the exchange interaction. The theory shows the appearance of an attractive interaction when the Fermi velocity in the non-magnetic band is larger than the magnon velocity in the magnetic band. The electron-magnon scattering is suppressed by the appearance of a spin gap simultaneously with the superconducting state. Although this model may well describe a general class of materials to be discovered, the possibility that this theory could describe superconducting cuprates is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Diffusion quantum Monte Carlo study of the equation of state and point defects in aluminum

The many-body diffusion quantum Monte Carlo (DMC) method with twist-averaged boundary conditions is used to calculate the ground-state equation of state and the energetics of point defects in fcc aluminum using supercells up to 1331 atoms. The DMC equilibrium lattice constant differs from experiment by 0.008 A, or 0.2%, while the cohesive energy using DMC with backflow wave functions with improved nodal surfaces differs by 27 meV. DMC-calculated defect formation and migration energies agree with available experimental data, except for the nearest-neighbor divacancy, which is found to be energetically unstable, in agreement with previous density functional theory (DFT) calculations. DMC and DFT calculations of vacancy defects are in reasonably close agreement. Self-interstitial formation energies have larger differences between DMC and DFT, of up to 0.33eV, at the tetrahedral site. We also computed formation energies of helium interstitial defects where energies differed by up to 0.34eV, also at the tetrahedral site. The close agreement with available experiments demonstrates that DMC can be used as a predictive method to obtain benchmark energetics of defects in metals.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Simple Impurity Embedded in a Spherical Jellium: Approximations of Density Functional Theory compared to Quantum Monte Carlo Benchmarks

We study the electronic structure of a spherical jellium in the presence of a central Gaussian impurity. We test how well the resulting inhomogeneity effects beyond spherical jellium are reproduced by several approximations of density functional theory (DFT). Four rungs of Perdew's ladder of DFT functionals, namely local density approximation (LDA), generalized gradient approximation (GGA), meta-GGA and orbital-dependent hybrid functionals are compared against our quantum Monte Carlo (QMC) benchmarks. We identify several distinct transitions in the ground state of the system as the electronic occupation changes between delocalized and localized states. We examine the parameter space of realistic densities ($1 \le r_s\le 5$) and moderate depths of the Gaussian impurity ($Z<7$). The selected 18 electron system (with closed-shell ground state) presents $1d \to 2s$ transitions while the 30 electron system (with open-shell ground state) exhibits $1f \to 2p$ transitions. For the former system, the accuracy for the transitions is clearly improving with increasing sophistication of functionals with meta-GGA and hybrid functionals having only small deviations from QMC. However, for the latter system, we find much larger differences for the underlying transitions between our pool of DFT functionals and QMC. We attribute this failure to treatment of the exact exchange within these functionals. Additionally, we amplify the inhomogeneity effects by creating the system with spherical shell which leads to even larger errors in DFT approximations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Quantum Monte Carlo Calculations of Dihydrogen Binding Energetics on Ca Cations: an Assessment of Errors in Density Functionals for Weakly Bonded Systems

We investigate the binding of single and quadruple hydrogen molecules on a positively charged Ca ion. By comparing with benchmark quantum Monte Carlo (QMC) calculations we demonstrate wide variability in other more approximate electronic structure methods including common density functionals. Single determinant QMC calculations find no binding at short range by approximately 0.1 eV for the quadruple hydrogen molecule case, for a fixed hydrogen bond length of 0.77 Angstrom. Density functional calculations using common functionals such a LDA and B3LYP differ substantially from the QMC binding curve. We show that use of full Hartree-Fock exchange and PBE correlation(HFX+PBEC) obtains close agreement with the QMC results, both qualitatively and quantitatively. These results both motivate the use and development of improved functionals and indicate that caution is required applying electronic structure methods to weakly bound systems such as hydrogen storage materials based on metal ion decorated nanostructures.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Real-time dynamics of particle-hole excitations in Mott insulator-metal junctions

Charge excitations in Mott insulators (MIs) are distinct from their band-insulator counterparts and they can provide a mechanism for energy harvesting in solar cells based on strongly correlated electronic materials. In this paper, we study the real-time dynamics of holon-doublon pairs in a MI connected to metallic leads using the time-dependent density matrix renormalization group method. The transfer of charge across the MI-metal interface is controlled by both the electron-electron interaction strength within the MI and the voltage difference between the leads. We find an overall enhancement of the charge transfer as compared to the case of a (noninteracting) band insulator-metal interface with a similar band gap. Moreover, the propagation of holon-doublon excitations within the MI dynamically changes the spin-spin correlations, introducing time-dependent phase shifts in the spin structure factor.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Systematic reduction of sign errors in many-body calculations of atoms and molecules

The self-healing diffusion Monte Carlo algorithm (SHDMC) [Phys. Rev. B {\bf 79}, 195117 (2009), {\it ibid.} {\bf 80}, 125110 (2009)] is shown to be an accurate and robust method for calculating the ground state of atoms and molecules. By direct comparison with accurate configuration interaction results for the oxygen atom we show that SHDMC converges systematically towards the ground-state wave function. We present results for the challenging N$_2$ molecule, where the binding energies obtained via both energy minimization and SHDMC are near chemical accuracy (1 kcal/mol). Moreover, we demonstrate that SHDMC is robust enough to find the nodal surface for systems at least as large as C$_{20}$ starting from random coefficients. SHDMC is a linear-scaling method, in the degrees of freedom of the nodes, that systematically reduces the fermion sign problem.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Controlling the gap of fullerene microcrystals by applying pressure: the role of many-body effects

We studied theoretically the optical properties of C$_{60}$ fullerene microcrystals as a function of hydrostatic pressure with first-principles many-body theories. Calculations of the electronic properties were done in the GW approximation. We computed electronic excited states in the crystal by diagonalizing the Bethe-Salpeter equation (BSE). Our results confirmed the existence of bound excitons in the crystal. Both the electronic gap and optical gap decrease continuously and non-linearly as pressure of up to 6 GPa is applied. As a result, the absorption spectrum shows strong redshift. We also obtained that "negative" pressure shows the opposite behavior: the gaps increase and the optical spectrum shifts toward the blue end of the spectrum. Negative pressure can be realized by adding cubane (C$_8$H$_8$) or other molecules with similar size to the interstitials of the microcrystal. For the moderate lattice distortions studied here, we found that the optical properties of fullerene microcrystals with intercalated cubane are similar to the ones of an expanded undoped microcrystal. Based on these findings, we propose doped C60 as active element in piezo-optical devices.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Self-healing diffusion quantum Monte Carlo algorithms: methods for direct reduction of the fermion sign error in electronic structure calculations

We develop a formalism and present an algorithm for optimization of the trial wave-function used in fixed-node diffusion quantum Monte Carlo (DMC) methods. We take advantage of a basic property of the walker configuration distribution generated in a DMC calculation, to (i) project-out a multi-determinant expansion of the fixed-node ground-state wave function and (ii) to define a cost function that relates the fixed-node ground-state and the non-interacting trial wave functions. We show that (a) locally smoothing out the kink of the fixed-node ground-state wave function at the node generates a new trial wave-function with better nodal structure and (b) we argue that the noise in the fixed-node wave-function resulting from finite sampling plays a beneficial role, allowing the nodes to adjust towards the ones of the exact many-body ground state in a simulated annealing-like process. We propose a method to improve both single determinant and multi-determinant expansions of the trial wave-function. We test the method in a model system where benchmark configuration interaction calculations can be performed. Comparing the DMC calculations with the exact solutions, we find that the trial wave-function is systematically improved. The overlap of the optimized trial wave function and the exact ground state converges to 100% even starting from wave-functions orthogonal to the exact ground state. In the optimization process we find an optimal non-interacting nodal potential of density-functional-like form whose existence was predicted earlier[Phys.Rev. B {\bf 77}, 245110 (2008)]. We obtain the exact Kohn-Sham effective potential from the DMC data.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Systematic reduction of sign errors in many-body problems: generalization of self-healing diffusion Monte Carlo to excited states

A recently developed self-healing diffusion Monte Carlo algorithm [PRB 79, 195117] is extended to the calculation of excited states. The formalism is based on an excited-state fixed-node approximation and the mixed estimator of the excited-state probability density. The fixed-node ground state wave-functions of inequivalent nodal pockets are found simultaneously using a recursive approach. The decay of the wave-function into lower energy states is prevented using two methods: i) The projection of the improved trial-wave function into previously calculated eigenstates is removed. ii) The reference energy for each nodal pocket is adjusted in order to create a kink in the global fixed-node wave-function which, when locally smoothed out, increases the volume of the higher energy pockets at the expense of the lower energy ones until the energies of every pocket become equal. This reference energy method is designed to find nodal structures that are local minima for arbitrary fluctuations of the nodes within a given nodal topology. We demonstrate in a model system that the algorithm converges to many-body eigenstates in bosonic-like and fermionic cases.

0 citations0 reviewsNo code linkedOpen access