Search

We investigate beyond-mean-field corrections to the energy of an elongated homogeneous Bose gas strongly confined in two directions, with dipoles aligned along the long axis of the system. When the dipolar interaction reaches its critical strength, the mean-field approach predicts instability. However, similar to the free-space case, beyond-mean-field effects significantly alter the ground state of the system, leading to the formation of a self-bound atomic cloud known as a quantum droplet. Our analysis demonstrates that the beyond-mean-field contribution to the energy in the quasi-1D region, in addition to the confinement induced shift of the mean field energy, is proportional to the third power of the density $\sim n^3$. Therefore, it can be interpreted as an effective three-body repulsion that stabilizes the gas, preventing collapse and leading to a finite-density solution. We also show that the same effect plays a crucial role in the binding of strongly elongated dipolar droplets under harmonic confinement.

We study a mixture of a repulsive dipolar condensate and a degenerate Fermi gas in a quasi-one-dimensional geometry. We demonstrate that the presence of fermions, which attract bosons, drastically changes the behavior of the dipolar condensate. For strong enough boson-fermion attraction, a dipolar Bose-Fermi droplet appears in the mixture, and as the attraction becomes stronger, a roton excitation develops in the Bogoliubov spectrum, leading to the formation of a supersolid, and eventually a crystal of isolated droplets. We describe the system by coupled extended Gross-Pitaevskii (bosons) and Hartree-Fock (fermions) equations. We study the excitation spectrum of the system and identify a number of Goldstone and Higgs modes in the supersolid regime., Comment: 10 pages, 11 figures

The spin squeezing protocols allow the dynamical generation of massively correlated quantum many-body states, which can be utilized in entanglement-enhanced metrology and technologies. We study a quantum simulator generating twisting dynamics realized in a two-component Bose-Hubbard model with dipolar interactions. We show that the interplay of contact and long-range dipolar interactions between atoms in the superfluid phase activates the anisotropic two-axis counter-twisting mechanism, accelerating the spin squeezing dynamics and allowing the Heisenberg-limited accuracy in spectroscopic measurements.

We describe a quantum droplet of a Bose-Bose mixture squeezed by an external harmonic forces in one spatial direction. Our approach is based on the self-consistent method formulated in [1]. The true spatial droplet profile in the direction of confinement is accounted for, however local density approximation is assumed in the free directions. We define a numerical approach to find the beyond-mean-field contribution to the chemical potential (Lee-Huang-Yang chemical potential) -- the quantity that determines the droplet's profile. In addition to the numerical approach, we find the Lee-Huang-Yang potential in the analytic form in two limiting cases: a perturbative result for a strong confinement and a semiclassical expression when confinement is very weak.

We define a formalism of a self-consistent description of the ground state of a weakly interacting Bose system, accounting for higher order terms in expansion of energy in the diluteness parameter. The approach is designed to be applied to a Bose-Bose mixture in a regime of weak collapse where quantum fluctuations lead to stabilization of the system and formation of quantum liquid droplets. The approach is based on the Generalized Gross -- Pitaevskii equation accounting for quantum depletion and anomalous density terms. The equation is self-consistently coupled to modified Bogoliubov equations. The modification we introduce resolves the longstanding issue of missing phonon-branch excitations when higher order terms are included. Our method ensures a gapless phononic low-energy excitation spectrum, crucial to correctly account for quantum fluctuations. We pay particular attention to the case of droplets harmonically confined in some directions. The method allows to determine the Lee-Huang-Yang-type contribution to the chemical potential of inhomogeneous droplets when the local density approximation fails.

We investigate a Bose-Einstein condensate of $F= 1$ $^{87}$Rb atoms in a 2D spin-dependent optical lattice generated by intersecting laser beams with a superposition of polarizations. For $^{87}$Rb the effective interaction of an atom with the electromagnetic field contains a scalar and a vector (called as fictitious magnetic field, $B_{fic}$) potentials. The Rb atoms behave as a quantum rotor (QR) with angular momentum given by the sum of the atomic rotational motion angular momentum and the hyperfine spin. The ground state of the QR is affected upon applying an external magnetic field, $B_{ext}$, perpendicular to the plane of QR motion and a sudden change of its topology occurs as the ratio $B_{ext}/B_{fic}$ exceeds critical value. It is shown that the change of topology of the QR ground state is a result of combined action of Zeeman and Einstein-de Haas effects. The first transfers atoms to the largest hyperfine component to polarize the sample along the field as the external magnetic field is increased. The second sweeps spin to rotational angular momentum, modifying the kinetic energy of the atoms., Comment: 5 pages + Supplemental material

We study coherent dynamics of two interacting Bose-Bose droplets by means of the extended Gross-Pitaevskii equation. The relative motion of the droplets couples to the phases of their components. The dynamics can be understood in terms of the evolution of zero-energy modes recovering symmetries spontaneously broken by the mean-field solution. These are translational symmetry and two U(1) symmetries, associated with the phases of the droplets' two components. A phase-dependent interaction potential and double Josephson-junction equations are introduced to explain the observed variety of different scenarios of collision. We show that the evolution of the droplets is a macroscopic manifestation of the hidden dynamics of their phases. The occurrence of nondissipative drag between the two supercurrents (Andreev-Bashkin effect) is mentioned.

Bose-Bose droplets are self-bound objects emerging from a mixture of two interacting Bose-Einstein condensates when their interactions are appropriately tuned. During droplet formation three continuous symmetries of the system's Hamiltonian are broken: translational symmetry and two U1 symmetries, allowing for arbitrary choice of phases of the mean-field wavefunctions describing the two components. Breaking of these symmetries must be accompanied by appearance of zero-energy excitations in the energy spectrum of the system recovering the broken symmetries. Normal modes corresponding to these excitations are the zero-energy modes. Here we find analytic expressions for these modes and introduce Hamitonians generating their time evolution -- dynamics of the droplet's centers of mass as well as dynamics of the phases of the two droplet's wavefunctions. When internal types of excitations (quasiparticles) are neglected then the very complex system of a quantum droplet is described using only few "global" degrees of freedom - the position of the center of mass of the droplet and two phases of two wave-functions, all these being quantum operators. This gives the possibility of describing in a relatively easy way processes of interaction of these quantum droplets, such as collisions., Comment: 17 pages, no figures, Category: cond-mat.quant-gas

We study the spin distillation of spinor gases of bosonic atoms and find two different mechanisms in ${}^{52}$Cr and $^{23}$Na atoms, both of which can cool effectively. The first mechanism involves dipolar scattering into initially unoccupied spin states and cools only above a threshold magnetic field. The second proceeds via equilibrium relaxation of the thermal cloud into empty spin states, reducing its proportion in the initial component. It cools only below a threshold magnetic field. The technique was initially demonstrated experimentally for a chromium dipolar gas [B. Naylor et al., Phys. Rev. Lett. 115, 243002 (2015)], whereas here we develop the concept further and provide an in-depth understanding of the required physics and limitations involved. Through numerical simulations, we reveal the mechanisms involved and demonstrate that the spin distillation cycle can be repeated several times, each time resulting in a significant additional reduction of the thermal atom fraction. Threshold values of magnetic field and predictions for the achievable temperature are also identified., Comment: 15 pages, 11 figures. Extended discussion and improved sectioning

Recently observed Pauli crystals are structures formed by trapped ultracold atoms with the Fermi statistics. Interactions between these atoms are switched off, so their relative positions are determined by joined action of the trapping potential and the Pauli exclusion principle. Numericalmodeling is used in this paper to find the Pauli crystals in a two-dimensional isotropic harmonic trap, three-dimensional harmonic trap, and a two-dimensional square well trap. The Pauli crystals do not have the symmetry of the trap -- the symmetry is broken by the measurement of positions and, in many cases, by the quantum state of atoms in the trap. Furthermore, the Pauli crystals are compared with the Coulomb crystals formed by electrically charged trapped particles. The~structure of the Pauli crystals differs from that of the Coulomb crystals, this provides evidence that the exclusion principle cannot be replaced by a two-body repulsive interaction but rather has to be considered to be a specifically quantum mechanism leading to many-particle correlations., Comment: This article belongs to the Special Issue "Symmetries and the Pauli Exclusion Principle"

We study the problem of the stability of a two-component droplet. The standard solution known from the literature is based on a particular form of the mean field energy functional, in particular on distinction of hard mode and soft mode contributions. By imposing the constraint on densities of the two species which minimizes the hard mode energy, the problem is reduced to a stability analysis of a one component system. As opposed to this, we address the issue in full generality. Our analysis is valid for arbitrary forms of energy density. We formulate constraints which correspond to the physically relevant situation of a system which has unconstrained volume and may evaporate particles. For the specific case of a two component Bose-Bose droplet we find approximate analytic solutions and compare them to the standard result. We show that the densities of both components of a stable droplet are limited to a range depending on interaction strength, in contrast to the original unique solution.

We study the properties of the Bose-Fermi mixture from the perspective of reaching a state of a self-bound quantum droplet. The variational analysis shows that the system exhibits bistability. For weak repulsion between bosons, one of the equilibrium states, smaller in size, spherically symmetric, and with negative energy, corresponds to quantum droplet, the other with always positive energy represents the elongated droplet-like state immersed in the sea of a fermionic cloud. For stronger repulsion between bosons the bifurcation is seized and only the former state is left. Now it represents an elongated object which, for strong enough boson-fermion attraction, gets negative energy. It becomes an excited Bose-Fermi droplet when the trap is released, what is demonstrated by solving the quantum hydrodynamics equations for the Bose-Fermi system. To depict our ideas we consider the $^{133}$Cs-$^6$Li mixture under ideal conditions, i.e. we assume no losses., Comment: 6 pages, 5 figures

We study the final stages of the evolution of a binary system consisted of a black hole and a white dwarf star. We implement the quantum hydrodynamic equations and carry out numerical simulations. As a model of a white dwarf star, we consider a zero temperature droplet of attractively interacting degenerate atomic bosons and spin-polarized atomic fermions. Such mixtures are investigated experimentally nowadays. We find that the white dwarf star is stripped off its mass while passing the periastron. Due to nonlinear effects, the accretion disk originated from the white dwarf becomes fragmented and the onset of a quantum turbulence with giant quantized vortices present in the bosonic component of the accretion disk is observed. The binary system ends its life in a spectacular way, revealing quantum features underlying the white dwarf star's structure. We find a charged mass, falling onto a black hole, could be responsible for recently discovered ultraluminous X-ray bursts. The simulations show that the final passage of a white dwarf near a black hole can cause a gamma-ray burst., Comment: 14 pages, 7 figures. Corresponding author -> MN email: m.nikolajuk@uwb.edu.pl; Printed in Scientific Reports

We study weakly interacting mixtures of ultracold atoms composed of bosonic and fermionic species in 2D and 1D. When interactions between particles are appropriately tuned, self-bound quantum liquids can be formed. We show that while formation of these droplets in 2D is due to the higher order correction terms contributing to the total energy and originating in quantum fluctuations, in 1D geometry the quantum fluctuations have a negligible role on formation of the self-bound systems. The leading mean-field interactions are then sufficient for droplet formation in 1D. We analyse stability conditions for 2D and 1D systems and predict values of equilibrium densities of droplets., Comment: Modified analysis in 2D and new results in 1D; 5 pages, 6 figures

We study a liquid quantum droplets in a mixture of two-component Bose-Einstein condensates under a variable confinement introduced along one or two spatial dimensions. Despite the atom-atom scattering has a three-dimensional character, discreetness of the available modes in the reduced dimension(s) strongly influences the zero-point energy -- the Lee-Huang-Yang term. In a weakly interaction limit, it is the leading correction to the mean-field energy at the crossover from three to two dimensions, or from three to one dimension. We analyze the properties of the droplets at the dimensional crossovers, and provide the demanding conditions for accessing quasi-low dimensions. We predict new kinds of droplets which are formed only due to the quantum fluctuations when the mean-field interaction vanishes. Our results pave the way for exploring new states of quantum matter, and are important for experiments with liquid quantum droplets in reduced dimensions., Comment: 5 pages + supplemental material

We study the stability of a zero temperature mixture of attractively interacting degenerate bosons and spin-polarized fermions in the absence of confinement. We demonstrate that higher order corrections to the standard mean-field energy can lead to a formation of Bose-Fermi liquid droplets -- self-bound systems in three-dimensional space. The stability analysis of the homogeneous case is supported by numerical simulations of finite systems by explicit inclusion of surface effects. We discuss the experimental feasibility of formation of quantum droplets and indicate the main obstacle -- inelastic three-body collisions., Comment: Resubmission to SciPost

The best known manifestation of the Fermi-Dirac statistics is the Pauli exclusion principle: no two identical fermions can occupy the same one-particle state. This principle enforces high order correlations in systems of many identical fermions and is responsible for a particular geometric arrangement of trapped particles even when all mutual interactions are absent [1]. These geometric structures, called Pauli crystals, are predicted for a system of $N$ identical atoms trapped in a harmonic potential. They emerge as the most frequent configurations in a collection of single-shot pictures of the system. Here we study how fragile Pauli crystals are when realistic experimental limitations are taken into account. The influence of the number of single-shots pictures available to analysis, thermal fluctuations and finite efficiency of detection are considered. The role of these sources of noise on the possibility of experimental observation of Pauli crystals is shown and conditions necessary for the detection of the geometrical arrangements of particles are identified., Comment: 11 pages, 6 figures

System of a two-flavor mixture of ultra-cold fermions confined in a one-dimensional harmonic trap is studied in the frame of the center of mass. We present a numerical method of obtaining energetic spectra in this frame for an arbitrary mass ratio of fermionic species. We identify a specific invariant encoded in many-body correlations which enable one to determine an eigenstate of the Hamiltonian and to label excitations of the center of mass. The tool presented may be particularly useful in experimental analysis of the interparticle interactions which do not affect the center of mass excitations in a harmonic potential.

We study the impact of spin-exchange collisions on the dynamics of Bose-Einstein condensation, by rapidly cooling a chromium multi-component Bose gas. Despite relatively strong spin-dependent interactions, the critical temperature for Bose-Einstein condensation is reached before the spin-degrees of freedom fully thermalize. The increase in density due to Bose-Einstein condensation then triggers spin dynamics, hampering the formation of condensates in spin excited states. Small metastable spinor condensates are nevertheless produced, and manifest strong spin fluctuations., Comment: 5 pages, 4 figures

Eigenstates of Bose particles with repulsive contact interactions in one-dimensional space with periodic boundary conditions can be found with the help of the Bethe ansatz. The type~II excitation spectrum identified by E. H. Lieb, reproduces the dispersion relation of dark solitons in the mean-field approach. The corresponding eigenstates possess translational symmetry which can be broken in measurements of positions of particles. We analyze emergence of single and double solitons in the course of the measurements and investigate dynamics of the system. In the weak interaction limit, the system follows the mean-field prediction for a short period of time. Long time evolution reveals many-body effects that are related to an increasing uncertainty of soliton positions. In the strong interaction regime particles behave like impenetrable bosons. Then, the probability densities in the configuration space become identical to the probabilities of non-interacting fermions but the wave-functions themselves remember the original Bose statistics. Especially, the phase flips that are key signatures of the solitons in the weak interaction limit, can be observed in the time evolution of the strongly interacting bosons., Comment: 11 pages, 9 figures

Dynamical properties of a few ultra-cold fermions confined in a double-well potential is studied. We show that the dynamics, which is governed by single-particle tunnelings for vanishing interactions, is completely different for strong interactions. Depending on the details of the configuration, for sufficiently strong interactions (repulsions or attractions) the particle flow through the barrier can be accelerated or slowed down. This effect cannot be explained with the single-particle picture. It is clarified with direct inspection to the spectrum of the few-body Hamiltonian., Comment: Authors' version accepted in EPL

Recently developed techniques allow for simultaneous measurements of the positions of all ultra cold atoms in a trap with high resolution. Each such single shot experiment detects one element of the quantum ensemble formed by the cloud of atoms. Repeated single shot measurements can be used to determine all correlations between particle positions as opposed to standard measurements that determine particle density or two-particle correlations only. In this paper we discuss the possible outcomes of such single shot measurements in case of cloud of ultra-cold non-interacting Fermi atoms. We show that the Pauli exclusion principle alone leads to correlations between particle positions that originate from unexpected spatial structures formed by the atoms.

A system of two species of fermions of different mass confined in a one-dimensional harmonic trap is studied with an exact diagonalization approach. It is shown independently on the number of particles that a mass difference between fermionic species induces a separation in the lighter flavor system. The mechanism of emerging of separated phases is explained phenomenologically and confirmed with the help of a direct inspection of the ground-state of the system. Finally, it is shown that the separation driven by a mass difference, in contrast to the separation induced by a difference of populations, is robust to the interactions with thermal environment.

In this work we present a very simple and efficient numerical scheme which can be applied to study the dynamics of bosonic systems like, for instance, spinor Bose-Einstein condensates with nonlocal interactions but equally well works for Fermi gases. The method we use is a modification of well known Split Operator Method (SOM). We carefully examine this algorithm in the case of $F=1$ spinor Bose-Einstein condensate without and with dipolar interactions and for strongly interacting two-component Fermi gas. Our extension of the SOM method has many advantages: it is fast, stable, and keeps constant all the physical constraints (constants of motion) at high level., Comment: 21 pages, 7 figures

We study, within a framework of the classical fields approximation, the density correlations of a weakly interacting expanding Bose gas for the whole range of temperatures across the Bose-Einstein condensation threshold. We focus on elongated quasi-one-dimensional systems where there is a huge discrepancy between the existing theory and experimental results (A. Perrin et al., Nature Phys. 8, 195 (2012)). We find that the density correlation function is not reduced for temperatures below the critical one as it is predicted for the ideal gas or for a weakly interacting system within the Bogoliubov approximation. This behavior of the density correlations agrees with the above mentioned experiment with the elongated system. Although the system was much larger then studied here we believe that the behavior of the density correlation function found there is quite generic. Our theoretical studies indicate also large density fluctuations in the trap in the quasicondensate regime where only phase fluctuations were expected. We argue that the enhanced density fluctuations can originate in the presence of interactions in the system, or more precisely in the existence of spontaneous dark solitons in the elongated gas at thermal equilibrium., Comment: Submited to PRA

We study a strongly attractive system of a few spin-1/2 fermions confined in a one-dimensional harmonic trap, interacting via two-body contact potential. Performing exact diagonalization of the Hamiltonian we analyze the ground state and the thermal state of the system in terms of one-- and two--particle reduced density matrices. We show how for strong attraction the correlated pairs emerge in the system. We find that the fraction of correlated pairs depends on temperature and we show that this dependence has universal properties analogous to the gap function known from the theory of superconductivity. In contrast to the standard approach based on the variational ansatz and/or perturbation theory, our predictions are exact and are valid also in a strong attraction limit. Our findings contribute to the understanding of strongly correlated few-body systems and can be verified in current experiments on ultra-cold atoms., Comment: 6 figures

Originally, the Hubbard model has been derived for describing the behaviour of strongly-correlated electrons in solids. However, since over a decade now, variations of it are also routinely being implemented with ultracold atoms in optical lattices. We review some of the rich literature on this subject, with a focus on more recent non-standard forms of the Hubbard model. After an introduction to standard (fermionic and bosonic) Hubbard models, we discuss briefly common models for mixtures, as well as the so called extended Bose-Hubbard models, that include interactions between neighboring sites, next-neighboring sites, and so on. The main part of the review discusses the importance of additional terms appearing when refining the tight-binding approximation on the original physical Hamiltonian. Even when restricting the models to the lowest Bloch band is justified, the standard approach neglects the density-induced tunneling (which has the same origin as the usual on-site interaction). The importance of these contributions is discussed for both contact and dipolar interactions. For sufficiently strong interactions, also the effects related to higher Bloch bands become important even for deep optical lattices. Different approaches that aim at incorporating these effects, mainly via dressing the basis Wannier functions with interactions, leading to effective, density-dependent Hubbard-type models, are reviewed. We discuss also examples of Hubbard-like models that explicitly involve higher $p$-orbitals, as well as models that couple dynamically spin and orbital degrees of freedom. Finally, we review mean-field nonlinear-Schr\"odinger models of the Salerno type that share with the non-standard Hubbard models the nonlinear coupling between the adjacent sites. In that part, discrete solitons are the main subject of the consideration. We conclude by listing some future open problems., Comment: expanded version 47pp, accepted in Rep. Prog. Phys

The main obstacle in experimental realization of the Einstein-de Haas effect in a Bose-Einstein condensate is necessity of a very precise control of the extremely small (of the order of tens of $\mu$G) external magnetic field. In this paper we numerically study the response of a rubidium condensate to an optimized time-dependent magnetic field. We find a significant transfer of atoms from the initial maximally polarized state to the next Zeeman component at magnetic fields of the order of tens of milligauss. We propose an experiment in which such an optimization scheme could enable the observation of the Einstein-de Haas effect in a rubidium atom condensate.

We numerically study the dynamics of a spinor chromium condensate in low magnetic fields. We show that the condensate evolution has a resonant character revealing rich structure of resonances similar to that already discussed in the case of alkali-atoms condensates. This indicates that dipolar resonances occur commonly in the systems of cold atoms. In fact, they have been already observed experimentally. We further simulate two recent experiments with chromium condensates, in which the threshold in spin relaxation and the spontaneous demagnetization phenomena were observed. We demonstrate that both these effects originate in resonant dynamics of chromium condensate., Comment: 8 pages, 9 figures

We study the ground-state properties of bosons loaded into the $p$-band of a one dimensional optical lattice. We show that the phase diagram of the system is substantially affected by the anharmonicity of the lattice potential. In particular, for a certain range of tunneling strength, the full many-body ground state of the system becomes degenerate. In this region, an additional symmetry of the system, namely the parity of the occupation number of the chosen orbital, is spontaneously broken. The state with nonvanishing staggered angular momentum, which breaks the time-reversal symmetry, becomes the true ground state of the system., Comment: 5 pages, 2 figures

We study a spin dynamics of two magnetic Chromium atoms trapped in a single site of a deep optical lattice in a resonant magnetic field. Dipole-dipole interactions couple spin degrees of freedom of the two particles to their motion in the site. The motion is quantized, therefore a trap geometry combined with two-body contact s-wave interactions influence a spin dynamics through the energy spectrum of the two atom system. Anharmonicity and anisotropy of the site results in a `fine' structure of two body eigenenergies. The structure can be easily resolved by a weak magnetic dipole-dipole interactions. As an example we examine the effect of anharmonicity and anisotropy of the binding potential on the Einstein-de Haas effect. We show that the weak dipolar interactions provide a perfect tool for a precision spectroscopy of the energy spectrum of the interacting few particle system., Comment: 14 pages, 8 figures

We study in detail the effect of quasicondensation. We show that this effect is strictly related to dimensionality of the system. It is present in one dimensional systems independently of interactions - exists in repulsive, attractive or in non-interacting Bose gas in some range of temperatures below characteristic temperature of the quantum degeneracy. Based on this observation we analyze the quasicondensation in terms of a ratio of the two largest eigenvalues of the single particle density matrix for the ideal gas. We show that in the thermodynamic limit in higher dimensions the second largest eigenvalue vanishes (as compared to the first one) with total number of particles as $\simeq N^{-\gamma}$ whereas goes to zero only logarithmically in one dimension. We also study the effect of quasicondensation for various geometries of the system: from quasi-1D elongated one, through spherically symmetric 3D case to quasi-2D pancake-like geometry.

Using classical field approximation we present the first study of statistical properties of one dimensional Bose gas with attractive interaction. The canonical probability distribution is generated with the help of a Monte Carlo method. This way we obtain not only the depletion of the condensate with growing temperature but also its fluctuations. The most important is our discovery of a reduced coherence length, the phenomenon observed earlier only for the repulsive gas, known as quasicondensation., Comment: 4 pages, 4 figures

We study a spinor condensate of Rb-87 atoms in F = 1 hyperfine state confined in an optical dipole trap. Putting initially all atoms in mF = 1 component we observe a significant transfer of atoms to other, initially empty Zeeman states exclusively due to dipolar forces. Because of conservation of a total angular momentum the atoms going to other Zeeman components acquire an orbital angular momentum and circulate around the center of the trap. This is a realization of Einstein-de Haas effect in a system of cold gases. We show that the transfer of atoms via dipolar interaction is possible only when the energies of the initial and the final sates are equal. This condition can be fulfilled utilizing a resonant external magnetic field, which tunes energies of involved states via the linear Zeeman effect. We found that there are many final states of different spatial density which can be tuned selectively to the initial state. We show a simple model explaining high selectivity and controllability of weak dipolar interactions in the condensate of Rb-87 atoms., Comment: 10 pages, 8 figures

We consider a simple experimental setup, based on a harmonic confinement, where a Bose-Einstein condensate and a thermal cloud of weakly interacting alkali atoms are trapped in two different vessels connected by a narrow channel. Using the classical field approximation, as described in J. Phys. B 40, R1 (2007) and optimized in Phys. Rev. A 81, 013629 (2010) for an arbitrary trapping potential, we theoretically investigate the analog of the celebrated superfluid helium fountain effect. We show that this thermo-mechanical effect might indeed be observed in this system. By analyzing the dynamics of the system, we are able to identify the superfluid and normal components of the flow as well as to distinguish the condensate fraction from the superfluid component. We show that the superfluid component can easily flow from the colder vessel to the hotter one while the normal component is practically blocked in the latter., Comment: 13 pages, 11 figures, 3 tables

Monte Carlo method within, so called, classical fields approximation is applied to one dimensional weakly interacting repulsive Bose gas trapped in a harmonic potential. Equilibrium statistical properties of the condensate are calculated within a canonical ensemble. We also calculate experimentally relevant low order correlation functions of the whole gas.

We study a square plaquette of four optical microtraps containing ultracold $^{87}$Rb atoms in F=1 hyperfine state. In a presence of external resonant magnetic field the dipolar interactions couple initial $m_F=1$ component to other Zeeman sublevels. This process is a generalization of the Einstein-de Haas effect to the case when the external potential has only $C_4$ point-symmetry. We observe that vortex structures appear in the initially empty $m_F=0$ state. Topological properties of this state are determined by competition between the local axial symmetry of the individual trap and the discrete symmetry of the plaquette. For deep microtraps vortices are localized at individual sites whereas for shallow traps only one discrete vortex appears in the plaquette. States created in these two opposite cases have different topological properties related to $C_4$ point-symmetry.

We study dynamics of two interacting ultra cold Bose atoms in a harmonic oscillator potential in one spatial dimension. Making use of the exact solution of the eigenvalue problem of a particle in the delta-like potential we study time evolution of initially separable state of two particles. The corresponding time dependent single particle density matrix is obtained and diagonalized and single particle orbitals are found. This allows to study decoherence as well as creation of entanglement during the dynamics. The evolution of the orbital corresponding to the largest eigenvalue is then compared to the evolution according to the Gross-Pitaevskii equation. We show that if initially the center of mass and relative degrees of freedom are entangled then the Gross-Pitaevskii equation fails to reproduce the exact dynamics and entanglement is produced dynamically. We stress that predictions of our study can be verified experimentally in an optical lattice in the low-tunneling limit., Comment: 9 figures, 5 movies available on-line

We investigate numerically the free-fall expansion of a $^{87}$Rb atoms condensate at nonzero temperatures. The classical field approximation is used to separate the condensate and the thermal cloud during the expansion. We calculate the radial and axial widths of the expanding condensate and find clear evidence that the thermal component changes the dynamics of the condensate. Our results are confronted against the experimental data.

We optimize the classical field approximation of the version described in J. Phys. B 40, R1 (2007) for the oscillations of a Bose gas trapped in a harmonic potential at nonzero temperatures, as experimentally investigated by Jin et al. [Phys. Rev. Lett. 78, 764 (1997)]. Similarly to experiment, the system response to external perturbations strongly depends on the initial temperature and on the symmetry of perturbation. While for lower temperatures the thermal cloud follows the condensed part, for higher temperatures the thermal atoms oscillate rather with their natural frequency, whereas the condensate exhibits a frequency shift toward the thermal cloud frequency (m=0 mode), or in the opposite direction (m=2 mode). In the latter case, for temperatures approaching critical, we find that the condensate begins to oscillate with the frequency of the thermal atoms, as in the m=0 mode. A broad range of frequencies of the perturbing potential is considered., Comment: 10 pages, 9 figures, added references, changed content, improved English