Density-Matrix Embedding Based Multi-Configurational Perturbation Theory Approach to Single-Ion Magnets

Abstract: Multi-configurational wave-function theory (MC-WFT) that combines complete active space self-consistent field (CASSCF) approach with subsequent state interaction (SI) treatment of spin-orbit coupling (SOC), abbreviated as CASSCF-SO, plays important roles in microscopic understanding of single-ion magnets (SIMs) with different central transition metal or lanthanide ions and various coordination environments, but its application to SIMs with complex structure is severely limited due to its highly demanding computational cost. Density-matrix embedding theory (DMET) provides a systematic and mathematically rigorous framework to combine low-level mean field approaches like Hartree-Fock and high-level MC-WFT methods like CASSCF-SO, which is particularly promising to SIMs. As a continuation of our previous work on DMET+CASSCF for $3d$ SIMs (Ai, Sun, and Jiang, J. Phys. Chem. Lett. 2022, 13, 10627), we extend the methodology by considering dynamic correlation on top of CASSCF using the second-order $n$-electron valence perturbation theory (NEVPT2) in the DMET framework, abbreviated as DMET+NEVPT2, and benchmark the accuracy of this approach to molecular magnetic anisotropy in a set of typical transition metal complexes. We found that DMET+NEVPT2 can give the results very close to all-electron treatment, and can be systematically improved for higher accuracy by expanding the region treated as the central cluster, while the computation cost is dramatically reduced due to the reduction of the number of orbitals by DMET construction. Our findings suggest that DMET is capable of accounting for most of the dynamic correlation that is important for magnetic anisotropy in typical SIMs, and can be useful for further high-accuracy spin-phonon study and high-throughput computations.