# Non-autonomous bright-dark solitons and Rabi oscillations in multi-component Bose-Einstein condensates

## Abstract

We study the dynamics of non-autonomous bright-dark matter-wave solitons in two- and three- component Bose-Einstein condensates. Our setting includes a time-dependent parabolic potential and scattering length, as well as Rabi coupling of the separate hyperfine states. By means of a similarity transformation, we transform the non-autonomous coupled Gross-Pitaevskii equations into the completely integrable Manakov model with defocusing nonlinearity, and construct the explicit form of the non-autonomous soliton solutions. The propagation characteristics for the one-soliton state, and collision scenarios for multiple soliton states are discussed in detail for two types of time-dependent nonlinearities: a kink-like one and a periodically modulated one, with appropriate time-dependence of the trapping potential. We find that in the two-component condensates the nature of soliton propagation is determined predominantly by the nature of the nonlinearity, as well as the temporal modulation of the harmonic potential; switching in this setting is essentially due to Rabi coupling. We also perform direct numerical simulation of the non-autonomous two-component coupled Gross-Pitaevskii equations to corroborate our analytical predictions. More interestingly, in the case of the three-component condensates, we find that the solitons can lead to collision-induced energy switching (energy-sharing collision), that can be profitably used to control Rabi switching or vice-versa. An interesting possibility of reversal of the nature of the constituent soliton, i.e., bright (dark) into dark (bright) due to Rabi coupling is demonstrated in the three-component setting.

###### pacs:

03.75.Mn, 03.75.Lm, 02.30.Ikkanna_phy@bhc.edu.in (corresponding author), dfrantz@phys.uoa.gr

Journal Reference: J. Phys. A: Math. Theor. 46 (2013) 475201

## 1 Introduction

Studies on atomic Bose-Einstein condensates (BECs) have received considerable attention in recent years [1]. In this context, it has been shown that nonlinear effects in matter waves can emerge, with this particular direction having attracted much interest both in theory and in experiments [2, 3]. From a theoretical viewpoint, this interest arises from the fact that many of such nonlinear effects can be understood in the framework of lowest-order mean-field theory, namely by means of the so-called Gross-Pitaevskii equation (GPE), which is nothing but the ubiquitous nonlinear Schrödinger (NLS) equation with external potential [1, 2, 3]. The GPE is a nonlinear evolution equation (with the nonlinearity originating from the interatomic interactions) and, as such, it permits the study of a variety of interesting purely nonlinear phenomena. The latter, have primarily been studied by treating the condensate as a purely nonlinear coherent matter-wave, i.e., from the viewpoint of the dynamics of nonlinear waves. Relevant studies have already been summarized in various books [2] and reviews; see, e.g., [4] for bright solitons, [5] for dark solitons, [6] for vortices, [7] for dynamical instabilities in BECs, and so on.

The formation and dynamical properties of solitons in BECs are determined by the nature of their two-body atomic interactions, i.e., the sign of the -wave scattering length which may be positive (negative) for repulsive (attractive) inter-atomic interactions. For instance, atomic dark (bright) solitons are formed in condensates with repulsive (attractive) interactions [2, 3, 4, 5] (note that bright gap solitons are also possible in repulsive BECs, but those are critically formed due to the presence of a periodic—optical lattice—potential [8]). On the other hand, it should be stressed that the above conditions can dramatically change in the case of multi-component BECs, composed by two or more different hyperfine states of the same atom species (e.g., Rb) [9, 10]: in such a case, e.g., binary condensates with repulsive inter- and intra-species interactions can support mixed bright-dark (BD) solitons. This type of vector soliton consists of a dark soliton in one component coupled to a bright soliton in the second component. These solitons are usually referred to as “symbiotic” solitons, because the bright component cannot be supported in a stand-alone fashion (it is only supported as such in attractive BECs [4]; see also [5]), and is only sustained because of the presence of its dark counterpart, which acts as an effective external trapping potential. Such atomic BD solitons were predicted in theory [11] and observed in experiments from different experimental groups [12, 13, 14]. Notice that solitons in three-component condensates (e.g., spinor BECs) have been studied too [15, 16, 17] and mixed, BD solitons were predicted to occur as well [18].

On the other hand, coming back to the role of the sign and magnitude of the scattering length, it is important to note that they can be adjusted over a relatively large range by employing external magnetic, electric or optical fields near Feshbach resonances [19]. This possibility has given rise to many theoretical and experimental studies, with a prominent example being the formation of bright matter-wave solitons and soliton trains in attractive condensates [20], by switching the interatomic interactions from repulsive to attractive. Many theoretical works studied the BEC dynamics under temporal modulation of the nonlinearity. In particular, the application of such a “Feshbach resonance management” (FRM) technique [21] can be used to stabilize attractive higher-dimensional BEC against collapse [22], to create robust matter-wave breathers in the effectively one-dimensional (1D) condensate [21, 23], or compress bright solitons in the presence of expulsive potentials [24]. Mathematically speaking, if the scattering length is time dependent then so does the nonlinearity coefficient in the GPE, and the latter becomes non-autonomous. In the pioneering works [25], it has been shown that the single-component non-autonomous GPE can be transformed into the standard integrable NLS equation. The procedure developed in Ref. [25] can be generalized to two-component condensates [26] and the corresponding non-autonomous matter-wave solitons were studied extensively with aid of the explicit soliton solutions reported in Refs. [27, 28, 29].

This work deals with the dynamics of multi-component BECs, particularly two- and three-component BECs, with temporally modulated scattering length and trapping potential, and in the presence of Rabi coupling; the latter is accounted for by a linear coupling between separate wave functions induced by a radio frequency [9] (see also Refs. [30, 31, 32] for a theoretical analysis and applications). Recently, coherent many body Rabi oscillations have been observed experimentally and theoretically discussed in Ref. [33]. Using two successive transformations, with one constraint equation being in the form of a Ricatti equation, we reduce the original multicomponent non-autonomous GPEs to the integrable defocusing vector NLS equation, known as the Manakov model [34]. These transformations can be viewed as a rotation followed by a similarity transformation. This way, we find non-autonomous mixed (BD) solitons and analyze their dynamics and interactions in detail, in both two- and three-component settings. We find that in the two-component system, soliton propagation is determined predominantly by the nature of the nonlinearity, as well as the temporal modulation of the harmonic potential; switching in this setting is essentially due to Rabi coupling. Our analytical predictions are corroborated with results from direct simulations; very good agreement between the two is found. Then by going one step further in the case of the three-component system, we find that the solitons can lead to collision-induced energy switching (energy-sharing collision) that can be profitably used to control Rabi switching or vice-versa. An interesting possibility of reversal of the nature of the constituent soliton, i.e., bright (dark) into dark (bright) induced by the Rabi coupling is demonstrated in the three-component setting.

The remaining part of this paper is arranged as follows. In Sec. II, we introduce the generalized non-autonomous coupled GPEs, for two-component condensates, with time-varying interaction strength and Rabi coupling. Also, we construct the one- and two-soliton solutions to the 1D coupled non-autonomous GPEs and analyze their propagation characteristics. In Sec. III, we consider the three-component non-autonomous GPEs with Rabi coupling, and explore the interesting energy exchange collisional features of non-autonomous matter-wave solitons. Finally, in Sec. IV, we present our conclusions.

## 2 Non-autonomous BD solitons in two-component BECs and Rabi oscillations

In the mean field approximation, the dynamics of two-component condensates can be well described by a pair of two coupled three-dimensional (3D) GPEs (see, e.g., Ref. [35]). In the physically relevant case of highly anisotropic (cigar-shaped) traps, these GPEs can be reduced to an effectively 1D system (see, e.g., Ref. [3]), which is of the following dimensionless form:

(1) |

Here, is the macroscopic wave function of the component, time and spatial coordinate are respectively measured in units of the inverse transversal frequency and the transverse harmonic oscillator length ( denotes atomic mass). In Eq. (1) the nonlinearity coefficients are proportional to the -wave scattering lengths , namely (where is the Bohr radius), and are assumed to be subject to the FRM technique (see below); furthermore are the external potentials [for () the potential is expulsive (confining)], and the linear cross coupling parameters represent the Rabi coupling: in fact these coefficients are proportional to the Rabi frequency, and are responsible for the transfer of atoms of component to component and vice-versa. In our case, considering two different hyperfine states of the same atom, the external time-varying potentials are equal, i.e., . Here, the strength of the parabolic trap is given by , where is the temporally modulated axial trap frequency. Finally, the macroscopic wave functions are normalized such that (), where is the number of atoms in the component.

### 2.1 Non-autonomous BD one-soliton solution

Based on experimental results pertaining to two-component Rb BECs (see, e.g., [13]), we can assume that the scattering length ratios are almost equal to one, and also they can be tuned through Feshbach resonance (as it was demonstrated in the experimental works [36, 37]). Hence, we consider the case of equal interaction strengths, i.e., in Eq. (1), and choose the linear coupling coefficients as , with .

In such a case, as shown in Appendix A, we can use two successive transformations to map the non-autonomous system (1) into an integrable autonomous nonlinear system. This is possible mainly due to the temporal dependence of the nonlinearity and external harmonic potential; in the absence of such dependence these transformations are not possible, in general. Notice, however, that in the absence of potential (homogeneous system), one can transform Eq. (1) into an integrable coupled nonlinear Schrödinger (CNLS) system even in the absence of inhomogeneities. The above mentioned transformations are: (i) a unitary transformation (see, e.g., Refs. [30, 32]) to reduce Eq. (1) to a non-autonomous system of equations without the Rabi coupling term; (ii) a similarity transformation, which reduces the aforementioned system to the following integrable defocusing CNLS equations [38, 39]:

(2) |

Equation (2) is the completely integrable Manakov system and admits -soliton solutions of dark-dark and mixed (BD) types. Here, we focus only on mixed type soliton solutions due to their recent experimental observations and for their special dynamical features [12, 13].

The mixed one-soliton solution of Eq. (2) can be obtained by Hirota’s direct method [41]. In standard form, the mixed one-soliton solution of the two-component defocusing CNLS equations with the bright (dark) part appearing in the () component (see Refs. [27, 39]) is given in Appendix B.1. Using the analytical form of this type of soliton, as well as the transformations of Appendix A, we find an exact analytical non-autonomous soliton solution of Eq. (1) in the form:

(3) |

where and are given by:

(4) | |||

(5) |

where , and the other parameters and functions involved in the above equations are defined in Appendices A and B.

The above mixed soliton bears resemblance to the recently observed, so-called “beating” dark-dark soliton [14, 42], where the bright and dark parts coexist in the same component with a constant asymptotic value as ; however, the form of this BD soliton still depends on the Rabi coupling and the time-dependent nonlinearity and parabolic potential. Notice that due to the presence of Rabi coupling, there will be a periodic switching between the two components along with oscillating background. Another important observation for the integrable two-component case is the dependence of the existence of e.g., the bright part of the mixed soliton on its dark counterpart (see Eqs. (6gh)-(6gi) and relevant discussion in Appendix B.1). Thus, it is impossible to make any one of the constituents (bright/dark) of the mixed soliton to vanish completely. In fact, this will result in singularities in the other soliton part. This holds even in the presence of the Rabi coupling for the two-component system.

To elucidate the dynamics of such special type of mixed soliton we consider the following two examples.

Example 1: Kink-like nonlinearity

First we note that a fast time modulated nonlinearity displays interesting dynamics in BECs [43]. A relatively sudden jump in the nonlinearity coefficient can be well represented by a kink-like nonlinearity of the form:

(6a) | |||

with the associated atomic scattering length being ; here, denotes the time scale characterizing the jump, and is an arbitrary constant. | |||

This is shown in Fig. 1(a) and the corresponding external harmonic potential is displayed in Fig. 1(b). The time-dependent nonlinearity (t) admitting the above form can be realized for the following scattering length [25], | |||

(6b) | |||

where is the value of scattering length far from the Feshbach resonance, is the applied time-varying magnetic field, is the resonant value of the magnetic field, and is the resonance width in the presence of the magnetic field. The corresponding strength of the time- dependent magnetic trap is determined by a Ricatti equation [see Eq. (6gv) in Appendix A] as follows: | |||

(6c) |

The propagation of the non-autonomous mixed soliton in the absence and presence of Rabi coupling are depicted in the top panels of Figs. 2 and 3, respectively. We have also solved the system (1) numerically by means of the split-step Fourier method, and the resulting numerical plots for the above choice of the time-dependent nonlinearity and potential strength (6c) are given in the bottom panel.

Comparing the top and bottom panels, one can see that there is a very good agreement between numerical and analytical results. We find that, in the absence of the Rabi coupling, there are no oscillations for the solitons and the background. However, the shape, width and velocity of solitons are modulated. Then, the introduction of Rabi term leads to an exchange of number of atoms between the components, which ultimately results in beating oscillations, as shown in Fig. 3. For small values of , the background just starts to oscillate (see Fig. 3). However, by increasing the parameter , one can observe rapid oscillations of the background also. To illustrate this, the contour plots of for and are shown in Fig. 4. Similar observations can be made for the second component, , as well. These oscillations of the background are due to the exchange of condensates between the soliton and the background.

Another important observation following from Figs. 3 and 4, as well as from expressions (3)-(5) is that the increase of the value of Rabi coupling leads a significant portion of the dark part of mixed soliton to appear in the first component (), along with the bright part accompanied by beating effects. Similar effects take place in component too. Thus, in a given component one can have co-existing oscillating BD soliton.

Enlarging Fig. 3 in the region to , we observe that the non-autonomous mixed soliton is absent in this region. This is due to the form of the nonlinearity , which becomes zero in the range = to and reaches a saturation for large positive values. Thus, even though the Rabi coupling leads to an exchange of atoms among the two components leading to oscillations, the nature of soliton propagation is predominantly determined by the nature of time-dependent nonlinearity.

Example 2: Periodically modulated nonlinearity

Next, we consider another physically interesting example, namely the case of a periodic time-varying nonlinearity of the form:

(6ga) | |||

where is an arbitrary real constant, is the characteristic frequency, and is a real constant parameter. In this case, the form of atomic scattering length is . Figure 5(a) represents the form of such a time-varying nonlinearity. The expression for the pertinent time-dependent external magnetic field required to achieve the above form of the nonlinearity can be determined from the formula (6b). The corresponding strength of the magnetic trap admits the form | |||

(6gb) |

in order to satisfy Eq. (6gv). Note that the potential is sign reversible and can support the same type of soliton for both positive and negative signs; see Fig. 5(b).

It is apparent from the periodic nature of the nonlinearity, and the time modulation of the strength of the external potential, that the condensates also execute periodic oscillations in both the components even in the absence of linear coupling as shown in Fig. 6. The numerical results, obtained by a direct integration of Eq. (1) are shown in the bottom panel of Fig. 6 in the absence of Rabi coupling.

Next, we introduce the Rabi coupling between the two components; then, as expected, there will be an exchange of atoms between the components, as well as between the soliton and the background. This exchange induces oscillations in the density of condensates which, in turn, modulate the oscillations due to the periodic nature of and and result in soliton beating in both components. This becomes clear in Fig. 7, from the plots showing the non-autonomous soliton (3) for the choice (6ga) (upper panel) and the corresponding numerical results (bottom panel).

As in the previous example, here also both parts of the mixed soliton co-exist in a given component due to Rabi switching.
The left panel of Fig. 8 displays the propagation of mixed soliton for a small value of (=). It is observed that the soliton only oscillates periodically, but the background is not oscillating. The right panel of Fig. 8 displays that for larger value of (say ), the background oscillates rapidly and significant switching of dark and bright parts among the components occurs, thereby resulting in the co-existence of both dark and bright parts in the same component. By increasing the values of , we also observe “creeping” soliton propagation with beating effects, which are not presented here. Similar type of creeping soliton appears in the presence of inhomogeneities too.

From the above two examples, we can conclude that in the non-autonomous two-component GPEs. (1) the nature of soliton propagation is determined predominantly by the temporal dependence of and while the switching of condensates is completely dependent on the Rabi coupling.

### 2.2 Two soliton solution of non-autonomous two-component GPEs

We now turn our attention to the case of mixed two-soliton solution of the integrable 2-CNLS system (2). The explicit form of this solution, given in Appendix B.1, contains all the information regarding the collision of two solitons. We will employ this autonomous soliton solution to construct the exact two soliton solution of the corresponding non-autonomous system (1) in detail.

The asymptotic analysis of the two-soliton solution of the two component GPEs system (2) given in Appendix C.1, shows that the two solitons undergo standard elastic collision as =, , where represents the amplitude of soliton in component. Here and in the following, the superscript (subscript) of A (or ) represents the number of soliton (component) while -(+) appearing in the corresponding quantities indicates their form before (after) collision. Thus the amplitude and speed of the bright and dark components are preserved after the collision, except for a phase-shift , where and are defined in Appendix . This elastic collision of solitons is shown in Fig. 9.

The exact dynamics of the non-autonomous two mixed solitons can be understood after transforming the two-soliton solution [see Eqs. (B.3)-(B.13)] by using the transformations given in Appendix A [see Eqs. (6gr) and (6go)]. In the following analysis, the collision dynamics for two interesting forms of the nonlinearity coefficient , namely kink-like and the periodically modulated nonlinearity are considered in detail.

### (i) Collision dynamics in the presence of kink-like modulated nonlinearity

Figure 10 shows the two-soliton collision in the non-autonomous two-component GPEs (1) with kink-like nonlinearity, whose form is given by Eq. (6a) in the absence of Rabi terms, and the corresponding strength of the magnetic trap is given by Eq. (6c). The figure shows that the two solitons undergo a collision around = and get well separated for larger values of .

Next, we include the Rabi coupling and plot the non-autonomous two soliton collision in Fig. 11. We notice from Figs. 10 and 11 that the nature of soliton collision in the presence of Rabi term looks alike the soliton collision in the absence of Rabi term. Additionally, in the present case, there is an oscillating exchange of atoms between the two components. The background oscillates in a periodic manner, in which the oscillations in a given component is maximum while it is minimum in the other.

Here we present the asymptotic analysis of the non-autonomous two-soliton solution of two-component GPEs with kink-like nonlinearity, obtained by making use of (B.3)-(B.13) and the transformations (6go) and (6gr).
For the choice , , the asymptotic forms of the non-autonomous solitons (, ) well before and after collision can be expressed as below:

Before collision

(6gh) | |||||

(6gi) | |||||

After collision

(6gj) | |||||

(6gk) | |||||

where , , , , , and , where , and all other quantities are given in Appendix B [cf. (B.6)-(B.13)]. Here, () is the amplitude of the soliton in the th component, , before (after) collision, in the absence of Rabi-term and the inhomogeneity, which are defined in Appendix C.1.

We notice that in the absence of Rabi term , there is no oscillatory terms as expected. By computing the densities before and after collision and noticing (since , see Eqs. (C.1)-(C.4)), we identify that the nature of collision is elastic except for a phase-shift. This phase-shift is same for both solitons and that can be found as . There will be oscillations in bright-dark solitons and also in the background along with a modulation due to nonlinearity and Rabi term.

### (ii) Collision dynamics in the presence of periodically modulated nonlinearity

The collision of two solitons in the two-component GPEs with periodically modulated nonlinearity in the absence and presence of Rabi coupling are depicted in Figs. 12 and 13 respectively.

For this case the strength of the parabolic trap is given by Eq. (6gb).

The role of time modulated scattering length is to introduce periodic modulations in the soliton profile before and after collision uniformly. Meanwhile, the Rabi term leads to a periodic exchange of condensates between the components along with oscillations in the background. As in the previous example, here also the Rabi coupling makes it feasible to have both dark and bright parts of the mixed soliton in first and second components. We find that the solitons exhibit elastic collision even in the presence of Rabi coupling, though there is an exchange of atoms between the components. This example shows that the Rabi coupling does not affect the elastic nature of the collision in the two-component GPEs. An asymptotic analysis similar to that of kink-like nonlinearity can be carried out for this case too.

## 3 Non-autonomous BD soliton in three-component BECs and Rabi oscillations

Following our considerations for the binary BEC case, we now proceed with the investigation of the three-component system. The dynamics of three-component BECs in 1D with equal time-dependent interaction strengths (i.e., , j,l=1,2,3), and in the presence of external time-dependent harmonic potential , is governed by the following dimensionless non-autonomous three-coupled GPEs (see, e.g., [3]):

(6gl) |

where (x,t) are the wave functions of the condensates. As before, coefficients account for the (linear) Rabi coupling and are chosen to be equal (i.e., , =1,2,3); furthermore, the strength of the nonlinear coupling is given by .

As in the two-component case, we employ a rotation transformation and a similarity transformation (see details in Appendix A) and reduce Eq. (6gl) in the form:

(6gm) |

The above model, three-coupled NLS system (6gm), is also a completely integrable system (the three-component generalization of Manakov system with defocusing nonlinearity) and the soliton solutions can be obtained by various methods, e.g., the Inverse Scattering Transform method, the Hirota’s direct method, the Bäcklund transformation method, etc.

### 3.1 Non-autonomous mixed (bright-bright-dark) one-soliton solution

The explicit form of the mixed one-soliton solution of Eq. (6gm), in the form of a bright-bright-dark soliton, can be obtained by means of the Hirota’s method [40, 41]. Here, we consider a mixed soliton solution, with the bright part appearing in and components, and with the dark part in component. The exact analytical form of this type of mixed one- and two-solitons are given in Appendix B.2. Then we make use of the transformations, corresponding to three-component case given in Appendix. A, for obtaining the exact non-autonomous soliton solutions of (6gl), as in the two-component GPE system.

We again consider the same two forms for which were discussed in the previous section during our analysis of two-component condensates.

(a) Kink-like nonlinearity

Figure 14 shows the mixed one soliton solution of the non-autonomous three-component GPEs. (6gl) with kink-like nonlinearity. In this case, the shape of the soliton is affected significantly by increasing the value of as in the two-component case.

Also, we notice that the oscillations of the background are rapid as compared to that of two component case, even for the same value of Rabi coupling . This can be inferred by comparing the bottom panels of Fig. 14 and Fig. 3.
Thus, an increase in the number of components can make the exchange of atoms between soliton and its background is too faster.

(b) Periodically modulated nonlinearity

The propagation of three-component non-autonomous mixed one-soliton in the presence of periodically modulated nonlinearity having the form of Mathieu function [see Eq. (6ga)] is shown in figure 15. As in the two component case, here also occurs oscillatory transfer of atoms among the components accompanied by oscillations of the background due to Rabi coupling.

The above two examples show that the role of Rabi coupling on three-component non-autonomous BECs is similar to that of two-component BECs except for an increase in the oscillations due to the additional component. Apart from this, we would like to give emphasis to a particular dynamical feature of the three-component system, in the presence of Rabi coupling, which is not possible in the two-component case. It can be recalled from the study on two-component non-autonomous case that it is impossible to make any component (bright/dark) of the mixed soliton to vanish completely, in the presence —as well as in the absence— of Rabi coupling. Contrary, in the three-component system, it is possible to make any one of the bright parts of the mixed soliton to vanish completely before the introduction of Rabi coupling term. The additional component brings an additional arbitrariness to the non-autonomous GPE system (6gl), which allows to make the density of any one of the bright parts of the mixed soliton to be zero before the incorporation of the Rabi term. The role of the Rabi term is to switch the condensates between all three-components. So, in the presence of Rabi term, some amount of field will be transformed to the component where there was no field before its introduction. Thus, in three-component condensates one can transfer a significant part of the condensates to a component which admits no field at all before the introduction of Rabi coupling from the other components. This is shown in Figs. 16 and 17, in which the condensate is completely absent in the second component in Fig. 16 but, after introducing the Rabi coupling, appreciable portion of condensate appears in the same component as shown in Fig. 17.

### 3.2 Non-autonomous mixed (bright-bright-dark) two-soliton collision

The explicit form of mixed two-soliton solution of (6gm) with the bright parts appearing in components and the dark part in , obtained by Hirota’s bilinearization method is given in Appendix B.2 and the corresponding asymptotic analysis is presented in Appendix C.2.

From the asymptotic expressions (C.5) and (C.7) one can find that the amplitudes (condensate densities) of the colliding solitons before and after collision can be related

where and , , where , , , and are defined in Appendix B.2. It is instructive to notice that , are not unimodular, in general, and also depend on dark soliton parameters and in addition to the bright soliton parameters. This will result in the energy sharing collision displaying suppression (enhancement) of condensate density in the bright part of a given soliton with commensurate changes in the bright part of the other colliding soliton. However, one can have standard elastic collision for the choice = , for which , , become unimodular. Additionally, the colliding solitons also experience a phase-shift = = .

It can also be inferred from the asymptotic expressions (C.6) and (C.8) that the dark solitons in the component always exhibit elastic collision as , . These dark solitons also undergo a phase-shift, same as that of bright solitons.

This study can be straightforwardly extended to integrable -component CNLS system with defocusing nonlinearity, for . It is interesting to note that the above kind of energy sharing collision in such -component CNLS system can take place only if the bright part of the mixed soliton appears at least in two-components.

Such fascinating energy sharing collision in the three-component GPE system (6gm), is shown in Fig. 18.

In Fig. 18, the bright part of the mixed soliton is enhanced in the component whereas it
is suppressed in the component after its collision with the soliton . The reverse scenario
takes place for the bright parts of the mixed soliton .
Notice that the dark parts of mixed solitons and are unaffected by the collision. The
energy sharing collision is characterized by an exchange of condensates among the bright parts of the
colliding mixed solitons and , leaving the dark part unaltered after collision.
This type of energy sharing collision in different integrable CNLS systems appearing in nonlinear
optics has been extensively studied in references [28, 40, 29, 39].

The non-autonomous two-soliton solution can be expressed as

(6gn) |

where are given in Eqs. (B.18-B.19) in which X and T are defined by equations (A.6) and (A.7), respectively. Now, it is of further interest to study whether the above discussed energy sharing collision still prevails in the presence of time-varying nonlinearity and Rabi coupling. For illustrative purpose, again we consider the two types of the time-dependent nonlinearities discussed in the two-component case.

### (i) Energy sharing collision of mixed solitons in non-autonomous three-component GPE.

### (a) Kink-like nonlinearity

The soliton collision for a kink-like modulated nonlinearity in the context of the non-autonomous three-component GPEs. (6gl), in the presence of Rabi term, is depicted in Fig. 19. The density of the bright part of the mixed soliton gets enhanced and is suppressed after the collision in the component, whereas the component experiences reverse effects for and , and the dark parts of mixed solitons in component remain unaltered. Thus, the switching nature of energy sharing collision in the autonomous system (6gm) is unaffected by the presence of Rabi term in the non-autonomous GP system, for this choice of time-dependent nonlinearity and external potential. This shows that the energy sharing collision of bright parts of the mixed solitons can exist for a wide range of time-dependence of nonlinearity, external potential and Rabi coupling, for which the non-autonomous system (6gl) is integrable. An asymptotic analysis of the non-autonomous two-soliton solution (6gn) can be carried out as in the two-component case. One can notice that the energy sharing collision still prevails as , , in general, which ultimately leads to . The phase-shift can be found to be the same for both the solitons and is given by , where , and are defined in Appendix B.

In the two-component case, we observe a gradual increase in the density of condensates in the bright parts in the region = to and it remains constant for , due to the nature of nonlinearity (Fig. 10) and onset of Rabi oscillations due to Rabi coupling (see Fig. 11). Note that, there the collision is elastic in the absence as well as in the presence of Rabi coupling and is not an energy sharing one. Thus, in two-component BECs the growth of condensate is purely due to the nature of .

However, in the present case, for the mixed soliton there is a suppression in the density of its bright part in after interaction due to the energy sharing collision. This suppression balances the enhancement in its density in this regime due to the nature of the time-varying nonlinearity, and will result in standard soliton with constant amplitude after interaction (see Fig. 19). Also, the Rabi term leads to beating effects in this soliton in , but with a smaller period as compared to that of the two-component case. However, the opposite effect takes place for in the second component . This shows that energy sharing collision can be profitably used in altering the exchange of condensates through Rabi switching.

### (b) Periodic modulated nonlinearity

Energy sharing collision of mixed solitons in the three-component GPEs with periodically modulated nonlinearity is shown in Fig. 20. We observe that in this case also the energy sharing collision for the non-autonomous system (1) takes place as that of the autonomous system (2) (see Fig. 18).

We observe another dramatic switching of condensates during collision of two solitons in the non-autonomous three-component GPEs. (6gl) with periodically modulated time-dependent nonlinearity, for smaller values of Rabi term . Here, before collision, we have an oscillating bright soliton part in the component, which completely transforms to an oscillating mixed (bright-dark) soliton in the same component after collision. In the second component, , the amplitude of bright soliton part gets enhanced after the collision. For the soliton , in component enhancement (suppression) with periodic oscillations takes place. In the component the solitons undergo elastic collision and there occurs only switching of condensates due to Rabi coupling. Notably, in component the pure bright part of is transformed to mixed (BD) part after collision (see component in the bottom panel of Fig. 20). This interesting collision is a consequence of the combined effects of time-dependent nonlinearity and external potential, Rabi term and the exchange of condensates due to energy sharing collision. Another noticeable effect arising due to the Rabi coupling, particularly for this type of periodic nonlinearity, is an increase in the relative separation distance between the solitons well before and after collision. Furthermore, the collision takes place faster in the presence of Rabi term.

## 4 Conclusions

In this work, we have studied the dynamics of non-autonomous mixed bright-dark matter-wave solitons in two- and three-component Bose-Einstein condensates. Our setting included a time-dependent parabolic potential and scattering length, as well as Rabi coupling between separate hyperfine states. We have transformed the non-autonomous two- and three-component Gross-Pitaevskii equations, into a set of integrable defocusing two- and three-component Manakov autonomous systems by means of two successive transformations. These transformations can be viewed as a rotation followed by a similarity transformation. Then, with the aid of the two- and three-component soliton solutions of defocusing Manakov systems, and by inverting the transformations, we have studied the dynamics of single- and multiple-non-autonomous matter-wave solitons. Our considerations involved two different choices of nonlinearities namely, a kink-like and Mathieu-like ones, together with their corresponding time-varying harmonic potentials.

Our study on the two-component case has shown that the nature of soliton propagation depends significantly on the nature of nonlinearity and the associated potential. The switching depends purely on the Rabi coupling. Here, the bright-dark matter-wave solitons exhibit elastic collisions along with oscillations due to Rabi coupling, and profile modulation due to the nonlinearity and strength of the potential. In this two-component non-autonomous case, the dark-bright solitons exist as “symbiotic” structures, i.e., the bright component exist only due to the presence of its dark counterpart. Hence it is impossible to make any part (either bright or dark) to be zero. Also, Rabi coupling makes it possible to have a special type of mixed soliton in which bright and dark parts can co-exist in the same component. Such structures are reminiscent to the so-called beating dark-dark solitons, which have recently been observed in autonomous binary BECs [14, 42]. The present study shows the possibility of observing such structures, with co-existing bright and dark parts in the same component, also in non-autonomous multi-component systems in the presence of Rabi coupling. Our analytical results of the two-component non-autonomous GPE system are found to be in very good agreement with results of direct numerical simulations.

We have also studied non-autonomous three-component BECs. An important observation in this case is that the bright part of the mixed soliton can be absent in any of the three-component condensate in the absence of Rabi coupling. The introduction of Rabi coupling makes it feasible to have some portion of mixed-soliton in that component. Also, in the three-component BECs the number of oscillations increases as compared to that of two-component case and ultimately the exchange of atoms between the soliton and background also increases in a periodic manner.

Finally, we considered the collision of non-autonomous solitons in the three-component condensate. In this case, the non-autonomous matter-wave solitons undergo interesting energy-sharing collision, leading to an exchange of atoms, in addition to Rabi coupling, in the bright ( and ) components; notice that in the third component () the solitons undergo elastic collision. The matter-wave mixed (bright-bright-dark) soliton collision in non-autonomous three-component condensates also displays another fascinating feature, namely the reversal of the type of mixed- soliton part (i.e., bright type to dark type or vice-versa).

Our study can be extended to other types of forms for nonlinearity modulation by identifying the corresponding modulation of the harmonic potential using the Riccati Eq. (6gv) and vice- versa. The agreement between the numerical and analytical results indicates that the presented exact non-autonomous solution can well be utilised for the purpose of studying non-integrable multi- component GPE type systems with Rabi coupling. It is of future interest to investigate the effect of spatial and spatio-temporal modulations of nonlinearity. Works along these directions are in progress. This study may have ramifications in the design of matter-wave devices (e.g., switches) and also in experiments involving growth and oscillations of condensates.

## Acknowledgments

The work of T.K. is supported by Department of Science and Technology (DST), Government of India, in the form of a major research project. R.B.M. acknowledges the financial support from DST in the form of Project Assistant. T.K. and R.B.M. also thank the principal and management of Bishop Heber College for constant support and encouragement. The work of F.T., H.E.N., and D.J.F. was partially supported by the Special Account for Research Grants of the University of Athens.

## Appendix A Mapping non-autonomous GPEs to integrable CNLS systems

In this Appendix, we explicitly show that the non-autonomous two- and three-component GP systems [cf. Eqs. (1) and (6gl)] with Rabi coupling can be transformed to integrable CNLS systems, i.e., the Manakov model and its three component generalization, with defocusing nonlinearities; this will be done by means of two successive transformations. In particular, in the case of the two-component BEC, first we apply the following unitary transformation to Eq. (1) with and :

(6go) |

In the case of the three-component BECs, a similar transformation is applied to Eq. (6gl); this transformation is of the form [32, 30]:

This way, we obtain the following set of non-autonomous equations without the Rabi coupling term:

(6gq) |

where and for Eq. (1), while and for Eq. (6gl). Then, by performing the similarity transformation

(6gr) |

where

(6gs) | |||

(6gt) | |||

(6gu) |

with and being arbitrary real constants, Eq. (6gq) can be transformed into the system of Eq. (2) for the two-component system, or to Eq. (6gm) for the three-component system, with the condition

(6gv) |

where . Equation (6gv) is nothing but a Riccati-type equation. Notice that similar type of transformation has been reported for two-component BECs [26] but in the absence of Rabi coupling.

## Appendix B One- and two-soliton solutions of the integrable two- and three- component defocusing CNLS systems

### b.1 Bright-dark solitons in the two-component system (2)

The bright-dark one-soliton solution of Eq. (2), with bright (dark) part appearing in the () component, obtained by Hirota’s bilinearization method [27, 39].

(6gh) | |||

(6gi) |

where , = -, , , and . The bright part of the mixed soliton is characterized by three complex parameters , , and , and one real parameter . Note that does not affect the amplitude of the soliton. The dark soliton part of the mixed soliton also influences the bright part through the parameters and . Also, the solution becomes non-singular only for the choice