Orbital Magnetization of Interacting Electrons

We derive an exact expression for the orbital magnetization of electrons with short-range interactions (such as density-density interactions) in terms of exact zero-frequency response functions of the zero-field system. The result applies to weakly and strongly correlated electrons at zero and finite temperature, provided that the local grand potential density only depends on local thermodynamic parameters. We benchmark the formula for non-interacting and weakly-coupled electrons. To zeroth and first orders in the interaction strength, it agrees with the modern theory of orbital magnetization and its recent generalization to self-consistent Hartree-Fock bands. Our work provides an exact framework of interacting orbital magnetization beyond mean-field treatments, and paves the way for quantitative studies of strongly correlated electrons in external magnetic fields.