Unusual discrete soliton and breather modes collective excitations in Bose-Einstein condensates in optical lattice

1. Unusual discrete soliton and breather modes collective excitations in Bose-Einstein condensates in optical lattice Bishwajyoti Dey Department of Physics University of Pune, Pune with Galal Al- Akhaly

2. Nonlinear localised excitations – solitons, breathers, compactons. These solutions are nonspreading – retain their shape in time. Solitons, breathers form if the peaking due to nonlinearity is balanced by the spreading due to linear dispersion. For compactons balance with nonlinear dispersion is required. For discrete systems the localization is due to the discreteness combined with the nonlinearity of the system. For linear systems, the discrete translational invariance have to be broken (adding impurity) to obtain spatially localized mode (Anderson Localization). For nonlinear systems one can retain discrete translational symmetry and still obtain localized excitations. Self localised solutions. Bright solitons have been observed in BEC where the linear spreading due to dispersion is compensated by the attractive nonlinear interactions between the atoms.

3. Solitons and BreathersNonlinear energy localization in continuous media is well studied • Small amplitude, long waves in a channel of slowly changing depths • Kortweg de Vries (KdV) equation Subscripts denote partial derivatives. q(x,t) wave shape. Has spatially localized solutions called “ Soliton” Two soliton pass through Each other without loss of identity

4. 7

5. Breathers in continuum systems Breather solutions are known for PDEs. Nongeneric • Well known sine-Gordon equation has breather solution

6. Breather Excitations in continuous systems (described by PDE’s) are non generic • Decade-long efforts to look for localized excitations like breathers in other partial differential equations (PDE’s) which describes evolution of physical fieldsproduced no result. • Reason is simple. Resonances of the nonlinear localized excitation with spectrum of plane waves (linear) of the system. • This resonance problem can be avoided if we consider a discrete system (LATTICE) instead of the continuum. • In discrete lattice the plane wave frequency have always finite upper bound. We can excite the nonlinear localized excitation above the upper bound frequency to avoid resonances.

7. Fermi-Pasta-Ulam lattice, V=0 and W is anharmonic. Klein-Gordon lattice, etc., V=nonlinear, W is harmonic. Discrete Nonlinear Schrodinger (DNLS) Lattice. Ablowitz-Ladik lattice Model. Salerno Lattice Model. More recently, dynamics of Bose-Einstein Condensate loaded on an deep optical lattice. Nonlinear Lattices

8. Nonlinear Schrodinger equation with cubic nonlinearity (PDE) Discrete nonlinear Schrodinger equation (nonlinear lattice model) Ablowitz-Ladik equation (integrable lattice model)

9. Salerno equation DNLS equation for the dynamics of Bose-Einstein condensate amplitude

10. The Lattice Problem : nonlinear lattice • Spatial discreteness and Nonlinearity For nonlinear lattice, onsite potential can be nonlinear, or W (intersite interaction) can be nonlinear (anharmonic) or both can be nonlinear. Linearize equation of motion around classical ground state

11. Dispersion relation for small amplitude plane wave

12. Page mode (right), Sievers-Takeno Mode (left) Two-site DB (right)

13. ENERGY LOCALIZATION : Local energy density If DBs are excited, initial local energy stay within the DB. The function should notdecay with time. Dashed line e5 vs time Solid line –Total energy of the chain e5

14. First unambiguous observation of BEC was reported by Eric Cornell, Carl Wieman (1995) in Colorado (US). BEC was observed cooling a gas of rubidium-87 to a temperature 170nK Fig. Velocity distribution. The axes are x and z velocities and third axis is number density of atoms. Macroscopic fraction (~10%) of the atoms are in the ground state.

15. Optical lattice: an artificial crystal of light – a periodic intensity pattern that is formed by the interference of two or more laser beams. More lasers give 3D spatial structure. Trapping atom in optical lattice – atoms can be trapped in the bright or dark regions of the optical lattice via Stark shift. Strength of the optical potential confining can be increased by increasing laser intensity. BEC mounted on a optical lattice is like electrons in a periodic potential of ions in conventional solid. Condensate atoms plays the role of electrons and optical lattice the role of ions.

16. Atoms trapped in an optical lattice move due to quantum tunneling even if the potential depth of the lattice point exceeds the kinetic energy. Strongly interacting limit. However when the well depth is large then the interaction energy of the atoms become more than the hopping energy, then the atoms will be trapped in potential minima and cannot move freely. This phase is called Mott insulator. Atoms in an optical lattice provide an ideal quantum system where all parameters can be controlled. This can be used to observe effects which are difficult to observe in real crystals. Examples: Bloch oscillation, Efimov effect, Superfluid to Mott insulator transition Superconductivity etc.

17. Cold atom experiments tread into the land of two-dimensional superconductivity M. Randeria Physics 5, 10 (2012). Sommer et al. Phys. Rev. Lett 108, 045302 (2012). A deeper understanding of strongly interacting two-dimensional superconductors and their normal states can give insights into high Tc superconductors. 2D attractive Fermi-gas has rich property even in normal state. Pair binding is enhanced in 2D. Pairing can occur without condensation over a large temperature gap leading to pseudogap effects. Analytical meanfield calculations showed that in 2D the binding energy of many-body systems is same as the binding energy of just a pair. This is verified by experiment showing great prediction power of mean-field theory.

18. Sommer et al. [1] use radio-frequency (rf) spectroscopy to determine the binding energy of paired lithium atoms in a cold gas. This cartoon shows a gas with two species of fermions, denoted by red and blue, which are analogous to the spin-up and spin-down electrons in a metal. The presence of an optical lattice potential (blue curve) tunes the dimensionality and forces the gas into stacks of quasi-two-dimensional layers. The interaction between the two types of fermions is tuned using a Feshbach resonance. Absorption of the rf photon (red wavy line) converts a fermion in the red state to a different internal state shown as green, and measures the pair binding energy.

19. Bloch Oscillation The quantum dynamics of accelerated particles in periodic potential leads to an oscillatory motion instead of a linear increase in velocity. This is termed as Bloch oscillation The periodicity of the potential implies eigenfunctions obey relation In presence of an accelerating force F, the quasimomentum evolves linearly in time In combination with the periodicity of the band structure, this causes an oscillatory motion, the Bloch oscillation. The oscillation period is .

20. Bloch Oscillation In solid state systems scattering due to impurity of the crystals structure leads to damping of Bloch oscillations on time scales much shorter than the oscillation period itself. Hence difficult to observe experimentally. Optical lattice on the other hand constitute a perfect optical crystal and BEC on optical lattice have enabled the first direct observation of Bloch oscillation. Due to interactions between atoms, Bloch oscillation decays from dynamical instabilities.

21. Efimov effect Quantum Mechanics of three-body systems : Efimov in 1970 predicted that there can exist bound states (Efimov states) of three particles even if the two-particle attraction is too weak to allow two particles to form a pair. The sequence of three-body bound states have universal properties, it is insensitive to the details of two-body potential at short distances. Efimov’s theoretical prediction could only be verified experimentally in 2006 in ultra cold gas of cesium atoms 36 years after its predictions ( Kraemer et al, Nature 440, 315 (2006)). Since then, Efimov effect have been observed in other BEC’s, Bose-Bose, Bose-Fermi, Fermi-Fermi and bosonic dipoles ( Ferlaino and Grimm, Physics 3, 9 (2010)).

22. Efimov’s prediction that the binding of few particles is universal has also been confirmed experimentally (Berninger et al, Phys. Rev. Lett 107, 120401 (2011)). It is this concept of universality which makes Efimov physics possible in various physical systems, such as, atomic physics, nuclear physics, strongly correlated system etc, where the relevant energy and length scales differ by many order of magnitude.

23. Dynamics of BEC in an optical lattice: order parameter and mean-field theory The many-body Hamiltonian describing N-interacting Bosons confined by an external potential is given by where are boson field operators, is the two-body interaction potential. The field operators can be written as where are the single-particle wave function and are the corresponding annihilation operators defined as, with commutation rules

24. Using the Heisenberg equation the time evolution of the field operator is given by Bogoliubov first order theory for the excitations of interacting Bose-gas where is a classical field, the order parameter or the wave function of the condensate. The condensate density . Assuming that only binary collisions at low energy are relevant and these collisions are characterized by a single parameter, the s-wave scattering length, independent of the details of the two-body potential, we replace The coupling constant where a is the scattering length.

25. This yields the equation for the order parameter, the Gross-Pitaevskii equation The GP equation can be written as where A system of atoms with attractive two-body interactions, is unstable against collapse above certain critical number of atoms Nc. An addition of a repulsive three-body interaction can overcome the collapse and region of stability for the condensate can be extended beyond Nc.

26. In presence of three-body interactions the Gross-Pitaevskii equation become which depend on the Hamiltonian of a single trapped atom as well as two- and three-body coupling constants g2 and g3. The three-body coupling constant g3 has been derived from a microscopic theory of three-body collisions in a BEC (Kohler, PRL (2002)). The spatial coordinates are chosen as the vector from atom 1 to atom 2 ( r12) and the vector from the center of mass of atoms 1 and 2 to atom 3

27. Dynamics of BEC: Gross- Pitaevskii (GP) equation – treating the condensate as classical field. GP equation is a variant of the Nonlinear Schrodinger equation (NLS) incorporating an external potential used to confine the condensate. Multicomponent GP equation for spinor condensate. Dimensionality reduction possible in the presence of external periodic potential generated by the optical lattices and in the discrete limit.

28. Deep periodic optical potential limit – tight binding model The linear Bloch waves exhibit strong localization in the deep potential limit. Condensate wave function is described with localized Wannier states associated with lowest band. where is the condensate wave function localized in trap n with the orthonormal conditions Using in GP equation and integrating using the orthonormal conditions above we get the dynamics of the condensate described by the discrete nonlinear Schrodinger equation as (DNLS)

29. DNLS (discrete nonlinear Schrödinger) Equation where

30. The DNLS equation is the equation of motion and can be derived from the Hamiltonian where and are conjugate variable. Both the Hamiltonian and thenorm are conserved quantity.

31. Variational Dynamics To study the dynamical regime of a high density BEC in an array, we consider dynamical evolution of a Gaussian profile wave packet and introduce the variational wave function where the variational parameters and are center and width respectively of the density and and are their associated momenta. The dynamical evolution of the variational wave packet can be obtained by a variational principle from the Lagrangian

32. Using Euler-Lagrange equation the variational equations of motion are The pairs and are conjugate dynamical variables w.r.t. the Hamiltonian The variational equations can be solved numerically to obtain the variational dynamics of the system.

33. The wave packet group velocity is given by and the inverse effective mass is given by where

34. Numerical solution of DNLS equation DNLS equation is also solved numerically to compare with the variational dynamics results and also to check stability of the dynamics and phase diagrams over long period of time. Write the amplitude of the order parameter in terms of two components and . DNLS then can be written as The coupled nonlinear equations are solved using Runge-Kutta method .

35. The variational wave function is used as initial condition and The Hamiltonian and the norm are checked at each steps of the integration to look for their constancy over time.

36. Dynamical instability of Bloch oscillation: In presence of an accelerating force (tilted wash board potential) the quasimomentum is , where . Linear regime: for zero condensate interactions , the center of condensate oscillate as (exact solution) Similarly, the width of the condensate density oscillates as (exact solution) No instability (decay) in linear regime.

37. Numerically, for Bloch oscillation, we calculate numerical average position defined as It is easy to show that , the average position of the center of density. Similarly, gives the numerical width of the wavepacket.

38. Bloch Oscillation: no instability in absence of interactions

39. Nonlinear regime: In this case the equation for the center of the density is given by Note : even though there is a damping term in the equation, the dynamics is fully Hamiltonian. The apparent damping is due to the divergence of the effective mass with time due to which the Bloch oscillation decays.

40. The Bloch oscillation decays as Decay of Bloch oscillation: effect of nonlinear Interactions.(GA,BD 2011) Anderson & Kasevich, Science (1998).

41. Phase diagram of the interacting BEC BEC with deep optical lattice potential supports many interesting phases. Phase diagrams can be obtained from the coupled variational equations and the corresponding Hamiltonian.

42. The trajectories in the plane can be obtained as The condition implies for . is obtained from the condition . For and which implies This gives and The wave packet stops as the effective mass goes to infinity. This corresponds to the self-trapped regime in the phase diagram.

43. On the other hand, for , But and the effective mass There is complete spreading of the wave packet giving rise to the diffusive regime. The critical line separating these two regime (the self-trapped and the diffusive) is obtained from the condition as

44. (GA,BD 2011) Self -trapping Diffusion

45. SOLITON For negative effective mass, i.e. for we get another interesting phase from the fixed point of the trajectory. This gives a regime in the phase diagram where soliton solutions are allowed. The center of mass moves with constant velocity and the shape of the wavepacket do not change with time. Soliton solutions are allowed for the parameter values For there are no soliton solutions, as in this case the trajectory do not have fixed points.

46. Soliton solution from direct numerical integration of the Gross-Pitaevskii equation. (GA,BD 2011)

47. DISCRETE BREATHER Another interesting phase is the discrete breather which is a spatially localized and time-periodic solution. In this case oscillate with time. The trajectories in the plane are closed. We have discrete breather solution with center of mass travelling with nearly constant velocity and with oscillating width. oscillate around constant value.

48. Phase space trajectories

49. Phase diagram of the scalar BEC in deep optical lattice (GA,BD 2011) Parameters: Inset:

50. Phase diagram in the Efimov region : Numerical results shows that the soliton exist only for large value of . For large value of , the soliton line approach the critical line .