N-boson time-dependent problem: A reformulation with stochastic wave functions
Abstract
We present a numerically tractable method to solve exactly the evolution of a N boson system with binary interactions. The density operator of the system ρ is obtained as the stochastic average of particular operators \|Ψ1><Ψ2\| of the system. The states \|Ψ1,2> are either Fock states \|N:φ1,2> or coherent states \|coh:φ1,2> with each particle in the state φ1,2(x). We determine the conditions on the evolution of φ1,2, which involves a stochastic element, under which we recover the exact evolution of ρ. We discuss various possible implementations of these conditions. The well known positive P representation arises as a particular case of the coherent state ansatz. We treat numerically two examples: a two-mode system and a one-dimensional harmonically confined gas. These examples, together with an analytical estimate of the noise, show that the Fock state ansatz is the most promising one in terms of precision and stability of the numerical solution.
- Publication:
-
Physical Review A
- Pub Date:
- February 2001
- DOI:
- 10.1103/PhysRevA.63.023606
- arXiv:
- arXiv:cond-mat/0003399
- Bibcode:
- 2001PhRvA..63b3606C
- Keywords:
-
- 03.75.Fi;
- 05.10.Gg;
- 42.50.-p;
- Stochastic analysis methods;
- Quantum optics;
- Condensed Matter - Soft Condensed Matter
- E-Print:
- 21 pages, 5 figures, submitted to Phys.Rev. A