bragg gratings ( bgs ) are light - controlling structures produced by a periodic variation of the refractive index along an optical fiber or waveguide .
devices based on fiber gratings , such as dispersion compensators , sensors , elements of laser cavities , etc . , are widely used in optical systems ka:99 .
gap solitons ( gss ) in fiber gratings , alias bg solitons stsi:94 , are supported through the balance between the bg - induced dispersion , which incorporates the bandgap in the system s linear spectrum , and the kerr nonlinearity of the fiber or waveguide .
analytical solutions for bg solitons in the standard fiber - grating model are well known acwa:89,chjo:89 .
slightly more than half of the analytical found family is stable , as first demonstrated , within the framework of the variational approximation , in ref .
@xcite , and then in a consistent numerical form in ref .
later , an accurate numerical technique for simulations of gss was elaborated in ref .
@xcite , and a rigorous analysis of their stability was developed in ref .
@xcite    following the theoretical prediction , bg solitons were created in the experiment , using a relatively short bg written in the cladding of an optical fiber @xcite and a virtual optically - induced grating coen:04 . in particular , an essential physical achievement was the creation of slow bg solitons moving at the velocity much lower than the velocity of light @xmath0 in vacuum ( namely , @xmath1 @xcite ) .
gss are one of fundamental species of solitons in optics , as well as in other nonlinear media .
in particular , the gs concept was extended to bose - einstein condensates ( becs ) , where such solitons were predicted in a condensate loaded into a periodic potential ( optical lattice ) bako:02,oski:03,br:04 , and then created in the experiment @xcite .
an issue of great significance to the fundamental studies and applications is the development of methods for control of bg solitons .
one of them is _ apodization _
@xcite , i.e. , the use of a grating with the bragg reflectivity gradually varying along the fiber ( or waveguide ) . in particular ,
a possibility was predicted to slow down the soliton and eventually bring it to a halt in a properly apodized bg mama:04 .
experimentally , it has been shown that the apodization helps to couple solitons into fiber gratings @xcite , and may also split them @xcite .
the retardation of the bg solitons to the above - mentioned low velocity , @xmath1 , was also performed with the help of the apodization @xcite .
technologies which make it possible to fabricate regular fiber gratings with periodic apodization , i.e. , an effective _ superlattice _ built on top of the bg , are well established @xcite .
analysis of the light propagation in the so developed _ supergratings _ was performed in @xcite .
it was shown that the supergrating gives rise to extra gaps in the system s spectrum (  rowland ghost gaps " ) .
solitons in the gaps were looked for in an approximation based on  supercoupled - mode equations " , which , essentially , assume that a soliton amounts to a slowly varying envelope of the supergrating s bloch function(s ) .
as is known , such an approximation ( the averaging method ) generally applies to the description of gss near bandgap edges @xcite .
a related problem was considered in @xcite , _
viz_. , gss in the bec loaded into a double - periodic optical lattice .
the period doubling opens an additional narrow  mini - gap " , where stable gss exist .
the objective of our work is to investigate the existence of gss in harmonic superlattices created on top of the usual bg in fibers or waveguides with the intrinsic kerr nonlinearity , extending the analysis initiated in yame:05 .
our model has periodic potentials in one spatial dimension .
analytical approaches for gss in two or more spatial dimensions were also developed recently in @xcite .
the paper is organized as follows .
we first give our mathematical model of the supergrating and analyze the bandgap structure in the model , using its dispersion relation , in section  [ s : prl ] .
this constitutes a non - trivial extension of the classical spectral theory for the mathieu equation smjo:99 , and predicts parameter values at which gss may be expected to exist in the model .
the sequence of new bandgaps not present in the unmodulated bg is produced . in sections  [
s : mel ] and [ s : av ] we establish the existence of gss in the newly opened bandgaps , using a combination of analytical and numerical techniques .
we find various families of gss in the supergrating , including those which do not exist in the ordinary gs model ( in particular , stable asymmetric solitons , and multi - soliton bound states ) .
using the melnikov method , we first prove that gss in the original bg model have counterparts in the extended model .
a similar treatment was used for a perturbed hamiltonian partial differential equations in @xcite .
then , the method of averaging is used to demonstrate the emergence of gss in the new bandgaps .
the most important result is that in each new bandgap , there is a single family of fundamental symmetric gss which _ entirely _ fills the gap .
such gss near band edges in nonlinear schrdinger equations with periodic potentials were also studied by different approaches in hwakya:11,pesuki:04 . in section  [
s : bif ] we discuss bifurcations of gss from a geometric point of view and illustrate this general approach in numerical computations for our problem .
the majority of the newly established gss are _ stable _ solutions of the model equations .
in particular , the family of fundamental gss in the new bandgaps is stable even for negative frequencies , when the ordinary gss are unstable @xcite .
we do not study the stability of gss in this paper in detail , referring to @xcite for the verification of the stability by means of direct simulations .
the standard model of the bragg grating ( bg ) in an optical fiber is based on the equations @xcite @xmath2 where @xmath3 and @xmath4 are amplitudes of the right- and left - propagating waves in the fiber . in these equations
, @xmath5 is the coordinate along the fiber , and @xmath6 is time . a family of exact soliton solutions of ( [ e : standard ] )
have been found in @xcite .    here
we will consider a situation where the bg is subjected to periodic modulation in @xmath5 , with a period @xmath7 , which gives rise to a _
supergrating_. the general form of the extended model reads @xmath8
v+\mu \cos \left ( kx+\delta \right ) ~u+\left ( |v|^{2}+\textstyle\frac{1}{2}|u|^{2}\right ) u = 0 , \\ & iv_{t}-iv_{x}+\left [ 1-\varepsilon \cos \left ( kx\right ) \right ] u+\mu \cos \left ( kx+\delta \right ) ~v+\left ( |u|^{2}+\textstyle\frac{1}{2}|v|^{2}\right ) v = 0.\end{split } \label{e : uv}\ ] ]    the real parameter @xmath9 accounts for the periodic modulation of the bg strength . in real fiber gratings , it may be typically @xmath10 cm ru:86,brst:97 , while the total length of the grating may be up to @xmath10 m. the modulation is implemented through the change imposed on the local variation of the refractive index in the fiber s cladding and may also affect the local index felt by each wave .
this is accounted for by the perturbation parameter @xmath11 and a phase shift @xmath12 . in practice , the harmonic superlattice corresponding to eqs .
( [ e : uv ] ) may be implemented if the optical interference pattern , which is used to burn the grating into the fiber , is created with spatial beatings .
we consider stationary gss , which are found via the ansatz @xmath13substituting this into eqs .
( [ e : uv ] ) , we obtain @xmath14 v+\mu \cos \left ( kx+\delta \right ) u+\left ( |v|^{2}+\textstyle\frac{1}{2}|u|^{2}\right ) u = 0 , \\ & \omega v - iv^{\prime } + \left [ 1-\varepsilon \cos ( kx)\right ] u+\mu \cos \left ( kx+\delta \right ) v+\left ( |u|^{2}+\textstyle\frac{1}{2}|v|^{2}\right ) v = 0,\end{split } \label{e : uvans}\]]where the prime stands for @xmath15 .
we first analyze the bandgap structure of ( [ e : uv ] ) by looking for solutions of the linearized system in the form of @xmath16 with @xmath17 being the propagation constant .
thus , we have to find periodic solutions of the linearization of ( [ e : uvans ] ) @xmath18 v+\mu \cos \left ( kx+\delta \right ) u = 0 , \\ & ( \omega + q)v - iv^{\prime } + \left [ 1-\varepsilon \cos ( kx)\right ] u+\mu \cos \left ( kx+\delta \right ) v = 0.\end{split } \label{e : uvlin}\]]note that in the unperturbed problem , with @xmath19 , the dispersion relation is given by @xmath20 , so that we find the well - known gap in the spectrum @xmath21 .
let us consider the case of @xmath22
. for @xmath23 non - trivial solutions of ( [ e : uvlin ] ) emerge due to parametric resonance , similar to the situation in the classical mathieu equation @xcite .
consequently , we expect new gaps to open up at points @xmath24 , @xmath25 .
we will denote the corresponding gaps by @xmath26 , while the central gap will be denoted by @xmath27 . in the central gap @xmath27 , gss for @xmath28
are found as perturbations of the soliton for @xmath29 , established in acwa:89
. we will use melnikov s method @xcite to prove the existence of gap solitons in section  [ s : mel ] .
for the non - central gaps , however , we will apply the averaging method to investigate the emergence of ( small ) gss for @xmath28 in section  [ s : av ] .
straightforward perturbation theory for parametrically excited systems yields approximation to the gaps for small @xmath9 .
for example , an approximation to the gap with @xmath30 close to @xmath31 , i.e. close to the unperturbed first gap at @xmath32 , can be obtained by looking for solutions of ( [ e : uvlin ] ) to lowest order approximation @xmath33 using this ansatz we derive at the solvability condition at order @xmath34 , @xmath35
\\ \qquad \qquad \qquad = \varepsilon ^{2}\left [ ( \omega ^{2}-q^{2})^{2}-k^{2}(\omega ^{2}+q^{2}+1)-1\right ] .\end{array } \label{e : solv}\]]observe that @xmath36 , since we are interested in the gap close to the unperturbed one with @xmath32 .
moreover , note that at the boundaries of the gaps periodic solutions with wavenumber @xmath37 or @xmath38 exist @xcite .
hence , being interested in solutions with @xmath39 we introduce @xmath40 and expand ( [ e : solv ] ) for small @xmath41 to obtain @xmath42 a new gap appears if this quadratic equation for @xmath43 has no real solutions .
thus we obtain as the condition for the gap , that @xmath44 in terms of the original parameter @xmath30 we find that the gaps open up at @xmath45 and their widths are @xmath46 similar expressions can be obtained for the higher - order gaps @xmath47 for @xmath48 showing that the gap width of the band corresponding to @xmath49 scales with @xmath50 .
we can treat the case of @xmath51 similarly .
these perturbative results can be extended numerically by detecting parameter values , for which eqs .
( [ e : uvlin ] ) possesses periodic solutions with wavenumber @xmath37 or @xmath38 .
it is well known that these parameters form the boundary of the gaps in the spectrum ( compare with above ) .
following @xcite we have computed these parameter values as solutions of the corresponding eigenvalue problem in matlab .
the results for several combinations of parameter values @xmath37 and @xmath52 are shown in fig .
f : linear . in this figure
the shaded and white areas stand for bands and gaps , respectively .
it is found that varying @xmath37 merely changes the width of bands and gaps with respect to @xmath53 , whereas changing @xmath52 leads to more complex effects . in particular , for @xmath54 gaps opening at @xmath29 can vanish again when @xmath55 or @xmath52 are increased .
in this and next sections we investigate the existence of gss in eq .
( e : uv ) , using both analytical and numerical techniques .
these solitons are described by homoclinic solutions to the origin @xmath56 in eqs .
( e : uvans ) .
before we present the computations let us take a closer look at properties of eqs .
( [ e : uvans ] ) . it will be important that these equations allow for a symmetry reduction by setting @xmath57 , where @xmath58 denotes complex conjugation . in the following
, we will thus only consider the reduced equation @xmath59 let @xmath60 , where @xmath61 .
( [ e : u ] ) is rewritten as @xmath62b -\frac{3}{2}(a^{2}+b^{2})b , \\ b^{\prime}=&[\omega-1+\varepsilon\cos kx+\mu\cos(kx+\delta)]a + \frac{3}{2}(a^{2}+b^{2})a . \end{split } \label{e : ab}\ ] ] we will consider as a dynamical system in the three - dimensional phase space @xmath63 , with @xmath64 as the circle of length @xmath65 .    note that if @xmath66 , then eqs .  ( [ e : ab ] ) are _ reversible _ , that is , invariant under compositions of time reversal and the ( linear ) involutions , @xmath67 we refer to @xcite for a general review of reversible systems
. of particular importance for us will be _ symmetric _ homoclinic orbits , which are mapped to itself under the action of @xmath68 or @xmath69 .
it is well known that an orbit is symmetric if and only if it intersects the invariant plane of the involution .
now we assume that @xmath70 and @xmath71 , and set @xmath72 for some @xmath73 .
the cases of @xmath74 and @xmath75 will be treated in section  [ s : av ] and a : bound , respectively .    for @xmath29 eqs .
( [ e : ab ] ) become a planar hamiltonian system , @xmath76 with a hamiltonian @xmath77+\frac{3}{8}(a^2+b^2)^2,\ ] ] where @xmath78 and @xmath79 represent the canonical momentum and coordinates , respectively .
the equilibrium at the origin is a hyperbolic saddle in ( e : ab0 ) and has a pair of homoclinic orbits @xmath80 which are symmetric under the involution @xmath81 . when @xmath55 and/or @xmath52 are nonzero but sufficiently small , the origin is still a hyperbolic saddle in and has two - dimensional stable and unstable manifolds which may intersect transversely .
such intersection yields transverse homoclinic orbits to the saddle at the origin , which persist under variation of the parameters @xcite . along certain curves in the parameter space ,
however , the intersection of the manifolds may become tangential , resulting in bifurcations of the homoclinic orbits .
see section [ s : bif ] for details .
here we aim to prove the existence of such transverse homoclinic orbits in eqs .
( [ e : ab ] ) by means of the melnikov method @xcite . from formula ( 4.5.6 ) in @xcite
, we derive the melnikov functions @xmath82 for @xmath83 as @xmath84\,a_{\pm } ( x)\,b_{\pm } ( x)\cos k(x+x_{0})dx   \notag \\ & + 2\mu \int_{-\infty } ^{\infty } a_{\pm } ( x)\,b_{\pm } ( x)\cos [ k(x+x_{0})+\delta ] dx   \notag \\ = & -\varepsilon \left ( \int_{-\infty } ^{\infty } [ 2\omega + 3(a_{\pm } ^{2}(x)+b_{\pm } ^{2}(x))]\,a_{\pm } ( x)\,b_{\pm } ( x)\sin kx\,dx\right ) \sin kx_{0 }   \notag \\ & -2\mu \left ( \int_{-\infty } ^{\infty } a_{\pm } ( x)\,b_{\pm } ( x)\sin kx\,dx\right ) \sin ( kx_{0}+\delta ) ,   \label{e : mel}\end{aligned}\]]where we used the fact that @xmath85 and @xmath86 are even and odd functions of @xmath5 , respectively . substituting ( [ e : homo0 ] ) into ( e : mel ) and using the method of residues ,
we compute @xmath87^{2}}\sin kx\,dx   \notag \\ & = \frac{\pi k}{3}\,\mathrm{cosech}\!\left ( \frac{k\pi } { 2\sin \theta } \right ) \sinh \!\left ( \frac{k\theta } { 2\sin \theta } \right )   \label{e : int1}\end{aligned}\]]and @xmath88a_{\pm } ( x)\,b_{\pm } ( x)\sin kx\,dx   \notag \\ & = \frac{8}{9}\sin ^{5}\theta \int_{-\infty } ^{\infty } \frac{\sinh ( 2x\sin \theta ) } { [ \cosh ( 2x\sin \theta ) + \cos \theta ] ^{3}}\sin kx\,dx   \notag \\ & = \frac{\pi k}{9}\,\mathrm{cosech}\!\left ( \frac{k\pi } { 2\sin \theta } \right ) \left [ -2\sinh \!\left ( \frac{k\theta } { 2\sin \theta } \right ) \cos \theta + k\cosh \!\left ( \frac{k\theta } { 2\sin \theta } \right ) \right ] .
\label{e : int2}\end{aligned}\]]see [ a : cal ] for derivations of ( [ e : int1 ] ) and ( [ e : int2 ] ) .
hence , the melnikov functions become @xmath89where we used the relation @xmath90 .
let us assume that @xmath91 then we easily see that @xmath92 has simple zeros at @xmath93 where @xmath94
this implies that there exist transverse homoclinic orbits near @xmath95 and @xmath96 . in particular , when @xmath66 , we have @xmath97 so that the first and second orbit are symmetric under @xmath81 and @xmath98 , respectively .
thus , we can prove that there exist gss with @xmath99 in ( [ e : uv ] ) .     for @xmath100 , @xmath101 and @xmath102 : ( a ) @xmath103 ; ( b ) @xmath104 .
, title="fig : " ] for @xmath100 , @xmath101 and @xmath102 : ( a ) @xmath103 ; ( b ) @xmath104 . , title="fig : " ]     of ( [ e : uv ] ) in the central gap with @xmath102 , @xmath105 and @xmath101 .
the solid curve shows the real part and the dashed curve the imaginary part of the solution .
, title="fig : " ] ( a ) @xmath106     of ( [ e : uv ] ) in the central gap with @xmath102 , @xmath105 and @xmath101 .
the solid curve shows the real part and the dashed curve the imaginary part of the solution .
, title="fig : " ] ( b ) @xmath107     of ( [ e : uv ] ) in the central gap with @xmath105 , @xmath102 , @xmath108 and @xmath109 .
, title="fig : " ] + ( a ) @xmath106     of ( [ e : uv ] ) in the central gap with @xmath105 , @xmath102 , @xmath108 and @xmath109 .
, title="fig : " ] + ( b ) @xmath107    to illustrate the above analysis we now introduce the poincar section @xmath110 in phase space , and consider the return map induced by the flow of .
we have computed the stable and unstable manifolds of the origin using the continuation tool ` auto ` @xcite with assistance of the ` hommap ` driver @xcite . in particular , to compute the stable and unstable manifolds , small segments ( @xmath111 ) in the stable and unstable subspaces of the origin were taken and solutions of ( [ e : ab ] ) with initial conditions on the segments were numerically integrated .
the same approach was also used for computation of the one - dimensional unstable manifolds for a three - dimensional poincar map in @xcite . as detected in the theory , the stable and unstable manifolds
are observed to intersect transversely and gss could be found .
the stable and unstable manifolds on the poincar section @xmath112 for @xmath113 , @xmath114 , @xmath115 and @xmath102 are drawn in fig .
[ f : im0 ] .
intersections of these manifolds give rise to homoclinic orbits or gss , and symmetric gss in ( [ e : uv ] ) for the same parameter values are plotted in fig .
[ f : big0 ] .
these gss are not visibly different from the unperturbed ones of . in fig .
[ f : big0 m ] asymmetric gss in are also plotted for @xmath116 , @xmath117 , @xmath114 , @xmath115 and @xmath102 . note that there are different asymmetric gss which have almost the same shapes as those in fig .
[ f : big0 m ] but have peaks at a position shifting by @xmath118
. plots of gss throughout this paper show the real part of the solutions as a solid curve , whereas the imaginary part is plotted as a dashed curve .
in order to investigate the existence of gss of ( [ e : uv ] ) with @xmath99 in the non - central gaps , we shall use the higher - order averaging method mu:88,mu:91 .
let us describe a general setup for the averaging analysis .
we seek homoclinic solutions to the origin of ( [ e : ab ] ) near @xmath119 in the limit of small @xmath120 .
note that for @xmath121 the origin is an equilibrium in ( [ e : ab ] ) , at which the jacobian matrix has a pair of purely imaginary eigenvalues @xmath122 when @xmath123 .
let @xmath124 and @xmath125 be non - negative integers which will be determined for each gap later . define a detuning parameter @xmath126 by @xmath127 let @xmath128 \displaystyle\frac{mk}{2(\omega _ { \pm m}+1)}\sin \frac{mk}{2}x & \displaystyle\cos \frac{mk}{2}x\end{pmatrix},\ ] ] which is the fundamental matrix of the linearization at the origin of the unperturbed system ( [ e : ab0 ] ) with @xmath129 . using the transformation @xmath130 \displaystyle\frac{m}{2(\omega _ { \pm m}+1)}\eta \end{pmatrix } \label{e : trans}\]]in ( [ e : ab ] ) , one obtains @xmath131b \\[1ex ] \lbrack \varepsilon ^{j}\omega + \cos kx+\bar{\mu}\cos ( kx+\delta ) + \frac{3}{2}\varepsilon ^{l}(a^{2}+b^{2})]a\end{pmatrix } ,   \label{e : xieta}\end{aligned}\]]where @xmath132 and @xmath78 and @xmath79 are represented by @xmath133 and @xmath134 via ( [ e : trans ] ) .
homoclinic solutions in ( [ e : xieta ] ) , equivalently in ( [ e : ab ] ) , can be approximately obtained by the averaging method .
these homoclinic solutions correspond to gps of ( [ e : uv ] ) with @xmath99 in the non - central gaps . while in the first gaps @xmath135
the standard ( i.e. , first - order ) averaging method is sufficient , we have to use the higher - order averaging method for the higher - order gaps .
this difficulty occurs because the potential considered here consists of the pure harmonic component @xmath136 .
but it is important to note that higher - order averaging is also needed for more general potentials with higher fourier components , if these higher components are small compared to the first one .
this is the reason for us to present the higher - order method here .
the necessary lengthy computations can be easily performed with computer algebra systems .
specifically , we have used the program ` haverage.m ` for the computer software ` mathematica ` @xcite , developed in @xcite ( see also @xcite ) .
we exemplarily present the results for the first three gaps @xmath135 , @xmath137 and @xmath138 .
similar second - order averaging analyses can also calculate the boundaries of the central gap ( see [ a : bound ] ) .      for the computations in the first gap at @xmath139 , we set @xmath140 and @xmath141 .
the standard averaging procedure for ( [ e : xieta ] ) yields the first - order averaged system @xmath142 , \\ & \eta^{\prime}=\varepsilon\left [ \frac{\omega_{\pm 1}(1 + 2\omega)+\bar{\mu}\cos\delta}{k}\xi + \frac{\bar{\mu}\sin\delta}{k}\eta + \gamma_{\pm 1}(\xi^2+\eta^2)\xi\right ] , \end{split } \label{e : av1}\ ] ] where @xmath143 eqs .
are easily obtained by averaging the right hand side of with respect to @xmath5 over @xmath144 $ ] and such computations can be done by the ` mathematica ` program ` haverage.m ` . in the averaged system ( [ e : av1 ] ) with @xmath145 , we perform a rotational transformation @xmath146 to obtain a system of a more convenient form @xmath147 , \\ & \tilde{\eta}^{\prime}=\varepsilon\left [ \frac{2\omega_{\pm 1}}{k}(\omega_{1\pm}+\omega)\tilde{\xi } + \gamma_{\pm 1}(\tilde{\xi}^2+\tilde{\eta}^2)\tilde{\xi}\right ] , \end{split } \label{e : av1a}\ ] ] where @xmath148 and @xmath149 with @xmath150.\ ] ] note that eq .
( [ e : av1a ] ) is still valid for @xmath151 if @xmath152 is taken when @xmath153 and if @xmath154 is taken when @xmath155 .
the system ( [ e : av1a ] ) is hamiltonian with a hamilton function @xmath156 + \frac{1}{4}\gamma_{\pm 1}(\tilde{\xi}^2+\tilde{\eta}^2)^2\right\},\end{aligned}\ ] ] for which the origin @xmath157 is a saddle if @xmath158 , i.e. , @xmath159 we easily see that if @xmath160 then @xmath161 ( resp .
@xmath162 ) and a closing of the gap region @xmath163 ( resp .
@xmath164 ) occurs .
see also fig .
f : linear(d ) .
if @xmath165 , then @xmath166 , and thus such a closing does not occur .
let us assume that @xmath166 , i.e. , condition  ( [ e : closing ] ) does not hold .
a straightforward analysis of the level sets of @xmath167 shows the existence of two homoclinic orbits in ( [ e : av1a ] ) and hence in ( e : av1 ) .
moreover , we can obtain analytical expressions of the homoclinic orbits as @xmath168 } { \cosh[(4\omega_{\pm 1}\omega_{1\pm}\sin\phi / k)\varepsilon x]+\cos\phi } , \right .
\notag \\ & \qquad \left .
\mp 2\sqrt{\frac{2\omega_{\pm 1}\omega_{1\pm}}{3k\gamma_{\pm 1}}}\sin\phi \cos\frac{1}{2}\phi\ , \frac{\cosh[(2\omega_{\pm 1}\omega_{1\pm } \sin\phi / k)\varepsilon x ] } { \cosh[(4\omega_{\pm 1}\omega_{1\pm } \sin\phi / k)\varepsilon x]+\cos\phi } \right ) ,   \label{e : homo}\end{aligned}\ ] ] where @xmath169 .
it is easy to see that , since the coefficient in front of @xmath170 in @xmath167 is negative , the homoclinic orbits intersect the @xmath171-axis and are symmetric about it , forming a figure - of - eight configuration . in general , it is very difficult to make precise statements about the symmetry of homoclinic orbits in ( [ e : av1 ] ) .
for some special cases , however , we have the following results :    1 .   when @xmath172 and @xmath173 , the homoclinic orbits in the averaged system ( [ e : av1 ] )
are symmetric about the @xmath174-axis if @xmath175 and about the @xmath176-axis if @xmath177 in the negative gap , while they are always symmetric about the @xmath174-axis in the positive gap .
2 .   when @xmath178 and @xmath179 , the homoclinic orbits in ( [ e : av1 ] )
are symmetric about the @xmath174-axis if @xmath180 and about the @xmath176-axis if @xmath181 in the positive gap , while they are always symmetric about the @xmath174-axis in the negative gap .
3 .   when @xmath22 , the homoclinic orbits in ( [ e : av1 ] )
are symmetric about the @xmath174-axis in both the positive and negative gaps .    from the existence of homoclinic orbits in the averaged system ,
we claim that there exist homoclinic orbits in and hence stationary gss in . for this consider ( [ e : xieta ] ) and ( e : av1 ) as dynamical systems defined in the phase space @xmath63 .
recall that we set @xmath140 and @xmath141 in ( [ e : xieta ] ) .
we introduce a cross section @xmath182 at @xmath183 and define the poincar maps @xmath184 for ( [ e : xieta ] ) and ( [ e : av1 ] ) , respectively . for both @xmath185 and @xmath186 the origin @xmath187 is a fixed point and has stable and unstable manifolds .
moreover , the stable and unstable manifolds , @xmath188 and @xmath189 , coincide along the homoclinic orbits in ( [ e : av1 ] ) for @xmath186 . by the averaging theorem ( see theorem  4.1.1 of @xcite ) , the stable and unstable manifolds , @xmath190 and @xmath191 , for @xmath185 are @xmath192-close to @xmath188 and @xmath193 . since eq .
( e : xieta ) as well as eq .
( [ e : av1 ] ) are hamiltonian , the poincar maps @xmath194 are area - preserving .
( solid lines ) which are close to @xmath195 ( broken lines ) but do not intersect .
the regions @xmath196 and @xmath197 with @xmath198 are hatched distinctly . ]
now assume that the stable and unstable manifolds @xmath199 and @xmath200 do not intersect , and go inward and outward , respectively , as shown in fig .  [
f : asu ] .
we want to show that this yields a contradiction and define a region @xmath196 , which is encircled by two parts of @xmath199 and @xmath200 and a line connecting two points @xmath201 and @xmath202 on @xmath199 and @xmath200 .
the region @xmath196 is mapped by the poincar map @xmath203 to another region @xmath204 which is encircled by two different parts of @xmath205 and @xmath200 and a curve connecting two points @xmath206 and @xmath207 on @xmath199 and @xmath200 .
it is obvious that @xmath208 is larger than @xmath196 and contradicts the fact that @xmath203 is area - preserving .
similarly , we can not assume that @xmath199 and @xmath200 go outward and inward , respectively .
therefore , we can prove that @xmath199 and @xmath200 intersect and there exists a homoclinic orbit in the full system ( e : xieta ) and hence in ( [ e : u ] ) .
on the other hand , eq .   can be averaged up to @xmath209 for any integer @xmath210 but has very rapid oscillations compared with the unperturbed homoclinic orbits which vary very slowly with @xmath211 . applying a result by neishtadt @xcite
( see also @xcite ) , we see that the splitting distance between @xmath199 and @xmath200 is exponentially small with respect to @xmath212 at most . thus , there exist stationary gs solutions to the system of equations ( [ e : uv ] ) in the gap regions .    moreover , assume that @xmath66 , i.e. , @xmath22 or @xmath213 , so that eqs .
( [ e : ab ] ) is reversible with respect to @xmath214 and @xmath215 , @xmath216 .
then it follows from the above arguments that the homoclinic orbit in the averaged system ( [ e : av1 ] ) is symmetric about the @xmath176-axis ( resp .
@xmath174-axis ) . by the persistence of symmetric orbits in reversible systems
, there exist symmetric homoclinic orbits about the @xmath176-axis ( resp .
@xmath174-axis ) in the full system ( [ e : xieta ] ) and with respect to @xmath81 and @xmath217 ( resp . @xmath98 and @xmath218 ) in ( [ e : ab ] ) and hence in ( [ e : u ] ) . in particular ,
the homoclinic orbits in ( [ e : u ] ) have an even ( resp .
odd ) real part @xmath78 and odd ( resp . even )
imaginary part @xmath219 .
we can use similar arguments to the central gap when condition  ( [ e : con ] ) does not hold , and prove that there must be homoclinic orbits .
however , the splitting distance between the stable and unstable manifolds is not exponentially small with respect to @xmath55 but @xmath220 at most since both of the perturbations in and unperturbed homoclinic orbits vary with @xmath221 .     for @xmath179 . ]     for the associated poincar map : ( a ) @xmath222 ; ( b ) @xmath223 ; ( c ) @xmath224.,title="fig : " ] for the associated poincar map : ( a ) @xmath222 ; ( b ) @xmath223 ; ( c ) @xmath224.,title="fig : " ] +   for the associated poincar map : ( a ) @xmath222 ; ( b ) @xmath223 ; ( c ) @xmath224.,title="fig : " ]    in addition , we can detect a non - transverse homoclinic orbit in ( [ e : ab ] ) when @xmath179 and @xmath225 as follows . fix @xmath179 . then the melnikov functions @xmath92 have graphs as shown in fig .
[ f : mel ] , depending on @xmath226 , @xmath227 or @xmath228 . since @xmath229 represent signed measures of their distance ( see , e.g. , section  4.5 of @xcite ) , noting that @xmath230 and eqs .
( [ e : ab ] ) is reversible , we can draw the behavior of the stable and unstable manifolds @xmath231 for the poincar map of ( [ e : ab ] ) as shown in figs .
[ f : su](a ) and ( b ) .
using the fact that the poincar map is area - preserving and applying a discussion similar to the above one for ( e : xieta ) , we see that @xmath190 and @xmath191 coincide for some @xmath55 near @xmath232 , say @xmath233 , as shown in fig .
f : su(c ) .
thus , the point @xmath234 is very degenerate in the parameter space .
this also suggests that new interesting behavior may occur at the point in the original pdes ( [ e : uv ] ) .     when @xmath108 , @xmath104 , @xmath235 and @xmath102 : ( a ) @xmath236 ; ( b ) @xmath237 ; ( c ) @xmath238 .
the red and blue curves represent the stable and unstable manifolds , respectively.,title="fig : " ] when @xmath108 , @xmath104 , @xmath235 and @xmath102 : ( a ) @xmath236 ; ( b ) @xmath237 ; ( c ) @xmath238 .
the red and blue curves represent the stable and unstable manifolds , respectively.,title="fig : " ] +   when @xmath108 , @xmath104 , @xmath235 and @xmath102 : ( a ) @xmath236 ; ( b ) @xmath237 ; ( c ) @xmath238 .
the red and blue curves represent the stable and unstable manifolds , respectively.,title="fig : " ]    .
the parameter values @xmath239 of fig .
[ f : im0](b ) with @xmath102 is plotted as @xmath240 . ]    in fig .
[ f : im0 ] we show the results of numerical computations of the stable and unstable manifolds by ` auto97 ` with ` hommap ` near the homoclinic bifurcation point @xmath241 for @xmath242 , @xmath107 and @xmath102 .
we see that as predicted by the theory , the invariant manifolds almost coincide at @xmath243 while they split when @xmath55 is different from @xmath232 .
thus , the theoretical prediction for the homoclinic bifurcations is very precise .
figure  [ f : hbc ] also shows the approximate homoclinic bifurcation curves given by eq . ,
on which the stable and unstable manifolds coincide , in the @xmath244-space for @xmath245 .
we set @xmath246 and @xmath247 , and perform a second - order averaging method to prove the existence of gap solitons of ( [ e : uvans ] ) with @xmath99 in the second gaps appearing at @xmath248 .
to avoid the complexity of expressions , we give our results only for @xmath22 in the gaps of higher orders .
the second - order averaged system has been obtained using the ` mathematica ` program ` haverage.m ` , and it reads @xmath249 \right .
\\ & \qquad + \beta_{\pm 2}\left[k^{2}\xi^{2}+ ( \omega_{\pm 2}+1)^{2 } \eta^{2}\right]\big\}\ ; \eta , \\
\eta^{\prime}= & \displaystyle{\frac{\varepsilon^{2}}{(\omega_{\pm 2}+1)^{2 } } } \left\{\alpha_{\pm 2 } \left[(6(\omega_{\pm 2}+1)\omega+1)k^{2}- 5(\omega_{\pm 2}+1)^{2}\right ] \right .
\\ & \qquad + \beta_{\pm 2}\left[k^{2}\xi^{2 } + ( \omega_{\pm 2}+1)^{2}\eta^{2}\right]\big\ } \ ; \xi , \end{split } \label{e : av2}\ ] ] with @xmath250 } { 16k^{2}(\omega_{\pm 2}+1)^{2}}.\end{aligned}\ ] ] note that @xmath251 and @xmath252 .
( e : av2 ) is a hamiltonian system with the hamilton function @xmath253k^{2}-5(\omega_{\pm 1/1}+1)^{2}\ } } { 2(\omega_{\pm 1/1}+1)^{2}}\xi^{2 }   \notag \\ &
\quad + \frac{\alpha_{\pm 1/1 } \{[6(\omega_{\pm 1/1}+1)\omega-5]k^{2}+(\omega_{\pm 1/1}+1)^{2}\ } } { 2k^{2}}\eta^{2 }   \notag \\ &
\quad + \frac{\beta_{\pm 1/1}}{4k^{2}(\omega_{\pm 2/1}+1)^{2 } } [ k^{2}\xi^{2}+(\omega_{\pm 1/1}+1)^{2}\eta^{2}]^{2}\biggr).\end{aligned}\ ] ]    let @xmath254 the origin is a saddle of ( [ e : av2 ] ) and has a pair of homoclinic orbits if @xmath255 and @xmath256 when the signs ` @xmath257 ' and ` @xmath258 ' are taken in the subscript of @xmath53 in ( [ e : av2 ] ) , respectively , i.e. , @xmath259 or @xmath260    again , we find the averaged system to be hamiltonian with a saddle equilibrium at the origin , and a study of the corresponding level set of the hamiltonian shows the existence of a figure - of - eight of homoclinic orbits . these orbits intersect the @xmath174-axis . by the same argument as above , this implies the existence of gap solitons for ( [ e : u ] ) in the gaps @xmath137 having an even @xmath78-component and an odd @xmath79-component .
we set @xmath261 , @xmath262 and @xmath263 , and perform the third - order averaging procedure for ( [ e : xieta ] ) using the ` mathematica ` program ` haverage.m ` to obtain the averaged system @xmath264k^{2 } -4(\omega_{\pm 3}+1)^{2}\right\}\eta \\ & + \frac{\varepsilon^{3}}{9k^{2 } } \left\{\alpha_{\pm 3}-\beta_{\pm 3}\left[9k^{2}\xi^{2 } + 4(\omega_{\pm 3}+1)^{2}\eta^{2}\right]\right\}\eta , \\
\eta^{\prime}= & \frac{\varepsilon^{2}}{4(\omega_{\pm 3}+1)^{2}}\rho_{\pm 3 } \left\{[32(\omega_{\pm 3}+1)\omega-9]k^{2 } -4(\omega_{\pm 3}+1)^{2}\right\}\xi \\ & + \frac{\varepsilon^{3}}{4(\omega_{\pm 3}+1)^{2 } } \left\{\alpha_{\pm 3}+\beta_{\pm 3}\left[9k^{2}\xi^{2 } + 4(\omega_{\pm 3}+1)^{2}\eta^{2}\right]\right\}\xi , \end{split } \label{e : av3o}\ ] ] where @xmath265 [ 81k^{4}-56(\omega_{\pm3}+1)^{2}k^{2}+16(\omega_{\pm3}+1)^{4 } ] } { 512k^{4}(\omega_{\pm 3}+1)^{2 } } ,
\\ \beta_{\pm 3 } = & \frac{81k^{4}+24(\omega_{\pm 3}+1)^{2}k^{2}+16(\omega_{\pm 3}+1)^{4 } } { 64k^{2}(\omega_{\pm 3}+1)^{2}}.\end{aligned}\ ] ] note that @xmath266 and @xmath267 since @xmath268^{2}+18(\omega_{\pm3}+1)^{2}>0.\end{aligned}\ ] ]    let @xmath269 where @xmath270 note that @xmath271 and @xmath272 .
( [ e : av3o ] ) becomes @xmath273\}\eta , \\
\eta^{\prime}= & \frac{\varepsilon^{3}}{4(\omega_{\pm 3}+1)^{2}}\ { ( \gamma_{\pm 3}\nu+\alpha_{\pm 3 } ) + \beta_{\pm 3}[9k^{2}\xi^{2}+4(\omega_{\pm 3}+1)^{2}\eta^{2}]\}\xi , \end{split } \label{e : av3}\ ] ] where @xmath274>0.\ ] ] eq .
( [ e : av3 ] ) is a hamiltonian system with the hamilton function @xmath275^{2}\biggr ) .
\label{e : h3}\end{aligned}\ ] ]    we easily see that @xmath276 and set @xmath277 the origin is a saddle and has a pair of homoclinic orbits if @xmath278 , i.e. , @xmath279 especially , in the regions ( [ e : region3 ] ) the coefficients of @xmath280 and @xmath281 in the hamiltonian ( [ e : h3 ] ) are positive and negative , respectively , so that the pair of homoclinic orbits draws a vertical figure - eight in the @xmath282-phase plane .
similar computations reveal the existence of gap solitons in the higher order gaps .
these solutions are odd in @xmath78 ( and even in @xmath79 ) or vice versa , depending on whether the gap number is odd or even .
a general proof of this is beyond the scope of this paper . in the following we illustrate and extend this perturbative results by numerical computations .       when @xmath101 and @xmath102 : ( a ) @xmath283 , @xmath284 ; ( b ) @xmath285 , @xmath284 ; ( c ) @xmath286 , @xmath287 ; ( d ) @xmath288 , @xmath287 ; ( e ) @xmath289 , @xmath290 ; ( f ) @xmath291 , @xmath290 .
the red and blue curves represent the stable and unstable manifolds , respectively.,title="fig : " ] when @xmath101 and @xmath102 : ( a ) @xmath283 , @xmath284 ; ( b ) @xmath285 , @xmath284 ; ( c ) @xmath286 , @xmath287 ; ( d ) @xmath288 , @xmath287 ; ( e ) @xmath289 , @xmath290 ; ( f ) @xmath291 , @xmath290 .
the red and blue curves represent the stable and unstable manifolds , respectively.,title="fig : " ] +   when @xmath101 and @xmath102 : ( a ) @xmath283 , @xmath284 ; ( b ) @xmath285 , @xmath284 ; ( c ) @xmath286 , @xmath287 ; ( d ) @xmath288 , @xmath287 ; ( e ) @xmath289 , @xmath290 ; ( f ) @xmath291 , @xmath290 .
the red and blue curves represent the stable and unstable manifolds , respectively.,title="fig : " ] when @xmath101 and @xmath102 : ( a ) @xmath283 , @xmath284 ; ( b ) @xmath285 , @xmath284 ; ( c ) @xmath286 , @xmath287 ; ( d ) @xmath288 , @xmath287 ; ( e ) @xmath289 , @xmath290 ; ( f ) @xmath291 , @xmath290 .
the red and blue curves represent the stable and unstable manifolds , respectively.,title="fig : " ] +   when @xmath101 and @xmath102 : ( a ) @xmath283 , @xmath284 ; ( b ) @xmath285 , @xmath284 ; ( c ) @xmath286 , @xmath287 ; ( d ) @xmath288 , @xmath287 ; ( e ) @xmath289 , @xmath290 ; ( f ) @xmath291 , @xmath290 .
the red and blue curves represent the stable and unstable manifolds , respectively.,title="fig : " ] when @xmath101 and @xmath102 : ( a ) @xmath283 , @xmath284 ; ( b ) @xmath285 , @xmath284 ; ( c ) @xmath286 , @xmath287 ; ( d ) @xmath288 , @xmath287 ; ( e ) @xmath289 , @xmath290 ; ( f ) @xmath291 , @xmath290 .
the red and blue curves represent the stable and unstable manifolds , respectively.,title="fig : " ]     + ( a ) @xmath292 , @xmath293 , @xmath22     of ( [ e : uvans ] ) in the gaps @xmath135 for @xmath102 , symmetric with respect to @xmath217 or @xmath218 .
gray curves in the panels correspond to solutions of the averaged system ( [ e : av1]).,title="fig : " ] + ( b ) @xmath294 , @xmath293 , @xmath22     of ( [ e : uvans ] ) in the gaps @xmath135 for @xmath102 .
gray curves in the panels correspond to solutions of the averaged system ( [ e : av1 ] ) .
, title="fig : " ] + ( a ) @xmath292 , @xmath293 , @xmath295 , @xmath296     of ( [ e : uvans ] ) in the gaps @xmath135 for @xmath102 .
gray curves in the panels correspond to solutions of the averaged system ( [ e : av1 ] ) .
, title="fig : " ] + ( b ) @xmath294 , @xmath293 , @xmath295 , @xmath296     of ( [ e : uvans ] ) in the gaps @xmath135 for @xmath297 , @xmath298 , @xmath299 and @xmath300 .
, title="fig : " ] + ( a ) @xmath301     of ( [ e : uvans ] ) in the gaps @xmath135 for @xmath297 , @xmath298 , @xmath299 and @xmath300 .
, title="fig : " ] + ( b ) @xmath302    we computed the stable and unstable manifolds of the origin and homoclinic orbits in ( [ e : ab ] ) , i.e. , gss with @xmath99 in ( [ e : uv ] ) , using the computer package , ` auto ` with ` hommap ` . in our computations we found gss due to the transverse intersection of the stable and unstable manifolds , thus confirming the analytical results . for several parameter values ,
the stable and unstable manifolds for the poincar section @xmath303 are shown in fig .
[ f : im123 ] and gss in ( [ e : uv ] ) are plotted in figs .
[ f : big123]-[f : big1 ] .
in addition to the solutions of ( [ e : ab ] ) , the figures include homoclinic solutions of the averaged systems ( [ e : av1 ] ) , ( [ e : av2 ] ) and ( [ e : av3 ] ) , rescaled according to the transformation  ( [ e : trans ] ) .
as can be seen , they provide a good match to the envelope of the solutions to ( [ e : ab ] ) . in all figures
presented , the parameter @xmath37 has been set to @xmath102 .    in closing this section , we return to the parameter values of fig .
f : linear(d ) . in this case a new bandgap opens up for small @xmath55 , but closes at @xmath304 .
see also formula  .
our computations reveal that gss exist in the gaps both above and below this parameter value . in fig .
[ f : big1a ] we show solution profiles at @xmath305 and @xmath306 .
we further discuss the behavior of gss under variation of the system parameters .
when the parameter values are varied , homoclinic tangencies , i.e. , homoclinic bifurcations , can occur , leading to the creation of additional gss and bifurcations of gss .
we have found such bifurcations in all gap regions including the central one .    ) in the gaps @xmath164 and @xmath27 for @xmath102 and @xmath109 : ( a ) and ( b ) @xmath283 , @xmath101 ; ( c ) @xmath283 , @xmath307 ; ( d ) @xmath104 , @xmath101 . in plates  ( b ) and ( c ) , the point `` @xmath240 '' represents pitchfork and saddle - node bifurcations , respectively , title="fig:",scaledwidth=46.0% ] ) in the gaps @xmath164 and @xmath27 for @xmath102 and @xmath109 : ( a ) and ( b ) @xmath283 , @xmath101 ; ( c ) @xmath283 , @xmath307 ; ( d ) @xmath104 , @xmath101 . in plates  ( b ) and ( c ) , the point `` @xmath240 '' represents pitchfork and saddle - node bifurcations , respectively , title="fig:",scaledwidth=46.0% ] + ) in the gaps @xmath164 and @xmath27 for @xmath102 and @xmath109 : ( a ) and ( b ) @xmath283 , @xmath101 ; ( c ) @xmath283 , @xmath307 ; ( d ) @xmath104 , @xmath101 . in plates  ( b ) and ( c ) , the point `` @xmath240 '' represents pitchfork and saddle - node bifurcations , respectively , title="fig:",scaledwidth=46.0% ] ) in the gaps @xmath164 and @xmath27 for @xmath102 and @xmath109 : ( a ) and ( b ) @xmath283 , @xmath101 ; ( c ) @xmath283 , @xmath307 ; ( d ) @xmath104 , @xmath101 . in plates  ( b ) and ( c ) , the point `` @xmath240 '' represents pitchfork and saddle - node bifurcations , respectively , title="fig:",scaledwidth=46.0% ] +    as an example , we plot bifurcation diagrams of gss for eqs .
( [ e : uvans ] ) in the gaps @xmath164 and @xmath27 when @xmath292 or 0.5 and @xmath296 , in fig .
[ f : bd ] . here
@xmath308 represents the point at which each gs crosses the section @xmath309 . in the diagrams we observe saddle - node and pitchfork bifurcations of gss
the occurrence of these bifurcations is explained from the behavior of the poincar map of ( [ e : ab ] ) for the section @xmath182 , as follows .
+     let us first assume that @xmath54 , such that the poincar map is not reversible . in this case ,
a tangency between the stable and unstable manifolds as shown in fig .
[ f : hbmap](a ) leads to a saddle - node bifurcation , in which two new gss are created .
see fig .
[ f : bd](c ) for a numerical evidence .    now suppose that @xmath22 and let us choose a cross - section , such that the poincar map is reversible with respect to @xmath81 and @xmath98 .
we can distinguish two different bifurcation scenarios .
if the stable and unstable manifolds have a cubic tangency on the invariant plane of @xmath81 or @xmath98 , but intersect this plane transversely as shown in fig .
[ f : hbmap](b ) , then a pitchfork bifurcation of homoclinic orbits takes place . when the parameters are varied , a pair of homoclinic orbits appear , while the original homoclinic orbit still exists due to the persistence of symmetric orbits .
see figs .
[ f : bd](a ) and ( b ) .
the original gs is symmetric , but the two new orbits are asymmetric .    on the other hand , if the stable and unstable manifolds are tangent to the invariant plane at their quadratic tangency as shown in fig .
[ f : hbmap](c ) , then a saddle - node bifurcation of symmetric gss occurs .
see fig .
f : bd(b ) . note that a saddle - node bifurcation of asymmetric gss can occur even though @xmath22 .     when @xmath102 and @xmath109 : ( a ) @xmath283 , @xmath307 , @xmath310 ; ( b ) @xmath283 , @xmath311 , @xmath101 ; ( c ) @xmath104 , @xmath312 , @xmath101 .
the red and blue curves represent the stable and unstable manifolds , respectively . in plate  (
a ) the invariant manifolds have quadratic tangencies inside the circles.,title="fig : " ]   when @xmath102 and @xmath109 : ( a ) @xmath283 , @xmath307 , @xmath310 ; ( b ) @xmath283 , @xmath311 , @xmath101 ; ( c ) @xmath104 , @xmath312 , @xmath101 .
the red and blue curves represent the stable and unstable manifolds , respectively . in plate  ( a )
the invariant manifolds have quadratic tangencies inside the circles.,title="fig : " ] +   when @xmath102 and @xmath109 : ( a ) @xmath283 , @xmath307 , @xmath310 ; ( b ) @xmath283 , @xmath311 , @xmath101 ; ( c ) @xmath104 , @xmath312 , @xmath101 .
the red and blue curves represent the stable and unstable manifolds , respectively . in plate
( a ) the invariant manifolds have quadratic tangencies inside the circles.,title="fig : " ]    we illustrate these general ideas with computations for equation .
figure  [ f : im ] shows the stable and unstable manifolds of the poincar map of on the section @xmath182 for the parameter values when saddle - node or pitchfork bifurcations occur in figs .
[ f : bd ] .    figure  [ f : im](a ) corresponds to fig .
[ f : hbmap](a ) and depicts the situation when the map is not reversible ( @xmath178 ) and the stable and unstable manifolds have a quadratic tangency .
similarly , in figs .
[ f : im](b ) and ( c ) the poincar map is reversible , and the manifolds have cubic and quadratic tangencies on the invariant plane @xmath313 as in the sketches in figs .
[ f : hbmap](b ) and ( c ) .
of ( [ e : uvans ] ) born at the bifurcations of fig .
[ f : bd](a)-(c ) for @xmath287 , @xmath283 , @xmath102 and @xmath101 .
the gss in plates  ( a ) and ( b ) are symmetric and born at the saddle - node bifurcation while the gss in plate  ( c ) is asymmetric and born at the pitchfork bifurcation . ]     of ( [ e : uvans ] ) born at the bifurcations of fig .
[ f : bd](a)-(c ) for @xmath287 , @xmath283 , @xmath102 and @xmath101 .
the gss in plates  ( a ) and ( b ) are symmetric and born at the saddle - node bifurcation while the gss in plate  ( c ) is asymmetric and born at the pitchfork bifurcation . ]     of ( [ e : uvans ] ) born at the bifurcations of fig .
[ f : bd](a)-(c ) for @xmath287 , @xmath283 , @xmath102 and @xmath101 .
the gss in plates  ( a ) and ( b ) are symmetric and born at the saddle - node bifurcation while the gss in plate  ( c ) is asymmetric and born at the pitchfork bifurcation . ]     of ( [ e : uvans ] ) born at the saddle - node bifurcation of fig .
[ f : bd](d ) for @xmath105 , @xmath104 , @xmath101 and @xmath102 . ]     of ( [ e : uvans ] ) born at the saddle - node bifurcation of fig .  [ f : bd](d ) for @xmath105 , @xmath104 , @xmath101 and @xmath102 . ]
figures  [ f : newgs-1 ] and [ f : newgs0 ] show the profiles of gss born at the bifurcations in figs .
[ f : bd](a)-(c ) and ( d ) , respectively . in particular , the gss in fig .
[ f : newgs0 ] are bound states with two humps .
standard results about transverse homoclinic orbits explain the existence of multi - pulse gss , see @xcite , and we also find such gss in the other gap regions .
we note that bifurcations of symmetric homoclinic orbits in reversible systems that are _ autonomous _ have been discussed in buchto:96,kn:97 .
our analysis of the poincar map also yields comparable bifurcation results in systems with periodic forcing like eqs .
( [ e : ab ] ) .
in this paper we have studied the existence of gss ( gap solitons ) in a model of the bg ( bragg grating ) subjected to periodic modulation along the fiber , thus creating the supergrating .
the corresponding mathematical model amounts to a parametrically forced extension of the standard bg model with the additional symmetry - breaking forcing term .
as was known previously @xcite , the supergrating opens up a set of new bandgaps . using the melnikov method and averaging techniques
, we have demonstrated that families of gss exist robustly in the new gaps .
we have also analyzed bifurcations of gss and showed that , in addition to the fundamental gss , the model also supports families of bound states of the solitons .
many of the gss found in the model are completely stable , as was verified by means of direct simulations for some of them @xcite .
the present paper extends the results of @xcite and gives detailed mathematical proofs for them , although we did not analyze the stability problem here in detail .
in particular , we have discussed the effects of the symmetry - breaking forcing term and demonstrated that it can lead to the closing and re - opening of the bandgaps . in the presence of the new term ,
gss exist robustly ( as generic solutions ) in all the bandgaps , but , as the symmetry of the equation is broken , they undergo many bifurcations .
this research was initiated and many important parts were carried out when the authors worked with alan champneys or visited him at the university of bristol .
we thank him for discussions and comments as well as his hospitality and support .
k.y . acknowledges support from the japan society for the promotion of science , grant - in - aid for scientific research ( c ) nos .  18560056 , 21540124 and 22540180 .
]    let @xmath314^j}\ ] ] and let @xmath315 for @xmath316 . to estimate ( [ e : ij ] ) , we consider the complex integral @xmath317 where @xmath318 is a rectangle contour , @xmath319 and @xmath320 , depicted in fig .
[ f : a1 ] , with @xmath321 .
the complex function @xmath322 is singular only at @xmath323 which are @xmath124th - order poles , inside of @xmath318 . by the residue theorem
, we have @xmath324 where @xmath325 is the residue of @xmath326 at @xmath327 .
we compute @xmath328 so that @xmath329 and @xmath330.\end{aligned}\ ] ]    on the other hand , we can write @xmath331 since by @xmath332 @xmath333 and the second and fourth integrals in ( [ e : int3 ] ) tend to zero as @xmath334 , we have @xmath335 where we used ( [ e : fj ] ) .
hence , we obtain @xmath336 and @xmath337,\ ] ] from which eqs .
( [ e : int1 ] ) and ( [ e : int2 ] ) follow .
applying the third - order averaging method , we calculate the boundaries of the central gap . again
, the ` mathematica ` program `` haverage.m '' is used to obtain the third - order averaged systems .
let @xmath338 and @xmath339 . using the transformation @xmath340 , we rewrite ( [ e : ab ] ) as @xmath341\eta+o(\varepsilon^{5/2 } ) , \\
\eta^{\prime}= & \varepsilon^{1/2}[\cos kx+\bar{\mu}\cos(kx+\delta)]\xi + \varepsilon^{3/2}\left(-\nu\xi+\frac{3}{2}\xi^{3}\right)+o(\varepsilon^{5/2 } ) .
\end{split } \label{e : abr}\ ] ] the third - order averaged system for ( [ e : abr ] ) is obtained as @xmath342 ,   \label{e : avr}\ ] ] where @xmath343 was taken as the small parameter . in the averaged system ( [ e : avr ] ) , the origin is a saddle if @xmath344 hence , we have the approximate right boundary of the central gap , @xmath345      let @xmath346 and @xmath339 . using the transformation @xmath347 , we rewrite ( [ e : ab ] ) as @xmath348\eta -\varepsilon^{3/2}\left(\nu\eta+\frac{3}{2}\eta^{3}\right)+o(\varepsilon^{5/2 } ) , \\
\eta^{\prime}= & -2\varepsilon^{1/2}\xi + \varepsilon^{3/2}[\cos kx+\bar{\mu}\cos(kx+\delta)]\xi+o(\varepsilon^{5/2 } ) , \end{split } \label{e : abl}\ ] ] the third - order averaged system for ( [ e : abl ] ) is obtained as @xmath349,\quad \eta^{\prime}=-2\varepsilon^{1/2}\xi,\ ] ] in which the origin is a saddle if @xmath350 thus , we have the approximate left boundary of the central gap , @xmath351                        e.  j. doedel , a.  r. champneys , t.  r. fairgrieve , yu .
a. kuznetsov , b.  sandstede , and x.  j. wang , continuation and bifurcation software for ordinary differential equations , 1997 ( an upgraded version is available at ` http://cmvl.cs.concordia.ca/auto/ ` ) .
t. dohnal , d.  e. pelinovsky and g. schneider , coupled - mode equations and gap solitons in a two - dimensional nonlinear elliptic problem with a separable periodic potential , _ j. nonlinear sci .
_ 19 ( 2009 ) 95131 .                                                          c.  sim , averaging under fast quasiperiodic forcing , in : j.  seimenis ( ed . ) , _ hamiltonian mechanics : integrability and chaotic behavior , proceedings of a nato advanced research workshop held in toru , june 28-july 2 , 1993 _ , plenum , new york , 1994 , pp .