Semiclassical analysis of Bose-Hubbard dynamics

In this work the two-site Bose-Hubbard model is studied analytically in the limit of weak coupling u and large number of particles N. In particular, the difference in the occupation between the two sites, where initially all particles are at one site, was calculated analytically. This quantity exhibits collapses and revivals that superimpose rapid oscillations. Excellent agreement with the exact numerical solution was found. The semiclassical approximation where 1/N plays the role of Planck's constant was used and perturbation theory to order u² was applied. The occupation difference was calculated also for the case where initially both sites are occupied provided that the difference in occupation is sufficiently large. This work provides an analytical description of results that were so far found only numerically. Similar behavior and analysis are expected for a large variety of physical situations.