We report on the findings of a blind challenge devoted to determining the frozencore, full configuration interaction (FCI) ground state energy of the benzene molecule in a standard correlation-consistent basis set of double-ζ quality. As a broad international endeavour, our suite of wave function-based correlation methods collectively represents a diverse view of the high-accuracy repertoire offered by modern electronic structure theory. In our assessment, the evaluated high-level methods are all found to qualitatively agree on a final correlation energy, with most methods yielding an estimate of the FCI value around −863 mE H. However, we find the root-mean-square deviation of the energies from the studied methods to be considerable (1.3 mE H), which in light of the acclaimed performance of each of the methods for smaller molecular systems clearly displays the challenges faced in extending reliable, near-exact correlation methods to larger systems. While the discrepancies exposed by our study thus emphasize the fact that the current state-of-the-art approaches leave room for improvement, we still expect the present assessment to provide a valuable community resource for benchmark and calibration purposes going forward.

The semistochastic heat-bath configuration interaction method is a selected configuration interaction plus perturbation theory method that has provided near-full configuration interaction (FCI) levels of accuracy for many systems with both single- and multi-reference character. However, obtaining accurate energies in the complete basis-set limit is hindered by the slow convergence of the FCI energy with respect to basis size. Here, we show that the recently developed basis-set correction method based on range-separated density functional theory can be used to significantly speed up basis-set convergence in SHCI calculations. In particular, we study two such schemes that differ in the functional used and apply them to transition metal atoms and monoxides to obtain total, ionization, and dissociation energies well converged to the complete-basis-set limit within chemical accuracy.

We present a version of the T-moves approach for treating nonlocal pseudopotentials in diffusion Monte Carlo, which has much smaller time-step errors than the existing T-moves approaches, while at the same time preserving desirable features such as the upper-bound property for the energy. In addition, we modify the reweighting factor of the projector used in diffusion Monte Carlo to reduce the time-step error. The latter is applicable not only to pseudopotential calculations but also to all-electron calculations.

No abstract