"Infinite boundary conditions for response functions and limit cycles in iDMRG, demonstrated for bilinear-biquadratic spin-1 chains"
M. Binder and T. Barthel arXiv:1804.09163,
Response functions <A(x,t)B(y,0)> for one-dimensional strongly correlated quantum many-body systems can be computed with matrix product state (MPS) techniques. Especially, when one is interested in spectral functions or dynamic structure factors of translation-invariant systems, the response for some range |x-y|<L is needed. We demonstrate how the number of required time-evolution runs can be reduced substantially: (a) If finite-system simulations are employed, the number of time-evolution runs can be reduced from L to 2sqrt(L). (b) To go beyond, one can employ infinite MPS (iMPS) such that two evolution runs suffice. To this purpose, iMPS that are heterogeneous only around the causal cone of the perturbation are evolved in time, i.e., the simulation is done with so-called infinite boundary conditions. Computing overlaps of these states, spatially shifted relative to each other, yields the response functions for all distances |x-y|. As a specific application, we compute the dynamic structure factor for ground states of bilinear-biquadratic spin-1 chains with very high resolution and explain the underlying low-energy physics. To determine the initial uniform iMPS for such simulations, infinite-system density-matrix renormalization group (iDMRG) can be employed. We discuss that, depending on the system and chosen bond dimension, iDMRG with a cell size n may converge to a non-trivial limit cycle of length m. This then corresponds to an iMPS with an enlarged unit cell of size mn.
"Fundamental limitations for measurements in quantum many-body systems"
T. Barthel and J. Lu arXiv:1802.04378,
Dynamical measurement schemes are an important tool for the investigation of quantum many-body systems, especially in the age of quantum simulation. Here, we address the question whether generic measurements can be implemented efficiently if we have access to a certain set of experimentally realizable measurements and can extend it through time evolution. For the latter, two scenarios are considered (a) evolution according to unitary circuits and (b) evolution due to Hamiltonians that we can control in a time-dependent fashion. We find that the time needed to realize a certain measurement to a predefined accuracy scales in general exponentially with the system size -- posing a fundamental limitation. The argument is based, on the construction of epsilon-packings for manifolds of observables with identical spectra and a comparison of their cardinalities to those of epsilon-coverings for quantum circuits and unitary time-evolution operators. The former is related to the study of Grassmann manifolds.
"Typical one-dimensional quantum systems at finite temperatures can be simulated efficiently on classical computers"
T. Barthel arXiv:1708.09349,
It is by now well-known that ground states of gapped one-dimensional (1d) quantum-many body systems with short-range interactions can be studied efficiently using classical computers and matrix product state techniques. A corresponding result for finite temperatures was missing. For 1d systems that can be described by a local 1+1d field theory, it is shown here that the cost for the classical simulation at finite temperatures grows in fact only polynomially with the inverse temperature and is system-size independent -- even for quantum critical systems. In particular, we show that the thermofield double state (TDS), a purification of the equilibrium density operator, can be obtained efficiently in matrix-product form. The argument is based on the scaling behavior of Rényi entanglement entropies in the TDS. At finite temperatures, they obey the area law. For quantum critical systems, the entanglement is found to grow only logarithmically with inverse temperature, S~log(beta). The field-theoretical results are confirmed by quasi-exact numerical simulations of quantum magnets and interacting bosons.
"Symmetric minimally entangled typical thermal states for canonical and grand-canonical ensembles"
M. Binder and T. Barthel arXiv:1701.03872,
pdf, Phys. Rev. B 95, 195148 (2017)
Based on the density matrix renormalization group (DMRG), strongly correlated quantum many-body systems at finite temperatures can be simulated by sampling over a certain class of pure matrix product states (MPS) called minimally entangled typical thermal states (METTS). When a system features symmetries, these can be utilized to substantially reduce MPS computation costs. It is conceptually straightforward to simulate canonical ensembles using symmetric METTS. In practice, it is important to alternate between different symmetric collapse bases to decrease autocorrelations in the Markov chain of METTS. To this purpose, we introduce symmetric Fourier and Haar-random block bases that are efficiently mixing. We also show how grand-canonical ensembles can be simulated efficiently with symmetric METTS. We demonstrate these approaches for spin-1/2 XXZ chains and discuss how the choice of the collapse bases influences autocorrelations as well as the distribution of measurement values and, hence, convergence speeds.
"Matrix product purifications for canonical ensembles and quantum number distributions"
T. Barthel arXiv:1607.01696,
pdf, Phys. Rev. B 94, 115157 (2016)
Matrix product purifications (MPPs) are a very efficient tool for the simulation of strongly correlated quantum many-body systems at finite temperatures. When a system features symmetries, these can be used to reduce computation costs substantially. It is straightforward to compute an MPP of a grand-canonical ensemble, also when symmetries are exploited. This paper provides and demonstrates methods for the efficient computation of MPPs of canonical ensembles under utilization of symmetries. Furthermore, we present a scheme for the evaluation of global quantum number distributions using matrix product density operators (MPDOs). We provide exact matrix product representations for canonical infinite-temperature states, and discuss how they can be constructed alternatively by applying matrix product operators to vacuum-type states or by using entangler Hamiltonians. A demonstration of the techniques for Heisenberg spin-1/2 chains explains why the difference in the energy densities of canonical and grand-canonical ensembles decays as 1/L.
"Finite-temperature effects on interacting bosonic one-dimensional systems in disordered lattices"
L. Gori, T. Barthel, A. Kumar, E. Lucioni, L. Tanzi, M. Inguscio, G. Modugno, T. Giamarchi, C. D'Errico, and G. Roux arXiv:1512.04238,
pdf, Phys. Rev. A 93, 033650 (2016)
We analyze the finite-temperature effects on the phase diagram describing the insulating properties of interacting 1D bosons in a quasi-periodic lattice. We examine thermal effects by comparing experimental results to exact diagonalization for small-sized systems and to density-matrix renormalization group (DMRG) computations. At weak interactions, we find short thermal correlation lengths, indicating a substantial impact of temperature on the system coherence. Conversely, at strong interactions, the obtained thermal correlation lengths are significantly larger than the localization length, and the quantum nature of the T=0 Bose glass phase is preserved up to a crossover temperature that depends on the disorder strength. Furthermore, in the absence of disorder, we show how quasi-exact finite-T DMRG computations, compared to experimental results, can be employed to estimate the temperature, which is not directly accessible in the experiment.
"A matrix product algorithm for stochastic dynamics on networks, applied to non-equilibrium Glauber dynamics"
T. Barthel, C. De Bacco, and S. Franz arXiv:1508.03295,
pdf, Phys. Rev. E 97, 010104(R) (2018)
We introduce and apply a novel efficient method for the precise simulation of stochastic dynamical processes on locally tree-like graphs. Networks with cycles are treated in the framework of the cavity method. Such models correspond, for example, to spin-glass systems, Boolean networks, neural networks, or other technological, biological, and social networks. Building upon ideas from quantum many-body theory, the new approach is based on a matrix product approximation of the so-called edge messages -- conditional probabilities of vertex variable trajectories. Computation costs and accuracy can be tuned by controlling the matrix dimensions of the matrix product edge messages (MPEM) in truncations. In contrast to Monte Carlo simulations, the algorithm has a better error scaling and works for both, single instances as well as the thermodynamic limit. We employ it to examine prototypical non-equilibrium Glauber dynamics in the kinetic Ising model. Because of the absence of cancellation effects, observables with small expectation values can be evaluated accurately, allowing for the study of decay processes and temporal correlations.
"Phase Diagram of an Extended Quantum Dimer Model on the Hexagonal Lattice"
T. M. Schlittler, T. Barthel, G. Misguich, J. Vidal, and R. Mosseri arXiv:1507.04643,
pdf, Phys. Rev. Lett. 115, 217202 (2015)
We introduce a quantum dimer model on the hexagonal lattice that, in addition to the standard three-dimer kinetic and potential terms, includes a competing potential part counting dimer-free hexagons. The zero-temperature phase diagram is studied by means of quantum Monte Carlo simulations, supplemented by variational arguments. It reveals some new crystalline phases and a cascade of transitions with rapidly changing flux (tilt in the height language). We analyze perturbatively the vicinity of the Rokhsar-Kivelson point, showing that this model has the microscopic ingredients needed for the devil's staircase scenario [Fradkin et al. Phys. Rev. B 69, 224415 (2004)], and is therefore expected to produce fractal variations of the ground-state flux.
"Phase diagram of the hexagonal lattice quantum dimer model: Order parameters, ground-state energy, and gaps"
T. M. Schlittler, R. Mosseri, and T. Barthel arXiv:1501.02242,
pdf, Phys. Rev. B 96, 195142 (2017)
The phase diagram of the quantum dimer model on the hexagonal (honeycomb) lattice is computed numerically, extending on earlier work by Moessner, Sondhi, and Chandra. The different ground state phases are studied in detail using several local and global observables. In addition, we analyze imaginary-time correlation functions to determine ground state energies as well as gaps to the first excited states. This leads in particular to a confirmation that the intermediary so-called plaquette phase is gapped -- a point which was previously advocated with general arguments and some data for an order parameter, but required a more direct proof. On the technical side, we describe an efficient world-line quantum Monte Carlo algorithm with improved cluster updates that increase acceptance probabilities by taking account of potential terms of the Hamiltonian during the cluster construction. The Monte Carlo simulations are supplemented with variational computations.
"Minimally entangled typical thermal states versus matrix product purifications for the simulation of equilibrium states and time evolution"
M. Binder and T. Barthel arXiv:1411.3033,
pdf, Phys. Rev. B 92, 125119 (2015)
For the simulation of equilibrium states and finite-temperature response functions of strongly-correlated quantum many-body systems, we compare the efficiencies of two different approaches in the framework of the density matrix renormalization group (DMRG). The first is based on matrix product purifications. The second, more recent one, is based on so-called minimally entangled typical thermal states (METTS). For the latter, we highlight the interplay of statistical and DMRG truncation errors, discuss the use of self-averaging effects, and describe schemes for the computation of response functions. For critical as well as gapped phases of the spin-1/2 XXZ chain and the one-dimensional Bose-Hubbard model, we assess computation costs and accuracies of the two methods at different temperatures. For almost all considered cases, we find that, for the same computation cost, purifications yield more accurate results than METTS -- often by orders of magnitude. The METTS algorithm becomes more efficient only for temperatures well below the system's energy gap. The exponential growth of the computation cost in the evaluation of response functions limits the attainable timescales in both methods and we find that in this regard, METTS do not outperform purifications.
"Bound states and entanglement in the excited states of quantum spin chains"
J. Mölter, T. Barthel, U. Schollwöck, and V. Alba arXiv:1407.0066,
pdf, J. Stat. Mech. P10029 (2014), Special Issue "Quantum Entanglement in Condensed Matter Physics"J. Stat. Mech. 2014
We investigate the 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 Smax for the entanglement entropy. This exhibits an intermediate behaviour 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 the 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 a finite particle density, the semiclassical bound Smax becomes inaccurate. For highly-entangled states the subsystem entropy is proprtional to the chord length, signalling 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.
"Multispinon continua at zero and finite temperature in a near-ideal Heisenberg chain"
B. Lake, D. A. Tennant, J.-S. Caux, T. Barthel, U. Schollwöck, S. E. Nagler, and C. D. Frost arXiv:1307.4071,
pdf, Phys. Rev. Lett. 111, 137205 (2013)
The space- and time-dependent response of many-body quantum systems is the most informative aspect of their emergent behaviour. The dynamical structure factor, experimentally measurable using neutron scattering, can map this response in wavevector and energy with great detail, allowing theories to be quantitatively tested to high accuracy. Here, we present a comparison between neutron scattering measurements on the one-dimensional spin-1/2 Heisenberg antiferromagnet KCuF3, and recent state-of-the-art theoretical methods based on integrability and density matrix renormalization group simulations. The unprecedented quantitative agreement shows that precise descriptions of strongly correlated states at all distance, time and temperature scales are now possible, and highlights the need to apply these novel techniques to other problems in low-dimensional magnetism.
"Domain-wall melting in ultracold-boson systems with hole and spin-flip defects"
J. C. Halimeh, A. Wöllert, I. P. McCulloch, U. Schollwöck, and T. Barthel arXiv:1307.0513,
pdf, Phys. Rev. A 89, 063603 (2014)
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.
"Algebraic versus exponential decoherence in dissipative many-particle systems"
Z. Cai and T. Barthel arXiv:1304.6890,
pdf, Phys. Rev. Lett. 111, 150403 (2013), covered by Phys.org "Quantum particles find safety in numbers"phys.org, 10/16/2013
The interplay between dissipation and internal interactions in quantum many-body systems gives rise to a wealth of novel phenomena. Here we investigate spin-1/2 chains with uniform local couplings to a Markovian environment using the time-dependent density matrix renormalization group (tDMRG). For the open XXZ model, we discover that the decoherence time diverges in the thermodynamic limit. The coherence decay is then algebraic instead of exponential. This is due to a vanishing gap in the spectrum of the corresponding Liouville superoperator and can be explained on the basis of a perturbative treatment. In contrast, decoherence in the open transverse-field Ising model is found to be always exponential. In this case, the internal interactions can both facilitate and impede the environment-induced decoherence.
"Precise evaluation of thermal response functions by optimized density matrix renormalization group schemes"
T. Barthel arXiv:1301.2246,
pdf, New J. Phys. 15, 073010 (2013)
This paper provides a study and discussion of earlier as well as novel more efficient schemes for the precise evaluation of finite-temperature response functions of strongly correlated quantum systems in the framework of the time-dependent density matrix renormalization group (tDMRG). The computational costs and bond dimensions as functions of time and temperature are examined at the example of the spin-1/2 XXZ Heisenberg chain in the critical XY phase and the gapped Néel phase. The matrix product state purifications occurring in the algorithms are in one-to-one relation with corresponding matrix product operators. This notational simplification elucidates implications of quasi-locality on the computational costs. Based on the observation that there is considerable freedom in designing efficient tDMRG schemes for the calculation of dynamical correlators at finite temperatures, a new class of optimizable schemes, as recently suggested in arXiv:1212.3570, is explained and analyzed numerically. A specific novel near-optimal scheme that requires no additional optimization reaches maximum times that are typically increased by a factor of two, when compared against earlier approaches. These increased reachable times make many more physical applications accessible. For each of the described tDMRG schemes, one can devise a corresponding transfer matrix renormalization group (TMRG) variant.
"Scaling of the thermal spectral function for quantum critical bosons in one dimension"
T. Barthel, U. Schollwöck, and S. Sachdev arXiv:1212.3570,
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.
"Quasi-locality and efficient simulation of Markovian quantum dynamics"
T. Barthel and M. Kliesch arXiv:1111.4210,
pdf, Phys. Rev. Lett. 108, 230504 (2012)
We consider open many-body systems governed by a time-dependent quantum master equation with short-range interactions. With a generalized Lieb-Robinson bound, we show that the evolution in this very generic framework is quasilocal; i.e., the evolution of observables can be approximated by implementing the dynamics only in a vicinity of the observables' support. The precision increases exponentially with the diameter of the considered subsystem. Hence, time evolution can be simulated on classical computers with a cost that is independent of the system size. Providing error bounds for Trotter decompositions, we conclude that the simulation on a quantum computer is additionally efficient in time. For experiments and simulations in the Schrödinger picture, our result can be used to rigorously bound finite-size effects.
"Solving condensed matter ground state problems by semidefinite relaxations"
T. Barthel and R. Hübener arXiv:1106.4966,
pdf, Phys. Rev. Lett. 108, 200404 (2012)
We present a new generic approach to the condensed-matter ground-state problem which is complementary to variational techniques and works directly in the thermodynamic limit. Relaxing the ground-state problem, we obtain semidefinite programs (SDP). These can be solved efficiently, yielding strict lower bounds to the ground-state energy and approximations to the few-particle Green's functions. As the method is applicable for all particle statistics, it represents in particular a novel route for the study of strongly correlated fermionic and frustrated spin systems in D>1 spatial dimensions. It is demonstrated for the XXZ model and the Hubbard model of spinless fermions. The results are compared against exact solutions, quantum Monte Carlo, and Anderson bounds, showing the competitiveness of the SDP method.
"Dissipative quantum Church-Turing theorem"
M. Kliesch, T. Barthel, C. Gogolin, M. Kastoryano, and J. Eisert arXiv:1105.3986,
pdf, Phys. Rev. Lett. 107, 120501 (2011), see also D. Browne "Viewpoint: Quantum simulation hits the open road"Physics 4, 72 (2011)
We show that the time evolution of an open quantum system, described by a possibly time dependent Liouvillian, can be simulated by a unitary quantum circuit of a size scaling polynomially in the simulation time and the size of the system. An immediate consequence is that dissipative quantum computing is no more powerful than the unitary circuit model. Our result can be seen as a dissipative Church-Turing theorem, since it implies that under natural assumptions, such as weak coupling to an environment, the dynamics of an open quantum system can be simulated efficiently on a quantum computer. Formally, we introduce a Trotter decomposition for Liouvillian dynamics and give explicit error bounds. This constitutes a practical tool for numerical simulations, e.g., using matrix-product operators. We also demonstrate that most quantum states cannot be prepared efficiently.
"Real-space renormalization yields finitely correlated states"
T. Barthel, M. Kliesch, and J. Eisert arXiv:1003.2319,
pdf, Phys. Rev. Lett. 105, 010502 (2010)
Real-space renormalization approaches for quantum lattice systems generate certain hierarchical classes of states that are subsumed by the multi-scale entanglement renormalization ansatz (MERA). It is shown that, with the exception of one dimension, MERA states can be efficiently mapped to finitely correlated states, also known as projected entangled pair states (PEPS), with a bond dimension independent of the system size. Hence, MERA states form an efficiently contractible class of PEPS and obey an area law for the entanglement entropy. It is shown further that there exist other efficiently contractible schemes violating the area law.
"Contraction of fermionic operator circuits and the simulation of strongly correlated fermions"
T. Barthel, C. Pineda, and J. Eisert arXiv:0907.3689,
pdf, Phys. Rev. A 80, 042333 (2009), also in Virtual Journal of Nanoscale Science and Technology20, Issue 20 (2009), and in Virtual Journal of Quantum Information9, Issue 11 (2009)
A fermionic operator circuit is a product of fermionic operators of usually different and partially overlapping support. Further elements of fermionic operator circuits (FOCs) are partial traces and partial projections. The presented framework allows for the introduction of fermionic versions of known qudit operator circuits (QUOC), important for the simulation of strongly correlated d-dimensional systems: the multiscale entanglement renormalization ansaetze (MERA), tree tensor networks (TTN), projected entangled pair states (PEPS), or their infinite-size versions (iPEPS etc.). After the definition of a FOC, we present a method to contract it with the same computation and memory requirements as a corresponding QUOC, for which all fermionic operators are replaced by qudit operators of identical dimension. A given scheme for contracting the QUOC relates to an analogous scheme for the corresponding fermionic circuit, where additional marginal computational costs arise only from reordering of modes for operators occurring in intermediate stages of the contraction. Our result hence generalizes efficient schemes for the simulation of d-dimensional spin systems, as MERA, TTN, or PEPS to the fermionic case.
"Unitary circuits for strongly correlated fermions"
C. Pineda, T. Barthel, and J. Eisert arXiv:0905.0669,
pdf, Phys. Rev. A 81, 050303(R) (2010)
We introduce a scheme for efficiently describing pure states of strongly correlated fermions in higher dimensions using unitary circuits. A local way of computing local expectation values is presented. We formulate a dynamical reordering scheme, corresponding to time-adaptive Jordan-Wigner transformation that avoids non-local string operators and only keeps suitably ordered the causal cone. Primitives of such a reordering scheme are highlighted. Fermionic unitary circuits can be contracted with the same complexity as in the spin case. The scheme gives rise to a variational description of fermionic models that does not suffer from a sign problem. We present a numerical example on a 9x9 fermionic lattice model to show the functioning of the approach.
"Magnetism, coherent many-particle dynamics, and relaxation with ultracold bosons in optical superlattices"
T. Barthel, C. Kasztelan, I. P. McCulloch, and U. Schollwöck arXiv:0809.5141,
pdf, Phys. Rev. A 79, 053627 (2009)
We study how well magnetic models can be implemented with ultracold bosonic atoms of two different hyperfine states in an optical superlattice. The system is captured by a two-species Bose-Hubbard model, but realizes in a certain parameter regime actually the physics of a spin-1/2 Heisenberg magnet, describing the second-order hopping processes. Tuning of the superlattice allows for controlling the effect of fast first-order processes versus the slower second-order ones. Using the density-matrix renormalization-group method, we provide the evolution of typical experimentally available observables. The validity of the description via the Heisenberg model, depending on the parameters of the Hubbard model, is studied numerically and analytically. The analysis is also motivated by recent experiments [S. Foelling et al., Nature (London) 448, 1029 (2007); S. Trotzky et al., Science 319, 295 (2008)] where coherent two-particle dynamics with ultracold bosonic atoms in isolated double wells were realized. We provide theoretical background for the next step, the observation of coherent many-particle dynamics after coupling the double wells. Contrary to the case of isolated double wells, relaxation of local observables can be observed. The tunability between the Bose-Hubbard model and the Heisenberg model in this setup could be used to study experimentally the differences in equilibration processes for nonintegrable and Bethe ansatz integrable models. We show that the relaxation in the Heisenberg model is connected to a phase averaging effect, which is in contrast to the typical scattering driven thermalization in nonintegrable models. We discuss the preparation of magnetic ground states by adiabatic tuning of the superlattice parameters.
"Quasiperiodic Bose-Hubbard model and localization in one-dimensional cold atomic gases"
G. Roux, T. Barthel, I. P. McCulloch, C. Kollath, U. Schollwöck, and T. Giamarchi arXiv:0802.3774,
pdf, Phys. Rev. A 78, 023628 (2008)
We compute the phase diagram of the one-dimensional Bose-Hubbard model with a quasiperiodic potential by means of the density-matrix renormalization group technique. This model describes the physics of cold atoms loaded in an optical lattice in the presence of a superlattice potential whose wavelength is incommensurate with the main lattice wavelength. After discussing the conditions under which the model can be realized experimentally, the study of the density vs the chemical potential curves for a nontrapped system unveils the existence of gapped phases at incommensurate densities interpreted as incommensurate charge-density-wave phases. Furthermore, a localization transition is known to occur above a critical value of the potential depth V_2 in the case of free and hard-core bosons. We extend these results to soft-core bosons for which the phase diagrams at fixed densities display new features compared with the phase diagrams known for random box distribution disorder. In particular, a direct transition from the superfluid phase to the Mott-insulating phase is found at finite V_2. Evidence for reentrances of the superfluid phase upon increasing interactions is presented. We finally comment on different ways to probe the emergent quantum phases and most importantly, the existence of a critical value for the localization transition. The latter feature can be investigated by looking at the expansion of the cloud after releasing the trap.
"Dephasing and the steady state in quantum many-particle systems"
T. Barthel and U. Schollwöck arXiv:0711.4896,
pdf, Phys. Rev. Lett. 100, 100601 (2008), also in Virtual Journal of Quantum Information8, Issue 3 (2008)
We discuss relaxation in bosonic and fermionic many-particle systems. For integrable systems, time evolution can cause a dephasing effect, leading for finite subsystems to steady states. We explicitly derive those steady subsystem states and devise sufficient prerequisites for the dephasing to occur. We also find simple scenarios, in which dephasing is ineffective and discuss the dependence on dimensionality and criticality. It follows further that, after a quench of system parameters, entanglement entropy will become extensive. This provides a way of creating strong entanglement in a controlled fashion.
"Entanglement entropy in collective models"
J. Vidal, S. Dusuel, and T. Barthel arXiv:cond-mat/0610833,
pdf, J. Stat. Mech. P01015 (2006)
We discuss the behaviour of the entanglement entropy of the ground state in various collective systems. Results for general quadratic two-mode boson models are given, yielding the relation between quantum phase transitions of the system (signalled by a divergence of the entanglement entropy) and the excitation energies. Such systems naturally arise when expanding collective spin Hamiltonians at leading order via the Holstein-Primakoff mapping. In a second step, we analyse several such models (the Dicke model, the two-level Bardeen-Cooper-Schrieffer model, the Lieb-Mattis model and the Lipkin-Meshkov-Glick model) and investigate the properties of the entanglement entropy over the whole parameter range. We show that when the system contains gapless excitations the entanglement entropy of the ground state diverges with increasing system size. We derive and classify the scaling behaviours that can be met.
"Entanglement entropy beyond the free case"
T. Barthel, S. Dusuel, and J. Vidal arXiv:cond-mat/0606436,
pdf, Phys. Rev. Lett. 97, 220402 (2006), also in Virtual Journal of Nanoscale Science and Technology14, Issue 24 (2006)
We present a perturbative method to compute the ground state entanglement entropy for interacting systems. We apply it to a collective model of mutually interacting spins in a magnetic field. At the quantum critical point, the entanglement entropy scales logarithmically with the subsystem size, the system size, and the anisotropy parameter. We determine the corresponding scaling prefactors and evaluate the leading finite-size correction to the entropy. Our analytical predictions are in perfect agreement with numerical results.
"Entanglement scaling in critical two-dimensional fermionic and bosonic systems"
T. Barthel, M.-C. Chung, and U. Schollwöck arXiv:cond-mat/0602077,
pdf, Phys. Rev. A 74, 022329 (2006)
We relate the reduced density matrices of quadratic fermionic and bosonic models to their Green's function matrices in a unified way and calculate the scaling of the entanglement entropy of finite systems in an infinite universe exactly. For critical fermionic two-dimensional (2D) systems at $T=0$, two regimes of scaling are identified: generically, we find a logarithmic correction to the area law with a prefactor dependence on the chemical potential that confirms earlier predictions based on the Widom conjecture. If, however, the Fermi surface of the critical system is zero-dimensional, then we find an area law with a sublogarithmic correction. For a critical bosonic 2D array of coupled oscillators at $T=0$, our results show that the entanglement entropy follows the area law without corrections.
"Entanglement and boundary critical phenomena"
H.-Q. Zhou, T. Barthel, J. O. Fjærestad, and U. Schollwöck arXiv:cond-mat/0511732,
pdf, Phys. Rev. A 74, 050305(R) (2006), also in Virtual Journal of Nanoscale Science and Technology14, Issue 23 (2006)
We investigate boundary critical phenomena from a quantum information perspective. Bipartite entanglement in the ground state of one-dimensional quantum systems is quantified using the Renyi entropy S_alpha, which includes the von Neumann entropy (alpha=1) and the single-copy entanglement (alpha=infinity) as special cases. We identify the contribution from the boundary entropy to the Renyi entropy, and show that there is an entanglement loss along boundary renormalization group (RG) flows. This property, which is intimately related to the Affleck-Ludwig g-theorem, can be regarded as a consequence of majorization relations between the spectra of the reduced density matrix along the boundary RG flows. We also point out that the bulk contribution to the single-copy entanglement is half of that to the von Neumann entropy, whereas the boundary contribution is the same.
Evolution of quantum many-body systems using tensor network states --
College Park MD 06/2018
Classical simulation of quantum matter and its limits --
Durham NC 05/2018
Fractal structures in the phase diagram of the honeycomb-lattice quantum dimer model --
Durham NC 12/2017
Typical 1d quantum systems at finite temperatures can be simulated efficiently on classical computers -- Minsk 11/2017, Cambridge MA 11/2017, Benasque 02/2018, Los Angeles 03/2018
Entanglement and computational complexity for 1D quantum systems at finite temperatures -- Munich 06/2017, Orsay 07/2017, Paris 07/2017
Entanglement and computational complexity for 1D quantum many-body systems -- Durham 03/2017, Wilmington 04/2017
A novel matrix product algorithm for stochastic dynamics on networks -- Baltimore 03/2016
Accurate T>0 response functions for strongly-correlated quasi-1D systems using DMRG -- Durham NC 09/2015
Devil's staircase in a quantum dimer model on the hexagonal lattice --
Strongly-correlated systems at finite T: METTS versus purifications --
Quantum dimers on the honeycomb lattice, phase diagram with devil's staircase --
Finite-temperature DMRG simulations of the Bose-Hubbard model --
Efficient parametrization and simulation of condensed matter systems using tensor network states --
Durham NC 04/2014
Simulating condensed matter systems with tensor network states and discovery of algebraic decoherence --
Charlottesville 01/2014, Lyon 03/2014, Saclay 03/2014
Tensor network state techniques for condensed matter systems and applications --
Algebraic versus exponential decoherence in dissipative many-particle systems -- Garching 10/2013, Orsay 11/2013, Paris 11/2013, Augsburg 01/2014, Houston 02/2014, Toulouse 04/2014, Durham NC 04/2014