Bose–Einstein condensate in a double-well potential: Feynman path integral variational approach

A Bose–Einstein condensate in a double-well potential is considered by using the Feynman path integral theory combined with a variational approach. The system consists of N interacting bosons confined in the double well which is taken as a harmonic potential with a Gaussian barrier. The calculation is carried out within the first cumulant approximation measured with respect to the harmonic action containing variational parameters. Assuming a separable expression for the trial action corresponds to the usual mean field approximation. Performing the variational calculations, we obtain analytical results for the ground state energy and its condensate wavefunction. Using a projection onto the odd contributions, the first excited state energy and condensate wavefunction are determined, too. We find very good agreement with a full numerical solution of the Gross–Pitaevskii equation over the full range of potential parameters. (paper)