The NewMag crystal-field code for f-element systems: Implementation for extended active spaces, second-order correlated energies, and generalisation to fn configurations

Abstract: We report a massive update of the NewMag program, which enables to extract Stevens crystal-field parameters (CFPs) after relativistic and multiconfigurational calculations are performed. The code can now post-treat ORCA and OpenMolcas outputs that contain state-average complete active space self-consistend field (SA-CASSCF) calculations with minimal or extended active spaces, second-order NEVPT2 or CASPT2 calculations, and spin-orbit configuration interaction (SOCI) calculations. Stevens parameters are extracted following the original Stevens convention with the so-called "extended" parameters. Rotation invariant indicators are also computed, after applying a correction of the parameter values to ensure normalization, following Rudowicz. At the SOCI level, the spin-orbit coupling (SOC) constant is also extracted within a spherical approximation of the SOC operator. For all the reference calculation levels, the model spectrum is reconstructed, allowing a direct assessment of its quality. The reported implementation is successfully applicable to f-element systems with a non-void, non-half-filled/empty or non-full f shell, that is for f^n configurations with n != 0, != 7 and != 14. Similarities and differences with the SINGLE_ANISO and AILFT codes are discussed. With selected examples, we showcase the interest of determining CFPs in f-element systems to understand their magnetic and optical properties in general, and more specifically the added value of NewMag. Finally, application of this approach to d-element systems is discussed, to reveal when it readily and successfully applies and when it may fail in reproducing satisfactorily the ab initio energies.