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 .