Auxiliary master equation approach within stochastic wave functions: Application to the interacting resonant level model.

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 to benchmark the new developments. In particular, our results are obtained by employing the stochastic wave functions method to solve the auxiliary open quantum system arising within AMEA. This development allows us to reach extremely low wall times 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, N_{B}, used within the auxiliary open quantum system formally to the limit of N_{B}→∞. 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. Further, we complement this benchmark by the presentation of spectral functions for higher temperatures where we find different behaviors around zero frequency depending on the hybridization strength.

A WENO-solver combined with adaptive momentum discretization for the Wigner transport equation and its application to resonant tunneling diodes.

We present a novel numerical scheme for the deterministic solution of the Wigner transport equation, especially suited to deal with situations in which strong quantum effects are present. The unique feature of the algorithm is the expansion of the Wigner function in local basis functions, similar to finite element or finite volume methods. This procedure yields a discretization of the pseudo-differential operator that conserves the particle density on arbitrarily chosen grids. The high flexibility in refining the grid spacing together with the weighted essentially non-oscillatory (WENO) scheme for the advection term allows for an accurate and well-resolved simulation of the phase space dynamics. A resonant tunneling diode is considered as test case and a detailed convergence study is given by comparing the results to a non-equilibrium Green's functions calculation. The impact of the considered domain size and of the grid spacing is analyzed. The obtained convergence of the results towards a quasi-exact agreement of the steady state Wigner and Green's functions computations demonstrates the accuracy of the scheme, as well as the high flexibility to adjust to different physical situations.