Jiehong Lao¹, Zheng Zhou^,2,∗, Xili Zhang^,1, Fuqiu Ye^,3, Honghua Zhong^,41Department of Physics and Key Laboratory of Low-dimensional Quantum Structures and Quantum Control of Ministry of Education, Hunan Normal University, Changsha 410081, China
²Department of Physics, Hunan Institute of Technology, Hengyang 421002, China
³Department of Physics, Jishou University, Jishou 416000, China
⁴Institute of Mathematics and Physics, Central South University of Forestry and Technology, Changsha 410004, China

First author contact: ∗Author to whom any correspondence should be addressed.
Received:2021-02-11Revised:2021-03-20Accepted:2021-03-22Online:2021-04-20


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


PDF (0KB)MetadataMetricsRelated articlesExportEndNote|Ris|BibtexFavorite
Cite this article
Jiehong Lao, Zheng Zhou, Xili Zhang, Fuqiu Ye, Honghua Zhong. Oscillatory stability of quantum droplets in $\text{}{ \mathcal P }{ \mathcal T }$-symmetric optical lattice. Communications in Theoretical Physics, 2021, 73(6): 065103- doi:10.1088/1572-9494/abf093

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 of ³⁹K with appropriate signs of the inter-component 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 of ${ \mathcal P }{ \mathcal T }$-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, inter-core coupling, and optical lattice potential. Slowly moving ${ \mathcal P }{ \mathcal T }$-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 of ${ \mathcal P }{ \mathcal T }$-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 guarantee ${ \mathcal P }{ \mathcal T }$ symmetry in the pioneering work [27]. The ${ \mathcal P }{ \mathcal T }$ 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 operator $\hat{P}$ is defined as $\vec{p}\to -\vec{p}$, $\vec{x}\to -\vec{x}$ ($\vec{p}$, $\vec{x}$ stand for momentum and position operators, respectively) and the time operator $\hat{T}$ is defined as $\vec{p}\to -\vec{p}$, $\vec{x}\to \vec{x}$, $\vec{{\rm{i}}}\to -\vec{{\rm{i}}}$. A ${ \mathcal P }{ \mathcal T }$-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 spontaneous ${ \mathcal P }{ \mathcal T }$ 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] in ${ \mathcal P }{ \mathcal T }$-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 in ${ \mathcal P }{ \mathcal T }$-symmetric optical lattices. This paper is structured as follows. In section 2, we introduce the model, and ${ \mathcal P }{ \mathcal T }$-symmetric phase transition is revealed. In section 3, the stability of the ${ \mathcal P }{ \mathcal T }$-symmetric QDs is studied. The ${ \mathcal P }{ \mathcal T }$-symmetric QDs exhibits oscillatory stability with the increase of N and g in the unbroken ${ \mathcal P }{ \mathcal T }$-symmetric phase. The ${ \mathcal P }{ \mathcal T }$-symmetric QDs are entirely unstable in the broken ${ \mathcal P }{ \mathcal T }$-symmetric phase. The collisions of QDs in ${ \mathcal P }{ \mathcal T }$-symmetric optical lattices are investigated in section 4. The collisions of ${ \mathcal P }{ \mathcal T }$-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 unstable ${ \mathcal P }{ \mathcal T }$-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 1D ${ \mathcal P }{ \mathcal T }$-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.
(1)$\begin{eqnarray}{\rm{i}}{\partial }_{t}\psi =-\displaystyle \frac{1}{2}{\partial }_{{xx}}\psi +g| \psi {| }^{2}\psi -| \psi | \psi +[V(x)+{\rm{i}}W(x)]\psi ,\end{eqnarray}$
where g > 0 represents the strength of the cubic self-repulsion, $V(x)={V}_{0}{\sin }^{2}(x)$ and $W(x)={W}_{0}\sin (2x)$ are the real part and imaginary part of complex lattice potential, respectively. V₀ and W₀ satisfy the ${ \mathcal P }{ \mathcal T }$-symmetry, i.e. V(−x) = V(x), W(−x) = −W(x). V₀ and W₀ are 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$\mu$t, where φ(x) is the complex function and $\mu$ is the chemical potential. In this case, φ(x) satisfies
(2)$\begin{eqnarray}\displaystyle \frac{1}{2}{\partial }_{{xx}}\phi -g| \phi {| }^{2}\phi +| \phi | \phi -[V(x)+{\rm{i}}W(x)]\phi =\mu \phi ,\end{eqnarray}$
to determine the linear stability properties of the QDs in the ${ \mathcal P }{ \mathcal T }$-symmetric optical lattices, we consider small perturbations on the solutions of equation (1) in the form,
(3)$\begin{eqnarray}\psi =(\phi +u{{\rm{e}}}^{-{\rm{i}}\lambda t}+{\nu }^{* }{{\rm{e}}}^{{\rm{i}}{\lambda }^{* }t}){{\rm{e}}}^{-{\rm{i}}\mu t},\end{eqnarray}$
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
(4)$\begin{eqnarray}\left(\begin{array}{cc}{\hat{L}}_{1} & {\hat{L}}_{2}\\ -{\hat{L}}_{2}^{* } & -{\hat{L}}_{1}^{* }\end{array}\right)\left(\begin{array}{c}u\\ \nu \end{array}\right)=\lambda \left(\begin{array}{c}u\\ \nu \end{array}\right),\end{eqnarray}$
where
(5)$\begin{eqnarray}\begin{array}{rcl}{\hat{L}}_{1} & = & -\displaystyle \frac{1}{2}{\partial }_{{xx}}+2g| \phi {| }^{2}-\displaystyle \frac{3}{2}| \phi | +V(x)+{\rm{i}}W(x)-\mu ,\\ {\hat{L}}_{2} & = & g{\phi }^{2}-\displaystyle \frac{{\phi }^{2}}{2| \phi | }.\end{array}\end{eqnarray}$
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. $\tfrac{1}{2}{\partial }_{{xx}}\phi -[V(x)+{\rm{i}}W(x)]\phi =\mu \phi $. 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)e^ikx, 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 of ${ \mathcal P }{ \mathcal T }$-symmetric optical lattices. It is worth noting that the eigenvalues of the complex potential band structure are usually complex. However, for ${ \mathcal P }{ \mathcal T }$-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 V₀ = 0.3, and imaginary part strength W₀ = 0.01, W₀ = 0.15, and W₀ = 0.3, respectively. Above the threshold value W₀ = 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 W₀ = 0.3. There exists complex eigenvalues and QDs are unstable. In the broken ${ \mathcal P }{ \mathcal T }$-symmetric phase (W₀ > 0.15), $\mu$ 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 ${ \mathcal P }{ \mathcal T }$-symmetric phase (W₀ < 0.15).

Figure 1.

New window|Download| PPT slide
Figure 1.(a) Typical profile of ${ \mathcal P }{ \mathcal T }$-symmetric optical lattices, where black solid line represents real part, and red imaginary line describes imaginary part. (b) Band-gap structure of ${ \mathcal P }{ \mathcal T }$-symmetric optical lattices. Red solid line: V₀ = 0.3, W₀ = 0.01. Blue imaginary line: V₀ = 0.3, W₀ = 0.15. Green dotted-dashed line: V₀ = 0.3, W₀ = 0.3. (c) Imaginary part of complex eigenvalues corresponding to (b).

3. Stability analysis for QDs in ${ \mathcal P }{ \mathcal T }$-symmetric optical lattices

In this section, we discuss the spatial profiles and stability of QDs in ${ \mathcal P }{ \mathcal T }$-symmetric optical lattices. A set of numerical solutions are obtained by using the the imaginary-time-integration method. Fixing the strength of ${ \mathcal P }{ \mathcal T }$-symmetric optical lattices as V₀ = 0.3 and W₀ = 0.05, we depict the profiles of static droplets for different values of condensate norm $N={\int }_{-\infty }^{\infty }| \psi {| }^{2}{\rm{d}}x$, which is proportional to the number of atoms in the condensate, in figures 2(a) and (b). Borrowing the concept of ${ \mathcal P }{ \mathcal T }$-symmetric solitons [39], we term the wave-function profiles of these self-bound states satisfying the property of φ^*(x) = φ(−x) as ${ \mathcal P }{ \mathcal T }$-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.

Figure 2.

New window|Download| PPT slide
Figure 2.The profiles of ${ \mathcal P }{ \mathcal T }$-symmetric QDs in ${ \mathcal P }{ \mathcal T }$-symmetric optical lattices for different values of the condensate norm N = 1 (a) and N = 10 (b). The real and imaginary parts of ${ \mathcal P }{ \mathcal T }$ symmetric optical lattices are V₀ = 0.3 and W₀ = 0.05, respectively. The other parameters are set as g = 1.


Below, the stability of ${ \mathcal P }{ \mathcal T }$-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 the ${ \mathcal P }{ \mathcal T }$-symmetric QDs through scanning a broad range of N and g. Figure 3(a) shows the linear stability analysis of ${ \mathcal P }{ \mathcal T }$-symmetric QDs with different values of condensate norm N and the strength of the cubic self-repulsion g in the unbroken ${ \mathcal P }{ \mathcal T }$-symmetric phase (W₀ = 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 ${ \mathcal P }{ \mathcal T }$-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 of ${ \mathcal P }{ \mathcal T }$-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 of ${ \mathcal P }{ \mathcal T }$-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 the ${ \mathcal P }{ \mathcal T }$-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 the ${ \mathcal P }{ \mathcal T }$-symmetric QDs become robustly stable at sufficiently large values of g. Figure 3(b) shows the linear stability analysis of ${ \mathcal P }{ \mathcal T }$-symmetric QDs with different values of condensate norm N and the strength of the cubic self-repulsion g in the broken ${ \mathcal P }{ \mathcal T }$-symmetric phase (W₀ = 0.2). For the broken ${ \mathcal P }{ \mathcal T }$-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 of ${ \mathcal P }{ \mathcal T }$-symmetric QDs, the phase diagram of the stability of ${ \mathcal P }{ \mathcal T }$-symmetric QDs with respect to the intensity of imaginary part W₀ and the cubic self-repulsion g in figure 3(c) for N = 4. As expected, the ${ \mathcal P }{ \mathcal T }$-symmetric QDs are entirely unstable in the broken ${ \mathcal P }{ \mathcal T }$-symmetric phase (W₀ > 0.15), and the ${ \mathcal P }{ \mathcal T }$-symmetric QDs exhibit oscillating stability behavior with the increase of g for fixing N and W₀ in the unbroken ${ \mathcal P }{ \mathcal T }$-symmetric phase (W₀ < 0.15). The physical origin of oscillatory stability with the increase of N and g in the unbroken ${ \mathcal P }{ \mathcal T }$-symmetric phase maybe is the competitive interplay between the QDs in optical lattice and its gain and loss, which must ensure ${ \mathcal P }{ \mathcal T }$ symmetry of the system satisfied. These results are novel with respect to that of QDs in ${ \mathcal P }{ \mathcal T }$-symmetric dual-core couplers [18].

Figure 3.

New window|Download| PPT slide
Figure 3.(a) A phase diagram of the stability of ${ \mathcal P }{ \mathcal T }$-symmetric QDs with respect to N and g in the unbroken ${ \mathcal P }{ \mathcal T }$-symmetric phase (W₀ = 0.05). (b) A phase diagram of the stability of ${ \mathcal P }{ \mathcal T }$-symmetric QDs with respect to N and g in the broken ${ \mathcal P }{ \mathcal T }$-symmetric phase (W₀ = 0.2). (c) Phase diagram of ${ \mathcal P }{ \mathcal T }$-symmetric QDs relative to W₀ and g by fixing N = 4. (d) The stability $\mathrm{Im}[\lambda ]$ 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 ${ \mathcal P }{ \mathcal T }$-symmetric QDs of pentagram markers in (d), respectively. (g) The stability $\mathrm{Im}[\lambda ]$ 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 ${ \mathcal P }{ \mathcal T }$-symmetric QDs of pentagram markers in (g), respectively.

4. Collisions of ${ \mathcal P }{ \mathcal T }$-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)
(6)$\begin{eqnarray}\psi (x,t=0)={{\rm{e}}}^{{\rm{i}}{kx}}{\psi }_{1}(x+{x}_{0})+{{\rm{e}}}^{-{\rm{i}}{kx}}{\psi }_{2}(x-{x}_{0}),\end{eqnarray}$
where ψ₁(x) and ψ₂(x) are the stationary shapes of QDs with normalization N₁ and N₂, respectively. k is a kick that determines the velocity of droplet. This ansatz approximates a solution comprising two initial QDs located at x₀ and −x₀. When W₀ is 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 multi-humped 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 multi-humped 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.

Figure 4.

New window|Download| PPT slide
Figure 4.Typical examples of density plots for collisions between two ${ \mathcal P }{ \mathcal T }$-symmetric QDs in ${ \mathcal P }{ \mathcal T }$-symmetric optical lattices, launched as per equation (6) with x₀ = 64. (a) N₁ = 2, N₂ = 8, k = 0.1, g = 1; (b) N₁ = 2, N₂ = 8, k = 2, g = 1; (c) N₁ = 4, N₂ = 5.5, k = 0.1, g = 1; (d) N₁ = 4, N₂ = 4, k = 0.1, g₁ = 0.57, g₂ = 0.8. The other parameters are set as V₀ = 0.3 and W₀ = 0.05.

5. Conclusion and discussion

We study the stability and collisions of ${ \mathcal P }{ \mathcal T }$-symmetric QDs forming in a binary bosonic condensate trapped in ${ \mathcal P }{ \mathcal T }$-symmetric optical lattices. It is found that the stability of ${ \mathcal P }{ \mathcal T }$-symmetric QDs in the ${ \mathcal P }{ \mathcal T }$-symmetric system depends strongly on the values of the imaginary part W₀ of the ${ \mathcal P }{ \mathcal T }$-symmetric optical lattices, self-repulsion strength g, and the condensate norm N. Due to band gap closure and complex eigenvalues occur, the ${ \mathcal P }{ \mathcal T }$-symmetric QDs are entirely unstable in the broken ${ \mathcal P }{ \mathcal T }$-symmetric phase (W₀ > 0.15). The ${ \mathcal P }{ \mathcal T }$-symmetric QDs exhibit oscillatory stability with the increase of N and g. Finally, collisions between stable ${ \mathcal P }{ \mathcal T }$-symmetric QDs with unequal norms are systematically studied in ${ \mathcal P }{ \mathcal T }$-symmetric optical lattices. The slowly moving QDs tend to merge into breather after the collision. The quasi-elastic 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 stable ${ \mathcal P }{ \mathcal T }$-symmetric QDs in the ${ \mathcal P }{ \mathcal T }$-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 unstable ${ \mathcal P }{ \mathcal T }$-symmetric QDs, tend to merge into breathers after the collision.

Significantly, compared with conservative systems, the appearance of ${ \mathcal P }{ \mathcal T }$ symmetry concept has brought a variety of influences on theoretical and experimental researches in the past decade [46-49]. The effect of ${ \mathcal P }{ \mathcal T }$-symmetry on the stability of QDs is still lacking. Therefore, extending the concept of ${ \mathcal P }{ \mathcal T }$ 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.

Reference View Option By original order
By published year
By cited within times
By Impact factor

[1]Petrov D S 2015Quantum mechanical stabilization of a collapsing Bose−Bose mixture
Phys. Rev. Lett 115 155302
DOI:10.1103/PhysRevLett.115.155302 [Cited within: 1]

[2]Chomaz L, Baier S, Petter D, Mark M J, Wächtler F, Santos L, Ferlaino F 2016Quantum-fluctuation-driven crossover from a dilute Bose−Einstein condensate to a macrodroplet in a dipolar quantum fluid
Phys. Rev. X 6 041039
DOI:10.1103/PhysRevX.6.041039 [Cited within: 1]

[3]Schmitt M, Wenzel M, Böttcher F, Ferrier-Barbut I, Pfau T 2016Self-bound droplets of a dilute magnetic quantum liquid
Nature 539 259
DOI:10.1038/nature20126 [Cited within: 1]

[4]Cabrera C R, Tanzi L, Sanz J, Naylor B, Thomas P, Cheiney P, Tarruell L 2018Quantum liquid droplets in a mixture of Bose−Einstein condensates
Science 359 301
DOI:10.1126/science.aao5686 [Cited within: 1]

[5]Lee D T, Huang K, Yang N C 1957Eigenvalues and eigenfunctions of a bose system of hard spheres and its low-temperature properties
Phys. Rev. 106 1135
DOI:10.1103/PhysRev.106.1135 [Cited within: 1]

[6]Wächtler F, Santos L 2016Ground-state properties and elementary excitations of quantum droplets in dipolar Bose−Einstein condensates
Phys. Rev. A 94 043618
DOI:10.1103/PhysRevA.94.043618 [Cited within: 1]

[7]Baillie D, Wilson R M, Bisset R N, Blakie P B 2016Self-bound dipolar droplet: a localized matter wave in free space
Phys. Rev. A 94 021602(R)
DOI:10.1103/PhysRevA.94.021602

[8]Ferrier-Barbut I, Kadau H, Schmitt M, Wenzel M, Pfau T 2016Observation of quantum droplets in a strongly dipolar bose gas
Phys. Rev. Lett. 116 215301
DOI:10.1103/PhysRevLett.116.215301

[9]Ferrier-Barbut I, Wenze M, Bötcher F, Langen T, Isoard M, Stringari S, Pfau T 2018Scissors mode of dipolar quantum droplets of dysprosium atoms
Phys. Rev. Lett. 120 160402
DOI:10.1103/PhysRevLett.120.160402

[10]Cidrim A, dos Santos F E A, Henn E A L, Macrì T 2018Vortices in self-bound dipolar droplets
Phys. Rev. A 98 023618
DOI:10.1103/PhysRevA.98.023618 [Cited within: 1]

[11]Staudinger C, Mazzanti F, Zillich R E 2018Self-bound Bose mixtures
Phys. Rev. A 98 023633
DOI:10.1103/PhysRevA.98.023633 [Cited within: 1]

[12]Cheiney P, Cabrera C R, Sanz J, Naylor B, Tanzi L, Tarruell L 2018Bright soliton to quantum droplet transition in a mixture of Bose−Einstein condensates
Phys. Rev. Lett. 120 135301
DOI:10.1103/PhysRevLett.120.135301

[13]Semeghini G, Ferioli G, Masi L, Mazzinghi C, Wolswijk L, Minardi F, Modugno M, Modugno G, Inguscio M, Fattori M 2018Self-bound quantum droplets of atomic mixtures in free space
Phys. Rev. Lett. 120 235301
DOI:10.1103/PhysRevLett.120.235301 [Cited within: 1]

[14]Astrakharchik G E, Malomed B A 2018Dynamics of one-dimensional quantum droplets
Phys. Rev. A 98 013631
DOI:10.1103/PhysRevA.98.013631 [Cited within: 3]

[15]Zhou Z, Yu X, Zou Y, Zhong H H 2019Dynamics of quantum droplets in a one-dimensional optical lattice
Commun. Nonlinear Sci. Numer. Simul. 78 104881
DOI:10.1016/j.cnsns.2019.104881 [Cited within: 2]

[16]Dong L W, Qi W, Peng P, Wang L X, Zhou H, Huang C M 2020Multi-stable quantum droplets in optical lattices
Nonlinear Dyn. 102 303 310
DOI:10.1007/s11071-020-05967-y [Cited within: 1]

[17]Liu B, Zhang H F, Zhong R X, Zhang X L, Qin X Z, Huang C Q, Li Y Y, Malomed B A 2019Symmetry breaking of quantum droplets in a dual-core trap
Phys. Rev. A 99 053602
DOI:10.1103/PhysRevA.99.053602 [Cited within: 1]

[18]Zhou Z, Zhu B, Wang H B, Zhong H H 2020Stability and collisions of quantum droplets in ${ \mathcal P }{ \mathcal T }$-symmetric dual-core couplers
Commun. Nonlinear Sci. Numer. Simul. 91 105424
DOI:10.1016/j.cnsns.2020.105424 [Cited within: 2]

[19]Petrov D S, Astrakharchik G 2016Ultradilute low-dimensional liquids
Phys. Rev. Lett. 117 100401
DOI:10.1103/PhysRevLett.117.100401 [Cited within: 1]

[20]Li Y Y, Chen Z P, Luo Z H, Huang C Q, Tan H S, Pang W, Malomed B A 2018Two-dimensional vortex quantum droplets
Phys. Rev. A 98 063602
DOI:10.1103/PhysRevA.98.063602

[21]Kartashov Y V, Malomed B A, Torner L 2019Metastability of quantum droplet clusters
Phys. Rev. Lett. 122 193902
DOI:10.1103/PhysRevLett.122.193902

[22]Ferioli G, Semeghini G, Masi L, Giusti G, Modugno G, Inguscio M, Gallemí A, Recati A, Fattori M 2019Collisions of self-bound quantum droplets
Phys. Rev. Lett. 122 090401
DOI:10.1103/PhysRevLett.122.090401 [Cited within: 1]

[23]Bullough R K, Wadati M 2004Optical solitons and quantum solitons
J. Opt. B: Quantum Semiclass. Opt. 6 S205
DOI:10.1088/1464-4266/6/5/007 [Cited within: 1]

[24]McDonald G D, Kuhn C C N, Hardman K S, Bennetts S, Everitt P J, Atlin P A, Debs J E, Close J D, Robins N P 2014Bright solitonic matter-wave interferometer
Phys. Rev. Lett. 113 013002
DOI:10.1103/PhysRevLett.113.013002 [Cited within: 1]

[25]Witthaut D, Trimborn F, Hennig H, Kordas G, Geisel T, Wimberger S 2011Beyond mean-field dynamics in open Bose−Hubbard chains
Phys. Rev. A 83 063608
DOI:10.1103/PhysRevA.83.063608 [Cited within: 1]

[26]Messiah A 1999 Quantum Mechanics New York Dover
[Cited within: 1]

[27]Bender C M, Boettcher S 1998Real spectra in non-Hermitian Hamiltonians having ${ \mathcal P }{ \mathcal T }$ symmetry
Phys. Rev. Lett. 80 5243
DOI:10.1103/PhysRevLett.80.5243 [Cited within: 1]

[28]Jin L, Song Z 2021Symmetry-protected scattering in non-hermitian linear systems
Chin. Phys. Lett. 38 024202
DOI:10.1088/0256-307X/38/2/024202 [Cited within: 1]

[29]Jin L, Song Z 2018Incident direction independent wave propagation and unidirectional lasing
Phys. Rev. Lett. 121 073901
DOI:10.1103/PhysRevLett.121.073901 [Cited within: 1]

[30]Bender C M, Dunne G V, Meisinger P N 1999Complex periodic potentials with real band spectra
Phys. Lett. A 252 272
DOI:10.1016/S0375-9601(98)00960-8 [Cited within: 1]

[31]El-Ganainy R, Makris K G, Christodoulides D N, Musslimani Z H 2007Theory of coupled optical ${ \mathcal P }{ \mathcal T }$ symmetric structures
Opt. Lett. 32 2632
DOI:10.1364/OL.32.002632

[32]Jin L, Song Z 2010Wave packet dynamics in a non-Hermitian ${ \mathcal P }{ \mathcal T }$ symmetric tight-binding chain
Commun. Theor. Phys. 54 73
DOI:10.1088/0253-6102/54/1/14

[33]Jin L, Song Z 2010Scaling behavior and phase diagram of a ${ \mathcal P }{ \mathcal T }$-symmetric non-hermitian Bose−Hubbard system
Ann. Phys. 330 013629


[34]Graefe E, Korsch H J, Niederle A E 2010Quantum-classical correspondence for a non-hermitian Bose−Hubbard dimer
Phys. Rev. A 82 013629
DOI:10.1103/PhysRevA.82.013629 [Cited within: 1]

[35]Makris K G, El-Ganainy R, Christodoulides D N, Musslimani Z H 2008Beam dynamics in ${ \mathcal P }{ \mathcal T }$ symmetric optical lattices
Phys. Rev. Lett. 100 103904
DOI:10.1103/PhysRevLett.100.103904 [Cited within: 1]

[36]Musslimani Z H, Makris K G, El-Ganainy R, Christodoulides D N 2008Optical solitons in ${ \mathcal P }{ \mathcal T }$ periodic potentials
Phys. Rev. Lett. 100 030402
DOI:10.1103/PhysRevLett.100.030402

[37]Guo A, Salamo G J, Duchesne D, Morandotti R, Ravat M V, Aimez V, Siviloglou G A, Christodoulides D N 2009Observation of ${ \mathcal P }{ \mathcal T }$-symmetry breaking in complex optical potentials
Phys. Rev. Lett. 103 093902
DOI:10.1103/PhysRevLett.103.093902

[38]Makris K G, El-ganainy R, Christodoulides D N, Musslimani Z H 2010${ \mathcal P }{ \mathcal T }$-symmetric optical lattices
Phys. Rev. A 81 063807
DOI:10.1103/PhysRevA.81.063807

[39]Nixon S, Ge L, Yang J 2012Stability analysis for solitons in ${ \mathcal P }{ \mathcal T }$-symmetric optical lattices
Phys. Rev. A 85 023822
DOI:10.1103/PhysRevA.85.023822 [Cited within: 1]

[40]Li K, Kevrekidis P G 2011${ \mathcal P }{ \mathcal T }$-symmetric oligomers: Analytical solutions, linear stability, and nonlinear dynamics
Phys. Rev. E 83 066608
DOI:10.1103/PhysRevE.83.066608 [Cited within: 1]

[41]Graefe E M 2012Stationary states of a ${ \mathcal P }{ \mathcal T }$ symmetric two-mode Bose−Einstein condensate
J. Phys. A: Math. Theor 45 444015
DOI:10.1088/1751-8113/45/44/444015 [Cited within: 1]

[42]Haag D, Dast D, Cartarius H, Wunner G 2018${ \mathcal P }{ \mathcal T }$-symmetric gain and loss in a rotating Bose-Einstein condensate
Phys. Rev. A 97 033607
DOI:10.1103/PhysRevA.97.033607

[43]Dizdarevic D, Dast D, Haag D, Main J, Cartarius H, Wunner G 2015Cusp bifurcation in the eigenvalue spectrum of ${ \mathcal P }{ \mathcal T }$-symmetric Bose−Einstein condensates
Phys. Rev. A 91 033636
DOI:10.1103/PhysRevA.91.033636

[44]Dast D, Haag D, Cartarius H, Main J, W unner G 2016Bose−Einstein condensates with balanced gain and loss beyond mean-field theory
Phys. Rev. A 94 053601
DOI:10.1103/PhysRevA.94.053601 [Cited within: 1]

[45]Barashenkov I V, Suchkov S V, Sukhorukov A A, Dmitriev S V, Kivshar Y S 2012Breathers in ${ \mathcal P }{ \mathcal T }$-symmetric optical couplers
Phys. Rev. A 86 053809
DOI:10.1103/PhysRevA.86.053809 [Cited within: 1]

[46]Miri M A, Aceves A B, Kottos T, Kovanis V, Christodoulidesr D N 2012Bragg solitons in nonlinear ${ \mathcal P }{ \mathcal T }$-symmetric periodic potentials
Phys. Rev. A 86 033801
DOI:10.1103/PhysRevA.86.033801 [Cited within: 1]

[47]Longhi S 2010${ \mathcal P }{ \mathcal T }$-symmetric laser absorber
Phys. Rev. A 82 031801
DOI:10.1103/PhysRevA.82.031801

[48]Bender N, Factor S, Bodyfelt J D, Ramezani H, Christodoulides D N, Ellis F M, Kottos T 2011Experimental study of active LRC circuits with ${ \mathcal P }{ \mathcal T }$-symmetries
Phys. Rev. A 84 040101
DOI:10.1103/PhysRevA.84.040101

[49]Bender N, Factor S, Bodyfelt J D, Ramezani H, Christodoulides D N, Ellis F M, Kottos T 2013Observation of asymmetric transport in structures with active nonlinearities
Phys. Rev. A 110 234101
DOI:10.1103/PhysRevLett.110.234101 [Cited within: 1]