Dynamics of interacting bosons using the Herman-Kluk semiclassical initial value representation

Research output: Contribution to journalResearch articleContributedpeer-review

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 (n -> infinity). An error is, however, incurred in next-to-next-to-leading order (small parameter 1/n), 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.

Details

Original languageEnglish
Article number165303
Number of pages18
JournalJournal of Physics A: Mathematical and Theoretical
Volume49
Issue number16
Publication statusPublished - 22 Apr 2016
Peer-reviewedYes

External IDs

WOS 000372195600015
Scopus 84961595312
ORCID /0000-0002-7806-3525/work/142234190

Keywords

Keywords

  • interacting bosons, semiclassical dynamics, semiclassial initial value representation, ultracold gases in optical lattices