A study of the fixed-node error in quantum Monte Carlo calculations of electronic transitions: The case of the singlet n→π∗ (CO) transition of the acrolein

Abstract

We report fixed-node diffusion Monte Carlo (FN-DMC) calculations of the singlet n-->pi( *) (CO) vertical transition of acrolein. The impact of the fixed-node approximation on the excitation energy is investigated. To do that, trial wave functions corresponding to various nodal patterns are used. They are constructed by using either a minimal complete-active-space self-consistent field (CASSCF) calculation involving an oxygen lone pair n and the pi( *) (CO) molecular orbitals or a more complete set involving all the molecular orbitals expected to play a significant role in the excitation process. Calculations of both states have been performed with molecular orbitals optimized separately for each state via standard "state specific" CASSCF calculations or by using a common set of optimized orbitals ["state averaged" CASSCF calculations] whose effect is to introduce some important correlation between the nodal patterns of the two electronic states. To investigate the role of the basis set three different basis of increasing size have been employed. The comparative study based on the use of all possible combinations of basis sets, active spaces, and type of optimized molecular orbitals shows that the nodal error on the difference of energies is small when chemically relevant active space and state-averaged-type CASSCF wave functions are used, although the fixed-node error on the individual total energies involved can vary substantially. This remarkable result obtained for the acrolein suggests that FN-DMC calculations based on a simple strategy (use of standard ab initio wave functions and no Monte Carlo optimization of molecular orbital parameters) could be a working computational tool for computing electronic transition energies for more general systems.