Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userComparison of MPS based real time evolution algorithms for Anderson Impurity Models
71
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We perform a detailed comparison of two Matrix Product States (MPS) based time evolution algorithms for Anderson Impurity Models. To describe the bath, we use both the star-geometry as well as the commonly employed Wilson chain geometry. For each bath geometry, we use either the Time Dependent Variational Principle (TDVP) or the Time Evolving Block Decimation (TEBD) to perform the time evolution. To apply TEBD for the star-geometry, we use a specially adapted algorithm that can deal with the long-range coupling terms. Analyzing the major sources of errors, one expects them to be proportional to the system size for all algorithms. Surprisingly, we find errors independent of system size except for TEBD in chain geometry. Additionally, we show that the right combination of bath representation and time evolution algorithm is important. While TDVP in chain geometry is a very precise approach, TEBD in star geometry is much faster, such that for a given accuracy it is superior to TDVP in chain geometry. This makes the adapted version of TEBD in star geometry the most efficient method to solve impurity problems.
rate research
Read More
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 (tDMRG) 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.
We study the single impurity Anderson model (SIAM) using the equations of motion method (EOM), the non-crossing approximation (NCA), the one-crossing approximation (OCA), and Wilsons numerical renormalization group (NRG). We calculate the density of states and the linear conductance focusing on their dependence on the chemical potential and on the temperature paying special attention to the Kondo and Coulomb blockade regimes for a large range of model parameters. We report that some standard approximations based on the EOM technique display a rather unexpected poor behavior in the Coulomb blockade regime even at high temperatures. Our study offers a critical comparison between the different methods as well as a detailed compilation of the shortcomings and limitations due the approximations involved in each technique, thus allowing for a cost-benefit analysis of the different solvers that considers both numerical precision and computational performance.
We present a general theory for the Fano resonance in Anderson impurity systems. It is shown that the broadening of the impurity level leads to an additional and important contribution to the Fano resonance around the Fermi surface, especially in the mixed valence regime. This contribution results from the interference between the Kondo resonance and the broadened impurity level. Being applied to the scanning tunnelling microscopic experiments, we find that our theory gives a consistent and quantitative account for the Fano resonance lineshapes for both Co and Ti impurities on Au or Ag surfaces. The Ti systems are found to be in the mixed valence regime.
We introduce a block Lanczos (BL) recursive technique to construct quasi-one-dimensional models, suitable for density-matrix renormalization group (DMRG) calculations, from single- as well as multiple-impurity Anderson models in any spatial dimensions. This new scheme, named BL-DMRG method, allows us to calculate not only local but also spatially dependent static and dynamical quantities of the ground state for general Anderson impurity models without losing elaborate geometrical information of the lattice. We show that the BL-DMRG method can be easily extended to treat a multi-orbital Anderson impurity model. We also show that the symmetry adapted BL bases can be utilized, when it is appropriate, to reduce the computational cost. As a demonstration, we apply the BL-DMRG method to three different models for graphene: (i) a single adatom on the honeycomb lattice, (ii) a substitutional impurity in the honeycomb lattice, and (iii) an effective model for a single carbon vacancy in graphene. Our analysis reveals that, for the particle-hole symmetric case at half filling of electron density, the ground state of model (i) behaves as an isolated magnetic impurity with no Kondo screening while the ground state of the other two models forms a spin singlet state. We also calculate the real-space dependence of the spin-spin correlation functions between the impurity site and the conduction sites for these three models. Our results clearly show that, reflecting the presence of absence of unscreened magnetic moment at the impurity site, the spin-spin correlation functions decay as $r^{-3}$, differently from the non-interacting limit ($r^{-2}$), for model (i) and as $ r^{-4}$, exactly the same as the non-interacting limit, for models (ii) and (iii) in the asymptotic $r$, where $r$ is the distance between the impurity site and the conduction site.
We present a continuous-time Monte Carlo method for quantum impurity models, which combines a weak-coupling expansion with an auxiliary-field decomposition. The method is considerably more efficient than Hirsch-Fye and free of time discretization errors, and is particularly useful as impurity solver in large cluster dynamical mean field theory (DMFT) calculations.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
