A key property of many-body localization, the localization of quantum particles in systems with both quenched disorder and interactions, is the area law entanglement of even highly excited eigenstates of many-body localized Hamiltonians. Matrix Product States (MPS) can be used to efficiently represent low entanglement (area law) wave functions in one dimension. An important application of MPS is the widely used Density Matrix Renormalization Group (DMRG) algorithm for finding ground states of one dimensional Hamiltonians. Here, we describe two algorithms, the Shift and Invert MPS (SIMPS) and excited state DMRG which finds highly-excited eigenstates of many-body localized Hamiltonians. Excited state DMRG uses a modified sweeping procedure to identify eigenstates whereas SIMPS is a shift-inverse procedure that applies the inverse of the shifted Hamiltonian to a MPS multiple times to project out the targeted eigenstate. To demonstrate the power of these methods we verify the breakdown of the Eigenstate Thermalization Hypothesis (ETH) in the manybody localized phase of the random field Heisenberg model, show the saturation of entanglement in the MBL phase and generate local excitations.

We introduce the cut averaged entanglement entropy in disordered periodic spin chains and prove it to be a concave function of subsystem size for individual eigenstates. This allows us to identify the entanglement scaling as a function of subsystem size for individual states in inhomogeneous systems. Using this quantity, we probe the critical region between the many-body localized (MBL) and ergodic phases in finite systems. In the middle of the spectrum, we show evidence for bimodality of the entanglement distribution in the MBL critical region, finding both volume law and area law eigenstates over disorder realizations as well as within single disorder realizations. The disorder averaged entanglement entropy in this region then scales as a volume law with a coefficient below its thermal value. We discover in the critical region, as we approach the thermodynamic limit, that the cut averaged entanglement entropy density falls on a one-parameter family of curves. Finally, we also show that without averaging over cuts the slope of the entanglement entropy vs. subsystem size can be negative at intermediate and strong disorder, caused by rare localized regions in the system.

Strong disorder in interacting quantum systems can give rise to the phenomenon of Many-Body Localization (MBL), which defies thermalization due to the formation of an extensive number of quasi local integrals of motion. The one particle operator content of these integrals of motion is related to the one particle orbitals of the one particle density matrix and shows a strong signature across the MBL transition as recently pointed out by Bera et al. [Phys. Rev. Lett. 115, 046603 (2015); Ann. Phys. 529, 1600356 (2017)]. We study the properties of the one particle orbitals of many-body eigenstates of an MBL system in one dimension. Using shift-and-invert MPS (SIMPS), a matrix product state method to target highly excited many-body eigenstates introduced in [Phys. Rev. Lett. 118, 017201 (2017)], we are able to obtain accurate results for large systems of sizes up to L = 64. We find that the one particle orbitals drawn from eigenstates at different energy densities have high overlap and their occupations are correlated with the energy of the eigenstates. Moreover, the standard deviation of the inverse participation ratio of these orbitals is maximal at the nose of the mobility edge. Also, the one particle orbitals decay exponentially in real space, with a correlation length that increases at low disorder. In addition, we find a "1/f " distribution of the coupling constants of a certain range of the number operators of the OPOs, which is related to their exponential decay.

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 Rényi entropies are non-Gaussian. We find that Rényi entropies with different Rényi index transition into the MBL phase at different points and have different scaling behavior, suggesting a multifractal behavior.

The twofold twist defects in the D(Z k ) quantum double model (abelian topological phase) carry non-abelian fractional Majorana-like characteristics. We align these twist defects in a line and construct a one dimensional Hamiltonian which only includes the pairwise interaction. For the defect chain with even number of twist defects, it is equivalent to the Z k clock model with periodic boundary condition (up to some phase factor for boundary term), while for odd number case, it maps to Z k clock model with duality twisted boundary condition. At critical point, for both cases, the twist defect chain enjoys an additional translation symmetry, which corresponds to the Kramers-Wannier duality symmetry in the Z k clock model and can be generated by a series of braiding operators for twist defects. We further numerically investigate the low energy excitation spectrum for k = 3, 4, 5 and 6. For even-defect chain, the critical points are the same as the Z k clock conformal field theories (CFTs), while for odd-defect chain, when k = 4, the critical points correspond to orbifolding a Z2 symmetry of CFTs of the even-defect chain. For k = 4 case, we numerically observe some similarity to the Z4 twist fields in SU (2)1/D4 orbifold CFT.