Response properties in phaseless auxiliary field quantum Monte Carlo

Abstract: We present a method for calculating first-order response properties in phaseless auxiliary field quantum Monte Carlo (AFQMC) through the application of automatic differentiation (AD). Biases and statistical efficiency of the resulting estimators are discussed. Our approach demonstrates that AD enables the calculation of reduced density matrices (RDMs) with the same computational cost scaling as energy calculations, accompanied by a cost prefactor of less than four in our numerical calculations. We investigate the role of self-consistency and trial orbital choice in property calculations.