Next Article in Journal
The Novel Scaling of Tsallis Parameters from the Transverse Momentum Spectra of Charged Particles in Heavy-Ion Collisions
Previous Article in Journal
Effect of Protoneutron Star Magnetized Envelops in Neutrino Energy Spectra
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Rabi Frequency Management of Collapsing Quasi-Two-Dimensional Bose-Einstein Condensates with Pseudospin-1/2

1
Ulugh Beg Astronomical Institute, Astronomy St. 33, Tashkent 100052, Uzbekistan
2
Tashkent Institute of Irrigation and Agricultural Mechanization Engineering, Kori Niyozi 39, Tashkent 100000, Uzbekistan
3
Department of Economics, Sustainable Agriculture and Digital Technology, Samarkand Branch of Tashkent State Agrarian University, A.Temur 7, Okdaryo, Samarkand 191200, Uzbekistan
4
Physics Faculty, National University of Uzbekistan, Tashkent 100174, Uzbekistan
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 31 March 2022 / Revised: 25 April 2022 / Accepted: 26 April 2022 / Published: 28 April 2022

Abstract

:
The collapse of quasi-two-dimensional pseudospin-1/2 Bose-Einstein condensate of attracting atoms with intra- and cross-spin interaction is studied in the presence of the Rabi coupling. The condensate dynamics is presented as a function of the self-interaction and Rabi frequency. The evolution of two components of the condensate by using the Gross-Pitaevskii equations is investigated. The initial Gaussian ansatz for two-component wave functions is selected for the better interpretation of the numerical results. The intra-spin-coupling modifies the critical number of atoms causing the collapse while the collapse is observed only in a single pseudospin component. It is demonstrated that for cross-spin-coupling only double spin-components collapse can occur.

1. Introduction

Understanding of behavior of nonlinear quantum systems, including solitons and instantons, is critically important for understanding of quantum field and particles physics [1]. Thus, studies of experimentally available nonlinear quantum condensed matter systems can shed light on variety of processes occuring on much smaller spatial scales and involving much higher energies. One of those nonlinear quantum condensed matter systems is Bose-Einstein condensate (BEC) of interacting atoms [2]. The BECs open new research avenues in low-energy quantum physics as well as in astrophysics and cosmology, where they are actively studied for the understanding of the nature of mysterious Dark Matter [3] and possible alternatives of the relativistic neutron stars as boson stars. The nonlinear self-interaction that can be “repulsive” or “attractive” due to the Feshbach resonance [4] or dipolar interactions [5,6] is critically important for the BECs properties. The attractive interaction in the BEC can produce solitons [7], quantum droplets [8,9], collapse processes [10], and many other phenomena.
The BEC collapse is a dynamical process, corresponding to the squeezing of the characteristic size of the condensate into a point and sequentially an explosion as occurs in the experiments [11,12,13,14]. The collapse process depends on the dimension of the system and interatomic interaction. For instance, in a one-dimensional system, a cubic interatomic attraction in the condensate characterizes soliton dynamics [15] rather than the collapse. However, in the presence of the stronger quintic nonlinearity, the condensate can collapse [16]. In three dimensional systems with cubic nonlinearity attractive interaction always leads to the BEC collapse while in two dimensions a critical number of atoms is needed for the collapse to occur. The stabilization of the collapse opens a wide range of possible features of self-attracting BEC with the studies carried out by introducing a variety of obstacles into the self-attraction of the condensate [17,18,19,20,21,22,23].
In this paper we study the quasi-two-dimensional BEC of atoms with pseudospin-1/2. Studying the cubic nonlinearity, the main focus will be on two types of interactions: intra- and cross-spin coupling of components [24,25,26]. To produce the spin dynamics, we use an external synthetic Rabi magnetic field which leads to redistribution of the condensate density between the spin components. As a result, we will demonstrate the collapse of the condensate depending on the type of interactions: (i) with the intra-spin interaction the collapse can be observed only in one spin component depending on the Rabi field strength; (ii) with the cross-spin interaction the collapse can be observed only in both spin components for a sufficient strength of the Rabi field.
Controlled collapse of two and three-dimensional condensate is demonstrated in [17] through changing scattering length of atoms, that is in the framework of time-dependent inter-atomic interaction g. As a result, the stability of the condensate is limited in time, eventually switching into collapse. Here, we apply a Rabi field to rotate the BEC pseudospin and the spin rotation leads to the oscillation of condensate density between spin components consequently modifying inter-atomic interaction. Thus we obtain the effect being similar to the oscillated interatomic scattering length used in Ref. [17]. Despite this similarity, the physics of the problem studied here is considerably richer since, as it will be shown, various forms of collapse and collapse inhibition can occur here.

2. The Model and Main Parameters

The collapse dynamics of quasi-two-dimensional BEC with pseudospin-1/2 coupled by the synthetic Rabi magnetic field is characterized by evolution of the wave functions ψ i ψ i ( r , t ) ( i = 1 , 2 ) with the total norm N, r ( x , y ) . The evolution of the wave functions is described by the Gross-Pitaevskii equations (GPEs) in the form:
i t ψ 1 = 2 2 m Δ ψ 1 g 1 + | κ | ( 1 κ ) | ψ 1 | 2 + ( 1 + κ ) | ψ 2 | 2 ψ 1 + Ω 2 ψ 2 , i t ψ 2 = 2 2 m Δ ψ 2 g 1 + | κ | ( 1 κ ) | ψ 2 | 2 + ( 1 + κ ) | ψ 1 | 2 ψ 2 + Ω 2 ψ 1 .
Here m is the atomic mass, the interaction constant g = 4 π 2 | a s | / m a z , a z is the condensate extension along the z axis, and a s < 0 is the s-wave scattering length. The Rabi frequency Ω corresponds to the effective Zeeman magnetic field directed along the x-axis. In the system of GPEs (1) κ = 1 , 0 , 1 describes the intra-, total and cross-spin interaction, correspondingly. The spin components are normalized by
N i = 2 π 0 ψ i ( r , t ) 2 r d r ,
where ( i = 1 , 2 ) , N i N i ( t ) and the total norm N = N 1 + N 2 is conserved.
Since the dynamics of the collapse is characterized by the density profile of the condensate components, we study the evolution of widths of the spin components and spin dynamics driven by Rabi field for given inter-atomic interactions. The evolution of the width of the spin components is characterized by time-dependent inverse participation ratio in the form:
a i ( t ) = N i 2 π 0 ψ i ( r , t ) 4 r d r 1 / 2 .
The nonlinear condensate dynamics has a strong impact on the evolution of spin components defined by
σ j ( t ) = 2 π N 0 { ψ 1 * , ψ 2 * } T σ j { ψ 1 , ψ 2 } r d r ,
where σ j ( j = x , y , z ) are corresponding Pauli matrices, “*” stands for the complex conjugate and “ T ” is the transposition. Consequently, the length of the spin vector P ( t ) [27]
P ( t ) = j σ j ( t ) 2 1 ,
corresponding to the decreased purity of the condensate in the spin subspace.
From Equation (1) the total energy of the system is given by
E T κ = E K + E I κ + E R .
Here E K , E I κ , and E R are the kinetic, interaction, and Rabi magnetic field energies, respectively, determined by
E K = 2 π 0 2 2 m ψ 1 * Δ ψ 1 + ψ 2 * Δ ψ 2 r d r ,
E I κ = π g 1 + | κ | 0 ( 1 κ ) | ψ 1 | 4 + | ψ 2 | 4 + 2 ( 1 + κ ) | ψ 1 | 2 | ψ 2 | 2 r d r ,
E R = π Ω 0 ψ 1 * ψ 2 + ψ 2 * ψ 1 r d r .
From Equations (4) and (9) follows, the Rabi magnetic field energy can be characterized with spin parameter σ x by E R = N σ x / 2 .
For the numerical solution of Equation (1) we assume that at t = 0 one has g = 0 and Ω = 0 , and the condensate is prepared in the ground state of harmonic potential. At t > 0 the harmonic potential is suddenly switched off and simultaneously the interaction and the Rabi magnetic field are switched on causing the BEC dynamics. Thus, we start with the Gaussian ground state wave function with the initial spin along the z-axis
ψ 1 ( r , t = 0 ) = N π a 0 2 exp r 2 2 a 0 2 ,
where a 0 is initial width of the condensate and ψ 2 ( r , t = 0 ) = 0 . Hereafter we use the units with m = 1 and consequently dimensionless g = 4 π | a s | / a z , unit of length is arbitrary and unit of time is 2 .
To have a reference point we assume that the condensate is prepared in the ground state form (10) and in Equation (1) Ω = 0 , κ = 0 . Thus, Equation (6) yields
E 0 = N 2 π Λ a 0 2 ,
where Λ = ( g N λ cr ) / 2 and λ cr = 2 π is the critical interaction parameter g N to obtain the collapse as characterized by condition E 0 for a 0 satisfied only if Λ > 0 [28]. In this case the corresponding time dependent width is given by [29]
a v ( t ) = a 0 1 Λ t 2 a 0 4 ,
with the collapse time being T c = a 0 2 / Λ .
For a qualitative explanation of the Rabi magnetic field effect, we use Gaussian ansatz for two-component wave functions in the form,
ψ 1 ( r , t ) = N π a ˜ 1 2 ( t ) 1 / 2 exp r 2 2 a ˜ 1 2 ( t ) cos Ω 2 t , ψ 2 ( r , t ) = N π a ˜ 2 2 ( t ) 1 / 2 exp r 2 2 a ˜ 2 2 ( t ) sin Ω 2 t .
Here, we assume the a ˜ i 2 ( t ) is a complex number. The real part Re [ a ˜ i 2 ( t ) ] > 0 in order to have a wave function of Gaussian form with a ˜ i 2 ( t ) = Re [ a ˜ i 2 ( t ) ] .
In the wave functions (13), Ω is the spin precession rate corresponding to oscillations of the condensate density between spin components. From the ansatz functions (13) we can characterize widths of spin components by
2 π N 0 r 3 | ψ 1 ( r , t ) | 2 d r 1 / 2 = | a ˜ 1 ( t ) | 2 a ˜ 1 ( t ) cos Ω 2 t , 2 π N 0 r 3 | ψ 2 ( r , t ) | 2 d r 1 / 2 = | a ˜ 2 ( t ) | 2 a ˜ 2 ( t ) sin Ω 2 t .
Here the widths of the spin-projected wave packets are characterized by | a ˜ 1 ( t ) | 2 / a ˜ 1 ( t ) where factors with cosine and sine functions are explicitly related to the spin rotation. If in Equation (14) Im [ a ˜ i 2 ( t ) ] = 0 , then the spin components widths are characterized by a ˜ i ( t ) .
Since the kinetic and Rabi-field related energies do not explicitly depend on the self-interaction type, then for ansatz functions (13) the Equations (7) and (9) become:
E K = N 2 1 a ˜ 1 2 ( t ) cos 2 Ω 2 t + 1 a ˜ 2 2 ( t ) sin 2 Ω 2 t ,
E R = Ω N a ˜ 1 ( t ) a ˜ 2 ( t ) | a ˜ 1 ( t ) a ˜ 2 ( t ) | 2 a ˜ 1 2 ( t ) | a ˜ 2 ( t ) | 4 + a ˜ 2 2 ( t ) | a ˜ 1 ( t ) | 4 | a ˜ 1 2 ( t ) + a ˜ 2 * 2 ( t ) | 2 sin Ω t .
Consequently, as can be seen from (16), the Rabi energy can be positive or negative, in order to compensate kinetic and self-interaction energies depending on the spin state.

3. Numerical Results and Discussion

3.1. Intra-Spin Interaction

In this section we discuss the dynamics of the intra-spin interacting condensate with κ = 1 in Equation (1). Since at time t = 0 , the initial spin (10) is directed along the z-axis, at time t > 0 the applied Rabi magnetic field starts to rotate the spin and the wave function of the second spin component became nonzero. Correspondingly, the density of the condensate included in (1) splits between two spin components, thus reducing the interaction between atoms because of the absence of spin cross-coupling term. As a result, although for a weak Rabi field, the attractive self-interaction will be reduced by reorientation of spin and density redistribution, the condensate eventually collapses. However, for a sufficiently strong Rabi field, spin rotation provides oscillations of interaction on the 2 π / Ω timescale and can prevent the collapse, similar to the effect studied in Ref. [17].
Now, for intra-spin interaction κ = 1 we apply the ansatz functions (13) to (8) and we obtain the interaction energy in the form
E I κ = 1 = g N 2 4 π cos 4 Ω 2 t a ˜ 1 2 ( t ) | a ˜ 1 ( t ) | 4 g N 2 4 π sin 4 Ω 2 t a ˜ 2 2 ( t ) | a ˜ 2 ( t ) | 4 .
On the right hand side of the Equation (17), the first and second terms characterize the interaction energies of the intra-spin components and the Rabi field. They are introduced to demonstrate oscillations of these energies between the spin components. As follows from Equation (14), in Equation (17) factors a ˜ i 2 ( t ) / | a ˜ i ( t ) | 4 characterize the inverse square of the widths of spin components. Then the condensate collapses at least in one spin component, if a ˜ i 2 ( t ) / | a ˜ i ( t ) | 4 0 and E I κ = 1 , respectively. On the other hand, it is difficult to collapse both spin components due to the cosine and sine functions in Equation (17), since due to the fourth power, these functions can remain close to zero for a significantly long time compared to the collapse time.
Since the spin rotation frequency is Ω , the characteristic spin flip time is T s f = π / Ω + 2 π k / Ω , k = 0 , 1 , 2 , . At the time π / 2 Ω , the occupations of the spin components become close, then the interaction energy (17) has decreased by a factor of 2 . As a result, T c = π / 2 Ω could, in general, be a probable moment of collapse of the double spin components under the condition of sufficiently large self-interaction, g N > 4 π . On the other hand, at t = π / 2 Ω it follows from Equation (16) that the Rabi field-related energy can be either positive or negative with corresponding kinetic and self-interaction energies being dependent on the spin state. Therefore, simultaneous collapse in both spin components is difficult to achieve. Eventually, π / 2 Ω is a border time of collapse due to the transfer of density between two spin components. As a result, the condition for the collapse of the first spin component before the first spin flip is T s f > 2 T c ( k = 0 ), i.e., Ω < π g N 2 π / ( 2 2 a 0 2 ) . This means that for a weak Rabi field, the first spin component collapses, and for a sufficiently large Ω , one of the spin components can collapse, as depends on the Rabi field.
The Figure 1 presents the stability and collapse of diagram of the two spin components in the ( g N , Ω ) -plane. In the figure, the panels (a) and (b) are plotted for two different initial widths presented in the caption, and the initial wave function is described by (10). The results show that when the Rabi coupling is switched on, the critical value of g N is shifted up from 2 π to 2 π × 1.33 . The transitions between different types of collapse show a complex domain structure. For g N < 4 π , the Rabi frequency can provide a stable BEC at Ω π g N 2 π / ( 2 2 a 0 2 ) . Also, comparing panels (a) and (b) of the Figure 1, it can be seen that decreasing the initial width of the condensate increases the Rabi collapse field of each spin component, requiring a stronger Rabi field.
The Figure 2 demonstrate dynamics of the condensate width and spin evolution, defined by Equations (3) and (4). The panels (a), (b), (c) a related to the values ( g N , Ω ) given by points (1,2,3) in Figure 1, correspondingly. The panel (a1) presents collapse of the first spin component because for given values of ( g N , Ω ) = ( 1.32 λ cr , 0.1 λ cr ) the condition Ω < π g N 2 π / ( 2 2 a 0 2 ) is satisfied. Nonetheless, by comparing a 1 ( t ) and a ( t ) one can see that this weak Rabi frequency changes the collapse time with respect to Ω = 0 . Also, for the same interaction and Rabi frequency, from the Figure 2(a2) it follows that the total spin is reoriented in another direction rather that follows for a rotation expected due to the Rabi field, because σ x 0 . Consequently, the length of the spin vector P ( t ) < 1 near the T c .
The panel (b1) of Figure 2 shows that a relatively strong Rabi frequency at ( g N , Ω ) = ( 1.32 λ cr , 0.42 λ cr ) can stabilize the BEC and, correspondingly, the spin (b2) is rotating perfectly in the spin subspace ( σ y , σ z ) , with the conserved length P ( t ) 1 . The panel (c1) of Figure 2 demonstrates collapse of the second spin component for the interaction g N > 1.33 λ cr (point (3), ( g N , Ω ) = ( 1.57 λ cr , 0.42 λ cr ) ) . One can see from the panel (c1) that at time t < T s f / 2 the widths a ( t ) and a 1 ( t ) are the same because the interaction is relatively strong therefore first spin component is squeezed. From panel (c2) it follows that at time T s f < T c the spin is flipped to spin-down and then the condensate collapses. Here also close to time T c the spin length is P ( t ) < 1 . Comparison of a ( t ) and collapsing a i ( t ) in panels (a1) and (c1) implies that switching the collapse between spin components increases the collapse time T c compared the realization without the spin-related effects. It may be noticed that in the panels (a2), (b2) and (c2) always σ x < 0 that is reorientation of the initial spin state. The sharp jump of the widths a i ( t ) in panels (b1) and (c1) corresponds to spin flip time T s f .
The Figure 3 presents density cross-section profiles of the condensate for the same values of the Rabi frequency and interaction parameters from panels (a) and (b) of the Figure 2 for collapse and stable, BEC correspondingly. The panels (a1), (a2) and (c1), (c2) of Figure 3 demonstrate the density evolution of the spin components. The panels (b1) and (d1) demonstrate spin rotation of σ z and panels (b2) and (d2) show spin reorientation in the σ x direction. From panels (a1) and (b1) it follows that the condensate collapses in the first spin component at a time T c < T s f because of a weak the Rabi frequency. It follows from panel (a2) that the density of the second spin component is almost zero since spin cannot rotate but it is reoriented because σ x 0 . The panels (c1) and (c2) clearly show the evolution of the density transition between spin-up and spin-down states and panel (d1) shows pure spin rotation due to a large Rabi frequency. It follows from the panels (b2) and (d2) that the initial spin is reoriented and σ x < 0 that can compensate the interaction energy for collapse or stable BEC. All presented density plots are in the same scales shown on the right side with arbitrary units.
As a result, we numerically demonstrated that in the case of intra-spin interaction, the Rabi frequency oscillations lead to different time dependence of the density profiles in the spin up to spin down components. Density redistribution between the components decrease self-interaction in each component. Thus, the Rabi rotation at frequency Ω can produce BEC stable against the collapse or collapse of one of the components. The intra-spin interaction can lead to the reorientation of the spin from that expected for the Rabi rotation by producing σ x < 0 state. In addition, the initial extension of the condensate strongly modifies the effect of the Rabi field on the BEC collapse.

3.2. Cross-Spin Interaction

In this section we discuss the cross-spin interaction with κ = 1 in Equation (1), that can provide collapse of the condensate by applying the Rabi field. Since we take the initial spin (10) along z-axis ( ψ 2 ( x , t = 0 ) = 0 ), then it follows from Equation (1) that to switch on the cross-spin interaction the Rabi frequency should be nonzero.
With the cross-spin coupling the self-interaction energy (8) is given by
E I κ = 1 = g N 2 4 π sin 2 [ Ω t ] a ˜ 12 1 ( t ) ,
where a ˜ 12 ( t ) is
a ˜ 12 ( t ) = | a ˜ 1 ( t ) | 2 a ˜ 1 ( t ) + | a ˜ 2 ( t ) | 2 a ˜ 2 ( t ) .
On the other hand, at the moment of time T k = π k / Ω the energy (18) also becomes zero due to sin 2 [ Ω t ] function. Thus, a sufficient self-interaction occurs only if the BEC spin | σ z ( t ) | < 1 . Eventually, the condensate can collapse only at a ˜ 12 ( t ) 0 , and it follows from Equation (19) that the collapse can be ensured only in two spin components simultaneously.
The Figure 4 demonstrates the stability and collapse diagram of two spin components in the ( g N , Ω ) -plane. It follows from the figure that collapse can occur as Ω exceeding a certain critical value of 0.4 π and sufficiently large g N > 6 π .
The Figure 5 demonstrates the collapse and stable dynamics of the BEC for given values of ( g N , Ω ) from fixed points (1) and (2) of Figure 4, correspondingly. The panels (a1) and (b1) correspond to the time evolution of the width, and the panels (a2) and (b2) correspond to the total spin rotation dynamics. The plot (a1) demonstrates the widths of both spin components tend to zero, as consistent with the interaction energy (18). Near the collapse time T c , the widths show some increase since, as it follows from plot (a2), at this time the spin is flipped to spin-down and σ x 0 leading to a decrease in the self-interaction. Thus, this behavior of the widths of the spin components can be explained as follows. Since E R σ x , it follows from the plot (a2) that the attraction interaction energy (18) is compensated by the Rabi magnetic field energy at σ x > 0 and a sharp change in this energy at σ x 0 leads to expanding of the condensate. Finally, due to a relatively fast spin rotation, one obtains σ x > 0 and consequently, the condensate collapses. Also, due to this interplay of the energies in the panel (a1) one can see the collapse time with cross-spin interaction is considerably larger than a ( t ) that corresponds to the spinless condensate. The panel (b1) of Figure 5 demonstrates the widths of spin components with their interrelated behavior for stable condensate with the values of ( g N , Ω ) from the fixed point (2) of Figure 4. The sharp changes in the widths correspond to transitions between the spin up and down states with the Rabi frequency. It follows from panel (b2) that the spin is rotating in the ( σ y , σ z ) plane with the length P ( t ) 1 . In this case, due to fast oscillations of σ x ( t ) , Rabi magnetic field cannot inhibit collapse, but it still keeps the condensate from expanding for a relatively long time. The panels (c1) and (c2) demonstrate cross-section density plot of spin components and panels (d1) and (d2) show corresponding spin rotation dynamics in ( x , t ) plane. In fact, the panels (c) and (d) present the density plot simulation of the panels (a1) and (a2), correspondingly. Here the condensate is squeezed and spin dynamics are clearly demonstrated.
To summarize, in this section we presented numerical evidence that the condensate with cross-spin interaction and initially non-coupled spin components can collapse in the presence of Rabi magnetic field. For the self-interaction exceeding some critical value, the collapse of the condensate can be produced by spin reorientation from the initial state even for a relatively small Ω . This implies that the self-interaction energy could be compensated by interaction with the Rabi field when the spin state with σ x ( t > 0 ) > 0 is formed.
Moreover, Figure 6 demonstrates for intra- and cross-spin interaction that the initially Gaussian wave functions become non-Gaussian with time. The plots of Figure 6 show the evolution of the expectation values of r 4 i and r 2 i 2 defined as:
r d i = 2 π N i 0 r d | ψ i ( r , t ) | 2 r d r ,
where ψ i ( r , t ) is numerical solution of Equation (1) and N i is defined by Equation (2). For two-dimensional Gaussians one has r 4 i 2 r 2 i 2 = 0 , and the nonzero values demonstrate that the BEC shape strongly deviates from the Gaussian profile.

4. Conclusions

In this paper, we investigated the collapse of quasi-two-dimensional BEC of atoms with pseudospin-1/2 and in a synthetic Rabi magnetic field. The collapse of condensate is explored assuming intra- or cross-spin attractive interactions. With the intra-spin interaction, the Rabi magnetic field affects the collapse of the condensate by modifying number of atoms in each spin component. In this case the spin rotation leads to oscillations in the self-interaction energy with the double Rabi frequency that can stabilize the condensate. With the cross-spin interaction, the condensate collapses by rotation of the spin, which enhances the self-interaction.
We obtained that for the intra- and cross-spin coupling both, the initially Gaussian wave function becomes non-Gaussian at time t > 0 . As a result, we determine collapse and stability diagram in the ( g N , Ω ) plane for both interactions. The intra-spin-coupling modifies the critical number of atoms causing the collapse while the collapse is observed only in a single pseudospin component. It is demonstrated that for cross-spin-coupling only double spin-components collapse can occur.
To connect the results with the realization of the BEC experiments, for s-wave scattering length | a s | 100 a B 5 × 10 3 μ m ( a B being the Bohr radius) and the condensate extension along the z axis a z 1 μ m , the dimensionless interaction constant g = 4 π | a s | / a z 0.05 . The velocity of the collapse is v c g N / M a ( 0 ) . For the BEC of 87 Rb atoms at a ( 0 ) 10   μ m and N 10 4 this estimate yields v c 0.1 cm/s and the corresponding time scale T c = a ( 0 ) / v c 0.01 s. Therefore, the Rabi frequency providing the cross-over condition Ω T c 1 is of the order of Ω 100 Hz. The conditions of weak and strong Rabi frequency against the BEC collapse, Ω T c 1 and Ω T c 1 , correspondingly, can be achieved by modifying initial a ( 0 ) , the number of particles N , and Rabi frequency Ω in the available experimental range [30].

Author Contributions

Conceptualization, S.N.M. and B.J.A.; methodology, B.J.A. and S.N.M.; software, S.N.M.; validation, B.J.A.; formal analysis, B.J.A. and S.N.M.; investigation, S.N.M. and B.J.A.; resources, S.N.M.; writing—original draft preparation, S.N.M.; writing—review and editing, B.J.A.; visualization, S.N.M.; supervision, S.N.M. and B.J.A.; project administration, B.J.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No data is available.

Acknowledgments

This research is supported by Grant F-FA-2021-432 and Sh.M. partly is supported by the Fund for Financing Science and Innovation Support of Ministry for Innovative Development of Uzbekistan.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BECBose-Einstein condensate
GPEGross-Pitaevskii equation

References

  1. Rajaraman, R. Solitons and Instantons: An Introduction to Solitons and Instantons in Quantum Field Theory; North-Holland Personal Library, North-Holland Publishing Company: Amsterdam, The Netherlands, 1982. [Google Scholar]
  2. Dalfovo, F.; Giorgini, S.; Pitaevskii, L.P.; Stringari, S. Theory of Bose-Einstein condensation in trapped gases. Rev. Mod. Phys. 1999, 71, 463. [Google Scholar] [CrossRef] [Green Version]
  3. Sikivie, P.; Yang, Q. Bose-Einstein Condensation of Dark Matter Axions. Phys. Rev. Lett. 2009, 103, 111301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Chin, C.; Grimm, R.; Julienne, P.; Tiesinga, E. Feshbach resonances in ultracold gases. Rev. Mod. Phys. 2010, 82, 1225. [Google Scholar] [CrossRef]
  5. Giovanazzi, S.; Görlitz, A.; Pfau, T. Tuning the Dipolar Interaction in Quantum Gases. Phys. Rev. Lett. 2002, 89, 130401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Kumar, R.K.; Tomio, L.; Gammal, A. Spatial separation of rotating binary Bose-Einstein condensates by tuning the dipolar interactions. Phys. Rev. A 2019, 99, 043606. [Google Scholar] [CrossRef] [Green Version]
  7. Kivshar, Y.S.; Malomed, B.A. Dynamics of solitons in nearly integrable systems. Rev. Mod. Phys. 1989, 61, 763. [Google Scholar] [CrossRef]
  8. Ferrier-Barbut, I.; Kadau, H.; Schmitt, M.; Wenzel, M.; Pfau, T. Observation of Quantum Droplets in a Strongly Dipolar Bose Gas. Phys. Rev. Lett. 2016, 116, 215301. [Google Scholar] [CrossRef]
  9. Semeghini, G.; Ferioli, G.; Masi, L.; Mazzinghi, C.; Wolswijk, L.; Minardi, F.; Modugno, M.; Modugno, G.; Inguscio, M.; Fattori, M. Self-Bound Quantum Droplets of Atomic Mixtures in Free Space. Phys. Rev. Lett. 2018, 120, 235301. [Google Scholar] [CrossRef] [Green Version]
  10. Bergé, L. Wave collapse in physics: Principles and applications to light and plasma waves. Phys. Rep. 1998, 303, 259. [Google Scholar] [CrossRef]
  11. Greiner, M.; Mandel, O.; Hänsch, T.W.; Bloch, I. Collapse and revival of the matter wave field of a Bose–Einstein condensate. Nature 2002, 419, 51. [Google Scholar] [CrossRef]
  12. Donley, E.A.; Claussen, N.R.; Cornish, S.L.; Roberts, J.L.; Cornell, E.A.; Wieman, C.E. Dynamics of collapsing and exploding Bose–Einstein condensates. Nature 2001, 412, 295. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Sackett, C.A.; Stoof, H.T.C.; Hulet, R.G. Growth and Collapse of a Bose-Einstein Condensate with Attractive Interactions. Phys. Rev. Lett. 1998, 80, 2031. [Google Scholar] [CrossRef] [Green Version]
  14. Cornish, S.L.; Claussen, N.R.; Roberts, J.L.; Cornell, E.A.; Wieman, C.E. Stable 85Rb Bose-Einstein Condensates with Widely Tunable Interactions. Phys. Rev. Lett. 2000, 85, 1795. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Mardonov, S.; Modugno, M.; Sherman, E.Y.; Malomed, B.A. Rabi-coupling-driven motion of a soliton in a Bose-Einstein condensate. Phys. Rev. A 2019, 99, 013611. [Google Scholar] [CrossRef] [Green Version]
  16. Abdullaev, F.K.; Salerno, M. Gap-Townes solitons and localized excitations in low-dimensional Bose-Einstein condensates in optical lattices. Phys. Rev. A 2005, 72, 033617. [Google Scholar] [CrossRef] [Green Version]
  17. Abdullaev, F.K.; Caputo, J.G.; Kraenkel, R.A.; Malomed, B.A. Controlling collapse in Bose-Einstein condensates by temporal modulation of the scattering length. Phys. Rev. A 2003, 67, 013605. [Google Scholar] [CrossRef] [Green Version]
  18. Kagan, Y.; Surkov, E.L.; Shlyapnikov, G.V. Evolution of a Bose-condensed gas under variations of the confining potential. Phys. Rev. A 1996, 54, R1753. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Petrov, D.S. Quantum Mechanical Stabilization of a Collapsing Bose-Bose Mixture. Phys. Rev. Lett. 2015, 115, 155302. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Yu, Z.F.; Zhang, A.X.; Tang, R.A.; Xu, H.P.; Gao, J.M.; Xue, J.K. Spin-orbit-coupling stabilization of a collapsing binary Bose-Einstein condensate. Phys. Rev. A 2017, 95, 033607. [Google Scholar] [CrossRef]
  21. Wang, L.X.; Dai, C.Q.; Wen, L.; Liu, T.; Jiang, H.F.; Saito, H.; Zhang, S.G.; Zhang, X.F. Dynamics of vortices followed by the collapse of ring dark solitons in a two-component Bose-Einstein condensate. Phys. Rev. A 2018, 97, 063607. [Google Scholar] [CrossRef]
  22. Kartashov, Y.V.; Sherman, E.Y.; Malomed, B.A.; Konotop, V.V. Stable two-dimensional soliton complexes in Bose–Einstein condensates with helicoidal spin–orbit coupling. New J. Phys. 2020, 22, 103014. [Google Scholar] [CrossRef]
  23. Zeng, L.; Zeng, J. Preventing critical collapse of higher-order solitons by tailoring unconventional optical diffraction and nonlinearities. Commun. Phys. 2020, 3, 26. [Google Scholar] [CrossRef]
  24. Sakaguchi, H.; Sherman, E.Y.; Malomed, B.A. Vortex solitons in two-dimensional spin-orbit coupled Bose-Einstein condensates: Effects of the Rashba-Dresselhaus coupling and Zeeman splitting. Phys. Rev. E 2016, 94, 032202. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Kartashov, Y.V.; Konotop, V.V.; Modugno, M.; Sherman, E.Y. Solitons in Inhomogeneous Gauge Potentials: Integrable and Nonintegrable Dynamics. Phys. Rev. Lett. 2019, 122, 064101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Chen, Z.; Li, Y.; Malomed, B.A.; Salasnich, L. Spontaneous symmetry breaking of fundamental states, vortices, and dipoles in two- and one-dimensional linearly coupled traps with cubic self-attraction. Phys. Rev. A 2017, 96, 033621. [Google Scholar] [CrossRef] [Green Version]
  27. Mardonov, S.; Modugno, M.; Sherman, E.Y. Dynamics of a macroscopic spin qubit in spin–orbit coupled Bose–Einstein condensates. J. Phys. B At. Mol. Opt. Phys. 2015, 48, 115302. [Google Scholar] [CrossRef]
  28. Mardonov, S.N. Collapse of quasi-two-dimensional symmetrized Dresselhaus spin-orbit coupled Bose–Einstein condensate. Arab. J. Math. 2021, 11. [Google Scholar] [CrossRef]
  29. Mardonov, S.; Sherman, E.Y.; Muga, J.G.; Wang, H.W.; Ban, Y.; Chen, X. Collapse of spin-orbit-coupled Bose-Einstein condensates. Phys. Rev. A 2015, 91, 043604. [Google Scholar] [CrossRef] [Green Version]
  30. Kartashov, Y.V.; Astrakharchik, G.E.; Malomed, B.A.; Torner, L. Frontiers in multidimensional self-trapping of nonlinear fields and matter. Nat. Rev. Phys. 2019, 1, 185. [Google Scholar] [CrossRef]
Figure 1. The diagram of stable and collapsing condensate is plotted in  ( g N , Ω ) plane for intra-spin interaction κ = 1 given in Equation (1). The left panel (a) and the right panel (b) are plotted for initial width of condensate a 0 = 1 and a 0 = 1 / 2 , correspondingly. Below in the next figures we will present characteristic dynamics of the condensate for the values of ( g N , Ω ) from the fixed points (1,2,3) in the panel (a).
Figure 1. The diagram of stable and collapsing condensate is plotted in  ( g N , Ω ) plane for intra-spin interaction κ = 1 given in Equation (1). The left panel (a) and the right panel (b) are plotted for initial width of condensate a 0 = 1 and a 0 = 1 / 2 , correspondingly. Below in the next figures we will present characteristic dynamics of the condensate for the values of ( g N , Ω ) from the fixed points (1,2,3) in the panel (a).
Particles 05 00012 g001
Figure 2. The plot of spin-components width (a1c1) and spin rotation dynamics (a2c2) versus time for intra-spin interaction κ = 1 in Equation (1). The panels (ac) correspond to the values of ( g N , Ω ) given as ( 1.32 λ cr , 0.1 λ cr ) , ( 1.32 λ cr , 0.42 λ cr ) , ( 1.57 λ cr , 0.42 λ cr ) from fixed points (1,2,3) in the Figure 1, correspondingly. In the panels (a1,c1) the width a ( t ) are plotted for Ω = 0 to compare the collapse times with non-spin effect.
Figure 2. The plot of spin-components width (a1c1) and spin rotation dynamics (a2c2) versus time for intra-spin interaction κ = 1 in Equation (1). The panels (ac) correspond to the values of ( g N , Ω ) given as ( 1.32 λ cr , 0.1 λ cr ) , ( 1.32 λ cr , 0.42 λ cr ) , ( 1.57 λ cr , 0.42 λ cr ) from fixed points (1,2,3) in the Figure 1, correspondingly. In the panels (a1,c1) the width a ( t ) are plotted for Ω = 0 to compare the collapse times with non-spin effect.
Particles 05 00012 g002
Figure 3. The density plot of the condensate parameters taken from panels (a,b) of Figure 2 is presented for the corresponding values of ( g N , Ω ) . The panels (a1,c1) are demonstrating | ψ 1 | 2 and (a2,c2) are demonstrating | ψ 2 | 2 . The panels (b1,d1) are demonstrating σ z and (b2,d2) are demonstrating σ x . The plots are presented in the plane ( x , t ) for intra-spin interaction κ = 1 in Equation (1). The values of ( g N , Ω ) are given from the points (1, 2) of the Figure 1, correspondingly, for collapsing and stable BEC. The density scales on the right side are the same for all panels with arbitrary units. Here the panels (a,b) are presenting collapse of the first spin component with corresponding spin dynamics and the panels (c,d) are presenting stable condensate with total density transition between spin components and corresponding spin rotation.
Figure 3. The density plot of the condensate parameters taken from panels (a,b) of Figure 2 is presented for the corresponding values of ( g N , Ω ) . The panels (a1,c1) are demonstrating | ψ 1 | 2 and (a2,c2) are demonstrating | ψ 2 | 2 . The panels (b1,d1) are demonstrating σ z and (b2,d2) are demonstrating σ x . The plots are presented in the plane ( x , t ) for intra-spin interaction κ = 1 in Equation (1). The values of ( g N , Ω ) are given from the points (1, 2) of the Figure 1, correspondingly, for collapsing and stable BEC. The density scales on the right side are the same for all panels with arbitrary units. Here the panels (a,b) are presenting collapse of the first spin component with corresponding spin dynamics and the panels (c,d) are presenting stable condensate with total density transition between spin components and corresponding spin rotation.
Particles 05 00012 g003
Figure 4. The diagram of the stable and collapsing BEC is presented in ( g N , Ω ) -plane for cross-spin interaction with κ = 1 in Equation (1). The plot corresponds to initial width a 0 = 1 . For the values of ( g N , Ω ) from fixed points (1,2) characteristic condensate dynamics will be presented below in the next figures.
Figure 4. The diagram of the stable and collapsing BEC is presented in ( g N , Ω ) -plane for cross-spin interaction with κ = 1 in Equation (1). The plot corresponds to initial width a 0 = 1 . For the values of ( g N , Ω ) from fixed points (1,2) characteristic condensate dynamics will be presented below in the next figures.
Particles 05 00012 g004
Figure 5. The panels (a1,b1) present the widths of the spin components defined by Equation (3) and the panels (a2,b2) present the spin rotation dynamics defined by Equation (4) versus time. The panels (a,b) are plotted for the values of ( g N , Ω ) given as ( 3.3 λ cr , 0.3 λ cr ) , ( 3.3 λ cr , 0.5 λ cr ) from the fixed points (1, 2) of Figure 4 for collapse and stable BEC, correspondingly. In the panel (a1) the width a ( t ) is plotted for Ω = 0 to compare the collapse time with that in the spinless condensate. The panels (c1,c2) are the cross section density plots of the collapse dynamics of | ψ 1 | 2 and | ψ 2 | 2 , and the panels (d1,d2) present the cross section density plot of spin rotation of σ z and σ x in the plane ( x , t ) , with correspondingly values of ( g N , Ω ) from the panel (a), correspondingly.
Figure 5. The panels (a1,b1) present the widths of the spin components defined by Equation (3) and the panels (a2,b2) present the spin rotation dynamics defined by Equation (4) versus time. The panels (a,b) are plotted for the values of ( g N , Ω ) given as ( 3.3 λ cr , 0.3 λ cr ) , ( 3.3 λ cr , 0.5 λ cr ) from the fixed points (1, 2) of Figure 4 for collapse and stable BEC, correspondingly. In the panel (a1) the width a ( t ) is plotted for Ω = 0 to compare the collapse time with that in the spinless condensate. The panels (c1,c2) are the cross section density plots of the collapse dynamics of | ψ 1 | 2 and | ψ 2 | 2 , and the panels (d1,d2) present the cross section density plot of spin rotation of σ z and σ x in the plane ( x , t ) , with correspondingly values of ( g N , Ω ) from the panel (a), correspondingly.
Particles 05 00012 g005
Figure 6. The plot of non-Gaussian behavior of the spin components versus time. Blue lines correspond to the spin up and red dashed lines corresponds to the spin down. The panels (a,b) correspond to intra-spin interaction with the value of ( g N , Ω ) from fixed points (1) and (2) of Figure 1, correspondingly. The panel (c) corresponds to cross-spin interaction with the value of ( g N , Ω ) from fixed points (1) of the Figure 4.
Figure 6. The plot of non-Gaussian behavior of the spin components versus time. Blue lines correspond to the spin up and red dashed lines corresponds to the spin down. The panels (a,b) correspond to intra-spin interaction with the value of ( g N , Ω ) from fixed points (1) and (2) of Figure 1, correspondingly. The panel (c) corresponds to cross-spin interaction with the value of ( g N , Ω ) from fixed points (1) of the Figure 4.
Particles 05 00012 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mardonov, S.N.; Ahmedov, B.J. Rabi Frequency Management of Collapsing Quasi-Two-Dimensional Bose-Einstein Condensates with Pseudospin-1/2. Particles 2022, 5, 135-145. https://0-doi-org.brum.beds.ac.uk/10.3390/particles5020012

AMA Style

Mardonov SN, Ahmedov BJ. Rabi Frequency Management of Collapsing Quasi-Two-Dimensional Bose-Einstein Condensates with Pseudospin-1/2. Particles. 2022; 5(2):135-145. https://0-doi-org.brum.beds.ac.uk/10.3390/particles5020012

Chicago/Turabian Style

Mardonov, Shukhrat N., and Bobomurat J. Ahmedov. 2022. "Rabi Frequency Management of Collapsing Quasi-Two-Dimensional Bose-Einstein Condensates with Pseudospin-1/2" Particles 5, no. 2: 135-145. https://0-doi-org.brum.beds.ac.uk/10.3390/particles5020012

Article Metrics

Back to TopTop