We investigate the distribution of the eigenvalues of the reduced density matrix (entanglement spectrum) after a global quantum quench. We show that in an appropriate scaling limit the lower part of the entanglement spectrum exhibits ``universality''. In the scaling limit and at asymptotically long times the distribution of the entanglement spectrum depends on two parameters that can be determined from the Rényi entropies. We show that two typical scenarios occur. In the first one, the distribution of the entanglement spectrum levels is similar to the one describing the ground-state entanglement spectrum in Conformal Field Theories. In the second scenario, the lower levels of the entanglement spectrum are highly degenerate and their distribution is given by a series of Dirac deltas. We benchmark our analytical results in free-fermion chains, such as the transverse field Ising chain and the XX chain, in the rule 54 chain, and in Bethe ansatz solvable spin models.
We discuss the entanglement spectrum of the ground state of a (1+1)-dimensional system in a gapped phase near a quantum phase transition. In particular, in proximity to a quantum phase transition described by a conformal field theory (CFT), the system is represented by a gapped Lorentz invariant field theory in the ``scaling limit'' (correlation length $\ensuremath{\xi}$ much larger than microscopic ``lattice'' scale ``
a''), and can be thought of as a CFT perturbed by a relevant perturbation. We show that for such (1+1) gapped Lorentz invariant field theories in infinite space, the low-lying entanglement spectrum obtained by tracing out, say, left half-infinite space, is precisely equal to the physical spectrum of the unperturbed gapless, i.e., conformal field theory defined on a finite interval of length ${L}_{\ensuremath{\xi}}=\phantom{\rule{0.16em}{0ex}}ln(\ensuremath{\xi}/a)$ with certain boundary conditions. In particular, the low-lying entanglement spectrum of the gapped theory is the finite-size spectrum of a boundary conformal field theory, and is always discrete and universal. Each relevant perturbation, and thus each gapped phase in proximity to the quantum phase transition, maps into a particular boundary condition. A similar property has been known to hold for Baxter's corner transfer matrices in a very special class of fine-tuned, namely, integrable off-critical lattice models, for the entire entanglement spectrum and independent of the scaling limit. In contrast, our result applies to completely general gapped Lorentz invariant theories in the scaling limit, without the requirement of integrability, for the low-lying entanglement spectrum. While the entanglement spectrum of the ground state of a gapped theory on a finite interval of length
2R with suitable boundary conditions, bipartitioned into two equal pieces, turns out to exhibit a crossover between the finite-size spectra of the same CFT with in general different boundary conditions as the system size
R crosses the correlation length from the ``critical regime'' $R\ensuremath{\ll}\ensuremath{\xi}$ to the ``gapped regime'' $R\ensuremath{\gg}\ensuremath{\xi}$, the physical spectrum on a finite interval of length
R with the same boundary conditions, on the other hand, is known to undergo a dramatic reorganization during the same crossover from being discrete to being continuous.
We study the scaling behavior of the eigenvalues of correlation matrices, which characterize the entanglement of a subsystem with its complement part of a total pure state. A most distinguishing feature of entanglement entropy is its logarithmic dependence on the subsystem size for the groundstate of one-dimensional critical systems. Despite its robust universal character and relevance to a wide range of topics, a thorough understanding of this result requires sophisticated mathematical physics techniques or conformal field theory. The aim of our work is to shed light on this from the underlying eigenvalue distribution perspective. The central object is the correlation matrix, which takes the form of Toeplitz or block-Toeplitz matrix. We develop a circulant matrix approximation in the large matrix dimension limit, thus allowing for the individual eigenvalues behavior to be analysed analytically. We find that for both free lattice fermions and transverse field Ising chain, eigenvalues in the bulk of the eigenspectrum scales as
1/LA with the subsystem size. Together with the extensivity of the entropy function, it explains the robust
log2LA scaling of entanglement at criticality. Perturbing from the entanglement-free limit of the Ising chain, we find a smooth crossover behavior to `non-critical' scaling that is characterized by a very slow logarithmic dependence rendering it seemingly a constant entanglement value expected of non-critical systems.
The following thesis focuses on the scaling of entanglement entropy in lower dimensions and is divided into three main parts. Chapter 2 studies the thermal reduced density matrices in fermion and spin systems on ladders. Chapter 3 studies the many-body localization phase transition in a Rokhsar-Kivelson type wave function. Chapter 4 studies the subleading correction term of entanglement entropy in 2+1 dimensional scale invariant systems. Chapter 5 studies the bulk-boundary correspondence in 3 + 1 dimensional topological phases and its entanglement entropy. In chapter 2, we investigate the reduced density matrices for a model of free fermions on a two-leg ladder (gapped by the inter-chain tunneling operator) and in 1/2 spin systems on a ladder with a gapped ground state using exact solutions and several controlled approximations. We calculate the reduced density matrix and the entanglement entropy for a leg of the ladder (i.e. cut made between the chains). In the fermionic system we find the exact form of the reduced density matrix for one of the chains and determine the entanglement spectrum explicitly. Here we find that in the weak tunneling limit of the ladder the entanglement entropy of one chain of the gapped ladder has a simple and universal form dictated by conformal invariance. In the case of the spin system, we consider the strong coupling limit by using perturbation theory and get the reduced density matrix by the Schmidt decomposition. The entanglement entropies of a general gapped system of two coupled conformal field theories (in 1+1 dimensions) is discussed using the replica trick and scaling arguments. We show that 1) for a system with a bulk gap, the reduced density matrix has the form of a thermal density matrix, 2) the long-wavelength modes of one subsystem (a chain) of a gapped coupled system are always thermal, 3) the von Neumann entropy equals to the thermodynamic entropy of one chain, and 4) the bulk gap plays the role of effective temperature. In chapter 3, we construct a family of many-body wave functions to study the many-body localization phase transition. The wave functions have a Rokhsar-Kivelson form, in which the weight for the configurations are chosen from the Gibbs weights of a classical spin glass model, known as the Random Energy Model, multiplied by a random sign structure to represent a highly excited state. These wave functions show a phase transition into an MBL phase. In addition, we see three regimes of entanglement scaling with subsystem size: scaling with entanglement corresponding to an infinite temperature thermal phase, constant scaling, and a sub-extensive scaling between these limits. Near the phase transition point, the fluctuations of the Renyi entropies are non-Gaussian. We find that Renyi entropies with different Renyi index transition into the MBL phase at different points and have different scaling behavior, suggesting a multifractal behavior. In chapter 4, we study the universal scaling behavior of the entanglement entropy of critical theories in 2+1 dimensions. We specially consider two fermionic scale-invariant models, free massless Dirac fermions and a model of fermions with quadratic band touching, and numerically study the two-cylinder entanglement entropy of the models on the torus. We find that in both cases the entanglement entropy satisfies the area law and has the subleading term which is a scaling function of the aspect ratios of the cylindrical regions. We test the scaling of entanglement in both the free fermion models using three possible scaling functions for the subleading term derived from a) the quasi-one-dimensional conformal field theory, b) the bosonic quantum Lifshitz model, and c) the holographic AdS/CFT correspondence. For the later case we construct an analytic scaling function using holography, appropriate for critical theories with a gravitational dual description. We find that the subleading term in the fermionic models is well described, for a range of aspect ratios, by the scaling form derived from the quantum Lifshitz model as well as that derived using the AdS/CFT correspondence (in this case only for the Dirac model). For the case where the fermionic models are placed on a square torus we find the fit to the different scaling forms is in agreement to surprisingly high precision. In chapter 5, we discuss (2+1)-dimensional gapless surface theories of bulk (3+1)-dimensional topological phases, such as the BF theory at level K, and its generalization. In particular, we put these theories on a flat (2+1) dimensional torus T3 parameterized by its modular parameters, and compute the partition functions obeying various twisted boundary conditions. We show the partition functions are transformed into each other under SL(3, Z) modular transformations, and furthermore establish the bulk-boundary correspondence in (3+1) dimensions by matching the modular S and T matrices computed from the boundary field theories with those computed in the bulk. We propose the three-loop braiding statistics can be studied by constructing the modular S and T matrices from an appropriate boundary field theory. We also study the EE for 3 + 1- dimensional topological phase with or without three-loop braiding phase.
We study the entanglement entropies in one-dimensional open critical systems, whose effective description is given by a conformal field theory with boundaries. We show that, for pure-state systems formed by the ground state or by the excited states associated to primary fields, the entanglement entropies have a finite-size behavior that depends on the correlation of the underlying field theory. The analytical results are checked numerically, finding excellent agreement for the quantum chains ruled by the theories with central charge
c=1/2 and 1.
Entanglement spectrum of finite-size correlated electron systems are investigated using the Gutzwiller projection technique. The product of largest eigenvalue and rank of the block reduced density matrix, which is a measure of distance of the state from the maximally entangled state of the corresponding rank, is seen to characterise the insulator to metal crossover in the state. The fraction of distinct eigenvalues exhibits a `chaotic' behaviour in the crossover region, and it shows a `integrable' behaviour at both insulating and metallic ends. The integrated entanglement spectrum obeys conformal field theory (CFT) prediction at the metal and insulator ends, but shows a noticeable deviation from CFT prediction in the crossover regime, thus it can also track a metal-insulator crossover. A modification of the CFT result for the entanglement spectrum for finite size is proposed which holds in the crossover regime also. The adjacent level spacing distribution of unfolded non-zero eigenvalues for intermediate values of Gutzwiller projection parameter
g is the same as that of an ensemble of random matrices obtained by replacing each block of reduced density matrix by a random real symmetric Toeplitz matrix. It is strongly peaked at zero, with an exponential tail proportional to
e−(n/R)s, where
s is the adjacent level spacing,
n is number of distinct eigenvalues and
R is the rank of the reduced density matrix.