The theory discussed in the previous section, including the extension to RASPT2, has been implemented in a development version of OpenMolcas. In the following, we present and discuss the numerical results obtained with it. All calculations employed the frozen core approximation in association with the atomic compact Cholesky decomposition to generate an on-the-fly auxiliary basis set for the resolution-of-identity treatment of the electron repulsion integrals. The ionization potential-electron affinity (IPEA) shift and symmetry constraints were not used, but an imaginary level shift parameter of 0.2i was selected for all calculations. For XDW, the parameter ζ is set to 50 for all calculations. The selection of the parameter affects the computed results, but the impact of the parameter has been examined in the past literature, and therefore our primary interest is not to examine this influence. In the reference SA-CASSCF calculation, the weight for each state is equal regardless of the number of averaged states.