Predicting the physical properties of strongly correlated electron systems (SCES) requires the combined Density Functional Theory and Dynamical Mean-Field Theory (DFT+DMFT) framework. However, the computational bottle-neck lies in the impurity solver. This work focuses on rigorously benchmarking the Multi-Orbital Iterated Perturbation Theory (MO-IPT) impurity solver against the numerically exact Continuous-Time Quantum Monte Carlo (CTQMC) method, specifically through a detailed study of the three-orbital Hubbard model. The two-orbital model serves as a canonical system to test performance under realistic multi-band conditions. Our systematically derived MO-IPT overcomes the critical limitations of CTQMC, including exponential scaling with orbital number, the need for analytic continuation, and the severe fermionic sign problem. By demonstrating high accuracy across varying interaction strengths and temperatures, we validate MO-IPT as an inexpensive, yet reliable, replacement solver essential for future material-specific investigations.
Although several impurity solvers in the dynamical mean field theory (DMFT) have been proposed, especially in multiband systems, there are practical difficulties arising from a trade-off between numerical costs and reliability. In this study, we reinterpret the iterative perturbation theory (IPT) as an approximation which captures the strong correlation effects by mimicking the particular frequency structures of the exact full vertex, and extend it such that it can have efficiency and reliability simultaneously by modifying IPT vertex using the parquet equations. We apply this method to several models to evaluate their validity. We confirm that our method shows good agreements with the numerically exact continuous-time quantum Monte Carlo method in the single-site DMFT calculation.
We extend the Kajueter-Kotliar [] iterated perturbation theory (KK-IPT) away from half-filling to nonequilibrium steady states. We benchmark the resulting nonequilibrium KK-IPT approach against the auxiliary master equation approach (AMEA), whose accuracy is controlled in and out of equilibrium. As expected, in equilibrium, KK-IPT reproduces the AMEA results for different fillings with high accuracy at the level of both spectral properties and electron densities. Out of equilibrium, we study quantum transport across a correlated impurity and compute the differential conductance and spectral functions. We find very good agreement between nonequilibrium KK-IPT and AMEA in the parameter regime where the latter is reliable, in particular at moderate temperatures and biases. Although a controlled benchmark is not available in the low-temperature, low-bias regime, where AMEA becomes less reliable, we show that this nonequilibrium KK-IPT impurity solver satisfies the exact spectral sum rules for the first and second moments to high accuracy throughout the entire parameter range studied. These results support nonequilibrium KK-IPT as an approximate description of nonequilibrium steady states away from half-filling. At the same time, comparing against AMEA the double occupancy obtained from the nonequilibrium Galitskii-Migdal expression for the interaction energy shows that the deviation from AMEA remains small near half-filling for moderate and large values of the bias, but grows markedly away from half-filling, delineating the regime in which the method can be trusted quantitatively rather than merely qualitatively.
The local moment approach (LMA) has presented itself as a powerful semianalytical quantum impurity solver (QIS) in the context of the dynamical mean-field theory (DMFT) for the periodic Anderson model and it correctly captures the low-energy Kondo scale for the single impurity model, having excellent agreement with the Bethe ansatz and numerical renormalization group (NRG) results. However, the most common correlated lattice model, the Hubbard model, has not been explored well within the LMA+DMFT framework beyond the insulating phase. Here in our work, within the framework we complete the filling-interaction phase diagram of the single band Hubbard model at zero temperature. Our formalism is generic to any particle filling and can be extended to finite temperature. We contrast our results with another QIS, namely the iterated perturbation theory (IPT) and show that the second spectral moment sum rule improves better as the Hubbard interaction strength grows stronger in LMA, whereas it severely breaks down after the Mott transition in IPT. For the metallic case, the Fermi liquid (FL) scaling agreement with the NRG spectral density supports the fact that the FL scale emerges from the inherent Kondo physics of the impurity model. We also show that, in the metallic phase, the FL scaling of the spectral density leads to universality which extends to infinite frequency range at infinite correlation strength (strong coupling). At large interaction strength, the off half-filling spectral density forms a pseudogap near the Fermi level and filling-controlled Mott transition occurs as one approaches the half-filling. As a response property, we finally study the zero temperature optical conductivity and find universal features such as absorption peak position governed by the FL scale and a doping independent crossing point, often dubbed the isosbestic point in experiments.
We study the phase diagram of the ionic Hubbard model (IHM) at half filling on a Bethe lattice of infinite connectivity using dynamical mean-field theory (DMFT), with two impurity solvers, namely, iterated perturbation theory (IPT) and continuous time quantum Monte Carlo (CTQMC). The physics of the IHM is governed by the competition between the staggered ionic potential $\mathrm{\ensuremath{\Delta}}$ and the on-site Hubbard
U. We find that for a finite $\mathrm{\ensuremath{\Delta}}$ and at zero temperature, long-range antiferromagnetic (AFM) order sets in beyond a threshold
U=UAF via a first-order phase transition. For
U smaller than
UAF the system is a correlated band insulator. Both methods show a clear evidence for a quantum transition to a half-metal (HM) phase just after the AFM order is turned on, followed by the formation of an AFM insulator on further increasing
U. We show that the results obtained within both methods have good qualitative and quantitative consistency in the intermediate-to-strong-coupling regime at zero temperature as well as at finite temperature. On increasing the temperature, the AFM order is lost via a first-order phase transition at a transition temperature ${T}_{AF}(U,\mathrm{\ensuremath{\Delta}})$ [or, equivalently, on decreasing
U below ${U}_{AF}(T,\mathrm{\ensuremath{\Delta}})$], within both methods, for weak to intermediate values of
U/t. In the strongly correlated regime, where the effective low-energy Hamiltonian is the Heisenberg model, IPT is unable to capture the thermal (Neel) transition from the AFM phase to the paramagnetic phase, but the CTQMC does. At a finite temperature
T,DMFT+CTQMC shows a second phase transition (not seen within
DMFT+IPT) on increasing
U beyond
UAF. At
UN>UAF, when the Neel temperature
TN for the effective Heisenberg model becomes lower than
T, the AFM order is lost via a second-order transition. For $U\ensuremath{\gg}\mathrm{\ensuremath{\Delta}}, {T}_{N}\ensuremath{\sim}{t}^{2}/U(1\ensuremath{-}{x}^{2})$, where $x=2\mathrm{\ensuremath{\Delta}}/U$ and thus
TN increases with increase in $\mathrm{\ensuremath{\Delta}}/U$. In the three-dimensional parameter space of $(U/t,\phantom{\rule{0.28em}{0ex}}T/t,\phantom{\rule{4.pt}{0ex}}\text{and}\phantom{\rule{4.pt}{0ex}}\mathrm{\ensuremath{\Delta}}/t)$, as
T increases, the surface of first-order transition at ${U}_{AF}(T,\mathrm{\ensuremath{\Delta}})$ and that of the second-order transition at ${U}_{N}(T,\mathrm{\ensuremath{\Delta}})$ approach each other, shrinking the range over which the AFM order is stable. There is a line of tricritical points that separates the surfaces of first- and second-order phase transitions.
We solve the impurity problem which arises within nonequilibrium dynamical mean-field theory for the Hubbard model by means of a self-consistent perturbation expansion around the atomic limit. While the lowest order, known as the noncrossing approximation (NCA), is reliable only when the interaction
U is much larger than the bandwidth, low-order corrections to the NCA turn out to be sufficient to reproduce numerically exact Monte Carlo results in a wide parameter range that covers the insulating phase and the metal-insulator crossover regime at not too low temperatures. As an application of the perturbative strong-coupling impurity solver we investigate the response of the double occupancy in the Mott insulating phase of the Hubbard model to a dynamical change in the interaction or the hopping, a technique which has been used as a probe of the Mott insulating state in ultracold fermionic gases.