Quantum ground state of self-organized atomic crystals in optical resonators

Cold atoms, driven by a laser and simultaneously coupled to the quantum field of an optical resonator, may self-organize in periodic structures. These structures are supported by the optical lattice, which emerges from the laser light they scatter into the cavity mode and form when the laser intensity exceeds a threshold value. We study theoretically the quantum ground state of these structures above the pump threshold of self-organization by mapping the atomic dynamics of the self-organized crystal to a Bose-Hubbard model. We find that the quantum ground state of the self-organized structure can be the one of a Mott insulator, depending on the pump strength of the driving laser. For very large pump strengths, where the intracavity-field intensity is maximum and one would expect a Mott-insulator state, we find intervals of parameters where the phase is compressible. These states could be realized in existing experimental setups.

Low-energy behavior of strongly interacting bosons on a flat-band lattice above the critical filling factor

Bosons interacting repulsively on a lattice with a flat lowest band energy dispersion may, at sufficiently small filling factors, enter into a Wigner-crystal-like phase. This phase is a consequence of the dispersionless nature of the system, which in turn implies the occurrence of single-particle localized eigenstates. We investigate one of these systems-the sawtooth lattice-filled with strongly repulsive bosons at filling factors infinitesimally above the critical point where the crystal phase is no longer the ground state. We find, in the hard-core limit, that the crystal retains its structure in all but one of its cells, where it is broken. The broken cell corresponds to an exotic kind of repulsively bound state, which becomes delocalized. We investigate the excitation spectrum of the system analytically and find that the bound state behaves as a single particle hopping on an effective lattice with reduced periodicity, and is therefore gapless. Thus, the addition of a single particle to a flat-band system at critical filling is found to be enough to make kinetic behavior manifest.

Bose-Einstein condensation dynamics in three dimensions by the pseudospectral and finite-difference methods

We suggest a pseudospectral method for solving the three-dimensional time-dependent Gross-Pitaevskii (GP) equation, and use it to study the resonance dynamics of a trapped Bose-Einstein condensate induced by a periodic variation in the atomic scattering length. When the frequency of oscillation of the scattering length is an even multiple of one of the trapping frequencies along the x, y or z direction, the corresponding size of the condensate executes resonant oscillation. Using the concept of the differentiation matrix, the partial-differential GP equation is reduced to a set of coupled ordinary differential equations, which is solved by a fourth-order adaptive step-size control Runge-Kutta method. The pseudospectral method is contrasted with the finite-difference method for the same problem, where the time evolution is performed by the Crank-Nicholson algorithm. The latter method is illustrated to be more suitable for a three-dimensional standing-wave optical-lattice trapping potential.

Tightly bound gap solitons in a Fermi gas

Within the framework of the mean-field hydrodynamic model of a degenerate Fermi gas ( DFG), we study, by means of numerical methods and variational approximation ( VA), the formation of fundamental gap solitons ( FGSs) in a DFG ( or in a BCS superfluid generated by weak interaction between spin- up and spin- down fermions), which is trapped in a periodic optical- lattice ( OL) potential. An effectively one- dimensional ( 1D) con. guration is considered, assuming strong transverse confinement; in parallel, a proper 1D model of the DFG ( which amounts to the known quintic equation for the Tonks- Girardeau gas in the OL) is considered too. The FGSs found in the first two bandgaps of the OL- induced spectrum ( unless they are very close to edges of the gaps) feature a ( tightly bound) shape, being essentially confined to a single cell of the OL. In the second bandgap, we also find antisymmetric tightly bound subfundamental solitons ( SFSs), with zero at the midpoint. The SFSs are also confined to a single cell of the OL, but, unlike the FGSs, they are unstable. The predicted solitons, consisting of similar to 10(4) - 10(5) atoms, can be created by available experimental techniques in the DFG of Li-6 atoms.

Mean-field model for the interference of matter-waves from a three-dimensional optical trap

Using the mean-field time-dependent Gross-Pitaevskii equation we study the formation of a repulsive Bose-Einstein condensate on a combined optical and harmonic traps in two and three dimensions and subsequent generation of the interference pattern upon the removal of the combined traps as in the experiment by, Greiner et al. [Nature (London 415 (2002) 39]. For optical traps of moderate strength, interference pattern of 27 (9) prominent bright spots is found to be formed in three. (two) dimensions on a cubic (square) lattice in agreement with experiment. Similar interference pattern can also be formed upon removal of the optical lattice trap only. The pattern so formed can oscillate for a long time in the harmonic trap which can be observed experimentally. (C) 2003 Elsevier B.V. B.V. All rights reserved.

Three-dimensional solitons in cross-combined linear and nonlinear optical lattices

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Josephson oscillation of a superfluid Fermi gas

Using the complete numerical solution of a time-dependent three-dimensional rnean-field model we study the Josephson oscillation of a superfluid Fermi gas (SFG) at zero temperature formed in a combined axially-symmetric harmonic plus one-dimensional periodic optical-lattice (OL) potentials after displacing the harmonic trap along the axial OL axis. We study the dependence of Josephson frequency on the strength of the OL potential. The Josephson frequency decreases with increasing strength as found in the experiment of Cataliotti et al. [Science 293, 843 (2001)] for a Bose-Einstein condensate and of the experiment of Pezze et al. [Phys. Rev. Lett. 93, 120401 (2004)] for an ideal Fermi gas. We demonstrate a breakdown of Josephson oscillation in the SFG for a large displacement of the harmonic trap. These features of Josephson oscillation of a SFG can be tested experimentally.

Gap solitons in fermion superfluids

We consider a dynamical model of a superfluid Fermi gas in the Bardeen-Cooper-Schrieffer regime trapped in a periodic optical lattice (OL) potential. The model is based on an equation for complex order parameter phi of the superfluid, which is derived from the relevant energy density and includes a self-repulsive term similar to phi(7/3). By means of the variational approximation (VA) and numerical simulations, we find families of stable one- and two-dimensional (I D and 2D) gap solitons (GSs) in this model. Chiefly, they are compact objects trapped in a single cell of the OL. Families of stable even and odd bound states of these GSs are also found in one dimension. A 3D GS family is constructed too, but solely within the framework of the VA. In the linear limit, the VA predicts an almost exact position of the left edge of the first band-gap in the OL-induced spectrum. The full VA provides an accurate description of families of I D and 2D fundamental GSs. We also demonstrate that a I D GS can be safely transported by an OL moving at a moderate velocity. (C) 2009 IMACS. Published by Elsevier B.V. All rights reserved.

Matter-wave two-dimensional solitons in crossed linear and nonlinear optical lattices

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Gap solitons in superfluid boson-fermion mixtures

Using coupled equations for the bosonic and fermionic order parameters, we construct families of gap solitons (GSs) in a nearly one-dimensional Bose-Fermi mixture trapped in a periodic optical-lattice (OL) potential, the boson and fermion components being in the states of the Bose-Einstein condensation and Bardeen-Cooper-Schrieffer superfluid, respectively. Fundamental GSs are compact states trapped, essentially, in a single cell of the lattice. Full families of such solutions are constructed in the first two band gaps of the OL-induced spectrum, by means of variational and numerical methods, which are found to be in good agreement. The families include both intragap and intergap solitons, with the chemical potentials of the boson and fermion components falling in the same or different band gaps, respectively. Nonfundamental states, extended over several lattice cells, are constructed too. The GSs are stable against strong perturbations.

The dynamics of bright matter wave solitons in a quasi one-dimensional Bose-Einstein condensate with a rapidly varying trap

The dynamics of a bright matter wave soliton in a quasi one-dimensional Bose-Einstein condensate (BEC) with a periodically rapidly varying time trap is considered. The governing equation is based on averaging the fast modulations of the Gross-Pitaevskii (GP) equation. This equation has the form of a GP equation with an effective potential of a more complicated structure than an unperturbed trap. In the case of an inverted (expulsive) quadratic trap corresponding to an unstable GP equation, the effective potential can be stable. For the bounded space trap potential it is showed that bifurcation exists, i.e. the single-well potential bifurcates to the triple-well effective potential. The stabilization of a BEC cloud on-site state in the temporary modulated optical lattice is found. This phenomenon is analogous to the Kapitza stabilization of an inverted pendulum. The analytical predictions of the averaged GP equation are confirmed by numerical simulations of the full GP equation with rapid perturbations.

Three-dimensional solitons in cross-combined linear and nonlinear optical lattices

The existence and stability of three-dimensional (3D) solitons, in cross-combined linear and nonlinear optical lattices, are investigated. In particular, with a starting optical lattice (OL) configuration such that it is linear in the x-direction and nonlinear in the y-direction, we consider the z-direction either unconstrained (quasi-2D OL case) or with another linear OL (full 3D case). We perform this study both analytically and numerically: analytically by a variational approach based on a Gaussian ansatz for the soliton wavefunction and numerically by relaxation methods and direct integrations of the corresponding Gross-Pitaevskii equation. We conclude that, while 3D solitons in the quasi-2D OL case are always unstable, the addition of another linear OL in the z-direction allows us to stabilize 3D solitons both for attractive and repulsive mean interactions. From our results, we suggest the possible use of spatial modulations of the nonlinearity in one of the directions as a tool for the management of stable 3D solitons.

Interacting fermionic atoms in optical lattices - A quantum simulator for condensed matter physics

This thesis reports on the creation and analysis of many-body states of interacting fermionic atoms in optical lattices. The realized system can be described by the Fermi-Hubbard hamiltonian, which is an important model for correlated electrons in modern condensed matter physics. In this way, ultra-cold atoms can be utilized as a quantum simulator to study solid state phenomena. The use of a Feshbach resonance in combination with a blue-detuned optical lattice and a red-detuned dipole trap enables an independent control over all relevant parameters in the many-body hamiltonian. By measuring the in-situ density distribution and doublon fraction it has been possible to identify both metallic and insulating phases in the repulsive Hubbard model, including the experimental observation of the fermionic Mott insulator. In the attractive case, the appearance of strong correlations has been detected via an anomalous expansion of the cloud that is caused by the formation of non-condensed pairs. By monitoring the in-situ density distribution of initially localized atoms during the free expansion in a homogeneous optical lattice, a strong influence of interactions on the out-of-equilibrium dynamics within the Hubbard model has been found. The reported experiments pave the way for future studies on magnetic order and fermionic superfluidity in a clean and well-controlled experimental system.

Quantum simulation of a topological Mott insulator with Rydberg atoms in a Lieb lattice

We propose a realistic scheme to quantum simulate the so-far experimentally unobserved topological Mott insulator phase-an interaction-driven topological insulator-using cold atoms in an optical Lieb lattice. To this end, we study a system of spinless fermions in a Lieb lattice, exhibiting repulsive nearest-and next-to-nearest-neighbor interactions and derive the associated zero-temperature phase diagram within mean-field approximation. In particular, we analyze how the interactions can dynamically generate a charge density wave ordered, a nematic, and a topologically nontrivial quantum anomalous Hall phase. We characterize the topology of the different phases by the Chern number and discuss the possibility of phase coexistence. Based on the identified phases, we propose a realistic implementation of this model using cold Rydberg-dressed atoms in an optical lattice. The scheme, which allows one to access, in particular, the topological Mott insulator phase, robustly and independently of its exact position in parameter space, merely requires global, always-on off-resonant laser coupling to Rydberg states and is feasible with state-of-the-art experimental techniques that have already been demonstrated in the laboratory.

Condensados de Bose-Einstein em redes óticas: a transição superfluido-isolante de Mott em redes hexagonais e a classe de universalidade superfluido-vidro de Bose em 3D

Estudamos transições de fases quânticas em gases bosônicos ultrafrios aprisionados em redes óticas. A física desses sistemas é capturada por um modelo do tipo Bose-Hubbard que, no caso de um sistema sem desordem, em que os átomos têm interação de curto alcance e o tunelamento é apenas entre sítios primeiros vizinhos, prevê a transição de fases quântica superfluido-isolante de Mott (SF-MI) quando a profundidade do potencial da rede ótica é variado. Num primeiro estudo, verificamos como o diagrama de fases dessa transição muda quando passamos de uma rede quadrada para uma hexagonal. Num segundo, investigamos como a desordem modifica essa transição. No estudo com rede hexagonal, apresentamos o diagrama de fases da transição SF-MI e uma estimativa para o ponto crítico do primeiro lobo de Mott. Esses resultados foram obtidos usando o algoritmo de Monte Carlo quântico denominado Worm. Comparamos nossos resultados com os obtidos a partir de uma aproximação de campo médio e com os de um sistema com uma rede ótica quadrada. Ao introduzir desordem no sistema, uma nova fase emerge no diagrama de fases do estado fundamental intermediando a fase superfluida e a isolante de Mott. Essa nova fase é conhecida como vidro de Bose (BG) e a transição de fases quântica SF-BG que ocorre nesse sistema gerou muitas controvérsias desde seus primeiros estudos iniciados no fim dos anos 80. Apesar dos avanços em direção ao entendimento completo desta transição, a caracterização básica das suas propriedades críticas ainda é debatida. O que motivou nosso estudo, foi a publicação de resultados experimentais e numéricos em sistemas tridimensionais [Yu et al. Nature 489, 379 (2012), Yu et al. PRB 86, 134421 (2012)] que violam a lei de escala $\\phi= u z$, em que $\\phi$ é o expoente da temperatura crítica, $z$ é o expoente crítico dinâmico e $ u$ é o expoente do comprimento de correlação. Abordamos essa controvérsia numericamente fazendo uma análise de escalonamento finito usando o algoritmo Worm nas suas versões quântica e clássica. Nossos resultados demonstram que trabalhos anteriores sobre a dependência da temperatura de transição superfluido-líquido normal com o potencial químico (ou campo magnético, em sistemas de spin), $T_c \\propto (\\mu-\\mu_c)^\\phi$, estavam equivocados na interpretação de um comportamento transiente na aproximação da região crítica genuína. Quando os parâmetros do modelo são modificados de maneira a ampliar a região crítica quântica, simulações com ambos os modelos clássico e quântico revelam que a lei de escala $\\phi= u z$ [com $\\phi=2.7(2)$, $z=3$ e $ u = 0.88(5)$] é válida. Também estimamos o expoente crítico do parâmetro de ordem, encontrando $\\beta=1.5(2)$.