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

Abstract

Recent experimental progress using 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 methods, with propagators expressed as integrals over phase space 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 the Herman-Kluk (HK) propagator conserves unitarity in the semiclassical limit ($n\to \infty$), but a decay of the norm is seen for low $n$. The frozen Gaussian approximation (FGA, i.e. 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 limit $n\to \infty$. An error is, however, incurred in next-to-next-to-leading order (small parameter $1/n$), as seen upon numerical evaluation of the integral or by considering finite $n$ corrections to steepest descent. 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.

arXiv comment18 pages, 4 figures; accepted for publication in present form by J. Phys. A: Math. Theor

References in corpus (12)

• Probing the relaxation towards equilibrium in an isolated strongly correlated 1D Bose gas
• The cold atom Hubbard toolbox
• Time-resolved Observation and Control of Superexchange Interactions with Ultracold Atoms in Optical Lattices
• Direct Observation of Second Order Atom Tunnelling
• Coherent Backscattering in Fock Space: a Signature of Quantum Many-Body Interference in Interacting Bosonic Systems
• Semiclassical quantization of an N-particle Bose-Hubbard model
• Semiclassical quantization of the Bogoliubov spectrum
• Analytical results for Josephson dynamics of ultracold Bosons
• The semiclassical propagator in fermionic Fock space
• Time-dependent Semiclassics for Ultracold Bosons
• Semiclassical Spectrum of Small Bose-Hubbard Chains: A Normal Form Approach
• Semiclassical quantisation for a bosonic atom-molecule conversion system

Cited by in corpus (14)

• The Importance of the Pre-exponential Factor in Semiclassical Molecular Dynamics
• Gouy's Phase Anomaly in Electron Waves Produced by Strong-Field Ionization
• Deep zoom-in simulation of a fuzzy dark matter galactic halo
• Chaos in the three-site Bose-Hubbard model -- classical vs quantum
• Quantum revival patterns from classical phase-space trajectories
• Eigenstate thermalization scaling in approaching the classical limit
• The semiclassical propagator in Fock space: dynamical echo and many-body interference
• Exact variational dynamics of the multimode Bose-Hubbard model based on SU(M) coherent states
• Simulation of the Quantum Dynamics of Indistinguishable Bosons with the Coupled Coherent States Method
• Quantum Enhancement of Thermalization
• The Dynamics of the Hubbard Model through Stochastic Calculus and Girsanov Transformation
• A semiclassical theory of phase-space dynamics of interacting bosons
• Computing Quantum Mean Values in the Deep Chaotic Regime
• Electronic energies from coupled fermionic 'Zombie' states imaginary time evolution