Numerical renormalization group calculations for magnetic impurity systems with spin-orbit coupling and crystal-field effects

Exploiting symmetries in the numerical renormalization group (NRG) method significantly enhances performance by improving the accuracy, increasing the computational speed, and optimizing the memory efficiency. Published codes focus on continuous rotations and unitary groups, which generally are not applicable to systems with strong crystal-field effects. The PointGroupNRG code implements symmetries related to discrete rotation groups, which are defined by the user in terms of Clebsch-Gordan coefficients, together with particle conservation and spin rotation symmetries. In this paper we present a new version of the code that extends the available finite groups, previously limited to simply reducible point groups, in a way that all point and double groups become accessible. It also includes the full spin-orbital rotation group. Moreover, to improve the code's flexibility for impurities with complex interactions, this new version allows to choose between a standard Anderson Hamiltonian for the impurity or, as another novel feature, an ionic model that requires only the spectrum and the impurity Lehmann amplitudes.