quantized vortices are among the hallmarks of superfluidity , and have been observed in the gaseous bose - einstein condensates ( becs ) of @xmath0  @xcite and @xmath1  @xcite , where repulsive interatomic interactions play a crucial role in nucleation and stabilization of vortices . in the case of attractive interactions , however , creating vortices by stirring the condensate in a harmonic trap is difficult because the vortex state is energetically unfavorable against decaying into the non - vortex ground state .
a possible experimental scheme for creating a vortex state with attractive interactions is to prepare it in the repulsive regime and then switch the sign of the interaction to attractive by using , e.g. , the feshbach resonance  @xcite .
another possible scheme to create vortices in an attractive bec is to imprint topological phases on the condensate using either controlled interconversion between two components  @xcite or spin rotation by adiabatic change of the magnetic field  @xcite .
the aim of the present paper is to study the dynamical instabilities and ensuing dynamics of vortex states with attractive interactions in 3d axisymmetric traps .
there are two types of instabilities that cause a vortex state to decay : dynamical and thermodynamic instabilities .
dynamical instabilities , or modulational instabilities , set in when the bogoliubov spectrum acquires an imaginary part which causes an exponential growth in the corresponding mode .
vortex nucleation  @xcite , vortex exchange between two components  @xcite , splitting of a vortex with attractive interactions  @xcite , and disintegration of a doubly - quantized vortex with repulsive interactions  @xcite are all attributed to dynamical instabilities .
thermodynamic instabilities are caused by interactions between the condensate and the thermal cloud accompanied by dissipation of the energy and angular momentum .
reductions in the number of vortices  @xcite , vortex bending  @xcite , and vortex - lattice formation  @xcite are associated with thermodynamic instabilities .
the time scale of the latter instabilities at low temperatures is typically @xmath2 s or longer  @xcite .
in contrast , the time scale of dynamical instabilities can be much faster than that of thermodynamic instabilities since the former grows exponentially in time .
we shall therefore focus on dynamical instabilities in this paper .    in ref .
@xcite , pu _
_ reported that singly- and doubly - quantized vortices in 2d exhibit a series of dynamical instabilities when the strength of interaction exceeds the corresponding critical values .
we numerically solved the time - dependent gross - pitaevskii ( gp ) equation for a 2d system , and showed that the vortex splits into two clusters which revolve around the center of the trap , and undergo split - merge cycles or collapse  @xcite .
in the present paper , we examine such dynamics in full 3d systems and also study the dynamical instability of the doubly - quantized vortex , that was recently realized by the mit group  @xcite .
our primary findings are the following : +  ( 1 ) our 3d analysis reveals that the dynamical instabilities depend crucially on the trap geometry . for a singly - quantized vortex
, split - merge cycles can occur in a pancake - shaped trap , whereas the split fragments immediately collapse in a spherical trap . in a cigar - shaped trap
, there is a parameter range in which the vortex can collapse into a single cluster due to the dipole instability .
+  ( 2 ) the doubly - quantized vortex with attractive interactions is always unstable against disintegration of the vortex core .
if the strength of interaction is suddenly changed to within a certain range , the vortex splits into three clusters which undergo split - merge cycles before one or possibly more of them collapse .
+  ( 3 ) even when the lifetime of the condensate is very short , the phenomenon of vortex split can be observed by suddenly changing the interaction to a large attractive value .
the split of a doubly - quantized vortex is shown to be observable using a current experimental setup of the mit group .
+  ( 4 ) using a variational method , we find that both the split instability and the split - merge cycles can be understood by using only three basis functions .
+  ( 5 ) in the immediate vicinity of the critical strength of interaction for the vortex split , quantum fluctuations like the ones found in a 1d case  @xcite become prominent .
this paper is organized as follows .
section  [ s : formulation ] briefly reviews the gp and bogoliubov theories and discusses dynamical instabilities .
section  [ s : single ] studies dynamical instabilities of a singly - quantized vortex state , and sec .  [ s : double ] those of a doubly - quantized vortex state .
section  [ s : variational ] examines a variational method to understand the dynamical instabilities and split - merge phenomenon .
section  [ s : lowly ] shows that the split - merge phenomenon is accompanied by the emergence of a symmetry - broken low - lying state .
section  [ s : conclusion ] concludes this paper .
we consider a dilute - gas bec at zero temperature confined in an axisymmetric harmonic trap .
the dynamics of the condensate are described by the gp equation @xmath3 where @xmath4 is the atomic mass , @xmath5 with @xmath6 and @xmath7 being the radial and axial trap frequencies , and @xmath8 is the s - wave scattering length .
we normalize the length , time , energy , and wave function by @xmath9 , @xmath10 , @xmath11 , and @xmath12 , respectively , where @xmath13 is the number of atoms in the bec , and the wave function is normalized as @xmath14 . the gp equation ( [ gp ] ) is then written as @xmath15 where @xmath16 , @xmath17 , and @xmath18 .
the bogoliubov spectrum for a stationary state @xmath19 of the gp equation is obtained from the bogoliubov - de gennes equations    [ bogo ] @xmath20    where @xmath21 and @xmath22 is the chemical potential . for the axisymmetric vortex states
@xmath23 , each bogoliubov mode @xmath24 couples only to @xmath25 , where @xmath26 and @xmath27 are integers and @xmath28 is the azimuthal angle . in the following discussions
, we shall refer to the modes with @xmath29 , @xmath30 , and @xmath31 as the dipole , quadrupole , and octupole modes , respectively .
the stability of the system is characterized by the bogoliubov spectrum @xmath32 .
when the system is in the ground state , or in a metastable state for which there is an energy barrier that prevents the system from decaying into a lower energy state , all bogoliubov eigenenergies are real and positive .
an example of such a metastable state is a bec with an attractive interaction below the critical strength of interaction for collapse  @xcite . when there are negative eigenvalues in the bogoliubov spectrum , the system can thermodynamically decay into a lower energy state by dissipating energy to the environment .
when complex eigenvalues emerge in the bogoliubov spectrum , the corresponding modes grow exponentially in time and the system becomes dynamically unstable .
the time scale of such dynamical instability is given by the inverse of the imaginary part .
the difference between thermodynamic and dynamical instabilities is that the former requires energy dissipation , whereas the latter can occur via an energy - conserving process .
thus , the dynamical instability can be a dominant decay process at sufficiently low temperature in a high - vacuum chamber in which dissipation is negligible over a short time .
the time evolution of a system that deviates slightly from the stationary state @xmath19 of the gp equation can be described in terms of the bogoliubov modes @xmath33 and @xmath34 as  @xcite @xmath35,\ ] ] where @xmath36 s ( @xmath37 ) are small coefficients that describe the deviation .
when the eigenvalue of a bogoliubov mode becomes complex , the mode grows exponentially in time as    @xmath38 \nonumber \\ & = & e^{i ( q \phi - \mu t ) } \bigl\ { \tilde \psi_0 + e^{{\rm im}e_n t }    re}e_n t + \theta_n ) + ( \tilde u_n- \tilde v_n ) \sin(\delta m \phi - { \rm re}e_n t + \theta_n ) \bigr ] \bigr\},\end{aligned}\ ] ]    where we leave only the bogoliubov mode with the largest imaginary part @xmath39 , and define @xmath40 , @xmath41 , @xmath42 , and @xmath43 . since both @xmath44 and @xmath45 are solutions of eq .
( [ bogo ] ) , one of them satisfies @xmath46 .
equation  ( [ modegrow ] ) indicates that the @xmath27-fold density modulation grows exponentially in time while rotating at frequency @xmath47 .
note , however , that the above discussion , which is based on bogoliubov analysis , can be applied only to an initial stage of the growth of the unstable mode , and we must integrate gp equation ( [ gpn ] ) to follow the time evolution over a longer time .
we briefly describe the numerical procedure which will be used to solve eqs .
( [ gpn ] ) and ( [ bogo ] ) .
the stationary state @xmath19 of eq .
( [ gpn ] ) is obtained by the newton - raphson method  @xcite or the imaginary time method  @xcite in the restricted functional subspace of @xmath48 .
the bogoliubov spectra are then obtained by diagonalizing eq .
( [ bogo ] ) in each subspace of the angular momentum @xmath27 .
time evolution of the wave function is obtained by integrating eq .
( [ gpn ] ) using the crank - nicholson scheme  @xcite with a spatial discretization taken typically to be @xmath49 . in studying dynamical instabilities ,
we add a symmetry - breaking perturbation to the initial state such as @xmath50 where @xmath51 is an integer .
while dynamically unstable modes can grow from infinitesimal numerical noise over a sufficiently long time even when we start from an axisymmetric state , we add the perturbation ( [ perturb ] ) in order to simulate experimental noise in a controllable manner so that the results are reproducible .
this type of perturbation may be realized experimentally , e.g. , using off - resonant laser beams propagating along the trap axis .
the system is allowed to evolve in time until just before the collapse sets in , since numerical simulations of violent collapses and explosions  @xcite in 3d space present a formidable computational challenge .
in this section , we study the dynamical instabilities in a singly - quantized vortex state @xmath52 .
figure  [ f : l1 ] shows imaginary parts of the bogoliubov spectra as functions of the interaction parameter @xmath53 for spherical ( @xmath54 ) , pancake - shaped ( @xmath55 ) , and cigar - shaped ( @xmath56  @xcite ) traps .     for ( a ) @xmath54 ( spherical trap ) , ( b ) @xmath55 ( pancake - shaped ) , and ( c ) @xmath57 ( cigar - shaped ) .
the dotted lines show critical values @xmath58 below which the collapse of the vortex state is caused by the monopole instability . , width=317 ]    as seen from figs .
[ f : l1 ] ( a ) and ( b ) , for the spherical and pancake - shaped traps , the quadrupole mode ( @xmath59 ) first becomes unstable in the sense that the imaginary part @xmath60 appears for values of @xmath61 smaller than those of other modes .
a similar result is obtained in ref .
@xcite for a 2d system , which corresponds to the pancake - shaped trap with a large asymmetry parameter @xmath62  @xcite . for the cigar - shaped trap [ fig .
[ f : l1 ] ( c ) ] , on the other hand , the dipole mode ( @xmath63 ) first becomes unstable , which occurs for @xmath64  @xcite .
this suggests that when the strength of interaction @xmath61 is adiabatically increased , the dipole instability plays a dominant role of decay mechanisms for the cigar - shaped trap with @xmath64 ; otherwise the quadrupole mode is dominant .
figure  [ f : panl1 ] shows the time evolution of integrated column densities and phase profiles for @xmath55 .
( lower - left insets ) for @xmath65 and @xmath55 .
the initial state is the stationary singly - quantized vortex state for @xmath66 disturbed by perturbation ( [ perturb ] ) with @xmath67 .
the lower - right insets show gray - scale images of the phase of the integrated wave function @xmath68 , where the phase from @xmath69 to @xmath70 is represented in gray scale from black to white .
( f ) the integrated density @xmath71 at @xmath72 , where the interaction is switched to @xmath73 and the trap is turned off at @xmath74 .
the sensitivity of the imaging in ( f ) is 12 times higher than that in ( a)-(e ) .
the sizes of the images are @xmath75 in @xmath71 and the phase profiles , and @xmath76 in @xmath77 in ( a)-(e ) , and @xmath78 in ( f ) in units of @xmath79 .
, width=317 ]    the integrated column density is defined by @xmath80 and the phase of the integrated wave function @xmath81 is used to draw the phase profiles .
the main panels in fig .
[ f : panl1 ] show top views @xmath71 and the lower - left insets show side views @xmath77 .
the interaction parameter is taken to be @xmath65 , which exceeds the critical value for the quadrupole instability @xmath82 , but is below the critical value for the dipole instability @xmath83 .
the bogoliubov eigenvalue of the quadrupole mode at @xmath65 is @xmath84 .
we add the perturbation  ( [ perturb ] ) with @xmath67 to the initial state to trigger the quadrupolar dynamical instability . in fig .
[ f : panl1 ] ( b ) , we see that a density fluctuation in the quadrupole mode starts to grow with the lyapunov exponent @xmath85 and rotates at frequency @xmath86 in accordance with eq .
( [ modegrow ] ) .
the vortex then splits into two clusters which revolve around the center of the trap [ fig .
[ f : panl1 ] ( c ) ] , and subsequently these clusters unite to recover the original vortex [ fig .
[ f : panl1 ] ( e ) ] .
after this , the split - merge cycle repeats .
figure  [ f : panl1 ] ( f ) shows an expanded image of @xmath71 at @xmath72 , in which the interaction was switched from @xmath65 to @xmath73 and the trap was turned off at @xmath74 .
the fringe pattern is due to interference between two matter waves emanating from the two clusters shown in fig .
[ f : panl1 ] ( c ) .
the above results confirm that our findings in ref .
@xcite , which were obtained in a 2d system , hold true in a pancake - shaped trap and are therefore realizable experimentally .    when the vortex splits into @xmath51 clusters , each of these subsequently collapses if the effective strength of interaction per each cluster , i.e. , @xmath87 , exceeds the critical value @xmath88 for the collapse of a condensate without vortices , while the split - merge cycle occurs if @xmath87 is sufficiently smaller than @xmath89 .
in fact , in the case of @xmath55 in fig .
[ f : panl1 ] , @xmath90 is smaller than @xmath91 , and therefore the split clusters do not collapse , allowing the system to undergo split - merge cycles .
however , for the spherical trap , @xmath92 and @xmath93 , and hence the value @xmath94 at which the split occurs is always larger than @xmath89 , indicating that split clusters always collapse .
we performed numerical calculations for the spherical trap with @xmath61 slightly larger than @xmath95 , and confirmed that two clusters immediately collapse without undergoing a split - merge cycle ( as shown in figs .
4 ( d)-(f ) of ref .
the phenomenon , in which a vortex first splits into fragments and each fragment subsequently collapses , may be called `` fragmented collapse '' .
figure  [ f : mitl1 ] shows time evolution for the cigar - shaped trap with @xmath57 .
( upper panels ) , @xmath96 ( lower panels ) , and phase profiles ( insets ) for @xmath97 and @xmath57 .
the initial state is the stationary singly - quantized vortex state for @xmath98 disturbed by perturbation ( [ perturb ] ) with @xmath99 .
the sizes of the images are @xmath75 for @xmath71 and @xmath100 for @xmath77 in units of @xmath79 .
, width=317 ]    the interaction parameter is taken to be @xmath97 , at which only the dipole mode is unstable [ see fig .  [
f : l1 ] ( c ) ] .
the initial state is the stationary state for @xmath97 plus the small perturbation ( [ perturb ] ) with @xmath99 .
we see that the atoms tend towards one side of the ring due to the dipole instability [ fig .
[ f : mitl1 ] ( b ) ] , after which collapse occurs [ fig .
[ f : mitl1 ] ( c ) ] .
since the dipole mode has angular momenta @xmath101 and @xmath102 , the phase defect that exists from the outset slightly deviates from the center associated with the growth of the @xmath103 component , and another phase defect enters due to the @xmath104 component as seen in the inset of fig .  [
f : mitl1 ] ( c ) . for @xmath57 ,
the critical value at which the dipole mode becomes unstable is @xmath105 , which exceeds @xmath106 , and therefore the dipole instability is always followed by collapse .
we note that the dipole instability in fig .
[ f : mitl1 ] is distinct from the dissipative vortex spiral - out phenomenon that arises from the dipole mode with negative eigenvalue  @xcite . in the former ,
the quasiparticles that change angular momentum by @xmath107 are excited in pairs , and hence the angular momentum is conserved .
the energy is also conserved , since the real parts of the excitation energies have the same magnitude with opposite signs . in the latter ,
the atoms are transferred only into the non - vortex ground state decreasing the angular momentum and dissipating the energy .
a doubly - quantized vortex state @xmath108 of @xmath109 bec has recently been realized by the mit group  @xcite using a topological phase - imprinting method  @xcite , where the interaction is repulsive .
the doubly - quantized vortex with repulsive interactions is predicted to be dynamically unstable against disintegration into two singly - quantized vortices  @xcite . in the case of attractive interactions
, we will show that not only the vortex disintegration but also various new phenomena , such as a three - fold split - merge process , occur due to dynamical instabilities .
the imaginary parts of the bogoliubov eigenvalues for the doubly - quantized vortex states are shown in fig .
[ f : l2 ] .     for ( a ) @xmath54 , ( b ) @xmath110 , and ( c ) @xmath57 .
the dotted lines show the critical value @xmath58 below which the collapse of the vortex state is caused by the monopole instability .
, width=317 ]    we note that unlike the case of the singly - quantized vortex state , the quadrupole mode ( @xmath59 ) is always dynamically unstable for @xmath111 , indicating that no dynamically stable state exists for the case of attractive interactions
. therefore , an adiabatic decrease in @xmath53 from positive to negative always gives rise to the quadrupole instability .
figure  [ f : sphl2 ] shows time evolutions of the integrated densities @xmath71 ( main panels ) and @xmath112 ( lower - left insets ) and the phase profiles ( lower - right insets ) for the doubly - quantized vortex state in a spherical trap , where the initial state is a stationary state of the gp equation with @xmath97 plus the perturbation given in eq .
( [ perturb ] ) with @xmath67 .
( main panels ) , @xmath113 ( lower - left insets ) , and phase profiles ( lower - right insets ) for @xmath54 .
the initial state is the stationary doubly - quantized vortex state for @xmath114 disturbed by perturbation ( [ perturb ] ) with @xmath67 , and the interaction is switched to @xmath97 at @xmath115 .
the sizes of the images are @xmath75 for @xmath71 and the phase profiles , and @xmath116 for @xmath112 in units of @xmath117 .
, width=317 ]    for this strength of interaction , only the quadrupole mode is dynamically unstable [ cf .
[ f : l2 ] ( a ) ] .
the vortex core disintegrates , as can be seen from the density and phase profiles .
the angular momenta of the quadrupole excitations are @xmath118 and @xmath102 .
vortex disintegration is attributed to the growth of the @xmath103 component , and the four phase defects in fig .
[ f : sphl2 ] ( d ) are due to the @xmath119 component .
the density distribution exhibits various changes such as high - density regions occurring in the opposite edge of the vortices [ fig .
[ f : sphl2 ] ( c ) ] , then between the two vortices [ fig .
[ f : sphl2 ] ( d ) ] , and then the density distribution becomes similar to the split vortex as in the case of a singly - quantized vortex [ fig .
[ f : sphl2 ] ( e ) ] . finally , in fig .
[ f : sphl2 ] ( f ) the central density exceeds critical density and collapses .
another important distinction of fig .
[ f : l2 ] from fig .
[ f : l1 ] is that there is a range of @xmath53 over which the imaginary part of the eigenvalue of the octupole mode ( @xmath120 ) becomes largest .
this implies that if @xmath53 is suddenly changed into that regime , the octupolar instability dominates until the quadrupole mode grows .
figure  [ f : panl2 ] shows a typical example of such situations with @xmath55 , where the initial state is the noninteracting stationary state with vorticity of two subject to perturbation  ( [ perturb ] ) with @xmath121 , and the interaction parameter is changed from @xmath122 to @xmath123 at @xmath115 .
( main panels ) , @xmath96 ( lower - left insets ) , and phase profiles ( lower - right insets ) for the pancake - shaped trap with @xmath55 .
the initial state is the stationary doubly - quantized vortex state for @xmath114 disturbed by perturbation ( [ perturb ] ) with @xmath121 , and the interaction is switched to @xmath124 at @xmath115 .
( i ) the integrated density @xmath71 at @xmath125 , where the interaction is switched to @xmath73 and the trap is turned off at @xmath126 .
the sensitivity of the imaging is 12 times higher than that in ( a)-(h ) .
the sizes of the images are @xmath75 in @xmath71 and the phase profiles , and @xmath76 in @xmath77 in ( a)-(h ) , and @xmath78 in ( i ) in units of @xmath79 .
, width=317 ]    as the octupole mode grows , the vortex splits into three clusters which revolve around the center of the trap [ fig .  [
f : panl2 ] ( c ) ] .
interestingly , the three clusters merge to recover the original vortex [ fig .  [
f : panl2 ] ( d ) ] as in the two - cluster case shown in fig .
[ f : panl1 ] .
this is because the effective strength of interaction per each cluster @xmath127 is smaller than the critical value for collapse of the non - vortex state @xmath128 .
subsequently , the vortex once again splits into three clusters , and the balance in atomic numbers between these is then broken by the quadrupole instability .
it follows that one cluster [ denoted by a in fig .
[ f : panl2 ] ( f ) ] grows , then another cluster [ b in fig .
[ f : panl2 ] ( g ) ] grows , and finally the remaining cluster [ c in fig .
[ f : panl2 ] ( h ) ] grows and eventually collapses .
this phenomenon is similar to the seesaw - like oscillation between two clusters due to the dipole instability  @xcite . at the initial stage of disintegration
in which the octupole mode plays a dominant role [ figs .
[ f : panl2 ] ( a)-(e ) ] , the split in the phase defects is not very prominent because the angular momenta of the @xmath129 modes are @xmath130 and @xmath131 . when the quadrupole mode becomes significant , the outward shifts of the phase defects manifest themselves [ figs .
[ f : panl2 ] ( g)-(h ) ] due to the @xmath103 component of the quadrupole mode .
figure  [ f : panl2 ] ( i ) shows the integrated density at @xmath125 , following the interaction being switched from @xmath124 to @xmath73 and the trap being turned off at @xmath126 .
we see the three - fold fringe pattern resulting from interference between the matter waves that emanate from the three clusters shown in fig .
[ f : panl2 ] ( c ) .    in an elongated cigar - shaped trap
, dynamics along the trap axis also plays an important role .
figure  [ f : mitl2 ] shows the integrated density and phase profiles for @xmath57 .
( upper panels ) , phase profiles ( insets ) , @xmath96 ( lower - left panels ) , and @xmath132 ( lower - right panels ) for the cigar - shaped trap with @xmath57 .
the initial state is the stationary doubly - quantized vortex state for @xmath114 disturbed by perturbation ( [ perturb ] ) with @xmath67 , and the interaction is switched to @xmath133 at @xmath115 .
the sizes of the images are @xmath75 for @xmath71 and the phase profiles , and @xmath134 for @xmath77 and @xmath112 in units of @xmath135 . , width=317 ]    the initial state is the noninteracting stationary state with vorticity of two plus perturbation ( [ perturb ] ) with @xmath67 , and the interaction is switched from @xmath122 to @xmath133 at @xmath115 .
the atoms first gather along the symmetry axis ( i.e. , the @xmath136 axis ) towards @xmath137 [ fig .  [
f : mitl2](b ) ] , then the vortex begins to disintegrate since the density is high enough .
subsequently the atoms rebound along the @xmath136 direction [ fig .
[ f : mitl2](c ) ] .
we see that the disintegrated vortex lines in fig .
[ f : mitl2 ] ( c ) are straight , unlike the repulsive case where twisting of vortex lines occurs  @xcite . in fig .
[ f : mitl2 ] ( d ) , the condensate splits into two high - density regions containing the disintegrated vortices as a consequence of the rebound .    in an experiment using @xmath1 ,
the lifetime of the condensate is as short as 1 or 2 ms near the feshbach resonance  @xcite , which indicates that the lifetime is @xmath138 ( in our dimensionless units of time ) for the optical dipole trap used in the mit group ( @xmath139 hz , @xmath140 hz , and @xmath56 )  @xcite . for dynamical instabilities to occur within this short lifetime
, the strength of attractive interaction @xmath61 must be much larger than those in the above cases .
figure  [ f : mitl2exp ] shows the density and phase profiles just before the collapse , where the initial state is the noninteracting stationary state with vorticity of two plus perturbation ( [ perturb ] ) with @xmath141 , @xmath30 , @xmath31 , and @xmath142 , and the interaction is switched from @xmath122 to @xmath143 at @xmath115 .    , phase profiles , and @xmath132 ( from top to bottom ) for the cigar - shaped trap with @xmath57 .
the initial state is the stationary doubly - quantized vortex state for @xmath114 disturbed by perturbation ( [ perturb ] ) with ( a ) @xmath99 ,
( b ) @xmath67 , ( c ) @xmath121 , and ( d ) @xmath144 .
the interaction is switched from @xmath122 to @xmath145 at @xmath115 .
the sizes of the images are @xmath75 for @xmath71 and the phase profiles , and @xmath146 for @xmath112 in units of @xmath117 .
, width=317 ]    we find that the collapsing dynamics depend on the initial perturbation @xmath51 , since the system has several multipole instabilities for @xmath145 .
we note that the vortex split occurs only around @xmath137 , where the density is high enough .
thus , the vortex - split phenomena can be observed for large values of @xmath61 , even when the lifetime of the condensate is very short .
to understand the split - merge phenomena shown in the preceding sections , we restrict ourselves to 2d space in this and next sections .
we employ a gaussian variational wave function for the condensate as @xmath147 where @xmath26 is the vorticity and @xmath148 is a variational parameter characterizing the size of the condensate . substituting eq .
( [ gaussian0 ] ) into the gp energy functional @xmath149 = \int d{\bf r } \left [ \frac{1}{2 } \psi^ * \left ( -\nabla^2 +   \omega^2 r^2 \right ) \psi + \frac{g}{2 } |\psi|^4 \right]\ ] ] and minimizing it with respect to @xmath148 gives @xmath150^{\frac{1}{4}}.\ ] ] for @xmath151 for @xmath152 and @xmath153 for @xmath154 agree reasonably well with the numerically obtained critical values @xmath155 and @xmath156 , respectively .    to study the time evolution from eq .
( [ gaussian0 ] ) , we use a variational wave function given by @xmath157 where @xmath158 and @xmath159 is given by eq .
( [ gaussian0 ] ) with @xmath26 replaced by @xmath4 . in eq .
( [ varpsi ] ) we use the same width @xmath148 as eq .
( [ d ] ) for all @xmath159 and neglect the radial degrees of freedom to simplify the following calculation and to obtain the results analytically .
substituting eq .
( [ varpsi ] ) into the action @xmath160 = \int dt e[\psi ] - i \int d{\bf r } dt \psi^ *   \frac{\partial}{\partial t } \psi\ ] ] yields @xmath161 where @xmath162 , and @xmath163 , which vanishes for @xmath164 . from @xmath165 , the variational equation of motion for @xmath166
becomes @xmath167 we assume @xmath168 and @xmath169 , and neglect terms of the order of @xmath170 .
equation  ( [ eomc ] ) for @xmath171 reduces to @xmath172 and hence we obtain @xmath173 . using eq .
( [ c1 ] ) and defining @xmath174 , we obtain a closed set of equations    @xmath175    from which the eigenfrequencies can be found .
the condition under which the eigenfrequencies are real is given by @xmath176 ^ 2 \nonumber \\ & & - 4g^2 \left ( g^{m , 2q - m}_{q , q } \right)^2 \geq 0.\end{aligned}\ ] ]    first we consider the singly - quantized vortex ( @xmath152 ) .
stability of the quadrupole mode is determined by @xmath177 which becomes negative for @xmath178 .
this value is in reasonable agreement with the numerically obtained critical value for the quadrupole mode in 2d , @xmath179 .
however , the discriminant @xmath180 is always positive , which contradicts the numerical result that dipole instability does occur for @xmath181 .
also , for @xmath182 or @xmath119 , the numerical calculation shows that dynamical instability occurs for @xmath183 , while @xmath184 is always positive for @xmath185 [ @xmath53 must be larger than @xmath186 as seen from eq .
( [ d ] ) ] . in the case of the doubly - quantized vortex ( @xmath154 )
, @xmath187 is always negative , which is consistent with the fact that the quadrupole mode is always dynamically unstable for @xmath188 .
however , @xmath189 fails to reproduce the behavior for @xmath190 , where stable and unstable regions appear alternatively  @xcite .
the discriminant @xmath191 becomes negative for @xmath192 , which is in reasonable agreement with the numerical result that the octupole instability sets in for @xmath193 .
however , @xmath194 contradicts the numerical result that dipole instability does occur for @xmath195 .
thus , the gaussian ansatz ( [ varpsi ] ) is viable for some cases , whereas more sophisticated trial functions are needed to completely describe the dynamical instabilities .
the variational results can be improved by using the different variational width @xmath148 for each @xmath159 in eq .
( [ varpsi ] ) , which may be determined by , e.g. , the method used in ref .  @xcite . however , near and above the critical point , this variational method fails , since the norm of the bogoliubov function vanishes when a complex eigenvalue emerges .
equation  ( [ d13 ] ) is expanded near the critical point as @xmath196 where @xmath197 . just above the critical point , therefore , the imaginary part of the complex eigenvalue is proportional to @xmath198  @xcite .
all the imaginary parts shown in figs .  [ f : l1 ] and [ f : l2 ] ( except @xmath59 in fig .
[ f : l2 ] ) have this @xmath53 dependence near the critical points [ so do the @xmath63 curves in figs .  [ f : l1 ] ( c ) and [ f : l2 ] ( c ) in the immediate vicinity of the critical points ] .
the @xmath59 lines for @xmath154 in fig .
[ f : l2 ] are proportional to @xmath61 near @xmath122 , which agrees with @xmath199 .    the split - merge cycles shown in figs .
[ f : panl1 ] and [ f : panl2 ] are recurrence phenomena arising from nonlinearity . to understand the behavior in fig .
[ f : panl1 ] qualitatively , we reduce the number of modes in eq .
( [ eomc ] ) to three , i.e. , @xmath200 , @xmath31 , and @xmath131 .
equation  ( [ eomc ] ) then reduces to simultaneous nonlinear differential equations with three variables , which can be solved numerically .
figure  [ f:3mode ] shows the time evolution of @xmath201 , @xmath202 , and @xmath203 for @xmath204 , where the initial condition is @xmath205 and @xmath206 .
( solid curve ) , @xmath207 ( dashed curve ) , and @xmath208 ( dotted curve ) in the three - mode approximation of eq .
( [ eomc ] ) for @xmath204 .
the initial condition is @xmath209 and @xmath210 .
the insets show the density and phase profiles at @xmath115 and @xmath211 , where the sizes of the images are @xmath75 in units of @xmath212 .
, width=317 ]    the vortex can be seen to split at @xmath213 , but subsequently recovers the original vortex during @xmath214 , and then splits again at @xmath215 .
thus , the split - merge phenomenon can be qualitatively reproduced by a simple three - mode approximation .
the split - merge dynamics shown in fig .  [ f : panl1 ] implies the existence of a low - lying state in which two clusters revolve around each other without changing their shapes . to find such a low - lying stationary state ,
we numerically seek a local minimum of the gp energy functional with the angular momentum held fixed at that of the initial vortex state @xmath216 , where @xmath217 .
the gp energy functional with lagrange multipliers is given by @xmath218 = e[\psi ] + \int d{\bf r } \left ( -\mu |\psi|^2 + \omega \psi^ *   \hat l_z \psi \right),\ ] ] where @xmath219 $ ] is given in eq .
( [ e ] ) , and the chemical potential @xmath22 and the angular frequency @xmath220 of rotation of the system are introduced as the lagrange multipliers in order to fix the normalization and the angular momentum .    using the newton - raphson method  @xcite starting from appropriate initial states ,
we obtain the stationary states in the frame of reference rotating at frequency @xmath220 .
figure  [ f : lowly ] shows various properties of the axisymmetric vortex state and the low - lying state .
( solid curves ) and chemical potential @xmath22 ( dashed curves ) of the axisymmetric singly - quantized vortex state ( upper branch ) and low - lying soliton state ( lower branch for @xmath221 ) .
gray - scale images show the density profile of the axisymmetric vortex state for @xmath222 and those of the low - lying states for @xmath223 and @xmath224 .
( b ) the bogoliubov excitation energies of the quadrupole modes above the axisymmetric vortex state ( @xmath225 ) and those of the corresponding modes above the low - lying soliton state ( @xmath221 ) .
( c ) the number of virtually excited atoms in the bogoliubov modes .
the inset is the logarithmic plot of the excitation energies and the number of excited atoms ( both of them refer to the common scale ) as functions of @xmath226 .
the vertical dotted lines in ( a)-(c ) show @xmath227 . , width=264 ]    when @xmath61 is increased , the stationary state bifurcates at the critical point @xmath228 of the quadrupole instability , at which the low - lying state emerges as represented by the lower branches of @xmath229 and @xmath22 in fig .
[ f : lowly ] ( a ) .
density profiles of the low - lying states for @xmath223 and @xmath124 are shown as gray - scale images in fig .
[ f : lowly ] ( a ) .
it is clear from them that the split - merge cycle is due to oscillations between the axisymmetric vortex state ( upper branch ) and the low - lying soliton state ( lower branch ) .
it is interesting to note that the energies of the axisymmetric and low - lying states are smoothly connected at the critical point , while the chemical potential has a kink , which implies that the nature of the phase transition between the vortex state and the soliton state is of second order .
the excitation energies are plotted in fig .
[ f : lowly ] ( b ) , where @xmath230 denotes the energy of the lowest excitation with angular momentum @xmath231 above the axisymmetric vortex state . as @xmath53 approaches @xmath228 in the axisymmetric regime , the excitation energy @xmath232 vanishes , despite the fact that both @xmath233 and @xmath233 are nonzero at @xmath234 .
this implies that quasiparticles with @xmath235 can be excited in pairs without changing energy and angular momentum .
hence , the two modes @xmath235 may be strongly coupled via the pair excitations . if we follow the axisymmetric branch [ the upper branch in fig .
[ f : lowly ] ( a ) ] , the eigenenergies of the two modes become complex .
if we follow the low - lying soliton branch , one of the two modes becomes the goldstone mode associated with the axisymmetry breaking , and the energy of the other mode @xmath236 increases with decreasing @xmath53 .
figure  [ f : lowly ] ( c ) shows the number of virtually excited atoms @xmath237 in the quadrupole modes @xmath238 for @xmath239 and that in the corresponding mode @xmath240 for @xmath241 . at @xmath234 the number of excited atoms diverges , indicating that the bogoliubov approximation breaks down .
the inset of fig .
[ f : lowly ] ( c ) shows the logarithmic plot of the excitation energies and the numbers of excited atoms .
close to the critical point @xmath242 , the excitation energies are found to be proportional to @xmath243 and the numbers of excited atoms to @xmath244 .
we note that these behaviors near the critical point ( the kink in @xmath22 , the softening of the excitation mode , the divergence in the number of excited atoms , and their power laws ) are very similar to those occurring in the transition between the uniform - density and soliton states in a 1d ring  @xcite . thus , in the immediate vicinity of the critical point at which the vortex split occurs , the mean - field and bogoliubov approximations break down due to large quantum fluctuations , and we must employ , e.g. , the exact diagonalization method  @xcite to study this regime
. however , the mean - field theory used in the present paper is still valid except for this narrow critical region .
we studied dynamical instabilities in singly- and doubly - quantized vortex states with attractive interactions , and found a rich variety of dynamics triggered by the dynamical instabilities .
the dynamics depend on the strength of interaction @xmath53 , the trap geometry @xmath62 , the initial perturbations , and the manner ( adiabatic or sudden ) in which @xmath53 changes .
a singly - quantized vortex in a trap with @xmath245 first exhibits the quadrupole instability as @xmath53 is adiabatically decreased .
the vortex then splits into two clusters revolving around the center of the trap . in a spherical trap the split fragments immediately collapse for any value of @xmath53 that gives rise to the split instability , whereas in a pancake - shaped trap with @xmath55
they unite to recover the original vortex and then the split - merge cycles repeat for @xmath61 slightly exceeding @xmath95 . in a cigar - shaped trap with @xmath64 ,
the dipole instability causes the atoms to concentrate in one cluster and then the system collapses .
the doubly - quantized vortex state is always dynamically unstable against disintegration of the vortex core for @xmath188 , and therefore the disintegration is unavoidable if we adiabatically decrease @xmath53 .
if @xmath53 is suddenly changed to some negative value with an appropriate initial perturbation to the condensate , the octupole instability can dominate the quadrupole instability , and the vortex splits into three clusters . in a pancake - shaped trap , these undergo split - merge cycles as in the case of a singly - quantized vortex , but eventually one of them collapses due to the quadrupole instability that always exists .
we have also demonstrated that the split phenomena can be observed using the current experimental setup of the mit group .
to understand the dynamical instability and the split - merge cycles , we performed a variational analysis .
the emergence of a complex eigenvalue in some modes was reproduced by the gaussian trial functions .
we were able to reproduce the split - merge cycles by assuming only three relevant modes with angular momenta @xmath246 and @xmath247 , suggesting that only a few modes are associated with the split - merge cycles .    in the immediate vicinity of the critical strength of interaction of the dynamical instability ,
the transition from the axisymmetric vortex state to the split state is accompanied by divergence in the number of virtually excited atoms and by softening of the excitation mode according to the bogoliubov theory , which is very similar to the 1d case discussed in ref .
the dynamical instabilities in vortex states of attractive becs therefore involve a fundamental generic issue of the decay of a many - particle quantum system .
vortex - split phenomena are attributed to the interplay between the attractive interaction and the topological phase defect , that is , the condensate is bound to break axisymmetry , since they can not gather at the phase defect .
the emergence of the dynamical instabilities and ensuing dynamics may be regarded as pattern formation due to symmetry breaking .
more complicated and interesting pattern might be formed in the vortex lattice , which will be published elsewhere .
dynamics after the collapse were not studied because 3d simulation of the collapse and explosion phenomena exceeds our present computational ability .
the burst atoms ejected from each cluster interfere with each other and exhibit some anisotropic pattern if the burst atomic cloud is coherent .
therefore experiments on vortex split and fragmented collapse will provide a test of whether or not the burst atoms are coherent .
we thank j. k. chin for discussions .
this work was supported by the special coordination funds for promoting science and technology and a grant - in - aid for scientific research ( grant no . 15340129 ) by the ministry of education , science , sports , and culture of japan , and by the yamada science foundation .