Oscillatory stability of quantum droplets in PT-symmetric optical lattice

Jiehong Lao,Zheng Zhou,Xili Zhang,Fuqiu Ye and Honghua Zhong

1 Department of Physics and Key Laboratory of Low-dimensional Quantum Structures and Quantum Control of Ministry of Education,Hunan Normal University,Changsha 410081,China

2 Department of Physics,Hunan Institute of Technology,Hengyang 421002,China

3 Department of Physics,Jishou University,Jishou 416000,China

4 Institute of Mathematics and Physics,Central South University of Forestry and Technology,Changsha 410004,China

Abstract We study stability and collisions of quantum droplets(QDs)forming in a binary bosonic condensate trapped in parity-time(PT)-symmetric optical lattices.It is found that the stability of QDs in the PT-symmetric system depends strongly on the values of the imaginary part W0 of the PT-symmetric optical lattices,self-repulsion strength g,and the condensate norm N.As expected,thePT-symmetric QDs are entirely unstable in the brokenPT-symmetric phase.However,the PT-symmetric QDs exhibit oscillatory stability with the increase of N and g in the unbroken PT-symmetric phase.Finally,collisions betweenPT-symmetric QDs are considered.The collisions of droplets with unequal norms are completely different from that in free space.Besides,a stable PT-symmetric QDs collides with an unstable ones tend to merge into breathers after the collision.

Keywords:quantum droplets,PT symmetry,optical lattice

1.Introduction

In recent years,a new type of self-bound quantum droplets(QDs)has attracted much attention in the field of ultracold atoms.Based on the theoretical suggestions described in the[1],QDs were experimentally created in dipolar bosonic gases of erbium[2],in dysprosium[3],and in mixtures of two atomic states of39K with appropriate signs of the intercomponent and intra-component interactions[4].Owing to the balance between the attraction generated by the dipole–dipole interaction and the Lee−Huang−Yang(LHY)repulsion caused by the quantum fluctuation[5],QDs has been observed in the experiment of the dipole boson gas[6–10].Another possibility for the creation of QDs was realized by means of the balance of the inter-component attraction and LHY repulsion in the binary boson mixture[11–13].

Theoretically,the LHY correction takes different forms in different spatial dimensions in the modified Gross–Pitaevskii equations.In the one-dimensional(1D)geometry,the competition of the quadratic LHY attraction and the effective cubic mean-field repulsion can form 1D QDs.QDs in 1D free space are studied and two different physical regimes are found,corresponding to small droplets of an approximately Gaussian shape and large puddles with a broad flat-top plateau,respectively[14].Small droplets collide quasielastically,while large colliding droplets may merge or suffer fragmentation,depending on their relative velocity.Optical lattice is an effective tool to deal with matter-wave solitons in Bose–Einstein condensates(BECs)and optical solitons in nonlinear media.Therefore,the introduction of optical lattice potential will influence the dynamics of QDs.It has been demonstrated that there are two kinds of physically regimes in optical lattices,namely,on-site QDs and off-site QDs[15].Both on-site QDs regardless of the value of norm and off-site QDs for the norm exceeding a critical value are completely stable.Unlike free space,the slowly moving small QDs in optical lattices merge after the collision.Multi-stable QDs are also studied in optical lattices[16].A continuous family of odd-symmetric and even-symmetric QDs exists in both first finite gap and semi-infinite gap.In addition,spontaneous symmetry breaking of QDs in a dual-core couplers was studied[17].QDs feature spatial density profiles of two different types:bell-shaped and flat-top ones,for relatively small and large values of N,respectively.The spontaneous symmetry breaking of QDs occurs and the symmetry is restored with the increase of the total norm.Collisions between moving QDs trend to merge into breathers.Very recently,the stability and collision ofPT-symmetric QDs have been studied in dual-core couplers with balanced gain and loss[18].It is found that the stability of symmetric QDs depends critically on the competition of gain and loss,intercore coupling,and optical lattice potential.Slowly moving PT-symmetric QDs tend to merge into breathers,while the fast-moving ones occur quasielastic collision and suffer fragmentation for small and large values of N,respectively.

Recently,QDs in 1D systems have been generalized to 2D and 3D by modifying the form of the LHY term[19–22].There are many potential values for stable QDs,such as manipulations of quantum information[23]and matter-wave interferometry[24].Therefore,it is of great significance to study the stability of QDs.On the other hand,the dynamics of BECs in optical lattices have been widely studied[25],and the study ofPT-symmetric optical lattices originates from quantum mechanics.For real eigenenergies and probability conservation,it is usually assumed that the Hamiltonian operator be Hermitian in quantum mechanics[26].However,it is proved that non-Hermitian Hamiltonian can also have entirely real spectra so long as they guaranteePT symmetry in the pioneering work[27].ThePT symmetry and non-Hermitian systems are in extensively investigated in the theory and experiments in recent decades.The asymmetric light transport was systematically investigated[28].The unidirectional absorber and unidirectional laser were proposed[29].In general,the parity operatoris defined asstand for momentum and position operators,respectively)and the time operatoris defined asAPT-symmetric optical lattice potential V(x)is realized by the condition V(x)=V*(−x).The most interesting feature of this pseudo-Hermitian Hamiltonian is the existence of a critical threshold beyond which a phase transition occurs due to spontaneousPT symmetry breaking[30–34].In this regime,the spectrum is no longer real.Because the equivalency between Schrödinger equation and the equation describing the propagation of light,optics[35–40]and BECs[41–44]inPT-symmetric systems have been reported in the past years.

In the present work,we study the stability and collisions of 1D QDs formed by Bose gas inPT-symmetric optical lattices.This paper is structured as follows.In section 2,we introduce the model,andPT-symmetric phase transition is revealed.In section 3,the stability of thePT-symmetric QDs is studied.ThePT-symmetric QDs exhibits oscillatory stability with the increase of N and g in the unbroken PT-symmetric phase.ThePT-symmetric QDs are entirely unstable in the brokenPT-symmetric phase.The collisions of QDs inPT-symmetric optical lattices are investigated in section 4.The collisions ofPT-symmetric QDs with equal norms are consistent with that of QDs in optical lattices.However,the collisions of droplets with unequal norms are completely different from that in free space.In addition,the collision between a stable and an unstablePT-symmetric QDs tends to merge into breathers after the collision.Finally,the paper is concluded by section 5.

2.The model

We investigate the system that QDs are formed in the binary condensate with mutually symmetric spinor components trapped in the 1DPT-symmetric optical lattices.The dynamics of the QDs in this system can be described by the modified Gross–Pitaevskii equations,including the LHY correction in the scaled form.

where g>0 represents the strength of the cubic self-repulsion,V(x)=V0sin2(x)andW(x)=W0sin (2x)are the real part and imaginary part of complex lattice potential,respectively.V0and W0satisfy thePT-symmetry,i.e.V(−x)=V(x),W(−x)=−W(x).V0and W0are the amplitudes of the real lattice and imaginary lattice,respectively.Stationary solutions of equation(1)are sought in the form ψ(x,t)=φ(x)e−iμt,where φ(x)is the complex function and μ is the chemical potential.In this case,φ(x)satisfies

to determine the linear stability properties of the QDs in the PT-symmetric optical lattices,we consider small perturbations on the solutions of equation(1)in the form,

where u and ν are the infinitesimal perturbations,the imaginary part of λ represents an unstable growth rate.The superscript*represents the complex conjugation.After substitution into equation(1)and linearizing,we obtain the following linear eigenvalue problem

where

Evidently,the QDs are linearly unstable if λ has an imaginary parts,but they are stable if λ is real.

The linear property of periodic potential can be understood by studying the linear problem of equation(2),i.e.Since the potentials V(x)and W(x)of equation(1)are π-periodic,according to Floquet Bloch theorem,the form of eigenfunction is φ=Φk(x)eikx,where Φk(x+π)=Φk(x)and k stands for the real Bloch momentum.In figure 1(a),we show the typical profiles of the real and imaginary parts ofPT-symmetric optical lattices.It is worth noting that the eigenvalues of the complex potential band structure are usually complex.However,forPT-symmetric periodic potential,as long as the parameter of system is lower than the phase transition point,the eigenvalues are completely real.Figure 1(b)shows the linear band gap structure with real part strength V0=0.3,and imaginary part strength W0=0.01,W0=0.15,and W0=0.3,respectively.Above the threshold value W0=0.15(phase transition point),the eigenvalue spectra become complex,and all bloch bands merge together concurrently.Figure 1(c)shows the imaginary parts of complex eigenvalues occurs at W0=0.3.There exists complex eigenvalues and QDs are unstable.In the broken PT-symmetric phase(W0>0.15),μ exists non-zero imaginary part.So,any QDs are unstable to perturbations and there is an overall instability effect on propagation.In the following,we discuss the dynamics of QDs in the unbroken PT-symmetric phase(W0＜0.15).

Figure 1.(a)Typical profile of PT-symmetric optical lattices,where black solid line represents real part,and red imaginary line describes imaginary part.(b)Band-gap structure of PT-symmetric optical lattices.Red solid line:V0=0.3,W0=0.01.Blue imaginary line:V0=0.3,W0=0.15.Green dotted–dashed line:V0=0.3,W0=0.3.(c)Imaginary part of complex eigenvalues corresponding to(b).

Figure 2.The profiles of PT-symmetric QDs in PT-symmetric optical lattices for different values of the condensate norm N=1(a)and N=10(b).The real and imaginary parts of PT symmetric optical lattices are V0=0.3 and W0=0.05,respectively.The other parameters are set as g=1.

3.Stability analysis for QDs in PT-symmetric optical lattices

In this section,we discuss the spatial profiles and stability of QDs inPT-symmetric optical lattices.A set of numerical solutions are obtained by using the the imaginary-time-integration method.Fixing the strength ofPT-symmetric optical lattices as V0=0.3 and W0=0.05,we depict the profiles of static droplets for different values of condensate normwhich is proportional to the number of atoms in the condensate,in figures 2(a)and(b).Borrowing the concept ofPT-symmetric solitons[39],we term the wave-function profiles of these self-bound states satisfying the property of φ*(x)=φ(−x)asPT-symmetric QDs.QDs in the system have two different types of spatial profiles,viz.,approximately Gaussian-shaped QDs[see figure 2(a)]and multi-humped QDs[see figure 2(b)]for relatively small and large values of N,respectively.

Below,the stability ofPT-symmetric QDs is considered.In reality,the system can not be completely isolated,and the dynamics of the system is expected to be robust to small fluctuations.Only if perturbed QDs can survive for a sufficiently long time,QDs can be observed in the experiment.Therefore,it is very important to study the stability of QDs.We can study stabilization of thePT-symmetric QDs through scanning a broad range of N and g.Figure 3(a)shows the linear stability analysis ofPT-symmetric QDs with different values of condensate norm N and the strength of the cubic self-repulsion g in the unbrokenPT-symmetric phase(W0=0.05).The QDs exhibit oscillating stability behavior for different strength of the cubic self-repulsion g and norm N.Fixing the parameter as g=1,we study stabilization of the PT-symmetric QDs through scanning a broad range of N in figure 3(d),corresponding to the dashed line of the longitudinal section in figure 3(a).The stability ofPT-symmetric QDs with different norm N have oscillating behavior.The linear stability analysis results may be also verified by direct numerical evolution simulations,which can be realized by using the split-step Fourier method with 3% random noise added into the initial conditions.The real-time dynamic evolution ofPT-symmetric QDs for N=1.2 and N=2.5 are depicted in figures 3(e)and(f),respectively,corresponding to the black pentagram marks of figure 3(d).Fixing the parameter as N=4 corresponding to the dashed line of the cross section of figure 3(a),we study the oscillating stability behavior of thePT-symmetric QDs through scanning a broad range of g in figure 3(g).The results can be also verified by the real-time dynamic evolution for g=0.57 and g=0.8,as shown in figures 3(h)and(i),corresponding to the black pentagram marks of figure 3(g).Moreover,it is worth noting that thePT-symmetric QDs become robustly stable at sufficiently large values of g.Figure 3(b)shows the linear stability analysis ofPT-symmetric QDs with different values of condensate norm N and the strength of the cubic self-repulsion g in the brokenPT-symmetric phase(W0=0.2).For the brokenPT-symmetric phase,the stable QDs can not be found regardless of the values of N and g.As predicted,QDs are unstable because the parameter of system exceeds the phase transition point and the band gap is closed and complex eigenvalues occur.In order to further clarify the effect of gain and loss on the stability ofPT-symmetric QDs,the phase diagram of the stability ofPT-symmetric QDs with respect to the intensity of imaginary part W0and the cubic self-repulsion g in figure 3(c)for N=4.As expected,thePT-symmetric QDs are entirely unstable in the brokenPT-symmetric phase(W0>0.15),and thePT-symmetric QDs exhibit oscillating stability behavior with the increase of g for fixing N and W0in the unbrokenPT-symmetric phase(W0＜0.15).The physical origin of oscillatory stability with the increase of N and g in the unbrokenPT-symmetric phase maybe is the competitive interplay between the QDs in optical lattice and its gain and loss,which must ensurePT symmetry of the system satisfied.These results are novel with respect to that of QDs in PT-symmetric dual-core couplers[18].

Figure 3.(a)A phase diagram of the stability of PT-symmetric QDs with respect to N and g in the unbroken PT-symmetric phase(W0=0.05).(b)A phase diagram of the stability of PT-symmetric QDs with respect to N and g in the broken PT-symmetric phase(W0=0.2).(c)Phase diagram of PT-symmetric QDs relative to W0 and g by fixing N=4.(d)The stabilityIm[λ]versus the condensate norm N corresponds to a longitudinal section of g=1 in(a).(e)and(f)are direct simulations of the perturbed evolution of PT-symmetric QDs of pentagram markers in(d),respectively.(g)The stabilityIm[λ]versus the self-repulsion strength g corresponds to cross section of N=4 in(a).(g)and(i)are direct simulations of the perturbed evolution of PT-symmetric QDs of pentagram markers in(g),respectively.

Figure 4.Typical examples of density plots for collisions between two PT-symmetric QDs in PT-symmetric optical lattices,launched as per equation(6)with x0=64.(a)N1=2,N2=8,k=0.1,g=1;(b)N1=2,N2=8,k=2,g=1;(c)N1=4,N2=5.5,k=0.1,g=1;(d)N1=4,N2=4,k=0.1,g1=0.57,g2=0.8.The other parameters are set as V0=0.3 and W0=0.05.

4.Collisions of PT-symmetric QDs

One of the characteristics of a soliton is that its shape does not change when it collides with another soliton.From this point of view,verifying a persistence of the shape of droplets involved in the pairwise collisions is very important problem.The collision dynamics of QDs in optical lattices have been proved that slow-moving ground state QDs tend to merger into breathers after collision,while fast-moving QDs may lead to quasi-elastic collision or splitting,which depends on the value of k.We simulated the collision and solved the initial condition of equation(1)

where ψ1(x)and ψ2(x)are the stationary shapes of QDs with normalization N1and N2,respectively.k is a kick that determines the velocity of droplet.This ansatz approximates a solution comprising two initial QDs located at x0and−x0.When W0is relatively small,the collision dynamics of QDs with equal norms is consistent with that in real optical lattices[15].For a relatively small value of k,the Gaussian-shaped QDs tend to merger into breathers after the collision.With the increase of k,the quasi-elastic collision is restored.For multihumped QDs,a newly merged quiescent breather occurs after the collision.For the fast-moving QDs,they undergo fragmentation and a majority of the particles keep in the moving ones.We also analyze collision dynamics of QDs with unequal norms in figures 4(a)and(b).In figure 4(a),a Gaussian-shaped droplet(N=2)collides with a multihumped droplet with N=8.It is found that the large droplet forms quiescent breather and the trajectory of the small droplet is deflected,and the two droplets become highly excited,exhibiting internal periodic vibrations.The situation is quite different for QDs in free space[14].With the increase of k,the two droplets undergo fragmentation after the collision.Unstable solitons can transform into breather-like objects[45].In figure 4(c),a stable droplet(N=5.5)collides with an unstable droplet with N=4.An interesting phenomenon occurs,two droplets tend to merger into breathers after the collision at a relatively small value of k.Besides,the collision between two droplets with N=4 corresponding to g=0.57 and g=0.8 is shown in figure 4(d),respectively.QDs form two breathers of unequal norms after the collision at a relatively small value of k.

5.Conclusion and discussion

We study the stability and collisions ofPT-symmetric QDs forming in a binary bosonic condensate trapped in PT-symmetric optical lattices.It is found that the stability of PT-symmetric QDs in thePT-symmetric system depends strongly on the values of the imaginary part W0of the PT-symmetric optical lattices,self-repulsion strength g,and the condensate norm N.Due to band gap closure and complex eigenvalues occur,thePT-symmetric QDs are entirely unstable in the brokenPT-symmetric phase(W0>0.15).The PT-symmetric QDs exhibit oscillatory stability with the increase of N and g.Finally,collisions between stable PT-symmetric QDs with unequal norms are systematically studied inPT-symmetric optical lattices.The slowly moving QDs tend to merge into breather after the collision.The quasielastic collision for fast-moving small QDs,while large droplets undergo fragmentation and the particles in colliding QDs are kept in the moving ones.However,for collisions of droplets with unequal norms,compared with a large puddle droplet collides with a small droplet,the large droplet becomes highly excited,while the small droplet remains essentially in an unperturbed shape in the free space,the moving slowly large and small stablePT-symmetric QDs in thePT-symmetric optical lattices are highly excited after the collision[14].With the increase of k,both of the droplets undergo fragmentation after the collision.Besides,the collision between a stable and an unstablePT-symmetric QDs,tend to merge into breathers after the collision.

Significantly,compared with conservative systems,the appearance ofPT symmetry concept has brought a variety of influences on theoretical and experimental researches in the past decade[46–49].The effect ofPT-symmetry on the stability of QDs is still lacking.Therefore,extending the concept ofPT symmetry to QDs can broaden more unknown secrets of QDs.

Acknowledgments

This work is supported by the National Natural Science Foundation of China under Grant No.11 805 283.The Hunan Provincial Natural Science Foundation under Grant No.2019JJ30044 and No.2019JJ40060.The Scientific Research Fund of Hunan Provincial Education Department of China under Grant No.19A510 and No.20B162.