A quantum Monte Carlo method that combines the second-order many-body perturbation theory and Monte Carlo (MC) integration has been developed for correlation and correlation-corrected (quasiparticle) energy bands of one-dimensional solids. The sum-of-product expressions of correlation energy and self-energy are transformed, with the aid of a Laplace transform, into high-dimensional integrals, which are subject to a highly scalable MC integration with the Metropolis algorithm for importance sampling. The method can compute correlation energies of polyacetylene and polyethylene within a few $\text{m}{E}_{\text{h}}$ and quasiparticle energy bands within a few tenths of an eV. It does not suffer from the fermion sign problem and its description can be systematically improved by raising the perturbation order.

• Received 31 August 2014
• Revised 4 November 2014

DOI:https://doi.org/10.1103/PhysRevB.90.201110