Source author record

Thomas Barthel

Thomas Barthel 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

15works
7topics
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

15 published item(s)

preprint2022arXiv

Eigenstate entanglement scaling for critical interacting spin chains

With increasing subsystem size and energy, bipartite entanglement entropies of energy eigenstates cross over from the groundstate scaling to a volume law. In previous work, we pointed out that, when strong or weak eigenstate thermalization (ETH) applies, the entanglement entropies of all or, respectively, almost all eigenstates follow a single crossover function. The crossover functions are determined by the subsystem entropy of thermal states and assume universal scaling forms in quantum-critical regimes. This was demonstrated by field-theoretical arguments and the analysis of large systems of non-interacting fermions and bosons. Here, we substantiate such scaling properties for integrable and non-integrable interacting spin-1/2 chains at criticality using exact diagonalization. In particular, we analyze XXZ and transverse-field Ising models with and without next-nearest-neighbor interactions. Indeed, the crossover of thermal subsystem entropies can be described by a universal scaling function following from conformal field theory. Furthermore, we analyze the validity of ETH for entanglement in these models. Even for the relatively small system sizes that can be simulated, the distributions of eigenstate entanglement entropies are sharply peaked around the subsystem entropies of the corresponding thermal ensembles.

preprint2022arXiv

On the closedness and geometry of tensor network state sets

Tensor network states (TNS) are a powerful approach for the study of strongly correlated quantum matter. The curse of dimensionality is addressed by parametrizing the many-body state in terms of a network of partially contracted tensors. These tensors form a substantially reduced set of effective degrees of freedom. In practical algorithms, functionals like energy expectation values or overlaps are optimized over certain sets of TNS. Concerning algorithmic stability, it is important whether the considered sets are closed because, otherwise, the algorithms may approach a boundary point that is outside the TNS set and tensor elements diverge. We discuss the closedness and geometries of TNS sets, and we propose regularizations for optimization problems on non-closed TNS sets. We show that sets of matrix product states (MPS) with open boundary conditions, tree tensor network states (TTNS), and the multiscale entanglement renormalization ansatz (MERA) are always closed, whereas sets of translation-invariant MPS with periodic boundary conditions (PBC), heterogeneous MPS with PBC, and projected entangled-pair states (PEPS) are generally not closed. The latter is done using explicit examples like the W state, states that we call two-domain states, and fine-grained versions thereof.

preprint2022arXiv

Tensor Network States with Low-Rank Tensors

Tensor networks are used to efficiently approximate states of strongly-correlated quantum many-body systems. More generally, tensor network approximations may allow to reduce the costs for operating on an order-$N$ tensor from exponential to polynomial in $N$, and this has become a popular approach for machine learning. We introduce the idea of imposing low-rank constraints on the tensors that compose the tensor network. With this modification, the time and space complexities for the network optimization can be substantially reduced while maintaining high accuracy. We detail this idea for tree tensor network states (TTNS) and projected entangled-pair states. Simulations of spin models on Cayley trees with low-rank TTNS exemplify the effect of rank constraints on the expressive power. We find that choosing the tensor rank $r$ to be on the order of the bond dimension $m$, is sufficient to obtain high-accuracy groundstate approximations and to substantially outperform standard TTNS computations. Thus low-rank tensor networks are a promising route for the simulation of quantum matter and machine learning on large data sets.

preprint2020arXiv

The matrix product approximation for the dynamic cavity method

Stochastic dynamics of classical degrees of freedom, defined on vertices of locally tree-like graphs, can be studied in the framework of the dynamic cavity method which is exact for tree graphs. Such models correspond for example to spin-glass systems, Boolean networks, neural networks, and other technical, biological, and social networks. The central objects in the cavity method are edge messages -- conditional probabilities of two vertex variable trajectories. In this paper, we discuss a rather pedagogical derivation for the dynamic cavity method, give a detailed account of the novel matrix product edge message (MPEM) algorithm for the solution of the dynamic cavity equation as introduced in Phys. Rev. E 97, 010104(R) (2018), and present optimizations and extensions. Matrix product approximations of the edge messages are constructed recursively in an iteration over time. Computation costs and precision can be tuned by controlling the matrix dimensions of the MPEM in truncations. Without truncations, the dynamics is exact. Data for Glauber-Ising dynamics shows a linear growth of computation costs in time. In contrast to Monte Carlo simulations, the approach has a much better error scaling. Hence, it gives for example access to low probability events and decaying observables like temporal correlations. We discuss optimized truncation schemes and an extension that allows to capture models which have a continuum time limit.

preprint2017arXiv

Minimally entangled typical thermal states versus matrix product purifications for the simulation of equilibrium states and time evolution

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 the 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.

preprint2015arXiv

Finite-temperature effects on interacting bosonic 1D systems in disordered lattices

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.

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

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.

preprint2013arXiv

Algebraic versus exponential decoherence in dissipative many-particle systems

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.

preprint2013arXiv

Precise evaluation of thermal response functions by optimized density matrix renormalization group schemes

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 for 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.

preprint2012arXiv

Quasi-locality and efficient simulation of Markovian quantum dynamics

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 quasi-local, 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, the 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, our result can be used to rigorously bound finite-size effects.

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.

preprint2012arXiv

Solving condensed-matter ground-state problems by semidefinite relaxations

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.

preprint2010arXiv

Contraction of fermionic operator circuits and the simulation of strongly correlated fermions

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 ansatz (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.

preprint2010arXiv

Real-space renormalization yields finite correlations

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 spatial dimension, MERA states are actually finitely correlated states, i.e., projected entangled pair states (PEPS) with a bond dimension independent of the system size. Hence, real-space renormalization generates states which can be encoded with local effective degrees of freedom, and MERA states form an efficiently contractible class of PEPS that obey the area law for the entanglement entropy. It is shown further that there exist other efficiently contractible schemes violating the area law.