BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Vincenzo Alba

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

ResearcherUnclaimed source record
cond-mat.stat-mechcond-mat.str-elhep-thquant-phcond-mat.quant-gascond-mat.dis-nnmath-phmath.MP

Catalog footprint

What is connected

31works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

31 published item(s)

preprint2022arXiv

Growth of Rényi Entropies in Interacting Integrable Models and the Breakdown of the Quasiparticle Picture

Rényi entropies are conceptually valuable and experimentally relevant generalisations of the celebrated von Neumann entanglement entropy. After a quantum quench in a clean quantum many-body system they generically display a universal linear growth in time followed by saturation. While a finite subsystem is essentially at local equilibrium when the entanglement saturates, it is genuinely out-of-equilibrium in the growth phase. In particular, the slope of the growth carries vital information on the nature of the system's dynamics, and its characterisation is a key objective of current research. Here we show that the slope of Rényi entropies can be determined by means of a spacetime duality transformation. In essence, we argue that the slope coincides with the stationary density of entropy of the model obtained by exchanging the roles of space and time. Therefore, very surprisingly, the slope of the entanglement is expressed as an equilibrium quantity. We use this observation to find an explicit exact formula for the slope of Rényi entropies in all integrable models treatable by thermodynamic Bethe ansatz and evolving from integrable initial states. Interestingly, this formula can be understood in terms of a quasiparticle picture only in the von Neumann limit.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Hydrodynamics of quantum entropies in Ising chains with linear dissipation

We study the dynamics of quantum information and of quantum correlations after a quantum quench, in transverse field Ising chains subject to generic linear dissipation. As we show, in the hydrodynamic limit of long times, large system sizes, and weak dissipation, entropy-related quantities -- such as the von Neumann entropy, the Rényi entropies, and the associated mutual information -- admit a simple description within the so-called quasiparticle picture. Specifically, we analytically derive a hydrodynamic formula, recently conjectured for generic noninteracting systems, which allows us to demonstrate a universal feature of the dynamics of correlations in such dissipative noninteracting system. For any possible dissipation, the mutual information grows up to a time scale that is proportional to the inverse dissipation rate, and then decreases, always vanishing in the long time limit. In passing, we provide analytic formulas describing the time-dependence of arbitrary functions of the fermionic covariance matrix, in the hydrodynamic limit.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Noninteracting fermionic systems with localized losses: Exact results in the hydrodynamic limit

We investigate the interplay between unitary dynamics after a quantum quench and localized dissipation in a noninteracting fermionic chain. In particular, we consider the effect of gain and loss processes, for which fermions are added and removed incoherently. We focus on the hydrodynamic limit of large distances from the source of dissipation and of long times, with their ratio being fixed. In this limit, the localized dissipation gives rise to an effective delta potential (dissipative impurity), and the time-evolution of the local correlation functions admits a simple hydrodynamic description in terms of the fermionic occupations in the initial state and the reflection and transmission amplitudes of the impurity. We derive this hydrodynamic framework from the ab initio calculation of the microscopic dynamics. This allows us to analytically characterize the effect of dissipation for several theoretically relevant initial states, such as a uniform Fermi sea, homogeneous product states, or the inhomogeneous state obtained by joining two Fermi seas. In this latter setting, when both gain and loss processes are present, we observe the emergence of exotic nonequilibrium steady states with stepwise uniform density profiles. In all instances, for strong dissipation the coherent dynamics of the system is arrested, which is a manifestation of the celebrated quantum Zeno effect.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Unbounded entanglement production via a dissipative impurity

We investigate the entanglement dynamics in a free-fermion chain initially prepared in a Fermi sea and subjected to localized losses (dissipative impurity). We derive a formula describing the dynamics of the entanglement entropies in the hydrodynamic limit of long times and large intervals. The result depends only on the absorption coefficient of the effective delta potential describing the impurity in the hydrodynamic limit. Genuine dissipation-induced entanglement is certified by the linear growth of the logarithmic negativity. Finally, in the quantum Zeno regime at strong dissipation the entanglement growth is arrested (Zeno entanglement death).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Closure of the entanglement gap at quantum criticality: The case of the Quantum Spherical Model

The study of entanglement spectra is a powerful tool to detect or elucidate universal behaviour in quantum many-body systems. We investigate the scaling of the entanglement (or Schmidt) gap $δξ$, i.e., the lowest laying gap of the entanglement spectrum, at a two-dimensional quantum critical point. We focus on the paradigmatic quantum spherical model, which exhibits a second-order transition, and is mappable to free bosons with an additional external constraint. We analytically show that the Schmidt gap vanishes at the critical point, although only logarithmically. For a system on a torus and the half-system bipartition, the entanglement gap vanishes as $π^2/\ln(L)$, with $L$ the linear system size. The entanglement gap is nonzero in the paramagnetic phase and exhibits a faster decay in the ordered phase. The rescaled gap $δξ\ln(L)$ exhibits a crossing for different system sizes at the transition, although logarithmic corrections prevent a precise verification of the finite-size scaling. Interestingly, the change of the entanglement gap across the phase diagram is reflected in the zero-mode eigenvector of the spin-spin correlator. At the transition quantum fluctuations give rise to a non-trivial structure of the eigenvector, whereas in the ordered phase it is flat. We also show that the vanishing of the entanglement gap at criticality can be qualitatively but not quantitatively captured by neglecting the structure of the zero-mode eigenvector.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Entanglement equipartition in critical random spin chains

The reduced density matrix of many-body systems possessing an additive conserved quantity can be decomposed in orthogonal sectors which can be independently analyzed. Recently, these have been proven to equally contribute to entanglement entropy for one dimensional conformal and integrable systems. In this paper, we extend this equipartition theorem to the disordered critical systems by studying the random singlet phase. We analytically compute the disorder averaged symmetry resolved Rényi entropies and show the leading orders are independent of the symmetry sector. Our findings are cross-checked with simulations within the numerical strong disorder renormalization group. We also identify the first subleading term breaking equipartition which is of the form $s^2/\ln\ell$ where $s$ is the magnetization of a subsystem of length $\ell$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Entanglement gap, corners, and symmetry breaking

We investigate the finite-size scaling of the lowest entanglement gap $δξ$ in the ordered phase of the two-dimensional quantum spherical model (QSM). The entanglement gap decays as $δξ=Ω/\sqrt{L\ln(L)}$. This is in contrast with the purely logarithmic behaviour as $δξ=π^2/\ln(L)$ at the critical point. The faster decay in the ordered phase reflects the presence of magnetic order. We analytically determine the constant $Ω$, which depends on the low-energy part of the model dispersion and on the geometry of the bipartition. In particular, we are able to compute the corner contribution to $Ω$, at least for the case of a square corner.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Entanglement revivals as a probe of scrambling in finite quantum systems

The entanglement evolution after a quantum quench became one of the tools to distinguish integrable versus chaotic (non-integrable) quantum many-body dynamics. Following this line of thoughts, here we propose that the revivals in the entanglement entropy provide a finite-size diagnostic benchmark for the purpose. Indeed, integrable models display periodic revivals manifested in a dip in the block entanglement entropy in a finite system. On the other hand, in chaotic systems, initial correlations get dispersed in the global degrees of freedom (information scrambling) and such a dip is suppressed. We show that while for integrable systems the height of the dip of the entanglement of an interval of fixed length decays as a power law with the total system size, upon breaking integrability a much faster decay is observed, signalling strong scrambling. Our results are checked by exact numerical techniques in free-fermion and free-boson theories, and by time-dependent density matrix renormalisation group in interacting integrable and chaotic models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Molecular dynamics simulation of entanglement spreading in generalized hydrodynamics

We consider a molecular dynamics method, the so-called flea gas for computing the evolution of entanglement after inhomogeneous quantum quenches in an integrable quantum system. In such systems the evolution of local observables is described at large space-time scales by the Generalized Hydrodynamics approach, which is based on the presence of stable, ballistically propagating quasiparticles. Recently it was shown that the GHD approach can be joined with the quasiparticle picture of entanglement evolution, providing results for entanglement growth after inhomogeneous quenches. Here we apply the flea gas simulation of GHD to obtain numerical results for entanglement growth. We implement the flea gas dynamics for the gapped anisotropic Heisenberg XXZ spin chain, considering quenches from globally homogeneous and piecewise homogeneous initial states. While the flea gas method applied to the XXZ chain is not exact even in the scaling limit (in contrast to the Lieb--Liniger model), it yields a very good approximation of analytical results for entanglement growth in the cases considered. Furthermore, we obtain the {\it full-time} dynamics of the mutual information after quenches from inhomogeneous settings, for which no analytical results are available.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Spreading of correlations in Markovian open quantum systems

Understanding the spreading of quantum correlations in out-of-equilibrium many-body systems is one of the major challenges in physics. For {\it isolated} systems, a hydrodynamic theory explains the origin and spreading of entanglement via the propagation of quasi-particle pairs. However, when systems interact with their surrounding much less has been established. Here we show that the quasi-particle picture remains valid for open quantum systems: while information is still spread by quasiparticles, the environment modifies their correlation and introduces incoherent and mixing effects. For free fermions with gain/loss dissipation we provide formulae fully describing incoherent and quasiparticle contributions in the spreading of entropy and mutual information. Importantly, the latter is not affected by incoherent correlations. The mutual information is exponentially damped at short times and eventually vanishes signalling the onset of a classical limit. The behaviour of the logarithmic negativity is similar and this scenario is common to other dissipations. For weak dissipation, the presence of quasiparticles underlies remarkable scaling behaviors.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Entanglement and classical fluctuations at finite-temperature critical points

We investigate several entanglement-related quantities at finite-temperature criticality in the three-dimensional quantum spherical model, both as a function of temperature $T$ and of the quantum parameter $g$, which measures the strength of quantum fluctuations. While the von Neumann and the Rényi entropies exhibit the volume-law for any $g$ and $T$, the mutual information obeys an area law. The prefactors of the volume-law and of the area-law are regular across the transition, reflecting that universal singular terms vanish at the transition. This implies that the mutual information is dominated by nonuniversal contributions. This hampers its use as a witness of criticality, at least in the spherical model. We also study the logarithmic negativity. For any value of $g,T$, the negativity exhibits an area-law. The negativity vanishes deep in the paramagnetic phase, it is larger at small temperature, and it decreases upon increasing the temperature. For any $g$, it exhibits the so-called sudden death, i.e., it is exactly zero for large enough $T$. The vanishing of the negativity defines a "death line", which we characterise analytically at small $g$. Importantly, for any finite $T$ the area-law prefactor is regular across the transition, whereas it develops a cusp-like singularity in the limit $T\to 0$. Finally, we consider the single-particle entanglement and negativity spectra. The vast majority of the levels are regular across the transition. Only the larger ones exhibit singularities. These are related to the presence of a zero mode, which reflects the symmetry breaking. This implies the presence of sub-leading singular terms in the entanglement entropies. Interestingly, since the larger levels do not contribute to the negativity, sub-leading singular corrections are expected to be suppressed for the negativity.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Negativity spectrum of one-dimensional conformal field theories

The partial transpose $ρ_A^{T_2}$ of the reduced density matrix $ρ_A$ is the key object to quantify the entanglement in mixed states, in particular through the presence of negative eigenvalues in its spectrum. Here we derive analytically the distribution of the eigenvalues of $ρ_A^{T_2}$, that we dub negativity spectrum, in the ground sate of gapless one-dimensional systems described by a Conformal Field Theory (CFT), focusing on the case of two adjacent intervals. We show that the negativity spectrum is universal and depends only on the central charge of the CFT, similarly to the entanglement spectrum. The precise form of the negativity spectrum depends on whether the two intervals are in a pure or mixed state, and in both cases, a dependence on the sign of the eigenvalues is found. This dependence is weak for bulk eigenvalues, whereas it is strong at the spectrum edges. We also investigate the scaling of the smallest (negative) and largest (positive) eigenvalues of $ρ_A^{T_2}$. We check our results against DMRG simulations for the critical Ising and Heisenberg chains, and against exact results for the harmonic chain, finding good agreement for the spectrum, but showing that the smallest eigenvalue is affected by very large scaling corrections.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The classical mutual information in mean-field spin glass models

We investigate the classical Renyi entropy S_n and the associated mutual information I_n in the Sherrington-Kirkpatrick (S-K) model, which is the paradigm model of mean-field spin glasses. Using classical Monte Carlo simulations and analytical tools we investigate the S-K model on the n-sheets booklet. This is obtained by gluing together n independent copies of the model, and it is the main ingredient to construct the Renyi entanglement-related quantities. We find a glassy phase at low temperature, whereas at high temperature the model exhibits paramagnetic behavior, consistent with the regular S-K model. The temperature of the paramagnetic-glassy transition depends non-trivially on the geometry of the booklet. At high-temperatures we provide the exact solution of the model by exploiting the replica symmetry. This is the permutation symmetry among the fictitious replicas that are used to perform disorder averages (via the replica trick). In the glassy phase the replica symmetry has to be broken. Using a generalization of the Parisi solution, we provide analytical results for S_n and I_n, and for standard thermodynamic quantities. Both S_n and I_n exhibit a volume law in the whole phase diagram. We characterize the behavior of the corresponding densities S_n/N, I_n/N, in the thermodynamic limit. Interestingly, at the critical point the mutual information does not exhibit any crossing for different system sizes, in contrast with local spin models.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The entanglement negativity in random spin chains

We investigate the logarithmic negativity in strongly-disordered spin chains in the random-singlet phase. We focus on the spin-1/2 random Heisenberg chain and the random XX chain. We find that for two arbitrary intervals the disorder-averaged negativity and the mutual information are proportional to the number of singlets shared between the two intervals. Using the strong-disorder renormalization group (SDRG), we prove that the negativity of two adjacent intervals grows logarithmically with the intervals length. In particular, the scaling behavior is the same as in conformal field theory, but with a different prefactor. For two disjoint intervals the negativity is given by a universal simple function of the cross ratio, reflecting scale invariance. As a function of the distance of the two intervals, the negativity decays algebraically in contrast with the exponential behavior in clean models. We confirm our predictions using a numerical implementation of the SDRG method. Finally, we also implement DMRG simulations for the negativity in open spin chains. The chains accessible in the presence of strong disorder are not sufficiently long to provide a reliable confirmation of the SDRG results.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The quench action approach in finite integrable spin chains

We consider the problem of constructing the stationary state following a quantum quench, using the exact overlaps for finite size integrable models. We focus on the isotropic Heisenberg spin chain with initial state Néel or Majumdar-Ghosh (dimer), although the proposed approach is valid for an arbitrary integrable model. We consider only eigenstates which do not contain zero-momentum strings because the latter are affected by fictitious singularities that are very difficult to take into account. We show that the fraction of eigenstates that do not contain zero-momentum strings is vanishing in the thermodynamic limit. Consequently, restricting to this part of the Hilbert space leads to vanishing expectation values of local observables. However, it is possible to reconstruct the asymptotic values by properly reweighting the expectations in the considered subspace, at the price of introducing finite-size corrections. We also develop a Monte Carlo sampling of the Hilbert space which allows us to study larger systems. We accurately reconstruct the expectation values of the conserved charged and the root distributions in the stationary state, which turn out to match the exact thermodynamic results. The proposed method can be implemented even in cases in which an analytic thermodynamic solution is not obtainable.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Eigenstate thermalization hypothesis (ETH) and integrability in quantum spin chains

We investigate the eigenstate thermalization hypothesis (ETH) in integrable models, focusing on the spin-1/2 isotropic Heisenberg (XXX) chain. We provide numerical evidence that ETH holds for typical eigenstates (weak ETH scenario). Specifically, using a numerical implementation of state-of-the-art Bethe ansatz results, we study the finite-size scaling of the eigenstate-to-eigenstate fluctuations of the reduced density matrix. We find that fluctuations are normally distributed, and their standard deviation decays in the thermodynamic limit as L^{-1/2}, with L the size of the chain. This is in contrast with the exponential decay that is found in generic non-integrable systems. Based on our results, it is natural to expect that this scenario holds in other integrable spin models and for typical local observables. Finally, we investigate the entanglement properties of the excited states of the XXX chain. We numerically verify that typical mid-spectrum eigenstates exhibit extensive entanglement entropy (i.e., volume-law scaling).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Overlap distributions for quantum quenches in the anisotropic Heisenberg chain

The dynamics after a quantum quench is determined by the weights of the initial state in the eigenspectrum of the final Hamiltonian, i.e., by the distribution of overlaps in the energy spectrum. We present an analysis of such overlap distributions for quenches of the anisotropy parameter in the one-dimensional anisotropic spin-1/2 Heisenberg model (XXZ chain). We provide an overview of the form of the overlap distribution for quenches from various initial anisotropies to various final ones, using numerical exact diagonalization. We show that if the system is prepared in the antiferromagnetic Néel state (infinite anisotropy) and released into a non-interacting setup (zero anisotropy, XX point) only a small fraction of the final eigenstates gives contributions to the post-quench dynamics, and that these eigenstates have identical overlap magnitudes. We derive expressions for the overlaps, and present the selection rules that determine the final eigenstates having nonzero overlap. We use these results to derive concise expressions for time-dependent quantities (Loschmidt echo, longitudinal and transverse correlators) after the quench. We use perturbative analyses to understand the overlap distribution for quenches from infinite to small nonzero anisotropies, and for quenches from large to zero anisotropy.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Simulating the Generalized Gibbs Ensemble (GGE): a Hilbert space Monte Carlo approach

By combining classical Monte Carlo and Bethe ansatz techniques we devise a numerical method to construct the Truncated Generalized Gibbs Ensemble (TGGE) for the spin-1/2 isotropic Heisenberg ($XXX$) chain. The key idea is to sample the Hilbert space of the model with the appropriate GGE probability measure. The method can be extended to other integrable systems, such as the Lieb-Liniger model. We benchmark the approach focusing on GGE expectation values of several local observables. As finite-size effects decay exponentially with system size, moderately large chains are sufficient to extract thermodynamic quantities. The Monte Carlo results are in agreement with both the Thermodynamic Bethe Ansatz (TBA) and the Quantum Transfer Matrix approach (QTM). Remarkably, it is possible to extract in a simple way the steady-state Bethe-Gaudin-Takahashi (BGT) roots distributions, which encode complete information about the GGE expectation values in the thermodynamic limit. Finally, it is straightforward to simulate extensions of the GGE, in which, besides the local integral of motion (local charges), one includes arbitrary functions of the BGT roots. As an example, we include in the GGE the first non-trivial quasi-local integral of motion.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Bethe ansatz description of edge-localization in the open-boundary XXZ spin chain

At large values of the anisotropy Δ, the open-boundary Heisenberg spin-1/2 chain has eigenstates displaying localization at the edges. We present a Bethe ansatz description of this `edge-locking' phenomenon in the entire Δ>1 region. We focus on the simplest spin sectors, namely the highly polarized sectors with only one or two overturned spins, i.e., one-particle and two-particle sectors. Edge-locking is associated with pure imaginary solutions of the Bethe equations, which are not commonly encountered in periodic chains. In the one-particle case, at large anisotropies there are two eigenstates with imaginary Bethe momenta, related to localization at the two edges. For any finite chain size, one of the two solutions become real as the anisotropy is lowered below a certain value. For two particles, a richer scenario is observed, with eigenstates having the possibility of both particles locked on the same or different edge, one locked and the other free, and both free either as single magnons or as bound composites corresponding to `string' solutions. For finite chains, some of the edge-locked spins get delocalized at certain values of the anisotropy (`exceptional points'), corresponding to imaginary solutions becoming real. We characterize these phenomena thoroughly by providing analytic expansions of the Bethe momenta for large chains, large anisotropy, and near the exceptional points. In the large-chain limit all the exceptional points coalesce at the isotropic point (Δ=1) and edge-locking becomes stable in the whole Δ>1 region.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Entanglement negativity and conformal field theory: a Monte Carlo study

We investigate the behavior of the moments of the partially transposed reduced density matrix ρ^{T_2}_A in critical quantum spin chains. Given subsystem A as union of two blocks, this is the (matrix) transposed of ρ_A with respect to the degrees of freedom of one of the two. This is also the main ingredient for constructing the logarithmic negativity. We provide a new numerical scheme for calculating efficiently all the moments of ρ_A^{T_2} using classical Monte Carlo simulations. In particular we study several combinations of the moments which are scale invariant at a critical point. Their behavior is fully characterized in both the critical Ising and the anisotropic Heisenberg XXZ chains. For two adjacent blocks we find, in both models, full agreement with recent CFT calculations. For disjoint ones, in the Ising chain finite size corrections are non negligible. We demonstrate that their exponent is the same governing the unusual scaling corrections of the mutual information between the two blocks. Monte Carlo data fully match the theoretical CFT prediction only in the asymptotic limit of infinite intervals. Oppositely, in the Heisenberg chain scaling corrections are smaller and, already at finite (moderately large) block sizes, Monte Carlo data are in excellent agreement with the asymptotic CFT result.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Entanglement negativity via replica trick: a Quantum Monte Carlo approach

Motivated by recent developments in conformal field theory (CFT), we devise a Quantum Monte Carlo (QMC) method to calculate the moments of the partially transposed reduced density matrix at finite temperature. These are used to construct scale invariant combinations that are related to the negativity, a true measure of entanglement for two intervals embedded in a chain. These quantities can serve as witnesses of criticality. In particular, we study several scale invariant combinations of the moments for the 1D hard-core boson model. For two adjacent intervals unusual finite size corrections are present, showing parity effects that oscillate with a filling dependent period. These are more pronounced in the presence of boundaries. For large chains we find perfect agreement with CFT. Oppositely, for disjoint intervals corrections are more severe and CFT is recovered only asymptotically. Furthermore, we provide evidence that their exponent is the same as that governing the corrections of the mutual information. Additionally we study the 1D Bose-Hubbard model in the superfluid phase. Remarkably, the finite-size effects are smaller and QMC data are already in impressive agreement with CFT at moderate large sizes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Entanglement spectrum of the two dimensional Bose-Hubbard model

We study the entanglement spectrum (ES) of the Bose-Hubbard model on the two dimensional square lattice at unit filling, both in the Mott insulating and in the superfluid phase. In the Mott phase, we demonstrate that the ES is dominated by the physics at the boundary between the two subsystems. On top of the boundary-local (perturbative) structure, the ES exhibits substructures arising from one-dimensional dispersions along the boundary. In the superfluid phase, the structure of the ES is qualitatively different, and reflects the spontaneously broken U(1) symmetry of the phase. We attribute the basic low-lying structure to a so-called "tower of states" (TOS) Hamiltonian of the model. We then discuss how these characteristic structures evolve across the superfluid to Mott insulator transition and their influence on the behavior of the entanglement entropies. Finally, we briefly outline the implications of the ES structure on the efficiency of matrix-product-state based algorithms in two dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Boundary-locality and perturbative structure of entanglement spectra in gapped systems

The entanglement between two parts of a many-body system can be characterized in detail by the entanglement spectrum. Focusing on gapped phases of one-dimensional systems, we show how this spectrum is dominated by contributions from the boundary between the parts. This contradicts the view of an "entanglement Hamiltonian" as a bulk entity. The boundary-local nature of the entanglement spectrum is clarified through its hierarchical level structure, through the combination of two single-boundary spectra to form a two-boundary spectrum, and finally through consideration of dominant eigenfunctions of the entanglement Hamiltonian. We use the boundary-locality to formulate a perturbative scheme for calculating entanglement spectra.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Entanglement spectrum of the Heisenberg XXZ chain near the ferromagnetic point

We study the entanglement spectrum (ES) of a finite XXZ spin 1/2 chain in the limit Δ-> -1^+ for both open and periodic boundary conditions. At Δ=-1 (ferromagnetic point) the model is equivalent to the Heisenberg ferromagnet and its degenerate ground state manifold is the SU(2) multiplet with maximal total spin. Any state in this so-called "symmetric sector" is an equal weight superposition of all possible spin configurations. In the gapless phase at Δ>-1 this property is progressively lost as one moves away from the Δ=-1 point. We investigate how the ES obtained from the states in this manifold reflects this change, using exact diagonalization and Bethe ansatz calculations. We find that in the limit Δ->-1^+ most of the ES levels show divergent behavior. Moreover, while at Δ=-1 the ES contains no information about the boundaries, for Δ>-1 it depends dramatically on the choice of boundary conditions. For both open and periodic boundary conditions the ES exhibits an elegant multiplicity structure for which we conjecture a combinatorial formula. We also study the entanglement eigenfunctions, i.e. the eigenfunctions of the reduced density matrix. We find that the eigenfunctions corresponding to the non diverging levels mimic the behavior of the state wavefunction, whereas the others show intriguing polynomial structures. Finally we analyze the distribution of the ES levels as the system is detuned away from Δ=-1.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Entanglement entropy of two disjoint intervals in c=1 theories

We study the scaling of the Renyi entanglement entropy of two disjoint blocks of critical lattice models described by conformal field theories with central charge c=1. We provide the analytic conformal field theory result for the second order Renyi entropy for a free boson compactified on an orbifold describing the scaling limit of the Ashkin-Teller (AT) model on the self-dual line. We have checked this prediction in cluster Monte Carlo simulations of the classical two dimensional AT model. We have also performed extensive numerical simulations of the anisotropic Heisenberg quantum spin-chain with tree-tensor network techniques that allowed to obtain the reduced density matrices of disjoint blocks of the spin-chain and to check the correctness of the predictions for Renyi and entanglement entropies from conformal field theory. In order to match these predictions, we have extrapolated the numerical results by properly taking into account the corrections induced by the finite length of the blocks to the leading scaling behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Entanglement entropy of two disjoint blocks in critical Ising models

We study the scaling of the Renyi and entanglement entropy of two disjoint blocks of critical Ising models, as function of their sizes and separations. We present analytic results based on conformal field theory that are quantitatively checked in numerical simulations of both the quantum spin chain and the classical two dimensional Ising model. Theoretical results match the ones obtained from numerical simulations only after taking properly into account the corrections induced by the finite length of the blocks to their leading scaling behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Magnetic and glassy transitions in the square-lattice XY model with random phase shifts

We investigate the magnetic and glassy transitions of the square-lattice XY model in the presence of random phase shifts. We consider two different random-shift distributions: the Gaussian distribution and a slightly different distribution (cosine distribution) which allows the exact determination of the Nishimori line where magnetic and overlap correlation functions are equal. We perform Monte Carlo simulations for several values of the temperature and of the variance of the disorder distribution, in the paramagnetic phase close to the magnetic and glassy transition lines. We find that, along the transition line separating the paramagnetic and the quasi-long-range order phases, magnetic correlation functions show a universal Kosterlitz-Thouless behavior as in the pure XY model, while overlap correlations show a disorder-dependent critical behavior. This behavior is observed up to a multicritical point which, in the cosine model, lies on the Nishimori line. Finally, for large values of the disorder variance, we observe a universal zero-temperature glassy critical transition, which is in the same universality class as that occurring in the gauge-glass model.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The three-dimensional gauge-glass model

We investigate the temperature-disorder (T-S) phase diagram of a three-dimensional gauge glass model, which is a cubic-lattice nearest-neighbor XY model with quenched random phase shifts A_xy at the bonds, by numerical Monte Carlo simulations. We consider the uncorrelated phase-shift distribution P(A_xy)\sim \exp[(cos A_xy)/S], which has the pure XY model and the uniform distribution of random shifts as extreme cases at S=0 and S->infty respectively, and which gives rise to equal magnetic and overlap correlation functions when T=S. While the high-temperature phase is always paramagnetic, at low temperatures there is a ferromagnetic phase for weak disorder (small S) and a glassy phase at large disorder (large S). These three phases are separated by transition lines with different magnetic and glassy critical behaviors. The disorder induced by the random shifts turns out to be irrelevant at the paramagnetic-ferromagnetic transition line, where the critical behavior belongs to the 3D XY universality class of pure systems; disorder gives only rise to very slowly decaying scaling corrections. The glassy critical behavior along the finite-temperature paramagnetic-glassy transition line belongs to the gauge-glass universality class, with a quite large critical exponent nu=3.2(4). These transition lines meet at a multicritical point M, located at T=S=0.7840(2). The low-temperature ferromagnetic and glassy phases are separated by a third transition line, from M down to the T=0 axis, which is slightly reentrant.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Quasi-long-range order in the 2D XY model with random phase shifts

We study the square-lattice XY model in the presence of random phase shifts. We consider two different disorder distributions with zero average shift and investigate the low-temperature quasi-long-range order phase which occurs for sufficiently low disorder. By means of Monte Carlo simulations we determine several universal quantities which are then compared with the analytic predictions of the random spin-wave theory. We observe a very good agreement which indicates that the universal long-distance behavior in the whole low-disorder low-temperature phase is fully described by the random spin-wave theory.

0 citations0 reviewsNo code linkedOpen access