Arxiv:1804.10133V3 [Cond-Mat.Quant-Gas] 17 Dec 2018

Rise and fall of a bright soliton in an optical lattice Piero Naldesi1 ∗, Juan Polo Gomez1, Boris Malomed2;3, Maxim Olshanii4, Anna Minguzzi1 and Luigi Amico5;6;7;8;9 1 UniversitéGrenoble-Alpes, LPMMC, F-38000 Grenoble, France and CNRS, LPMMC, F-38000 Grenoble, France 2 Department of Physical Electronics, School of Electrical Engineering, Faculty of Engineering, Tel Aviv University, P.O.B. 39040, Ramat Aviv, Tel Aviv, Israel 3 Center for Light-Matter Interaction, Tel Aviv University, P.O.B. 39040, Ramat Aviv, Tel Aviv, Israel 4 Department of Physics, University of Massachusetts Boston, Boston, MA 02125, USA 5 Dipartimento di Fisica e Astronomia, Via S. Sofia 64, 95127 Catania, Italy 6 Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, Singapore 117543, Singapore 7 MajuLab, CNRS-UNS-NUS-NTU International Joint Research Unit, UMI 3654, Singapore 8 CNR-MATIS-IMM & INFN-Sezione di Catania, Via S. Sofia 64, 95127 Catania, Italy and 9 LANEF 'Chaire d'excellence', UniversitèGrenoble-Alpes & CNRS, F-38000 Grenoble, France (Dated: December 19, 2018) We study an ultracold atomic gas with attractive interactions in a one-dimensional optical lattice. We find that its excitation spectrum displays a quantum soliton band, corresponding to N-particle bound states, and a continuum band of other, mostly extended, states. For a system of a finite size, the two branches are degenerate in energy for weak interactions, while a gap opens above a threshold value of the interaction strength. We find that the interplay between degenerate extended and bound states has important consequences for both static and dynamical properties of the system. In particular, the solitonic states turn out to be protected from spatial perturbations and random disorder. We discuss how such dynamics implies that our system effectively provides an example of a quantum many-body system that, with the variation of the bosonic lattice filling, crosses over from integrable non-ergodic to non-integrable ergodic dynamics, through non-integrable non-ergodic regimes. Introduction and summary of the results. Solitary was demonstrated that the ground state energy may dis- waves in classical fluids may arise when wave disper- play instabilities that can be nonetheless cured by a suit- sion effects are compensated by non-linear interactions able choice of interactions and density [22]. Indeed, the [1]. The study of their mathematical properties has de- limits of vanishing interaction with a finite density, lead- fined important domains in mathematical research [2{4], ing to mean-field results for large number of bosons [23], with far reaching implications for pure and applied mod- or of vanishing density with finite interactions [24], have ern science [5{9]. Although quantum mechanics is an been thoroughly explored. In particular, by analysing intrinsically linear theory, solitons may emerge also in the solutions of the mean-field Gross-Pitaevskii equation quantum fluids: in this case, non-linearities do arise as it was found that [25, 26] a critical value of the attrac- a result of effective cooperative phenomena occurring in tion exists for which the ground state density undergoes quantum-many-particle systems. Indeed, solitons were to a transition from a uniform profile to a bright-soliton- demonstrated to emerge in different quantum mechani- type one (with the logic implied by the onset of modula- cal contexts, ranging from quantum material science to tional instabilities in the condensate). In the same limit, particle physics [10{15]. it was exactly demonstrated using the Bethe-Ansatz so- lution that density-density and higher-order correlation In this work, we are primarily motivated by the recent functions display a qualitative change of behaviour in cor- investigations in quantum fluids as provided by ultracold respondence to the critical value predicted by the mean- atoms trapped in one dimensional optical potentials. In field theory. these systems, bright solitons may emerge for attractive atom-atom interactions [3] (see also [16]) and they have In this work we focus on a bosonic system described been observed in several experiments [17{20]. From a by the Bose-Hubbard Model (BHM) [27, 28] (1) confined conceptual point of view, however, bosonic systems with in a one dimensional lattice, where density and inter- attractive interactions must be treated with care: be- actions strength are both finite. We perform a numeri- cause of the Bose statistics, the lowest energy state can be cal study of energy bands and quantum correlations us- arXiv:1804.10133v3 [cond-mat.quant-gas] 17 Dec 2018 macroscopically occupied with a density that is magnified ing the Density Matrix Renormalization Group (DMRG) by interactions. Important progress has been achieved method [29{33]. Among different aspects implied by the describing the bosonic fluid through a famous integrable lattice, here we exploit the energy-band structure of the theory as proposed by the Lieb-Liniger model, which is system, featuring characteristic bendings, foldings and amenable to an exact analysis [21]. Relying on that, it energy gaps. Such effects can indeed define new physical regimes in our system with peculiar bound states of solitonic type (see [34{36]). For finite attractive interactions we find that N-body ∗Electronic address: [email protected] bound states are formed; for weak interactions, however, 2 these bound states are degenerate with a second band dimensional lattice. The Hamiltonian reads of other states, mostly extended or involving lower-order L L bound states. By increasing the interaction strength, the X y jUj X H^ = −J a a +1 + h.c. − n^ (^n − 1) (1) bound states get more and more energetically favourable, j j 2 j j j=1 j=1 until a critical interaction strength Uc for which the band of bound states is completely separated by an energy gap y from the rest of the spectrum. where the operators ni = ai ai count the number of y bosons at the site i; the operators ai, ai obey the canon- The calculation of the dynamical structure factor (2) y ical commutation relations [a ; a ] = δ , and L is the provides the portrait of the band structure of the sys- i j ij tem, displaying a two-branch dispersion at low ener- number of sites. The parameters J, U in (1) are the gies as shown in Fig. 1. Such quantity is experimen- hopping amplitude and the strength of the on-site in- tally detectable by means in ultracold-atoms experiments teraction, respectively. Throughout this article, we will by Bragg scattering methods [37, 38]. We characterize consider only attractive interactions; both energies and the nature of ground and excited states in the spectrum the parameter U are in units of J. Times are in units of by monitoring the density-density correlations functions, ~=J. Equation (1) describes a closed system and there- which display a different spatial behaviour for extended fore it neglects 3-body losses. We will also always (un- less stated otherwise) consider open boundary conditions: and N-particle bound states, see Fig. 2. y y a1aL = aLa1 = 0. The BHM (1) is not solvable by the The changing of band structure and the opening of a coordinate Bethe ansatz. The failure results because of gap varying the interaction, has important consequences finite probabilities that a given site is occupied by more for the dynamics of the system. We devise a protocol than two particles, whose interaction cannot be factorized in which we prepare a quantum soliton in the middle of in 2-body scattering [40{42]. Nevertheless the plane-wave the chain and then we let it expand under the guidance ansatz of coordinate Bethe Ansatz works well in the so of Hamiltonian (1). For U = 0 only extended states are called two-particle sector, for which such probabilities are available in the dynamics; for 0 < U < Uc the extended vanishing[39, 43]. Despite the fact that the BHM is not and bound states are available; for U > Uc, at low ener- integrable, its continuous limit is the Bose-gas integrable gies, solely bound states exist. As a striking feature, in field theory [44] (see Supplemental Material). We note this latter regime the density keeps more and more the that a similar logic works also for discretization for the localized shape of the initial state and only a small frac- classical non-linear-Schroedinger equation [45]. tion of the state spreads over the lattice. The crossover Bound versus scattering states. Information on the between the two regimes is purely mesoscopic since Uc available excitations in the system as a function of their scales like the inverse of the number of particles (see sup- momentum k and energy ! is provided by the dynamical plemental material). We note that the expansion velocity structure factor S(k; !): experiences a crossover from a large value for U < Uc to a smaller value for for U >Uc. These features, that have X X −ikr 2 S(k; !) = j hαj e n^r j0i j δ(! − !α): (2) been predicted for two particles [39], clearly emerge in α6=0 r Fig. 4. Even though no quantum phase transition oc- curs in the system, we find that the expansion velocity wheren ^r is the number operator acting on the site r, close to Uc displays a universal scaling. Such a feature j0i is the ground state and α labels the states with in- illustrates the subtle interplay between interactions and creasing energy (ie α = 1 is the first excited state). The particle number in the dynamics of attractive bosons in peaks of S(k; !) reconstruct the energy bands of the sys- a lattice. tem [46, 47]. We observe that for small U a low-energy We also find that the occurrence of degenerate scatter- band separates from the rest of the spectrum.

Load more