Results are presented for the dynamics arising due to a sudden quench of a boson interaction parameter with the simultaneous switching on of a commensurate periodic potential, the latter providing a source of nonlinearity that can cause inelastic scattering. A quantum kinetic equation is derived perturbatively in the periodic potential and solved within the leading order gradient expansion. A two-particle irreducible formalism is employed to construct the stress-momentum tensor and hence the conserved energy. The dynamics is studied in detail in the phase where the boson spectrum remains gapless. The periodic potential is found to give rise to multiparticle scattering processes that relaxes the boson distribution function. At long times the system is found to thermalize with a thermalization time that depends in a nonmonotonic way on the amount of energy injected into the system due to the quantum quench. This nonmonotonic behavior arises due to the competing effect of an increase of phase space for scattering on the one hand, and an enhancement of the orthogonality catastrophe on the other hand as the quench amplitude is increased. The approach to equilibrium is found to be purely exponential for large quench amplitudes, and more complex for smaller quench amplitudes.
Results are presented for the dynamics arising due to a sudden quench of a boson interaction parameter with the simultaneous switching on of a commensurate periodic potential, the latter providing a source of nonlinearity that can cause inelastic scattering. A quantum kinetic equation is derived perturbatively in the periodic potential and solved within the leading order gradient expansion. A two-particle irreducible formalism is employed to construct the stress-momentum tensor and hence the conserved energy. The dynamics is studied in detail in the phase where the boson spectrum remains gapless. The periodic potential is found to give rise to multiparticle scattering processes that relaxes the boson distribution function. At long times the system is found to thermalize with a thermalization time that depends in a nonmonotonic way on the amount of energy injected into the system due to the quantum quench. This nonmonotonic behavior arises due to the competing effect of an increase of phase space for scattering on the one hand, and an enhancement of the orthogonality catastrophe on the other hand as the quench amplitude is increased. The approach to equilibrium is found to be purely exponential for large quench amplitudes, and more complex for smaller quench amplitudes.
The dynamics of interacting bosons in one dimension following the sudden switching on of a weak disordered potential is investigated. On time scales before quasiparticles scatter (prethermalized regime), the dephasing from random elastic forward scattering causes all correlations to decay exponentially fast, but the system remains far from thermal equilibrium. For longer times, the combined effect of disorder and interactions gives rise to inelastic scattering and to thermalization. A novel quantum kinetic equation accounting for both disorder and interactions is employed to study the dynamics. Thermalization turns out to be most effective close to the superfluid-Bose-glass critical point where nonlinearities become more and more important. The numerically obtained thermalization times are found to agree well with analytic estimates.
We investigate the dynamics of one-dimensional interacting bosons in an optical lattice after a sudden quench in the weakly interacting (Bose-Hubbard) and strongly interacting (sine-Gordon) regimes. While in a higher dimension, the Mott-superfluid phase transition is observed for weakly interacting bosons in deep lattices, in one dimension an instability is generated also for shallow lattices with a commensurate periodic potential pinning the atoms to the Mott state through a transition described by the sine-Gordon model. The present work aims at identifying the quench dynamics in both the Bose-Hubbard and sine-Gordon interaction regimes. We numerically exactly solve the time-dependent Schr\"odinger equation for a small number of atoms and obtain dynamical measures of several key quantities. We investigate the correlation dynamics of first and second order; both exhibit rich many-body features in the dynamics. We conclude that in both cases, dynamics exhibits collapse-revival phenomena, though with different timescales. We argue that the dynamical fragmentation is a convenient quantity to distinguish the dynamics especially near the pinning zone. To understand the relaxation process we measure the many-body information entropy. Bose-Hubbard dynamics clearly establishes the possible relaxation to the maximum entropy state. In contrast, the sine-Gordon dynamics is so fast that it does not exhibit any signature of relaxation in the present timescale of computation.
Abstract We study the tunneling dynamics of N = 10 one-dimensional interacting bosons confined in a temporally driven double well potential that imitates a quantum seesaw and how we can manipulate these dynamics by changing the drive of the seesaw potential. We emulate the seesaw with a driven double well potential and consider two driving protocols: an harmonic constant-frequency drive and a chirped drive with linearly increasing frequency. We consider the time-dependent many-body Schrödinger equation of a repulsively interacting quasi-one-dimensional few-boson system. We solve it by using the multiconfigurational time-dependent Hartree method for bosons (MCTDHB) as implemented in the MCTDH-X software. For an harmonic drive and at small values of the driving amplitude, the dynamics of the particles become very slow rendering a stationary-like state. In a phase-space picture the population imbalance between the wells follows a trajectory which is restricted to a comparatively small region of space. For an harmonic drive at intermediate amplitudes, the dynamics become periodic in nature, implying that the bosons populate each of wells periodically. At comparatively large amplitudes of the harmonic drive, the dynamics show features of chaos in phase-space representation. For the chirped drive with a driving frequency increasing linearly in time, the imbalance of the atoms in the seesaw, however, has a temporal evolution that is faster for certain frequency ranges. The tunneling dynamics in such cases, for small amplitudes, show the appearance of quasi-periodicity with simultaneously present slow and fast oscillations. Increasing the amplitude of the chirped drive, we observe that the dynamics, although being periodic, become severely damped in their amplitude. Our study establishes that by tuning the temporal evolution of the quantum seesaw, a precise control of tunneling dynamics of the correlated bosons can be achieved. Since harmonic driving and chirp frequency modulation of the seesaw are experimentally achievable, our simulations can be experimentally realized in laboratories dealing with cold atomic gases.
A perturbative renormalization group approach is employed to study the effect of a periodic potential on a system of one-dimensional bosons in a nonequilibrium steady state due to an initial interaction quench. The renormalization group flows are modified significantly from the well-known equilibrium Berezinski-Kosterlitz-Thouless form. They show several new features such as a generation of an effective temperature, generation of dissipation, as well as a change in the location of the quantum critical point separating the weak- and strong-coupling phases. Detailed results on the weak-coupling side of the phase diagram are presented, such as the renormalization of the parameters and the asymptotic behavior of the correlation functions. The physical origin of the generated temperature and friction is discussed.
We investigate the quantum dynamics of two bosons, trapped in a two-dimensional harmonic trap, upon quenching arbitrarily their interaction strength and thereby covering the entire energy spectrum. Utilizing the exact analytical solution of the stationary system, we derive a closed analytical form of the expansion coefficients of the time-evolved two-body wave function, whose dynamics is determined by an expansion over the postquench eigenstates. The emergent dynamical response of the system is analyzed in detail by inspecting several observables such as the fidelity, the reduced one-body densities, the radial probability density of the relative wave function in both real and momentum space, and the Tan contact, which unveils the existence of short range two-body correlations. When the system is initialized in its bound state, it is perturbed in the most efficient manner as compared to any other initial configuration. Moreover, starting from an interacting ground state, the two-boson response is enhanced for quenches toward the noninteracting limit.
Motivated by recent experiments, we study the relaxation dynamics and thermalization in the one-dimensional Bose-Hubbard model induced by a global interaction quench. Specifically, we start from an initial state that has exactly one boson per site and is the ground state of a system with infinitely strong repulsive interactions at unit filling. Using exact diagonalization and the density-matrix renormalization-group method, we compute the time dependence of such observables as the multiple occupancy and the momentum distribution function. Typically, the relaxation to stationary values occurs over just a few tunneling times. The stationary values are identical to the so-called diagonal ensemble on the system sizes accessible to our numerical methods and we further observe that the microcanonical ensemble describes the time averages of many observables reasonably well for small and intermediate interaction strength. The expectation values of observables in the canonical ensemble agree quantitatively with the time averages obtained from the quench at small interaction strengths, and qualitatively provide a good description even in parameter regimes where the microcanonical ensemble is not applicable due to finite-size effects. We discuss our numerical results in the framework of the eigenstate thermalization hypothesis. Moreover, we also observe that the diagonal and the canonical ensembles are practically identical for our initial conditions already on the level of their respective energy distributions for small interaction strengths. Finally, we discuss implications of our results for the interpretation of a recent sudden expansion experiment [Phys. Rev. Lett. 110, 205301 (2013)], in which the same interaction quench was realized.