We study the effects of integrability-breaking perturbations on the nonequilibrium evolution of many-particle quantum systems. We focus on a class of spinless fermion models with weak interactions. We employ equation of motion techniques that can be viewed as generalizations of quantum Boltzmann equations. We benchmark our method against time-dependent density matrix renormalization group computations and find it to be very accurate as long as interactions are weak. For small integrability breaking, we observe robust prethermalization plateaux for local observables on all accessible time scales. Increasing the strength of the integrability-breaking term induces a "drift" away from the prethermalization plateaux towards thermal behavior. We identify a time scale characterizing this crossover.

We consider quantum quenches in an integrable quantum chain with tuneable-integrabilitybreaking interactions. In the case where these interactions are weak, we demonstrate that at intermediate times after the quench local observables relax to a prethermalized regime, which can be described by a density matrix that can be viewed as a deformation of a generalized Gibbs ensemble. We present explicit expressions for the approximately conserved charges characterizing this ensemble. We do not find evidence for a crossover from the prethermalized to a thermalized regime on the time scales accessible to us. Increasing the integrability-breaking interactions leads to a behavior that is compatible with eventual thermalization.

We present exact analytical solutions for the zero-energy modes of two-dimensional massless Dirac fermions fully confined within a smooth one-dimensional potential V͑x͒ =−␣ / cosh͑␤x͒, which provides a good fit for potential profiles of existing top-gated graphene structures. We show that there is a threshold value of the characteristic potential strength ␣ / ␤ for which the first mode appears, in striking contrast to the nonrelativistic case. A simple relationship between the characteristic strength and the number of modes within the potential is found. An experimental setup is proposed for the observation of these modes. The proposed geometry could be utilized in future graphene-based devices with high on/off current ratios.

We review two important non-perturbative approaches for extracting the physics of low- dimensional strongly correlated quantum systems. Firstly, we start by providing a comprehensive review of non-Abelian bosonization. This includes an introduction to the basic elements of confor- mal field theory as applied to systems with a current algebra, and we orient the reader by presenting a number of applications of non-Abelian bosonization to models with large symmetries. We then tie this technique into recent advances in the ability of cold atomic systems to realize complex symme- tries. Secondly, we discuss truncated spectrum methods for the numerical study of systems in one and two dimensions. For one-dimensional systems we provide the reader with considerable insight into the methodology by reviewing canonical applications of the technique to the Ising model (and its variants) and the sine-Gordon model. Following this we review recent work on the development of renormalization groups, both numerical and analytical, that alleviate the effects of truncating the spectrum. Using these technologies, we consider a number of applications to one-dimensional systems: properties of carbon nanotubes, quenches in the Lieb-Liniger model, 1+1D quantum chro- modynamics, as well as Landau-Ginzburg theories. In the final part we move our attention to consider truncated spectrum methods applied to two-dimensional systems. This involves combining truncated spectrum methods with matrix product state algorithms. We describe applications of this method to two-dimensional systems of free fermions and the quantum Ising model, including their non-equilibrium dynamics.

We employ equation of motion techniques to study the non-equilibrium dynamics in a lattice model of weakly interacting spinless fermions. Our model provides a simple setting for analyzing the effects of weak integrability breaking perturbations on the time evolution after a quantum quench. We establish the accuracy of the method by comparing results at short and intermediate times to timedependent density matrix renormalization group computations. For sufficiently weak integrabilitybreaking interactions we always observe prethermalization plateaux, where local observables relax to non-thermal values at intermediate time scales. At later times a crossover towards thermal behaviour sets in. We determine the associated time scale, which depends on the initial state, the band structure of the non-interacting theory, and the strength of the integrability breaking perturbation. Our method allows us to analyze in some detail the spreading of correlations and in particular the structure of the associated light cones in our model. We find that the interior and exterior of the light cone are separated by an intermediate region, the temporal width of which appears to scale with a universal power-law t 1/3 .

We show that confinement in the quantum Ising model leads to nonthermal eigenstates, in both continuum and lattice theories, in both one (1D) and two dimensions (2D). In the ordered phase, the presence of a confining longitudinal field leads to a profound restructuring of the excitation spectrum, with the low-energy two-particle continuum being replaced by discrete 'meson' modes (linearly confined pairs of domain walls). These modes exist far into the spectrum and are atypical, in the sense that expectation values in the state with energy E do not agree with the microcanonical (thermal) ensemble prediction. Single meson states persist above the two meson threshold, due to a surprising lack of hybridization with the (n ≥ 4)-domain wall continuum, a result that survives into the thermodynamic limit and that can be understood from analytical calculations. The presence of such states is revealed in anomalous post-quench dynamics, such as the lack of a light cone, the suppression of the growth of entanglement entropy, and the absence of thermalization for some initial states. The nonthermal states are confined to the ordered phase -the disordered (paramagnetic) phase exhibits typical thermalization patterns in both 1D and 2D in the absence of integrability.

There is a dichotomy in the nonequilibrium dynamics of quantum many body systems. In the presence of integrability, expectation values of local operators equilibrate to values described by a generalized Gibbs ensemble, which retains extensive memory about the initial state of the system. On the other hand, in generic systems such expectation values relax to stationary values described by the thermal ensemble, fixed solely by the energy of the state. At the heart of understanding this dichotomy is the eigenstate thermalization hypothesis (ETH): individual eigenstates in nonintegrable systems are thermal, in the sense that expectation values agree with the thermal prediction at a temperature set by the energy of the eigenstate. In systems where ETH is violated, thermalization can be avoided. Thus establishing the range of validity of ETH is crucial in understanding whether a given quantum system thermalizes. Here we study a simple model with confinement, the quantum Ising chain with a longitudinal field, in which ETH is violated. Despite an absence of integrability, there exist rare (nonthermal) states that persist far into the spectrum. These arise as a direct consequence of confinement: pairs of particles are confined, forming new 'meson' excitations whose energy can be extensive in the system size. We show that such states are nonthermal in both the continuum and in the low-energy spectrum of the corresponding lattice model. We highlight that the presence of such states within the spectrum has important consequences, with certain quenches leading to an absence of thermalization and local observables evolving anomalously. arXiv:1808.10782v2 [cond-mat.str-el]

Over the past two decades, advances in computational algorithms have revealed a curious property of the two-dimensional Hubbard model (and related theories) with hole doping: the presence of close-in-energy competing ground states that display very different physical properties. On the one hand, there is a complicated state exhibiting intertwined spin, charge, and pair density wave orders. We call this 'type A'. On the other hand, there is a uniform d-wave superconducting state that we denote as 'type B'. We advocate, with the support of both microscopic theoretical calculations and experimental data, dividing the high-temperature cuprate superconductors into two corresponding families, whose properties reflect either the type A or type B ground states at low temperatures. We review the anomalous properties of the pseudogap phase that led us to this picture, and present a modern perspective on the role that umklapp scattering plays in these phenomena in the type B materials. This reflects a consistent framework that has emerged over the last decade, in which Mott correlations at weak coupling drive the formation of the pseudogap. We discuss this development, recent theory and experiments, and open issues. * If one did not refermionise the spin sector, the non-linear interaction term within the bosonic theory reads −4g a =b cos Θ a cos Θ b , which is suggestive of a high symmetry. * This is easy to see in the purely bosonic language, as the interaction term −4g cos Θ sf cos Θ f at strong coupling pins the bosons to (Θ sf , Θ f ) = (2mπ, 2nπ or (Θ sf , Θ f ) = [2m + 1]π, [2n + 1]π with n, m ∈ Z. * Interestingly, recent gauge theory calculations that are proposed to apply to the pseudogap phase of the cuprates find such a deformation of the Fermi surface. See figure 7(a) of (Sachdev et al. 2018).