The decay rate of quasiparticles in quantum dots is studied through the real time calculation of the single-particle Green function in the self-consistent approximation. The method avoids exact diagonalization, transforming the problem into a system of coupled non-linear integral equations which may be solved iteratively. That allows us to study systems larger than previously treated in the literature. Our results for the inverse participation ratio of many-body states show that the threshold energy for the quasiparticle disintegration is ${E}^{*}\ensuremath{\sim}\sqrt{g}\ensuremath{\Delta}.$ The delocalization transition is soft rather than sharp. Three different regimes as function of the effective interaction strength may be clearly identified at high energies.
The decay rate of quasiparticles in quantum dots is studied through the real time calculation of the single-particle Green function in the self-consistent approximation. The method avoids exact diagonalization, transforming the problem into a system of coupled non-linear integral equations which may be solved iteratively. That allows us to study systems larger than previously treated in the literature. Our results for the inverse participation ratio of many-body states show that the threshold energy for the quasiparticle disintegration is E∗\ensuremath∼g\ensuremathΔ. The delocalization transition is soft rather than sharp. Three different regimes as function of the effective interaction strength may be clearly identified at high energies.
In interacting quantum systems, the single-particle Green's function is expected to decay in time due to the interaction induced decoherence of quasiparticles. In the limit of weak interaction strengths (Δ), a naive application of Fermi's Golden Rule (FGR) predicts an O(Δ2) quasiparticle decay rate. However, for 1d fermions on the lattice at T>0, this calculation gives a divergent result and the scaling of the quasiparticle lifetime with interaction strength remains an open question. In this work we propose a solution to this question: combining numerical simulations using the recently introduced dissipation-assisted operator evolution (DAOE) method, with non-perturbative diagrammatic re-summations, we predict a logarithmic enhancement of the quasiparticle decay rate τ−1∼Δ2logΔ−2. We argue that this effect is present in a wide variety of well-known weakly interacting quantum fermionic and bosonic systems, and even in some classical systems, provided the non-interacting limit has quasiparticles with a generic dispersion.
We study the fermionic quasiparticle spectrum in a zero-temperature superfluid Fermi gas, and in particular how it is modified by different disintegration processes. On top of the disintegration by emission of a collective boson ( 1\to2 1 → 2 , subject of a previous study, PRL , 073404), we consider here disintegration events where three quasiparticles are emitted ( 1\to3 1 → 3 ). We show that both disintegration processes are described by a t t -matrix self-energy (as well as some highly off-resonant vacuum processes), and we characterize the associated disintegration continua. At strong coupling, we show that the quasiparticle spectrum is heavily distorted near the 1\to3 1 → 3 disintegration threshold. Near the dispersion minimum, where the quasiparticles remain well-defined, the main effect of the off-shell disintegration processes is to shift the location of the minimum by a value that corresponds to the Hartree shift in the BCS limit. With our approximation of the self-energy, the correction to the energy gap with respect to the mean-field result however remains small, in contrast with experimental measurements.
We develop a general approach to the nonequilibrium dynamics of quantum-impurity systems for arbitrary coupling strength. The numerical renormalization group is used to generate a complete basis set necessary for the correct description of the time evolution. We benchmark our method with the exact analytical solution for the resonant-level model. As a first application, we investigate the equilibration of an ultrasmall quantum dot subject to a sudden change of gate voltage and external magnetic field. Two distinct relaxation times are identified for the spin and charge dynamics.
The question of which non-interacting Green's function 'best' describes an interacting many-body electronic system is both of fundamental interest as well as of practical importance in describing electronic properties of materials in a realistic manner. Here, we study this question within the framework of Baym-Kadanoff theory, an approach where one locates the stationary point of a total energy functional of the one-particle Green's function in order to find the total ground-state energy as well as all one-particle properties such as the density matrix, chemical potential, or the quasiparticle energy spectrum and quasiparticle wave functions. For the case of the Klein functional, our basic finding is that minimizing the length of the gradient of the total energy functional over non-interacting Green's functions yields a set of self-consistent equations for quasiparticles that is identical to those of the quasiparticle self-consistent GW (QSGW) (van Schilfgaarde et al 2006 Phys. Rev. Lett. 96 226402-4) approach, thereby providing an a priori justification for such an approach to electronic structure calculations. In fact, this result is general, applies to any self-energy operator, and is not restricted to any particular approximation, e.g., the GW approximation for the self-energy. The approach also shows that, when working in the basis of quasiparticle states, solving the diagonal part of the self-consistent Dyson equation is of primary importance while the off-diagonals are of secondary importance, a common observation in the electronic structure literature of self-energy calculations. Finally, numerical tests and analytical arguments show that when the Dyson equation produces multiple quasiparticle solutions corresponding to a single non-interacting state, minimizing the length of the gradient translates into choosing the solution with largest quasiparticle weight.
In this thesis, I analyze properties of mesoscopic systems that exhibit strong quantum correlations. To this purpose, I adapt and use the density-matrix renormalization group (DMRG), a numerical method that has been developed specifically for analyzing such systems. The following systems are discussed: - Josephson effect between superconducting nanograins: I investigate the Josephson effect, i.e. the dependence of the energy of two weakly coupled superconductors on the difference between their superconducting phase, in the regime that the level spacing d is comparable to the bulk superconducting gap Delta. Because the BCS solution is inapplicable in this regime, I use the DMRG to calculate the ground state of the two coupled superconductors and extract the Josephson energy E_J. The standard BCS result for E_J is reproduced only in the continuum limit of vanishing d. As d is increased, however, E_J turns out to display a nonmonotonic behaviour. A tight-binding approximation for weak Josephson coupling explains the physical mechanism underlying this reentrance in a transparent way. - Well-defined quasiparticles in small metallic grains: I analyze zero-temperature spectral functions of mesoscopic systems such as quantum dots and metallic grains, in the limit of large conductance, where they can be described by a "universal Hamiltonian" model. I show that within this model, an important class of spectral functions is dominated by one single energy eigenstate only. For an interacting system this is a very peculiar property, which implies an infinite lifetime of the quasiparticles. Moreover, I show that the dominating eigenstate contains only a limited subclass of all excitations, which I characterize as the "No-Gaudino" excitations. Hence, these are sufficient to explain many of the properties of the systems under consideration. Besides its own physical significance, the dominance of the "No-Gaudino" states has also high practical value, because it permits the calculation of zero-temperature spectral functions with high accuracy using the DMRG. I illustrate the use of this method by calculating the tunneling density of states of metallic grains and the magnetic response of mesoscopic rings. - Real-time dynamics in spin-1/2 chains I study the nonequilibrium transport properties in spin-1/2 chains by solving the time evolution of a non-stationary initial state, and investigate the influence of different interaction strength and dimerization on the magnetization transport. To this purpose, I use the "adaptive time-dependent DMRG", a DMRG variant that was recently developed to solve the many-body Schroedinger equation with high accuracy. I find that the magnetization possesses a well-defined long-time limit, whose nature does not depend on the dimerization, but only on the strength J_z of the S^z-S^z-interaction: For |J_z| < 1 I find ballistic magnetization transport, and for J_z >1 almost no transport, with a sharp crossover at |J^z|=1. I explain this crossover as a subtle consequence of a quantum phase transition which occurs at the precise value J_z| = 1. - Many-body scattering states: I present a general method for calculating many-body scattering states that does not rely on the assumptions of perturbation theory or near-equilibrium. Due to the complexity of the problem, I limit the discussion to a description of the algorithm itself and a few proof-of-principle calculations only. The strategy is to calculate the many-body scattering state that results when a scatterer (e.g. a quantum dot) is connected to two leads, to which a voltage bias V is applied. This state is obtained by solving the many-body Lippmann-Schwinger equation using the DMRG.