Analytical results on the correlation functions of strongly correlated many-body systems are rare in the literature and their importance cannot be overstated. We present determinant representations for the space-, time-, and temperature-dependent correlation functions of the strongly interacting one-dimensional Hubbard model in the presence of an external trapping potential. These representations are exact and valid in both equilibrium and nonequilibrium scenarios like the ones initiated by a sudden change of the confinement potential. In addition, they can be implemented numerically very easily, significantly outperforming other numerical approaches. As applications of our results we investigate the single-particle spectral functions of systems with harmonic trapping and show that dynamical quasicondensation occurs for both fermionic and bosonic spin-$1/2$ systems released from a Mott insulator state.
Analytical results on the correlation functions of strongly correlated many-body systems are rare in the literature and their importance cannot be overstated. We present determinant representations for the space-, time-, and temperature-dependent correlation functions of the strongly interacting one-dimensional Hubbard model in the presence of an external trapping potential. These representations are exact and valid in both equilibrium and nonequilibrium scenarios like the ones initiated by a sudden change of the confinement potential. In addition, they can be implemented numerically very easily, significantly outperforming other numerical approaches. As applications of our results we investigate the single-particle spectral functions of systems with harmonic trapping and show that dynamical quasicondensation occurs for both fermionic and bosonic spin-1/2 systems released from a Mott insulator state.
We report on the derivation of determinant representations for the Green's functions and spectral function of the trapped Tonks-Girardeau gas on the lattice and in the continuum. Our results are valid for any type of statistics of the constituent particles, at zero and finite temperature and arbitrary confining potentials, including nonequilibrium scenarios induced by sudden changes of the external potential. In addition, they are also extremely efficient and easy to implement numerically with the main computational effort being represented by the calculation of partial overlaps of the dynamically evolved single-particle wave functions. In the lattice case we show that the spectral function of a system with a strong harmonic potential presents only two singular lines compared with three singular lines in the case of a homogeneous system.
Understanding the effect of correlations in interacting many-body systems is one of the main challenges in quantum mechanics. While the general problem can only be addressed by approximate methods and numerical simulations, in some limiting cases, it is amenable to exact solutions. I will present a family of exact solutions for strongly interacting 1D quantum gases in arbitrary external potential, valid for bosons, multicomponent fermions and mixtures. It allows us to obtain experimental observables such as the density profiles and momentum distribution at all momentum scales, beyond the Luttinger liquid approach. It also predicts the exact quantum dynamics at all the times, including the small oscillation regime yielding the collective modes of the system and the large quench regime where the system parameters are changed considerably. The solution can be extended to describe finite-temperature conditions, spin, and magnetization effects. In particular I will present large-size calculations of the frequency- and momentum-resolved spectral function of strongly repulsive lattice bosons and a symmetry-breaking exact solution for a Bose-Bose mixture. Our approach directly describes the experiments with trapped ultracold atoms where the strongly correlated regime in one dimension has been achieved.
In the first part of the Thesis we mostly concentrate on spectral properties of strongly correlated systems and on their equilibrium properties. This is accomplished by the general concept of imaginary-time dynamics which we apply to a number of different problems in which different strengths of this approach emerge. In Chapter 1 we introduce the formalism that allows for a connection between the quantum and the classical worlds. The connection is established by means of the imaginary-time quantum evolution which, under certain circumstances, is shown to be equivalent to a classical stochastic process. It is further shown that exact static and spectral properties of correlated systems can be obtained when this mapping is feasible. The relationship between the imaginary-time dynamics in different frameworks such as the path-integral and the perturbative one is also underlined. In Chapter 2 we present a specific implementation of the general ideas previously presented. In particular we introduced an extension to lattice systems of the Reptation Monte Carlo algorithm [30] which benefits of a sampling scheme based on directed updates. Specific improvements over the existing methodologies consist in the unbiased evaluation of the imaginary-time path integrals for bosons and a systematic scheme to improve over the Fixed-node approximation for fermions. Applications to the Hubbard and the Heisenberg models are presented. In Chapter 3 we demonstrate the application of the imaginary-time dynamics to the exact study of spectral properties. Subject of our attention is a highly anharmonic and correlated quantum crystal such as Helium 4 at zero temperature.[33] Concerning this system, we have obtained the first ab-initio complete phonon dispersion in good agreement with neutron spectroscopy experiments. Moreover, we have also studied the density excitations of solid helium in a region of wave-vectors in between the collective (phonon) and the single-particle regimes, where the presence of residual coherence in the dynamics shows analogies between the highly anharmonic crystal and the superfluid phase. In Chapter 4 we introduce a novel method, based on the imaginary-time dynamics, to obtain unbiased estimates of fermionic properties.[34] By means of this method and of a very accurate variational state, we provide strong evidence for the stability of a saturated ferromagnetic phase in the high-density regime of the two-dimensional infinite-U Hubbard model. By decreasing the electron density, we observe a discontinuous transition to a paramagnetic phase, accompanied by a divergence of the susceptibility on the paramagnetic side. This behavior, resulting from a high degeneracy among different spin sectors, is consistent with an infinite-order phase transition scenario. In Chapter 5 the use of imaginary-time dynamics in the context of finite-temperature response functions is highlighted. As an application, we study an intriguing quantum phase featuring both glassy order and Bose-Einstein condensation. [35] We introduce and validate a model for the role of geometrical frustration in the coexistence of off-diagonal long range order with an amorphous density profile. The exact characterization of the response of the system to an external density perturbation is what allows here to establish the existence of a spin-glass phase. The differences between such a phase and the otherwise insulating "Bose glasses" are further elucidated in the Chapter. In the second part of the Thesis we focus our attention on the dynamics of closed systems out of equilibrium. This is accomplished by both non-stochastic exact methods for the dynamics and the introduction of a novel time-dependent Variational Monte Carlo scheme. In Chapter 6 exact diagonalization schemes and renormalization-based methods for one-dimensional systems are introduced. We identify key phenomenological traits resulting from the many-body correlation in closed systems driven sufficiently away from equilibrium.[31] We provide evidences that the dynamics of interacting lattice bosons away from equilibrium can be trapped into extremely long-lived inhomogeneous metastable states. The slowing down of incoherent density excitations above a threshold energy, much reminiscent of a dynamical arrest on the verge of a glass transition, is identified as the key feature of this phenomenon. In Chapter 7 we present an extension to dynamical properties of the Variational Quantum Monte Carlo method.[32] This is accomplished by introducing a general class of time-dependent variational states which is based on the mapping of the many-body dynamics onto an instantaneous ground-state problem. The application of the method to the experimentally relevant quantum quenches of interacting bosons reveals the accuracy and the reliability of the introduced numerical scheme. We indeed obtain for the first time a consistent variational description of the approach to the equilibrium of local observables and underline the origin of the metastability and glassy behavior previously identified. In the very last part we draw our conclusions and show some possible paths for stimulating future research.
Although charge-spin separation has important consequences for the properties of one-dimensional (1D) Mott-Hubbard insulators, matrix elements in some of its dynamical correlation functions involve the coupling of spin and charge degrees of freedom. The corresponding interplay of the 1D Mott-Hubbard insulator's charge and spin degrees of freedom is an issue of both fundamental and technological interest, for instance, concerning their dynamics at subpicosecond timescales. On the other hand, the up- and down-spin one-particle spectral functions are the simplest dynamical correlation functions that involve excitation of both the charge and spin degrees of freedom. They are thus suitable to extract basic useful physical information on the above interplay at finite magnetic fields. Here the line shape of such functions at and in the $(k,\ensuremath{\omega})$ plane's vicinity of their cusp singularities is studied for the Mott-Hubbard insulator described by the 1D Hubbard model with one fermion per site at zero and finite magnetic fields. At zero field, they can be accessed in terms of electrons by photoemission experiments. At finite field, such functions and corresponding interplay of correlations and magnetic-field effects refer to an involved nonperturbative many-particle problem that is poorly understood. The Mott-Hubbard gap that separates the addition and removal spectral functions is calculated for all spin densities and interaction values. The qualitative differences in the one-particle properties of the Mott-Hubbard insulator and corresponding doped insulator are also investigated. The relation of our theoretical results and predictions to both condensed-matter and ultracold spin-1/2 atom systems is discussed.
We use nonequilibrium dynamical mean-field theory in combination with a recently developed Quantum Monte Carlo impurity solver to study the real-time dynamics of a Hubbard model which is driven out of equilibrium by a sudden increase in the on-site repulsion U. We discuss the implementation of the self-consistency procedure and some important technical improvements of the QMC method. The exact numerical solution is compared to iterated perturbation theory, which is found to produce accurate results only for weak interaction or short times. Furthermore, we calculate the spectral functions and the optical conductivity from a Fourier transform on the finite Keldysh contour, for which the numerically accessible time scales allow to resolve the formation of Hubbard bands and a gap in the strongly interacting regime. The spectral function, and all one-particle quantities that can be calculated from it, thermalize rapidly at the transition between qualitatively different weak- and strong-coupling relaxation regimes.
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.