Granulation of quantum matter -- the formation of persistent small-scale patterns -- is realized in the images of quasi-one-dimensional Bose-Einstein condensates perturbed by a periodically modulated interaction. Our present analysis of a mean-field approximation suggests that granulation is caused by the gradual transformation of phase undulations into density undulations. This is achieved by a suitably large modulation frequency, while for low enough frequencies the system exhibits a quasi-adiabatic regime. We show that the persistence of granulation is a result of the irregular evolution of the phase of the wavefunction representing an irreversible process. Our model predictions agree with numerical solutions of the Schrödinger equation and experimental observations. The numerical computations reveal the emergent many-body correlations behind these phenomena via the multi-configurational time-dependent Hartree theory for bosons (MCTDHB).
Granulation of quantum matter -- the formation of persistent small-scale patterns -- is realized in the images of quasi-one-dimensional Bose-Einstein condensates perturbed by a periodically modulated interaction. Our present analysis of a mean-field approximation suggests that granulation is caused by the gradual transformation of phase undulations into density undulations. This is achieved by a suitably large modulation frequency, while for low enough frequencies the system exhibits a quasi-adiabatic regime. We show that the persistence of granulation is a result of the irregular evolution of the phase of the wavefunction representing an irreversible process. Our model predictions agree with numerical solutions of the Schrödinger equation and experimental observations. The numerical computations reveal the emergent many-body correlations behind these phenomena via the multi-configurational time-dependent Hartree theory for bosons (MCTDHB).
We explore, both experimentally and theoretically, the response of an elongated Bose-Einstein condensate to modulated interactions. We identify two distinct regimes differing in modulation frequency and modulation strength. Longitudinal surface waves are generated either resonantly or parametrically for modulation frequencies near the radial trap frequency or twice the trap frequency, respectively. The dispersion of these waves, the latter being a Faraday wave, is well-reproduced by a mean-field theory that accounts for the 3D nature of the elongated condensate. In contrast, in the regime of lower modulation frequencies we find that no clear resonances occur, but with increased modulation strength, the condensate forms an irregular granulated distribution that is outside the scope of a mean-field approach. We find that the granulated condensate is characterized by large quantum fluctuations and correlations, which are well-described with single-shot simulations obtained from wavefunctions computed by a beyond mean-field theory at zero temperature, the multiconfigurational time-dependent Hartree for bosons method.
A mechanism leading to spontaneous breaking of the space symmetry, and to the appearance of spatial patterns in the most recently created state of matter, the Bose-Einstein condensate (BEC), is proposed. The pattern formation are due to a periodic modulation in time of the scattering length of the condensed atoms. Experimentally this can be achieved via the Feshbach resonance, or, for low dimensional (1D or 2D) condensates, also by modulation of the trapping potential along a spatial coordinate normal to the "working" space. The parametric excitation of patterns in condensates by considering the mean-field Gross-Pitaevskii equation, and its dissipative generalization are investigated.
The mean-field dynamics of Bose-Einstein condensates loaded in an optical lattice, confined by a parabolic potential, and subjected to a change of scattering length by means of the Feshbach resonance is considered. The system is described by the Gross-Pitaevskii (GP) equation with varying nonlinearity, which in a number of cases is reduced to one-dimensional (1D) perturbed nonlinear Schr\"odinger (NLS) equations, the particular form of which depends on the relation among the parameters of the problem. We analytically describe the adiabatic dynamics of periodic solutions of the respective NLS equations, provide a numerical study of 1D models confined by a parabolic trap, and carry out numerical simulations of the matter-wave dynamics within the framework of the radially symmetric 3D GP equation. Special attention is paid to processes of generation of trains of bright and dark matter solitons from initially periodic waves. The results of the 1D approximation are compared with direct numerical simulation of the original 3D GP eqation, showing remarkable coincidence for definite regions of parameters.
An experimentally realizable scheme of periodic sign-changing modulation of the scattering length is proposed for Bose-Einstein condensates similar to dispersion-management schemes in fiber optics. Because of controlling the scattering length via the Feshbach resonance, the scheme is named Feshbach-resonance management. The modulational-instability analysis of the quasiuniform condensate driven by this scheme leads to an analog of the Kronig-Penney model. The ensuing stable localized structures are found. These include breathers, which oscillate between the Thomas-Fermi and Gaussian configuration, or may be similar to the 2-soliton state of the nonlinear Schrödinger equation, and a nearly static state ("odd soliton") with a nested dark soliton. An overall phase diagram for breathers is constructed, and full stability of the odd solitons is numerically established.
We present semianalytical many-body results for energies and excitation frequencies for an inhomogeneous Bose-Einstein condensate over a wide range of atom numbers N for both small s-wave scattering lengths, typical of most laboratory experiments, and large scattering lengths, achieved by tuning through a Feshbach resonance. Our dimensional perturbation treatment includes two-body correlations at all orders and yields analytical results through first order by taking advantage of the high degree of symmetry of the condensate at the zeroth-order limit. Because N remains a parameter in our analytical results, the challenge of calculating energies and excitation frequencies does not rise with the number of condensate atoms. In this proof-of-concept paper the atoms are confined in a spherical trap and are treated as hard spheres. Our many-body calculations compare well to Gross-Pitaevskii results in the weakly interacting regime and depart from the mean-field approximation as the density approaches the strongly interacting regime. The excitation frequencies provide a particularly sensitive test of beyond-mean-field corrections. For example, for N=2000 atoms and an experimentally realized large scattering length of $a=0.433{a}_{ho}\phantom{\rule{0.3em}{0ex}}({a}_{ho}=\sqrt{\ensuremath{\hbar}∕m{\ensuremath{\omega}}_{ho}})$ we predict a 75 shift from the mean-field breathing mode frequency.
When a gas of bosonic atoms is cooled sufficiently near to absolute zero, it undergoes a startling transition from a collection of randomly-behaved individual atoms to a phase in which thousands, even millions, of atoms act in collective unison. In this new phase, known as a Bose-Einstein condensate (BEC), the atoms all occupy the exact same quantum mechanical ground state and can be collectively described by a single phase-coherent wavefunction of macroscopic size. But how exactly does this harmonious behaviour arise within a collection of incoherent atoms as their temperature decreases? Does the condensate wavefunction grow from a single central nucleus, or does it emerge simultaneously throughout the system? And how does the situation differ when the gas is cooled suddenly, away from thermodynamic equilibrium? Can non-equilibrium behaviour be inherited by the condensate that forms? In this thesis we aim to answer these and other questions by studying the formation dynamics of Bose-Einstein condensates and the development of phase coherence. The research described herein makes use of several theoretical tools, which we lay out on the workbench in Chapter 2: Semi-classical Hartree-Fock theory provides a mean-field description of a condensate and thermal cloud at equilibrium; Quantum kinetic theory is used to calculate the dynamics of condensate growth into an equilibrium state, from a non-equilibrium thermal cloud; The stochastic projected Gross-Pitaevskii equation (SPGPE) employs a classical-field description to encapsulate the non-equilibrium dynamics of the condensate and low-energy region of the atom cloud, coupled to an equilibrium thermal cloud. This thesis is composed of four principal studies, each constituting a separate chapter.In Chapter 3 we study the formation of a Bose-Einstein condensate in a cigar-shaped 3D harmonic trap, induced by the controlled addition of an attractive ``dimple'' potential along the weak axis in a recent experiment [Garrett et al., Phys. Rev. A 83, 013630 (2011)]. In this manner condensation is induced without cooling due to a localized increase in the phase space density. We perform a quantitative analysis of the thermodynamic transformation in both the sudden and adiabatic regimes for a range of dimple widths and depths. We find good agreement with equilibrium calculations based on self-consistent semi-classical Hartree-Fock theory describing the condensate and thermal cloud. We observe there is an optimal dimple depth that results in a maximum in the condensate fraction. We also study the non-equilibrium dynamics of condensate formation in the sudden turn-on regime, finding good agreement for the observed time dependence of the condensate fraction with calculations based on quantum kinetic theory. In Chapter 4 we study the onset of phase coherence during the formation of a BEC following a sudden quench to below the transition temperature. Using the SPGPE, we simulate a recent experiment performed at ANU in which a sudden RF quench is applied to a dilute gas of metastable Helium atoms confined in an elongated 3D harmonic magnetic trap. At regularly-spaced intervals following the quench, broadband RF pulses are applied to homogeneously outcouple small (<5%) pulses of atoms, which expand and fall to a microchannel plate detector. The second-order correlation function at the detector, postulated to provide a measure of phase coherence of the in-situ cloud from which they were outcoupled, is measured versus time and compared with the condensate growth curve. A delay was observed in the relaxation of the second-order correlation function, relative to the peak in the number of condensate atoms. In our simulations we incorporate the effects of periodic outcoupling both the condensate and thermal cloud, and monitor the time-dependence of first- and second-order correlation functions in both position and momentum space. We observe long-range phase coherence in the final equilibrium state, though mild persistent phase fluctuations are present. We also observe the curious phenomenon of ``superbunching'' - an enhancement of density fluctuations - at specific values of the momentum, which we relate to anomalous correlations of fluctuations. Despite obtaining a good quantitative fit to the experimentally-measured condensate growth curve in our SPGPE simulations, we do not observe any evidence of the delay in coherence observed in the experiment, which remains unexplained. In Chapter 5 we use the SPGPE to study the crossover from a true BEC with long-range phase coherence to a quasicondensate dominated by phase fluctuations, for the case of a Bose gas in an highly-elongated 3D harmonic trap. We use parameters corresponding to an experiment in which a sudden quench resulted in the formation of a non-equilibrium condensate, the dynamics of which we study in the chapter that follows. Sweeping through a range of temperatures and chemical potentials within the crossover region, we reveal the behaviour of phase fluctuations through analysis of correlation functions. We compare the ideal-gas approximation for the transition temperature with the more accurate value obtained via the Binder cumulant. By comparing the phase coherence length with the spatial extent across which density fluctuations are suppressed, we identify the BEC/quasicondensate crossover point. We compare directly with predictions of mean-field theory [Petrov et al., Phys. Rev. Lett. 87, 050404 (2001)] and find that phase fluctuations persist to lower temperatures in our SPGPE simulations. We calculate anomalous correlations of the fluctuation field and show that the resulting correction to the mean-field condensate self-interaction energy is significant within the quasicondensate regime. In Chapter 6 we study the formation of a BEC into a non-equilibrium state by simulating an experiment in which a dilute gas of Rubidium-87 atoms in a highly-elongated 3D harmonic magnetic trap is rapidly quenched to below the transition temperature. Following the quench, a condensate is observed to form with an excessively elongated shape, and subsequently undergo damped quadrupole shape oscillations [Shvarchuck et al., Phys. Rev. Lett. 89, 270404 (2002)]. We simulate this experiment by modifying the SPGPE to incorporate a hydrodynamic thermal cloud, parametrized by axially and temporally-dependent temperature and chemical potential. Comparison against simulations without this axial dependence suggests the hydrodynamic behaviour of the thermal cloud plays a crucial role in the formation of a non-equilibrium condensate. Analyses of the first-order correlation function and the condensate wavefunction reveal the spontaneous formation of quantum vortices, as well as strong phase fluctuations at equilibrium indicating quasicondensate behaviour.