matter - wave solitons , that can be created in bose - einstein condensates ( becs ) , are a research subject of great interest , both as nonlinear collective excitations in macroscopic quantum matter , and due to the potential they offer in applications to high - precision interferometry , atomic - wave soliton lasers , quantum information processing , and other emerging technologies .
dark and bright matter - wave solitons were experimentally created in nearly one - dimensional ( 1d ) trap configurations @xcite , filled with @xmath1rb and @xmath2li condensate , respectively . in these settings ,
strong radial confinement freezes transverse dynamics of the condensate , keeping atoms ( with mass @xmath3 ) in the ground state of the corresponding 2d harmonic - oscillator potential , @xmath4 , while a weak axial parabolic trap , @xmath5 , allows quasi - free motion of the soliton in the axial direction , @xmath6 .
new recent experiments with @xmath7rb @xcite and related numerical simulations @xcite have showed the existence of stable solitary waves in a weakly elongated ( @xmath8 ) trap . these essentially three - dimensional ( 3d ) bright solitons exhibit
rich and complex behavior , not present in their nearly-1d counterparts .
notable among the new features , are inelastic collisions between solitons , which are strongly sensitive to phase relations between them and their relative speed .    adding a periodic optical - lattice ( ol ) potential in the axial direction of the quasi-1d trap
makes it possible to create bright matter - wave solitons of the gap type in repulsive condensates @xcite , which is a particular manifestation of very rich physics of bec trapped in ols , as reviewed in refs .
@xcite .
recently , 2d and 3d localized states in the gross - pitaevskii equation ( gpe , which is the basic mean - field bec model ) , that may be sustained and stabilized by cellular potentials , i.e. , multidimensional ols , have attracted a great deal of interest , see review @xcite .
a great challenge to the experiment is creation of multidimensional matter - wave solitons , as well as the making of spatiotemporal solitons in nonlinear optics @xcite . in particular , multi - dimensional solitons trapped in a low - dimensional ol ( i.e. , 1d lattice in the 2d space @xcite , and 2d lattice in the 3d space @xcite ) , have been predicted @xcite .
as these solitons keep their mobility in the free direction , the latter settings can be used to test head - on and tangential collisions between solitons .
note that the quasi-1d lattice potential can not stabilize 3d solitons @xcite ( this becomes possible if the quasi-1d lattice is combined with the periodic time modulation of the nonlinearity provided by the feshbach - resonance - management technique @xcite ; a general account of the technique of periodic management for solitons can be found in book @xcite ) .
the low - dimensional ol can also be implemented as a radial ( axially symmetric ) lattice in both 2d @xcite-@xcite and 3d @xcite settings . in both cases , stable solitons in the self - attractive medium
were predicted , in the form of a spot trapped either at the central potential well , or ( in the 2d case ) in a radial potential trough ( additionally , the so - called azimuthons were predicted as azimuthally periodic deformations of vortices in the uniform 2d medium @xcite ; however , they are unstable in the case of the cubic nonlinearity ) . in the latter case , the spot - shaped soliton may run at a constant angular velocity in the trough ( thus being a _ rotary soliton _
the circular motion may be more convenient for the experimental study of mobility and collisions of matter - wave solitons ( and soliton trains ) than the already realized soliton - supporting settings in the cigar - shaped traps @xcite , as in the circular geometry the motion is not affected by the weak longitudinal trap .
this setting is also convenient for studies of patterns in self - repulsive becs . in particular ,
dark solitons were , thus far , created in cigar - shaped traps @xcite , where the background density is essentially non - uniform , as it vanishes at edges of the trap . in a circular trough ,
pairs of dark solitons can be created without this complication .
in fact , both axisymmetric ( azimuthally uniform ) vortices @xcite and dipole and quadrupole states , that may be regarded as stable complexes of two or four dark solitons @xcite , were predicted in the two - dimensional gpe with the repulsive nonlinearity and bessel - lattice radial potential .
the stable dipole and quadrupole patterns may also rotate at a constant angular velocity @xcite .
radial potentials in the form of @xmath9 or @xmath10 , adopted in refs . @xcite and @xcite , respectively ( @xmath11 and @xmath12 are positive constants , and @xmath13 and @xmath14 are the bessel functions ) , can support patterns of most types considered in those works in the linear limit of the model . indeed , axisymmetric states trapped in the central potential well , as well as azimuthally uniform vortices and dipolar and quadrupolar modes trapped in circular potential troughs , may be found as 2d eigenstates of the quantum - mechanical linear schrdinger equation in such radial potentials .
accordingly , nonlinear states with the same symmetries found in refs .
@xcite may be considered as continuations of their linear counterparts .
an exception is the spot - shaped soliton in a circular trough ( the one that can perform rotary motion along the trough , as mentioned above ) . on the contrary to that , in this paper we report results for the _ radially periodic _ ol potential , @xmath15 , which , obviously , can not support any radially localized state in the linear limit .
therefore , all localized patterns reported in this paper are essentially nonlinear objects , and they may be different from their counterparts found in the bessel lattices , as concerns their existence and stability conditions . the solitons that we report in this work may be categorized as solitons of the ordinary type in the model with self - attraction , and _ radial gap solitons _
( a new species of solitary waves ) , in the case of self - repulsion .
note that solitons of the gap type can not exist in the above - mentioned bessel - lattice potentials .
as for the attractive model , a somewhat similar problem was considered in refs . @xcite and
@xcite , where it was demonstrated that the self - attractive nonlinearity in the 1d and 2d geometry ( the latter setting is axisymmetric ) can turn a quasi - bound state in a potential well dug on top of a potential hill into a true bound state ( including 2d vortices @xcite ) . in the 3d setting , the situation is more complicated , but a similar effect is also possible @xcite .
it is relevant to mention that , besides the attractive or repulsive cubic nonlinearity , the radial - lattice potential can also be combined with the self - focusing saturable nonlinearity , characteristic of photorefractive media . it was predicted @xcite and very recently experimentally demonstrated @xcite that stable solitons are possible in this case too ( in fact , ref .
@xcite presents the first ever experimental realization of solitons in any radial lattice ) .
another very recent experimental result is the observation of light localization , in the form of spot - like and ring - like patterns , in a photorefractive medium equipped with a third - order bessel lattice subjected to azimuthal modulation , which corresponds to a combination of saturable nonlinearity and an effective potential @xmath16 , where @xmath17 is the angular coordinate @xcite .
a different experimental setup that also leads to a ring - shaped trap for the bec has been recently reported , based on a specially designed configuration of the external magnetic field that creates a guide for matter waves in the form of a torus @xcite . in parallel , nearly-1d solitons in toroidal traps loaded with self - attractive bec were studied theoretically @xcite ( in the experiment of ref .
@xcite , the condensate was self - repulsive ) .
the objective of this paper is to investigate various types of matter - wave solitons supported by the periodic radial lattice potential .
it should be mentioned that , besides the bec context , the model and results obtained in it may also be relevant to photonic - crystal fibers with a concentric ( rather than usual hexagonal ) structure , similar to that in fibers with the multilayer cladding @xcite . in that case
, the soliton will be a self - trapped beam propagating in the fiber .
the paper is organized as follows . in section
ii we formulate the model based on the respective two - dimensional gpe .
section iii is dealing with solitons trapped in the center ( including the radial gap solitons , in the case of the model with repulsion ) , for which results are obtained by means of the variational approximation ( va ; both static and dynamical variants of the va are presented ) , and in direct simulations .
section iv treats solitons trapped in annular channels ( potential troughs ) . by means of the va ,
we derive an effective 1d equation for the attractive model , which predicts a possibility of the existence of stable axisymmetric ring - shaped states , and the onset of their instability against azimuthal modulations at a finite density .
exact modulated ( cnoidal - wave ) solutions , which appear above the instability threshold of the uniform state , are found too , as well as exact solutions for azimuthal solitons .
direct numerical simulations corroborate the existence of stable uniform and modulated ring - shaped patterns , and reveal the radial gap solitons whose central love is trapped in the annular trough ( in the repulsive version of the model ) .
also reported are stable ring - shaped gap solitons carrying a pair of dark solitons on their crests .
section v addresses dynamics of moving two - dimensional solitons , trapped in the circular troughs of the attractive model .
it is shown that the slow solitons are stable , while fast ones are destroyed due to leakage of matter into the outer trough , under the action of the centrifugal force .
head - on collisions between in - phase solitons lead to collapse ( intrinsic blow - up ) , while @xmath0-out - of - phase solitons collide many times , and eventually merge into a single soliton ( without collapsing ) ; the same happens with in - phase solitons colliding tangentially , as they are trapped in adjacent troughs .
finally , section vi concludes the paper .
the starting point is the standard gpe in a normalized form @xcite , @xmath18 written for the single - atom wave function @xmath19 in polar coordinates @xmath20 . here ,
@xmath21 and @xmath22 correspond to the self - focusing and defocusing nonlinearity , respectively ( alias negative and positive scattering length of interatomic collisions ) , and , as said above , @xmath23 is the radially periodic potential with strength @xmath11 and wavenumber @xmath24 ( by dint of obvious rescaling , we set @xmath25 ) . in terms of ols ,
this potential can be created by a cylindrical beam whose amplitude is modulated as @xmath26 ( see further discussion below ) . note that the potential with @xmath27 gives rise to a local minimum at @xmath28 , see an example in fig .
[ fig1 ] .    together with eq .
( [ gpe2d ] ) , we will use its variational representation .
obviously , the equation can be derived from the lagrangian @xmath29 ,   \label{lagr}\end{aligned}\ ] ] where @xmath30 stands for the complex conjugate .
.,width=302,height=151 ]    stationary axially symmetric states with chemical potential @xmath31 are looked for as solutions to eq .
( [ gpe2d ] ) in the form of @xmath32 , with a real function @xmath33 obeying the equation@xmath34 for given potential @xmath35 , solutions to eq .
( [ u ] ) form families parameterized by @xmath31 or , alternatively , by the norm,@xmath36    solutions to stationary equation ( [ u ] ) with @xmath21 ( self - attraction ) were obtained by means of a recently developed spectral - renormalization method for finding self - localized states of nonlinear - schrdinger ( nls ) type equations @xcite .
independently , the same stationary solutions were also generated by dint of the well - known method which uses the integration of gpe ( [ gpe2d ] ) in imaginary time @xcite .
stability of the solitary waves was then tested by adding small perturbations to them and integrating the gpe in real time .
a stable perturbed soliton would shed off some radiation , which was absorbed at boundaries of the integration domain , and relax into a slightly different stationary form , corresponding to the smaller norm .    in the case of repulsion ( @xmath37 ) ,
self - localized states were obtained by solving eq .
( [ gpe2d ] ) in real time , with an absorber placed at the boundary of the integration domain . in that case , a localized waveform ( e.g. , gaussian ) with a suitable norm was chosen as the initial condition .
in the course of the evolution , excess norm is radiated away with linear waves , which are absorbed at the boundary , and a self - localized state emerges ( see an example in fig .
[ fig4 ] below ) .
it is obvious that solitons found this way are stable .
the above - mentioned radial lattices of the bessel type can be naturally generated by a nondiffracting linear optical beam with the cylindrical symmetry ( bessel beams themselves are usually generated by means of axicoms ) @xcite .
on the other hand , the radial potential of the form @xmath38 can be induced , as said above , by a cylindrical beam whose amplitude is modulated as @xmath26 . the modulation can be provided by passing the cylindrical laser beam through a properly shaped plate @xcite .
unlike the bessel beam , the one with the @xmath26 transverse modulation will suffer conical diffraction , but it is not a problem for bec  experiments , as a tight optical trap created in the transverse direction may readily contain a 2d pancake - shaped configuration of the condensate , of the thickness @xmath39 @xmath31 m @xcite .
the diffraction of the paraxial beam , with waist diameter @xmath40 @xmath31 m , relevant to the experiment , on such a short propagation distance is completely negligible , as the respective diffraction length ( alias rayleigh range ) is estimated to be @xmath41 mm , where @xmath42 @xmath31 m is the beam s carrier wavelength .
it is known that a radial potential structure can easily trap a soliton in the central potential well @xcite .
we start the consideration of the present model with solutions of this type .
typical examples of the solitonic shape in the model with self - attraction ( @xmath21 ) are displayed in fig .
these solutions were obtained by means of both the above - mentioned spectral - renormalization method @xcite and integration in imaginary time @xcite .    , trapped ( in the model with self - attraction ) at the center of the radial structure , @xmath43 , with @xmath44 ( in this case , the norm is @xmath45 ; solid line ) , and @xmath46 ( @xmath47 , dashed line).,width=302,height=226 ]    in the repulsive model ( @xmath37 ) , numerical solutions were obtained , as said above , by solving eq .
( [ gpe2d ] ) in real time , in the presence of the boundary absorber .
an example of the relaxation of an initial gaussian pulse into a soliton is given in fig .
full views of the typical stable solitons trapped at the center of the radial structure in the model with self - attraction and self - repulsion are displayed in fig .
[ fig4 ] .
( in the presence of absorbers at the boundary of the integration domain ) from an initial gaussian with norm @xmath48 .
the norm of the eventually established state is @xmath49 , i.e. , @xmath50 of the initial value.,width=302,height=151 ]     ( with @xmath45 , @xmath51 , and @xmath52 ) .
right panel : a localized state in the repulsive bec , with @xmath49 , @xmath53 , amplitude @xmath54 , and @xmath55 , which was generated as shown in fig .
[ fig3].,title="fig:",width=151,height=113 ]   ( with @xmath45 , @xmath51 , and @xmath52 ) .
right panel : a localized state in the repulsive bec , with @xmath49 , @xmath53 , amplitude @xmath54 , and @xmath55 , which was generated as shown in fig .
[ fig3].,title="fig:",width=151,height=113 ]    as concerns the solitons in the repulsive model ( @xmath37 ) , a typical example of which is displayed in the right panel of fig .
[ fig4 ] , they may be clearly classified as solitons of the gap type , since they exist solely due to the interplay of the lattice potential and self - defocusing nonlinearity . the ordinary setting that gives rise to gap solitons is quite similar , with the ol in the 1d , 2d or 3d cartesian coordinates @xcite . slowly decaying fringes , attached to the central core of the soliton and evident in the right panel of fig .
[ fig4 ] ( and absent in the ordinary soliton displayed in the left panel ) , is a characteristic feature of gap solitons .
thus , this self - trapped state may be naturally called a _ radial gap soliton_.      to describe a stationary soliton trapped in the central potential well in an analytical approximation , we adopt a gaussian ansatz , @xmath56 comparison with fig .
[ fig2 ] shows that , in the model with attraction ( @xmath21 ) , this ansatz is quite adequate for solitons with a larger norm , trapped in a relatively weak lattice , but it may be inaccurate as an approximation for solitons with smaller @xmath57 trapped in a stronger lattice .
the application of the variational approximation ( va ) to stationary equation ( [ u ] ) , with the help of the accordingly modified lagrangian ( [ lagr ] ) , is straightforward ( cf .
@xcite , where the same gaussian ansatz was used to predict 2d solitons trapped in quasi-1d and square ols ) .
thus we arrive at a set of va - generated equations that relate the soliton s size @xmath58 to its norm and chemical potential : @xmath59 -\frac{1}{a^{2 } } ,   \label{va1}\ ] ] where @xmath60 , @xmath61 , and @xmath62is the standard error function .    the analysis of dependence @xmath63 following from eq .
( [ va1 ] ) with @xmath21 ( self - focusing ) predicts the existence of a family of solitons trapped at the center of the radial lattice , which are stable according to the vakhitov - kolokolov ( vk ) stability criterion @xcite , @xmath64 , as shown in fig .
although the vk criterion is only necessary for the stability , as it ignores a possibility of oscillatory instability with complex eigenvalues , direct simulations of eq .
( [ gpe2d ] ) confirm the stability of these solutions ( see below ) .     following from variational equations ( [ va1 ] ) , for the self - attraction case , @xmath65 .
according to the vk criterion , @xmath66 , stable solutions belong to upper branches of the curves.,width=302,height=226 ]    additional evidence for the existence of self - trapped localized states of bec in the present model can be provided by a more general , time - dependent , version of the va . to this end
, we define a generalized ansatz [ cf .
( [ ansatz - zentrum ] ) ] @xmath67 ,   \label{ansatz}\ ] ] where @xmath68is a real radial chirp . applying the standard va procedure @xcite ,
one arrives at the following evolution equation for the width of the localized state , @xmath69 ,   \label{att}\ ] ] where an effective nonlinearity strength is @xmath70 and the amplitude is given by @xmath71 equation ( [ att ] ) describes motion of a unit - mass particle with the coordinate @xmath72 in an effective potential @xmath73 in fig . [ fig6 ] the effective potential ( [ pot ] ) is depicted for different values of @xmath74 , for the attractive and repulsive versions of the model .    ) for different values of @xmath75 , at fixed strength @xmath76 of the radial optical lattice.,title="fig:",width=151,height=151 ] ) for different values of @xmath75 , at fixed strength @xmath76 of the radial optical lattice.,title="fig:",width=151,height=151 ]    solitons in the attractive bec ( @xmath21 ) exist when potential ( [ pot ] ) possesses a local minimum . in particular , for @xmath77 ( see fig .
[ fig6 ] ) , the potential minimum exists provided that @xmath78 . for the norm smaller than that corresponding to @xmath79
[ see eq . ( [ chi ] ) ] , the self - attraction is too weak to form a soliton , and the matter - wave pulse spreads out through barriers of the radial potential . on the other hand , in the case of attraction , @xmath74 is also limited from above , by the collapse occurring at a critical value , @xmath80 . in the absence of the lattice ( @xmath81 ) , a numerically exact critical norm is @xmath82 ( it is equal to the norm of the _ townes soliton _ ) @xcite . for any @xmath11 ,
the va predicts @xmath83 , as the value of coefficient @xmath74 at which the first term in effective potential ( [ pot ] ) vanishes .
according to eq .
( [ chi ] ) , this gives @xmath84 @xcite .
thus , the va predicts a finite existence region for stable self - trapped localized states in the attractive bec ; for example , it is @xmath85 for @xmath77 .
while the upper edge of this region is fixed at @xmath86 , numerical analysis of expression ( [ pot ] ) demonstrates that the lower critical norm necessary for the existence of the self - trapped state almost linearly decreases with the increase of the ol strength @xmath11 .
numerically exact solution of the @xmath87-independent version of full gpe ( [ gpe2d ] ) demonstrates that , although the collapse sets in at a critical value of the norm which is slightly higher than @xmath84 , similar to the case of the 2d soliton stabilized by the quasi-1d or square ol @xcite , the dependence of @xmath88 on @xmath11 is weak .
the physical mechanism behind the rise of the collapse threshold in a trapping potential is that the confinement allows the soliton to increase the internal  quantum pressure " , generated by the dispersion term in gpe ( [ gpe2d ] ) , which counteracts the nonlinear self - focusing of the wave packet .
predictions of the va for dynamical regimes can be compared to direct simulations in terms of the time dependence of the amplitude , @xmath89 , in a soliton which was excited by a sudden perturbation [ the variational prediction for @xmath89 is taken as per eq .
( [ ampl ] ) ] . to display an example , we take the stationary solution from fig .
[ fig2 ] , corresponding to @xmath77 and @xmath47 , and suddenly increase the lattice strength to @xmath90 .
comparison of the variational results with simulations of eq .
( [ gpe2d ] ) is presented in fig .
as observed , agreement for absolute values of the amplitude is , at best , qualitative , while the oscillation frequencies coincide much better .    )
(  gpe " ) with @xmath65 , in terms of the soliton s amplitude as a function of time .
intrinsic vibrations of the soliton were triggered by sudden increase of the lattice strength from @xmath46 to @xmath91 , for the self - trapped state with @xmath47 in fig .
[ fig2].,width=226,height=151 ]    in the case of the repulsive bec ( @xmath37 ) , the behavior is opposite , as seen in the right panel of fig .
[ fig6 ] . in this case , a larger norm for given ol strength @xmath11 causes spreading of the localized state , because the self - defocusing nonlinearity overcomes the confining effect of the trapping potential
. simulations of eq .
( [ gpe2d ] ) with @xmath37 show that the wave packet sheds off excessive norm in the form of linear waves , which are absorbed at domain boundaries , and relaxes into a stable shape after the norm becomes adequate .
for example , at @xmath92 stable localized states in the repulsive bec exist at @xmath93 @xmath94 [ i.e. , the norm is limited to values @xmath57 @xmath95 , see eq .
( [ chi ] ) ] .
the most essential property of radially periodic potentials is that localized states can self - trap not only in the central potential well , but also in remote circular troughs @xcite .
if the trough s curvature is not significant , these localized states resemble 2d solitons in the quasi-1d ol reported in ref .
@xcite . here , we aim to derive an asymptotic 1d equation for patterns trapped in a circular trough , and obtain its relevant solutions . to this end , we adopt the following ansatz for the solution , @xmath96 where @xmath97 is the coordinate running along the circumference of radius @xmath98 , @xmath99 is a slowly varying complex amplitude , and the transverse variable is @xmath100 , with @xmath101 being a radial potential minimum in a vicinity of which the ring - shaped pattern is trapped . assuming @xmath3 to be large enough , we consider the annular pattern placed in a remote circular trough , which may be treated as a quasi - rectilinear potential trap , with negligible curvature . obviously , the latter condition amounts to @xmath102 [ recall @xmath103 is the radial size in ansatz ( [ ansatz1 ] ) ] .    neglecting the curvature in eq .
( [ gpe2d ] ) , we rewrite it as a 2d equation with @xmath6 and @xmath104 treated as local cartesian coordinates , @xmath105 the next step is to reduce eq .
( [ no - curv ] ) to an effectively 1d equation along azimuthal direction @xmath104 . to this end
, we note that eq .
( [ no - curv ] ) corresponds to the lagrangian [ cf .
( [ lagr ] ) ] @xmath106.\end{aligned}\ ] ] in this lagrangian , we substitute ansatz ( [ ansatz1 ] ) and integrate the resulting expression in transverse direction @xmath6 . after a straightforward calculation , this yields an effective ( averaged ) lagrangian , @xmath107 , \label{leff}\end{aligned}\ ] ] where @xmath108 .
the application of the standard variational procedure to averaged lagrangian ( [ leff ] ) gives rise to a cumbersome system of coupled equations for the complex amplitude @xmath99 and real width @xmath109 .
the equations strongly simply if we assume that @xmath110 and @xmath103 are slowly varying functions of @xmath104 , which , in particular , corresponds to long - wave perturbations that account for the onset of the modulational instability ( mi ) , see below .
in particular , in the attractive model ( @xmath21 ) one may first adopt the lowest - order approximation , completely neglecting the @xmath104-dependence and assumes time - independent @xmath111 and @xmath103 , which gives rise to a result previously known from the application of the va to the nls equation in 1d with potential @xmath112 @xcite@xmath113 :    @xmath114 } .
\label{a^2}\ ] ]    this relation is cumbersome too , but straightforward analysis demonstrates that it takes a very simple form , @xmath115 , when @xmath11 is both small and large [ in the latter case , scaling @xmath116 is assumed , to lend the last three term in the lagrangian density in eq .
( [ leff ] ) the same order of magnitude ] .    to obtain a closed - form lagrangian for the amplitude field @xmath99 in the attractive model
, one should express @xmath103 in eq .
( [ leff ] ) in terms of @xmath117 , using relation ( [ a^2 ] ) .
this will lead , in the general case , to a messy result ; however , the above - mentioned simplest approximation , @xmath118 , which replaces eq .
( [ a^2 ] ) , yields a tractable expression , @xmath119 \right .
\nonumber \\ & & \left .
-\left ( 2+i\right ) \frac{\left\vert a_{y}\right\vert ^{2 } } { \left\vert a\right\vert } -2\varepsilon |a| \right\ } \label{lfinal}\end{aligned}\ ] ] [ the last term in the lagrangian density should be dropped if @xmath11 is small ; if @xmath11 is large , the term is obtained within the framework of the above scaling , that assumes small @xmath103 , hence @xmath120 .
finally , the action functional , @xmath121 , as defined by lagrangian ( [ lfinal ] ) , gives rise to the euler - lagrange equation , @xmath122 , in the following form : @xmath123 + 2\left\vert \tilde{a}\right\vert ^{2}\tilde{a}=0 ,   \label{new}\end{aligned}\ ] ] where @xmath124 [ this transformation eliminates a linear term in the equation generated by the last term in lagrangian ( [ lfinal ] ) ] .    formally setting @xmath125 in eq .
( [ new ] ) recovers an equation known as a model for the propagation of surface waves on a plasma layer with a sharp boundary @xcite . a class of nonpolynomial generalized nls equations similar to eq .
( [ new ] ) , and some of their fundamental solutions , such as single - soliton ones , were introduced in ref .
@xcite ( see also ref .
@xcite ) .    stationary solutions to eq . (
[ new ] ) with chemical potential @xmath31 are looked for in the ordinary form , @xmath126 , with @xmath127 satisfying equation @xmath128 + 2b^{3}+\mu b=0 $ ] , where @xmath129.@xmath113 this equation has its hamiltonian , @xmath130 setting @xmath131 , one obtains a family of soliton solutions , for any @xmath132 : @xmath133 for @xmath134 and @xmath132 , a family of periodic cnoidal - wave solutions is obtained in the form @xmath135 where @xmath136 is the jacobi s elliptic cosine , with modulus @xmath137 and @xmath138 are three roots of equation @xmath139 [ it is obtained by setting @xmath140 in eq .
( [ h ] ) ] .
solutions ( [ cn ] ) exist , for given @xmath132 , in a region of @xmath141 ; as said above , the limit of @xmath131 corresponds to soliton ( [ soliton ] ) , and the opposite limit , @xmath142    cnoidal solutions on the ring of radius @xmath98 must satisfy the corresponding periodic boundary conditions ( b.c . ) , @xmath143 ( for the ordinary 1d nls equations with repulsion and attraction , analytical solutions satisfying the periodic b.c .
are collected in refs . @xcite and @xcite , respectively ) . accordingly
, parameter @xmath144 in solutions ( [ cn ] ) , ( [ modulus ] ) is not continuous , but rather takes discrete values selected by matching the b.c . to the periodicity of @xmath145 , @xmath146 where @xmath147 is the complete elliptic integral , and
@xmath148 is an arbitrary integer .
note that eq .
( [ new ] ) features the galilean invariance : if @xmath149 is a solution , then its counterpart in the form of a solution moving at an arbitrary velocity @xmath150 is @xmath151 hence one can boost solitons ( [ soliton ] ) and cnoidal waves ( [ cn ] ) to a ( formally ) arbitrary velocity by means of transformation ( [ galileo ] ) .
in fact , the velocity is not arbitrary because the factor @xmath152 also must satisfy the periodic b.c . , which leads to a quantization condition , @xmath153 , with @xmath154 .
in the model with attraction ( @xmath21 ) , the existence of uniform ring - shaped solitons trapped in annular troughs of large radii is obvious , while a nontrivial issue is a possibility of finding a stability region for such solutions . in terms of eq .
( [ new ] ) , the ring soliton is represented by a cw ( continuous - wave ) solution with a constant amplitude @xmath155 , which was mentioned above as a limit case of the cnoidal solution family , with @xmath156 .
it is well known that , while all cw solutions of nls - like equations with attraction are modulationally unstable [ including nonpolynomial equations similar to eq .
( [ new ] ) @xcite ] , periodic b.c . may stabilize them if the cw amplitude is smaller than a certain critical value .
this argument suggests a possibility to find a stability region for the uniform ring solitons in the present model with attraction .
a standard analysis of the mi assumes a perturbation of the cw solution , @xmath157 \exp \left ( 2ib_{0}^{2}t+i\chi _ { 1}(y , t)\right ) , \ ] ] with infinitesimal amplitude and phase perturbations @xmath158 and @xmath159 .
eigenmodes of the perturbations may be looked for assuming that they are proportional to @xmath160 , where @xmath161 is an arbitrary real wavenumber , and @xmath162 is the corresponding instability growth rate .
straightforward calculations demonstrate that the cw solution to eq .
( [ new ] ) is subject to the mi under the condition @xmath163 > 0 .
\label{mi}\ ] ]    in the annular system that we are dealing with , @xmath161 is subjected to the geometric quantization , the same way as above , i.e. , @xmath164 with integer @xmath148 .
therefore , condition ( [ mi ] ) with lowest @xmath165 is not satisfied , making the cw solution _  modulationally stable _ , under the condition @xmath166 in a final form , this stability condition may be expressed in terms of the full norm of the axisymmetric ring soliton . substituting ansatz ( [ ansatz1 ] ) in eq .
( [ norm ] ) , and making use of the assumptions adopted above , one obtains @xmath167 .
further , substituting here the above approximation , @xmath168 , yields @xmath169 with regard to this relation , stability condition ( [ nomi ] ) amounts to a limitation on @xmath57 : @xmath170 ( note that this result does not contain any parameter ) .
a natural way to look at these stability conditions is to consider a situation with gradually increasing @xmath57 ( or @xmath155 ) .
the azimuthally uniform ring solution will loose its stability when the amplitude ( or norm ) attains the critical value defined by eq .
( [ nomi ] ) [ or by eq .
( [ thr ] ) ] .
one may expect that the onset of the mi gives rise to a _ bifurcation _ , which creates stable azimuthally modulated solutions , with a modulation depth scaling as @xmath171 for @xmath172 .
alternatively , one may fix the cw amplitude , @xmath155 , and gradually increase the trough s radius @xmath98 ; the loss of the cw stability and emergence of the modulated solutions should occur when @xmath98 attains a critical value following from eq .
( [ nomi ] ) , @xmath173    in fact , the modulated states generated by the mi onset should be nothing else but the stationary cnoidal wave ( [ cn ] ) , ( [ modulus ] ) .
indeed , with regard to eq .
( [ cw ] ) and the fact that @xmath174 , it follows form matching condition ( [ matching ] ) that , with the increase of @xmath98 for fixed @xmath155 , the cnoidal - wave solution appears , with an infinitely small modulation depth , at @xmath175 , where the critical value @xmath176 is _ precisely the same _ as defined by the onset of the mi [ i.e. , given by eq .
( [ cr ] ) ] .
a defect of the above approximate analysis is that , with regard to relation @xmath177 , stability condition ( [ thr ] ) may be cast in the form of @xmath178 , which _ does not _ comply with the underlying low - curvature assumption , @xmath179 .
this fact makes the existence of modulationally stable ring - shaped solitons in the present model with attraction doubtful . in direct simulations of eq .
( [ gpe2d ] ) with @xmath21 , we were unable to find modulationally stable annular uniform nonlinear states . in this connection
, it is relevant to recall that the linear limit of eq .
( [ gpe2d ] ) with potential @xmath180 does not allow any stationary solution localized in the radial direction , therefore the weak nonlinearity is essential in stabilizing the axisymmetric rings .
note that localized ring - shaped states may be possible in the linear limit of the bessel - lattice models considered in refs . @xcite and @xcite , with potentials @xmath181 and @xmath182 , respectively
; however , the azimuthal ( in)stability of such states was not studied in the presence of attraction . on the other hand ,
_ stable _ annular patterns ( both static ones and persistent ring - shaped breathers ) were recently found in an attractive model with a single circular trough , which was created not by the linear potential , but rather by modulation of the nonlinear coefficient in the radial direction @xcite . in that model ,
annular states may only be stable if the annulus is _ wide enough _ , with the ratio of its outer and inner radii exceeding @xmath183 .
vorticity - carrying rings were found to be unstable in the same setting .      in this section
we investigate the gap soliton nature and the band structure associated with ring shaped states in terms of the following 2d gpe @xmath184 , with @xmath185 the squared 2d angular momentum .
this equation corresponds to the @xmath186 limit ( @xmath187 ) of the spherically symmetric 3d gpe , with @xmath188 denoting spherical coordinates .
notice that apart for the factor of two in front of the term @xmath189 , due to spherical rather than cylindrical symmetry , eq .
( [ gpespher ] ) is the same as eq .
( [ gpe2d ] ) . from a physical point of view
( [ gpespher ] ) is appropriate to describe a condensate trapped in the outer part of the conical region created by a cylindrical beam of light passed through a very short focal distance lens . in the following
we shall refer to this setting as the `` pizza '' geometry since , in contrast with the pancake geometry , the thickness of the condensate in the axial direction is not uniform and reduces to zero as @xmath190 .
this geometry could be implemented by combining a cylindrical beam creating the radial lattice with a blue - detuned conical beam ( the pizza - shaper )    ) with @xmath191 , @xmath192 and @xmath193 .
right panel : the same as in the left panel for the repulsive case @xmath194 with @xmath195 and @xmath196.,title="fig:",width=151,height=151 ] ) with @xmath191 , @xmath192 and @xmath193 .
right panel : the same as in the left panel for the repulsive case @xmath194 with @xmath195 and @xmath196.,title="fig:",width=151,height=151 ]    launched in the opposite direction .
notice that at large radial distances ( i.e. in the outer part of the pizza ) the differences between pancake and pizza geometry become negligible especially if the thickness of the border is very small compared to the radius .
ring shaped states of eq .
( [ gpespher ] ) located in potential wells far from the center are expected , therefore , to be qualitatively similar to the ones of the pancake setting eq .
( [ gpe2d ] ) .    in the following
we look for ring gap solitons of eq .
( [ gpespher ] ) in the form of stationary states @xmath197 of the form @xmath198 by substituting eq .
( [ stat ] ) into eq.([gpespher ] ) we obtain the following equation for the radial function @xmath199 @xmath200 where @xmath201 are integer numbers which allow to satisfy the periodic boundary condition in @xmath202 : @xmath203 .
equation ( [ radial ] ) can be mapped into a 1d eigenvalue problem by means of the transformation @xmath204 leading to the following reduced radial equation @xmath205\xi = 0 \label{reducedradial}\ ] ] which has the form of a 1d schrdinger equation with periodic potential and nonlinear effective centrifugal barrier : @xmath206 ( here @xmath207 ) . for sufficiently large distances @xmath208 and fixed angular momentum and nonlinearity , this equation reduces to the mathieu equation with the well known band gap structure .
the part of the spectrum associated to excitations localized nearby the origin , however , is expected to be strongly influenced by the centrifugal barrier which at short distances enhances the nonlinearity and favors the formation of localized states .    to investigate the band structure and the gap solitons of eq .
( [ reducedradial ] ) we use a self - consistent method developed in ref .
the gap solitons obtained with the self - consistent method will then be used as initial conditions for time propagation of the radial 2d gpe in eq .
( [ gpe2d ] ) .     and corresponding to the energy levels @xmath209 of the left panel of fig .
[ fig8 ] , respectively .
panels c , d .
repulsive radial gap solitons ( thick lines ) trapped in circular troughs of the potential at @xmath210 and corresponding to the energy levels @xmath211 of the right panel of fig .
[ fig8 ] , respectively .
the thin lines corresponds to the radial potential @xmath212 scaled by a factor @xmath213 in the left panel and by a factor @xmath214 in the right panel , for graphical convenience.,title="fig:",width=132,height=132 ]   and corresponding to the energy levels @xmath209 of the left panel of fig .
[ fig8 ] , respectively .
panels c , d .
repulsive radial gap solitons ( thick lines ) trapped in circular troughs of the potential at @xmath210 and corresponding to the energy levels @xmath211 of the right panel of fig .
[ fig8 ] , respectively .
the thin lines corresponds to the radial potential @xmath212 scaled by a factor @xmath213 in the left panel and by a factor @xmath214 in the right panel , for graphical convenience.,title="fig:",width=132,height=128 ]     and corresponding to the energy levels @xmath209 of the left panel of fig .
[ fig8 ] , respectively .
panels c , d .
repulsive radial gap solitons ( thick lines ) trapped in circular troughs of the potential at @xmath210 and corresponding to the energy levels @xmath211 of the right panel of fig .
[ fig8 ] , respectively .
the thin lines corresponds to the radial potential @xmath212 scaled by a factor @xmath213 in the left panel and by a factor @xmath214 in the right panel , for graphical convenience.,title="fig:",width=132,height=132 ]   and corresponding to the energy levels @xmath209 of the left panel of fig .
[ fig8 ] , respectively .
panels c , d .
repulsive radial gap solitons ( thick lines ) trapped in circular troughs of the potential at @xmath210 and corresponding to the energy levels @xmath211 of the right panel of fig .
[ fig8 ] , respectively .
the thin lines corresponds to the radial potential @xmath212 scaled by a factor @xmath213 in the left panel and by a factor @xmath214 in the right panel , for graphical convenience.,title="fig:",width=132,height=132 ]    .
bottom panels .
3d view of the gap solitons corresponding to the bound states c ( left panel ) and d ( right panel ) of fig .
[ fig8].,title="fig:",width=151,height=113 ] .
bottom panels .
3d view of the gap solitons corresponding to the bound states c ( left panel ) and d ( right panel ) of fig .
[ fig8].,title="fig:",width=151,height=113 ]    .
bottom panels .
3d view of the gap solitons corresponding to the bound states c ( left panel ) and d ( right panel ) of fig .
[ fig8].,title="fig:",width=151,height=113 ] .
bottom panels .
3d view of the gap solitons corresponding to the bound states c ( left panel ) and d ( right panel ) of fig .
[ fig8].,title="fig:",width=151,height=113 ]    as a result we show the existence of radially symmetric gap solitons which correspond to the ring shaped states discussed before .    )
section of @xmath215 obtained from eq .
( [ gpe2d ] ) by taking as initial conditions the localized states depicted in panels a - d of fig .
the norm of the attractive states in panels @xmath209 are @xmath216 , @xmath217 , respectively , and correspond to bound states a , b , of fig .
repulsive states in panels @xmath211 have norm @xmath218 and @xmath219 , respectively , and correspond to the bound states c , d , of fig .
parameters are the same as in corresponding panels of fig .
[ fig8 ] . , title="fig:",width=132,height=132 ] ) section of @xmath215 obtained from eq .
( [ gpe2d ] ) by taking as initial conditions the localized states depicted in panels a - d of fig .
the norm of the attractive states in panels @xmath209 are @xmath216 , @xmath217 , respectively , and correspond to bound states a , b , of fig .
repulsive states in panels @xmath211 have norm @xmath218 and @xmath219 , respectively , and correspond to the bound states c , d , of fig .
parameters are the same as in corresponding panels of fig .
[ fig8 ] . , title="fig:",width=132,height=132 ]    ) section of @xmath215 obtained from eq .
( [ gpe2d ] ) by taking as initial conditions the localized states depicted in panels a - d of fig .
the norm of the attractive states in panels @xmath209 are @xmath216 , @xmath217 , respectively , and correspond to bound states a , b , of fig .
repulsive states in panels @xmath211 have norm @xmath218 and @xmath219 , respectively , and correspond to the bound states c , d , of fig .
parameters are the same as in corresponding panels of fig .
[ fig8 ] . , title="fig:",width=132,height=132 ] ) section of @xmath215 obtained from eq .
( [ gpe2d ] ) by taking as initial conditions the localized states depicted in panels a - d of fig .
the norm of the attractive states in panels @xmath209 are @xmath216 , @xmath217 , respectively , and correspond to bound states a , b , of fig .
repulsive states in panels @xmath211 have norm @xmath218 and @xmath219 , respectively , and correspond to the bound states c , d , of fig .
parameters are the same as in corresponding panels of fig .
, title="fig:",width=132,height=132 ]    in fig .
[ fig8 ] we depict the lower energy part of the spectrum obtained from eq .
( [ reducedradial ] ) in the case of attractive and repulsive interactions .
we see in both cases the existence of a series of discrete levels located in the gaps between the bands which correspond to radially symmetric localized states .
the wave functions of the energy levels labeled a - d in fig .
[ fig8 ] are shown in fig .
[ fig9 ] .
the symmetries of these localized states is similar to the ones of intrinsic localized modes of nonlinear lattices and gap solitons of continuous 1d gpe with a periodic potential . adopting the same terminology introduced for these cases we shall refer to them as the onsite symmetric and the onsite asymmetric ring gap solitons .
a 3d plot of these states is reported in fig .
[ fig10 ] ( notice that the asymmetric state gives two rings of matter located in the same trough of the potential in the attractive case and many ring oscillations in the repulsive case ) .
the stability of these gap solitons has been checked by time propagation under the 2d gpe in eq .
( [ gpe2d ] ) .
the results are reported in fig .
[ fig11 ] for the @xmath220 section of of @xmath221 as a function of time .
we see that for the attractive case the gap solitons are always unstable , the onsite symmetric state becomes modulationally unstable at @xmath222 while the asymmetric state decays into a state localized in the center and a series of ring solitons which expand trough the lattice keeping the same symmetry as the initial state . for the repulsive case
we see that while the asymmetric state in the second band gap ( panel d ) is unstable , the onsite symmetric ring gap soliton in the first band gap is quite stable under time evolution .
the stability of the onsite symmetric repulsive ring gap solitons and the modulational instability of the repulsive ones has been found also for other parameter values and seems to be a general property of the model .
the existence of other types of ring gap solitons ( which resemble intersite symmetric and asymmetric intrinsic localized modes of nonlinear lattices ) is also possible but they look even more unstable under time evolution than the onsite asymmetric ones for both attractive and repulsive cases . in the following section we shall investigate onsite symmetric gap soliton states in more details .
generic examples of a ring - shaped soliton in the attractive model , and its counterpart in the repulsive one , found from direct solution of gpe ( [ gpe2d ] ) in , respectively , imaginary and real time , are displayed in fig .
[ fig12 ] ( in the latter case , real - time simulations were run with absorbers set at @xmath223 and @xmath28 ) .
the latter example represents _ radial gap solitons _ with the core trapped in a remote circular trough ( cf .
the right panel in fig .
[ fig4 ] that shows a radial gap soliton with the core trapped in the center ) .
the fact that the soliton belongs to the gap type is clearly attested to by conspicuous inner and outer radial fringes attached to the core ( the fringes are absent in the ordinary soliton in the left panel of the figure ) .    ,
trapped in a circular trough at @xmath224 .
the amplitude of this soliton is @xmath225 , and its chemical potential and norm are @xmath226 and @xmath227 .
note that analytical approximation ( [ nfinal ] ) yields , in this case , @xmath228 , which is very close to the numerically found value , and its modulational instability is definitely predicted by eq .
( [ thr ] ) .
right panel : a stable ring - shaped soliton ( actually , it is a _ radial gap soliton _ ) in the repulsive model with @xmath55 , trapped at @xmath229 . in the latter case , @xmath230 and @xmath231.,title="fig:",width=151,height=113 ] , trapped in a circular trough at @xmath224 .
the amplitude of this soliton is @xmath225 , and its chemical potential and norm are @xmath226 and @xmath227 .
note that analytical approximation ( [ nfinal ] ) yields , in this case , @xmath228 , which is very close to the numerically found value , and its modulational instability is definitely predicted by eq .
( [ thr ] ) .
right panel : a stable ring - shaped soliton ( actually , it is a _ radial gap soliton _ ) in the repulsive model with @xmath55 , trapped at @xmath229 . in the latter case , @xmath230 and @xmath231.,title="fig:",width=151,height=113 ]    the evolution of the unstable axisymmetric ring from fig .
[ fig12 ] ( left panel ) , and outcome of the evolution are displayed in fig .
[ fig13 ] .
it is observed that the system tends to form a stationary necklace - like pattern composed of six peaks towering above the remaining quasi - uniform background .
( left panel ) , triggered by a small azimuthal perturbation , @xmath232 .
the evolution is shown in terms of the amplitude , @xmath233 , taken along circumference @xmath234 , with @xmath224 , @xmath235 running from @xmath236 to @xmath237 ( due to the symmetry , only half of the circle is shown ) .
( b ) the field profile formed by @xmath238 , demonstrating the emergence of a regular necklace - shaped pattern.,title="fig:",width=302,height=226 ] * ( a ) *     ( left panel ) , triggered by a small azimuthal perturbation , @xmath232 .
the evolution is shown in terms of the amplitude , @xmath233 , taken along circumference @xmath234 , with @xmath224 , @xmath235 running from @xmath236 to @xmath237 ( due to the symmetry , only half of the circle is shown ) .
( b ) the field profile formed by @xmath238 , demonstrating the emergence of a regular necklace - shaped pattern.,title="fig:",width=302,height=226]*(b ) *    in the repulsive model , a stable ring ( i.e. , the radial gap soliton ) may additionally carry , on its crests , pairs of dark solitons .
an example of such a pattern is displayed in fig .
[ fig14 ] ( it resembles the dipole - mode ring soliton recently found in the model with repulsion and bessel lattice in ref .
@xcite ) .    , trapped at @xmath229 .
this state has norm @xmath239 and chemical potential @xmath240.,width=302,height=226 ]    stable static rings with a density modulation in the azimuthal direction , corresponding to analytical solutions ( [ cn ] ) for the cnoidal waves , have been found too in the attractive model .
a typical example of such a pattern is displayed in fig .
[ fig15 ] ; in fact , it was generated as an outcome of the evolution of an unstable axisymmetric ring with the same norm .    ) .
the solution s norm and chemical potential are @xmath241 and @xmath242.,width=302,height=226 ]    increase of the norm leads to a transition from weakly modulated stable patterns to deeply modulated ones .
an example of such a stable configuration , generated from four gaussian pulses , with the phase shifts @xmath0 between them , by integration of eq .
( [ gpe2d ] ) in imaginary time , is displayed in fig .
[ fig16 ] ( note that the total norm of this pattern exceeds that corresponding to fig .
[ fig15 ] by a factor of @xmath243 ) .     and @xmath244 .
the norm of each soliton is @xmath245 , with the phase difference @xmath237 between adjacent ones .
the lattice is embedded in an integration domain of size @xmath246 .
, width=302,height=226 ]
strongly localized solitary waves of the attractive bec , self - trapped in a large - radius annular potential channel , like those shown in fig . [ fig16 ] , can freely move along in the channel . if set in motion with a sufficiently small velocity
, the soliton can circulate in the channel indefinitely long , preserving its integrity , as shown in fig .
[ fig17 ] .    , trapped in the circular trough at @xmath229 , was set in motion in the clockwise direction with velocity @xmath247 .
as illustrated by the sequence of snapshots , the circular motion in the trough continues indefinitely long without any visible loss.,width=302,height=302 ]    however , when the circulation speed is too large , the centrifugal force acting on the soliton can cause the matter to tunnel into adjacent radial troughs .
the loss of the norm ( number of atoms ) caused by the centrifugal tunneling eventually brings the soliton s norm below the threshold value necessary for its existence , which leads to disintegration of the soliton , see fig .
[ fig18 ] .
( with initial norm @xmath248 ) is set in motion with a greater velocity , @xmath249 , the centrifugal force gives rise to the underbarrier leakage of atoms into the outer circular troughs .
as illustrated by the sequence of snapshots , the decaying soliton can not preserve its shape .
the norm remaining within the given radial channel at @xmath250 is @xmath251 , i.e. , @xmath252 of the initial value.,width=302,height=302 ]    if the lattice is weak , and the soliton is strongly self - trapped ( with a sufficiently large norm ) , the centrifugal force can also drive it away from the center , across the lattice . in that case ( not shown here ) , the soliton also suffers the radiation loss , and eventually disintegrates .
the present model offers unique possibilities for exploring interactions and collisions between matter - wave solitons , as the dynamics in circular channels is free of external perturbations , which are always present in previously reported settings ( in the form of the weak longitudinal trap intended to keep the solitons within a finite spatial domain ) .
the same advantage is offered by setups with the toroidal trap @xcite
. however , sideline ( tangential ) collisions , which are possible between solitons moving in adjacent potential troughs of the radial lattice , have no counterparts in the toroidal traps .    in head - on collisions of solitons in the attractive model ( in the same circular channel ) ,
two competing processes play major roles .
the first is the nonlinear self - focusing , which may give rise to intrinsic collapse of a large - norm  lump " temporarily formed by the colliding solitons , and the other is the interference due to coherence of the matter - wave solitons @xcite .
the interference pattern , featuring alternating regions of high and low density , can suppress the collapse .
the predominance of either mechanism depends on the time of interaction .
if it is small ( in the case of large collision speed ) , solitons can pass trough each other without collapse ( although with conspicuous loss , see below , hence the collision can not be termed quasi - elastic ) , even if their total norm exceeds the collapse threshold . on the contrary ,
slow collisions lead to the onset of collapse .
these two possible scenarios of the head - on collisions between in - phase ( mutually attracting ) solitons are shown in figs . [ fig19 ] and [ fig20 ] .
figure [ fig19 ] presents an example of a relatively fast collision . at the moment of the full overlap between the solitons , @xmath253 ,
the interference fringes are evident , and are also seen in the cross section along the circumference of the channel ( lower right panel ) . although the total norm of the colliding solitons is overcritical ( @xmath254 ) , and the collision time ( @xmath255 ) is larger than the collapse time ( which is estimated as @xmath256 ) , the solitons separate without collapse .
nevertheless , the collision is essentially inelastic , the loss being additionally enhanced by the matter leakage under the action of the centrifugal force .
on the other hand , fig .
[ fig20 ] shows that the solitons colliding at a sufficiently small velocity merge together and indeed blow up due to the collapse .     and zero phase difference , in the circular trough at @xmath229 .
one of the solitons is set in motion ( at @xmath257 ) at @xmath258 with speed @xmath249 in the clockwise direction , while the second one ( at @xmath259 ) stays quiescent .
the norm remaining within the given circular channel at @xmath260 is @xmath261 .
the soliton which was originally at rest loses less norm than the moving one .
, title="fig:",width=151,height=151 ]   and zero phase difference , in the circular trough at @xmath229 .
one of the solitons is set in motion ( at @xmath257 ) at @xmath258 with speed @xmath249 in the clockwise direction , while the second one ( at @xmath259 ) stays quiescent .
the norm remaining within the given circular channel at @xmath260 is @xmath261 .
the soliton which was originally at rest loses less norm than the moving one . , title="fig:",width=151,height=151 ]     and zero phase difference , in the circular trough at @xmath229 .
one of the solitons is set in motion ( at @xmath257 ) at @xmath258 with speed @xmath249 in the clockwise direction , while the second one ( at @xmath259 ) stays quiescent .
the norm remaining within the given circular channel at @xmath260 is @xmath261 .
the soliton which was originally at rest loses less norm than the moving one .
, title="fig:",width=151,height=151 ]   and zero phase difference , in the circular trough at @xmath229 .
one of the solitons is set in motion ( at @xmath257 ) at @xmath258 with speed @xmath249 in the clockwise direction , while the second one ( at @xmath259 ) stays quiescent .
the norm remaining within the given circular channel at @xmath260 is @xmath261 .
the soliton which was originally at rest loses less norm than the moving one .
, title="fig:",width=151,height=151 ]       , in one circular channel ( @xmath229 ) , ends up with collapse .
upper panels : density plots at @xmath258 and the moment of collapse ( @xmath262 ) .
lower panel : the evolution of the density profiles in the cross section along the circumference of the potential trough . as in fig .
[ fig19 ] , the initial norm of each soliton is @xmath263 .
, title="fig:",width=151,height=151 ] , in one circular channel ( @xmath229 ) , ends up with collapse .
upper panels : density plots at @xmath258 and the moment of collapse ( @xmath262 ) .
lower panel : the evolution of the density profiles in the cross section along the circumference of the potential trough . as in fig .
[ fig19 ] , the initial norm of each soliton is @xmath263 . , title="fig:",width=151,height=151 ]    , in one circular channel ( @xmath229 ) , ends up with collapse .
upper panels : density plots at @xmath258 and the moment of collapse ( @xmath262 ) .
lower panel : the evolution of the density profiles in the cross section along the circumference of the potential trough . as in fig .
[ fig19 ] , the initial norm of each soliton is @xmath263 . ,
width=302,height=226 ]    two @xmath0 out - of - phase solitons colliding in the same potential trough bounce from each other due to repulsion between them . then , repeated bouncing collisions are observed , which is similar to the dynamics in the bessel lattice @xcite , as well as in a quasi-1d ol , equipped with a weak longitudinal trap @xcite .
however , in contrast to those works , we have observed that many collisions in the circular trough gradually wash out the phase coherence between the solitons . as a result ,
the solitons cease to repel each other , and a clear interference pattern disappears .
eventually , the originally repelling solitons merge together , shedding the excessive norm away , as shown in fig .
[ fig21 ] .
the intrinsic collapse _ does not _ occur in this case , as the loss makes the norm of the finally established single soliton smaller than the critical value at the collapse threshold .    , but
the phase difference of @xmath264 between them .
solitons bounce at @xmath265 and @xmath266 .
finally , the solitons merge into a single one , without collapse.,width=226,height=226 ]    collisions of matter - wave solitons in the attractive bec is accompanied by strong exchange of matter between them .
the present setting makes it possible to investigate related phenomena for solitons moving in adjacent channels of the radial potential , and thus colliding tangentially .
the result is that two in - phase solitons are unstable against the flow of matter in this setting , and tend to merge ( a similar trend was observed in the model with a quasi-1d lattice in ref .
nevertheless , the solitons may avoid collapse and separate , if the collision time is small due to a large relative velocity , as shown in fig .
[ fig22 ] .    , moving at velocities @xmath267 in adjacent troughs ( at @xmath229 and @xmath268 ) of the radial lattice with @xmath244 .
the solitons survive the collision , and remain localized , although the collision is inelastic , giving rise to radiation loss . , title="fig:",width=113,height=113 ] , moving at velocities @xmath267 in adjacent troughs ( at @xmath229 and @xmath268 ) of the radial lattice with @xmath244 .
the solitons survive the collision , and remain localized , although the collision is inelastic , giving rise to radiation loss . ,
title="fig:",width=113,height=113 ] , moving at velocities @xmath267 in adjacent troughs ( at @xmath229 and @xmath268 ) of the radial lattice with @xmath244 .
the solitons survive the collision , and remain localized , although the collision is inelastic , giving rise to radiation loss . , title="fig:",width=113,height=113 ]    on the other hand , if the solitons are loosely bound ( i.e. their norms are small and/or the radial lattice is weak ) , the slow lateral collision gives rise to merger of the solitons into one , excessive norm being shed off with linear waves , as demonstrated in fig .
[ fig23 ] .    ,
but for solitons with smaller norms , @xmath269 , launched to move at velocities @xmath270.,title="fig:",width=113,height=113 ] , but for solitons with smaller norms , @xmath269 , launched to move at velocities @xmath270.,title="fig:",width=113,height=113 ] , but for solitons with smaller norms , @xmath269 , launched to move at velocities @xmath270.,title="fig:",width=113,height=113 ]
in this work , we have proposed a new setting to explore two - dimensional ( 2d ) localized states in self - attractive and repulsive becs , in the form of a periodic radial optical lattice .
a crucial difference from previously studied 2d models with bessel lattices is the fact that no localized state may exist in the linear limit of the present model , hence all solitary states in it are  true solitons " , impossible without the nonlinearity . besides the bec
, the model may also apply to spatial solitons ( beams ) in photonic crystal fibers with a circular intrinsic structure .
the existence of such objects was demonstrated by means of different variants of the variational approximation , and in direct numerical simulations .
we have investigated the localized states trapped in the central potential well and in remote circular potential troughs . in both cases ,
a new species , namely , stable _ radial gap solitons _ , have been identified in the model with self - repulsion ( in addition , the radial gap soliton trapped in a circular trough can carry stable pairs of dark solitons on its crest ) . in remote troughs ,
we have investigated ring - shaped patterns delocalized in the azimuthal direction , as well as strongly localized azimuthal solitons ( including moving ones ) . in the attractive model , solutions of both types
are described by an effective 1d equation supplemented by the periodic boundary conditions ( it is a nonpolynomial nls equation with the coordinate running along the circumference ) . in particular ,
a finite threshold of the azimuthal modulational instability ( mi )  of axisymmetric ring - shaped states was predicted , and exact cnoidal - wave solutions , generated by the mi - induced bifurcation were found in the framework of the 1d equation .
azimuthal solitons were also found as solutions to the same equation .
the existence of stable azimuthally uniform and weakly modulated ring - shaped states was corroborated by direct simulations .
dynamics of completely localized solitons circulating in the annular potential troughs was investigated in the attractive model by means of direct simulations .
in particular , the solitons with sufficiently small velocities remain stable indefinitely long , while high velocities give rise to leakage of matter into the adjacent ( more remote ) trough under the action of the centrifugal force , which eventually destroys the soliton .
collisions between solitary waves running in the same or adjacent circular channels were investigated too .
head - on collisions of the in - phase solitons in one trough lead to the collapse ; @xmath0-out of phase solitons bounce from each other many times , but gradually loose their mutual coherence , and eventually merge into a single soliton without collapsing , shedding off excess norm .
in - phase solitons colliding in adjacent channels may also merge into a single soliton .
b.b.b . is partially supported by the grant no .
20 - 06 from the fund for fundamental research of the uzbek academy of sciences . the work of b.a.m .
was supported , in a part , by the center - of - excellence grant no .
8006/03 from the israel science foundation .
ms acknowledges financial support from the miur , through the inter - university project:_dynamical properties of bose - einstein condensates in optical lattices _ ,
prin-2005 .
99 j. denschlag , j. e. simsarian , d. l. feder , c. w. clark , l. a. collins , j. cubizolles , l. deng , e. w. hagley , k. helmerson , w. p. reinhardt , s. l. rolston , b. i. schneider , and w. d. phillips , science * 287 * , 97 ( 2000 ) ; s. burger , k. bongs , s. dettmer , w. ertmer , k. sengstock , a. sanpera , g. v. shlyapnikov , and m. lewenstein , phys .
lett . * 83 * , 5198 ( 1999 ) ;                                                                  b. b. baizakov , m. salerno , and b. a. malomed , in : _ nonlinear waves : classical and quantum aspects _ , ed . by f. kh .
abdullaev and v. v. konotop , p. 61
( kluwer academic publishers : dordrecht , 2004 ) ; also available at http://rsphy2.anu.edu.au/asd124/baizakov_2004_61_nonlinear waves.pdf .