Track synthesis, materials and catalysis as it is published.
Search papers, labs, and topics across Lattice.
The authors updated the NewMag program to post-process relativistic multiconfigurational electronic structure outputs (SA-CASSCF, CASPT2, NEVPT2, and SOCI) from ORCA and OpenMolcas, extracting Stevens crystal-field parameters and spin-orbit coupling constants across arbitrary open-shell $f^n$ configurations ($n \neq 0, 7, 14$). Accurately mapping ab initio correlated wavefunctions to effective crystal-field Hamiltonians is critical for predicting magnetic anisotropy, optical transitions, and spin-relaxation dynamics in heavy-element complexes. The generalized implementation incorporates normalized rotation-invariant metrics, spherical spin-orbit approximations, and automated model spectrum reconstruction to directly quantify the fidelity of the parameterized Hamiltonian against reference electronic energies.
Effective spin Hamiltonians can now be extracted directly and reliably from second-order correlated relativistic wavefunctions across arbitrary $f^n$ configurations, eliminating manual active-space workarounds in magnetic property modeling.
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.