Expansion dynamics of interacting fermions in a lattice is simulated within the one-dimensional (1D) Hubbard model, using the essentially exact time-evolving block decimation (TEBD) method. In particular, the expansion of an initial band-insulator state is considered. We analyze the simulation results based on the dynamics of a two-site two-particle system, the so-called Hubbard dimer. Our findings describe essential features of a recent experiment on the expansion of a Fermi gas in a two-dimensional lattice. We show that the Hubbard-dimer dynamics, combined with a two-fluid model for the paired and nonpaired components of the gas, gives an efficient description of the full dynamics. This should be useful for describing dynamical phenomena of strongly interacting Fermions in a lattice in general.
Expansion dynamics of interacting fermions in a lattice is simulated within the one-dimensional (1D) Hubbard model, using the essentially exact time-evolving block decimation (TEBD) method. In particular, the expansion of an initial band-insulator state is considered. We analyze the simulation results based on the dynamics of a two-site two-particle system, the so-called Hubbard dimer. Our findings describe essential features of a recent experiment on the expansion of a Fermi gas in a two-dimensional lattice. We show that the Hubbard-dimer dynamics, combined with a two-fluid model for the paired and nonpaired components of the gas, gives an efficient description of the full dynamics. This should be useful for describing dynamical phenomena of strongly interacting Fermions in a lattice in general.
In this work we study dynamics of fermionic quantum gases. Radio frequency spectroscopy, Josephson oscillations, and time-of-flight dynamics are studied theoretically, and with the help of numerical simulations using the time-evolving block decimation algorithm. \n The simulations of radio frequency spectroscopy shed light on the nature of excitations in a fermionic superfluid confined in a lattice. The Josephson oscillations in the presence of different potentials for up and down spins show that the asymmetric Josephson effect cannot be described by the traditionally established view of coherent tunnelling of Cooper pairs. The expansion dynamics of a band-insulating state reveals that we can explain the time development of two-component systems in a lattice using a simple two-site model. Finally, the time-of-flight expansion of an imbalanced superfluid restricted to move in one dimension provides a remarkably simple way for detecting the long-sought-for Fulde-Ferrell-Larkin-Ovchinnikov superfluid state.
We investigate the expansion of bosons and fermions in a homogeneous lattice after a sudden removal of the trapping potential using exact numerical methods. As a main result, we show that in one dimension, both bosonic and fermionic Mott insulators expand with the same velocity, irrespective of the interaction strength, provided the expansion starts from the ground state of the trapped gas. Furthermore, their density profiles become identical during the expansion; the asymptotic density dynamics is identical to that of initially localized, noninteracting particles, and the asymptotic velocity distribution is flat. The expansion velocity for initial correlated Mott insulating states is therefore independent of the interaction strength and particle statistics. Interestingly, this nonequilibrium dynamics is sensitive to the interaction driven quantum phase transition in the Bose-Hubbard model; while being constant in the Mott phase, the expansion velocity decreases in the superfluid phase and vanishes for large systems in the noninteracting limit. These results are compared to the setup of a recent experiment [Ronzheimer et al., Phys. Rev. Lett. 110, 205301 (2013)], where the trap opening was combined with an interaction quench from infinitely strong interactions to finite values. In the latter case, the interaction quench breaks the universal dynamics in the asymptotic regime and the expansion depends on the interaction strength. We carry out an analogous analysis for a two-component Fermi gas, with similar observations. In addition, we study the effect of breaking the integrability of hard-core bosons in different ways; while the fast ballistic expansion from the ground state of Mott insulators in one dimension remains unchanged for finite interactions, we observe strong deviations from this behavior on a two-leg ladder even in the hard-core case. This change in dynamics bares similarities with the dynamics in the dimensional crossover from one to two dimensions observed in the aformentioned experimental study.
We simulate a balanced attractively interacting two-component Fermi gas in a one-dimensional lattice perturbed with a moving potential well or barrier. Using the time-evolving block decimation (TEBD) method, we study different velocities of the perturbation and distinguish two velocity regimes based on clear differences in the time evolution of particle densities and the pair correlation function. We show that, in the slow regime, the densities deform as particles are either attracted by the potential well or repelled by the barrier, and a wave front of hole or particle excitations propagates at the maximum group velocity. Simultaneously, the initial pair correlations are broken and coherence over different sites is lost. In contrast, in the fast regime, the densities are not considerably deformed and the pair correlations are preserved.
We study dynamical (quasi)condensation in the Fermi-Hubbard model starting from a completely uncorrelated initial state of adjacent doubly occupied sites. We show that upon expansion of the system in one dimension, dynamical (quasi)condensation occurs not only for large interactions via the condensation of doublons, but also for small interactions. The behavior of the system is distinctly different in the two parameter regimes, underlining a different mechanism at work. We address the question of whether the dynamical (quasi)condensation effect persists in the thermodynamic limit. For this purpose, we use the time-dependent two-particle reduced density matrix method, which allows the extension to large system sizes, long propagation times, and two-dimensional (2D) systems. Our results indicate that the effect vanishes in the thermodynamic limit. However, especially in 2D, further investigation beyond numerically tractable system sizes calls for the use of quantum simulators, for which we show that the described effect can be investigated by probing density fluctuations.
In the last years the impressive progress on the experimental side led to a variety of new experiments allowing to address systems out of equilibrium.In this way the behavior of such systems far from equilibrium is no longer a purely theoretical issue but indeed observable.New experimental techniques, like particles trapped in optical lattices, render a realization of quantum systems with nearly arbitrary system parameters possible and provide a possibility to study their time evolution.Systems out of equilibrium are characterized by the fact, that these systems are in highly excited states giving rise to totally new fascinating properties.In the present thesis one-and two-dimensional fermionic Hubbard models out of equilibrium are discussed.The system is taken out of equilibrium by a so-called interaction quench.At the beginning the system is prepared in the groundstate of the non-interacting Hamiltonian.At a time t the interaction between the fermions is suddenly turned on so that the time evolution is governed by the whole, interacting Hamiltonian.Hence the system is prepared in the groundstate of one Hamiltonian but evolves according to a di_erent Hamiltonian.Consequently the system ends up in a highly excited state.To describe such a system a method based on an expansion of the Heisenberg equations of motion to highest order possible is developed in this thesis.This method provides an exact description of the time evolution on short and intermediate time scales after the quench.As the method reveal exact results and does not rely on any perturbative assumption, a study of arbitrarily large interaction strengths is possible.Besides, the method is one of the few methods capable of two-dimensional systems.In the following the method used in this thesis is explained and advantages and disadvantages of the approach are thematized.For this purpose the results of the developed iterated equation of motion approach are compared to results obtained in calculation up to second order in the interaction.Then the dynamics of the onedimensional model is discussed with a focus on the relation of the results derived by the iterated equations of motion approach to results obtained by bosonization theory, the behavior for strong interactions and the dynamical transition from the weak to the strong quench regime.Furthermore, the Hubbard model is studied on a two-dimensional square lattice.This model is fundamentally di_erent from the one-dimensional model: In contrast to the one-dimensional model the two-dimensional model is not integrable allowing a true relaxation of the system.For this system a calculation up to second order in the interaction is performed and compared to the results of the iterated equation of motion approach.Besides, the time evolution of the momentum distribution and the influence of doping on the dynamics is studied.Moreover, a first estimate for relaxation times is provided without relying on the assumption of a mixed state.
The single-band Hubbard Hamiltonian [1] is the simplest model of interacting electrons in a lattice.It was developed by John Hubbard, who introduced the model to take into account the local electronic correlations in narrow energy bands.The model has a kinetic term that allows hopping of electrons between atomic sites and a term that considers the on-site interaction between electrons.This model is largely used in solid state physics to study magnetic properties of solids, insulator-metal transition (Mott transition) and high temperature superconductors.However, this is a problem that has exact analytical solution in a very few limiting cases and must be treated approximately or numerically.In this work, we propose a new methodology to analytically solve the Hubbard Hamiltonian, mapping it into a two-site model (the Hubbard dimer).To obtain the Green's function for the lattice, we employed the cumulant expansion technique, using as a 'seed' the exact two-site Green's function (the atomic method) [2].From the Green's function, we obtained the density of states [DOS] and the occupation numbers as functions of the external parameters of the model, and compared the results with the Hubbard I approximation [1].