Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:0807.0581 · Physical Review B

Transport properties and Kondo correlations in nanostructures: Time-dependent DMRG method applied to quantum dots coupled to Wilson chains

Luis G. G. V. Dias da Silva, Fabian Heidrich‐Meisner, Adrian E. Feiguin, C. A. Büsser, George Balster Martins, E. V. Anda, Elbio R. Dagotto

Published Nov 19, 2008• 88 citations
arXivHTMLPDF

Abstract

We apply the adaptive time-dependent density-matrix renormalization-group method (tDMRG) to the study of transport properties of quantum-dot systems connected to metallic leads. Finite-size effects make the usual tDMRG description of the Kondo regime a numerically demanding task. We show that such effects can be attenuated by describing the leads by ``Wilson chains,'' in which the hopping matrix elements decay exponentially away from the impurity $({t}_{n}\ensuremath{\propto}{\ensuremath{\Lambda}}^{\ensuremath{-}n/2})$. For a given system size and in the linear-response regime, results for $\ensuremath{\Lambda}>1$ show several improvements over the undamped $\ensuremath{\Lambda}=1$ case: perfect conductance is obtained deeper in the strongly interacting regime and current plateaus remain well defined for longer time scales. Similar improvements were obtained in the finite-bias regime up to bias voltages of the order of the Kondo temperature. These results show that with the proposed modification, the tDMRG characterization of Kondo correlations in the transport properties can be substantially improved, while it turns out to be sufficient to work with much smaller system sizes. We discuss the numerical cost of this approach with respect to the necessary system sizes and the entanglement growth during the time evolution.

DOI
10.1103/physrevb.78.195317

Related papers

Adaptive time-dependent density-matrix renormalization-group technique for calculating the conductance of strongly correlated nanostructures

K. A. Al-Hassanieh, Adrian E. Feiguin, José A. Riera +2 more

A procedure based on the recently developed ``adaptive'' time-dependent density-matrix-renormalization-group (DMRG) technique is presented to calculate the zero temperature conductance of nanostructures, such as quantum dots (QDs) or molecular conductors, when represented by a small number of active levels. The leads are modeled using noninteracting tight-binding Hamiltonians. The ground state at time zero is calculated at zero bias. Then, a small bias is applied between the two leads, the wave function is DMRG evolved in time, and currents are measured as a function of time. Typically, the current is expected to present periodicities over long times, involving intermediate well-defined plateaus that resemble steady states. The conductance can be obtained from those steady-state-like currents. To test this approach, several cases of interacting and noninteracting systems have been studied. Our results show excellent agreement with exact results in the noninteracting case. More importantly, the technique also reproduces quantitatively well-established results for the conductance and local density of states in both the cases of one and two coupled interacting QDs. The technique also works at finite bias voltages, and it can be extended to include interactions in the leads.

arXiv:cond-mat/0601411•Physical Review B•May 3, 2006•87 citations

Conductance of correlated systems: real-time dynamics in finite systems

Alexander Branschädel, Guenter Schneider, Peter Schmitteckert

Numerical time evolution of transport states using time dependent Density Matrix Renormalization Group (td-DMRG) methods has turned out to be a powerful tool to calculate the linear and finite bias conductance of interacting impurity systems coupled to non-interacting one-dimensional leads. Several models, including the Interacting Resonant Level Model (IRLM), the Single Impurity Anderson Model (SIAM), as well as models with different multi site structures, have been subject of investigations in this context. In this work we give an overview of the different numerical approaches that have been successfully applied to the problem and go into considerable detail when we comment on the techniques that have been used to obtain the full I--V-characteristics for the IRLM. Using a model of spinless fermions consisting of an extended interacting nanostructure attached to non-interacting leads, we explain the method we use to obtain the current--voltage characteristics and discuss the finite size effects that have to be taken into account. We report results for the linear and finite bias conductance through a seven site structure with weak and strong nearest-neighbor interactions. Comparison with exact diagonalisation results in the non-interacting limit serve as a verification of the accuracy of our approach. Finally we discuss the possibility of effectively enlarging the finite system by applying damped boundaries and give an estimate of the effective system size and accuracy that can be expected in this case.

arXiv:1004.4178•Apr 23, 2010•0 citations

Kondo effect out of equilibrium

Thomas Korb, Herbert Schoeller

Correlations are a central ingredient of modern condensed matter physics. With the advent of mesoscopic physics, it became possible to study correlated systems in non-equilibrium situations, which can be reached by e.g. driving a constant current through a quantum dot. A basic question is how correlations on the one hand and a non-equilibrium current on the other hand influence each other. Here this interplay is studied for the model system of a quantum dot in the Kondo regime using the real time renormalization group (RTRG) method. This method is a functional renormalization group (RG) method allowing a non-perturbative description of correlated systems in non-equilibrium. For the Kondo quantum dot the method is generalized to the case of double vertices. The RTRG method not only provides equations for the calculation of the renormalized coupling between the quantum dot spin and the reservoir electrons, but also equations for the calculation of the non-equilibrium decoherence rate and the non-equilibrium current.The RG equations are derived to second order. This approximation is expected to be valid as long as the quantum dot does not enter the strong coupling regime, i.e. the effective coupling strength between the quantum dot spin and the electrons does not become too large. It is found that this is the case if the bias voltage applied to the dot or the temperature is much larger than the Kondo temperature.It is shown that the RG equations lead to coupling functions depending on the energy of the incident and the scattered electron. This dependence is simplified to the dependence on the average energy of the incident and the scattered electron in the weak coupling regime.A closer inspection of the RG equations for the couplings shows that they not only contain terms leading to a logarithmic increase of the coupling strength but also subleading terms which are unimportant in the weak coupling regime. The RG equations are solved analytically for the case where the subleading terms are neglected. It is found that the coupling functions develop resonances at energies associated with the chemical potentials of the reservoirs. The height of these resonances is determined by the decoherence rate. However, for bias voltages much larger than the Kondo temperature these resonances do not affect the conductance for a two terminal setup, so the conductance does not provide a measure for the decoherence rate in this case. It is found that the conductance for the two terminal setup is a universal function of voltage/Kondo temperature and temperature/Kondo temperature in the weak coupling regime. This function is calculated analytically in the limits of vanishing voltage or temperature. Furthermore, it is found that for zero temperature the decoherence rate is proportional to the non-equilibrium current flowing through the system. This supports the physical picture that a non-equilibrium current is associated with a decoherence rate, which cuts the logarithmic scaling of the coupling strength.Finally, the RG equations are solved numerically for zero temperature including the subleading corrections found before. As expected, the subleading corrections have no impact on the decoherence rate and the conductance in the weak coupling regime. On the other hand, it is found that they lead to strong deviations from the results of the analytic leading logarithmic approximation if the bias voltage approaches the Kondo temperature. Generally, it is found that the subleading corrections lead to a saturation of the decoherence rate and the conductance to some finite value. However, for the numerical solution a numerical regularization parameter has to be introduced. It turns out that the saturation value found depends on its value, so no statements on the precise value of the decoherence rate and the current can be made in this regime.

•Jan 1, 2006•0 citations

Real-time simulations of nonequilibrium transport in the single-impurity Anderson model

Fabian Heidrich‐Meisner, Adrian E. Feiguin, Elbio R. Dagotto

One of the main open problems in the field of transport in strongly interacting nanostructures is the understanding of currents beyond the linear response regime. In this work, we consider the single-impurity Anderson model and use the adaptive time-dependent density matrix renormalization group method to compute real-time currents out of equilibrium. We first focus on the particle-hole symmetric point where Kondo correlations are the strongest and then extend the study of the nonequilibrium transport to the mixed-valence regime. As a main result, we present accurate data for the current-voltage characteristics of this model.

arXiv:0903.2414•Physical Review B•Jun 29, 2009•179 citations

Open Wilson chain numerical renormalization group approach to steady-state nonequilibrium quantum transport

Anand Manaparambil, Frithjof B. Anders

The numerical renormalization group (NRG) approach was developed to identify and quantify different equilibrium regimes of quantum impurity systems (QISs) with unprecedented accuracy by a tailored finite size representation. Out of equilibrium, the steady-state density operator is not of the Boltzmannian form but one that is determined by the imposed boundary conditions. We extend the NRG to the nonequilibrium setting by augmenting each Wilson site with a reservoir, whose coupling functions are calculated via a continuous fraction expansion in order to recover the continuum limit exactly. The nonequilibrium parameters such as a finite bias as well as a finite temperature gradient enters through the Bloch-Redfield tensor (BRT), whose zero eigenvector gives the steady-state density operator. We used the resulting open chain full density matrix (OC-FDM) approach with an effective single lead description to investigate the charge and spin transport through a quantum dot (QD) under finite bias and temperature gradient. The influence of lead asymmetry and an external magnetic field on the transport properties are also studied. For completeness, we have also investigated local properties of the QD, such as charge fluctuations and find excellent agreement with real-time quantum Monte Carlo (RT-QMC) data. The OC-FDM approach was able to explore Kondo energy scales as low as $T_K/D \approx 10^{-8}$ in the non-equilibrium regime, as well as show convergence with the established equilibrium benchmarks, such as a quantitative agreement with full density matrix numerical renormalization group (FDM-NRG) and Fermi-liquid scaling at small bias and temperature. Owing to the effective single lead description, a single OC-FDM data point takes orders of magnitude less time on a standard laptop, compared to other state-of-the-art numerical methods.

arXiv:2606.14635•Physical review. B./Physical review. B•Oct 23, 2026•0 citations

Time-Dependent DMRG Study on Quantum Dot under a Finite Bias Voltage

Shunsuke Kirino, Tatsuya Fujii, Jize Zhao +1 more

Resonant tunneling through quantum dot under a finite bias voltage at zero temperature is investigated by using the adaptive time-dependent density matrix renormalization group(TdDMRG) method. Quantum dot is modeled by the Anderson Hamiltonian with the 1-D nearest-neighbor tight-binding leads. Initially the ground state wave function is calculated with the usual DMRG method. Then the time evolution of the wave function due to the slowly changing bias voltage between the two leads is calculated by using the TdDMRG technique. Even though the system size is finite, the expectation values of current operator show steady-like behavior for a finite time interval, in which the system is expected to resemble the real nonequilibrium steady state of the infinitely long system. We show that from the time intervals one can obtain quantitatively correct results for differential conductance in a wide range of bias voltage. Finally we observe an anomalous behavior in the expectation value of the double occupation operator at the dot $$ as a function of bias voltage.

arXiv:0805.0218•Journal of the Physical Society of Japan•Jul 25, 2008•48 citations
We apply the adaptive time-dependent density-matrix renormalization-group method (tDMRG) to the study of transport properties of quantum-dot systems connected to metallic leads. Finite-size effects make the usual tDMRG description of the Kondo regime a numerically demanding task. We show that such effects can be attenuated by describing the leads by ``Wilson chains,'' in which the hopping matrix elements decay exponentially away from the impurity (tn\ensuremath∝\ensuremathΛ\ensuremath−n/2)({t}_{n}\ensuremath{\propto}{\ensuremath{\Lambda}}^{\ensuremath{-}n/2})(tn​\ensuremath∝\ensuremathΛ\ensuremath−n/2). For a given system size and in the linear-response regime, results for \ensuremathΛ>1\ensuremath{\Lambda}>1\ensuremathΛ>1 show several improvements over the undamped \ensuremathΛ=1\ensuremath{\Lambda}=1\ensuremathΛ=1 case: perfect conductance is obtained deeper in the strongly interacting regime and current plateaus remain well defined for longer time scales. Similar improvements were obtained in the finite-bias regime up to bias voltages of the order of the Kondo temperature. These results show that with the proposed modification, the tDMRG characterization of Kondo correlations in the transport properties can be substantially improved, while it turns out to be sufficient to work with much smaller system sizes. We discuss the numerical cost of this approach with respect to the necessary system sizes and the entanglement growth during the time evolution.
A procedure based on the recently developed ``adaptive'' time-dependent density-matrix-renormalization-group (DMRG) technique is presented to calculate the zero temperature conductance of nanostructures, such as quantum dots (QDs) or molecular conductors, when represented by a small number of active levels. The leads are modeled using noninteracting tight-binding Hamiltonians. The ground state at time zero is calculated at zero bias. Then, a small bias is applied between the two leads, the wave function is DMRG evolved in time, and currents are measured as a function of time. Typically, the current is expected to present periodicities over long times, involving intermediate well-defined plateaus that resemble steady states. The conductance can be obtained from those steady-state-like currents. To test this approach, several cases of interacting and noninteracting systems have been studied. Our results show excellent agreement with exact results in the noninteracting case. More importantly, the technique also reproduces quantitatively well-established results for the conductance and local density of states in both the cases of one and two coupled interacting QDs. The technique also works at finite bias voltages, and it can be extended to include interactions in the leads.
Numerical time evolution of transport states using time dependent Density Matrix Renormalization Group (td-DMRG) methods has turned out to be a powerful tool to calculate the linear and finite bias conductance of interacting impurity systems coupled to non-interacting one-dimensional leads. Several models, including the Interacting Resonant Level Model (IRLM), the Single Impurity Anderson Model (SIAM), as well as models with different multi site structures, have been subject of investigations in this context. In this work we give an overview of the different numerical approaches that have been successfully applied to the problem and go into considerable detail when we comment on the techniques that have been used to obtain the full I--V-characteristics for the IRLM. Using a model of spinless fermions consisting of an extended interacting nanostructure attached to non-interacting leads, we explain the method we use to obtain the current--voltage characteristics and discuss the finite size effects that have to be taken into account. We report results for the linear and finite bias conductance through a seven site structure with weak and strong nearest-neighbor interactions. Comparison with exact diagonalisation results in the non-interacting limit serve as a verification of the accuracy of our approach. Finally we discuss the possibility of effectively enlarging the finite system by applying damped boundaries and give an estimate of the effective system size and accuracy that can be expected in this case.
Correlations are a central ingredient of modern condensed matter physics. With the advent of mesoscopic physics, it became possible to study correlated systems in non-equilibrium situations, which can be reached by e.g. driving a constant current through a quantum dot. A basic question is how correlations on the one hand and a non-equilibrium current on the other hand influence each other. Here this interplay is studied for the model system of a quantum dot in the Kondo regime using the real time renormalization group (RTRG) method. This method is a functional renormalization group (RG) method allowing a non-perturbative description of correlated systems in non-equilibrium. For the Kondo quantum dot the method is generalized to the case of double vertices. The RTRG method not only provides equations for the calculation of the renormalized coupling between the quantum dot spin and the reservoir electrons, but also equations for the calculation of the non-equilibrium decoherence rate and the non-equilibrium current.The RG equations are derived to second order. This approximation is expected to be valid as long as the quantum dot does not enter the strong coupling regime, i.e. the effective coupling strength between the quantum dot spin and the electrons does not become too large. It is found that this is the case if the bias voltage applied to the dot or the temperature is much larger than the Kondo temperature.It is shown that the RG equations lead to coupling functions depending on the energy of the incident and the scattered electron. This dependence is simplified to the dependence on the average energy of the incident and the scattered electron in the weak coupling regime.A closer inspection of the RG equations for the couplings shows that they not only contain terms leading to a logarithmic increase of the coupling strength but also subleading terms which are unimportant in the weak coupling regime. The RG equations are solved analytically for the case where the subleading terms are neglected. It is found that the coupling functions develop resonances at energies associated with the chemical potentials of the reservoirs. The height of these resonances is determined by the decoherence rate. However, for bias voltages much larger than the Kondo temperature these resonances do not affect the conductance for a two terminal setup, so the conductance does not provide a measure for the decoherence rate in this case. It is found that the conductance for the two terminal setup is a universal function of voltage/Kondo temperature and temperature/Kondo temperature in the weak coupling regime. This function is calculated analytically in the limits of vanishing voltage or temperature. Furthermore, it is found that for zero temperature the decoherence rate is proportional to the non-equilibrium current flowing through the system. This supports the physical picture that a non-equilibrium current is associated with a decoherence rate, which cuts the logarithmic scaling of the coupling strength.Finally, the RG equations are solved numerically for zero temperature including the subleading corrections found before. As expected, the subleading corrections have no impact on the decoherence rate and the conductance in the weak coupling regime. On the other hand, it is found that they lead to strong deviations from the results of the analytic leading logarithmic approximation if the bias voltage approaches the Kondo temperature. Generally, it is found that the subleading corrections lead to a saturation of the decoherence rate and the conductance to some finite value. However, for the numerical solution a numerical regularization parameter has to be introduced. It turns out that the saturation value found depends on its value, so no statements on the precise value of the decoherence rate and the current can be made in this regime.
One of the main open problems in the field of transport in strongly interacting nanostructures is the understanding of currents beyond the linear response regime. In this work, we consider the single-impurity Anderson model and use the adaptive time-dependent density matrix renormalization group method to compute real-time currents out of equilibrium. We first focus on the particle-hole symmetric point where Kondo correlations are the strongest and then extend the study of the nonequilibrium transport to the mixed-valence regime. As a main result, we present accurate data for the current-voltage characteristics of this model.
The numerical renormalization group (NRG) approach was developed to identify and quantify different equilibrium regimes of quantum impurity systems (QISs) with unprecedented accuracy by a tailored finite size representation. Out of equilibrium, the steady-state density operator is not of the Boltzmannian form but one that is determined by the imposed boundary conditions. We extend the NRG to the nonequilibrium setting by augmenting each Wilson site with a reservoir, whose coupling functions are calculated via a continuous fraction expansion in order to recover the continuum limit exactly. The nonequilibrium parameters such as a finite bias as well as a finite temperature gradient enters through the Bloch-Redfield tensor (BRT), whose zero eigenvector gives the steady-state density operator. We used the resulting open chain full density matrix (OC-FDM) approach with an effective single lead description to investigate the charge and spin transport through a quantum dot (QD) under finite bias and temperature gradient. The influence of lead asymmetry and an external magnetic field on the transport properties are also studied. For completeness, we have also investigated local properties of the QD, such as charge fluctuations and find excellent agreement with real-time quantum Monte Carlo (RT-QMC) data. The OC-FDM approach was able to explore Kondo energy scales as low as TK/D≈10−8T_K/D \approx 10^{-8}TK​/D≈10−8 in the non-equilibrium regime, as well as show convergence with the established equilibrium benchmarks, such as a quantitative agreement with full density matrix numerical renormalization group (FDM-NRG) and Fermi-liquid scaling at small bias and temperature. Owing to the effective single lead description, a single OC-FDM data point takes orders of magnitude less time on a standard laptop, compared to other state-of-the-art numerical methods.
Resonant tunneling through quantum dot under a finite bias voltage at zero temperature is investigated by using the adaptive time-dependent density matrix renormalization group(TdDMRG) method. Quantum dot is modeled by the Anderson Hamiltonian with the 1-D nearest-neighbor tight-binding leads. Initially the ground state wave function is calculated with the usual DMRG method. Then the time evolution of the wave function due to the slowly changing bias voltage between the two leads is calculated by using the TdDMRG technique. Even though the system size is finite, the expectation values of current operator show steady-like behavior for a finite time interval, in which the system is expected to resemble the real nonequilibrium steady state of the infinitely long system. We show that from the time intervals one can obtain quantitatively correct results for differential conductance in a wide range of bias voltage. Finally we observe an anomalous behavior in the expectation value of the double occupation operator at the dot $$ as a function of bias voltage.