Source author record

Ulrich Schollwöck

Ulrich Schollwöck 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

27works
9topics
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

27 published item(s)

preprint2022arXiv

Snapshot-based detection of $\frac{1}{2}$-Laughlin states: coupled chains and central charge

Experimental realizations of topologically ordered states of matter, such as fractional quantum Hall states, with cold atoms are now within reach. In particular, optical lattices provide a promising platform for the realization and characterization of such states, where novel detection schemes enable an unprecedented microscopic understanding. Here we show that the central charge can be directly measured in current cold atom experiments using the number entropy as a proxy for the entanglement entropy. We perform density-matrix renormalization-group simulations of Hubbard-interacting bosons on coupled chains subject to a magnetic field with $α=\frac{1}{4}$ flux quanta per plaquette. Tuning the inter-chain hopping, we find a transition from a trivial quasi-one dimensional phase to the topologically ordered Laughlin state at magnetic filling factor $ν=\frac{1}{2}$ for systems of three or more chains. We resolve the transition using the central charge, on-site correlations, momentum distributions and the many-body Chern number. Additionally, we propose a scheme to experimentally estimate the central charge from Fock basis snapshots. The model studied here is experimentally realizable with existing cold atom techniques and the proposed observables pave the way for the detection and classification of a larger class of interacting topological states of matter.

preprint2021arXiv

Concept of orbital entanglement and correlation in quantum chemistry

A recent development in quantum chemistry has established the quantum mutual information between orbitals as a major descriptor of electronic structure. This has already facilitated remarkable improvements of numerical methods and may lead to a more comprehensive foundation for chemical bonding theory. Building on this promising development, our work provides a refined discussion of quantum information theoretical concepts by introducing the physical correlation and its separation into classical and quantum parts as distinctive quantifiers of electronic structure. In particular, we succeed in quantifying the entanglement. Intriguingly, our results for different molecules reveal that the total correlation between orbitals is mainly classical, raising questions about the general significance of entanglement in chemical bonding. Our work also shows that implementing the fundamental particle number superselection rule, so far not accounted for in quantum chemistry, removes a major part of correlation and entanglement previously seen. In that respect, realizing quantum information processing tasks with molecular systems might be more challenging than anticipated.

preprint2021arXiv

Probing the Hall Voltage in Synthetic Quantum Systems

In the context of experimental advances in the realization of artificial magnetic fields in quantum gases, we discuss feasible schemes to extend measurements of the Hall polarization to a study of the Hall voltage, allowing for direct comparison with solid state systems. Specifically, for the paradigmatic example of interacting flux ladders, we report on characteristic zero crossings and a remarkable robustness of the Hall voltage with respect to interaction strengths, particle fillings, and ladder geometries, which is unobservable in the Hall polarization. Moreover, we investigate the site-resolved Hall response in spatially inhomogeneous quantum phases.

preprint2021arXiv

Snapshot based characterization of particle currents and the Hall response in synthetic flux lattices

Quantum simulators are attracting great interest because they promise insight into the behavior of quantum many-body systems that are prohibitive for classical simulations. The generic output of quantum simulators are snapshots, obtained by means of projective measurements. A central goal of theoretical efforts must be to predict the exact same quantities that can be measured in experiments. Here, we report on the snapshot based calculation of particle currents in quantum lattice models with a conserved number of particles. It is shown how the full probability distribution of locally resolved particle currents can be obtained from suitable snapshot data. Moreover, we investigate the Hall response of interacting bosonic flux ladders, exploiting snapshots drawn from matrix-product states. Flux ladders are minimal lattice models, which enable microscopic studies of the Hall response in correlated quantum phases and they are successfully realized in current quantum-gas experiments. Using a specific pattern of unitary two-site transformations, it is shown that the Hall polarization and the Hall voltage can be faithfully computed from the snapshots obtained in experimentally feasible quench and finite-bias simulations.

preprint2021arXiv

Stochastic Adaptive Single-Site Time-Dependent Variational Principle

In recent years, the time-dependent variational principle (TDVP) method based on the matrix product state (MPS) wave function formulation has shown its great power in performing large-scale quantum dynamics simulations for realistic chemical systems with strong electron-vibration interactions. In this work, we propose a new stochastic adaptive single-site TDVP (SA-1TDVP) scheme to evolve the bond-dimension adaptively, which can integrate the tra-ditional advantages of both the high efficiency of single-site TDVP (1TDVP) variant and the high accuracy of the two-site TDVP (2TDVP) variant. Based on the assumption that the level statistics of entanglement Hamiltonians, which originate from the reduced density matrices of the MPS method, follows a Poisson or Wigner distribution, as generically predicted by random matrix theory, addi-tional random singular values are generated to expand the bond-dimension automatically. Tests on simulating the vibrationally-resolved quantum dynamics and absorption spectra in the pyrazine molecule and perylene bisimide (PBI) J-aggregate trimer as well as a spin-1/2 Heisenberg chain show that it can be automatic and as accurate as 2TDVP but reduce the computational time remarkably.

preprint2020arXiv

Absence of superconductivity in the pure two-dimensional Hubbard model

We study the superconducting pairing correlations in the ground state of the doped Hubbard model -- in its original form without hopping beyond nearest neighbor or other perturbing parameters -- in two dimensions at intermediate to strong coupling and near optimal doping. The nature of such correlations has been a central question ever since the discovery of cuprate high-temperature superconductors. Despite unprecedented effort and tremendous progress in understanding the properties of this fundamental model, a definitive answer to whether the ground state is superconducting in the parameter regime most relevant to cuprates has proved exceedingly difficult to establish. In this work, we employ two complementary, state-of-the-art many-body computational methods, constrained path (CP) auxiliary-field quantum Monte Carlo (AFQMC) and density matrix renormalization group (DMRG) methods, deploying the most recent algorithmic advances in each. Systematic and detailed comparisons between the two methods are performed. The DMRG is extremely reliable on small width cylinders, where we use it to validate the AFQMC. The AFQMC is then used to study wide systems as well as fully periodic systems, to establish that we have reached the thermodynamic limit. The ground state is found to be non-superconducting in the moderate to strong coupling regime in the vicinity of optimal hole doping.

preprint2020arXiv

Imaginary-time matrix product state impurity solver in a real material calculation: Spin-orbit coupling in Sr$_2$RuO$_4$

Using an imaginary-time matrix-product state (MPS) based quantum impurity solver we perform a realistic dynamical mean-field theory (DMFT) calculation combined with density functional theory (DFT) for Sr$_2$RuO$_4$. We take the full Hubbard-Kanamori interactions and spin-orbit coupling (SOC) into account. The MPS impurity solver works at essentially zero temperature in the presence of SOC, a regime of parameters currently inaccessible to continuous-time quantum Monte Carlo (CTQMC) methods, due to a severe sign problem. We show that earlier results obtained at high temperature, namely that the diagonal self-energies are nearly unaffected by SOC and that interactions lead to an effective enhancement of the SOC, hold even at low temperature. We observe that realism makes the numerical solution of the impurity model with MPS much more demanding in comparison to earlier works on Bethe lattice models, requiring several algorithmic improvements.

preprint2020arXiv

Plaquette versus ordinary $d$-wave pairing in the $t'$-Hubbard model on a width 4 cylinder

The Hubbard model and its extensions are important microscopic models for understanding high- $T_c$ superconductivity in cuprates. In the model with next-nearest-neighbor hopping $t'$ (the $t'$- Hubbard model), pairing is strongly influenced by $t'$ . In particular, a recent study on a width-4 cylinder observed quasi-long-rage superconducting order, associated with a negative $t'$ , which was taken to imply superconductivity in the two-dimensional (2D) limit. In this work we study more carefully pairing in the width-4 $t'$-Hubbard model. We show that in this specific system, the pairing symmetry with $t'<0$ is not the ordinary $d$-wave one would expect in the 2D limit. Instead we observe a so-called plaquette d-wave pairing. The plaquette d-wave exists only on a width-4 cylinder, and so is not representative of the 2D limit. We find that a negative $t'$ suppresses the conventional d-wave, leading to plaquette pairing. In contrast, a different $t''$ coupling acting diagonally on the plaquettes suppresses plaquette pairing, leading to conventional $d$-wave pairing.

preprint2015arXiv

A Strictly Single-Site DMRG Algorithm with Subspace Expansion

We introduce a strictly single-site DMRG algorithm based on the subspace expansion of the Alternating Minimal Energy (AMEn) method. The proposed new MPS basis enrichment method is sufficient to avoid local minima during the optimisation, similarly to the density matrix perturbation method, but computationally cheaper. Each application of $\hat H$ to $|Ψ\rangle$ in the central eigensolver is reduced in cost for a speed-up of $\approx (d + 1)/2$, with $d$ the physical site dimension. Further speed-ups result from cheaper auxiliary calculations and an often greatly improved convergence behaviour. Runtime to convergence improves by up to a factor of 2.5 on the Fermi-Hubbard model compared to the previous single-site method and by up to a factor of 3.9 compared to two-site DMRG. The method is compatible with real-space parallelisation and non-abelian symmetries.

preprint2015arXiv

How to discretize a quantum bath for real-time evolution

Many numerical techniques for the description of quantum systems that are coupled to a continuous bath require the discretization of the latter. To this end, a wealth of methods has been developed in the literature, which we classify as (i) direct discretization, (ii) orthogonal polynomial, and (iii) numerical optimization strategies. We recapitulate strategies (i) and (ii) to clarify their relation. For quadratic Hamiltonians, we show that (ii) is the best strategy in the sense that it gives the numerically exact time evolution up to a maximum time $t_\text{max}$, for which we give a simple expression. For non-quadratic Hamiltonians, we show that no such best strategy exists. We present numerical examples relevant to open quantum systems and strongly correlated systems, as treated by dynamical mean-field theory (DMFT).

preprint2015arXiv

Imaginary-time matrix product state impurity solver for dynamical mean-field theory

We present a new impurity solver for dynamical mean-field theory based on imaginary-time evolution of matrix product states. This converges the self-consistency loop on the imaginary-frequency axis and obtains real-frequency information in a final real-time evolution. Relative to computations on the real-frequency axis, required bath sizes are much smaller and less entanglement is generated, so much larger systems can be studied. The power of the method is demonstrated by solutions of a three band model in the single and two-site dynamical mean-field approximation. Technical issues are discussed, including details of the method, efficiency as compared to other matrix product state based impurity solvers, bath construction and its relation to real-frequency computations and the analytic continuation problem of quantum Monte Carlo, the choice of basis in dynamical cluster approximation, and perspectives for off-diagonal hybridization functions.

preprint2015arXiv

Lattice assisted spectroscopy: a generalized scanning tunnelling microscope for ultra-cold atoms

We show that the possibility to address and image single sites of an optical lattice, now an experimental reality, allows to measure the frequency-resolved local particle and hole spectra of a wide variety of one- and two-dimensional systems of lattice-confined strongly correlated ultracold atoms. Combining perturbation theory and time-dependent DMRG, we validate this scheme of lattice-assisted spectroscopy (LAS) on several example systems, such as the 1D superfluid and Mott insulator, with and without a parabolic trap, and finally on edge states of the bosonic Su-Schrieffer-Heeger model. We also highlight extensions of our basic scheme to obtain an even wider variety of interesting and important frequency resolved spectra.

preprint2015arXiv

Phase diagram of the J1-J2 Heisenberg model on the kagome lattice

We perform an extensive density matrix renormalization group (DMRG) study of the ground-state phase diagram of the spin-1/2 J_1-J_2 Heisenberg model on the kagome lattice. We focus on the region of the phase diagram around the kagome Heisenberg antiferromagnet, i.e., at J_2=0. We investigate the static spin structure factor, the magnetic correlation lengths, and the spin gaps. Our results are consistent with the absence of magnetic order in a narrow region around J_2\approx 0, although strong finite-size effects do not allow us to accurately determine the phase boundaries. This result is in agreement with the presence of an extended spin-liquid region, as it has been proposed recently. Outside the disordered region, we find that for ferromagnetic and antiferromagnetic J_2 the ground state displays signatures of the magnetic order of the \sqrt{3}\times\sqrt{3} and the q=0 type, respectively. Finally, we focus on the structure of the entanglement spectrum (ES) in the q=0 ordered phase. We discuss the importance of the choice of the bipartition on the finite-size structure of the ES.

preprint2015arXiv

Spectral functions and time evolution from the Chebyshev recursion

We link linear prediction of Chebyshev and Fourier expansions to analytic continuation. We push the resolution in the Chebyshev-based computation of $T=0$ many-body spectral functions to a much higher precision by deriving a modified Chebyshev series expansion that allows to reduce the expansion order by a factor $\sim\frac{1}{6}$. We show that in a certain limit the Chebyshev technique becomes equivalent to computing spectral functions via time evolution and subsequent Fourier transform. This introduces a novel recursive time evolution algorithm that instead of the group operator $e^{-iHt}$ only involves the action of the generator $H$. For quantum impurity problems, we introduce an adapted discretization scheme for the bath spectral function. We discuss the relevance of these results for matrix product state (MPS) based DMRG-type algorithms, and their use within dynamical mean-field theory (DMFT). We present strong evidence that the Chebyshev recursion extracts less spectral information from $H$ than time evolution algorithms when fixing a given amount of created entanglement.

preprint2014arXiv

Bound states and entanglement in the excited states of quantum spin chains

We investigate entanglement properties of the excited states of the spin-1/2 Heisenberg (XXX) chain with isotropic antiferromagnetic interactions, by exploiting the Bethe ansatz solution of the model. We consider eigenstates obtained from both real and complex solutions ("strings") of the Bethe equations. Physically, the former are states of interacting magnons, whereas the latter contain bound states of groups of particles. We first focus on the situation with few particles in the chain. Using exact results and semiclassical arguments, we derive an upper bound S_MAX for the entanglement entropy. This exhibits an intermediate behavior between logarithmic and extensive, and it is saturated for highly-entangled states. As a function of the eigenstate energy, the entanglement entropy is organized in bands. Their number depends on the number of blocks of contiguous Bethe-Takahashi quantum numbers. In presence of bound states a significant reduction in the entanglement entropy occurs, reflecting that a group of bound particles behaves effectively as a single particle. Interestingly, the associated entanglement spectrum shows edge-related levels. At finite particle density, the semiclassical bound S_MAX becomes inaccurate. For highly-entangled states S_A\propto L_c, with L_c the chord length, signaling the crossover to extensive entanglement. Finally, we consider eigenstates containing a single pair of bound particles. No significant entanglement reduction occurs, in contrast with the few-particle case.

preprint2014arXiv

Chebyshev Matrix Product State Impurity Solver for the Dynamical Mean-Field Theory

We compute the spectral functions for the two-site dynamical cluster theory and for the two-orbital dynamical mean-field theory in the density-matrix renormalization group (DMRG) framework using Chebyshev expansions represented with matrix product states (MPS). We obtain quantitatively precise results at modest computational effort through technical improvements regarding the truncation scheme and the Chebyshev rescaling procedure. We furthermore establish the relation of the Chebyshev iteration to real-time evolution, and discuss technical aspects as computation time and implementation in detail.

preprint2014arXiv

Competing regimes of motion in 1D mobile impurities

We show that a distinguishable mobile impurity inside a one-dimensional many-body state at zero temperature generally does not behave like a quasiparticle (QP). Instead, both the impurities dynamics as well as the ground state of the bath are fundamentally transformed by a diverging number of zero-energy excitations being generated, leading to what we call infrared-dominated (ID) dynamics. Combining analytics and DMRG numerics we provide a general formula for the power law governing ID dynamics at zero momentum, discuss a threshold beyond which quasiparticle dynamics may occur again, and study the competition between the ID and quasiparticle universality classes at larger impurity momenta.

preprint2014arXiv

Domain-wall melting in ultracold boson systems with holes and spin-flip defects

Quantum magnetism is a fundamental phenomenon of nature. As of late, it has garnered a lot of interest because experiments with ultracold atomic gases in optical lattices could be used as a simulator for phenomena of magnetic systems. A paradigmatic example is the time evolution of a domain-wall state of a spin-1/2 Heisenberg chain, the so-called domain-wall melting. The model can be implemented by having two species of bosonic atoms with unity filling and strong on-site repulsion U in an optical lattice. In this paper, we study the domain-wall melting in such a setup on the basis of the time-dependent density matrix renormalization group (tDMRG). We are particularly interested in the effects of defects that originate from an imperfect preparation of the initial state. Typical defects are holes (empty sites) and flipped spins. We show that the dominating effects of holes on observables like the spatially resolved magnetization can be taken account of by a linear combination of spatially shifted observables from the clean case. For sufficiently large U, further effects due to holes become negligible. In contrast, the effects of spin flips are more severe as their dynamics occur on the same time scale as that of the domain-wall melting itself. It is hence advisable to avoid preparation schemes that are based on spin-flips.

preprint2014arXiv

Quantum magnetism of bosons with synthetic gauge fields in one-dimensional optical lattices: a Density Matrix Renormalization Group study

In this paper, we provide a comprehensive study of the quantum magnetism in the Mott insulating phases of the 1D Bose-Hubbard model with abelian or non-abelian synthetic gauge fields, using the Density Matrix Renormalization Group (DMRG) method. We focus on the interplay between the synthetic gauge field and the asymmetry of the interactions, which give rise to a very general effective magnetic model: a XYZ model with various Dzyaloshinskii-Moriya (DM) interactions. The properties of the different quantum magnetic phases and phases transitions of this model are investigated.

preprint2014arXiv

Solving nonequilibrium dynamical mean-field theory using matrix product states

We solve the nonequilibrium dynamical mean-field theory (DMFT) using matrix product states (MPS). This allows us to treat much larger bath sizes and by that reach substantially longer times (factor $\sim$ 2 -- 3) than with exact diagonalization. We show that the star geometry of the underlying impurity problem can have substantially better entanglement properties than the previously favoured chain geometry. This has immense consequences for the efficiency of an MPS-based description of general impurity problems: in the case of equilibrium DMFT, it leads to an orders-of-magnitude speedup. We introduce an approximation for the two-time hybridization function that uses time-translational invariance, which can be observed after a certain relaxation time after a quench to a time-independent Hamiltonian.

preprint2012arXiv

Fractional Excitations in Cold Atomic Gases

We study the behavior of excitations in the tilted one-dimensional Bose-Hubbard model. In the phase with broken symmetry, fundamental excitations are domain-walls which show fractional statistics. Using perturbation theory, we derive an analytic model for the time evolution of these fractional excitations, and demonstrate the existence of a repulsively bound state above a critical center of mass momentum. The validity of the perturbative analysis is confirmed by the use of t- DMRG simulations. These findings open the path for experimental detection of fractional particles in cold atomic gases.

preprint2012arXiv

Quantum dynamics of a single, mobile spin impurity

Quantum magnetism describes the properties of many materials such as transition metal oxides and cuprate superconductors. One of its elementary processes is the propagation of spin excitations. Here we study the quantum dynamics of a deterministically created spin-impurity atom, as it propagates in a one-dimensional lattice system. We probe the full spatial probability distribution of the impurity at different times using single-site-resolved imaging of bosonic atoms in an optical lattice. In the Mott-insulating regime, a post-selection of the data allows to reduce the effect of temperature, giving access to a space- and time-resolved measurement of the quantum-coherent propagation of a magnetic excitation in the Heisenberg model. Extending the study to the bath's superfluid regime, we determine quantitatively how the bath strongly affects the motion of the impurity. The experimental data shows a remarkable agreement with theoretical predictions allowing us to determine the effect of temperature on the coherence and velocity of impurity motion. Our results pave the way for a new approach to study quantum magnetism, mobile impurities in quantum fluids, and polarons in lattice systems.

preprint2012arXiv

Scaling of the thermal spectral function for quantum critical bosons in one dimension

We present an improved scheme for the precise evaluation of finite-temperature response functions of strongly correlated systems in the framework of the time-dependent density matrix renormalization group. The maximum times that we can reach at finite temperatures T are typically increased by a factor of two, when compared against the earlier approaches. This novel scheme, complemented with linear prediction, allows us now to evaluate dynamic correlators for interacting bosons in one dimension. We demonstrate that the considered spectral function in the quantum critical regime with dynamic critical exponent z=2 is captured by the universal scaling form S(k,omega)=(1/T)*Phi(k/sqrt(T),omega/T) and calculate the scaling function precisely.

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

Landau-Zener sweeps and sudden quenches in coupled Bose-Hubbard chains

We simulate numerically the dynamics of strongly correlated bosons in a two-leg ladder subject to a time-dependent energy bias between the two chains. When all atoms are initially in the leg with higher energy, we find a drastic reduction of the inter-chain particle transfer for slow linear sweeps, in quantitative agreement with recent experiments. This effect is preceded by a rapid broadening of the quasi-momentum distribution of atoms, signaling the presence of a bath of low-energy excitations in the chains. We further investigate the scenario of quantum quenches to fixed values of the energy bias. We find that for large enough density the momentum distribution relaxes to that of an equilibrium thermal state with the same energy.

preprint2011arXiv

Probing the relaxation towards equilibrium in an isolated strongly correlated 1D Bose gas

The problem of how complex quantum systems eventually come to rest lies at the heart of statistical mechanics. The maximum entropy principle put forward in 1957 by E. T. Jaynes suggests what quantum states one should expect in equilibrium but does not hint as to how closed quantum many-body systems dynamically equilibrate. A number of theoretical and numerical studies accumulate evidence that under specific conditions quantum many-body models can relax to a situation that locally or with respect to certain observables appears as if the entire system had relaxed to a maximum entropy state. In this work, we report the experimental observation of the non-equilibrium dynamics of a density wave of ultracold bosonic atoms in an optical lattice in the regime of strong correlations. Using an optical superlattice, we are able to prepare the system in a well-known initial state with high fidelity. We then follow the dynamical evolution of the system in terms of quasi-local densities, currents, and coherences. Numerical studies based on the time-dependent density-matrix renormalization group method are in an excellent quantitative agreement with the experimental data. For very long times, all three local observables show a fast relaxation to equilibrium values compatible with those expected for a global maximum entropy state. We find this relaxation of the quasi-local densities and currents to initially follow a power-law with an exponent being significantly larger than for free or hardcore bosons. For intermediate times the system fulfills the promise of being a dynamical quantum simulator, in that the controlled dynamics runs for longer times than present classical algorithms based on matrix product states can efficiently keep track of.

preprint2010arXiv

Coulomb interaction effects and electron spin relaxation in the 1d Kondo lattice model

We study the effects of the Coulomb interaction in the one dimensional Kondo lattice model on the phase diagram, the static magnetic susceptibility and electron spin relaxation. We show that onsite Coulomb interaction supports ferromagnetic order and nearest neighbor Coulomb interaction drives, depending on the electron filling, either a paramagnetic or ferromagnetic order. Furthermore we calculate electron quasiparticle life times, which can be related to electron spin relaxation and decoherence times, and explain their dependence on the strength of interactions and the electron filling in order to find the sweet spot of parameters where the relaxation time is maximized. We find that effective exchange processes between the electrons dominate the spin relaxation and decoherence rate.