Historically, the RKKY or indirect exchange, interaction has been accepted as being able to be described by second order perturbation theory. A typical universal expression is usually given in this context. This approach, however, fails to incorporate many body ef- fects, quantum fluctuations, and other important details. In Chapter 2, a novel numerical approach is developed to tackle these problems in a quasi-exact, non-perturbative manner. Behind the method lies the main concept of being able to exactly map an n-dimensional lattice problem onto a 1-dimensional chain. The density matrix renormalization group algorithm is then employed to solve the newly cast Hamiltonian.
Historically, the RKKY or indirect exchange, interaction has been accepted as being able to be described by second order perturbation theory. A typical universal expression is usually given in this context. This approach, however, fails to incorporate many body ef- fects, quantum fluctuations, and other important details. In Chapter 2, a novel numerical approach is developed to tackle these problems in a quasi-exact, non-perturbative manner. Behind the method lies the main concept of being able to exactly map an n-dimensional lattice problem onto a 1-dimensional chain. The density matrix renormalization group algorithm is then employed to solve the newly cast Hamiltonian.
We review a completely unbiased and controlled numerical approach to quantum impurity problems consisting of exactly mapping the complete lattice Hamiltonian onto an equivalent one dimensional system through a unitary transformation. The resulting dimensional and entanglement reduction allows one to study single particle problems on arbitrary d-dimensional lattices using the density matrix renormalization group (DMRG) method. The real-space resolution allows one to position the impurities at the boundary or bulk of the sample and to study screening effects due to edge or surface modes. We describe how to generalize this approach to multi-impurity problems, discuss applications and possible extensions.
Strong correlations arising from interactions between particles are responsible for many phenomena in condensed matter physics, from superconductivity to heavy fermions. An accurate theoretical description of such complex systems requires to solve so-called many-body problems, which are inherently difficult to tackle due to the exponential growth of the underlying Fock space with the size of the system. Usually, single-particle descriptions or perturbation theory at first orders fail to provide satisfactory results in correlated regimes, and pertinent Ansätze or cutting edge numerical methods are required to solvethese problems.Among the wide variety of strongly correlated electronic systems, quantum impurity systems are of central interest. These systems, where a few localized degrees of freedom that undergo strong interactions hybridize with a larger system of free particles, have been widely studied for decades. Beyond their relevance to study dilute dynamic impurities in metals, quantum dots, or more complex nanostructures, they are exploited in a broader class of problems through dynamical mean field theory (DMFT), that is designed for fully interacting lattice problems. Despite the effort devoted to studying of these problems,some questions remain open and new approaches are required to address them.In this thesis, quantum impurity problems are tackled with the use of natural orbitals, a representation rooted in the framework of quantum chemistry. The numerical renormalization group, the standard non-perturbative method to treat these systems, is employed to demonstrate the hierarchical structure of correlations unveiled by natural orbitals. An ansatz for the wavefunction of the ground state is then proposed, separating explicitely a small set of correlated orbitals to a large uncorrelated one, leaving the problem to a few-body complexity. This simplification brought new insights for simulating quantumimpurity problems, and a new algorithm is developed to generate an optimized set of correlated orbitals independently of existing methods, going beyond their usual limitations.In the first place, the algorithm is used to compute the screening cloud in large one-dimensional chains and in a realistic 2D square lattice coupled to a localized adatom. This observable is of great interest to understand the screening mechanism of the magnetic moment of the impurity that plays a central role in transport properties, for unidimensional systems or in bulk metals experiencing RKKY interactions. Spatial correlations were found to be highly non-symmetric, breaking the s-wave picture. In one direction, corre-lations vanish way faster than the usual isotropic correlation length, while on the other one they spread further than this same expected length. This anisotropy suggests that impurities in bulk systems may interact differently depending on their relative positions with respect to the lattice.Taking further advantage of the efficiency of this new algorithm allowed us to perform a statistical study on large systems to investigate correlations around a quantum impurity in a disordered host, averaging over a large amount of realizations, which was not possible with previous approaches. The question motivating this study was the fate of the correlation cloud spreading into the disordered electronic bath. Although spatial correlations are expected to vanish when the disorder increases, long range correlationssurvive to some configurations of the potential, inducing a bimodal structure of probability distribution functions of relevant scales.Finally, these disordered problems are investigated in the non-interacting situation with random matrix theory, which reproduced qualitatively the results obtained before. The joint distribution of one-particle orbitals energies and amplitudes is calculated, which allows to compute any disordered averaged local correlation functions, and calculations in the large-N limit are presented.
A formulation of numerical real-space renormalization groups for quantum many-body problems is presented and several algorithms utilizing this formulation are outlined. The methods are presented and demonstrated using S=1/2 and S=1 Heisenberg chains as test cases. The key idea of the formulation is that rather than keep the lowest-lying eigenstates of the Hamiltonian in forming a new effective Hamiltonian of a block of sites, one should keep the most significant eigenstates of the block density matrix, obtained from diagonalizing the Hamiltonian of a larger section of the lattice which includes the block. This approach is much more accurate than the standard approach; for example, energies for the S=1 Heisenberg chain can be obtained to an accuracy of at least ${10}^{\mathrm{\ensuremath{-}}9}$. The method can be applied to almost any one-dimensional quantum lattice system, and can provide a wide variety of static properties.
In the past few years there has been a resurgence of interest in dynamical impurity problems, as a result of developments in the theory of correlated electron systems. The general dynamical impurity problem is a set of conduction electrons interacting with an impurity which has internal degrees of freedom. The simplest and earliest example, the Kondo problem, has attracted interest since the mid-sixties not only because of its physical importance but also as an example of a model displaying logarithmic divergences order by order in perturbation theory. It provided one of the earliest applications of the renormalization group method, which is designed to deal with just such a situation. As we shall see, the antiferromagnetic Kondo model is controlled by a strong-coupling fixed point, and the essence of the renormalization group solution is to carry out the global renormalization numerically starting from the original (weak-coupling) Hamiltonian. In these lectures, we shall describe an alternative route in which we identify an exactly solvable model which renormalizes to the same fixed point as the original dynamical impurity problem. This approach is akin to determining the critical behavior at a second order phase transition point by solving any model in a given universality class.
For the treatment of interacting electrons in crystal lattices approximations based on the picture of effective sites, coupled in a self-consistent fashion, have proven very useful. Particularly in the presence of strong local correlations, a local approach to the problem, combining a powerful method for the short ranged interactions with the lattice propagation part of the dynamics, determines the quality of results to a large extent. For a considerable time the non crossing approximation (NCA) in direct perturbation theory, an approach originally developed by Keiter for the Anderson impurity model, built a standard for the description of the local dynamics of interacting electrons. In the last couple of years exact methods like the numerical renormalization group (NRG) as pioneered by Wilson, have surpassed this approximation as regarding the description of the low energy regime. We present an improved approximation level of direct perturbation theory for finite Coulomb repulsion U, the crossing approximation one (CA1) and discuss its connections with other generalizations of NCA. CA1 incorporates all processes up to fourth order in the hybridization strength V in a self-consistent skeleton expansion, retaining the full energy dependence of the vertex functions. We reconstruct the local approach to the lattice problem from the point of view of cumulant perturbation theory in a very general way and discuss the proper use of impurity solvers for this purpose. Their reliability can be tested in applications to e.g. the Hubbard model and the Anderson-lattice model. We point out shortcomings of existing impurity solvers and improvements gained with CA1 in this context. This paper is dedicated to the memory of Hellmut Keiter.
This thesis contributes to the field of strongly correlated electron systems with studies in two distinct fields thereof: the specific nature of correlations between electrons in one dimension and quantum quenches in quantum impurity problems. In general, strongly correlated systems are characterized in that their physical behaviour needs to be described in terms of a many-body description, i.e. interactions correlate all particles in a complex way. The challenge is that the Hilbert space in a many-body theory is exponentially large in the number of particles. Thus, when no analytic solution is available - which is typically the case - it is necessary to find a way to somehow circumvent the problem of such huge Hilbert spaces. Therefore, the connection between the two studies comes from our numerical treatment: they are tackled by the density matrix renormalization group (DMRG) [1] and the numerical renormalization group (NRG) [2], respectively, both based on matrix product states. The first project presented in this thesis addresses the problem of numerically finding the dominant correlations in quantum lattice models in an unbiased way, i.e. without using prior knowledge of the model at hand. A useful concept for this task is the correlation density matrix (CDM) [3] which contains all correlations between two clusters of lattice sites. We show how to extract from the CDM, a survey of the relative strengths of the system’s correlations in different symmetry sectors as well as detailed information on the operators carrying long-range correlations and the spatial dependence of their correlation functions. We demonstrate this by a DMRG study of a one-dimensional spinless extended Hubbard model [4], while emphasizing that the proposed analysis of the CDM is not restricted to one dimension. The second project presented in this thesis is motivated by two phenomena under ongoing experimental and theoretical investigation in the context of quantum impurity models: optical absorption involving a Kondo exciton [5, 6, 7] and population switching in quantum dots [8, 9, 10, 11, 12, 13, 14, 15]. It turns out that both phenomena rely on the various manifestations of Anderson orthogonality (AO) [16], which describes the fact that the response of the Fermi sea to a quantum quench (i.e. an abrupt change of some property of the impurity or quantum dot) is a change of the scattering phase shifts of all the single-particle wave functions, therefore drastically changing the system. In this context, we demonstrate that NRG, a highly accurate method for quantum impurity models, allows for the calculation of all static and dynamic quantities related to AO and present an extensive NRG study for population switching in quantum dots.