Source author record

Jan von Delft

Jan von Delft appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

50works
14topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

50 published item(s)

preprint2022arXiv

Benchmark Calculations of Multiloop Pseudofermion fRG

The pseudofermion functional renormalization group (pffRG) is a computational method for determining zero-temperature phase diagrams of frustrated quantum magnets. In a recent methodological advance, the commonly employed Katanin truncation of the flow equations was extended to include multiloop corrections, thereby capturing additional contributions from the three-particle vertex [ arXiv:2011.01268 , arXiv:2011.01269 ]. This development has also stimulated significant progress in the numerical implementation of pffRG, allowing one to track the evolution of pseudofermion vertices under the renormalization group flow with unprecedented accuracy. However, cutting-edge solvers differ in their integration algorithms, heuristics to discretize Matsubara frequency grids, and more. To lend confidence in the numerical robustness of state-of-the-art multiloop pffRG codes, we present and compare results produced with two independently developed and algorithmically distinct solvers for Heisenberg models on three-dimensional lattice geometries. Using the cubic lattice Heisenberg (anti)ferromagnet with nearest and next-nearest neighbor interactions as a generic benchmark model, we find the two codes to quantitatively agree, often up to several orders of magnitude in digital precision, both on the level of spin-spin correlation functions and renormalized fermionic vertices for varying loop orders. These benchmark calculations further substantiate the usage of multiloop pffRG solvers to tackle unconventional forms of quantum magnetism.

preprint2022arXiv

Fulfillment of sum rules and Ward identities in the multiloop functional renormalization group solution of the Anderson impurity model

We investigate several fundamental characteristics of the multiloop functional renormalization group (mfRG) flow by hands of its application to a prototypical many-electron system: the Anderson impurity model (AIM). We first analyze the convergence of the algorithm in the different parameter regions of the AIM. As no additional approximation is made, the multiloop series for the local self-energy and response functions converge perfectly to the corresponding results of the parquet approximation (PA) in the weak- to intermediate-coupling regime. Small oscillations of the mfRG solution as a function of the loop order gradually increase with the interaction, hindering a full convergence to the PA in the strong-coupling regime, where perturbative resummation schemes are no longer reliable. By exploiting the converged results, we inspect the fulfillment of (i) sum rules associated to the Pauli principle and (ii) Ward identities related to conservation laws. For the Pauli principle, we observe a systematic improvement by increasing the loop order and including the multiloop corrections to the self-energy. This is consistent with the preservation of crossing symmetries and two-particle self-consistency in the PA. For the Ward identities, we numerically confirm a visible improvement by means of the Katanin substitution. At weak coupling, violations of the Ward identity are further reduced by increasing the loop order in mfRG. In this regime, we also determine the precise scaling of the deviations of the Ward identity as a function of the electronic interaction. For larger interaction values, the overall behavior becomes more complex, and the benefits of the higher-loop terms are mostly present in the contributions at large frequencies.

preprint2022arXiv

Multiloop flow equations for single-boson exchange fRG

The recently introduced single-boson exchange (SBE) decomposition of the four-point vertex of interacting fermionic many-body systems is a conceptually and computationally appealing parametrization of the vertex. It relies on the notion of reducibility of vertex diagrams with respect to the bare interaction $U$, instead of a classification based on two-particle reducibility within the widely-used parquet decomposition. Here, we re-derive the SBE decomposition in a generalized framework (suitable for extensions to, e.g., inhomogeneous systems or real-frequency treatments) following from the parquet equations. We then derive multiloop functional renormalization group (mfRG) flow equations for the ingredients of this SBE decomposition, both in the parquet approximation, where the fully two-particle irreducible vertex is treated as an input, and in the more restrictive SBE approximation, where this role is taken by the fully $U$-irreducible vertex. Moreover, we give mfRG flow equations for the popular parametrization of the vertex in terms of asymptotic classes of the two-particle reducible vertices. Since the parquet and SBE decompositions are closely related, their mfRG flow equations are very similar in structure.

preprint2021arXiv

A beginner's guide to non-abelian iPEPS for correlated fermions

Infinite projected entangled pair states (iPEPS) have emerged as a powerful tool for studying interacting two-dimensional fermionic systems. In this review, we discuss the iPEPS construction and some basic properties of this tensor network (TN) ansatz. Special focus is put on (i) a gentle introduction of the diagrammatic TN representations forming the basis for deriving the complex numerical algorithm, and (ii) the technical advance of fully exploiting non-abelian symmetries for fermionic iPEPS treatments of multi-band lattice models. The exploitation of non-abelian symmetries substantially increases the performance of the algorithm, enabling the treatment of fermionic systems up to a bond dimension $D=24$ on a square lattice. A variety of complex two-dimensional (2D) models thus become numerically accessible. Here, we present first promising results for two types of multi-band Hubbard models, one with $2$ bands of spinful fermions of $\mathrm{SU}(2)_\mathrm{spin} \otimes \mathrm{SU}(2)_\mathrm{orb}$ symmetry, the other with $3$ flavors of spinless fermions of $\mathrm{SU}(3)_\mathrm{flavor}$ symmetry.

preprint2021arXiv

Keldysh Functional Renormalization Group Treatment of Finite-Ranged Interactions in Quantum Point Contacts

We combine two recently established methods, the extended Coupled-Ladder Approximation (eCLA) [Phys. Rev. B 95, 035122 (2017)] and a dynamic Keldysh functional Renormalization Group (fRG) approach for inhomogeneous systems [Phys. Rev. Lett. 119, 196401 (2017)] to tackle the problem of finite-ranged interactions in quantum point contacts (QPCs) at finite temperature. Working in the Keldysh formalism, we develop an eCLA framework, proceeding from a static to a fully dynamic description. Finally, we apply our new Keldysh eCLA method to a QPC model with finite-ranged interactions and show evidence that an interaction range comparable to the length of the QPC might be an essential ingredient for the development of a pronounced 0.7-shoulder in the linear conductance. We also discuss problems arising from a violation of a Ward identity in second-order fRG.

preprint2020arXiv

iPEPS study of spin symmetry in the doped $t$-$J$ model

We study the two-dimensional $t$-$J$ model on a square lattice using infinite projected entangled pair states (iPEPS). At small doping, multiple orders, such as antiferromagnetic order, stripe order and superconducting order, are intertwined or compete with each other. We demonstrate the role of spin symmetry at small doping by either imposing SU(2) spin symmetry or its U(1) subgroup in the iPEPS ansatz, thereby excluding or allowing spontaneous spin-symmetry breaking respectively in the thermodynamic limit. From a detailed comparison of our simulations, we provide evidence that stripe order is pinned by long-range antiferromagnetic order. We also find that SU(2) iPEPS, which enforces a spin-singlet state, yields a uniform charge distribution and favors d-wave singlet pairing.

preprint2020arXiv

Non-Fermi-liquid Kondo screening under Rabi driving

We investigate a Rabi-Kondo model describing an optically driven two-channel quantum dot device featuring a non-Fermi-liquid Kondo effect. Optically induced Rabi oscillation between the valence and conduction levels of the dot gives rise to a two-stage Kondo effect: Primary screening of the local spin is followed by secondary nonequilibrium screening of the local orbital degree of freedom. Using bosonization arguments and the numerical renormalization group, we compute the dot emission spectrum and residual entropy. Remarkably, both exhibit two-stage Kondo screening with non-Fermi-liquid properties at both stages.

preprint2020arXiv

Orbital differentiation in Hund metals

Orbital differentiation is a common theme in multiorbital systems, yet a complete understanding of it is still missing. Here, we consider a minimal model for orbital differentiation in Hund metals with a highly accurate method: We use the numerical renormalization group as a real-frequency impurity solver for a dynamical mean-field study of three-orbital Hubbard models, where a crystal field shifts one orbital in energy. The individual phases are characterized with dynamic correlation functions and their relation to diverse Kondo temperatures. Upon approaching the orbital-selective Mott transition, we find a strongly suppressed spin coherence scale and uncover the emergence of a singular Fermi liquid and interband doublon-holon excitations. Our theory describes the diverse polarization-driven phenomena in the $t_{2g}$ bands of materials such as ruthenates and iron-based superconductors, and our methodological advances pave the way towards real-frequency analyses of strongly correlated materials.

preprint2020arXiv

Quantum Many-Body Simulations of the 2D Fermi-Hubbard Model in Ultracold Optical Lattices

Understanding quantum many-body states of correlated electrons is one main theme in modern condensed matter physics. Given that the Fermi-Hubbard model, the prototype of correlated electrons, has been recently realized in ultracold optical lattices, it is highly desirable to have controlled numerical methodology to provide precise finite-temperature results upon doping, to directly compare with experiments. Here, we demonstrate the exponential tensor renormalization group (XTRG) algorithm [Phys. Rev. X 8, 031082 (2018)], complemented with independent determinant quantum Monte Carlo (DQMC) offer a powerful combination of tools for this purpose. XTRG provides full and accurate access to the density matrix and thus various spin and charge correlations, down to unprecedented low temperature of few percents of the fermion tunneling energy scale. We observe excellent agreement with ultracold fermion measurements at both half-filling and finite-doping, including the sign-reversal behavior in spin correlations due to formation of magnetic polarons, and the attractive hole-doublon and repulsive hole-hole pairs that are responsible for the peculiar bunching and antibunching behavior of the antimoments.

preprint2020arXiv

Strongly correlated materials from a numerical renormalization group perspective: How the Fermi-liquid state of Sr$_2$RuO$_4$ emerges

The crossover from fluctuating atomic constituents to a collective state as one lowers temperature or energy is at the heart of the dynamical mean-field theory description of the solid state. We demonstrate that the numerical renormalization group is a viable tool to monitor this crossover in a real-materials setting. The renormalization group flow from high to arbitrarily small energy scales clearly reveals the emergence of the Fermi-liquid state of Sr$_2$RuO$_4$. We find a two-stage screening process, where orbital fluctuations are screened at much higher energies than spin fluctuations, and Fermi-liquid behavior, concomitant with spin coherence, below a temperature of 25 K. By computing real-frequency correlation functions, we directly observe this spin--orbital scale separation and show that the van Hove singularity drives strong orbital differentiation. We extract quasiparticle interaction parameters from the low-energy spectrum and find an effective attraction in the spin-triplet sector.

preprint2019arXiv

Signatures of Mottness and Hundness in archetypal correlated metals

Physical properties of multi-orbital materials depend not only on the strength of the effective interactions among the valence electrons but also on their type. Strong correlations are caused by either Mott physics that captures the Coulomb repulsion among charges, or Hund physics that aligns the spins in different orbitals. We identify four energy scales marking the onset and the completion of screening in orbital and spin channels. The differences in these scales, which are manifest in the temperature dependence of the local spectrum and of the charge, spin and orbital susceptibilities, provide clear signatures distinguishing Mott and Hund physics. We illustrate these concepts with realistic studies of two archetypal strongly correlated materials, and corroborate the generality of our conclusions with a model Hamiltonian study.

preprint2016arXiv

Dynamic structure factor of the spin-1/2 XXZ chain in a transverse field

The spin-$\frac{1}{2}$ XXZ chain with easy-plane anisotropy in a transverse field describes well the thermodynamic properties of the material ${\rm Cs_2CoCl_4}$ in a wide range of temperatures and fields including the region close to the spin-flop Ising quantum phase transition. For a comparison with prospective inelastic neutron scattering experiments on this compound, we present results of an extensive numerical study of its dynamic structure factor $\mathcal{S}^{αβ}(k,ω)$ using matrix-product-state (MPS) techniques. Close to criticality, the dynamic part of the correlator $\mathcal{S}^{xx}$ longitudinal to the applied field is incoherent and possesses a small total weight as the ground state is already close to saturation. The transverse correlator $\mathcal{S}^{zz}$, on the other hand, is dominated by a coherent single-particle excitation with additional spectral weight at higher energies that we tentatively attribute to a repulsively bound pair of particles. With increasing temperature, the latter quickly fades and spectral weight instead accumulates close to zero wavevector just above the single-particle energy. On a technical level, we compare the numerical efficiency of real-time evolution to an MPS-based Chebyshev expansion in the present context, finding that both methods yield results of similar quality at comparable numerical costs.

preprint2016arXiv

Lindblad-Driven Discretized Leads for Non-Equilibrium Steady-State Transport in Quantum Impurity Models: Recovering the Continuum Limit

The description of interacting quantum impurity models in steady-state nonequilibrium is an open challenge for computational many-particle methods: the numerical requirement of using a finite number of lead levels and the physical requirement of describing a truly open quantum system are seemingly incompatible. One possibility to bridge this gap is the use of Lindblad-driven discretized leads (LDDL): one couples auxiliary continuous reservoirs to the discretized lead levels and represents these additional reservoirs by Lindblad terms in the Liouville equation. For quadratic models governed by Lindbladian dynamics, we present an elementary approach for obtaining correlation functions analytically. In a second part, we use this approach to explicitly discuss the conditions under which the continuum limit of the LDDL approach recovers the correct representation of thermal reservoirs. As an analytically solvable example, the nonequilibrium resonant level model is studied in greater detail. Lastly, we present ideas towards a numerical evaluation of the suggested Lindblad equation for interacting impurities based on matrix product states. In particular, we present a reformulation of the Lindblad equation, which has the useful property that the leads can be mapped onto a chain where both the Hamiltonian dynamics and the Lindblad driving are local at the same time. Moreover, we discuss the possibility to combine the Lindblad approach with a logarithmic discretization needed for the exploration of exponentially small energy scales.

preprint2015arXiv

Effect of spin-orbit interactions on the 0.7 anomaly in quantum point contacts

We study how the conductance of a quantum point contact is affected by spin-orbit interactions, for systems at zero temperature both with and without electron-electron interactions. In the presence of spin-orbit coupling, tuning the strength and direction of an external magnetic field can change the dispersion relation and hence the local density of states in the point contact region. This modifies the effect of electron-electron interactions, implying striking changes in the shape of the 0.7-anomaly and introducing additional distinctive features in the first conductance step.

preprint2015arXiv

Fermi-liquid theory for the single-impurity Anderson model

We generalize Nozières' Fermi-liquid theory for the low-energy behavior of the Kondo model to that of the single-impurity Anderson model. In addition to the electrons' phase shift at the Fermi energy, the low-energy Fermi-liquid theory is characterized by four Fermi-liquid parameters: the two given by Nozières that enter to first order in the excitation energy, and two additional ones that enter to second order and are needed away from particle-hole symmetry. We express all four parameters in terms of zero-temperature physical observables, namely the local charge and spin susceptibilities and their derivatives with respect to the local level position. We determine these in terms of the bare parameters of the Anderson model using Bethe Ansatz and Numerical Renormalization Group (NRG) calculations. Our low-energy Fermi-liquid theory applies throughout the crossover from the strong-coupling Kondo regime via the mixed-valence regime to the empty-orbital regime. From the Fermi-liquid theory, we determine the conductance through a quantum dot symmetrically coupled to two leads in the regime of small magnetic field, low temperature and small bias voltage, and compute the coefficients of the $\sim B^2$, $\sim T^2$, and $\sim V^2$ terms \textit{exactly} in terms of the Fermi-liquid parameters. The coefficients of $T^2$, $V^2$ and $B^2$ are found to change sign during the Kondo to empty-orbital crossover. The crossover becomes universal in the limit that the local interaction is much larger than the level width. For completeness, we also compute the shot noise and discuss the resulting Fano factor.

preprint2015arXiv

Hexagon-singlet solid ansatz for the spin-1 kagome antiferromagnet

We perform a systematic investigation on the hexagon-singlet solid (HSS) states, which are a class of spin liquid candidates for the spin-1 kagome antiferromagnet. With the Schwinger boson representation, we show that all HSS states have exponentially decaying correlations and can be interpreted as a (special) subset of the resonating Affleck-Kennedy-Lieb-Tasaki (AKLT) loop states. We provide a compact tensor network representation of the HSS states, with which we are able to calculate physical observables efficiently. We find that the HSS states have vanishing topological entanglement entropy, suggesting the absence of intrinsic topological order. We also employ the HSS states to perform a variational study of the spin-1 kagome Heisenberg antiferromagnetic model. When we use a restricted HSS ansatz preserving lattice symmetry, the best variational energy per site is found to be $e_0 = -1.3600$. In contrast, when allowing lattice symmetry breaking, we find a trimerized simplex valence bond crystal with a lower energy, $e_0=-1.3871$.

preprint2015arXiv

Simplex valence-bond crystal in the spin-1 kagome Heisenberg antiferromagnet

We investigate the ground state properties of a spin-1 kagome antiferromagnetic Heisenberg model using tensor-network (TN) methods. We obtain the energy per site {$e_0=-1.41090(2)$ with $D^*=8$ multiplets retained (i.e., a bond dimension of $D=24$), and $e_0=-1.4116(4)$ from large-$D$ extrapolation,} by accurate TN calculations directly in the thermodynamic limit. The symmetry between the two kinds of triangles is spontaneously broken, with a relative energy difference of $δ\approx$ 19\%, i.e, there is a trimerization (simplex) valence-bond order in the ground state. The spin-spin, dimer-dimer, and chirality-chirality correlation functions are found to decay exponentially with a rather short correlation length, showing that the ground state is gapped. We thus identify the ground state be a simplex valence-bond crystal (SVBC). We also discuss the spin-1 bilinear-biquadratic Heisenberg model on a kagome lattice, and determine its ground state phase diagram. Moreover, we implement non-abelian symmetries, here spin SU(2), in the TN algorithm, which improves the efficiency greatly and provides insight into the tensor structures.

preprint2015arXiv

Symmetric Minimally Entangled Typical Thermal States

We extend White's minimally entangled typically thermal states approach (METTS) to allow Abelian and non-Ablian symmetries to be exploited when computing finite-temperature response functions in one-dimensional (1D) quantum systems. Our approach, called SYMETTS, starts from a METTS sample of states that are not symmetry eigenstates, and generates from each a symmetry eigenstate. These symmetry states are then used to calculate dynamic response functions. SYMETTS is ideally suited to determine the low-temperature spectra of 1D quantum systems with high resolution. We employ this method to study a generalized diamond chain model for the natural mineral azurite Cu${}_3$(CO${}_3)_2$(OH$)_2$, which features a plateau at $\frac{1}{3}$ in the magnetization curve at low temperatures. Our calculations provide new insight into the effects of temperature on magnetization and excitation spectra in the plateau phase, which can be fully understood in terms of the microscopic model.

preprint2014arXiv

Equilibrium Fermi-liquid coefficients for the fully screened N-channel Kondo model

We analytically and numerically compute three equilibrium Fermi-liquid coefficients of the fully screened $N$-channel Kondo model, namely $c_B$, $c_T$ and $c_\varepsilon$, characterizing the magnetic field and temperature dependence of the resisitivity, and the curvature of the equilibrium Kondo resonance, respectively. We present a compact, unified derivation of the $N$-dependence of these coefficients, combining elements from various previous treatments of this model. We numerically compute these coefficients using the numerical renormalization group, with non-Abelian symmetries implemented explicitly, finding agreement with Fermi-liquid predictions on the order of 5% or better.

preprint2014arXiv

Functional Renormalization Group Approach for Inhomogeneous Interacting Fermi-Systems

The functional renormalization group (fRG) approach has the property that, in general, the flow equation for the two-particle vertex generates $\mathcal{O}(N^4)$ independent variables, where $N$ is the number of interacting states (e.g. sites of a real-space discretization). In order to include the flow equation for the two-particle vertex one needs to make further approximations if $N$ becomes too large. We present such an approximation scheme, called the coupled-ladder approximation, for the special case of onsite interaction. Like the generic third-order-truncated fRG, the coupled-ladder approximation is exact to second order and is closely related to a simultaneous treatment of the random phase approximation in all channels, i.e. summing up parquet-type diagrams. The scheme is applied to a one-dimensional model describing a quantum point contact.

preprint2014arXiv

Identifying Symmetry-Protected Topological Order by Entanglement Entropy

According to the classification using projective representations of the SO(3) group, there exist two topologically distinct gapped phases in spin-1 chains. The symmetry-protected topological (SPT) phase possesses half-integer projective representations of the SO(3) group, while the trivial phase possesses integer linear representations. In the present work, we implement non-Abelian symmetries in the density matrix renormalization group (DMRG) method, allowing us to keep track of (and also control) the virtual bond representations, and to readily distinguish the SPT phase from the trivial one by evaluating the multiplet entanglement spectrum. In particular, using the entropies $S^I$ ($S^H$) of integer (half-integer) representations, we can define an entanglement gap $G = S^I - S^H$, which equals 1 in the SPT phase, and $-1$ in the trivial phase. As application of our proposal, we study the spin-1 models on various 1D and quasi-1D lattices, including the bilinear-biquadratic model on the single chain, and the Heisenberg model on a two-leg ladder and a three-leg tube. Among others, we confirm the existence of an SPT phase in the spin-1 tube model, and reveal that the phase transition between the SPT and the trivial phase is a continuous one. The transition point is found to be critical, with conformal central charge $c=3$ determined by fits to the block entanglement entropy.

preprint2014arXiv

On the relation between the 0.7-anomaly and the Kondo effect: Geometric Crossover between a Quantum Point Contact and a Kondo Quantum Dot

Quantum point contacts (QPCs) and quantum dots (QDs), two elementary building blocks of semiconducting nanodevices, both exhibit famously anomalous conductance features: the 0.7-anomaly in the former case, the Kondo effect in the latter. For both the 0.7-anomaly and the Kondo effect, the conductance shows a remarkably similar low-energy dependence on temperature $T$, source-drain voltage $V_{\rm sd}$ and magnetic field $B$. In a recent publication [F. Bauer et al., Nature, 501, 73 (2013)], we argued that the reason for these similarities is that both a QPC and a KQD feature spin fluctuations that are induced by the sample geometry, confined in a small spatial regime, and enhanced by interactions. Here we further explore this notion experimentally and theoretically by studying the geometric crossover between a QD and a QPC, focussing on the $B$-field dependence of the conductance. We introduce a one-dimensional model that reproduces the essential features of the experiments, including a smooth transition between a Kondo QD and a QPC with 0.7-anomaly. We find that in both cases the anomalously strong negative magnetoconductance goes hand in hand with strongly enhanced local spin fluctuations. Our experimental observations include, in addition to the Kondo effect in a QD and the 0.7-anomaly in a QPC, Fano interference effects in a regime of coexistence between QD and QPC physics, and Fabry-Perot-type resonances on the conductance plateaus of a clean QPC. We argue that Fabry-Perot-type resonances occur generically if the electrostatic potential of the QPC generates a flatter-than-parabolic barrier top.

preprint2014arXiv

Two-bath spin-boson model: Phase diagram and critical properties

The spin-boson model, describing a two-level system coupled to a bath of harmonic oscillators, is a generic model for quantum dissipation, with manifold applications. It has also been studied as a simple example for an impurity quantum phase transition. Here we present a detailed study of a U(1)-symmetric two-bath spin-boson model, where two different components of an SU(2) spin 1/2 are coupled to separate dissipative baths. Non-trivial physics arises from the competition of the two dissipation channels, resulting in a variety of phases and quantum phase transitions. We employ a combination of analytical and numerical techniques to determine the properties of both the stable phases and the quantum critical points. In particular, we find a critical intermediate-coupling phase which is bounded by a continuous quantum phase transition which violates the quantum-to-classical correspondence.

preprint2013arXiv

Observing the Nonequilibrium Dynamics of the Quantum Transverse-Field Ising Chain in Circuit QED

We show how a quantum Ising spin chain in a time-dependent transverse magnetic field can be simulated and experimentally probed in the framework of circuit QED with current technology. The proposed setup provides a new platform for observing the nonequilibrium dynamics of interacting many-body systems. We calculate its spectrum to offer a guideline for its initial experimental characterization. We demonstrate that quench dynamics and the propagation of localized excitations can be observed with the proposed setup and discuss further possible applications and modifications of this circuit QED quantum simulator.

preprint2013arXiv

The quantum transverse-field Ising chain in circuit QED: effects of disorder on the nonequilibrium dynamics

We study several dynamical properties of a recently proposed implementation of the quantum transverse-field Ising chain in the framework of circuit QED. Particular emphasis is placed on the effects of disorder on the nonequilibrium behavior of the system. We show that small amounts of fabrication-induced disorder in the system parameters do not jeopardize the observation of previously-predicted phenomena. Based on a numerical extraction of the mean free path of the system, we also provide a simple quantitative estimate for certain disorder effects on the nonequilibrium dynamics of the circuit QED quantum simulator. We discuss the transition from weak to strong disorder, characterized by the onset of Anderson localization of the system's wave functions, and the qualitatively different dynamics it leads to.

preprint2012arXiv

Critical and strong-coupling phases in one- and two-bath spin-boson models

For phase transitions in dissipative quantum impurity models, the existence of a quantum-to-classical correspondence has been discussed extensively. We introduce a variational matrix product state approach involving an optimized boson basis, rendering possible high-accuracy numerical studies across the entire phase diagram. For the sub-ohmic spin-boson model with a power-law bath spectrum $\propto \w^s$, we confirm classical mean-field behavior for $s<1/2$, correcting earlier numerical renormalization-group results. We also provide the first results for an XY-symmetric model of a spin coupled to two competing bosonic baths, where we find a rich phase diagram, including both critical and strong-coupling phases for $s<1$, different from that of classical spin chains. This illustrates that symmetries are decisive for whether or not a quantum-to-classical correspondence exists.

preprint2012arXiv

Efficient simulation of infinite tree tensor network states on the Bethe lattice

We show that the simple update approach proposed by Jiang et. al. [H.C. Jiang, Z.Y. Weng, and T. Xiang, Phys. Rev. Lett. 101, 090603 (2008)] is an efficient and accurate method for determining the infinite tree tensor network states on the Bethe lattice. Ground state properties of the quantum transverse Ising model and the Heisenberg XXZ model on the Bethe lattice are studied. The transverse Ising model is found to undergo a second-order quantum phase transition with a diverging magnetic susceptibility but a finite correlation length which is upper-bounded by 1/ln(q-1) even at the transition point (q is the coordinate number of the Bethe lattice). An intuitive explanation on this peculiar "critical" phenomenon is given. The XXZ model on the Bethe lattice undergoes a first-order quantum phase transition at the isotropic point. Furthermore, the simple update scheme is found to be related with the Bethe approximation. Finally, by applying the simple update to various tree tensor clusters, we can obtain rather nice and scalable approximations for two-dimensional lattices.

preprint2012arXiv

Spin-1/2 Kondo effect in an InAs nanowire quantum dot: the Unitary limit, conductance scaling and Zeeman splitting

We report on a comprehensive study of spin-1/2 Kondo effect in a strongly-coupled quantum dot realized in a high-quality InAs nanowire. The nanowire quantum dot is relatively symmetrically coupled to its two leads, so the Kondo effect reaches the Unitary limit. The measured Kondo conductance demonstrates scaling with temperature, Zeeman magnetic field, and out-of-equilibrium bias. The suppression of the Kondo conductance with magnetic field is much stronger than would be expected based on a g-factor extracted from Zeeman splitting of the Kondo peak. This may be related to strong spin-orbit coupling in InAs.

preprint2012arXiv

Stroboscopic observation of quantum many-body dynamics

Recent experiments have demonstrated single-site resolved observation of cold atoms in optical lattices. Thus, in the future it may be possible to take repeated snapshots of an interacting quantum many-body system during the course of its evolution. Here we address the impact of the resulting Quantum (anti-)Zeno physics on the many-body dynamics. We use the time-dependent density-matrix renormalization group to obtain the time evolution of the full many-body wave function, which is then periodically projected in order to simulate realizations of stroboscopic measurements. For the example of a one-dimensional lattice of spin-polarized fermions with nearest-neighbor interactions, we find regimes for which many-particle configurations are stabilized and destabilized depending on the interaction strength and the time between observations.

preprint2012arXiv

SU(3) Anderson impurity model: A numerical renormalization group approach exploiting non-Abelian symmetries

We show how the density-matrix numerical renormalization group (DM-NRG) method can be used in combination with non-Abelian symmetries such as SU(N), where the decomposition of the direct product of two irreducible representations requires the use of a so-called outer multiplicity label. We apply this scheme to the SU(3) symmetrical Anderson model, for which we analyze the finite size spectrum, determine local fermionic, spin, superconducting, and trion spectral functions, and also compute the temperature dependence of the conductance. Our calculations reveal a rich Fermi liquid structure.

preprint2011arXiv

A numerical algorithm for the explicit calculation of SU(N) and SL(N,C) Clebsch-Gordan coefficients

We present an algorithm for the explicit numerical calculation of SU(N) and SL(N,C) Clebsch-Gordan coefficients, based on the Gelfand-Tsetlin pattern calculus. Our algorithm is well-suited for numerical implementation; we include a computer code in an appendix. Our exposition presumes only familiarity with the representation theory of SU(2).

preprint2011arXiv

Anderson Orthogonality and the Numerical Renormalization Group

Anderson Orthogonality (AO) refers to the fact that the ground states of two Fermi seas that experience different local scattering potentials, say |G_I> and |G_F>, become orthogonal in the thermodynamic limit of large particle number N, in that |<G_I|G_F>| ~ N^(- Delta_AO^2 /2) for N->infinity. We show that the numerical renormalization group (NRG) offers a simple and precise way to calculate the exponent Delta_AO: the overlap, calculated as function of Wilson chain length k, decays exponentially, ~ exp(-k alpha), and Delta_AO can be extracted directly from the exponent alpha. The results for Delta_AO so obtained are consistent (with relative errors typically smaller than 1%) with two other related quantities that compare how ground state properties change upon switching from |G_I> to |G_F>: the difference in scattering phase shifts at the Fermi energy, and the displaced charge flowing in from infinity. We illustrate this for several nontrivial interacting models, including systems that exhibit population switching.

preprint2011arXiv

Anderson Orthogonality in the Dynamics After a Local Quantum Quench

We present a systematic study of the role of Anderson orthogonality for the dynamics after a quantum quench in quantum impurity models, using the numerical renormalization group. As shown by Anderson in 1967, the scattering phase shifts of the single-particle wave functions constituting the Fermi sea have to adjust in response to the sudden change in the local parameters of the Hamiltonian, causing the initial and final ground states to be orthogonal. This so-called Anderson orthogonality catastrophe also influences dynamical properties, such as spectral functions. Their low-frequency behaviour shows nontrivial power laws, with exponents that can be understood using a generalization of simple arguments introduced by Hopfield and others for the X-ray edge singularity problem. The goal of this work is to formulate these generalized rules, as well as to numerically illustrate them for quantum quenches in impurity models involving local interactions. As a simple yet instructive example, we use the interacting resonant level model as testing ground for our generalized Hopfield rule. We then analyse a model exhibiting population switching between two dot levels as a function of gate voltage, probed by a local Coulomb interaction with an additional lead serving as charge sensor. We confirm a recent prediction that charge sensing can induce a quantum phase transition for this system, causing the population switch to become abrupt. We elucidate the role of Anderson orthogonality for this effect by explicitly calculating the relevant orthogonality exponents.

preprint2011arXiv

Chebyshev matrix product state approach for spectral functions

We show that recursively generated Chebyshev expansions offer numerically efficient representations for calculating zero-temperature spectral functions of one-dimensional lattice models using matrix product state (MPS) methods. The main features of this Chebychev matrix product state (CheMPS) approach are: (i) it achieves uniform resolution over the spectral function's entire spectral width; (ii) it can exploit the fact that the latter can be much smaller than the model's many-body bandwidth; (iii) it offers a well-controlled broadening scheme; (iv) it is based on a succession of Chebychev vectors |t_n>, (v) whose entanglement entropies were found to remain bounded with increasing recursion order n for all cases analyzed here; (vi) it distributes the total entanglement entropy that accumulates with increasing n over the set of Chebyshev vectors |t_n>. We present zero-temperature CheMPS results for the structure factor of spin-1/2 antiferromagnetic Heisenberg chains and perform a detailed finite-size analysis. Making comparisons to three benchmark methods, we find that CheMPS (1) yields results comparable in quality to those of correction vector DMRG, at dramatically reduced numerical cost; (2) agrees well with Bethe Ansatz results for an infinite system, within the limitations expected for numerics on finite systems; (3) can also be applied in the time domain, where it has potential to serve as a viable alternative to time-dependent DMRG (in particular at finite temperatures). Finally, we present a detailed error analysis of CheMPS for the case of the noninteracting resonant level model.

preprint2011arXiv

Quantum quench of Kondo correlations in optical absorption

The interaction between a single confined spin and the spins of a Fermionic reservoir leads to one of the most spectacular phenomena of many body physics -- the Kondo effect. Here we report the observation of Kondo correlations in optical absorption measurements on a single semiconductor quantum dot tunnel-coupled to a degenerate electron gas. In stark contrast to transport experiments, absorption of a single photon leads to an abrupt change in the system Hamiltonian and a quantum quench of Kondo correlations. By inferring the characteristic power law exponents from the experimental absorption line-shapes, we find a unique signature of the quench in the form of an Anderson orthogonality catastrophe, originating from a vanishing overlap between the initial and final many-body wave-functions. We also show that the power-law exponents that determine the degree of orthogonality can be tuned by applying an external magnetic field which gradually turns the Kondo correlations off. Our experiments demonstrate that optical measurements on single artificial atoms offer new perspectives on many-body phenomena previously studied exclusively using transport spectroscopy. Moreover, they initiate a new paradigm for quantum optics where many-body physics influences electric field and intensity correlations.

preprint2011arXiv

Superradiant Phase Transitions and the Standard Description of Circuit QED

We investigate the equilibrium behaviour of a superconducting circuit QED system containing a large number of artificial atoms. It is shown that the currently accepted standard description of circuit QED via an effective model fails in an important aspect: it predicts the possibility of a superradiant quantum phase transition, even though a full microscopic treatment reveals that a no-go theorem for such phase transitions known from cavity QED applies to circuit QED systems as well. We generalize the no-go theorem to the case of (artificial) atoms with many energy levels and thus make it more applicable for realistic cavity or circuit QED systems.

preprint2011arXiv

Transport and Dephasing in a Quantum Dot: Multiply Connected Graph Model

Using the theory of diffusion in graphs, we propose a model to study mesoscopic transport through a diffusive quantum dot. The graph consists of three quasi-1D regions: a central region describing the dot, and two identical left- and right- wires connected to leads, which mimic contacts of a real system. We find the exact solution of the diffusion equation for this graph and evaluate the conductance including quantum corrections. Our model is complementary to the RMT-models describing quantum dots. Firstly, it reproduces the universal limit at zero temperature. But the main advantage compared to RMT-models is that it allows one to take into account interaction-induced dephasing at finite temperatures. Besides, the crossovers from open to almost closed quantum dots and between different regimes of dephasing can be described within a single framework. We present results for the temperature dependence of the weak localization correction to the conductance for the experimentally relevant parameter range and discuss the possibility to observe the elusive 0D-regime of dephasing in different mesoscopic systems.

preprint2010arXiv

Matrix product state approach for a two-lead, multi-level Anderson impurity model

We exploit the common mathematical structure of the numerical renormalization group and the density matrix renormalization group, namely, matrix product states, to implement an efficient numerical treatment of a two-lead, multi-level Anderson impurity model. By adopting a star-like geometry, where each species (spin and lead) of conduction electrons is described by its own Wilson chain, instead of using a single Wilson chain for all species together, we achieve a very significant reduction in the numerical resources required to obtain reliable results. We illustrate the power of this approach by calculating ground state properties of a four-level quantum dot coupled to two leads. The success of this proof-of-principle calculation suggests that the star geometry constitutes a promising strategy for future calculations the ground state properties of multi-band, multi-level quantum impurity models. Moreover, we show that it is possible to find an "optimal" chain basis, obtained via a unitary transformation (acting only on the index distinguishing different Wilson chains), in which degrees of freedom on different Wilson chains become effectively decoupled from each other further out on the Wilson chains. This basis turns out to also diagonalize the model's chain-to-chain scattering matrix. We demonstrate this for a spinless two-lead model, presenting DMRG-results for the mutual information between two sites located far apart on different Wilson chains, and NRG results with respect to the scattering matrix.

preprint2009arXiv

Constrained optimization of sequentially generated entangled multiqubit states

We demonstrate how the matrix-product state formalism provides a flexible structure to solve the constrained optimization problem associated with the sequential generation of entangled multiqubit states under experimental restrictions. We consider a realistic scenario in which an ancillary system with a limited number of levels performs restricted sequential interactions with qubits in a row. The proposed method relies on a suitable local optimization procedure, yielding an efficient recipe for the realistic and approximate sequential generation of any entangled multiqubit state. We give paradigmatic examples that may be of interest for theoretical and experimental developments.

preprint2009arXiv

The dephasing rate formula in the many body context

We suggest a straightforward approach to the calculation of the dephasing rate in a fermionic system, which correctly keeps track of the crucial physics of Pauli blocking. Starting from Fermi's golden rule, the dephasing rate can be written as an integral over the frequency transferred between system and environment, weighted by their respective spectral densities. We show that treating the full many-fermion system instead of a single particle automatically enforces the Pauli principle. Furthermore, we explain the relation to diagrammatics. Finally, we show how to treat the more involved strong-coupling case when interactions appreciably modify the spectra. This is relevant for the situation in disordered metals, where screening is important.

preprint2008arXiv

Influence Functional for Decoherence of Interacting Electrons in Disordered Conductors

We have rederived the controversial influence functional approach of Golubev and Zaikin (GZ) for interacting electrons in disordered metals in a way that allows us to show its equivalence, before disorder averaging, to diagrammatic Keldysh perturbation theory. By representing a certain Pauli factor (1-2 rho) occuring in GZ's effective action in the frequency domain (instead of the time domain, as GZ do), we also achieve a more accurate treatment of recoil effects. With this change, GZ's approach reproduces, in a remarkably simple way, the standard, generally accepted result for the decoherence rate. -- The main text and appendix A.1 to A.3 of the present paper have already been published previously; for convenience, they are included here again, together with five additional, lengthy appendices containing relevant technical details.

preprint2008arXiv

Interplay of mesoscopic and Kondo effects for transmission amplitude of few-level quantum dots

The magnitude and phase of the transmission amplitude of a multi-level quantum dot is calculated for the mesoscopic regime of level spacing large compared to level width. The interplay between Kondo correlations and the influence by neighboring levels is discussed. As in the single-level case, the Kondo plateaus of magnitude and phase disappear with increasing temperature. At certain gate voltages, ``stationary'' points are found at which the transmission phase is independent of temperature. Depending on the mesoscopic parameters of the adjacent levels (like relative sign and magnitude of tunneling matrix elements), the stationary points are shifted to or repelled by the neighboring level.

preprint2008arXiv

Numerical renormalization group calculation of near-gap peaks in spectral functions of the Anderson model with superconducting leads

We use the numerical renormalization group method (NRG) to investigate a single-impurity Anderson model with a coupling of the impurity to a superconducting host. Analysis of the energy flow shows, in contrast to previous belief, that NRG iterations can be performed up to a large number of sites, corresponding to energy differences far below the superconducting gap. This allows us to calculate the impurity spectral function very accurately for frequencies near the gap edge, and to resolve, in a certain parameter regime, sharp peaks in the spectral function close to the gap edge.

preprint2003arXiv

Decoherence without dissipation?

In a recent article, Ford, Lewis and O'Connell (PRA 64, 032101 (2001)) discuss a thought experiment in which a Brownian particle is subjected to a double-slit measurement. Analyzing the decay of the emerging interference pattern, they derive a decoherence rate that is much faster than previous results and even persists in the limit of vanishing dissipation. This result is based on the definition of a certain attenuation factor, which they analyze for short times. In this note, we point out that this attenuation factor captures the physics of decoherence only for times larger than a certain time t_mix, which is the time it takes until the two emerging wave packets begin to overlap. Therefore, the strategy of Ford et al of extracting the decoherence time from the regime t < t_mix is in our opinion not meaningful. If one analyzes the attenuation factor for t > t_mix, one recovers familiar behaviour for the decoherence time; in particular, no decoherence is seen in the absence of dissipation. The latter conclusion is confirmed with a simple calculation of the off-diagonal elements of the reduced density matrix.

preprint2003arXiv

Flow equation renormalization of a spin-boson model with a structured bath

We discuss the dynamics of a spin coupled to a damped harmonic oscillator. This system can be mapped to a spin-boson model with a structured bath, i.e. the spectral function of the bath has a resonance peak. We diagonalize the model by means of infinitesimal unitary transformations (flow equations), thereby decoupling the small quantum system from its environment, and calculate spin-spin correlation functions.

preprint2002arXiv

Decoherence of interacting electrons in disordered conductors: on the relation between influence functional and diagrammatic approaches

We establish a connection between the influence functional approach of Golubev and Zaikin (GZ) and Keldysh diagrammatic perturbation theory for calculating the decoherence time of interacting electrons in disordered metals; we show how the standard diagrams for the Cooperon self energy can be recovered from GZ's influence functional $e^{- (i S_R + S_I)}$. This allows us to shed light on GZ's claim that $S_R$ is irrelevant for decoherence: $S_R$ generates as many important self energy diagrams as $S_I$; GZ's neglect of $S_R$ is permissible only at high temperatures ($T > \hbar / τ_{el}$).

preprint1998arXiv

Bosonization for Beginners --- Refermionization for Experts

This tutorial gives an elementary and self-contained review of abelian bosonization in 1 dimension in a system of finite size $L$, following and simplifying Haldane's constructive approach. As a non-trivial application, we rigorously resolve (following Furusaki) a recent controversy regarding the tunneling density of states, $ρ_{dos} (ω)$, at the site of an impurity in a Tomonaga-Luttinger liquid: we use finite-size refermionization to show exactly that for g=1/2 its asymptotic low-energy behavior is $ρ_{dos}(ω) \sim ω$. This agrees with the results of Fabrizio & Gogolin and of Furusaki, but not with those of Oreg and Finkel'stein (probably because we capture effects not included in their mean-field treatment of the Coulomb gas that they obtained by an exact mapping; their treatment of anti-commutation relations in this mapping is correct, however, contrary to recent suggestions in the literature). --- The tutorial is addressed to readers unfamiliar with bosonization, or for those interested in seeing ``all the details'' explicitly; it requires knowledge of second quantization only, not of field theory. At the same time, we hope that experts too might find useful our explicit treatment of certain subtleties -- these include the proper treatment of the so-called Klein factors that act as fermion-number ladder operators (and also ensure the anti-commutation of different species of fermion fields), the retention of terms of order 1/L, and a novel, rigorous formulation of finite-size refermionization of both $e^{-i Φ(x)}$ and the boson field $Φ(x)$ itself.

preprint1998arXiv

The 2-Channel Kondo Model I: Review of Experimental Evidence for its Realization in Metal Nanoconstrictions

Certain zero-bias anomalies (ZBAs) in the voltage, temperature and magnetic field dependence of the conductance $G(V,T,H)$ of quenched Cu point contacts have previously been interpreted to be due to non-magnetic 2-channel Kondo (2CK) scattering from near-degenerate atomic two-level tunneling systems (Ralph and Buhrman, 1992; Ralph et al. 1994), and hence to represent an experimental realization of the non-Fermi-liquid physics of the T=0 fixed point of the 2-channel Kondo model. In this, the first in a series of three papers (I,II,III) devoted to 2-channel Kondo physics, we present a comprehensive review of the quenched Cu ZBA experiments and their 2CK interpretation, including new results on ZBAs in constrictions made from Ti or from metallic glasses. We first review the evidence that the ZBAs are due to electron scattering from stuctural defects that are not static, but possess internal dynamics. In order to distinguish between several mechanisms proposed to explain the experiments, we then analyze the scaling properties of the conductance at low temperature and voltage and extract from the data a universal scaling function $Γ(v)$. The theoretical calculation of the corresponding scaling function within the 2CK model is the subject of papers II and III. The main conclusion of our work is that the properties of the ZBAs, and most notably their scaling behavior, are in good agreement with the 2CK model and clearly different from several other proposed mechanisms.