We introduce a numerical algorithm to stochastically sample the dual fermion perturbation series around the dynamical mean field theory, generating all topologies of two-particle interaction vertices. We show results in the weak and strong coupling regime of the half-filled Hubbard model in two dimensions, illustrating that the method converges quickly where dynamical mean field theory is a good approximation, and show that corrections are large in the strong correlation regime at intermediate interaction. The fast convergence of dual corrections to dynamical mean field results illustrates the power of the approach and opens a practical avenue towards the systematic inclusion of non-local correlations in correlated materials simulations. An analysis of the frequency scale shows that only low-frequency propagators contribute substantially to the diagrams, putting the inclusion of higher order vertices within reach.
The open source ALPS (Algorithms and Libraries for Physics Simulations) project provides a collection of physics libraries and applications, with a focus on simulations of lattice models and strongly correlated systems. The libraries provide a convenient set of well-documented and reusable components for developing condensed matter physics simulation code, and the applications strive to make commonly used and proven computational algorithms available to a non-expert community. In this paper we present an updated and refactored version of the core ALPS libraries geared at the computational physics software development community, rewritten with focus on documentation, ease of installation, and software maintainability.
PROGRAM SUMMARY
Problems of finite temperature quantum statistical mechanics can be formulated in terms of imaginary (Euclidean) time Green's functions and self-energies. In the context of realistic Hamiltonians, the large energy scale of the Hamiltonian (as compared to temperature) necessitates a very precise representation of these functions. In this paper, we explore the representation of Green's functions and self-energies in terms of series of Chebyshev polynomials. We show that many operations, including convolutions, Fourier transforms, and the solution of the Dyson equation, can straightforwardly be expressed in terms of the series expansion coefficients. We then compare the accuracy of the Chebyshev representation for realistic systems with the uniform-power grid representation, which is most commonly used in this context. arXiv:1805.03521v2 [cond-mat.stat-mech]
We present an implementation of the self-energy embedding theory (SEET) for periodic systems and provide a fully self-consistent embedding solution for a simple realistic periodic problem -1D crystalline hydrogen -that displays many of the features present in complex real materials. For this system, we observe a remarkable agreement between our finite temperature periodic implementation results and well established and accurate zero temperature auxiliary quantum Monte Carlo data extrapolated to thermodynamic limit. We discuss differences and similarities with other Green's function embedding methods and provide the detailed algorithmic steps crucial for highly accurate and reproducible results.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.