We study the single-impurity Anderson model out of equilibrium under the influence of a bias voltage f and a magnetic field B. We investigate the interplay between the shift ( B w ) of the Kondo peak in the spin-resolved density of states (DOS) and the one ( B f ) of the conductance anomaly. In agreement with experiments and previous theoretical calculations we find that, while the latter displays a rather linear behavior with an almost constant slope as a function of B down to the Kondo scale, the DOS shift first features a slower increase reaching the same behavior as. Our auxiliary master equation approach yields highly accurate nonequilibrium results for the DOS and for the conductance all the way from within the Kondo up to the charge fluctuation regime, showing excellent agreement with a recently introduced scheme based on a combination of numerical renormalization group with time-dependent density matrix renormalization group.

We study equilibrium and nonequilibrium properties of the single-impurity Anderson model with a power-law pseudogap in the density of states. In equilibrium, the model is known to display a quantum phase transition from a generalized Kondo to a local moment phase. In the present work, we focus on the extension of these phases beyond equilibrium, i.e. under the influence of a bias voltage. Within the auxiliary master equation approach combined with a scheme based on matrix product states (MPS) we are able to directly address the current-carrying steady state. Starting with the equilibrium situation, we first corroborate our results by comparing with a direct numerical evaluation of ground state spectral properties of the system by MPS. Here, a scheme to locate the phase boundary by extrapolating the power-law exponent of the self energy produces a very good agreement with previous results obtained by the numerical renormalization group. Our nonequilibrium study as a function of the applied bias voltage is then carried out for two points on either side of the phase boundary. In the Kondo regime the resonance in the spectral function is splitted as a function of the increasing bias voltage. The local moment regime, instead, displays a dip in the spectrum near the position of the chemical potentials. Similar features are observed in the corresponding self energies. The Kondo split peaks approximately obey a power-law behavior as a function of frequency, whose exponents depend only slightly on voltage. Finally, the differential conductance in the Kondo regime shows a peculiar maximum at finite voltages, whose height, however, is below the accuracy level.

We present further developments of the auxiliary master equation approach (AMEA), a numerical method to simulate many-body quantum systems in as well as out of equilibrium, and apply it to the Interacting Resonant Level Model (IRLM) to benchmark the new developments. In particular, our results are obtained by employing the stochastic wave functions (SWF) method to solve the auxiliary open quantum system arising within AMEA. This development allows to reach extremely low walltimes for the calculation of correlation functions with respect to previous implementations of AMEA. An additional significant improvement is obtained by extrapolating a series of results obtained by increasing the number of auxiliary bath sites, NB, used within the auxiliary open quantum system formally to the limit of NB → ∞. Results for the current-voltage characteristics and for equilibrium correlation functions are compared with the one obtained by exact and matrix-product states based approaches.