Source author record

Andreas Weichselbaum

Andreas Weichselbaum 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

30works
6topics
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

30 published item(s)

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.

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

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

X-Symbols for Non-Abelian Symmetries in Tensor Networks

The full exploitation of non-abelian symmetries in tensor network states (TNS) derived from a given lattice Hamiltonian is highly attractive in various aspects. From a theoretical perspective, it can offer deep insights into the entanglement structure and quantum information content of strongly correlated quantum many-body states. From a practical perspective, it allows one to push numerical efficiency by orders of magnitude. Physical expectation values based on TNS require the full contraction of a given tensor network, with the elementary ingredient being a pairwise contraction. While well-established for no or just abelian symmetries, this can become quickly extremely involved and cumbersome for general non-abelian symmetries. As shown in this work, however, the latter can be tackled in a transparent and efficient manner by introducing so-called X-symbols which deal with the underlying pairwise contraction of generalized Clebsch-Gordan tensors (CGTs). These X-symbols can be computed deterministically once and for all, and hence also be tabulated. Akin to 6j-symbols, X-symbols are generally much smaller than their constituting CGTs. In applications, they solely affect the tensors of reduced matrix elements, and therefore, once tabulated, allow one to completely sidestep the explicit usage of CGTs, and thus to greatly increase numerical efficiency.

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

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

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

Local susceptibility and Kondo scaling in the presence of finite bandwidth

The Kondo scale TK for impurity systems is expected to guarantee universal scaling of physical quantities. However, in practice, not every definition of TK necessarily supports this notion away from the strict scaling limit. Specifically, this paper addresses the role of finite bandwidth D in the strongly-correlated Kondo regime. For this, various theoretical definitions of TK are analyzed based on the inverse magnetic impurity susceptibility at zero temperature. While conventional definitions in that respect quickly fail to ensure universal Kondo scaling for all D, this paper proposes an altered definition of TK_{sc} that allows universal scaling of dynamical or thermal quantities for a given fixed Hamiltonian. If the scaling is performed with respect to an external parameter which directly enters the Hamiltonian, such as magnetic field, the corresponding TK_{sc;B} for universal scaling differs, yet becomes equivalent to TK_{sc} in the scaling limit. The only requirement for universal scaling in the full Kondo parameter regime with a residual error of less than 1% is a well-defined isolated Kondo feature with TK < 0.01D. By varying D over a wide range relative to the bare energies of the impurity, for example, this allows a smooth transition from the Anderson to the Kondo model.

preprint2014arXiv

Non-Fermi liquid behaviour in non-equilibrium transport through Co doped Au chains connected to four-fold symmetric leads

We calculate the differential conductance as a function of temperature and bias voltage, $G(T,V)$, through Au monoatomic chains with a substitutional Co atom as a magnetic impurity, connected to a four-fold symmetric lead. The system was recently proposed as a possible scenario for observation of the overscreened Kondo physics. Stretching the chain, the system could be tuned through a quantum critical point (QCP) with three different regimes, overscreened, underscreened and non Kondo phases. We present calculations of the impurity spectral function by using the numerical renormalization group (NRG) for the three different regimes characterizing the QCP. Non trivial behaviour of the spectral function is reported near the QCP. Comparison with results using the non crossing approximation (NCA), shows that the latter is reliable in the overscreened regime, when the anisotropy is larger than the Kondo temperature. For these parameters, which correspond to realistic previous estimates, $G(T,V)$ calculated within NCA exhibits clear signatures of the non-Fermi liquid behaviour within the overscreened regime.

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.

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

Non-abelian symmetries in tensor networks: a quantum symmetry space approach

A general framework for non-abelian symmetries is presented for matrix-product and tensor-network states in the presence of orthonormal local as well as effective basis sets. The two crucial ingredients, the Clebsch-Gordan algebra for multiplet spaces as well as the Wigner-Eckart theorem for operators, are accounted for in a natural, well-organized, and computationally straightforward way. The unifying tensor-representation for quantum symmetry spaces, dubbed QSpace, is particularly suitable to deal with standard renormalization group algorithms such as the numerical renormalization group (NRG), the density matrix renormalization group (DMRG), or also more general tensor networks such as the multi-scale entanglement renormalization ansatz (MERA). In this paper, the focus is on the application of the non-abelian framework within the NRG. A detailed analysis is given for a fully screened spin-3/2 three-channel Anderson impurity model in the presence of conservation of total spin, particle-hole symmetry, and SU(3) channel symmetry. The same system is analyzed using several alternative symmetry scenarios. This includes the more traditional symmetry setting SU(2)^4, the larger symmetry SU(2)*U(1)*SU(3), and their much larger enveloping symplectic symmetry SU(2)*Sp(6). These are compared in detail, including their respective dramatic gain in numerical efficiency. In the appendix, finally, an extensive introduction to non-abelian symmetries is given for practical applications, together with simple self-contained numerical procedures to obtain Clebsch-Gordan coefficients and irreducible operators sets. The resulting QSpace tensors can deal with any set of abelian symmetries together with arbitrary non-abelian symmetries with compact, i.e. finite-dimensional, semi-simple Lie algebras.

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

Tensor networks and the numerical renormalization group

The full-density-matrix numerical renormalization group (NRG) has evolved as a systematic and transparent setting for the cal- culation of thermodynamical quantities at arbitrary temperatures within the NRG framework. It directly evaluates the relevant Lehmann representations based on the complete basis sets intro- duced by Anders and Schiller (2005). In addition, specific attention is given to the possible feedback from low energy physics to high energies by the explicit and careful construction of the full thermal density matrix, naturally generated over a distribution of energy shells. Specific examples are given in terms of spectral functions (fdmNRG), time-dependent NRG (tdmNRG), Fermi-Golden rule calculations (fgrNRG), as well as the calculation of plain thermodynamic expectation values. Furthermore, based on the very fact that, by its iterative nature, the NRG eigenstates are naturally described in terms of matrix product states, the language of tensor networks has proven enormously convenient in the description of the underlying algorithmic procedures. This paper therefore also provides a detailed introduction and discussion of the prototypical NRG calculations in terms of their corresponding tensor networks.

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

Discarded weight and entanglement spectra in the Numerical Renormalization Group

A quantitative criterion to prove and analyze convergence within the numerical renormalization group (NRG) is introduced. By tracing out a few further NRG shells, the resulting reduced density matrices carry relevant information on numerical accuracy as well as entanglement. Their spectra can thus be analyzed twofold. The smallest eigenvalues provide a sensitive estimate of how much weight is discarded in the low energy description of later iterations. As such, the discarded weight is a clear indicator of the accuracy of a specific NRG calculation. In particular, it indicates in a site-specific manner whether sufficiently many states have been kept within a single NRG run. The largest eigenvalues of the reduced density matrices, on the other hand, lend themselves to a straightforward analysis in terms of entanglement spectra, which can be combined into entanglement flow diagrams. The latter show strong similarities with the well-known standard energy flow diagram of the NRG, supporting the prevalent usage of entanglement spectra to characterize different physical regimes.

preprint2011arXiv

Incommensurate correlations in the anisotropic triangular Heisenberg lattice

We study the anisotropic spin-half antiferromagnetic triangular Heisenberg lattice in two dimensions, seen as a set of chains with couplings J (J') along (in between) chains, respectively. Our focus is on the incommensurate correlation that emerges in this system in a wide parameter range due to the intrinsic frustration of the spins. We study this system with traditional DMRG using cylindrical boundary conditions to least constrain possible incommensurate order. Despite that the limit of essentially decoupled chains J'/J < 0.5 is not very accessible numerically, it appears that the spin-spin correlations remain incommensurate for any finite 0 < J' < Jc', where Jc'/J > 1. The incommensurate wave vector q_J, however, approaches the commensurate value corresponding to the antiferromagnetic correlation of a single chain very rapidly with decreasing J'/J, roughly as q_J ~ pi - c_1 (J'/J)^n exp(-c_2 J/J').

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.

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.

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.