We present a phase-space method for the Bose-Hubbard model based on the Q-function representation. In particular, we consider two model Hamiltonians in the mean-field approximation; the first is the standard "one site" model where quantum tunneling is approximated entirely using mean-field terms; the second "two site" model explicitly includes tunneling between two adjacent sites while treating tunneling with other neighbouring sites using the mean-field approximation. The ground state is determined by minimizing the classical energy functional subject to quantum mechanical constraints, which take the form of uncertainty relations. For each model Hamiltonian we compare the ground state results from the Q-function method with the exact numerical solution. The results from the Q-function method, which are easy to compute, give a good qualitative description of the main features of the Bose-Hubbard model including the superfluid to Mott insulator. We find the quantum mechanical constraints dominate the problem and show there are some limitations of the method particularly in the weak lattice regime.
We present a phase-space method for the Bose-Hubbard model based on the Q-function representation. In particular, we consider two model Hamiltonians in the mean-field approximation; the first is the standard "one site" model where quantum tunneling is approximated entirely using mean-field terms; the second "two site" model explicitly includes tunneling between two adjacent sites while treating tunneling with other neighbouring sites using the mean-field approximation. The ground state is determined by minimizing the classical energy functional subject to quantum mechanical constraints, which take the form of uncertainty relations. For each model Hamiltonian we compare the ground state results from the Q-function method with the exact numerical solution. The results from the Q-function method, which are easy to compute, give a good qualitative description of the main features of the Bose-Hubbard model including the superfluid to Mott insulator. We find the quantum mechanical constraints dominate the problem and show there are some limitations of the method particularly in the weak lattice regime.
We present a phase-space method for the Bose–Hubbard model based on the Q-function representation. In particular, we consider two model Hamiltonians in the mean-field approximation; the first is the standard 'one-site' model where quantum tunnelling is approximated entirely using mean-field terms; the second 'two-site' model explicitly includes tunnelling between two adjacent sites while treating tunnelling with other neighbouring sites using the mean-field approximation. The ground state is determined by minimizing the classical energy functional subject to quantum mechanical constraints, which take the form of uncertainty relations. For each model Hamiltonian, we compare the ground state results from the Q-function method with the exact numerical solution. The results from the Q-function method, which are easy to compute, give a good qualitative description of the main features of the Bose–Hubbard model including the superfluid to Mott insulator. We find the quantum mechanical constraints dominate the problem and show there are some limitations of the method particularly in the weak lattice regime.
In this thesis we consider the application of phase-space methods to Bose-Einstein condensates; the work comprises of three main parts: Part I: A phase-space method for the Bose-Hubbard model; Part II: An analogue model of an expanding universe in Bose-Einstein condensates. and Part III: An analogue model of an acoustic Black Hole in Bose-Einstein condensates. In part I we present a phase-space method for the Bose-Hubbard model based on the Qfunction representation. In particular, we consider two model Hamiltonians in the mean-field approximation; the first is the standard “one site” model where quantum tunneling is approximated entirely using mean-field terms; the second “two site” model explicitly includes tunneling between two adjacent sites while treating tunneling with other neighbouring sites using the meanfield approximation. The ground state is determined by minimising the classical energy functional subject to quantum mechanical constraints, which take the form of uncertainty relations. For each model Hamiltonian we compare the ground state results from the Q-function method with the exact numerical solution. The results from the Q-function method, which are easy to compute, give a good qualitative description of the main features of the Bose-Hubbard model including the superfluid to Mott insulator. We find the quantum mechanical constraints dominate the problem and show there are some limitations of the method particularly in the weak lattice regime. Analogue models of gravity have been motivated by the possibility of investigating phenomena not readily accessible in their cosmological counterparts. In particular, the prediction of quasiparticle creation in ultra-cold Bose gases in specific configurations can be viewed as an analogue to either cosmological particle creation or the Hawking effect. In part II of this thesis we investigate the analogue of cosmological particle creation in an expanding universe by numerically simulating a Bose-Einstein condensate with a time-dependent scattering length. In particular, we simulate a 2D homogeneous condensate using the classical field method via the truncated Wigner approximation. We show that for several different expansion scenarios the calculated particle production is consistent with the underlying theory. For inflationary models we find the particle production for long wavelength modes coincides with the analytic theory within the acoustic approximation, whereas the particle production is suppressed for short wavelength (ie. free-particle like) modes. Moreover, particle production is enhanced for faster expansions, approaching the analytic result for the sudden expansion in the limit of a very fast expansion. For the case of a cyclic expansion, particle production peaks for a mode frequency that is approximately half of the driving frequency as expected for parametric resonance. In part III of this thesis we investigate an acoustic black hole in a Bose-Einstein condensate, formed by two de Laval nozzles in a ring configuration — a system we refer to as the quantum de Laval nozzle. Our model is formulated in one dimension with a sinusoidal potential. For nonzero superfluid flow, this system can exhibit stable transonic flow with both black and white hole sonic horizons. Stationary states are found by solving the time-independent Gross-Pitaevskii equation subject to a phase quantisation constraint. By solving the projected Bogoliubov-de Gennes equations for the system, we also find the discrete spectrum and quasiparticle modes. There are dynamical instabilities for certain values of winding number and potential depth, for which it is possible to construct pairs of normalisable modes. We further investigate the dynamics of the system using a classical field method based on the truncated Wigner approximation. For a low winding number and unstable configuration, we find exponential growth for the pair of unstable modes, whereas there is no growth in these modes for a stable configuration. This can be interpreted as non-degenerate parametric amplification, valid for short times. In contrast, for a large winding number, there is significant growth in modes for both stable and unstable configurations. This is indicative of higher order processes neglected in the quasiparticle picture, which is further reinforced by that fact that large winding number solutions require large nonlinearities. Finally, we consider the connection of our results with the usual semi-classical prediction of the Hawking effect. For an unstable configuration, the normalised unstable modes couple equal and opposite real frequencies, so that the growth in these modes represents the closest analogy with the Hawking effect for our quantum system.
The one-dimensional Bose-Hubbard model at a unit filling factor is studied by means of a very high-order symbolic perturbative expansion. Analytical expressions are derived for the ground-state quantities such as energy per site, variance of on-site occupation, and correlation functions: $⟨{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{a}}_{j}^{\ifmmode\dagger\else\textdagger\fi{}}{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{a}}_{j+r}⟩$ and $⟨{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{n}}_{j}{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{n}}_{j+r}⟩$. These findings are compared to numerics and good agreement is found in the Mott insulator phase. Our results provide analytical approximations to important observables in the Mott phase, and are also of direct relevance to future experiments with ultracold atomic gases placed in optical lattices. We also discuss the symmetry of the Bose-Hubbard model associated with the sign change of the tunneling coupling.
The Bose-Hubbard model describes the physics of a system of bosonic ultracold atoms in an optical lattice, in which a phase transition is present between a superfluid phase and a Mott insulator one. \nThe exact solution of this Hamiltonian is only feasible to find the ground-state of small systems, while other techniques (as mean-field schemes or quantum Monte Carlo) are necessary to study systems of larger size. \n \nAs a first application, we study the trapped model - relevant for the comparison with current experiments - through an inhomogeneous mean-field scheme. We describe some signatures of the phase crossover between superfluid and Mott insulator. In particular, the visibility of the quasimomentum distribution shows some kinks as a function of the lattice depth; we describe these features and we link them with the ones observed in other works in the literature. \n \nAs a second application, we use quantum Monte Carlo techniques to study the one-dimensional Bose-Hubbard model with long-range interactions and we focus on the appearance of the Haldane insulating phase, distinguishable from the Mott one through the presence of non-local hidden order. \nNon-local correlation functions are also used to describe the difference between the superfluid phase and the Mott insulator one.
The exact solution of a quantum Bethe lattice model in the thermodynamic limit amounts to solve a functional self-consistent equation. In this paper we obtain this equation for the Bose-Hubbard model on the Bethe lattice, under two equivalent forms. The first one, based on a coherent-state path integral, leads in the large connectivity limit to the mean-field treatment of Fisher et al. [Phys. Rev. B 40, 546 (1989)] at the leading order, and to the bosonic dynamical mean field theory as a first correction, as recently derived by Byczuk and Vollhardt [Phys. Rev. B 77, 235106 (2008)]. We obtain an alternative form of the equation using the occupation number representation, which can be easily solved with an arbitrary numerical precision, for any finite connectivity. We thus compute the transition line between the superfluid and Mott insulator phases of the model, along with thermodynamic observables and the space and imaginary-time dependence of correlation functions. The finite connectivity of the Bethe lattice induces a richer physical content with respect to its infinitely connected counterpart: a notion of distance between sites of the lattice is preserved, and the bosons are still weakly mobile in the Mott insulator phase. The Bethe lattice construction can be viewed as an approximation to the finite-dimensional version of the model. We show indeed a quantitatively reasonable agreement between our predictions and the results of Quantum Monte Carlo simulations in two and three dimensions.
Abstract This paper is a first step toward rigorously establishing the validity of mean-field approximations for the dynamics of quantum systems in a large dimensional limit. Our result is a derivation of the large d mean-field limit of the Bose–Hubbard model on a square lattice. The Bose–Hubbard model effectively describes bosons on a lattice with on-site interactions and nearest-neighbour hopping, serving as a foundational framework for understanding strong particle interactions and the superfluid to Mott insulator phase transition. We prove a trace norm estimate between the one-lattice-site reduced density of the Schrödinger dynamics and the mean-field dynamics in the limit of large dimension. Here, the mean-field approximation is in the hopping amplitude and not in the interaction, leading to a very rich and non-trivial mean-field equation. This mean-field equation does not only describe the condensate, as is the case when the mean-field description comes from a large particle number limit averaging out the interaction, but it allows for a phase transition since it contains the full non-trivial interaction. Our work is a rigorous justification of a simple case of the highly successful dynamical mean-field theory (DMFT) for bosons, which somewhat surprisingly yields many qualitatively correct results in three dimensions.