اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدAuxiliary Hamiltonian representation of the nonequilibrium Dyson equation
177
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
The nonequilibrium Dyson (or Kadanoff-Baym) equation, which is an equation of motion with long-range memory kernel for real-time Green functions, underlies many numerical approaches based on the Keldysh formalism. In this paper we map the problem of solving the Dyson equation in real-time onto a noninteracting auxiliary Hamiltonian with additional bath degrees of freedom. The solution of the auxiliary model does not require the evaluation of a memory kernel and can thus be implemented in a very memory efficient way. The mapping is derived for a self-energy which is local in space and is thus directly applicable within nonequilibrium dynamical mean-field theory (DMFT). We apply the method to study the interaction quench in the Hubbard model for an optical lattice with a narrow confinement, using inhomogeneous DMFT in combination with second-order weak-coupling perturbation theory. We find that, although the quench excites pronounced density oscillations, signatures of the two-stage relaxation similar to the homogeneous system can be observed by looking at the time-dependent occupations of natural orbitals.
قيم البحث
اقرأ أيضاً
We study the single-impurity Anderson model out of equilibrium under the influence of a bias voltage $phi$ and a magnetic field $B$. We investigate the interplay between the shift ($omega_B$) of the Kondo peak in the spin-resolved density of states (
DOS) and the one ($phi_B$) 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 $phi_B$ only for $|g| mu_B B gg k_B T_K$. 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 present a general scheme to map correlated nonequilibrium quantum impurity problems onto an auxiliary open quantum system of small size. The infinite fermionic reservoirs of the original system are thereby replaced by a small number $N_B$ of nonin
teracting auxiliary bath sites whose dynamics is described by a Lindblad equation. Due to the presence of the intermediate bath sites, the overall dynamics acting on the impurity site is non-Markovian. With the help of an optimization scheme for the auxiliary Lindblad parameters, an accurate mapping is achieved, which becomes exponentially exact upon increasing $N_B$. The basic idea for this scheme was presented previously in the context of nonequilibrium dynamical mean field theory. In successive works on improved manybody solution strategies for the auxiliary Lindblad equation, such as Lanczos exact diagonalization or matrix product states, we applied the approach to study the nonequilibrium Kondo regime. In the present paper, we address in detail the mapping procedure itself, rather than the many-body solution. In particular, we investigate the effects of the geometry of the auxiliary system on the accuracy of the mapping for given $N_B$. Specifically, we present a detailed convergence study for five different geometries which, besides being of practical utility, reveals important insights into the underlying mechanisms of the mapping. For setups with onsite or nearest-neighbor Lindblad parameters we find that a representation adopting two separate bath chains is by far more accurate with respect to other choices based on a single chain or a commonly used star geometry. A significant improvement is obtained by allowing for long-ranged and complex Lindblad parameters. These results can be of great value when studying Lindblad-type approaches to correlated systems.
The recently introduced auxiliary Hamiltonian approach [Balzer K and Eckstein M 2014 Phys. Rev. B 89 035148] maps the problem of solving the two-time Kadanoff-Baym equations onto a noninteracting auxiliary system with additional bath degrees of freed
om. While the original paper restricts the discussion to spatially local self-energies, we show that there exists a rather straightforward generalization to treat also non-local correlation effects. The only drawback is the loss of time causality due to a combined singular value and eigen decomposition of the two-time self-energy, the application of which inhibits one to establish the self-consistency directly on the time step. For derivation and illustration of the method, we consider the Hubbard model in one dimension and study the decay of the Neel state in the weak-coupling regime, using the local and non-local second-order Born approximation.
We derive an exact mapping from the action of nonequilibrium dynamical mean-field theory (DMFT) to a single-impurity Anderson model (SIAM) with time-dependent parameters, which can be solved numerically by exact diagonalization. The representability
of the nonequilibrium DMFT action by a SIAM is established as a rather general property of nonequilibrium Green functions. We also obtain the nonequilibrium DMFT equations using the cavity method alone. We show how to numerically obtain the SIAM parameters using Cholesky or eigenvector matrix decompositions. As an application, we use a Krylov-based time propagation method to investigate the Hubbard model in which the hopping is switched on, starting from the atomic limit. Possible future developments are discussed.
Perturbation theory using self-consistent Greens functions is one of the most widely used approaches to study many-body effects in condensed matter. On the basis of general considerations and by performing analytical calculations for the specific exa
mple of the Hubbard atom, we discuss some key features of this approach. We show that when the domain of the functionals that are used to realize the map between the non-interacting and the interacting Greens functions is properly defined, there exists a class of self-energy functionals for which the self-consistent Dyson equation has only one solution, which is the physical one. We also show that manipulation of the perturbative expansion of the interacting Greens function may lead to a wrong self-energy as functional of the interacting Greens function, at least for some regions of the parameter space. These findings confirm and explain numerical results of Kozik et al. for the widely used skeleton series of Luttinger and Ward [Phys. Rev. Lett. 114, 156402]. Our study shows that it is important to distinguish between the maps between sets of functions and the functionals that realize those maps. We demonstrate that the self-consistent Greens functions approach itself is not problematic, whereas the functionals that are widely used may have a limited range of validity.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
