Search

We report an extensive high-field study of atacamite Cu$_2$Cl(OH)$_3$, a material realization of quantum sawtooth chains with weak interchain couplings, in continuous and pulsed magnetic fields up to 58 T. In particular, we have mapped the entropy landscape for fields as high as 35 T and have identified a field-induced quantum critical point at 21.9(1) T for $\mathbf{H} \parallel c$ axis. The quantum critical point separates field regions with and without magnetic order, evidenced by our thermodynamic study and $^1$H nuclear magnetic resonance spectroscopy, but lies far below full saturation of the magnetization. Corroborated by numerical results using density-matrix renormalization group (DMRG) calculations, we find this behavior associated with a decoupling of the sawtooth chain into two magnetic subsystems consisting of a spin-$1/2$ antiferromagnetic Heisenberg chain (backbone of the sawtooth chain) in the presence of an exchange field produced by the remaining field-polarized spins.

We document the ground state phase diagram of the one-dimensional Kitaev chain with quasi-periodic disorder in the presence of two-body interactions. Our data was obtained for systems of $L=1000$ sites using large-scale density matrix renormalization group numerics and is benchmarked against known results for the clean system. We demonstrate that moderate quasi-periodic disorder stabilizes the topological phase both for repulsive and attractive interactions. For larger disorder strengths, the system features re-entrance behavior and multiple phase transitions.

The antiferromagnetic sawtooth chain is a prototypical example of a frustrated spin system with vertex-sharing triangles, giving rise to complex quantum states. Depending on the interaction parameters, this system has three phases, of which the gapless non-collinear phase (for strongly coupled basal spins and loosely attached apical spins) has received little theoretical attention so far. In this work, we comprehensively investigate the properties of the non-collinear phase using large-scale tensor network computations which exploit the full SU(2) symmetry of the underlying Heisenberg model. We study the ground state both for finite systems using the density-matrix renormalization group (DMRG) as well as for infinite chains via the variational uniform matrix-product state (VUMPS) formalism. Finite temperatures and correlation functions are tackled via imaginary- or real time evolutions, which we implement using the time-dependent variational principle (TDVP). We find that the non-collinear phase is characterized by a double-Q structure for the apex-apex correlations. Deep into the phase, two peaks merge into a single one indicating a 90-degree spiral. The apical spins are soft and highly susceptible to external perturbations; they form a large number of gapless magnetic states that are polarized by weak fields and cause a long low-temperature tail in the specific heat. The dynamic spin-structure factor exhibits additive contributions from a two-spinon continuum (excitations of the basal chain) and a gapless peak at $k=\pi/2$ (excitations of the apical spins). Small temperatures excite the gapless states and smear the spectral weight of the $k=\pi/2$ peak out into a homogeneous flat-band structure. Our results are relevant, e.g., for the material atacamite Cu$_2$Cl(OH)$_3$ in high magnetic fields.

We study the effects of disorder in a one-dimensional model of $\mathbb{Z}_{3}$ Fock parafermions which can be viewed as a generalization of the prototypical Kitaev chain. Exact diagonalization is employed to determine level statistics, participation ratios, and the dynamics of domain walls. This allows us to identify ergodic as well as finite-size localized phases. In order to distinguish Anderson from many-body localization, we calculate the time evolution of the entanglement entropy in random initial states using tensor networks. We demonstrate that a purely quadratic parafermion model does not feature Anderson but many-body localization due to the nontrivial statistics of the particles.

We recapitulate recent developments of the functional renormalization group (FRG) approach to the steady state of systems out of thermal equilibrium. In particular, we discuss second-order truncation schemes which account for the frequency-dependence of the two particle vertex and which incorporate inelastic processes. Our focus is on two different types of one-dimensional fermion chains: i) infinite, open systems which feature a translation symmetry, and ii) finite systems coupled to left and right reservoirs. In addition to giving a detailed and unified review of the technical derivation of the FRG schemes, we briefly summarize some of the key physical results. In particular, we compute the non-equilibrium phase diagram and analyze the fate of the Berezinskii-Kosterlitz-Thouless transition in the infinite, open system., Comment: Appearing as part of a EPJB collection "Developments in the functional renormalization group approach to correlated electron systems"

The last decade has witnessed an impressive progress in the theoretical understanding of transport properties of clean, one-dimensional quantum lattice systems. Many physically relevant models in one dimension are Bethe-ansatz integrable, including the anisotropic spin-1/2 Heisenberg (also called spin-1/2 XXZ chain) and the Fermi-Hubbard model. Nevertheless, practical computations of, for instance, correlation functions and transport coefficients pose hard problems from both the conceptual and technical point of view. Only due to recent progress in the theory of integrable systems on the one hand and due to the development of numerical methods on the other hand has it become possible to compute their finite temperature and nonequilibrium transport properties quantitatively. Most importantly, due to the discovery of a novel class of quasilocal conserved quantities, there is now a qualitative understanding of the origin of ballistic finite-temperature transport, and even diffusive or super-diffusive subleading corrections, in integrable lattice models. We shall review the current understanding of transport in one-dimensional lattice models, in particular, in the paradigmatic example of the spin-1/2 XXZ and Fermi-Hubbard models, and we elaborate on state-of-the-art theoretical methods, including both analytical and computational approaches. Among other novel techniques, we discuss matrix-product-states based simulation methods, dynamical typicality, and, in particular, generalized hydrodynamics. We will discuss the close and fruitful connection between theoretical models and recent experiments, with examples from both the realm of quantum magnets and ultracold quantum gases in optical lattices., Comment: Review article (78 pages, 34 figures). V3: Revised version, additional references

Many-body localized systems in which interactions and disorder come together defy the expectations of quantum statistical mechanics: In contrast to ergodic systems, they do not thermalize when undergoing nonequilibrium dynamics. What is less clear, however, is how topological features interplay with many-body localized phases as well as the nature of the transition between a topological and a trivial state within the latter. In this work, we numerically address these questions, using a combination of extensive tensor network calculations, specifically DMRG-X, as well as exact diagonalization, leading to a comprehensive characterization of Hamiltonian spectra and eigenstate entanglement properties., Comment: 7 pages, 8 figures, submitted to PRB

We introduce the Loschmidt amplitude as a powerful tool to perform spectroscopy of generic many-body wave functions and use it to interrogate the wave function obtained after ramping the transverse field quantum Ising model through its quantum critical point. Previous results are confirmed and a more complete understanding of the population of defects and of the effects of magnon-magnon interaction or finite-size corrections is obtained. The influence of quantum coherence is clarified., Comment: 4.5 pages + SM, 5 + 3 figures, accepted version

We study the dynamics arising from a double quantum quench where the parameters of a given Hamiltonian are abruptly changed from being in an equilibrium phase A to a different phase B and back (A$\to$B$\to$A). As prototype models, we consider the (integrable) transverse field Ising as well as the (non-integrable) ANNNI model. The return amplitude features non-analyticities after the first quench through the equilibrium quantum critical point (A$\to$B), which is routinely taken as a signature of passing through a so-called dynamical quantum phase transition. We demonstrate that non-analyticities after the second quench (B$\to$A) can be avoided and reestablished in a recurring manner upon increasing the time $T$ spent in phase B. The system retains an infinite memory of its past state, and one has the intriguing opportunity to control at will whether or not dynamical quantum phase transitions appear after the second quench., Comment: 6 pages, 5 figures

We study the finite-energy density phase diagram of spinless fermions with attractive interactions in one dimension in the presence of uncorrelated diagonal disorder. Unlike the case of repulsive interactions, a delocalized Luttinger-liquid phase persists at weak disorder in the ground state, which is a well-known result. We revisit the ground-state phase diagram and show that the recently introduced occupation-spectrum discontinuity computed from the eigenspectrum of one-particle density matrices is noticeably smaller in the Luttinger liquid compared to the localized regions. Moreover, we use the functional renormalization scheme to study the finite-size dependence of the conductance, which resolves the existence of the Luttinger liquid as well and is computationally cheap. Our main results concern the finite-energy density case. Using exact diagonalization and by computing various established measures of the many-body localization-delocalization transition, we argue that the zero-temperature Luttinger liquid smoothly evolves into a finite-energy density ergodic phase without any intermediate phase transition.

We analyze the dynamics of the return amplitude following a sudden quench in the three-state quantum Potts chain. For quenches crossing the quantum critical point from the paramagnetic to the ferromagnetic phase, the corresponding rate function is non-analytic at critical times and behaves linearly in their vicinity. In particular, we find no indication of a link between the time evolution close to the critical times and the scaling properties of the quantum critical point in the Potts chain., Comment: submitted to PRB

One of the most straightforward ways to study thermal properties beyond linear response is to monitor the relaxation of an arbitrarily large left-right temperature gradient $T_L-T_R$. In one-dimensional systems which support ballistic thermal transport, the local energy currents $\langle j(t)\rangle$ acquire a non-zero value at long times, and it was recently investigated whether or not this steady state fulfills a simple additive relation $\langle j(t\to\infty)\rangle=f(T_L)-f(T_R)$ in integrable models. In this paper, we probe the non-equilibrium dynamics of the Hubbard chain using density matrix renormalization group (DMRG) numerics. We show that the above form provides an effective description of thermal transport in this model; violations are below the finite-time accuracy of the DMRG. As a second setup, we study how an initially equilibrated system radiates into different non-thermal states (such as the vacuum).

Integrable models such as the spin-1/2 Heisenberg chain, the Lieb-Liniger or the one-dimensional Hubbard model are known to avoid thermalization, which was also demonstrated in several quantum-quench experiments. Another dramatic consequence of integrability is the zero-frequency anomaly in transport coefficients, which results in ballistic finite-temperature transport, despite the presence of strong interactions. While this aspect of nonergodic dynamics has been known for a long time, there has so far not been any unambiguous experimental realization thereof. We make a concrete proposal for the observation ballistic transport via local quantum quench experiments in fermionic quantum-gas microscopes. Such an experiment would also unveil the coexistence of ballistic and diffusive transport channels in one and the same system and provide a means of measuring finite-temperature Drude weights. The connection between local quenches and linear-response functions is established via time-dependent Einstein relations.

We illustrate how finite-temperature charge and thermal Drude weights of one-dimensional systems can be obtained from the relaxation of initial states featuring global (left-right) gradients in the chemical potential or temperature. The approach is tested for spinless interacting fermions as well as for the Fermi-Hubbard model, and the behaviour in the vicinity of singular points (such as half filling or isotropic chains) is discussed. We present technical details on how to implement the calculation in practice using the density matrix renormalization group and show that the non-equilibrium dynamics is often less demanding to simulate numerically and features simpler finite-time transients than the corresponding linear response current correlators; thus, new parameter regimes can become accessible. As an application, we determine the thermal Drude weight of the Hubbard model for temperatures T which are an order of magnitude smaller than those reached in the equilibrium approach. This allows us to demonstrate that at low T and half filling, thermal transport is successively governed by spin excitations and described quantitatively by the Bethe ansatz Drude weight of the Heisenberg chain.

We study the expectation values of observables and correlation functions at long times after a global quantum quench. Our focus is on metallic (`gapless') fermionic many-body models and small quenches. The system is prepared in an eigenstate of an initial Hamiltonian, and the time evolution is performed with a final Hamiltonian which differs from the initial one in the value of one global parameter. We first derive general relations between time-averaged expectation values of observables as well as correlation functions and those obtained in an eigenstate of the final Hamiltonian. Our results are valid to linear and quadratic order in the quench parameter g and generalize prior insights in several essential ways. This allows us to develop a phenomenology for the thermalization of local quantities up to a given order in g. Our phenomenology is put to a test in several case studies of one-dimensional models representative of four distinct classes of Hamiltonians: quadratic ones, effectively quadratic ones, those characterized by an extensive set of (quasi-) local integrals of motion, and those for which no such set is known (and believed to be nonexistent). We show that for each of these models, all observables and correlation functions thermalize to linear order in g. The more local a given quantity, the longer the linear behavior prevails when increasing g. Typical local correlation functions and observables for which the term O(g) vanishes thermalize even to order g^2. Our results show that lowest order thermalization of local observables is an ubiquitous phenomenon even in models with extensive sets of integrals of motion., Comment: version as published

We study the properties of excited states in one-dimensional many-body localized (MBL) systems using a matrix product state algorithm. First, the method is tested for a large disordered non-interacting system, where for comparison we compute a quasi-exact reference solution via a Monte Carlo sampling of the single-particle levels. Thereafter, we present extensive data obtained for large interacting systems of L~100 sites and large bond dimensions chi~1700, which allows us to quantitatively analyze the scaling behavior of the entanglement S in the system. The MBL phase is characterized by a logarithmic growth (L)~log(L) over a large scale separating the regimes where volume and area laws hold. We check the validity of the eigenstate thermalization hypothesis. Our results are consistent with the existence of a mobility edge.

We study the thermal conductivity of the one-dimensional Fermi-Hubbard model at finite temperature using a density matrix renormalization group approach. The integrability of this model gives rise to ballistic thermal transport. We calculate the temperature dependence of the thermal Drude weight at half filling for various interactions and moreover, we compute its filling dependence at infinite temperature. The finite-frequency contributions originating from the fact that the energy current is not a conserved quantity are investigated as well. We report evidence that breaking the integrability through a nearest-neighbor interaction leads to vanishing Drude weights and diffusive energy transport. Moreover, we demonstrate that energy spreads ballistically in local quenches with initially inhomogeneous energy density profiles in the integrable case. We discuss the relevance of our results for thermalization in ultra-cold quantum gas experiments and for transport measurements with quasi-one dimensional materials.

We study the interplay of interactions and disorder in a one-dimensional fermion lattice coupled adiabatically to infinite reservoirs. We employ both the functional renormalization group (FRG) as well as matrix product state techniques, which serve as an accurate benchmark for small systems. Using the FRG, we compute the length- and temperature-dependence of the conductance averaged over $10^4$ samples for lattices as large as $10^{5}$ sites. We identify regimes in which non-ohmic power law behavior can be observed and demonstrate that the corresponding exponents can be understood by adapting earlier predictions obtained perturbatively for disordered Luttinger liquids. In presence of both disorder and isolated impurities, the conductance has a universal single-parameter scaling form. This lays the groundwork for an application of the functional renormalization group to the realm of many-body localization.

We study the spin and thermal conductivity of spin-1/2 ladders at finite temperature. This is relevant for experiments with quantum magnets. Using a state-of-the-art density matrix renormalization group algorithm, we compute the current autocorrelation functions on the real-time axis and then carry out a Fourier integral to extract the frequency dependence of the corresponding conductivities. The finite-time error is analyzed carefully. We first investigate the limiting case of spin-1/2 XXZ chains, for which our analysis suggests non-zero dc-conductivities in all interacting cases irrespective of the presence or absence of spin Drude weights. For ladders, we observe that all models studied are normal conductors with no ballistic contribution. Nonetheless, only the high-temperature spin conductivity of XX ladders has a simple diffusive, Drude-like form, while Heisenberg ladders exhibit a more complicated low-frequency behavior. We compute the dc spin conductivity down to temperatures of the order of T~0.5J, where J is the exchange coupling along the legs of the ladder. We further extract mean-free paths and discuss our results in relation to thermal conductivity measurements on quantum magnets.

We generalize nonlinear Luttinger liquid theory to describe the dynamics of one-dimensional quantum critical systems at low temperatures. Analyzing density-matrix renormalization group results for the spin autocorrelation function in the XXZ chain we provide, in particular, direct evidence for spin diffusion in sharp contrast to the exponential decay in time predicted by conventional Luttinger liquid theory. Furthermore, we discuss how the frequencies and exponents of the oscillatory contributions from the band edges are renormalized by irrelevant interactions and obtain excellent agreement between our finite temperature nonlinear Luttinger liquid theory and the numerical data.

We study finite-temperature transport properties of the one-dimensional Hubbard model using the density matrix renormalization group. Our aim is two-fold: First, we compute both the charge and the spin current correlation function of the integrable model at half filling. The former decays rapidly, implying that the corresponding Drude weight is either zero or very small. Second, we calculate the optical charge conductivity sigma(omega) in presence of small integrability-breaking next-nearest neighbor interactions (the extended Hubbard model). The DC conductivity is finite and diverges as the temperature is decreased below the gap. Our results thus suggest that the half-filled, gapped Hubbard model is a normal charge conductor at finite temperatures. As a testbed for our numerics, we compute sigma(omega) for the integrable XXZ spin chain in its gapped phase.

We discuss a few simple modifications to time-dependent density matrix renormalization group (DMRG) algorithms which allow to access larger time scales. We specifically aim at beginners and present practical aspects of how to implement these modifications almost effortlessly within any standard matrix product state (MPS) based formulation of the method. Most importantly, we show how to 'combine' the Schroedinger and Heisenberg time evolutions of arbitrary pure states |psi> and operators A in the evaluation of _psi(t)=. This includes quantum quenches. The generalization (non-)thermal mixed state dynamics _rho(t)=Tr[rho A(t)] induced by an initial density matrix rho is straightforward. In the context of equilibrium (ground state or finite temperature T>0) correlation functions, one can extend the simulation time by a factor of two by 'exploiting time translation invariance', which is efficiently implementable within MPS DMRG. We present a simple analytic argument for why a recently-introduced disentangler succeeds in reducing the effort of time-dependent simulations at T>0. Finally, we advocate the python programming language as an elegant option for beginners to set up a DMRG code.

We study the spin- and energy dynamics in one-dimensional spin-1/2 systems induced by local quantum quenches at finite temperatures using a time-dependent density matrix renormalization group method. System sizes are chosen large enough to ensure that the time-dependent data for the accesible time scales represents the behavior in the thermodynamic limit. As a main result, we observe a ballistic spreading of perturbations of the energy density in the integrable spin-1/2 XXZ chain for all temperatures and exchange anisotropies, related to the divergent thermal conductivity in this model and the exact conservation of the energy current. In contrast, the spin dynamics is ballistic in the massless phase, but shows a diffusive behavior at high temperatures in the easy-axis phase in the case of a vanishing background spin density. We extract a quantitative estimate for the spin diffusion constant from the time dependence of the spatial variance of the spin density, which agrees well with values obtained from current-current correlation functions using an Einstein relation. As an example for non-integrable models, we consider two-leg ladders, for which we observe indications of diffusive energy and spin dynamics. The relevance of our results for recent experiments with quantum magnets and bosons in optical lattices is discussed.

Finite-temperature transport properties of one-dimensional systems can be studied using the time dependent density matrix renormalization group via the introduction of auxiliary degrees of freedom which purify the thermal statistical operator. We demonstrate how the numerical effort of such calculations is reduced when the physical time evolution is augmented by an additional time evolution within the auxiliary Hilbert space. Specifically, we explore a variety of integrable and non-integrable, gapless and gapped models at temperatures ranging from T=infty down to T/bandwidth=0.05 and study both (i) linear response where (heat and charge) transport coefficients are determined by the current-current correlation function and (ii) non-equilibrium driven by arbitrary large temperature gradients. The modified DMRG algorithm removes an 'artificial' build-up of entanglement between the auxiliary and physical degrees of freedom. Thus, longer time scales can be reached.

We investigate the dynamics following sudden quenches across quantum critical points belonging to different universality classes. Specifically, we use matrix product state methods to study the quantum Ising chain in the presence of two additional terms which break integrability. We find that in all models the rate function for the return probability to the initial state becomes a non-analytic function of time in the thermodynamic limit. This so-called `dynamical phase transition' was first observed in a recent work by Heyl, Polkovnikov, and Kehrein [Phys. Rev. Lett. 110, 135704 (2013)] for the exactly-solvable quantum Ising chain, which can be mapped to free fermions. Our results for `interacting theories' indicate that non-analytic dynamics is a generic feature of sudden quenches across quantum critical points. We discuss potential connections to the dynamics of the order parameter.

We revisit the problem of the spin Drude weight D of the integrable spin-1/2 XXZ chain using two complementrary approaches, exact diagonalization (ED) and the time-dependent density-matrix renormalization group (tDMRG). We pursue two main goals. First, we present extensive results for the temperature dependence of D. By exploiting time translation invariance within tDMRG, one can extract D for significantly lower temperatures than in previous tDMRG studies. Second, we discuss the numerical quality of the tDMRG data and elaborate on details of the finite-size scaling of the ED results, comparing calculations carried out in the canonical and grand-canonical ensembles. Furthermore, we analyze the behavior of the Drude weight as the point with SU(2)-symmetric exchange is approached and discuss the relative contribution of the Drude weight to the sum rule as a function of temperature., Comment: Revised version, additional ED and DMRG data for negative Delta, DMRG data for lower temperatures, discussion of Benz et al. 2005 corrected

Many low-dimensional materials are well described by integrable one-dimensional models such as the Hubbard model of electrons or the Heisenberg model of spins. However, the small perturbations to these models required to describe real materials are expected to have singular effects on transport quantities: integrable models often support dissipationless transport, while weak non-integrable terms lead to finite conductivities. We use matrix-product-state methods to obtain quantitative values of spin/electrical and thermal conductivities in an almost integrable gapless chain (an XXZ spin chain with staggered fields, or equivalently a spinless fermion chain with staggered on-site potentials). The results at low temperatures validate a scaling theory based on bosonization.

We study the real-time dynamics of spin chains driven out of thermal equilibrium by an initial temperature gradient T_L \neq T_R using density matrix renormalization group methods. We demonstrate that the nonequilibrium energy current saturates fast to a finite value if the linear-response thermal conductivity is infinite, i.e. if the Drude weight D is nonzero. Our data suggests that a nonintegrable dimerized chain might support such dissipationless transport (D>0). We show that the steady-state value J_E of the current for arbitrary T_L \neq T_R is of the functional form J_E=f(T_L)-f(T_R), i.e. it is completely determined by the linear conductance. We argue for this functional form, which is essentially a Stefan-Boltzmann law in this integrable model; for the XXX ferromagnet, f can be computed via thermodynamic Bethe ansatz in good agreement with the numerics. Inhomogeneous systems exhibiting different bulk parameters as well as Luttinger liquid boundary physics induced by single impurities are discussed briefly.

We study one-dimensional spinless fermions at zero and finite temperature T using the density matrix renormalization group. We consider nearest as well as next-nearest neighbor interactions; the latter render the system inaccessible by a Bethe ansatz treatment. Using an infinite-system alogrithm we demonstrate the emergence of Luttinger liquid physics at low energies for a variety of static correlation functions as well as for thermodynamic properties. The characteristic power law suppression of the momentum distribution n(k) function at T=0 can be directly observed over several orders of magnitude. At finite temperature, we show that n(k) obeys a scaling relation. The Luttinger liquid parameter and the renormalized Fermi velocity can be extracted from the density response function, the specific heat, and/or the susceptibility without the need to carry out any finite-size analysis. We illustrate that the energy scale below which Luttinger liquid power laws manifest vanishes as the half-filled system is driven into a gapped phase by large interactions.

We provide strong evidence that the relaxation dynamics of one-dimensional, metallic Fermi systems resulting out of an abrupt amplitude change of the two-particle interaction has aspects which are universal in the Luttinger liquid sense: The leading long-time behavior of certain observables is described by universal functions of the equilibrium Luttinger liquid parameter and the renormalized velocity. We analytically derive those functions for the Tomonaga-Luttinger model and verify our hypothesis of universality by considering spinless lattice fermions within the framework of the density matrix renormalization group.

We introduce a real time version of the functional renormalization group which allows to study correlation effects on nonequilibrium transport through quantum dots. Our method is equally capable to address (i) the relaxation out of a nonequilibrium initial state into a (potentially) steady state driven by a bias voltage and (ii) the dynamics governed by an explicitly time-dependent Hamiltonian. All time regimes from transient to asymptotic can be tackled; the only approximation is the consistent truncation of the flow equations at a given order. As an application we investigate the relaxation dynamics of the interacting resonant level model which describes a fermionic quantum dot dominated by charge fluctuations. Moreover, we study decoherence and relaxation phenomena within the ohmic spin-boson model by mapping the latter to the interacting resonant level model.

We propose an easily implemented approach to study time-dependent correlation functions of one dimensional systems at finite temperature T using the density matrix renormalization group. The entanglement growth inherent to any time-dependent calculation is significantly reduced if the auxiliary degrees of freedom which purify the statistical operator are time evolved with the physical Hamiltonian but reversed time. We exploit this to investigate the long time behavior of current correlation functions of the XXZ spin-1/2 Heisenberg chain. This allows a direct extraction of the Drude weight D at intermediate to large T. We find that D is nonzero -- and thus transport is dissipationless -- everywhere in the gapless phase. At low temperatures we establish an upper bound to D by comparing with bosonization.

We investigate two serially-aligned quantum dots in the molecular regime of large tunnel couplings t. A Zeeman field B is used to tune the energy difference of singlet and triplet spin configurations. Attaching this geometry to BCS source and drain leads with gap Delta and phase difference phi gives rise to an equilibrium supercurrent J. To compute J in presence of Coulomb interactions U between the dot electrons, we employ the functional renormalization group (FRG). For B\simt -- where the singlet and lowest-lying triplet spin states are equal in energy -- the current exhibits characteristics of a 0-pi transition similar to a single impurity. Its magnitude in the pi phase, however, jumps discontinuously at B=t, being smaller on the triplet side. Exploiting the flexibility of the FRG, we demonstrate that this effect is generic and calculate J for realistic experimental parameters Delta, U, and gate voltages epsilon. To obtain a more thorough understanding of the discontinuity, we analytically treat the limit Delta=infty where one can access the exact many-particle states. Finally, carrying out perturbation theory in the dot-lead couplings substantiates the intuitive picture that Cooper pair tunneling is favored by a singlet spin configuration while inhibited by a triplet one.

Using low-temperature scanning tunneling spectroscopy applied to the Cs-induced two-dimensional electron system (2DES) on p-type InSb(110), we probe electron-electron interaction effects in the quantum Hall regime. The 2DES is decoupled from p-doped bulk states and exhibits spreading resistance within the insulating quantum Hall phases. In quantitative agreement with calculations we find an exchange enhancement of the spin splitting. Moreover, we observe that both the spatially averaged as well as the local density of states feature a characteristic Coulomb gap at the Fermi level. These results show that electron-electron interaction effects can be probed down to a resolution below all relevant length scales., Comment: supplementary movie in ancillary files

We study transport properties of quantum impurity systems using the functional renormalization group. The latter is an RG-based diagrammatic tool to treat Coulomb interactions in a fast and flexible way. Prior applications, which employed a simple first-order (Hartree-Fock-like) scheme to truncate the FRG flow equations within the Matsubara formalism, succeeded in accurately describing linear transport of various quantum dot geometries at zero temperature T=0. In a nutshell, advance in this Thesis is three-fold. First, we introduce a frequency-dependent second-order approximation in order to eventually compute finite-energy properties such as the conductance at T>0 (mainly focusing on the single impurity Anderson model). Second, a generalisation of the Hartree-Fock-like approach to Keldysh space allows for addressing the non-equilibrium steady-state dynamics of the interacting resonant level model. Third, we investigate the physics of a quantum dot Josephson junction as well as the charging of a single narrow level using the first-order scheme., Comment: Summary only - for the complete Thesis please go to http://www.theorie.physik.uni-goettingen.de/~karrasch/publications/thesis_karrasch.pdf

We illustrate how to calculate the finite-temperature linear-response conductance of quantum impurity models from the Matsubara Green function. A continued fraction expansion of the Fermi distribution is employed which was recently introduced by Ozaki [Phys. Rev. B 75, 035123 (2007)] and converges much faster than the usual Matsubara representation. We give a simplified derivation of Ozaki's idea using concepts from many-body condensed matter theory and present results for the rate of convergence. In case that the Green function of some model of interest is only known numerically, interpolating between Matsubara frequencies is much more stable than carrying out an analytic continuation to the real axis. We demonstrate this explicitly by considering an infinite tight-binding chain with a single site impurity as an exactly-solvable test system, showing that it is advantageous to calculate transport properties directly on the imaginary axis. The formalism is applied to the single impurity Anderson model, and the linear conductance at finite temperatures is calculated reliably at small to intermediate Coulomb interactions by virtue of the Matsubara functional renormalization group. Thus, this quantum many-body method combined with the continued fraction expansion of the Fermi function constitutes a promising tool to address more complex quantum dot geometries at finite temperatures., Comment: version accepted by Phys. Rev. B

We study the steady-state current in a minimal model for a quantum dot dominated by charge fluctuations and analytically describe the time evolution into this state. The current is driven by a finite bias voltage V across the dot, and two different renormalization group methods are used to treat small to intermediate local Coulomb interactions. The corresponding flow equations can be solved analytically which allows to identify all microscopic cutoff scales. Exploring the entire parameter space we find rich non-equilibrium physics which cannot be understood by simply considering the bias voltage as an infrared cutoff. For the experimentally relevant case of left-right asymmetric couplings, the current generically shows a power-law suppression for large V. The relaxation dynamics towards the steady state features characteristic oscillations as well as an interplay of exponential and power-law decay.

We investigate equilibrium and steady-state non-equilibrium transport properties of a spinless resonant level locally coupled to two conduction bands of width ~\Gamma via a Coulomb interaction U and a hybridization t'. In order to study the effects of finite bias voltages beyond linear response, a generalization of the functional renormalization group to Keldysh frequency space is employed. Being mostly unexplored in the context of quantum impurity systems out of equilibrium, we benchmark this method against recently-published time-dependent density matrix renormalization group data. We thoroughly investigate the scaling limit \Gamma\to\infty characterized by the appearance of power laws. Most importantly, at the particle-hole symmetric point the steady-state current decays like J ~ V^{-\alpha_J} as a function of the bias voltage V>>t', with an exponent \alpha_J(U) that we calculate to leading order in the Coulomb interaction strength. In contrast, we do not observe a pure power-law (but more complex) current-voltage-relation if the energy \epsilon of the resonant level is pinned close to either one of the chemical potentials \pm V/2.

We study a quantum dot Josephson junction inside an Aharonov-Bohm environment. The geometry is modeled by an Anderson impurity coupled to two directly-linked BCS leads. We illustrate that the well-established picture of the low-energy physics being governed by an interplay of two distinct (singlet and doublet) phases is still valid for this interferometric setup. The phase boundary depends, however, non-monotonically on the coupling strength between the superconductors, causing the system to exhibit re-entrance behavior and multiple phase transitions. We compute the zero-temperature Josephson current and demonstrate that it can become negative in the singlet phase by virtue of the Coulomb interaction U. As a starting point, the limit of large superconducting energy gaps \Delta=\infty is solved analytically. In order to tackle arbitrary \Delta<\infty and U>0, we employ a truncated functional renormalization group scheme which was previously demonstrated to give quantitatively reliable results for the quantum dot Josephson problem.

Understanding the charging of exceptionally narrow levels in quantum dots in the presence of interactions remains a challenge within mesoscopic physics. We address this fundamental question in the generic model of a narrow level capacitively coupled to a broad one. Using bosonization we show that for arbitrary capacitive coupling charging can be described by an analogy to the magnetization in the anisotropic Kondo model, featuring a low-energy crossover scale that depends in a power-law fashion on the tunneling amplitude to the level. Explicit analytical expressions for the exponent are derived and confirmed by detailed numerical and functional renormalization-group calculations., Comment: 4+ pages, 3 figures, PRL version

We investigate the Josephson current in a single wall carbon nanotube connected to superconducting electrodes. We focus on the parameter regime in which transport is dominated by Kondo physics. A sizeable supercurrent is observed for odd number of electrons on the nanotube when the Kondo temperature Tk is sufficiently large compared to the superconducting gap. On the other hand when, in the center of the Kondo ridge, Tk is slightly smaller than the superconducting gap, the supercurrent is found to be extremely sensitive to the gate voltage Vbg. Whereas it is largely suppressed at the center of the ridge, it shows a sharp increase at a finite value of Vbg. This increase can be attributed to a doublet-singlet transition of the spin state of the nanotube island leading to a pi shift in the current phase relation. This transition is very sensitive to the asymmetry of the contacts and is in good agreement with theoretical predictions., Comment: 5 pages, 4 figures

We use the Matsubara functional renormalization group (FRG) to describe electronic correlations within the single impurity Anderson model. In contrast to standard FRG calculations, we account for the frequency-dependence of the two-particle vertex in order to address finite-energy properties (e.g, spectral functions). By comparing with data obtained from the numerical renormalization group (NRG) framework, the FRG approximation is shown to work well for arbitrary parameters (particularly finite temperatures) provided that the electron-electron interaction U is not too large. We demonstrate that aspects of (large U) Kondo physics which are described well by a simpler frequency-independent truncation scheme are no longer captured by the 'higher-order' frequency-dependent approximation. In contrast, at small to intermediate U the results obtained by the more elaborate scheme agree better with NRG data. We suggest to parametrize the two-particle vertex not by three independent energy variables but by introducing three functions each of a single frequency. This considerably reduces the numerical effort to integrate the FRG flow equations., Comment: accepted by J. Phys.: Condensed Matter

We investigate the Josephson current J(\phi) through a quantum dot embedded between two superconductors showing a phase difference \phi. The system is modeled as a single Anderson impurity coupled to BCS leads, and the functional and the numerical renormalization group frameworks are employed to treat the local Coulomb interaction U. We reestablish the picture of a quantum phase transition occurring if the ratio between the Kondo temperature T_K and the superconducting energy gap \Delta or, at appropriate T_K/\Delta, the phase difference \phi or the impurity energy is varied. We present accurate zero- as well as finite-temperature T data for the current itself, thereby settling a dispute raised about its magnitude. For small to intermediate U and at T=0 the truncated functional renormalization group is demonstrated to produce reliable results without the need to implement demanding numerics. It thus provides a tool to extract characteristics from experimental current-voltage measurements., Comment: version accepted for publication in PRB

We investigate the appearance of pi lapses in the transmission phase theta of a two-level quantum dot with Coulomb interaction U. Using the numerical and functional renormalization group methods we study the entire parameter space for spin-polarized as well as spin-degenerate dots, modeled by spinless or spinful electrons, respectively. We investigate the effect of finite temperatures T. For small T and sufficiently small single-particle spacings delta of the dot levels we find pi phase lapses between two transmission peaks in an overwhelming part of the parameter space of the level-lead couplings. For large delta the appearance or not of a phase lapse between resonances depends on the relative sign of the level-lead couplings in analogy to the U=0 case. We show that this generic scenario is the same for spin-polarized and spin-degenerate dots. We emphasize that in contrast to dots with more levels, for a two-level dot with small delta and generic dot-lead couplings (that is up to cases with special symmetry) the "universal" phase lapse behavior is already established at U=0. The most important effect of the Coulomb interaction is to increase the separation of the transmission resonances. The relation of the appearance of phase lapses to the inversion of the population of the dot levels is discussed. For the spin-polarized case and low temperatures we compare our results to recent mean-field studies. For small delta correlations are found to strongly alter the mean-field picture., Comment: submitted to NJP

We present a recently-developed renormalization group scheme, the functional renormalization group (fRG), as a many-particle method suited to account for the two-particle interactions between the electrons in complex quantum dot geometries. A detailed derivation of a truncated set of RG flow equations that requires only basic knowledge of the functional integral approach to many-particle physics is given. Using fRG we study linear-response transport properties of a variety of quantum dots in the zero-temperature limit. For parallel geometries, we observe a generic crossover from mesoscopic to universal behaviour of the transmission phase if the single-particle level spacing is gradually decreased. We tackle a couple of well-known systems (the SIAM, the side-coupled geometry and short Hubbard chains), showing that the fRG correctly captures Kondo physics. We make use of the little computational resources required to scan the whole parameter space of these systems, partly uncovering new physics. The high accuracy of the fRG is demonstrated by extensive comparisons with NRG data., Comment: Diploma thesis, 143 pages; for printer-friendly ps version see http://www.theorie.physik.uni-goettingen.de/~karrasch

Transmission phase \alpha measurements of many-electron quantum dots (small mean level spacing \delta) revealed universal phase lapses by \pi between consecutive resonances. In contrast, for dots with only a few electrons (large \delta), the appearance or not of a phase lapse depends on the dot parameters. We show that a model of a multi-level quantum dot with local Coulomb interactions and arbitrary level-lead couplings reproduces the generic features of the observed behavior. The universal behavior of \alpha for small \delta follows from Fano-type antiresonances of the renormalized single-particle levels., Comment: 4 pages, version accepted for publication in PRL

We investigate the effect of local Coulomb correlations on electronic transport through a variety of coupled quantum dot systems connected to Fermi liquid leads. We use a newly developed functional renormalization group scheme to compute the gate voltage dependence of the linear conductance, the transmission phase, and the dot occupancies. A detailed derivation of the flow equations for the dot level positions, the inter-dot hybridizations, and the effective interaction is presented. For specific setups and parameter sets we compare the results to existing accurate numerical renormalization group data. This shows that our approach covers the essential physics and is quantitatively correct up to fairly large Coulomb interactions while being much faster, very flexible, and simple to implement. We then demonstrate the power of our method to uncover interesting new physics. In several dots coupled in series the combined effect of correlations and asymmetry leads to a vanishing of transmission resonances. In contrast, for a parallel double-dot we find parameter regimes in which the two-particle interaction generates additional resonances., Comment: 18 pages, 18 figures included, version accepted for publication in PRB