Theoretical descriptions of the non-equilibrium dynamics of quantum many-body systems essentially employ either (i) explicit treatments, relying on the truncation of the expansion of the many-body wave function, (ii) compressed representations of the many-body wave function, or (iii) evolution of an effective (downfolded) representation through Green's functions. In this work, we select representative cases of each of the methods and address how these complementary approaches capture the dynamics driven by intense field perturbations to non-equilibrium states. Under strong driving, the systems are characterized by strong entanglement of the single-particle density matrix and natural populations approaching those of a strongly interacting equilibrium system. We generate a representative set of results that are numerically exact and form a basis for a critical comparison of the distinct families of methods. We demonstrate that the compressed formulation based on similarity-transformed Hamiltonians (coupled-cluster approach) is practically exact in weak fields and, hence, weakly or moderately correlated systems. Coupled cluster, however, struggles for strong driving fields, under which the system exhibits strongly correlated behavior, as measured by the von Neumann entropy of the single-particle density matrix. The dynamics predicted by Green's functions in the (widely popular) G W approximation are less accurate, but improve significantly upon the mean-field results in the strongly driven regime.
Theoretical descriptions of the non-equilibrium dynamics of quantum many-body systems essentially employ either (i) explicit treatments, relying on the truncation of the expansion of the many-body wave function, (ii) compressed representations of the many-body wave function, or (iii) evolution of an effective (downfolded) representation through Green's functions. In this work, we select representative cases of each of the methods and address how these complementary approaches capture the dynamics driven by intense field perturbations to non-equilibrium states. Under strong driving, the systems are characterized by strong entanglement of the single-particle density matrix and natural populations approaching those of a strongly interacting equilibrium system. We generate a representative set of results that are numerically exact and form a basis for a critical comparison of the distinct families of methods. We demonstrate that the compressed formulation based on similarity-transformed Hamiltonians (coupled-cluster approach) is practically exact in weak fields and, hence, weakly or moderately correlated systems. Coupled cluster, however, struggles for strong driving fields, under which the system exhibits strongly correlated behavior, as measured by the von Neumann entropy of the single-particle density matrix. The dynamics predicted by Green's functions in the (widely popular) G W approximation are less accurate, but improve significantly upon the mean-field results in the strongly driven regime.
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.
A numerical approach is presented that allows to compute nonequilibrium steady-state properties of strongly correlated quantum many-body systems. The method is imbedded in the Keldysh Green's function formalism and is based upon the idea of the variational cluster approach as far as the treatment of strong correlations is concerned. It appears that the variational aspect is crucial as it allows for a suitable optimization of a ``reference'' system to the nonequilibrium target state. The approach is neither perturbative in the many-body interaction nor in the field, that drives the system out of equilibrium, and it allows to study strong perturbations and nonlinear responses of systems in which also the correlated region is spatially extended. We apply the presented approach to nonlinear transport across a strongly correlated quantum wire described by the fermionic Hubbard model. We illustrate how the method bridges to cluster dynamical mean-field theory upon coupling two baths containing an increasing number of uncorrelated sites.
Ab-initio calculations are dominated by density functional theory (DFT), which provides an efficient and accurate description of the electronic structure for most materials. For materials with strong correlations, however, many of the available density functionals yield poor results. Most well known is the case of transition metal oxides, for which most density functionals produce a qualitatively incorrect description. There is a quest to improve the description by borrowing from methods specifically designed for strongly correlated materials. One of the method which is widely used to study strongly correlated electron system in the thermodynamic limit is the dynamical mean-field theory (DMFT). For calculations pertaining to real materials, hybrid approaches like DFT+DMFT are commonly used. We consider reduced density-matrix functional theory (rDMFT) to be a useful framework for a rigorous formulation of such hybrid theories. The link from rDMFT to many-particle wave functions has been established by Levy's constrained-search algorithm on the one hand. The link to many-body perturbation theory and Green's function, on the other hand, has been provided via the Luttinger-Ward functional. The development in the field of rDMFT proceeded analogously to the development of DFT. In order to avoid the full complexity of an explicit many-body description, most density-matrix functionals are not extracted from the exact expressions i.e. via Levy's constrained-search or Green's function based methods . Rather, one proceeds analogously to the development of density functionals, namely by searching models for the density-matrix functional, that capture the most essential physical effects while having an algebraic dependence on the density matrix. Among such models, the most prominent ones are the Müller functional, Sharma functional and Marques-Lathiotakis functional. One of the major arguments in favor of density-matrix functionals is that one of the most simple functionals, the Müller functional, seems to provide a correct description of the bond-dissociation problem, for which common density functionals fail. Despite the successes, the common density-matrix functionals also reproduce a number of features in a qualitatively incorrect manner. Several such discrepancies of the standard state of the art density-matrix functionals are discussed in details in this thesis. The objective of this thesis is to study the performance of the class of commonly used model density-matrix functionals and to set the stage for the desire to develop the framework for an approximate scheme to evaluate new class of functionals from exact Green's function based method.
In the previous chapter we studied the kinematics of many-body systems, and the form of operators representing the physical properties of a system, all of which were embodied by the quantum field. In this chapter we shall study the quantum dynamics of many-body systems, which can also be embodied by the quantum fields. We shall employ the fact that the quantum dynamics of a system, instead of being described in terms of the dynamics of the states or the evolution operator, i.e. as previously done through the Schrödinger equation, can instead be carried by the quantum fields. The quantum dynamics is therefore expressed in terms of the correlation functions or Green's functions of the quantum fields evaluated with respect to some state of the system. In particular we shall consider the general case of quantum dynamics for arbitrary non-equilibrium states. After introducing various types of Green's functions and relating them to measurable quantities, we will discuss the simplifications reigning for the special case of equilibrium states. Quantum dynamics Quantum dynamics can be described in different ways since quantum mechanics is a linear theory and the dynamics described by a unitary transformation of states. This will come in handy in the next chapter when we study a quantum theory in terms of its perturbative expansion using the so-called interaction picture. Here we first discuss the Schrödinger and Heisenberg pictures. The Schrödinger picture Having the Hamiltonian on the multi-particle space at hand we can consider the dynamics described in the multi-particle space.
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.
In this thesis we study three examples of interacting many-body systems undergoing a non equilibrium time evolution. Firstly we consider the time evolution in an integrable system: the sine-Gordon field theory in the repulsive regime. We will focus on the one point function of the semi-local vertex operator eiβφ(x)/2 on a specific class of initial states. By analytical means we show that the expectation value considered decays exponentially to zero at late times and we determine the decay time. The method employed is based on a form-factor expansion and uses the "Representative Eigenstate Approach" of Ref. [73] (a.k.a. "Quench Action"). In a second example we study the time evolution in models close to "special" integrable points characterised by hidden symmetries generating infinitely many local conservation laws that do not commute with one another, in addition to the infinite commuting family implied by integrability. We observe that both in the case where the perturbation breaks the integrability and when it breaks only the additional symmetries maintaining integrability, the local observables show a crossover behaviour from an initial to a final quasi stationary plateau. We investigate a weak coupling limit, identify a time window in which the effects of the perturbations become significant and solve the time evolution through a mean-field mapping. As an explicit example we study the XYZ spin-1/2 chain with additional perturbations that break integrability. Finally, we study the effects of integrability breaking perturbations on the non-equilibrium evolution of more general many-particle quantum systems, where the unperturbed integrable model is generic. We focus on a class of spinless fermion models with weak interactions. We employ equation of motion techniques that can be viewed as generalisations of quantum Boltzmann equations. We benchmark our method against time dependent density matrix renormalisation group computations and find it to be very accurate as long as interactions are weak. For small integrability breaking, we observe robust prethermalisation plateaux for local observables on all accessible time scales. Increasing the strength of the integrability breaking term induces a "drift" away from the prethermalisation plateaux towards thermal behaviour. We identify a time scale characterising this crossover.