In this work, we study the non-equilibrium dynamics of low-dimensional strongly correlated quantum systems, which in our case are generated by sudden and time-dependent quantum quenches. Such non-equilibrium dynamics are a particularly difficult problem, because the time evolution of the initial many-body quantum state, which is governed by the post-quench Hamiltonian, is highly non-trivial. On top of that, applying an operator to this many-body quantum state in order to calculate its expectation values, typically creates an intricate set of excitations. Throughout the thesis, we furthermore contrast the results that we have obtained using exact numerical methods for two microscopic models, namely the Heisenberg and Fermi-Hubbard chains, with results derived from the Luttinger liquid theory, which, in equilibrium, describes the universal behavior of the two models.In the second part of the thesis, we reduce the dimensionality even further and employ complementary numerical and analytical tools to examine the unusual non-equilibrium behavior of a ring-shaped quantum impurity hosting interacting spin-less fermions.
In this work, we study the non-equilibrium dynamics of low-dimensional strongly correlated quantum systems, which in our case are generated by sudden and time-dependent quantum quenches. Such non-equilibrium dynamics are a particularly difficult problem, because the time evolution of the initial many-body quantum state, which is governed by the post-quench Hamiltonian, is highly non-trivial. On top of that, applying an operator to this many-body quantum state in order to calculate its expectation values, typically creates an intricate set of excitations. Throughout the thesis, we furthermore contrast the results that we have obtained using exact numerical methods for two microscopic models, namely the Heisenberg and Fermi-Hubbard chains, with results derived from the Luttinger liquid theory, which, in equilibrium, describes the universal behavior of the two models.In the second part of the thesis, we reduce the dimensionality even further and employ complementary numerical and analytical tools to examine the unusual non-equilibrium behavior of a ring-shaped quantum impurity hosting interacting spin-less fermions.
In this thesis, strongly correlated quantum many-body systems in equilibrium and in out-of-equilibrium situations are investigated. This is done by applying and developing well established numerical methods. The focus of the thesis lies in the development and application of the density matrix renormalization group method (DMRG) to quantum many-body systems out of equilbrium. In this thesis, in addition to the DMRG, methods for the exact diagonalization of the Hamiltonian of the system, like the Lanczos- or the Jacobi-Davidson method, are treated. An extension of the Lanczos method makes it possible to treat the time evolution of strongly correlated quantum systems with an accuracy comparable to machine precision. This method is the basis for a possible extension of the DMRG for the treatment of systems out of equilibrium, the so-called adaptive time-dependent DMRG ("adaptive t-DMRG"). A second variant of the adaptive t-DMRG uses the Suzuki-Trotter decomposition of the time-evolution operator. An error analysis demonstrates that both methods, for a suitable choice of control parameters, have errors < 1% at the end of the time evolution in complicated quantities like, e.g., the momentum distribution function. We apply these numerical methods to investigate the quantum critical behavior of a variant of the Hubbard model and to treat two non-equilibrium situations of current interest. Extensive use of the DMRG makes it possible to clarify the quantum-critical behavior of the ionic Hubbard model; in particular, the numerical results demonstrate quite clearly that the correct scenario has two critical points; at the first critical point, only the charge degrees of freedom are critical, while at the second one only the spin degrees of freedom are critical. Quite surprisingly, the Mott-insulator phase in the strong-coupling regime shows a divergent electrical susceptibility. Next, the dynamics of a system of so-called soft-core bosons is treated. The particles initially are trapped in a deep box-potential and the dynamics is investigated after releasing them from this trapping potential. Similar to exact results obtained for so-called hard-core bosons, (quasi-)coherent matter waves emerge, demonstrating the possibility of realizing an atom laser in experiments on optical lattices. We find that the wave vector of the emerging matter wave can be tuned by changing the strength of the interaction between the particles. The second non-equilibrium situation is a so-called "quantum-quench" of a system of strongly correlated fermions, i.e., the evolution of the system after suddenly changing an intrinsic parameter like, e.g., the interaction strength between the particles, is investigated. In this case, we are mainly interested in the long-time behavior; in particular, the question arises whether, relying on Boltzmann's ergodic hypothesis, relaxation to a thermal state is obtained. We find that, in general, non-thermal final states are reached, which can be described by generalized Gibbs-Boltzmann ensembles.
We study the dynamics in low dimensional correlated quantum systems. Advances in this thesis are three-fold. First, as a methodical point we extend the real-time formulation of the functional renormalization group approach introduced in [1]. Describing the transient time dependence of a given microscopic model prepared in an initial non-equilibrium state poses a formidable challenge in the presence of interaction. We successfully extend the functional renormalization group approach of [1] to finite temperature as well as extended interacting geometries. Furthermore, for both extensions we test the validity of the perturbatively motivated treatment by explicitly benchmarking to complementary methods and find that for not too large interactions the functional renormalization group proves to be quantitatively accurate. Second, we examine the non-equilibrium dynamics of open as well as closed quantum systems in the presence of correlations. We study a prominent prototype model for coherence in small quantum systems coupled to dissipative environments, namely the spin-boson model. We reveal that prevailing literature does not cover the full complexity of the dynamics encountered. The relaxation dynamics formerly classified as coherent and incoherent should be refined according to the results presented in this thesis: we propose a subdivision of the coherent regime into partially and asymptotically coherent. For the former the coherence signature of the dynamics is restricted only to intermediate times, while in the latter it extends in the asymptotic regime. Interestingly, we explicitly show that elevating temperature might be conducive to coherence, which is opposite to the expected behavior. For those closed quantum systems known to fall into the so called Luttinger Liquid universality class in equilibrium, we provide evidence that the equilibrium concept of universality seemingly carries over to the non-equilibrium steady-state. This is rather counter-intuitive, as universality is usually attributed to the low energy degrees of freedom only. Third, we study the thermoelectric transport properties of quantum dots. Quantum dots are promising candidates for highly efficient thermoelectric devices as they are not bounded by the Wiedemann-Franz law which strongly limits the efficiency of ordinary bulk systems. For this we study a charge fluctuating model and a model which additionally allows for spin fluctuations and phononic degrees of freedom. We find that in certain regimes quantum dots show very high efficiency when operated as power engines. [1] D. M. Kennes, ‘A Renormalization Group Approach to the Time Evolution of Correlated Quantum Dots’ Master’s thesis, RWTH Aachen (2011).
Recent experimental achievements in controlling ultracold gases in optical lattices open a new perspective on quantum many-body physics. In these experimental setups, it is possible to study coherent time evolution of isolated quantum systems. These dynamics reveal new physics beyond the low-energy properties that are usually relevant in solid-state many-body systems. In this paper, we study the time evolution of antiferromagnetic order in the Heisenberg chain after a sudden change of the anisotropy parameter, using various numerical and analytical methods. As a generic result, we find that the order parameter, which can show oscillatory or non-oscillatory dynamics, decays exponentially except for the effectively non-interacting case of the XX limit. For weakly ordered initial states, we also find evidence for an algebraic correction to the exponential law. The study is based on numerical simulations using a numerical matrix product method for infinite system sizes (iMPS), for which we provide a detailed description and an error analysis. Additionally, we investigate in detail the exactly solvable XX limit. These results are compared to approximative analytical approaches including an effective description by the XZ model as well as by mean-field, Luttinger-liquid and sine-Gordon theories. The comparison reveals which aspects of non-equilibrium dynamics can, as in equilibrium, be described by low-energy theories and which are the novel phenomena specific to quantum quench dynamics. The relevance of the energetically high part of the spectrum is illustrated by means of a full numerical diagonalization of the Hamiltonian.
The study of quantum dynamics is an increasingly important aspect of atomic, molecular and optical physics, as well as condensed-matter physics and other disciplines. It is essential in order to understand non-equilibrium phenomena which do not fit within the framework of equilibrium thermodynamics, such as many-body localization, dynamical quantum phase transitions, and non-equilibrium steady states in driven-dissipative systems. However, the exact calculation of dynamical observables is extremely difficult in general for several reasons, which makes it challenging for theoretical predictions and analyses to keep pace with rapidly advancing experiments. In this thesis I use and develop several computational methods to investigate the dynamics of quantum many-body systems, with the two goals of better understanding the evolution of local observables and correlation functions in spin and Fermion systems as well as creating and testing more effective numerical methods for solving quantum dynamics problems. I explore quantum dynamics by using Hamiltonian parameter quenches in several commonly-used models in AMO physics and elsewhere: the Fermi-Hubbard model, the XXZ model, and the transverse-Ising model. The first represents a system of short-range interacting Fermions, while the latter two represent spin systems with different spin-spin couplings and magnetic field terms. All of these models can be realized in current experiments. First, I derive an exact solution for an extreme interaction quench in the Fermi-Hubbard model, from a strongly-interacting initial state to a noninteracting final Hamiltonian. From this result I show that there are nontrivial transient connected correlations, despite the fact that the system has very high temperature and a noninteracting Hamiltonian. Those features would suppress such correlations in equilibrium, but the dynamical correlations are still significant. Then, I develop and examine two numerical methods for solving dynamics of the spin expectation values and correlation functions in the transverse-Ising and XXZ models. Both methods rely on approximating a thermodynamically large system with a much smaller and more computationally tractable cluster of lattice sites, or a series of clusters. The first method is a dynamical numerical linked cluster expansion (NLCE). NLCEs have already been used successfully on systems in thermodynamic equilibrium, and I extended the technique to directly compute the dynamics of an observable after a quantum quench. I review the general NLCE method, describe my dynamical variant of it, and compute the quench dynamics of spin and correlation observables in the transverse-Ising and XXZ models. The second method, which I term the adaptive-boundary cluster method, relies on recently refined bounds on information propagation in quantum systems to choose a cluster of lattice sites that minimizes the finite size error bound of a dynamical computation. I briefly review these bounds, and describe a new algorithm which uses them to construct a cluster whose boundary is chosen according to a related error estimate. Such an adaptive-boundary cluster has the property that there exists a computable upper bound on the finite size error of observables measured at the cluster center, and this error bound is nearly minimal when compared to that of other clusters of the same size. I then compute the quench dynamics of spin and correlation observables in the transverse-Ising model using the adaptive boundary clusters. I apply this method to both uniform and strongly disordered systems. To evaluate the effectiveness of both dynamical NLCE and adaptive-boundary clusters, I compare the results obtained using them to results from a standard and general purpose numerical method: exact solvers applied a rectangular cluster with periodic boundary conditions. Both of the new methods usually converge faster and are more accurate than the standard method for comparable cluster sizes. This thesis illustrates that choosing numerical methods which take advantage of the geometric and entanglement structure of the initial state and Hamiltonian can substantially enhance our ability to study quantum dynamics, by reducing the computational cost of numerically exact solutions. I conclude with a discussion of various avenues for future research to build on this work.
Recent experimental advances on ultracold atomic gases and trapped ions have made it possible to simulate exactly solvable quantum systems of interacting particles. In particular, the feasibility of making rapid changes, so-called quantum quenches, to such set-ups has allowed experimentalists to probe non-equilibrium phenomena in closed interacting quantum systems. This, in turn, has spurred a considerable theoretical interest in quantum many-body systems out of equilibrium. In this thesis, we study non-equilibrium properties of quantum many-body systems in the framework of exactly solvable quantum field theory in one spatial dimension. Specific systems include interacting fermions described by the Luttinger model and effective descriptions of spin chains using conformal field theory (CFT). Special emphasis is placed on heat and charge transport, studied from the point of view of quench dynamics, and, in particular, the effects of breaking conformal symmetries on transport properties. Examples include the Luttinger model with non-local interactions, breaking Lorentz and scale invariance, and inhomogeneous CFT, which generalizes standard CFT in that the usual propagation velocity v is replaced by a function v(x) that depends smoothly on the position x, breaking translation invariance. The quench dynamics studied here is for quantum quenches between, in general, different smooth inhomogeneous systems. An example of this is the so-called smooth-profile protocol, in which the initial state is defined by, e.g., smooth inhomogeneous profiles of inverse temperature and chemical potential, and the time evolution is governed by a homogeneous Hamiltonian. Using this protocol, we compute exact analytical results for the full time evolution of the systems mentioned above. In particular, we derive finite-time results that are universal in the sense that the same relations between the non-equilibrium dynamics and the initial profiles hold for any unitary CFT. These results also make clear that heat and charge transport in standard CFT are purely ballistic. Finally, we propose and study an inhomogeneous CFT model with v(x) given by a random function. We argue that this model naturally emerges as an effective description of one-dimensional quantum many-body systems with certain static random impurities. Using tools from wave propagation in random media, we show that such impurities lead to normal and anomalous diffusive contributions to heat transport on top of the ballistic one known from standard CFT.
We study non-equilibrium and strongly-correlated dynamics in two contexts. We begin by analysing quantum many-body systems out of equilibrium through the lens of cold atomic impurities in Bose gases. Such highly-imbalanced mixtures provide a controlled arena for the study of interactions, dissipation, decoherence and transport in a many-body quantum environment. Specifically we investigate the oscillatory dynamics of a trapped and initially highly-localised impurity interacting with a weakly-interacting trapped quasi low-dimensional Bose gas. This relates to and goes beyond a recent experiment by the Inguscio group in Florence. We witness a delicate interplay between the self-trapping of the impurity and the inhomogeneity of the Bose gas, and describe the dissipation of the energy of the impurity through phononic excitations of the Bose gas. We then study the transport of a driven, periodically-trapped impurity through a quasi one-dimensional Bose gas. We show that placing the weakly-interacting Bose gas in a separate periodic potential leads to a phononic excitation spectrum that closely mimics those in solid state systems. As a result we show that the impurity-Bose gas system exhibits phonon-induced resonances in the impurity current that were predicted to occur in solids decades ago but never clearly observed. Following this, allowing the bosons to interact strongly, we predict the effect of different strongly-correlated phases of the Bose gas on the motion of the impurity. We show that, by observing the impurity, properties of the excitation spectrum of the Bose gas, e.g., gap and bandwidth, may be inferred along with the filling of the bosonic lattice. In other words the impurity acts as a probe of its environment. To describe the dynamics of such a strongly-correlated system we use the powerful and near-exact time-evolving block decimation (TEBD) method, which we describe in detail. The second part of this thesis then analyses, for the first time, the performance of this method when applied to simulate non-equilibrium classical stochastic processes. We study its efficacy for a well-understood model of transport, the totally-asymmetric exclusion process, and find it to be accurate. Next, motivated by the inefficiency of sampling-based numerical methods for high variance observables we adapt and apply TEBD to simulate a path-dependent observable whose variance increases exponentially with system size. Specifically we calculate the expected value of the exponential of the work done by a varying magnetic field on a one-dimensional Ising model undergoing Glauber dynamics. We confirm using Jarzynski's equality that the TEBD method remains accurate and efficient. Therefore TEBD and related methods complement and challenge the usual Monte Carlo-based simulators of non-equilibrium stochastic processes.