Dynamics of the Creation of a Rotating Bose–Einstein Condensation by Two Photon Raman Transition Using a Laguerre–Gaussian Laser Pulse

: We present numerical simulations to unravel the dynamics associated with the creation of a vortex in a Bose–Einstein condensate (BEC), from another nonrotating BEC using two-photon Raman transition with Gaussian (G) and Laguerre–Gaussian (LG) laser pulses. In particular, we consider BEC of Rb atoms at their hyperﬁne ground states conﬁned in a quasi two dimensional harmonic trap. Optical dipole potentials created by G and LG laser pulses modify the harmonic trap in such a way that density patterns of the condensates during the Raman transition process depend on the sign of the generated vortex. We investigate the role played by the Raman coupling parameter manifested through dimensionless peak Rabi frequency and intercomponent interaction on the dynamics during the population transfer process and on the ﬁnal population of the rotating condensate. During the Raman transition process, the two BECs tend to have larger overlap with each other for stronger intercomponent interaction strength.

The helical phase front of Laguerre-Gaussian (LG) laser beams has been associated with its orbital angular momentum (OAM) in the paraxial regime [55]. A photon of such an LG laser modes has phase profile e ilφ , and carries lh unit OAM in the transverse plane, where φ is the angular coordinate and l is an integer, known as the winding number of the beam. Such LG modes are known to transfer OAM from an optical field to the Rydberg atom [56], BECs [57][58][59][60], and to create a mechanical rotation of particles [61,62]. It was shown that a coherent coupling between the ground state of condensate with a rotating condensate in vortex state, can be achieved by the transfer of OAM of photons to the condensed atoms through Raman transitions [37]. Quantum dynamics of such vortex coupler using LG beam was studied, and an off-axis motion of the quantized vortex cores was interpreted as the collapse and revival of the atoms of the condensate [63]. Besides, a pair of LG laser modes with unequal phase windings couple internal atomic states of BEC through Raman transitions, and thus giving rise to spin and orbital angular momentum coupling in the ground states of a spinor BEC [64,65]. Moreover, it has been shown that almost all the atoms in the non-rotating BEC can be transferred to the BEC with vortex, by employing LG beams [66,67].
Although an impressive volume of literature has been devoted to this subject, few of its vital aspects remain further to be explored. One such aspect constitutes the role played by the interaction between two BEC components on the population transfer. Indeed, during the transfer process, atoms of two condensates are present in two different hyperfine states, one with vorticity and another without vorticity. Thus, not only the atom-laser coupling, but also the atom-atom interaction between two different components is expected to influence the population transfer process. Note that the focus of the majority of the previous studies has been on the complete particle transfer from one quantum state to another. However, it is expected that by maneuvering atom-light coupling and intercomponent interaction one could achieve a population transfer of any desired value. In this way it equips us to realize a binary-mixture where one component contains a vortex, and the other does not, thus emulating the so-called vortex-bright-soliton structure [16,68]. Additionally, it is also desirable to know, through the miscibility parameter [15,69,70], how atoms in the condensate with a vortex penetrate into atoms of the condensate without any vortex during the transfer process [15,69,70]. Therefore, motivated by experimental accessibility [71,72] and theoretical novelty of the problem, we theoretically address these important aspects of the transfer mechanism in this paper.
We investigate the dynamics of population transfer from a nonrotating BEC to a Raman coupled rotating BEC by employing LG and Gaussian (G) pulses. In this process, the atoms in rotating BEC gain angular momentum from the LG laser pulse. We consider pulsed G and LG beams as the pump and Stokes beams , respectively, to transfer the atoms from one hyperfine level to another. In particular, we choose the temporal width of the pulses to be in the same time scale determined by the trap frequency. This consideration provides us the framework to understand the dynamics during the transfer process. Numerically integrating the Raman coupled multicomponent Gross-Pitaevskii equations, we point out the following key points: (i) the sign of the vorticity of the condensate as well as the initial growth region of the vortex state, captured within the density patterns, depend upon which laser mode is chosen as pump or Stokes beam and (ii) the repulsive inter-component atomic interaction and peak Rabi frequency of laser beams determine the number of atoms transferred to the non-rotating BEC. By calculating the overlap integral between the two condensates we also quantify how two condensates penetrate into each other during the transfer process.
We have organized the remainder of this paper as follows. In Section 2 we describe the theory of transfer mechanism. Section 3 provides a brief description of the numerical schemes used. In Section 4 we present results of the complete particle transfer and effects of the inter-particle interactions and the Raman coupling parameter on the final population of the rotating BEC. In Section 5, we discuss the implication and possible future extensions associated with the results presented. Appendix A presents the effects of trap frequencies and the time-delay between the pulses on the population transfer. Finally, in Appendix B we briefly outline the Hamiltonian and the derivation of the equations of motions.

Theoretical Methods
In our study, we consider BEC of alkali atoms trapped in a quasi two-dimensional harmonic trap confined in the x − y plane with z axis being the quantization axis. In order to transfer OAM from the optical beam to the BEC, we consider three electronic levels of the alkali atoms are coupled by a pair of laser pulses in Λ-type configuration as shown in Figure 1. Atoms of initially prepared BEC are at the state |1 , one of the hyperfine levels of the electronic ground state of atoms. The state |3 is an intermediate non-resonant excited state. The final state is considered to be |2 , another hyperfine level of the electronic ground state of the atoms. The atoms are irradiated by two laser pulses propagating collinearly parallel to the quantization axis [73]. We remark that with the dipole approximation of the atomic transitions, the changes in the internal spin states of atoms are dictated by polarizations of two light fields. However, the changes in external orbital motion of the atoms of BEC around the quantization axis are determined by the difference of the orbital angular momentum (OAM) of two light fields [57]. Let us consider that the OAM of the twisted laser pulses for the transition from state |1 to state |3 is l 1 and for |3 to |2 transition is l 2 . Then, the electric field vectors involved in this absorption or emission transitions can be written as (for i = 1 and 2) where E i (t),ˆ i , k i and ω i are the corresponding time dependent amplitude profile, polarization vector, wave number and frequency of the i-th pulse, respectively. We consider the temporal amplitude profiles of both pulses have the same form [74]: where τ 1(2) is the temporal position of the peak value of electric field E 1 (2) . Maximum amplitude E max and pulse duration T are the same for both pulses. The optical absorptionemission cycle imparts OAM onto the atoms in final state |2 and creates a vortex in the BEC with charge (l 1 − l 2 ) unit. Because of collinearity of the E 1 and E 2 pulses, no additional linear momentum is generated in the final state. In addition to such two-photon transitions in atomic BEC, these lasers also create extra confining potential, namely optical dipole potentials for the atoms in the states |1 and |2 [75]. In practice the value of detuning Δ is large, which ensures the negligible populations in state |3 . This allows us to eliminate the state |3 adiabatically. During the transfer process, atoms are present in both the hyperfine states, |1 and |2 . Therefore, coherent evolution of the condensates of atoms in these two states, characterized by wavefunctions Ψ 1 (x, y, t) and Ψ 2 (x, y, t) respectively, are governed by two Raman coupled Gross-Pitaevskii equations (see Appendix B for the derivation) and, where is maximum light intensity of both pulses and d is the atomic transition dipole moment and w 1(2) is the beam waist of the corresponding laser pulse. Therefore, the effective trap potentials felt by atoms of the condensates are We derive Equations (3) and (4) by nondimensionalizing Equations (A11) and (A12) respectively. For this, we scale the spatial coordinates by oscillator length a osc = √h /mω, time by 1/ω and condensate wavefunctions by N/a 3 osc . Here, m is the mass of the atoms and N is the total number of atoms in the system, and ω is the trapping frequencies along x and y directions of the harmonic trap. We denote N 1 and N 2 as the number of atoms in condensates Ψ 1 and Ψ 2 respectively, and consider the total number, N = N 1 + N 2 , is conserved during and after the transfer process. We point out that initially N 1 = N and N 2 = 0. Note that, the parameter associated with the peak Rabi frequency, V max , contains parameters from the considered atomic transition, laser pulses and the trap of the condensate. The quasi-2D configuration of the trap is achieved by ensuring large trapping frequency in z direction, that is, ω z ω. The intra and inter-component coupling strengths are G ii = 2N √ 2πλa ii /a osc and G ij = G ji = 2N √ 2πλa ij /a osc , respectively, and λ = ω z /ω is the anisotropy parameter. The intra-component and inter-component scattering lengths are denoted by a ii and a ij , respectively. Initially, only the condensate Ψ 1 is present within the trap. With two photon Raman transitions, the condensate Ψ 2 grows by gaining atoms from the condensate Ψ 1 . During this process atoms in Ψ 2 gain (l 1 − l 2 ) unit orbital angular momentum, which is manifested as a phase factor e i(l 1 −l 2 )φ in the condensate wavefunction Ψ 2 . The phase factor e −i(l 1 −l 2 )φ in the coupling term of Equation (3) ensures that no angular momentum is transferred back to the atoms in condensate Ψ 1 . Transfer of this angular momentum to the condensate Ψ 2 results in generating quantized vortex in the condensate. A quantized vortex in a BEC is point like topological defect which is manifested in the phase profile of the condensate wavefunction Ψ 2 . Around the vortex the phase of the condensate wavefunction changes by κ × 2π, where κ is an integer, which is referred to as the winding number or charge of the vortex.
A system of two component BECs can exhibit two phases, miscible or immiscible, depending on the the strengths of intracomponent and intercomponent interactions. At zero temperature, two defect free condensates in a homogeneous trap are miscible when a 2 12 ≤ a 11 a 22 , and immiscible for a 2 12 ≥ a 11 a 22 [76]. However, these conditions are modified when the condensates are considered in inhomogeneous trap [77]. Effects of finite temperature [70] and topological defects [15] on the miscible-immiscible transition have been reported. A well known measure to characterize these two phases is the overlap integral defined as [15,69,70] where n 1(2) (x, y) = Ψ 1(2) (x, y) 2 are the densities of the condensates. Λ = 0 corresponds that the two condensates are spatially separated, that is, the system is in immiscible phase. Whereas, Λ = 1 implies maximal spatial overlap between the condensates, that is, the system is in complete miscible phase.
To this end, by utilizing two-photon Raman transition, we transfer the atoms from one initially populated quantum state to another unpopulated state via an intermediate state, see Figure 1. A pump field links state |1 to electronically excited state |3 , and Stokes field links state |3 to another low energy state |2 . We perform a one-way controlled particle transfer from one hyperfine state to another, where atoms of the daughter state carries one unit vortcity, either positive or negative, thus enabling us to create a rotating BEC. In this context, coherent population transfer is possible if the Stokes field precedes, but temporally overlaps with, the pump field, and the pulses are applied adiabatically.

Numerical Methods
We start with a BEC of N atoms at state |1 , in the absence of laser pulses. Therefore, we set terms associated with laser pulses in Equation (3) to be zero to obtain the initial solution. Then, the wavefunction of the initial BEC, Ψ 1 , is generated by solving Equation (3) in imaginary time using split-time Crank-Nicolson method [78]. The initial wavefunction of BEC of the atoms in state |2 , Ψ 2 , is considered to be zero. Using these two initial wave functions, we evolve the system in presence of laser pulses. For this, we solve the coupled GP equations in Equations (3) and (4) in real time. The phase imprinting in the Ψ 2 occurs dynamically due to the two photon Raman transitions, which is obtained by considering, and Since Ψ 2 is zero at the initial time t 0 , l 1 − l 2 unit vortex is imprinted on Ψ 2 at t 1 = t 0 + δt and vorticity of Ψ 1 remains zero. This transfer of angular momentum continues, as long as both pulses are present. However, since the process is one-way, it stops when all the atoms in condensate Ψ 1 are transferred to the rotating condensate. For simulations, we choose a square grid of 300 × 300 grid points with a grid spacing δ x = δ y = 0.05a osc and time step Δ t = 0.0001ω −1 . In our study, we consider hyperfine states of 87 Rb with |1, −1 as |1 and |2, +1 as |2 . The intracomponent scattering lengths a 11 and a 22 of these two states are 100.4a 0 and 95.44a 0 [79] respectively, where a 0 is the Bohr radius. The trap frequency ω = 2π × 30.832 Hz [80] and the anisotropy parameter λ = 40 are the same for both condensates. For this system the oscillator length a osc = 1.94 μm. Furthermore, the relation μ 1 (2) hω z holds throughout the time evolution indicating that quasi-2D configuration is maintained always. Total number of atoms in the system is N = 10 4 . To create a BEC with a vortex of charge −1 unit, we use G pulse as "pump" of which l 1 = 0, and LG Pulse as "Stokes" with l 2 = 1. If we interchange the "pump" and "Stokes" laser pulses, a vortex of charge +1 unit will be created in the BEC. For simulations, we use the pulses with same temporal duration of T = 4.9 ms.

Creation of Vortex in the BEC
In G-LG pulse sequence, we employ G pulse as pump and LG pulse as Stokes, for which l 1 = 0 and l 2 = 1 respectively. For this arrangement, we consider τ 1 = 1.4 and τ 2 = 1.0 in the units of 1/ω. During the Raman transitions of atoms from state |1 to |2 an amount of −1 unit OAM is transferred to the atoms in state |2 . Here, we describe the transfer process. First, a photon from the G laser pulse which has zero OAM is absorbed by the atom in |1 . As a result, the atom is excited to an intermediate excited state |3 . Then, a photon with +1 unit OAM is emitted by the atom at the state |3 onto the LG beam. After this emission process the atom comes back to another ground state |2 . The conservation of the total angular momentum of the system, that is, the total angular momentum of atom plus light pulses, ensures that atom at the state |2 gains −1 unit OAM. Thus, −1 unit vorticity is created in the condensate Ψ 2 . Similarly, +1 unit vorticity can be created in the condensate Ψ 2 through LG-G pulse sequence, where we use LG pulse as pump and G pulse as Stokes of which l 1 = 1 and l 2 = 0 respectively.

Density Evolution of the Condensates
We have discussed that the sign of the vorticity in condensate Ψ 2 depends on the laser modes chosen as pump and Stokes beam. Here, we point out how the sign of the vorticity can be inferred from the changes of density patterns of the condensates during the transfer process. Furthermore, these density patterns serve as promising candidates to elucidate the residual excitations created during the light-matter interaction, since these excitations leave their foot-prints on the density profiles creating additional humps and dips [81,82]. Figure 2(a 1 -a 6 ,b 1 -b 6 ) illustrate the density profiles of the condensates during the Raman transitions, when the vortex of charge −1 unit is generated in the condensate Ψ 2 . Whereas, Figure 2(c 1 -c 6 ,d 1 -d 6 ) illustrate the density profiles when +1 unit vortex is created. In the lower-left corner of each density profile we mention the fraction of atoms in the condensate with respect total number of atoms in the system.
From the comparison between the Figure 2(a 1 -a 6 ,b 1 -b 6 ,c 1 -c 6 ,d 1 -d 6 ), it is evident that density patterns of the condensates during the creation of −1 unit vortex are different from the case of creation of +1 unit vortex. During the initial growth of the condensate Ψ 2 , the atoms occupy the central region of the trap when the vortex of charge −1 unit is created, whereas the atoms occupy the peripheral region of the trap when the +1 unit vortex created. At t = 0, the laser pulses are absent and the condensate Ψ 1 is populated by all the atoms in the system, hence, the condensate Ψ 2 is empty. It is important to mention that for the coherent population transfer, we apply the Stokes beam first. Therefore, in the early stage of the dynamics population of condensate Ψ 2 remains zero. Once the pump beam is applied, condensate Ψ 2 starts growing at the expense of atoms being transferred from the condensate Ψ 1 . At the same time, a vortex of either −1 or +1 unit is imprinted on condensate Ψ 2 depending on the angular momenta of the pump and Stokes beams. For the case of LG-G pulse sequence, which is illustrated in Figure 2(c 1 -c 6 ,d 1 -d 6 ), we observe that 11% of atoms has been transferred in the first 10.84 ms, but 68% of atoms are transferred in the next 1.96 ms. In contrast to this, we observe that fewer numbers of atoms are transferred to condensate Ψ 2 at the same time instants when compared to the G-LG pulse sequence, which is also evident from Figure 2(b 1 -b 6 ). In both the cases, the generated vortex appears with core, that is, zero density region at the center of condensate Ψ 2 , which is visible in the density profiles of Ψ 2 shown in Figure 2(b 1 -b 6 ,d 1 -d 6 ). It is worth noting that density depleted region at the center of the trap is also observed in the density profiles of condensate Ψ 1 during the creation of −1 unit vortex in Ψ 2 (Figure 2(a 1 -a 6 )). However, such a hole is absent in the condensate Ψ 1 , when +1 unit vortex is created. To understand the nature of the density depleted regions, we study the phase profiles of the condensates. We confirm the presence of phase discontinuity at the center of condensate Ψ 2 for both the cases. It is mentioned earlier that the phase of the condensate wavefunction changes by κ × 2π around a quantized vortex, where κ is the winding-number or charge of the vortex. We compute the winding number κ to be −1 when we use G as pulse and LG as Stokes beam, whereas κ = +1 when we consider LG-G pulse sequence. On the other hand, the phase profile of the condensate Ψ 1 does not possess phase discontinuity during the transfer process for both cases. Thus, the hole in condensate Ψ 1 which is generated during the application of G-LG pulse sequence, is not a vortex. Focusing our discussion on G-LG pulse sequence, we ascribe the presence of hole in condensate Ψ 1 to the distortion of harmonic trap potential by the optical dipole potential. In this case, the optical dipole potential is induced by the G laser pulse for the condensate Ψ 1 and by the LG laser pulse for the condensate Ψ 2 . Note that, at t = 0 ms the laser pulses are absent and the minimum of the harmonic oscillator occurs at the center of the trap. Hence, we obtain a pancake-shaped density profile of the condensate Ψ 1 , which has maximum density at the trap center to minimize trap potential energy. Then, during the application of laser pulse, the G-pulse gradually creates a rotationally symmetric "hump" at the center of the trap, which increases the potential energy at the trap center. Therefore, the minimum of the effective trap potential V eff,1 gets shifted radially away from the center, resulting in a rotationally symmetric annular region as the new minimum of the potential. It is important to mention that the density profile of a condensate in a binary mixture depends on the effective trap potential in conjunction with the number of atoms in the condensate, intra and intercomponent scattering length. Therefore, the atoms of the condensate Ψ 1 move away from center of the trap and settle at the annular region to minimize the trap potential energy. This creates a hole at the center of the density profile of the condensate Ψ 1 . Since the optical dipole potential induced by LG pulse has parabolic form around the center of the trap, the position of the minimum of the effective potential V eff,2 does not change over time. However, the steepness of this effective potential changes with time. It increases up to time t = τ 1 and then gradually decreases back to its initial value which is determined by the considered harmonic potential. Therefore, the atoms in the condensate Ψ 2 are always pushed towards the center of the trap to minimize trap potential energy. As a result, during the growth of Ψ 2 , the central region of the trap is occupied by the transferred atoms first, and then rest of the region is occupied.
For LG-G pulse sequence, laser modes of pump and Stokes beam are interchanged. Now the optical dipole potential is induced by the LG laser pulse for the condensate Ψ 1 and by the G laser pulse for the condensate Ψ 2 . Therefore, with the increase of the steepness of the parabolic potential, which is generated by the LG pulse, the atoms in the condensate Ψ 1 are pushed towards the central region of the trap. The atoms that are transferred to condensate Ψ 2 experience the "hump" in the trap potential at the center, which is created by the G pulse. Thus, the atoms in condensate Ψ 2 are pushed towards an annular minimum region of the effective trap potential. This results in larger core of the vortex in condensate Ψ 2 during the transfer process, which is to be contrasted with the previous case.

Root-Mean-Square Radius of the Condensates
The growth rate of condensate Ψ 2 can be inferred from the rate of change of rms radii of the condensates. In Figure 3 we illustrate the evolution of the r rms of both condensates during the transfer process for the cases when G-LG and LG-G pulse sequences are considered. From the comparison between the considered cases, we can infer that the growth rate of the condensate Ψ 2 is faster in the case of LG-G pulse sequence than the case of G-LG pulse sequence. Note that, for the chosen pulses, the strength of the Raman interaction term V is always maximum, at the boundary of the trap. However, atoms in the condensate try to occupy the minimum of the trap potential to minimize the trap potential energy. In particular, the effective trap potential V eff,2 of condensate Ψ 2 has a minimum at the center of trap for G-LG pulse sequence, but at a distance r = w 0 ln 4V 2 (t)/w 2 0 /2 from the center, for LG-G pulse sequence. Therefore, in the later case, the minimum of the effective trap potential is closer to the trap boundary where the Raman coupling V term is maximum.
This suggests that the growth rate of the condensate Ψ 2 depends on the distance between the position of the minimum of effective trap potential and the position of maximum Raman coupling. After the transfer process, the rms radius of Ψ 2 oscillates around a mean value. The frequency of such residual radial oscillations, as can be seen from Figure 3, is approximately ω = ω/3 for both pulse sequences. The amplitude of oscillation is much smaller than the mean radius of condensate. Most importantly, such a small amplitude of oscillation indicates that negligible amount residual excitations have been activated during the population transfer. However, detail analysis of such excitations is out of the scope of this work.

Effects of Intercomponent Interaction
We now discuss the effects of intercomponent interaction between the two condensates, during the transfer process and the final population of the condensate Ψ 2 . We consider the LG-G pulse sequence as the representative example.
The scattering length a 12 , which quantifies interactions between the atoms of the two different components, plays an important role in determining spatial wavefunctions and the energy of the condensates. Indeed, for certain temporal duration of pulses and intercomponent scattering length, the strength of the atom-light interaction V max have to be monitored to get the desired population of atoms in the state |2 . In the Figure 4, we present the number of atoms in condensate Ψ 2 at the end of the transfer process as a function of a 12 and V max . We vary peak Rabi frequency V max from 1 to 100 and intercomponent atomic scattering length a 12 from 70a 0 to 110a 0 . Peak Rabi frequency can be controlled either by changing peak light intensity of the pulse or by changing the detuning. Whereas, the scattering length can be varied through the magnetic Feshbach resonance [83]. We observe complete population transfer from condensate Ψ 1 to condensate Ψ 2 when V max is greater than 100 (not shown in the diagram). Intercomponent interactions merely affect the transfer process. In this region, atom-light interaction is strong enough to affect any density distribution determined by a 12 . However, this situation does not hold for intermediate values of V max , predominantly between 100 and 10. In this region, stronger is the intercomponent interaction, larger is the number of atoms transferred to condensate Ψ 2 . For small values of V max , larger values of a 12 suppresses the transfer process,which is evident from Figure 4. It is important to mention that in this limit, we observe the growth of condensate Ψ 2 is different for different values of intercomponent atomic scattering length. That is, depending on the strength of the atom-light interaction, a 12 affects the population transfer in different manner. For example, for V max = 1, the final population of Ψ 2 is suppressed for larger a 12 , whereas, for V max = 10, strong interaction increases the population in Ψ 2 (see Figure 5a,b). The colorbar shows the fraction of atoms in condensate Ψ 2 with respect to the total number of atoms N = 10 4 in the system at the end of the transfer process. The population transfer from Ψ 1 to Ψ 2 is carried-out using LG(pump)-G(stokes) pulse sequence. The system is confined in a harmonic trapping potential with frequency ω = 2π × 30.832 Hz and the intraspecies scattering lengths are a 11 = 100.04a 0 and a 22 = 95.44a 0 . In order to gain further intuition regarding the combined effect of the atom-light coupling and inter-component interaction we resort to the density evolution. Figure 6 presents few representative snapshots of the density profiles for V max = 1 and a 12 = 80a 0 . Note that this particular parameter set corresponds to 45% particle transfer to Ψ 2 (see Figure 5), therefore, enabling the creation of binary mixture of almost equal particles. This, in par-ticular, is better comprehended from the density profiles depicted in Figure 6. Preparing the initial state characterized by all the atoms residing at the state |1 ( Figure 6(a 1 ,b 1 )), the light-matter interaction is initiated by employing LG-and G beams as pump and stokes beam, respectively. Note that the large vortex core of Ψ 2 during the early stage of population transfer is caused by the presence of the Gaussian potential barrier at its center (see Figure 6(b 2 )). However, the same vortex core gradually shrinks as more number of particles are transferred and the Gaussian barrier gradually diminishes (see Figure 6(b 3 ,b 4 )). Finally, the BEC at |2 possess 45% of the total number of particles. A close inspection of Figure 6(a 4 -a 6 ,b 4 -b 6 ) reveals that a breathing motion characterized by expansion and contraction of the density profiles has been triggered in both Ψ 1 and Ψ 2 . Another observation is that Ψ 2 , after the population transfer, exhibits larger vortex core when compared to the case of complete particle transfer, see Figure 2(d 6 ) and Figure 6(b 6 ). In addition, we observe the peak Rabi frequency plays an important role in determining the miscibility between the two components during light-matter interaction. This is in contrast to the case when the light field is absent, that is, miscibility of two condensates is determined by the intra and intercomponent interactions. To illustrate this, we have considered Rabi frequencies, V max = 1 and V max = 10, for which both the condensates Ψ 1 and Ψ 2 have finite number of atoms N 1 and N 2 , even after the light-matter interaction. For these two cases, we show the variation of the miscibility parameter Λ with time in Figure 7. Note that just after the initiation of the transfer process, condensate Ψ 2 grows within the condensate Ψ 1 , resulting in gradual increase of Λ. However, when a sufficient number of atoms have been transferred to condensate Ψ 2 and both the pulses have significant temporal overlap, mutual repulsion between the condensates and the optical dipole potential tend to push the two condensates away from each other. This results in decrease of Λ. Again, the overlap between the condensates and hence Λ increases as pulses gradually die down. It is important to notice that during the light-matter interaction we obtain larger values of Λ for larger values of a 12 . This indicates, the stronger the intercomponent repulsion between the two condensates, the larger the overlap between them is. This is to be contrasted with the case when light-matter interaction is absent, in which, larger intercomponent repulsion separates the condensates spatially. After the light-matter interaction, that is, when the optical dipole potentials disappear, the miscibility between the condensates is determined by intra-and inter-component interactions.

Conclusions
In conclusion, we have shown that how two-photon Raman transition can be used to generate a rotating BEC with vorticity of either sign, by transferring atoms from another condensate. In this transition, atoms gain angular momentum from the LG laser pulse before being transferred to a rotating condensate. Density patterns of the condensates during the light-matter interaction depend on sign of the vorticity of the rotating condensate.
In particular, we have show that how a specific choice of pump and stokes beams can alter the locations, within the trap, where the transferred particles start accumulating. Notably, this result stems form the optical potentials felt by the atoms of each individual state. Most importantly, the interchange of pump and stokes laser beams not only changes the sign of vorticity, but also gives rise different dipole potentials influencing both the growth rate and growth region of the new Bose-Eintein condensate. The growth of the new condensate with −1 unit vorticity is started from the central region of the trap, but a condensate with +1 unit vorticity starts to grow from the peripheral region of the same. Furthermore, the "smoothness" of the density profiles and the temporal evolution of rms radius imply that there are very few excitations emerging during the light-matter interaction.
Moreover, we have shown that the number of transferred atoms can be monitored by tuning the intercomponent interaction, if the peak Rabi frequency of light-matter interaction is low and in particular, large intercomponent interaction subdues this transfer process. In this way, by maneuvering the atom-light interaction strength and the intercomponent scattering length one can create binary mixture of condensates.
Another major finding from our investigation is that intercomponent interaction kind of plays an opposite role in the process of phase separation during the Raman transition process, in contrast to literature [84] when such dynamical perturbation is absent. We find that a stronger intercomponent interaction favors greater miscibility between the condensates during the light-matter interaction.
Finally, we point out that the storage of a photon pair entangled in OAM space through Raman transition in the cold atomic ensemble has served as a sandbox to study information processing [85]. Besides, because atoms can have higher spin manifolds than light, the extension of our work to the spinor BEC would be an important study. Various topological properties can be developed in the ground state depending on Rabi frequency and atom-atom interaction strength, for example, a Mermin-Ho vortex or a meron pair phase [86], and might lead to the exhibition of non-Abelian braiding statistics [87] which is particularly interesting for topological quantum computing protocols [88]. We expect our study will shed light for further research in this direction.

Data Availability Statement:
All the data that support the plots within this paper are available from the corresponding author upon reasonable request.
Acknowledgments: K.M. is thankful to Subrata Das for technical assistance. K.M. is also grateful to Physical Research Laboratory, Ahmedabad for the hospitality during the initial stages of this work. S.B. and D.A. gratefully thank Arko Roy and Pekko Kuopanportti for insightful discussions.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A. Effects of Harmonic Trap and Time-Delay between the Pulses
Here we discuss the effects of the trapping frequencies of the harmonic potentials and the time-delay between two pulses during the particle transfer and on the final population of the component Ψ 2 . In particular, we are interested in the weak atom-light coupling regime where both components contain finite number of particles. To that purpose we set V max = 1, also the interspecies scattering length is a 12 = 80a 0 . Figure A1a presents time evolution of the fraction of total number of particles in Ψ 2 for different trapping frequencies ω. It is evident that consideration of a different trapping frequency only changes the time scale of the relevant phenomenology. Remarkably, it does not alter the final population of the Ψ 2 to a great extent. The onset and the growth rate of the population transfer are influenced by the trapping frequency. For instance, for ω/(2π) = 30.832 Hz, the transfer starts at earlier time instant and occurs at a faster rate when compared to the others.
Next we inspect the effect of varying time-delay between the two pulses on the population transfer in the weak light-matter coupling regime (see Figure A1b). We fix the peak location (τ 2 ) of the stokes pulse at τ 2 = 10.33 ms and vary the peak location (τ 1 ) of the pump pulse. We notice that maximum particle transfer is approximately 45% of the total population, which is achieved when τ 1 lies in the range 11.36-12.39 ms. However, when the time delay is too large (τ 1 = 15.49 ms) or too small (τ 1 = 10.36 ms) the number of particles in Ψ 2 decreases. Besides, when the time delay is very small, in other words temporal peaks of the two pulses are very close to each other, the growth of the population shows a very different behavior, see the blue curve in Figure A1b. The noticeable swell, during the light-matter interaction, in the curve corresponding to τ 1 = 10.36 ms can be related to the creation of significant amount of excitations in both components. These excitations stemming from the interactions with laser pulses also remain in both components after the interaction is over. We remark that such residual excitations leave their finger-prints onto the density profiles of the components (not shown here for the brevity), and bear the signatures of the break-down of coherent population transfer. ih Using (1) and (4) x 2 +y 2 w 2 1 (2) ) (A14) and Here the BEC is considered to be confined at z = 0 plane and ω 1 ≈ ω 2 .