  The LDA+U method is widely used to study the properties of realistic solids
with strong electron correlations. One of its main shortcomings is that it does
not provide direct access to the temperature dependence of material properties
such as the Dzyaloshinskii-Moriya interaction, the anomalous Hall conductivity,
and the spin-orbit torque. While the method of spectral moments allows us in
principle to compute these quantities directly at finite temperatures, the
standard two-pole approximation can be applied only to Hamiltonians that are
effectively of single-band type. We do a first step to explore if the method of
spectral moments may replace the LDA+U method in first-principles calculations
of correlated solids with many bands in cases where the direct assessment of
the temperature dependence of equilibrium and response functions is desired:
The spectral moments of many-band Hamiltonians of correlated electrons do not
commute and therefore they do not possess a system of common eigenvectors. We
show that nevertheless the spectral function may be constructed from the
spectral moments by solving a system of coupled non-linear equations.
Additionally, we show how to compute the anomalous Hall conductivity of
correlated electrons from this spectral function. We demonstrate the method for
the Hubbard-Rashba model, where the standard two-pole approximation cannot be
applied, because spin-orbit interaction (SOI)couples the spin-up and the
spin-down bands. In the quest for new quantum states that arise from the
combination of SOI and correlation effects, the Hartree-Fock approximation is
frequently used to obtain a first approximation for the phase diagram. We
propose that using the many-band generalization of the selfconsistent moment
method instead of Hartree-Fock in such exploratory model calculations may
improve the accuracy significantly, while keeping the computational burden low.
