Linear-scaling fixed-node diffusion quantum Monte Carlo: accounting for the nodal information in a density matrix-based scheme.

A reformulation of the fixed-node diffusion quantum Monte Carlo method (FN-DQMC) in terms of the N-particle density matrix is presented, which allows us to reduce the computational effort to linear for the evaluation of the local energy. The reformulation is based on our recently introduced density matrix-based approach for a linear-scaling variational QMC method [J. Kussmann et al., Phys. Rev. B. 75, 165107 (2007)]. However, within the latter approach of using the positive semi-definite N-particle trial density (rhoN T(R)=mid R:Psi(T)(R)mid R:(2)), the nodal information of the trial function is lost. Therefore, a straightforward application to the FN-DQMC method is not possible, in which the sign of the trial function is usually traced in order to confine the random walkers to their nodal pockets. As a solution, we reformulate the FN-DQMC approach in terms of off-diagonal elements of the N-particle density matrix rhoN T(R;R'), so that the nodal information of the trial density matrix is obtained. Besides all-electron moves, a scheme to perform single-electron moves within N-PDM QMC is described in detail. The efficiency of our method is illustrated for exemplary calculations.