This article aims to compare the so‐called iterated perturbation theory (IPT) and auxiliary master equation approach (AMEA) impurity solvers for a Mott insulating system driven out of equilibrium by a static electric field. Electronic heat bath and optical phonons are employed as dissipation mechanism of the current‐induced Joule heat that the excited electrons of the lattice experience as the result of the field's driving. Despite its simplicity, the IPT approach yields results that qualitatively are in good agreement with those obtained within the AMEA impurity solver, although fails to reproduce some correlation effects.
This article aims to compare the so‐called iterated perturbation theory (IPT) and auxiliary master equation approach (AMEA) impurity solvers for a Mott insulating system driven out of equilibrium by a static electric field. Electronic heat bath and optical phonons are employed as dissipation mechanism of the current‐induced Joule heat that the excited electrons of the lattice experience as the result of the field's driving. Despite its simplicity, the IPT approach yields results that qualitatively are in good agreement with those obtained within the AMEA impurity solver, although fails to reproduce some correlation effects.
We characterize the current-carrying nonequilibrium steady state (NESS) in a single-band Hubbard model confronted with a static electric field in the presence of quenched disorder. Beyond the linear response regime, the electric field amplitude must compensate for at least half of the band gap to have a non-negligible stationary current. As the disorder is not expected to dissipate the extra energy injected by the field, optical phonons assisted by a fermionic heat bath serve as dissipation channels for the current-induced Joule heat generated by the accelerated electrons. The NESS of the system is addressed employing the dynamical mean-field theory using the so-called auxiliary master equation approach as impurity solver. Disorder effects are treated locally via the coherent potential approximation (CPA) and the self-consistent Born (SCB) approach. In the regime in which the two schemes yield similar results, we employ the SCB as it is computationally cheaper than the CPA. We show that, in a purely electronic setup, the disorder-induced dephasing cannot contribute states within the gap but only smear out the edges of the Hubbard bands. When phonons are taken into account, the different nature of disorder-induced dephasing and phonon-related dissipation becomes clear. We show that although both disorder and electron-phonon interaction enhance the current at off-resonant fields, disorder effects play a marginal role since they cannot provide in-gap states which are instead brought about by phonons and represent the privileged relaxation pathway for excited electrons.
We characterize the response of a Mott insulating system to a static electric field in terms of its conducting and spectral properties. Dissipation is included by a coupling to fermionic baths and to either optical or acoustic phonons. This paper extends and completes the analysis made in a previous work by the authors [Mazzocchi et al., Phys. Rev. B 106, 125123 (2022)]. In the present work, phonons are included diagrammatically within the Migdal approximation by also including self-consistency from the electronic feedback. The nonequilibrium steady state is addressed by means of the dynamical mean-field theory based on the nonequilibrium Green's function approach, while the so-called auxiliary master equation approach is employed as impurity solver. With optical phonons the self-consistency suppresses the steady-state current for field strengths comparable to the band gap with respect to the non-self-consistent case. This is due to the interaction of phonons with the hot electrons of the lattice which increases their temperature, thus providing a less effective relaxation channel for the current-induced Joule heat. In addition, in the case of optical phonons, the results are essentially independent of the temperature of the fermionic baths, as the latter is sensibly smaller than their characteristic frequency. On the other hand, with acoustic phonons the steady-state current is slightly suppressed by the self-consistent treatment only at field strengths close to half of the gap and especially at very small phonon frequency. Also, in this case, the results seem to slightly depend on the temperature of the fermionic baths.
We study the spectral and transport properties of a Mott insulator driven by a static electric field into a nonequilibrium steady state. For the dissipation, we consider two mechanisms: wide band fermion reservoirs and phonons included within the Migdal approximation. The electron correlations are treated via nonequilibrium dynamical mean field theory with an impurity solver suitable for strong correlations. We find that dissipation via phonons is limited to restricted ranges of field values around Wannier-Stark resonances. To cover the full range of field strengths, we allow for a small coupling with fermionic baths, which stabilizes the solution. When considering both dissipation mechanisms, we find that phonons enhance the current for field strengths close to half of the gap while lowering it at the gap resonance as compared to the purely electronic dissipation used by Murakami and Werner [Murakami and Werner, Phys. Rev. B 98, 075102 (2018)]. Once phonons are the only dissipation mechanism, the current in the metallic phase is almost one order of magnitude smaller than the typical values obtained by coupling to a fermionic bath. In this case, the transport regime is characterized by an accumulation of charge in the upper Hubbard band described by two effective chemical potentials.
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.
In this work, we develop a non-equilibrium steady-state non-crossing approximation (NESS-NCA) impurity solver applicable to general impurity problems. The choice of the NCA as the impurity solver enables both a more accurate description of correlation effects with larger Coulomb interaction and scalability to multi-orbital systems. Based on this development, we investigate strongly correlated non-equilibrium states of a dissipative lattice system under constant electric fields. Both the electronic Coulomb interaction and the electric field are treated non-perturbatively using dynamical mean-field theory in its non-equilibrium steady-state form (NESS-DMFT) with the NESS-NCA impurity solver. We validate our implementation using a half-filled single-band Hubbard model attached to a fictitious free Fermion reservoir, which prevents temperature divergence. As a result, we identify metallic and insulating phases as functions of the electric field and the Coulomb interaction along with a phase coexistence region amid the metal-to-insulator transition (MIT). We find that the MIT driven by the electric field is qualitatively similar to the equilibrium MIT as a function of temperature, differing from results in previous studies using the iterative perturbation theory (IPT) impurity solver. Finally, we highlight the importance of the morphology of a correlated system under the influence of an electric field.
We study a model for photovoltaic energy collection consisting of a Mott-insulating layer in presence of acoustic phonons, coupled to two wide-band fermion leads at different chemical potentials and driven into a nonequilibrium steady state by a periodic electric field. We treat electron correlations with nonequilibrium dynamical mean-field theory using the so-called auxiliary master-equation approach as impurity solver and include dissipation by acoustic phonons via the Migdal approximation. For a small hybridization to the leads, we obtain a peak in the photocurrent as a function of the driving frequency which can be associated with impact ionization processes. For larger hybridizations the shallow peak suggests a suppression of impact ionization with respect to direct photovoltaic excitations. Acoustic phonons slightly enhance the photocurrent for small driving frequencies and suppress it at frequencies around the main peak at all considered hybridization strengths.