Abstract:We present a method to construct pseudo-BCS wave functions from the one-body density matrix. The resulting many-body wave function, which can be produced for any fermion systems, including those with purely repulsive interactions, has the form of a number-projected BCS form, or antisymmetrized geminal power (AGP). Such wave functions provide a better ansatz for correlated fermion systems than a single Slater determinant, and often better than a linear combination of Slater determinants (for example from a truncated active space calculation). We describe a procedure to build such a wave function conveniently from a given reduced density matrix of the system, rather than from a mean-field solution (which gives a Slater determinant for repulsive interactions). The pseudo-BCS wave function thus obtained reproduces the density matrix or minimizes the difference between the input and resulting density matrices. One application of the pseudo-BCS wave function is in auxiliary-field quantum Monte Carlo (AFQMC) calculations as the trial wave function to control the sign/phase problem. AFQMC is often among the most accurate general methods for correlated fermion systems. We show that the pseudo-BCS form further reduces the constraint bias and leads to improved accuracy compared to the usual Slater determinant trial wave functions, using the two-dimensional Hubbard model as an example. Furthermore, the pseudo-BCS trial wave function allows a new systematically improvable self-consistent approach, with pseudo-BCS trial wave function iteratively generated by AFQMC via the one-body density matrix. https://doi.org/10.1103/PhysRevResearch.3.013065 Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI. Published by the American Physical Society

Self-consistent optimization of the trial wave-function in constrained path auxiliary field Quantum Monte Carlo using mixed estimators

Self-consistent Optimization of the Trial Wave Function Within the Constrained Path Auxiliary Field Quantum Monte Carlo Method Using Mixed Estimators

Unbiasing Fermionic Auxiliary-Field Quantum Monte Carlo with Matrix Product State Trial Wavefunctions

Self-Consistent Solution Of Kohn-Sham Equation By Real-Space Finite-Difference Method

The Design of New Practical Constraints in Auxiliary-Field Quantum Monte Carlo

Coupling Quantum Monte Carlo and Independent-Particle Calculations: Self-consistent Constraint for the Sign Problem Based on the Density or the Density Matrix

Scalable Quantum Monte Carlo with Direct-Product Trial Wave Functions

Benchmark Study of the Two-Dimensional Hubbard Model with Auxiliary-Field Quantum Monte Carlo Method

Pseudo-BCS wave function from density matrix decomposition: Application in auxiliary-field quantum Monte Carlo

Auxiliary-field quantum Monte Carlo calculations with multiple-projector pseudopotentials

Re-anchoring Quantum Monte Carlo with Tensor-Train Sketching

Quantum Monte Carlo with reoptimized perturbatively selected configuration-interaction wave functions

Accelerating the Convergence of Auxiliary-Field Quantum Monte Carlo inSolids with Optimized Gaussian Basis Sets

Contextual Subspace Auxiliary-Field Quantum Monte Carlo: Improved bias with reduced quantum resources

Stochastic Multiconfigurational Self-Consistent Field Theory.

Some recent developments in auxiliary-field quantum Monte Carlo for real materials

An Auxiliary-Field Quantum Monte Carlo Perspective on the Ground State of the Dense Uniform Electron Gas: An Investigation with Hartree-Fock Trial Wavefunctions

Finite-Temperature Auxiliary-Field Quantum Monte Carlo: Self-Consistent Constraint And Systematic Approach To Low Temperatures

Benchmark Phaseless Auxiliary-Field Quantum Monte Carlo Method for Small Molecules

Auxiliary Field Quantum Monte Carlo for Dilute Neutrons on the Lattice