We present an algorithm to evaluate Matsubara sums for Feynman diagrams comprised of bare Green's functions with single-band dispersions with local U Hubbard interaction vertices. The algorithm provides an exact construction of the analytic result for the frequency integrals of a diagram that can then be evaluated for all parameters U , temperature T , chemical potential µ, external frequencies and internal/external momenta. This method allows for symbolic analytic continuation of results to the real frequency axis, avoiding any ill-posed numerical procedure. When combined with diagrammatic Monte-Carlo, this method can be used to simultaneously evaluate diagrams throughout the entire T -U -µ phase space of Hubbard-like models at minimal computational expense.

We present a general method to optimize the evaluation of Feynman diagrammatic expansions, which requires the automated symbolic assignment of momentum/energy conserving variables to each diagram. With this symbolic representation, we utilize the pole structure of each diagram to automatically sort the Feynman diagrams into groups that are likely to contain nearly equal or nearly cancelling diagrams, and we show that for some systems this cancellation is exact. This allows for a potentially massive cancellation during the numerical integration of internal momenta variables, leading to an optimal suppression of the 'sign problem' and hence reducing the computational cost. Although we define these groups using a frequency space representation, the equality or cancellation of diagrams within the group remains valid in other representations such as imaginary time used in standard diagrammatic Monte Carlo. As an application of the approach we apply this method, combined with algorithmic Matsubara integration (AMI) [Phys. Rev. B 99, 035120 (2019)] and Monte Carlo methods, to the Hubbard model self-energy expansion on a 2D square lattice up to sixth order which we evaluate and compare with existing benchmarks. arXiv:1911.11129v1 [cond-mat.str-el]

We report a scalable Fortran implementation of the phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) and demonstrate its excellent performance and beneficial scaling with respect to system size. Furthermore, we investigate modifications of the phaseless approximation that can help to reduce the overcorrelation problems common to the ph-AFQMC. We apply the method to the 26 molecules in the HEAT set, the benzene molecule, and water clusters. We observe a mean absolute deviation of the total energy of 1.15 kcal/mol for the molecules in the HEAT set, close to chemical accuracy. For the benzene molecule, the modified algorithm despite using a single-Slater-determinant trial wavefunction yields the same accuracy as the original phaseless scheme with 400 Slater determinants. Despite these improvements, we find systematic errors for the CN, CO 2 , and O 2 molecules that need to be addressed with more accurate trial wavefunctions. For water clusters, we find that the ph-AFQMC yields excellent binding energies that differ from CCSD(T) by typically less than 0.5 kcal/mol.

We implement the phaseless auxiliary field quantum Monte Carlo method using the plane-wave based projector augmented wave method and explore the accuracy and the feasibility of applying our implementation to solids. We use a singular value decomposition to compress the two-body Hamiltonian and, thus, reduce the computational cost. Consistent correlation energies from the primitive-cell sampling and the corresponding supercell calculations numerically verify our implementation. We calculate the equation of state for diamond and the correlation energies for a range of prototypical solid materials. A down-sampling technique along with natural orbitals accelerates the convergence with respect to the number of orbitals and crystal momentum points. We illustrate the competitiveness of our implementation in accuracy and computational cost for dense crystal momentum point meshes compared to a well-established quantum-chemistry approach, the coupled-cluster ansatz including singles, doubles, and perturbative triple particle–hole excitation operators.

No abstract