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

DOI10.1088/1751-8113/49/16/165303MaRDI QIDQ3186631zbMATH OpenOpenAlexFDO

Authors Shouryya Ray, Paula Ostmann, Lena Simon, Frank Grossmann, Walter T. Strunz

Publication date 10 August 2016

Published in Journal of Physics A: Mathematical and Theoretical (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1512.07194

zbMATH Keywords

Herman-Klukinitial value representationinteracting ultracold Bose gasestime-dependent semiclassics

Mathematics Subject Classification ID

Numerical computation of eigenvalues and eigenvectors of matrices (65F15) Dynamic lattice systems (kinetic Ising, etc.) and systems on graphs in time-dependent statistical mechanics (82C20) Interacting particle systems in time-dependent statistical mechanics (82C22)

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 ($noinfty$), 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 $noinfty$. 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.

Recommendations

• From the coherent state path integral to a semiclassical initial value representation of the quantum mechanical propagator
• Comment on  Semiclassical approximations in phase space with coherent states
• Bose-Hubbard Hamiltonian: quantum chaos approach
• Discretising the Herman-Kluk propagator
• Spectra of harmonium in a magnetic field using an initial value representation of the semiclassical propagator

Cited in

(6)

• On the reduced dynamics of a subset of interacting bosonic particles
• The semiclassical propagator in Fock space: dynamical echo and many-body interference
• Describing many-body bosonic waveguide scattering with the truncated Wigner method
• Comment on  Semiclassical approximations in phase space with coherent states
• Symmetry-induced many-body quantum interference in chaotic bosonic systems: an augmented truncated Wigner method
• The dynamics of the Hubbard model through stochastic calculus and Girsanov transformation

This page was built for publication: Dynamics of interacting bosons using the Herman-Kluk semiclassical initial value representation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q3186631)