We study the out-of-equilibrium dynamics of ultracold bosons in a double- and triple-well potential within the Bose-Hubbard model by means of the semiclassical Herman-Kluk propagator and compare the results to the frequently applied ``classical dynamics'' calculation in terms of the truncated Wigner approximation (TWA). For the double-well system we find the semiclassical results in excellent agreement with the numerically exact ones, while the TWA is not able to reproduce any revivals of the wave function. The triple-well system turns out to be more difficult to handle due to the irregularity of the corresponding classical phase space. Here, deviations of the TWA from the exact dynamics appear even for short times, while better agreement is obtained using the semiclassical approach presented in this article.
We study the out-of-equilibrium dynamics of ultracold bosons in a double- and triple-well potential within the Bose-Hubbard model by means of the semiclassical Herman-Kluk propagator and compare the results to the frequently applied ``classical dynamics'' calculation in terms of the truncated Wigner approximation (TWA). For the double-well system we find the semiclassical results in excellent agreement with the numerically exact ones, while the TWA is not able to reproduce any revivals of the wave function. The triple-well system turns out to be more difficult to handle due to the irregularity of the corresponding classical phase space. Here, deviations of the TWA from the exact dynamics appear even for short times, while better agreement is obtained using the semiclassical approach presented in this article.
Abstract Recent experimental progress in monitoring the dynamics of ultracold gases in optical lattices necessitates a quantitative theoretical description for a significant number of bosons. In the present paper, we investigate if time-dependent semiclassical initial value methodology, with propagators expressed as integrals over phase space and using classical trajectories, is suitable to describe interacting bosons, concentrating on a single mode. Despite the nonlinear contribution from the self-interaction, the corresponding classical dynamics allows for a largely analytical treatment of the semiclassical propagator. We find that application of the Herman–Kluk (HK) propagator conserves unitarity in the semiclassical limit, but a decay of the norm is seen for low particle numbers n . The frozen Gaussian approximation (FGA) (HK with unit prefactor) is explicitly shown to violate unitarity in the present system for non-vanishing interaction strength, even in the semiclassical limit. Furthermore, we show by evaluating the phase space integral in steepest descent approximation, that the HK propagator reproduces the exact spectrum correctly in the semiclassical limit ( ). An error is, however, incurred in next-to-next-to-leading order (small parameter ), as seen upon numerical evaluation of the integral and confirmed analytically by considering finite n corrections to the steepest descent calculations. The FGA, in contrast, is only accurate to lowest order, and an erroneous next-to-leading order term in the energy spectrum was found analytically. Finally, as an example application, we study the dynamics of wave packets by computing the time evolution of the Wigner function. While the often-used truncated Wigner approximation cannot capture any interferences present in the exact quantum mechanical solution (known analytically), we find that the HK approach, despite also using classical information only, reproduces the salient features of the exact solution correctly.
Abstract We study the phase-space representation of dynamics of bosons in the semiclassical regime where the occupation number of the modes is large. To this end, we employ the van Vleck-Gutzwiller propagator to obtain an approximation for the Green’s function of the Wigner distribution. The semiclassical analysis incorporates interference of classical paths and reduces to the truncated Wigner approximation (TWA) when the interference is ignored. Furthermore, we identify the Ehrenfest time after which the TWA fails. As a case study, we consider a single-mode quantum nonlinear oscillator, which displays collapse and revival of observables. We analytically show that the interference of classical paths leads to revivals, an effect that is not reproduced by the TWA or a perturbative analysis.
We study the phase-space representation of dynamics of bosons in the semiclassical regime where the occupation number of the modes is large. To this end, we employ the van Vleck-Gutzwiller propagator to obtain an approximation for the Green's function of the Wigner distribution. The semiclassical analysis incorporates interference of classical paths and reduces to the truncated Wigner approximation (TWA) when the interference is ignored. Furthermore, we identify the Ehrenfest time after which the TWA fails. As a case study, we consider a single-mode quantum nonlinear oscillator, which displays collapse and revival of observables. We analytically show that the interference of classical paths leads to revivals, an effect that is not reproduced by the TWA or a perturbative analysis.
We study the dynamics of ultracold bosons in a double-well potential within the two-mode Bose-Hubbard model by means of semiclassical methods. By applying a Wentzel-Kramers-Brillouin (WKB) quantization we find analytical results for the energy spectrum, which are in excellent agreement with numerically exact results. They are valid in the energy range of plasma oscillations, both in the Rabi and the Josephson regime. Adopting the reflection principle and the Poisson summation formula we derive an analytical expression for the dynamics of the population imbalance depending only on the few relevant parameters of the system. This allows us to discuss its characteristic dynamics, especially the oscillation frequency and the collapse and revival time, as a function of the model parameters, leading to a deeper understanding of Josephson physics. We find that our formulas match previous experimental observations.
We develop and utilize the SU(3) truncated Wigner approximation (TWA) in order to analyze far-from-equilibrium quantum dynamics of strongly interacting Bose gases in an optical lattice. Specifically, we explicitly represent the corresponding Bose-Hubbard model at an arbitrary filling factor with restricted local Hilbert spaces in terms of SU(3) matrices. Moreover, we introduce a discrete Wigner sampling technique for the SU(3) TWA and examine its performance as well as that of the SU(3) TWA with the Gaussian approximation for the continuous Wigner function. We directly compare outputs of these two approaches with exact computations regarding dynamics of the Bose-Hubbard model at unit filling with a small size and that of a fully connected spin-1 model with a large size. We show that both approaches can quantitatively capture quantum dynamics on a timescale of $\ensuremath{\hbar}/(Jz)$, where J and z denote the hopping energy and the coordination number. We apply the two kinds of SU(3) TWA to dynamical spreading of a two-point correlation function of the Bose-Hubbard model on a square lattice with a large system size, which has been measured in recent experiments. Noticeable deviations between the theories and experiments indicate that proper inclusion of effects of the spatial inhomogeneity, which is not straightforward in our formulation of the SU(3) TWA, may be necessary.
We analyze the time evolution of the Bose-Hubbard model after a sudden quantum quench to a weakly interacting regime. Specifically, motivated by a recent experiment at Kyoto University, we numerically simulate redistribution of the kinetic and on-site interaction energies at an early time, which was observed in nonequilibrium dynamics of ultracold Bose gases in a cubic optical lattice starting with a singly occupied Mott-insulator state. In order to compute the short-time dynamics corresponding to the experimental situation, we apply the truncated-Wigner approximation to the Bose-Hubbard model on a cubic lattice. We show that our semiclassical approach quantitatively reproduces the fast redistribution dynamics. We further analyze spatial spreading of density-density correlations at equal time in the Bose-Hubbard model on a square lattice with a large filling factor. When the system is initially prepared in a coherent state, we find that a propagation velocity of the correlation wave packet in the correlation function strongly depends on the final interaction strength, and it is bounded by twice the maximum group velocity of the elementary excitations. In contrast, when the system is initially in a Mott-insulator state, the propagation velocity of the wave packet is approximately independent of the final interaction strength.