Search

A recent paper on Fe2 [A. Kalemos, J. Chem. Phys. 142, 244304 (2015)] critiqued our previous work on the system [Hoyer et al., J. Chem. Phys. 141, 204309 (2014)]. In this comment, we explain the nature of our previously reported potential energy curve for Fe2 and we discuss our computed properties for Fe2. Additionally, we fix a labeling error that was present in our previous work, although this error is unrelated to the main point of discussion.

A multiconfigurational self-consistent field method based on the concept of generalized active space (GAS) is presented. GAS wave functions are obtained by defining an arbitrary number of active spaces with arbitrary occupation constraints. By a suitable choice of the GAS spaces, numerous ineffective configurations present in a large complete active space (CAS) can be removed, while keeping the important ones in the CI space. As a consequence, the GAS self-consistent field approach retains the accuracy of the CAS self-consistent field (CASSCF) ansatz and, at the same time, can deal with larger active spaces, which would be unaffordable at the CASSCF level. Test calculations on the Gd atom, Gd(2) molecule, and oxoMn(salen) complex are presented. They show that GAS wave functions achieve the same accuracy as CAS wave functions on systems that would be prohibitive at the CAS level., (© 2011 American Institute of Physics)

We propose a new approach to determine a suitable zeroth-order wavefunction for multiconfigurational perturbation theory. The same ansatz as in complete active space (CAS) wavefunction optimization is used but it is split in two parts, a principal space (A) and a much larger extended space (B). Löwdin's partitioning technique is employed to map the initial eigenvalue problem to a dimensionality equal to that of (A) only. Combined with a simplified expression for the (B) portion of the wavefunction, we are able to drastically reduce the storage and computational demands of the wavefunction optimization. This scheme is used to produce reference wavefunctions and energies for subsequent second-order perturbation theory (PT2) corrections. Releasing the constraint of computing the exact CAS energy and wavefunction prior to the PT2 treatment introduces a nonstandard paradigm for multiconfigurational methods. Based on the results of test calculations, we argue that principal parts with only few percents of the total number of CAS configurations could provide final multiconfigurational PT2 energies of the same accuracy as in the standard paradigm. In the future, algorithmic improvements for this scheme will bring into reach active spaces much beyond the present limit of CAS-based methods, therefore allowing for accurate studies of systems featuring strong correlation.