We use the time dependent variational matrix product state (tVMPS) approach to investigate the dynamical properties of the single impurity Anderson model (SIAM). Under the Jordan-Wigner transformation, the SIAM is reformulated into two spin-1/2 XY chains with local magnetic fields along the z-axis. The chains are connected by the longitudinal Ising coupling at the end points. The ground state of the model is searched variationally within the space spanned by the matrix product state (MPS). The temporal Green's functions are calculated both by the imaginary and real time evolutions, from which the spectral information can be extracted. The possibility of using the tVMPS approach as an impurity solver for the dynamical mean field theory is also addressed. Finite temperature density operator is obtained by the ancilla approach. The results are compared to those from the Lanczos and the Hirsch-Fye quantum Monte-Carlo methods.
We use the time dependent variational matrix product state (tVMPS) approach to investigate the dynamical properties of the single impurity Anderson model (SIAM). Under the Jordan-Wigner transformation, the SIAM is reformulated into two spin-1/2 XY chains with local magnetic fields along the z-axis. The chains are connected by the longitudinal Ising coupling at the end points. The ground state of the model is searched variationally within the space spanned by the matrix product state (MPS). The temporal Green's functions are calculated both by the imaginary and real time evolutions, from which the spectral information can be extracted. The possibility of using the tVMPS approach as an impurity solver for the dynamical mean field theory is also addressed. Finite temperature density operator is obtained by the ancilla approach. The results are compared to those from the Lanczos and the Hirsch-Fye quantum Monte-Carlo methods.
An emergent and promising tensor-network-based impurity solver is to represent the Feynman-Vernon influence functional as a matrix product state, where the bath is integrated out analytically. Here we present an approach to calculate the equilibrium impurity spectral function based on the recently proposed Grassmann time-evolving matrix product operators method. The central idea is to perform a quench from a separable impurity-bath initial state as in the nonequilibrium scenario. The retarded Green's function $G(t+{t}_{0},{t}^{\ensuremath{'}}+{t}_{0})$ is then calculated after an equilibration time t0 such that the impurity and bath are approximately in thermal equilibrium. There are two major advantages of this method. First, since we focus on real-time dynamics, we do not need to perform the numerically ill-posed analytic continuation as in imaginary- time evolution-based methods. Second, the required bond dimension of the matrix product state in real-time calculations is observed to be much smaller than that in imaginary-time calculations, leading to a significant improvement in numerical efficiency. The accuracy of this method is demonstrated using the single-orbital Anderson impurity model and benchmarked against the continuous-time quantum Monte Carlo method.
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.
An essential role for the description of correlated quantum many-particle systems is played by impurity problems. These consist of a small number of localized orbitals, where the electrons are subject to a Coulomb-interaction, coupled to noninteracting bath degrees of freedom. Impurity models are used for the description of magnetic impurities in metals, of quantum points in nanostructures, and in the context of the dynamical mean field theory (DMFT). For their solution a large number of numerical methods exists. We mention here quantum Monte Carlo methods, the numerical renormalization group, or the exact diagonalization method. In this thesis we introduce a new approach to impurity problems, which is based on the functional renormalization group (fRG). In contrast to the methods mentioned above this approach is approximative, but it can be applied with comparatively lower numerical effort. The aim of this thesis is to examine the advantages and limits of this method. An important role in this connection is played by the single impurity Anderson model (SIAM), which is already understood very well with other methods and is seen as a prototype for models with strong local correlations. The spin fluctuations of this model are governed by the so-called Kondo energy scale and a central question of this thesis is whether this scale can be reproduced with our approach. We begin with a short introduction to the Green's function formalism for correlated quantum many particle systems. Then the fRG flow equations for the one-particle irreducible vertex functions are deduced, and several approximation schemes are discussed. In chapter 4 we introduce the fRG scheme used in this thesis. At first the SIAM is mapped to a semi-infinite chain, in which the interacting orbital is given by the first site of this chain. The system is then subdivided into two parts. The first part, which is called 'core', consists of the interacting orbital and the first L bath sites. The remaining bath sites form the second part. At the beginning both parts are decoupled from each other such that the core can be solved exactly. Starting with this exact solution the coupling between the core and the remaining bath sites is switched on slowly in a renormalization group flow. This way it is possible to calculate the one-particle and two-particle correlation functions on the correlated site. We call this flow scheme 'hybridization flow'. The flow equations are formulated in an effective theory on the first bath site outside the core, which turns out to be advantageous compared to other implementations of the hybridization flow. In chapter 5 this hybridization flow method is applied to the SIAM with semi-elliptic bath density of states. We discuss the differences which arise for different core-sizes (L = 0,1,2,3) and for different truncations of the fRG flow equations. In the cases L=0 and L=2 the local Fermi liquid properties of the SIAM and the associated Kondo scale can not be reproduced in both approximation schemes. For the core sizes L=1 and L=3 the ground state of the core is already a spin singlet state in the beginning of the flow and the Fermi liquid properties are reproduced. The Kondo energy scale, which determines the width of the central quasi-particle peak and the size of the local spin-susceptibility, is however not accurately resolved. This is partly due to the finite temperature used in our implementation such that for larger interaction strength the temperature is above the Kondo energy scale. The application of the hybridization flow method to quantum impurity problems in the context of the DMFT is discussed in chapter 6. In this approach lattice models with a local interaction, like the Hubbard model, are mapped to an effective impurity problem for a single lattice site embedded in a dynamical mean field representing the influence of the other electrons. This scheme is exact in infinitely many dimensions and a nonperturbative approximation method for finite-dimensional systems. In particular the Mott metal-insulator transition in the Hubbard model is successfully described by this method. One class of quantities, which is easily accessible from the hybridization flow scheme, is given by the local two-particle correlation functions. These play an essential role in non-local extensions of the DMFT. Therefore we mainly focus on these quantities. At first we apply our flow scheme to the Hubbard model on the Bethe lattice in infinite dimensions and calculate the local one-particle irreducible vertex functions in the insulating and the metallic phase in good agreement with previous calculations, which use the exact diagonalization method as impurity solver. After this we extend our flow scheme to a cluster DMFT method, which includes short-range antiferromagnetic correlations and calculate again the local and non-local vertex functions for the two-dimensional Hubbard model. We compare these to the results of the single-site DMFT.
We provide a detailed formulation of the recently proposed variational approach [Y. Ashida, T. Shi, M.-C. Ba\~nuls, J. I. Cirac, and E. Demler, Phys. Rev. Lett. 121, 026805 (2018)] to study ground-state properties and out-of-equilibrium dynamics for generic quantum spin-impurity systems. Motivated by the original ideas of Tomonaga, Lee, Low, and Pines, we construct a canonical transformation that completely decouples the impurity from the bath degrees of freedom. By combining this transformation with a Gaussian ansatz for the fermionic bath, we obtain a family of variational many-body states that can efficiently encode the strong entanglement between the impurity and fermions of the bath. We give a detailed derivation of equations of motions in the imaginary- and real-time evolutions on the variational manifold. We benchmark our approach by applying it to investigate ground-state and dynamical properties of the anisotropic Kondo model and compare results with those obtained using the matrix-product state (MPS) ansatz. We show that our approach can achieve an accuracy comparable to MPS-based methods with several orders of magnitude fewer variational parameters than the corresponding MPS ansatz. Comparisons to the Yosida ansatz and the exact solution from the Bethe ansatz are also discussed. We use our approach to investigate the two-lead Kondo model and analyze its long-time spatiotemporal behavior and the conductance behavior at finite bias and magnetic fields. The obtained results are consistent with the previous findings in the Anderson model and the exact solutions at the Toulouse point.
Abstract Tensor-network-based methods are promising candidates to solve quantum impurity problems (QIP). They are free of sampling noises and the sign problem compared to state-of-the-art continuous-time quantum Monte Carlo methods. Recent progress made in tensor-network-based impurity solvers is to use the Feynman–Vernon influence functional to integrate out the bath analytically, retaining only the impurity dynamics and representing it compactly as a matrix product state. The recently proposed Grassmann time-evolving matrix product operator (GTEMPO) method is one of the representative methods in this direction. In this work, we systematically study the performance of GTEMPO in solving equilibrium QIPs at a finite temperature with a semicircular spectrum density of the bath. Our results show that its computational cost would generally increase as the temperature goes down and scale exponentially with the number of orbitals. In particular, the single-orbital Anderson impurity model can be efficiently solved with this method, for two orbitals we estimate that one could possibly reach inverse temperature β ≈ 20 if high-performance computing techniques are utilized, while beyond that only very high-temperature regimes can be reached in the current formalism. Our work paves the way to apply GTEMPO as an imaginary-time impurity solver.
We present an infinite Grassmann time-evolving matrix product operator method for quantum impurity problems, which directly works in the steady state. The method embraces the well-established infinite matrix product state algorithms with the recently developed GTEMPO method, and benefits from both sides: it obtains real-time Green's functions without sampling noises and bath discretization error, it is applicable for any temperature without the sign problem, and its computational cost is independent of the transient dynamics and does not scale with the number of baths. We benchmark the method on the finite-temperature equilibrium Green's function in the noninteracting limit against exact solutions and in the single-orbital Anderson impurity model against GTEMPO calculations. We also study the zero-temperature nonequilibrium steady state of an impurity coupled to two baths with a voltage bias, obtaining consistent particle currents with existing calculations. The method is ideal for studying steady-state quantum transport, and can be readily used as an efficient real-time impurity solver in the dynamical mean-field theory and its nonequilibrium extension.