Reconciling the Classical-Field Method with the Beliaev Broken Symmetry Approach

We present our views on the issues raised in the chapter by Griffin and Zaremba [A. Griffin and E. Zaremba, in Quantum Gases: Finite Temperature and Non-Equilibrium Dynamics, N. P. Proukakis, S. A. Gardiner, M. J. Davis, and M. H. Szymanska, eds., Imperial College Press, London (in press)]. We review some of the strengths and limitations of the Bose symmetry-breaking assumption, and explain how such an approach precludes the description of many important phenomena in degenerate Bose gases. We discuss the theoretical justification for the classical-field (c-field) methods, their relation to other non-perturbative methods for similar systems, and their utility in the description of beyond-mean-field physics. Although it is true that present implementations of c-field methods cannot accurately describe certain collective oscillations of the partially condensed Bose gas, there is no fundamental reason why these methods cannot be extended to treat such scenarios. By contrast, many regimes of non-equilibrium dynamics that can be described with c-field methods are beyond the reach of generalised mean-field kinetic approaches based on symmetry-breaking, such as the ZNG formalism.

Excitation spectrum of bosons in a finite one-dimensional circular waveguide via the Bethe ansatz

The exactly solvable Lieb-Liniger model of interacting bosons in one-dimension has attracted renewed interest as current experiments with ultra-cold atoms begin to probe this regime. Here we numerically solve the equations arising from the Bethe ansatz solution for the exact many-body wave function in a finite-size system of up to twenty particles for attractive interactions. We discuss the novel features of the solutions, and how they deviate from the well-known string solutions [H. B. Thacker, Rev. Mod. Phys. textbf{53}, 253 (1981)] at finite densities. We present excited state string solutions in the limit of strong interactions and discuss their physical interpretation, as well as the characteristics of the quantum phase transition that occurs as a function of interaction strength in the mean-field limit. Finally we compare our results to those of exact diagonalization of the many-body Hamiltonian in a truncated basis. We also present excited state solutions and the excitation spectrum for the repulsive 1D Bose gas on a ring.