it is believed that the direct detection of gravitational waves ( gws ) will bring the era of gravitational wave astronomy .
the interferometer detectors are now under operation and awaiting the first signal of gws  @xcite .
it is also known that pulsar timing arrays ( ptas ) can be used as a detector for gws @xcite .
these detectors are used to search for very low frequency ( @xmath0 ) gravitational waves , where the lower limit of the observable frequencies is determined by the inverse of total observation time @xmath1 .
indeed , the total observation time has a crucial role in ptas , because ptas are most sensitive near the lower edge of observable frequencies @xcite . taking into account its sensitivity ,
the first direct detection of the gravitational waves might be achieved by ptas .
the main target of ptas is the stochastic gravitational wave background ( sgwb ) generated by a large number of unresolved sources with the astrophysical origin or the cosmological origin in the early universe .
the promising sources are super massive black hole binaries  @xcite , cosmic ( super)string  @xcite , and inflation  @xcite .
previous studies have assumed that the sgwb is isotropic and unpolarized  @xcite .
these assumptions are reasonable for the primary detection of the sgwb , but the deviation from the isotropy and the polarizations should have rich information of sources of gravitational waves .
recently , the cross - correlation formalism has been generalized to deal with anisotropy in the sgwb @xcite .
result of this work enables us to consider arbitrary levels of anisotropy , and a bayesian approach was performed by using this formalism @xcite . on the other hand , for the anisotropy of the sgwb , the cross - correlation formalism has been also developed in the case of interferometer detectors  @xcite .
as to the polarization , there are works including the ones motivated by the modified gravity  @xcite
. we can envisage supermassive black hole binaries emit circularly polarized sgwb due to the chern - simons term  @xcite .
there may also exist cosmological sgwb with circular polarization in the presence of parity violating term in gravity sector  @xcite .    in this paper
, we investigate the detectability of circular polarization in the sgwb by ptas .
we characterize sgwb by the so called stokes @xmath2 parameter  @xcite and calculate generalized overlap reduction functions ( orfs ) so that we can probe the circular polarization of the sgwb .
we also discuss a method to separate the intensity ( @xmath3 mode ) and circular polarization ( @xmath2 mode ) of the sgwb .
the paper is organized as follows . in section [ sec :
stokes parameters for a plane gravitational wave ] , we introduce the stokes parameters for monochromatic plane gravitational waves , and clarify the physical meaning of the stokes parameters @xmath3 and @xmath2 . in section [ sec : formulation ] , we formulate the cross - correlation formalism for anisotropic circularly polarized sgwb with ptas .
the basic framework is essentially a combination of the formalism of @xcite , and the polarization decomposition formula of the sgwb derived in @xcite . in section [ sec : the generalized overlap reduction function for circular polarization ] , we calculate the generalized orfs for the @xmath2 mode .
the results for @xmath3 mode are consistent with the previous work  @xcite . in section [ sec : separation method ] , we give a method for separation between the @xmath3 mode and @xmath2 mode of the sgwb .
the final section is devoted to the conclusion . in appendixes , we present analytic results for the generalized overlap reduction functions . in this paper
, we will use the gravitational units @xmath4 .
let us consider the stokes parameters for plane waves traveling in the direction @xmath5 , which can be described by @xmath6 \
, \\ & & h_{xy}(t , z)=h_{yx}(t , z)={\rm re}[b_{\times}\mathrm{e}^{-iw(t - z ) } ] \ .\end{aligned}\ ] ] for an idealized monochromatic plane wave , complex amplitudes @xmath7 and @xmath8 are constants .
polarization of the plane gws is characterized by the tensor , ( see @xcite and also electromagnetic case @xcite ) @xmath9 where @xmath10 take @xmath11 .
any @xmath12 hermitian matrix can be expanded by the pauli and the unit matrices with real coefficients .
hence , the @xmath13 hermitian matrix @xmath14 can be written as @xmath15 where @xmath16 by analogy with electromagnetic cases , @xmath17 and @xmath2 are called stokes parameters . comparing with , we can read off the stokes parameters as @xmath18= b_{+}^{\ast}b_{\times}+ b_{\times}^{\ast}b_{+},\\ v&=&-2{\rm i m } [ b_{+}^{\ast}b_{\times}]=i ( b_{+}^{\ast}b_{\times}- b_{\times}^{\ast}b_{+}).\label{stv}\end{aligned}\ ] ] apparently , the real parameter @xmath3 is the intensity of gws . in order to reveal the physical meaning of the real parameter @xmath2 , we define the circular polarization bases @xcite @xmath19 from the relation @xmath20 we see @xmath21
thus , we can rewrite the stokes parameters - as @xmath22 from the above expression , we see that the real parameter @xmath2 characterizes the asymmetry of circular polarization amplitudes .
the other parameters @xmath23 and @xmath24 have additional information about linear polarizations by analogy with the electromagnetic cases .
alternatively , we can also define the tensor @xmath25 in circular polarization bases @xmath26 where @xmath27 .
note that the stokes parameters satisfy a relation @xmath28    next , we consider the transformation of the stokes parameters under rotations around the @xmath5 axis . the rotation around the @xmath5 axis is given by @xmath29 where @xmath30 is the angle of the rotation .
the gws traveling in the direction @xmath5 @xmath31 transform as @xmath32 where we took the transverse traceless gauge @xmath33 after a short calculation , we obtain @xmath34 using and , the four stokes parameters ( [ sti])-([stv ] ) transform as @xmath35 as you can see , the parameters @xmath23 and @xmath24 depend on the rotation angle @xmath30 .
this reflects the fact that @xmath23 and @xmath24 parameters characterize linear polarizations .
note that this transformation is similar to the transformation of electromagnetic case except for the angle @xmath36 and can be rewritten as @xmath37
in this section , we study anisotropic distribution of sgwb and focus on the detectability of circular polarizations with pulsar timing arrays .
we combine the analysis of @xcite and that of @xcite . in sec.[subsec : the spectral ] , we derive the power spectral density for anisotropic circularly polarized sgwb @xmath38 .
then we also derive the dimensionless density parameter @xmath39 which is expressed by the frequency spectrum of intensity @xmath40  @xcite . in sec.[subsec : the signal ] , we extend the generalized orfs to cases with circular polarizations characterized by the parameter @xmath2 . for simplicity ,
we consider specific anisotropic patterns with @xmath41 expressed by the spherical harmonics @xmath42 .      in the transverse traceless gauge , metric perturbations @xmath43 with a given propagation direction @xmath44
can be expanded as @xcite @xmath45 where the fourier amplitude satisfies @xmath46 as a consequence of the reality of @xmath43 , @xmath47 , @xmath48 is the frequency of the gws , @xmath49 are spatial indices , @xmath50 label polarizations .
note that the fourier amplitude @xmath51 satisfies the relation @xmath52 where @xmath53 was defined by .
the polarized tensors @xmath54 are defined by @xmath55 where @xmath56 and @xmath57 are unit orthogonal vectors perpendicular to @xmath58 .
the polarization tensors satisfy @xmath59 with polar coordinates , the direction @xmath44 can be represented by @xmath60 and the polarization basis vectors read @xmath61    we assume the fourier amplitudes @xmath62 are random variables , which is stationary and gaussian .
however , they are not isotropic and unpolarized .
the ensemble average of fourier amplitudes can be written as @xcite @xmath63 where @xmath64 here , the bracket @xmath65 represents an ensemble average , and @xmath66 is the dirac delta function on the two - sphere .
the gw power spectral density @xmath38 is a hermitian matrix , and satisfies @xmath67 because of the relation @xmath46 .
therefore , we have the relations @xmath68 note that the stokes parameters are not exactly the same as the expression of , but they have the relation and characterize the same polarization .
we further assume that the sgwbs satisfy @xmath69 we also assume the directional dependence of the sgwb is frequency independent @xcite .
this implies the gw power spectral density is factorized into two parts , one of which depends on the direction while the other depends on the frequency .
because of the transformations - , the parameters @xmath3 and @xmath2 have spin 0 and the parameters @xmath70 have spin @xmath71  @xcite . to analyze the sgwb on the sky , it is convenient to expand the stokes parameters by spherical harmonics @xmath72
. however , since @xmath70 parameters have spin @xmath71 , they have to be expanded by the spin - weighted harmonics @xmath73 @xcite .
thus , we obtain @xmath74 in this paper , we study specific anisotropic patterns with @xmath41 for simplicity .
therefore , we can neglect @xmath23 and @xmath24 from now on .
thus , the gw power spectral density becomes @xmath75 where @xmath76 so , we focus on the parameters @xmath3 and @xmath2 . in what follows , we will use the following shorthand notation @xmath77    next , we consider the dimensionless density parameter  @xcite @xmath78 where @xmath79 is the critical density , @xmath80 is the present value of the hubble parameter , @xmath81 is the energy density of gravitational waves , and @xmath82 is the energy density in the frequency range @xmath48 to @xmath83 .
the bracket @xmath65 represents the ensemble average .
however , actually , we take a spatial average over the wave lengths @xmath84 of gws or a temporal average over the periods @xmath85 of gws . here
, we assumed the ergodicity , namely , the ensemble average can be replaced by the temporal average .
using , , , as well as @xmath46 and @xmath86 , we get @xmath87 then we define @xmath88 hence , the dimensionless quantity @xmath39 in is given by @xmath89 where the spherical harmonics are orthogonal and normalized as @xmath90 using @xmath91 , we obtain @xmath92 without loss of generality , we normalize the monopole moment as @xmath93 so , becomes @xmath94      the time of arrival of radio pulses from the pulsar is affected by gws .
consider a pulsar with frequency @xmath95 located in the direction @xmath96 . to detect the sgwb ,
let us consider the redshift of the pulse from a pulsar @xcite @xmath97 where @xmath98 is a frequency detected at the earth and @xmath96 is the direction to the pulsar .
the unit vector @xmath44 represents the direction of propagation of gravitational plane waves .
we also defined the difference between the metric perturbations at the pulsar @xmath99 and at the earth @xmath100 as @xmath101 the gravitational plane waves at each point is defined as @xmath102 for the sgwb , the redshift have to be integrated over the direction of propagation of the gravitational waves @xmath44 : @xmath103 we choose a coordinate system @xmath104 and assume that the amplitudes of the metric perturbation at the pulsar and the earth are the same .
then becomes @xmath105 and therefore , reads @xmath106 where we have defined the pattern functions for pulsars @xmath107 note that our convention for the fourier transformation is @xmath108 therefore , the fourier transformation of can be written as @xmath109    in the actual signals from a pulsar , there exist noises .
hence , we need to use the correlation analysis .
we consider the signals from two pulsars @xmath110 where @xmath111 labels the pulsar . here
, @xmath112 denotes the signal from the pulsar and @xmath113 denotes the noise intrinsic to the measurement .
we assume the noises are stationary , gaussian and are not correlated between the two pulsars .
to correlate the signals of two measurements , we define @xmath114 where @xmath1 is the total observation time and @xmath115 is a real filter function which should be optimal to maximize signal - to - noise ratio . in the case of interferometer
, the optimal filter function falls to zero for large @xmath116 compered to the travel time of the light between the detecters .
since the signals of two detectors are expected to correlate due to the same effect of the gravitational waves , the optimal filter function should behave this way .
then , typically one of the detectors is very close to the other compared to the total observation time @xmath1 .
therefore , the total observation time @xmath1 can be extended to @xmath117 @xcite .
in contrast , in the case of pta , it is invalid that @xmath1 is very large compered to the travel time of the light between the pulsars .
nevertheless , we can assume that one of the two @xmath1 can be expanded to @xmath117 , because in situations @xmath118 and @xmath119 it is known that we can ignore the effect of the distance @xmath120 of pulsars .
in this case , it is clear that any locations of the pulsars are optimal and optimal filter function should behave like as the interferometer case @xcite .    using these assumptions @xmath118 and @xmath119 , we can rewrite as @xmath121 where @xmath122 note that @xmath123 satisfies @xmath124 , because @xmath125 is real .
moreover , to deal with the unphysical region @xmath126 we require @xmath127 .
thus , @xmath123 becomes real .
taking the ensemble average , using @xmath128 , @xmath118 , and assuming the noises in the two measurements are not correlated , we get @xmath129\ , \label{s2}\end{aligned}\ ] ] where we have defined @xmath130 the functions @xmath131 and @xmath132 are called the generalized orfs , which describe the angular sensitivity of the pulsars for the sgwb . note that , as we already mentioned , we consider the cases of @xmath41 for simplicity
. then we have assumed @xmath118 and @xmath119 , this assumption implies that approximately becomes @xmath133 due to the rapid oscillation of the phase factor .
therefore , the distance @xmath120 of the pulsars does not appear in the generalized orfs , and hence the generalized orfs do not depend on the frequency .    as you can see from ,
the correlation of the two measurements involve both the total intensity and the circular polarization .
however , the degeneracy can be disentangled by using separation method , which will be discussed in the section [ sec : separation method ] .
in this section , we consider the generalized orfs for circular polarizations : @xmath134 where we defined @xmath135 in the above , we have used and the fact that the generalized orfs do not depend on frequency . for computation of the generalized orfs for circular polarizations ,
it is convenient to use the computational frame @xcite defined by @xmath136 where @xmath137 is the angular separation between the two pulsars . using - , , and
, one can easily show that @xmath138 we therefore get @xmath139 the explicit form of the spherical harmonics reads @xmath140 where @xmath141 is the normalization factor .
the associated legendre functions are given by @xmath142 and @xmath143 with the legendre functions @xmath144\ .\label{pl}\end{aligned}\ ] ] using the spherical harmonics , becomes @xmath145 where we have used the fact that the function of @xmath146 is odd parity in the case of @xmath147 and is even parity in the case of @xmath148 .
note that the generalized orfs for circular polarizations are real functions . in the case of @xmath149 and/or @xmath150 ,
the integrand in vanishes .
therefore , we can not detect circular polarizations for these cases .
this fact for @xmath151 implies that we do not need to consider auto - correlation for a single pulsar .
this is the reason why we neglected auto - correlation term in .    integrating ,
we get the following form for @xmath152 : @xmath153 for @xmath154 , we have obtained @xmath155 \ , \\ \gamma^{v}_{1 - 1}&=&\gamma^{v}_{11 } \ , \end{aligned}\ ] ] recall that @xmath156 .
the derivation of this formula for @xmath154 can be found in appendix [ sec : angular integral of the generalized overlap reduction function for dipole circular polarization ] .    for @xmath157 , we derived the following : @xmath158\ , \\ \gamma^{v}_{2 - 1}&=&\gamma^{v}_{21}\ , \\
\gamma^{v}_{22}&=&-\frac{\sqrt{30\pi}}{6}(1-\cos\xi)\left[2-\cos\xi+6\left(\frac{1-\cos\xi}{1+\cos\xi}\right)\log\left(\sin\frac{\xi}{2}\right)\right]\ , \\ \gamma^{v}_{2 - 2}&=&-\gamma^{v}_{22}\ , \end{aligned}\ ] ] for @xmath159 , the results are @xmath160\ , \\ \gamma^{v}_{3 - 1}&=&\gamma^{v}_{31}\ , \\ \gamma^{v}_{32}&=&\frac{\sqrt{210\pi}}{24}(1-\cos\xi)\left[8 - 5\cos\xi-\cos^2\xi+24\left(\frac{1-\cos\xi}{1+\cos\xi}\right)\log\left(\sin\frac{\xi}{2}\right)\right]\ , \\ \gamma^{v}_{3 - 2}&=&-\gamma^{v}_{3 - 2}\ , \\ \gamma^{v}_{33}&=&-\frac{\sqrt{35\pi}}{16}\sin\xi\left(\frac{1-\cos\xi}{1+\cos\xi}\right)\left[11 - 6\cos\xi-\cos^2\xi+32\left(\frac{1-\cos\xi}{1+\cos\xi}\right)\log\left(\sin\frac{\xi}{2}\right)\right]\ , \\ \gamma^{v}_{3 - 3}&=&\gamma^{v}_{33}\ .\end{aligned}\ ] ] in fig .
[ gv ] , we plotted these generalized orfs as a function of the angular separation between the two pulsars @xmath137 .
it is apparent that considering the @xmath2 mode does not make sense when we only consider the isotropic ( @xmath152 ) orf . on the other hand ,
when we consider anisotropic ( @xmath161 ) orfs , it is worth taking into account polarizations .
the polarizations of the sgwb would give us rich information both of super massive black hole binaries and of inflation in the early universe .
as a function of the angular separation between the two pulsars @xmath137 . in fig .
[ gv](a ) , we find the orf for the monopole ( l=0 ) is trivial . in fig .
[ gv](b ) , the orfs for the dipole ( l=1 ) are shown . in fig .
[ gv](c ) , the orfs for the quadrupole ( l=2 ) are depicted . in fig .
[ gv](d ) , the orfs for the octupole ( l=3 ) are plotted .
the black solid curve , the blue dashed curve , the red dotted curve , the dark - red space - dotted curve , and the green long - dashed curve represent @xmath149 , @xmath162 , @xmath163 , @xmath164 , @xmath165 , respectively.,title="fig:",width=340 ] ( a ) @xmath152     as a function of the angular separation between the two pulsars @xmath137 . in fig . [ gv](a ) , we find the orf for the monopole ( l=0 ) is trivial . in fig .
[ gv](b ) , the orfs for the dipole ( l=1 ) are shown . in fig .
[ gv](c ) , the orfs for the quadrupole ( l=2 ) are depicted . in fig .
[ gv](d ) , the orfs for the octupole ( l=3 ) are plotted .
the black solid curve , the blue dashed curve , the red dotted curve , the dark - red space - dotted curve , and the green long - dashed curve represent @xmath149 , @xmath162 , @xmath163 , @xmath164 , @xmath165 , respectively.,title="fig:",width=340 ] ( b ) @xmath154     as a function of the angular separation between the two pulsars @xmath137 . in fig . [ gv](a ) , we find the orf for the monopole ( l=0 ) is trivial . in fig .
[ gv](b ) , the orfs for the dipole ( l=1 ) are shown . in fig .
[ gv](c ) , the orfs for the quadrupole ( l=2 ) are depicted . in fig .
[ gv](d ) , the orfs for the octupole ( l=3 ) are plotted .
the black solid curve , the blue dashed curve , the red dotted curve , the dark - red space - dotted curve , and the green long - dashed curve represent @xmath149 , @xmath162 , @xmath163 , @xmath164 , @xmath165 , respectively.,title="fig:",width=340 ] ( c ) @xmath157     as a function of the angular separation between the two pulsars @xmath137 . in fig .
[ gv](a ) , we find the orf for the monopole ( l=0 ) is trivial . in fig .
[ gv](b ) , the orfs for the dipole ( l=1 ) are shown . in fig .
[ gv](c ) , the orfs for the quadrupole ( l=2 ) are depicted . in fig .
[ gv](d ) , the orfs for the octupole ( l=3 ) are plotted .
the black solid curve , the blue dashed curve , the red dotted curve , the dark - red space - dotted curve , and the green long - dashed curve represent @xmath149 , @xmath162 , @xmath163 , @xmath164 , @xmath165 , respectively.,title="fig:",width=340 ] ( d ) @xmath159    using the same procedure described in the above to derive the generalized orfs for circular polarizations , we can also derive the generalized orfs for the intensity @xmath166 where @xmath167 the angular integral in this case was performed in @xcite .
the results are summarized in appendix [ sec : the generalized overlap reduction function for intensity ] .
in this section , we separate the @xmath3 mode and @xmath2 mode of the sgwb with correlation analysis @xcite . to this aim ,
we use four pulsars ( actually we need at least three pulsars ) , and define correlations of @xmath168 @xmath169 where @xmath170 label the pulsars . comparing with , we obtain @xmath171 \ ,
\label{1c12}\\ & & c_{34}(f)=\sum_{lm}^{l=3}\left[c_{lm}^{i}i(f)\gamma_{lm,34}^{i}+c_{lm}^{v}v(f)\gamma_{lm,34}^{v}\right ] \ .\label{1c34}\end{aligned}\ ] ] if the @xmath3 mode and @xmath2 mode of the sgwb are dominated by a certain @xmath172 and @xmath173 , and become @xmath174 \ , \label{2c12 } \\ & & c_{34}(f)=\left[c _ { l m}^{i}i(f)\gamma _ { l m,34}^{i}+c _ { l ' m'}^{v}v(f)\gamma _ { l ' m',34}^{v}\right ] \ .\label{2c34}\end{aligned}\ ] ] to separate the intensity and the circular polarization , we take the following linear combinations @xmath175 where we defined coefficients @xmath176 as you can see , @xmath177 contains only @xmath40 , and @xmath178 contains only @xmath179 .
for the signal @xmath180 , the formulas corresponding to and are given by @xmath181 \ , \label{sp}\end{aligned}\ ] ] where @xmath182 denotes @xmath3 and @xmath2 .
we assume @xmath183 and that the noise in the four pulsars are not correlated .
we also assume that the ensemble average of fourier amplitudes of the noises @xmath184 is of the form @xmath185 where @xmath186 is the noise power spectral density .
the reality of @xmath187 gives rise to @xmath188 and therefore we obtain @xmath189 . without loss of generality
, we can assume @xmath190 then we obtain corresponding noises @xmath191 : @xmath192\ , \label{np}\end{aligned}\ ] ] where @xmath193^{1/2 } \label{sn12 } \ , \quad s_{n,34}(f ) \equiv [ s_{n,3}(f)s_{n,4}(f)]^{1/2 } \label{sn34 } \ .\end{aligned}\ ] ] using the inner product @xmath194 \ , \end{aligned}\ ] ] we can rewrite , as @xmath195 therefore , the optimal filter function can be chosen as @xmath196 using , we get optimal signal - to - noise ratio @xmath197^{1/2}\ .\label{snr}\end{aligned}\ ] ] plugging , , and into , we obtain @xmath198^{1/2}\ , \\
{ \rm snr}_{v}&=&\left[t\int_{-\infty}^{\infty}df\,\,\frac{\left(c^{v}_{{l}'{m}'}\right)^{2}v^{2}(f)\left(\gamma_{{l}'{m}',34}^{v}\gamma^{i}_{{l}{m},12}-\gamma_{{l}'{m}',12}^{v}\gamma^{i}_{{l}{m},34}\right)^2}{\left(\gamma^{i}_{{l}{m},12}\right)^2s^{2}_{n,34}(f)+\left(\gamma^{i}_{{l}{m},34}\right)^2s^{2}_{n,12}(f)}\right]^{1/2}\ .\end{aligned}\ ] ] if we assume all of the noise power spectral densities are the same , becomes @xmath199 thus , the compiled orfs can be defined as @xmath200^{1/2}}\ , \\ \gamma_{12:34}^{v}&\equiv&\frac{\gamma_{{l}'{m}',34}^{v}\gamma^{i}_{{l}{m},12}-\gamma_{{l}'{m}',12}^{v}\gamma^{i}_{{l}{m},34}}{\left[\left(\gamma^{i}_{{l}{m},12}\right)^2+\left(\gamma^{i}_{{l}{m},34}\right)^2\right]^{1/2}}\ .\end{aligned}\ ] ] this compiled orfs @xmath201 and @xmath202 describe the angular sensitivity of the four pulsars for the pure @xmath3 and @xmath2 mode of the sgwb , respectively .
note that , to do this separation , we must know a priori the coefficients @xmath203 and @xmath204 .
if we do not assume , the generalized orfs depend on the frequency . in this case
, it seems difficult to calculate these coefficients .
we next consider the case that @xmath3 mode and/or @xmath2 mode dominant in two or more @xmath205 . in this case , if we have a priori knowledge of the values of @xmath206 in each of @xmath205 for coefficients
@xmath203 and @xmath204 , we can separate @xmath3 mode and @xmath2 mode . for example , assume that @xmath3 mode is dominated by @xmath207 , while @xmath2 mode is dominated by @xmath208 , then and become @xmath209\ , \label{3c12}\\ & & c_{34}(f)=\left[c^{i}_{00}i(f)\left(\gamma_{00,34}^{i}+\frac{c^{i}_{11}}{c^{i}_{00}}\gamma_{11,34}^{i}\right)+c_{11}^{v}v(f)\gamma_{11,34}^{v}\right]\ .\label{3c34}\end{aligned}\ ] ] thus , we can separate @xmath3 mode and @xmath2 mode by using linear combinations @xmath210\ , \\
d_{v}&\equiv&a_{v}c_{34}(f)+b_{v}c_{12}(f ) \nonumber\\ & = & c_{11}^{v}v(f)\left[\gamma_{11,34}^{v}\left(\gamma_{00,12}^{i}+\frac{c^{i}_{11}}{c^{i}_{00}}\gamma_{11,12}^{i}\right)-\gamma_{11,12}^{v}\left(\gamma_{00,34}^{i}+\frac{c^{i}_{11}}{c^{i}_{00}}\gamma_{11,34}^{i}\right)\right]\ , \end{aligned}\ ] ] where @xmath211 as in the previous calculations , we can get the compiled orfs @xmath212^{1/2}}\ , \label{gi1234}\\ \gamma_{12:34}^{v}&\equiv&\frac{\gamma_{11,34}^{v}\left(\gamma_{00,12}^{i}+\displaystyle\frac{c^{i}_{11}}{c^{i}_{00}}\gamma_{11,12}^{i}\right)-\gamma_{11,12}^{v}\left(\gamma_{00,34}^{i}+\displaystyle\frac{c^{i}_{11}}{c^{i}_{00}}\gamma_{11,34}^{i}\right)}{\left[\left(\gamma_{00,12}^{i}+\displaystyle\frac{c^{i}_{11}}{c^{i}_{00}}\gamma_{11,12}^{i}\right)^2+\left(\gamma_{00,34}^{i}+\displaystyle\frac{c^{i}_{11}}{c^{i}_{00}}\gamma_{11,34}^{i}\right)^2\right]^{1/2}}\ .\label{gv1234}\end{aligned}\ ] ]    [ cols="^,^ " , ]     in fig .
[ cg ] we show some compiled orfs @xmath213 ( left panels ) and @xmath214 ( right panels ) as a function of the two angular separations @xmath137 and @xmath215 for two pulsar pairs , respectively .
we used the expressions of @xmath2 mode and @xmath3 mode ( see appendix [ sec : the generalized overlap reduction function for intensity ] ) , and we assumed @xmath216 for simplicity . in fig .
[ cg](a ) and [ cg](b ) , the @xmath3 mode is dominated by @xmath217 and @xmath2 mode is dominated by @xmath218 . in fig .
[ cg](c ) and [ cg](d ) , the @xmath3 mode is dominated by @xmath219 and @xmath2 mode is dominated by @xmath218 . in fig .
[ cg](e ) and [ cg](f ) , the @xmath3 mode is dominated by @xmath207 and @xmath2 mode is dominated by @xmath218 . in fig .
[ cg](e ) and [ cg](f ) , the @xmath3 mode is dominated by @xmath220 and @xmath2 mode is dominated by @xmath218 . by definition , in the case of @xmath221 ,
the compiled orfs are zero .
we have studied the detectability of the stochastic gravitational waves with ptas . in most of the previous works ,
the isotropy of sgwb has been assumed for the analysis .
recently , however , a stochastic gravitational wave background with anisotropy have been considered .
the information of the anisotropic pattern of the distribution should contain important information of the sources such as supermassive black hole binaries and the sources in the early universe .
it is also intriguing to take into account the polarization of sgwb in the pta analysis .
therefore , we extended the correlation analysis to circularly polarized sgwb and calculated generalized overlap reduction functions for them .
it turned out that the circular polarization can not be detected for an isotropic background .
however , when the distribution has anisotropy , we have shown that there is a chance to observe circular polarizations in the sgwb .
we also discussed how to separate polarized modes from unpolarized modes of gravitational waves .
if we have a priori knowledge of the abundance ratio for each mode in each of @xmath205 , we can separate @xmath3 mode and @xmath2 mode in general .
this would be possible if we start from fundamental theory and calculate the spectrum of sgwb .
in particular , in the case that the signal of lowest @xmath222 is dominant , we performed the separation of @xmath3 mode and @xmath2 mode explicitly .
this work was supported by grants - in - aid for scientific research ( c ) no.25400251 and " mext grant - in - aid for scientific research on innovative areas no.26104708 and `` cosmic acceleration''(no.15h05895 ) .
in this appendix , we perform angular integration of the generalized orf for dipole ( @xmath154 ) circular polarization ( see @xcite ) : @xmath223 where we have defined @xmath224 .
it is obvious that in the case of @xmath225 , integrand of the generalized orf is zero , because of @xmath226 , then we obtain @xmath227 then , using - , we calculate @xmath228 and we find @xmath229    therefore we only have to consider the dipole generalized orf in the case of @xmath154 , @xmath230 : @xmath231 where @xmath232 first , to calculate @xmath233 , we use contour integral in the complex plane . defining @xmath234 and substituting @xmath235 into , we can rewrite @xmath233 as @xmath236 } \ , \end{aligned}\ ] ] where @xmath237 denotes a unit circle .
we can factorize the denominator of the integrand and get @xmath238 where @xmath239 hereafter , the upper sign applies when @xmath240 and the lower one applies when @xmath241 .
note that we only consider the region @xmath242 , so we have used the relation @xmath243 in above expression . in the region
@xmath244 , @xmath245 is inside the unit circle @xmath237 except for @xmath246 and @xmath247 is outside the unit circle @xmath237 .
now , we can perform the integral using the residue theorem @xmath248 where @xmath249 the residues inside the unit circle @xmath237 can be evaluated as @xmath250\right\ }             = \frac{i(z_{+}+z_{-})}{2\sqrt{1-x^2}\sin\xi } \ , \end{aligned}\ ] ] @xmath251 thus , we obtain @xmath252 next , we consider @xmath253 defined in
. using , we can calculate @xmath253 as @xmath254    similarly , we can evaluate @xmath255 given in . to calculate @xmath255 in the complex plane , we again substitute into and obtain @xmath256 we use the residue theorem @xmath257 where @xmath258 the residues inside the unit circle @xmath237 can be calculated as @xmath259\right\ }        = \frac{i(z_{+}^2+z_{-}^2)}{4\sqrt{1-x^2}\sin\xi } \ , \end{aligned}\ ] ] @xmath260 therefore , @xmath255 becomes @xmath261 substituting to , we can calculate @xmath262 : @xmath263    finally , substituting and into , we get the generalized orf for @xmath264 @xmath265\ .\end{aligned}\ ] ]     as a function of the angular separation between the two pulsars @xmath137 .
[ gi](a ) shows monopole ( l=0 ) , fig .
[ gi](b ) shows dipole ( l=1 ) , fig .
[ gi](c ) shows quadrupole ( l=2 ) and fig .
[ gi](d ) shows octupole ( l=3 ) .
the black solid curve , the blue dashed curve , the dark - blue dash - dotted curve , the red dotted curve , the green long - dashed curve , the dark - green space - dashed curve represent @xmath149 , @xmath266 , @xmath267 , @xmath268 , @xmath269 , @xmath270 , respectively.,title="fig:",width=340 ] ( a ) @xmath152     as a function of the angular separation between the two pulsars @xmath137 .
[ gi](a ) shows monopole ( l=0 ) , fig .
[ gi](b ) shows dipole ( l=1 ) , fig .
[ gi](c ) shows quadrupole ( l=2 ) and fig .
[ gi](d ) shows octupole ( l=3 ) .
the black solid curve , the blue dashed curve , the dark - blue dash - dotted curve , the red dotted curve , the green long - dashed curve , the dark - green space - dashed curve represent @xmath149 , @xmath266 , @xmath267 , @xmath268 , @xmath269 , @xmath270 , respectively.,title="fig:",width=340 ] ( b ) @xmath154     as a function of the angular separation between the two pulsars @xmath137 . fig . [ gi](a ) shows monopole ( l=0 ) , fig .
[ gi](b ) shows dipole ( l=1 ) , fig .
[ gi](c ) shows quadrupole ( l=2 ) and fig .
[ gi](d ) shows octupole ( l=3 ) .
the black solid curve , the blue dashed curve , the dark - blue dash - dotted curve , the red dotted curve , the green long - dashed curve , the dark - green space - dashed curve represent @xmath149 , @xmath266 , @xmath267 , @xmath268 , @xmath269 , @xmath270 , respectively.,title="fig:",width=340 ] ( c ) @xmath157     as a function of the angular separation between the two pulsars @xmath137 .
[ gi](a ) shows monopole ( l=0 ) , fig .
[ gi](b ) shows dipole ( l=1 ) , fig .
[ gi](c ) shows quadrupole ( l=2 ) and fig .
[ gi](d ) shows octupole ( l=3 ) . the black solid curve , the blue dashed curve , the dark - blue dash - dotted curve , the red dotted curve , the green long - dashed curve , the dark - green space - dashed curve represent @xmath149 , @xmath266 , @xmath267 , @xmath268 , @xmath269 , @xmath270 , respectively.,title="fig:",width=340 ] ( d ) @xmath159
in this appendix , we show orfs for the intensity @xcite .
the following form for @xmath152 was derived in @xcite , and our expressions are identical to their expressions : @xmath271\ , \end{aligned}\ ] ] for , @xmath154 , we calculated as @xmath272\ , \\ \gamma^{i}_{11}&=&\frac{\sqrt{6\pi}}{12}\sin\xi\left[1 + 3(1-\cos\xi)\left\{1+\frac{4}{1+\cos\xi}\log\left(\sin\frac{\xi}{2}\right)\right\}\right]\ , \\ \gamma^{i}_{1 - 1}&=&-\gamma^{i}_{11}\ , \end{aligned}\ ] ] for @xmath157 , we obtain @xmath273\ , \\ \gamma^{i}_{21}&=&-\frac{\sqrt{30\pi}}{60}\sin\xi\left[21 - 15\cos\xi-5\cos^2\xi+60\left(\frac{1-\cos\xi}{1+\cos\xi}\right)\log\left(\sin\frac{\xi}{2}\right)\right]\ , \\ \gamma^{i}_{2 - 1}&=&-\gamma^{i}_{2 - 1}\ , \\ \gamma^{i}_{22}&=&\frac{\sqrt{30\pi}}{24}(1-\cos\xi)\left[9 - 4\cos\xi-\cos^2\xi+24\left(\frac{1-\cos\xi}{1+\cos\xi}\right)\log\left(\sin\frac{\xi}{2}\right)\right]\ , \\ \gamma^{i}_{2 - 2}&=&\gamma^{i}_{22}\ ,
\end{aligned}\ ] ] for @xmath159 , it is straightforward to reach the following @xmath274\ , \\ \gamma^{i}_{31}&=&\frac{\sqrt{21\pi}}{48}\sin\xi(1-\cos\xi)\left[34 + 15\cos\xi+5\cos^2\xi+\frac{96}{1+\cos\xi}\log\left(\sin\frac{\xi}{2}\right)\right]\ , \\ \gamma^{i}_{3 - 1}&=&-\gamma^{i}_{31}\ , \\ \gamma^{i}_{32}&=&-\frac{\sqrt{210\pi}}{48}(1-\cos\xi)\left[17 - 9\cos\xi-3\cos^2\xi-\cos^3\xi+48\left(\frac{1-\cos\xi}{1+\cos\xi}\right)\log\left(\sin\frac{\xi}{2}\right)\right]\ , \\ \gamma^{i}_{3 - 2}&=&\gamma^{i}_{32}\ , \\ \gamma^{i}_{33}&=&\frac{\sqrt{35\pi}}{48}\frac{(1-\cos\xi)^2}{\sin\xi}\left[34 - 17\cos\xi-4\cos^2\xi-\cos^3\xi+96\left(\frac{1-\cos\xi}{1+\cos\xi}\right)\log\left(\sin\frac{\xi}{2}\right)\right]\ , \\ \gamma^{i}_{3 - 3}&=&-\gamma^{i}_{33}\ .\end{aligned}\ ] ] these are plotted in fig .
the generalized orfs of total intensity are different from that of circular polarization in that the value for @xmath149 is non - trivial .
then the @xmath3 mode orfs for @xmath275 have value even in the case of @xmath151 .
this implies that we can consider auto - correlation for a single pulsar .
99 j.  aasi _ et al .
_ [ ligo scientific collaboration ] , class .
grav .   * 32 * , 074001 ( 2015 ) doi:10.1088/0264 - 9381/32/7/074001 [ arxiv:1411.4547 [ gr - qc ] ] .
f.  acernese _ et al . _
[ virgo collaboration ] , class .
grav .   * 32 * , no . 2 , 024001 ( 2015 ) doi:10.1088/0264 - 9381/32/2/024001 [ arxiv:1408.3978 [ gr - qc ] ] . k.  somiya [ kagra collaboration ] , class .
* 29 * , 124007 ( 2012 ) doi:10.1088/0264 - 9381/29/12/124007 [ arxiv:1111.7185 [ gr - qc ] ] .
s.  l.  detweiler , astrophys .
j.   * 234 * , 1100 ( 1979 ) .
doi:10.1086/157593 romani , r. w. ( 1989 ) .
timing a millisecond pulsar array . timing neutron stars , 113 - 117 .
l.  lentati _ et al .
_ , mon .  not .
soc .   * 453 * , 2576 ( 2015 ) [ arxiv:1504.03692 [ astro-ph.co ] ] .
z.  arzoumanian _ et al .
_ [ nanograv collaboration ] , arxiv:1508.03024 [ astro-ph.ga ] . e.  s.  phinney , astro - ph/0108028 . a.  vilenkin , phys .
rept .   * 121 * , 263 ( 1985 ) .
s.  kuroyanagi , k.  miyamoto , t.  sekiguchi , k.  takahashi and j.  silk , phys .
d * 87 * , no . 2 , 023522 ( 2013 ) [ phys .
d * 87 * , no .
6 , 069903 ( 2013 ) ] [ arxiv:1210.2829 [ astro-ph.co ] ] .
m.  maggiore , gr - qc/0008027 .
l.  p.  grishchuk , phys .
* 48 * , 1235 ( 2005 ) doi:10.1070/pu2005v048n12abeh005795 [ gr - qc/0504018 ] .
b.  allen and j.  d.  romano , phys .
d * 59 * , 102001 ( 1999 ) [ gr - qc/9710117 ] .
c.  m.  f.  mingarelli , t.  sidery , i.  mandel and a.  vecchio , phys .
d * 88 * , no . 6 , 062005 ( 2013 ) [ arxiv:1306.5394 [ astro-ph.he ] ] .
s.  r.  taylor and j.  r.  gair , phys .
d * 88 * , 084001 ( 2013 ) [ arxiv:1306.5395 [ gr - qc ] ] .
n.  seto and a.  taruya , phys .
lett .   * 99 * , 121101 ( 2007 ) doi:10.1103/physrevlett.99.121101 [ arxiv:0707.0535 [ astro - ph ] ] . n.  seto and a.  taruya , phys .  rev .
d * 77 * , 103001 ( 2008 ) [ arxiv:0801.4185 [ astro - ph ] ] . s.  j.  chamberlin and x.  siemens , phys .
d * 85 * , 082001 ( 2012 ) doi:10.1103/physrevd.85.082001 [ arxiv:1111.5661 [ astro-ph.he ] ] . j.  r.  gair , j.  d.  romano and s.  r.  taylor , phys .  rev .
d * 92 * , no .
10 , 102003 ( 2015 ) doi:10.1103/physrevd.92.102003 [ arxiv:1506.08668 [ gr - qc ] ] .
r.  jackiw and s.  y.  pi , phys .
d * 68 * , 104012 ( 2003 ) doi:10.1103/physrevd.68.104012 [ gr - qc/0308071 ] . m.  satoh , s.  kanno and j.  soda , phys .
d * 77 * , 023526 ( 2008 ) doi:10.1103/physrevd.77.023526 [ arxiv:0706.3585 [ astro - ph ] ] . c.  r.  contaldi , j.  magueijo and l.  smolin , phys .
* 101 * ( 2008 ) 141101 doi:10.1103/physrevlett.101.141101 [ arxiv:0806.3082 [ astro - ph ] ] .
t.  takahashi and j.  soda , phys .
lett .   * 102 * , 231301 ( 2009 ) doi:10.1103/physrevlett.102.231301 [ arxiv:0904.0554 [ hep - th ] ] .
j.  l.  cook and l.  sorbo , phys .
d * 85 * , 023534 ( 2012 ) [ phys .
d * 86 * , 069901 ( 2012 ) ] doi:10.1103/physrevd.86.069901 , 10.1103/physrevd.85.023534 [ arxiv:1109.0022 [ astro-ph.co ] ] .
i.  obata , t.  miura and j.  soda , phys .
d * 92 * , no .
6 , 063516 ( 2015 ) doi:10.1103/physrevd.92.063516 [ arxiv:1412.7620 [ hep - ph ] ] .
a. p. lightman , w. h. press , r. h. price , and s. a. teukolski , problem book in relativity and gravitation , 2nd ed .
( princeton university press , 1979 ) .
m. maggiore , gravitational waves , vol . 1 : theory and experiments ( oxford university press , 2008 )
. g. b. rybicki , a. p. lightman , radiative processes in astrophysics ( wiley - interscience , 1979 ) .
l. d. landau and e. m. lifshitz , the classical theory ( pergamon press , 1975 ) . c. misner , k. thorne and j. wheeler , gravitation , ( freeman 1973 ) .
b.  allen and a.  c.  ottewill , phys .
d * 56 * , 545 ( 1997 ) doi:10.1103/physrevd.56.545 [ gr - qc/9607068 ] .
u.  seljak and m.  zaldarriaga , phys .
lett .   * 78 * , 2054 ( 1997 ) doi:10.1103/physrevlett.78.2054 [ astro - ph/9609169 ] .
j. n. goldberg et al .
, journal of mathematical physics * 8 * , 2155 , 1967 . m.  anholm , s.  ballmer , j.  d.  e.  creighton , l.  r.  price and x.  siemens , phys .
d * 79 * , 084030 ( 2009 ) doi:10.1103/physrevd.79.084030 [ arxiv:0809.0701 [ gr - qc ] ] . l.  g.  book and e.  e.  flanagan , phys .
d * 83 * , 024024 ( 2011 ) doi:10.1103/physrevd.83.024024 [ arxiv:1009.4192 [ astro-ph.co ] ] . f.  a.  jenet and j.  d.  romano , am .
j.  phys .
* 83 * , 635 ( 2015 ) doi:10.1119/1.4916358 [ arxiv:1412.1142 [ gr - qc ] ] .
r.  w.  hellings and g.  s.  downs , astrophys .
j.   * 265 * , l39 ( 1983 ) . doi:10.1086/183954


rapid progress in the design and manufacture of optical fiber systems is a result of worldwide demand for ultra - high bit - rate optical communications .
this explains the growing interest of the soliton community in soliton - based optical fiber communication systems .
this area of research was considerably advanced in recent years  @xcite .
the most remarkable results include the application of the concept of the dispersion management to _ temporal optical solitons _ and soliton - based optical transmission systems , and the discovery of the so - called _ dispersion managed soliton_. high - speed optical communications require effective components such as high - performance broadband computer networks that can be developed by employing the concept of the bit - parallel - wavelength ( bpw ) pulse transmission that offers many of the advantages of both parallel fiber ribbon cable and conventional wavelength - division - multiplexing ( wdm ) systems  @xcite .
expanding development in the study of the soliton fiber systems has been observed in parallel with impressive research on their spatial counterparts , optical self - trapped beams or _
spatial optical solitons_. one of the key concepts in this field came from the theory of multi - frequency wave mixing and cascaded nonlinearities where a nonlinear phase shift is produced as a result of the parametric wave interaction  @xcite .
in all such systems , the nonlinear interaction between the waves of two ( or more ) frequencies is the major physical effect that can support coupled - mode multi - frequency solitary waves .
the examples of temporal and spatial solitons mentioned above have one common feature : they involve the study of solitary waves in multi - component nonlinear models .
the main purpose of this paper is to overview several different physical examples of multi - mode and/or multi - frequency solitary waves that occur for the pulse or beam propagation in nonlinear optical fibers and waveguides . for these purposes , we select three different cases : multi - wavelength solitary waves in bit - parallel - wavelength optical fiber links , multi - colour spatial solitons due to multistep cascading in optical waveguides with quadratic nonlinearities , and quasiperiodic solitons in the fibonacci superlattices .
we believe these examples display both the diversity and richness of the multi - mode soliton systems , and they will allow further progress to be made in the study of nonlinear waves in multi - component nonintegrable physical models .
because the phenomenon of the long - distance propagation of _ temporal optical solitons _ in optical fibers  @xcite is known to a much broader community of researchers in optics and nonlinear physics , first we emphasize _ the difference between temporal and spatial solitary waves_. indeed , for a long time stationary beam propagation in planar waveguides has been considered somewhat similar to the pulse propagation in fibers .
this approach is based on the so - called _ spatio - temporal analogy _ in wave propagation , meaning that the propagation coordinate @xmath0 is treated as the evolution variable and the spatial beam profile along the transverse direction in waveguides , is similar to the temporal pulse profile in fibers .
this analogy is based on a simple notion that both beam evolution and pulse propagation can be described by the cubic nonlinear schrdinger ( nls ) equation .
however , contrary to the widely accepted opinion , there is a crucial difference between temporal and spatial solitons . indeed , in the case of the nonstationary pulse propagation in fibers , the operation wavelength is usually selected near the zero point of the group - velocity dispersion .
this means that the absolute value of the fiber dispersion is small enough to be compensated by a weak nonlinearity such as that produced by the ( very weak ) kerr effect in optical fibers which leads to a relative nonlinearity - induced change in the refractive index .
therefore , nonlinearity in such systems is always weak and it should be well modeled by a cubic nls equation which is known to be integrable by means of the inverse - scattering technique . however , for very short ( e.g. , fs ) pulses the cubic nls equation describing the long - distance propagation of pulses should be corrected to include additional terms that would account for such effects as higher - order dispersion , raman scattering , etc .
all such corrections can be taken into account with the help of the perturbation theory  @xcite .
thus , in fibers nonlinear effects are weak and they become important only when dispersion is small ( near the zero - dispersion point ) affecting the pulse propagation over large distances ( of order of hundreds of meters or even kilometers ) .
the situation changes dramatically when we consider the propagation of multi - wavelength pulses with almost equal group velocities .
the corresponding model is described by a nonintegrable and rather complicated system of coupled nls equations , which we briefly discuss below .    in contrary to the pulse propagation in optical fibers ,
the physics underlying the stationary beam propagation in planar waveguides and bulk media is different . in this case
an optical beam is generated by a continuous wave ( cw ) source and it is time independent .
however , when the beam evolves with the propagation distance @xmath0 , it diffracts in the transverse spatial directions .
then , a nonlinear change in the refractive index should compensate for the beam spreading caused by diffraction _ which is not a small effect_. that is why to observe spatial solitons as self - trapped optical beams , much larger nonlinearities are usually required , and very often such nonlinearities are not of the kerr type ( e.g. they saturate at higher intensities ) .
this leads to the models of generalized nonlinearities with the properties of solitary waves different from those described by the integrable cubic nls equation .
propagation distances involved in the phenomenon of the beam self - focusing and spatial soliton propagation are of the order of millimeters or centimeters .
to achieve such large nonlinearities , one needs to use the optical materials with large nonlinearity - induced refractive index .
one of the possible way to overcome this difficulty is to use the so - called _ cascaded nonlinearities _ of noncentrosymmetric optical materials where nonlinear effects are accumulated due to parametric wave interaction under the condition of the wave phase matching .
such parametric wave - mixing effects generate novel classes of spatial optical solitons where resonant parametric coupling between the envelopes of two ( or more ) beams of different frequencies supports stable spatially localised waves even in a bulk medium ( see details in ref .
it is this kind of multi - component solitary waves that we discuss below .
a growing demand for high - speed computer communications requires an effective and inexpensive computer interconnection .
one attractive alternative to the conventional wdm systems is bpw single - mode fiber optics links for very high bandwidth computer communications  @xcite .
they differ from the wdm schemes in that no parallel to serial conversion is necessary , and parallel pulses are launched simultaneously on different wavelengths .
when the pulses of different wavelengths are transmitted simultaneously , the cross - phase modulation can produce an interesting _ pulse shepherding effect _
@xcite , when a strong ( `` shepherd '' ) pulse enables the manipulation and control of pulses co - propagating on different wavelengths in a multi - channel optical fiber link .    to describe the simultaneous transmission of @xmath1 different wavelengths in a nonlinear optical fiber , we follow the standard derivation  @xcite and obtain a system of @xmath1 coupled nonlinear schrdinger ( nls ) equations @xmath2 ) : @xmath3 { \displaystyle \qquad      + \chi_j \left ( |a_j|^2 + 2 \sum_{m \neq j } |a_m|^2
\right ) a_j      = 0 ,    } \end{array}\ ] ] where , for the @xmath4th wave , @xmath5 is the slowly varying envelope , @xmath6 and @xmath7 are the group velocity and group - velocity dispersion , respectively , and the nonlinear coefficients @xmath8 characterize the kerr effect
. equations  ( [ eq : nls_dim ] ) do not include the fiber loss , since the fiber lengths involved in bit - parallel links are only a small fraction of the attenuation length .
we measure the variables in the units of the central wavelength channel ( say , @xmath9 ) , and obtain the following normalized system of the @xmath1 coupled nls equations , @xmath10 { \displaystyle \qquad      + \gamma_j \left(|u_j|^2 + 2 \sum_{m\neq j } |u_m|^2\right ) u_j      = 0 ,    } \end{array}\ ] ] where @xmath11 , @xmath12 is the incident optical power in the central channel , @xmath13 , @xmath14 , so that @xmath15 . for the operating wavelengths spaced @xmath16 nm apart within the band @xmath17 nm , the coefficients @xmath18 and @xmath19 are different but close to @xmath20 . initially , in eq .
( [ eq : nls ] ) , we omit the mode walk - off effect described by the parameters @xmath21 ( so that @xmath22 ) .
this effect will be analysed later in this section .    to analyze the nonlinear modes , i.e. localized states of the bpw model ( [ eq : nls ] )
, we look for stationary solutions in the form , @xmath23 and therefore obtain the system of equations for the normalized mode amplitudes , @xmath24 { \displaystyle     \frac{1}{2 } \alpha_n \frac{d^2u_n}{dt^2 }      + \gamma_n \left ( |u_n|^2      + 2\sum_{m \neq n } |u_m|^2 \right ) u_n      = \lambda_n u_n ,   } \end{array}\ ] ] where @xmath25 , the amplitudes and time are measured in the units of @xmath26 and @xmath27 , respectively , and @xmath28 .    system ( [ eq : nls_nn ] ) has _ exact analytical solutions _ for @xmath1 coupled components , the so - called _
bpw solitons_. indeed , looking for solutions in the form @xmath29 , @xmath30 , we obtain the constraint @xmath31 , and a system of @xmath1 coupled algebraic equations for the wave amplitudes , @xmath32 in a special symmetric case , we take @xmath33 , and the solution of those equations is simple  @xcite : + @xmath34^{-1/2}$ ] .
analytical solutions can also be obtained in the _ linear limit _
, when the central frequency pulse ( at @xmath35 ) is large . then , linearizing eqs .
( [ eq : nls_nn ] ) for small @xmath36 , we obtain a decoupled nonlinear equation for @xmath37 and @xmath38 decoupled linear equations for @xmath39 . each of the latter possess a localized solution provided @xmath40 , where @xmath41 ^ 2 $ ] . in this limit
the central soliton pulse @xmath37 ( `` shepherd pulse '' ) can be considered as inducing an effective waveguide that supports a fundamental mode @xmath39 with the corresponding cutoff @xmath42 .
since , by definition , the parameters @xmath43 and @xmath44 are close to @xmath20 , we can verify that the soliton - induced waveguide supports maximum of two modes ( fundamental and the first excited one ) .
this is an important physical result that explains the effective robustness of the pulse guidance by the shepherding pulse .    to demonstrate a number of unique properties of the multi - channel bpw solitons , we consider the case @xmath45 in more details .
a comprehensive discussion of the case @xmath46 can be found in the preprint  @xcite .
we select the following set of the normalized parameters : @xmath47 , @xmath48 , and @xmath49 .
solitary waves of this four - wavelength bpw system can be found numerically as localized solutions of eqs .
( [ eq : nls_nn ] ) . figure  [ fig : bpw1 ] presents the lowest - order families of such localized solutions . in general
, they are characterized by @xmath38 parameters , but we can capture the characteristic features by presenting power dependencies along the line @xmath50 in the parameter space @xmath51 .
the power of the central - wavelength component ( @xmath35 ) does not depend on @xmath52 ( straight line @xmath53 ) .
thin dashed , dotted , and dash - dotted curves correspond to the three separate single - mode solitons of the multi - channel bpw system , ( 1 ) , ( 2 ) , and ( 3 ) , respectively , shown with the corresponding branches of ( 0 + 1 ) , ( 0 + 2 ) , and ( 0 + 3 ) two - mode solitons .
the latter curves start off from the bifurcation points on the @xmath37 branch at @xmath54 , @xmath55 , and @xmath56 , respectively .
close separation of those curves is the result of closeness of the parameters @xmath43 and @xmath44 for @xmath57 .
thick solid curves in fig .
[ fig : bpw1 ] correspond to the two- ( 1 + 2 ) and three - mode ( 0 + 1 + 2 ) localized solutions .
the latter solutions bifurcate and give birth to four - wavelength solitons ( 0 + 1 + 2 + 3 ) ( branch a - b ) .
two examples of such four - wave composite solitons are shown in fig .
[ fig : bpw1 ] ( bottom row ) .
the solution b is close to an exact sech - type solution at @xmath58 ( described above ) for @xmath45 , whereas the solution a is close to that approximately described in the linear limit in the vicinity of a bifurcation point .
importantly , for different values of the parameters @xmath59 , the uppermost bifurcation point for this branch ( open circle in fig . [
fig : bpw1 ] ) is not predicted by a simple linear theory and , due to the nonlinear mode coupling , it gets shifted from the branch of the central - wavelength soliton ( straight line ) to a two - mode branch ( 0 + 1 + 2 ) ( thick solid curve ) .    as a result , if we start on the right end of the horizontal branch and follow the lowest branches of the total power @xmath60 in fig .
[ fig : bpw1 ] , we pass the following sequence of the soliton families and bifurcation points : @xmath61 . if the modal parameters are selected closer to each other , the first two links of _ the bifurcation cascade _ disappear ( i.e. the corresponding bifurcation points merge ) , and the four - mode soliton bifurcates directly from the central - wavelength pulse , as predicted by the linear theory .
note however that the sequence and location of the bifurcation points is a function of the cross - section of the parameter space @xmath51 , and the results presented above correspond to the choice @xmath62 .    the qualitative picture of the cascading bifurcations preserves for other values of @xmath1 .
in particular , near the bifurcation point a mixed - mode soliton corresponds to the localized modes guided by the central - wavelength soliton ( shepherd ) pulse .
the existence of such soliton solutions is a key concept of bpw transmission when the data are launched in parallel carrying a desirable set of bits of information , all guided by the shepherd pulse at a selected wavelength .
effects of the walk - off on the multi - channel bpw solitons seems to be most dangerous for the pulse alignment in the parallel links . for nearly equal soliton components ,
it was shown long time ago @xcite that nonlinearity can provide an effective trapping mechanism to keep the pulses together . for the shepherding effect , the corresponding numerical simulations are presented in figs .
[ fig : bpw2](a - d ) for the four - channel bpw system .
initially , we launch a composite four - mode soliton as an unperturbed solution a [ see fig .  [
fig : bpw1 ] ] of eqs .
( [ eq : nls ] ) , without walk - off and centered at @xmath63 .
when this solution evolves along the propagation direction @xmath0 in the presence of small to moderate relative walk - off ( @xmath64 for @xmath65 ) , its components remain strongly localized and mutually trapped [ fig .
[ fig : bpw2](a , b ) ] , whereas it loses some energy into radiation for much larger values of the relative mode walk - off [ fig .
[ fig : bpw2](c , d ) ] .
recent progress in the study of cascading effects in optical materials with quadratic ( second - order or @xmath66 ) nonlinear response has offered new opportunities for all - optical processing , optical communications , and optical solitons  @xcite .
most of the studies of cascading effects employ parametric wave mixing processes with a single phase - matching and , as a result , two - step cascading  @xcite .
for example , the two - step cascading associated with type i second - harmonic generation ( shg ) includes the generation of the second harmonic ( @xmath67 ) followed by reconstruction of the fundamental wave through the down - conversion frequency mixing ( dfm ) process ( @xmath68 ) .
these two processes are governed by one phase - matched interaction and they differ only in the direction of power conversion .    the idea to explore more than one simultaneous nearly phase - matched process , or _
double - phase - matched ( dpm ) wave interaction _ , became attractive only recently  @xcite , for the purposes of all - optical transistors , enhanced nonlinearity - induced phase shifts , and polarization switching .
in particular , it was shown  @xcite that multistep cascading can be achieved by two second - order nonlinear cascading processes , shg and sum - frequency mixing ( sfm ) , and these two processes can also support a novel class of multi - colour parametric solitons  @xcite , briefly discussed below .      to introduce the simplest model of multistep cascading ,
we consider the fundamental beam with frequency @xmath69 entering a noncentrosymmetric nonlinear medium with a quadratic response . as a first step ,
the second - harmonic wave with frequency @xmath70 is generated via the shg process .
as a second step , we expect the generation of higher order harmonics due to sfm , for example , a third harmonic ( @xmath71 ) or even fourth harmonic ( @xmath72 )  @xcite .
when both such processes are nearly phase matched , they can lead , via down - conversion , to a large nonlinear phase shift of the fundamental wave  @xcite .
additionally , the multistep cascading can support _ a novel type of three - wave spatial solitary waves _ in a diffractive @xmath66 nonlinear medium , _ multistep cascading solitons_.    we start our analysis with the reduced amplitude equations derived in the slowly varying envelope approximation with the assumption of zero absorption of all interacting waves ( see , e.g. , ref .
@xcite ) . introducing the effect of diffraction in a slab waveguide geometry
, we obtain @xmath73 { \displaystyle   \qquad\qquad\qquad\qquad      + \chi_{2}a_{2 } a_{1}^{\ast}e^{-i\delta k_{2}z } = 0 ,     } \\*[9pt ] { \displaystyle       4
i k_{1 } \frac{\partial a_{2}}{\partial z }       + \frac{\partial^{2 } a_{2}}{\partial x^{2 } }       + \chi_{4 } a_{3 } a_{1}^{\ast } e^{-i\delta k_{3}z }     } \\*[9pt ] { \displaystyle   \qquad\qquad\qquad\qquad      + \chi_{5 } a_{1}^{2 } e^{i\delta k_{2}z }       = 0 ,     } \\*[9pt ] { \displaystyle       6 i k_{1}\frac{\partial a_{3}}{\partial z }       + \frac{\partial^{2 } a_{3}}{\partial x^{2 } }       + \chi_{3}a_{2}a_{1}e^{i\delta k_{3}z }       = 0 ,     } \end{array}\ ] ] where @xmath74 , @xmath75 , and @xmath76 , and the nonlinear coupling coefficients @xmath77 are proportional to the elements of the second - order susceptibility tensor which we assume to satisfy the following relations ( no dispersion ) , @xmath78 , @xmath79 , and @xmath80 .    in eqs .
( [ physeqns ] ) , @xmath81,@xmath82 and @xmath83 are the complex electric field envelopes of the fundamental harmonic ( fh ) , second harmonic ( sh ) , and third harmonic ( th ) , respectively , @xmath84 is the wavevector mismatch for the shg process , and @xmath85 is the wavevector mismatch for the sfm process .
the subscripts ` 1 ' denote the fh wave , the subscripts ` 2 ' denote the sh wave , and the subscripts ` 3 ' , the th wave .
following the technique earlier employed in refs .
@xcite , we look for stationary solutions of eq .
( [ physeqns ] ) and introduce the normalised envelope @xmath86 , @xmath87 , and @xmath88 according to the relations , @xmath89 { \displaystyle       a_{2 } = \frac{2 \beta k_{1}}{\chi_{2}}e^{2i\beta z + i\delta k_{2 } z } v ,     } \\*[9pt ] { \displaystyle       a_{3 } = \frac{\sqrt{2\chi_{2}}\beta      k_{1}}{\chi_{1}\sqrt{\chi_{5}}}e^{3i\beta z + i\delta k z } u ,    } \end{array}\ ] ] where @xmath90 . renormalising the variables as @xmath91 and @xmath92 , we finally obtain a system of coupled equations , @xmath93 { \displaystyle 2i \frac{\partial v}{\partial z } + \frac{\partial^{2 } v}{\partial x^{2 } } - \alpha v + \frac{1}{2 } w^{2 } + w^{\ast}u = 0 , } \\*[9pt ] { \displaystyle 3i\frac{\partial u}{\partial z } + \frac{\partial^{2 } u}{\partial x^{2 } } - \alpha_{1}u + \chi vw = 0 , } \\*[9pt ] \end{array}\ ] ] where @xmath94 and @xmath95 are two dimensionless parameters that characterise the nonlinear phase matching between the parametrically interacting waves .
dimensionless material parameter @xmath96 depends on the type of phase matching , and it can take different values of order of one .
for example , when both shg and sfm are due to quasi - phase matching ( qpm ) , we have @xmath97 $ ] , where @xmath98 .
then , for the first - order @xmath99 qpm processes ( see , e.g. , ref .
@xcite ) , we have @xmath100 , and therefore @xmath101 . when sfm is due to the third - order qpm process ( see , e.g. , ref .
@xcite ) , we should take @xmath102 , and therefore @xmath103 . at
last , when sfm is the fifth - order qpm process , we have @xmath104 and @xmath105 .
dimensionless equations ( [ normal ] ) present a fundamental model for three - wave multistep cascading solitons in the absence of walk - off . additionally to the type
i shg solitons ( see , e.g. , refs  @xcite ) , the multistep cascading solitons involve the phase - matched sfm interaction ( @xmath106 ) that generates a third harmonic wave .
two - parameter family of localised solutions consists of three mutually coupled waves .
it is interesting to note that , similar to the case of nondegenerate three - wave mixing  @xcite , eqs .
( [ normal ] ) possess an exact solution . to find it , we make a substitution @xmath107 , @xmath108 and @xmath109 , and obtain unknown parameters from the following algebraic equations @xmath110 valid for @xmath111 and @xmath112 .
equations ( [ exactsol ] ) have two solutions corresponding to _ positive _ and _ negative _ values of the amplitude ( @xmath113 ) .
this indicates a possibility of multi - valued solutions , even within the class of exact solutions .    in general , three - wave solitons of eqs .
( [ normal ] ) can be found only numerically .
figures [ fig : tr1](a ) and [ fig : tr1](b ) present two examples of solitary waves for different sets of the mismatch parameters @xmath114 and @xmath115 .
when @xmath116 [ see fig .
[ fig : tr1](a ) ] , which corresponds to an unmatched sfm process , the amplitude of the third harmonic is small , and it vanishes for @xmath117 .    to summarise different types of three - wave solitary waves , in fig .
[ fig : tr2 ] we plot the dependence of the total soliton power defined as @xmath118 .
it is clearly seen that for some values of @xmath119 ( including the exact solution at @xmath120 shown by two filled circles ) , there exist _ two different branches _ of three - wave solitary waves , and only one of those branches approaches , for large values of @xmath119 , a family of two - wave solitons of the cascading limit ( fig .  [ fig : tr2 ] , dashed ) . the slope of the branches changes from negative ( for small @xmath119 ) to positive ( for large @xmath119 ) , indicating a possible change of the soliton stability .
however , the detailed analysis of the soliton stability is beyond the scope of this paper ( see , e.g. , refs .
@xcite ) .
another type of multistep cascading parametric processes which involve only two frequencies , i.e. _ two - colour multistep cascading _
, can occur due to the vectorial interaction of waves with different polarization .
we denote two orthogonal polarization components of the fundamental harmonic ( fh ) wave ( @xmath121 ) as a and b , and two orthogonal polarizations of the second harmonic ( sh ) wave ( @xmath122 ) , as s and t. then , a simple multistep cascading process consists of the following steps .
first , the fh wave a generates the sh wave s via type i shg process .
then , by down - conversion
sa - b , the orthogonal fh wave b is generated .
at last , the initial fh wave a is reconstructed by the processes sb - a or ab - s , sa - a .
two principal second - order processes aa - s and ab - s correspond to _ two different components _ of the @xmath66 susceptibility tensor , thus introducing additional degrees of freedom into the parametric interaction .
different cases of such type of multistep cascading processes are summarized in table  [ tab : dpm ] .    to demonstrate some of the unique properties of the multistep cascading
, we discuss here how it can be employed for soliton - induced waveguiding effects in quadratic media . for this purpose
, we consider a model of two - frequency multistep cascading described by the principal dpm process ( c ) ( see table  [ tab : dpm ] above ) in the planar slab - waveguide geometry .
using the slowly varying envelope approximation with the assumption of zero absorption of all interacting waves , we obtain @xmath123 { \displaystyle       2 i k_{1}\frac{\partial b}{\partial z } + \frac{\partial^{2 } b }       { \partial x^{2 } } + \chi_2 s b^{\ast}e^{-i\delta k_2 z } = 0 ,    } \\*[9pt ] { \displaystyle       4 i k_{1 } \frac{\partial s}{\partial z }       + \frac{\partial^{2 } s}{\partial x^{2 } }       + 2 \chi_1 a^2 e^{i\delta k_1 z }       + 2 \chi_2 b^2 e^{i\delta k_2 z }       = 0 ,    } \end{array}\end{aligned}\ ] ] where @xmath74 , the nonlinear coupling coefficients @xmath77 are proportional to the elements of the second - order susceptibility tensor , and @xmath124 and @xmath125 are the corresponding wave - vector mismatch parameters .    to simplify the system ( [ eq_1 ] ) , we look for its stationary solutions and introduce the normalized envelopes @xmath126 , @xmath127 , and @xmath128 according to the following relations , @xmath129 , @xmath130 , and @xmath131 , where @xmath132 , @xmath133 , and @xmath134 , and the longitudinal and transverse coordinates are measured in the units of @xmath135 and @xmath136 , respectively .
then , we obtain a system of normalized equations , @xmath137    { \displaystyle      i \frac{\partial v}{\partial z } +       \frac{\partial^{2 } v}{\partial x^{2 } } - \alpha_1 v + \chi v^{\ast}w= 0 , }                        \\*[9pt ]    { \displaystyle      2i \frac{\partial w}{\partial z } +       \frac{\partial^{2}w}{\partial x^{2 } } - \alpha w+\frac{1}{2}(u^2+v^2)= 0 , }    \end{array}\ ] ] where @xmath138 , @xmath139 , and @xmath140 .
first of all , we notice that for @xmath141 ( or , similarly , @xmath142 ) , the dimensionless eqs .
( [ eq_n ] ) reduce to the corresponding model for the two - step cascading due to type i shg discussed earlier @xcite , and its stationary solutions are defined by the equations for real @xmath126 and @xmath128 , @xmath143    { \displaystyle     \frac{d^2 w}{d x^{2 } } - \alpha w + \frac{1}{2 } u^2 = 0 , }    \end{array}\ ] ] that possess a one - parameter family of two - wave localized solutions @xmath144 found earlier numerically for any @xmath145 , and also known analytically for @xmath146 , @xmath147 ( see ref .
@xcite ) .
then , in the small - amplitude approximation , the equation for real orthogonally polarized fh wave @xmath127 can be treated as an eigenvalue problem for an effective waveguide created by the sh field @xmath148 , @xmath149 v = 0.\ ] ] therefore , an additional parametric process allows to propagate a probe beam of one polarization in _ an effective waveguide _ created by a two - wave spatial soliton in a quadratic medium with fh component of another polarization .
however , this type of waveguide is different from what has been studied for kerr - like solitons because it is _ coupled parametrically _ to the guided modes and , as a result , the physical picture of the guided modes is valid , rigorously speaking , only in the case of stationary phase - matched beams . as a result , the stability of the corresponding waveguide and localized modes of the orthogonal polarization it guides is a key issue . in particular , the waveguide itself ( i.e. _ two - wave parametric soliton _ ) becomes unstable for @xmath150  @xcite .    in order to find the guided modes of the parametric waveguide created by a two - wave quadratic soliton
, we have to solve eq .
( [ eq_eigen ] ) where the exact solution @xmath148 is to be found numerically .
then , to address this problem analytically , approximate solutions can be used , such as those found with the help of the variational method  @xcite .
however , the different types of the variational ansatz used do not provide a very good approximation for the soliton profile at all @xmath114 .
for our eigenvalue problem ( [ eq_eigen ] ) , the function @xmath148 defines parameters of the guided modes and , in order to obtain accurate results , it should be calculated as close as possible to the exact solutions found numerically . to resolve this difficulty , below we suggest a novel `` almost exact '' solution that _ would allow to solve analytically many of the problems involving quadratic solitons _
, including the eigenvalue problem  ( [ eq_eigen ] ) .
first , we notice that from the exact solution at @xmath146 and the asymptotic result for large @xmath114 , @xmath151 , it follows that the sh component @xmath148 of eqs .
( [ eq_2 ] ) remains almost self - similar for @xmath152 .
thus , we look for the sh field in the form @xmath153 , where @xmath154 and @xmath155 are to be defined .
the solution for @xmath156 should be consistent with this choice of the shape for sh , and it is defined by the first ( linear for @xmath126 ) equation of the system ( [ eq_2 ] ) .
therefore , we can take @xmath126 in the form of the lowest guided mode , @xmath157 , that corresponds to an effective waveguide @xmath148 . by matching the asymptotics of these trial functions with those defined directly from eqs .
( [ eq_2 ] ) at small and large @xmath158 , we obtain the following solution , @xmath159 @xmath160 here , the third relation allows us to find @xmath154 for arbitrary @xmath114 as a solution of a cubic equation , and then to find all other parameters as functions of @xmath114 . for mismatches in the interval @xmath161
, the parameter values change monotonically in the regions : @xmath162 , @xmath163 , and @xmath164 .
it is really amazing that the analytical solution  ( [ eq_s]),([eq_p ] ) provides _ an excellent approximation _ for the profiles of the two - wave parametric solitons found numerically , with the relative errors not exceeding 1%3% for stable solitons
( e.g. when @xmath165 ) . as a matter of fact
, we can treat eqs .
( [ eq_s ] ) and  ( [ eq_p ] ) as an _
approximate scaling transformation _ of the family of two - wave bright solitons .
moreover , this solution allows us to capture some remarkable internal similarities and distinctions between the solitons existing in different types of nonlinear media .
in particular , as follows from eqs .
( [ eq_s ] ) and  ( [ eq_p ] ) , the fh component and the self - consistent effective waveguide ( created by the sh field ) have approximately the same stationary transverse profiles as for one - component solitons in a kerr - like medium with power - law nonlinear response  @xcite . for @xmath166 ( @xmath167 ) and @xmath168 ( @xmath169 )
our general expressions reduce to the known analytical solutions , and the fh profile is exactly the same as that for solitons in quadratic and cubic kerr media , respectively .
on the other hand , the strength of self - action for quadratic solitons depends on the normalized phase mismatch @xmath114 and , in general , the beam dynamics for parametric wave mixing can be very different from that observed in kerr - type media .
now , the eigenvalue problem ( [ eq_eigen ] ) can be readily solved analytically .
the eigenmode cutoff values are defined by the parameter @xmath119 that takes one of the discrete values , @xmath170 , where @xmath171^{1/2}$ ] .
number @xmath172 stands for the mode order @xmath173 , and the localized solutions are possible provided @xmath174 .
the profiles of the corresponding guided modes are @xmath175 where @xmath176 , @xmath177 is the hypergeometric function , and @xmath178 is the mode s amplitude which can not be determined within the framework of the linear analysis .    according to these results , a two - wave parametric soliton creates , a multi - mode waveguide and larger number of the guided modes can be observed for smaller @xmath114 .
figures  [ fig : al1](a , b ) show the dependence of the mode cutoff values @xmath179 for a fixed @xmath180 , and @xmath181 for a fixed @xmath114 , respectively . for the case
@xmath103 , the dependence has a simple form : @xmath182 ^ 2 $ ] .
because a two - wave soliton creates an induced waveguide parametrically coupled to its guided modes of the orthogonal polarization , the dynamics of the guided modes _ may differ drastically _ from that of conventional waveguides based on the kerr - type nonlinearities .
figures   show two examples of the evolution of guided modes . in the first example
[ see fig .
[ fig : wave_w](a - c ) ] , a weak fundamental mode is amplified via parametric interaction with a soliton waveguide , and the mode experiences a strong power exchange with the orthogonally polarized fh component through the sh field .
this process is accompanied by only a weak deformation of the induced waveguide [ see fig .  [
fig : wave_w](a )  dotted curve ] .
the resulting effect can be interpreted as a power exchange between two guided modes of orthogonal polarizations in a waveguide created by the sh field . in the second example , the propagation is stable [ see fig .  [
fig : wave_w](d ) ] .    when all the fields in eq .
( [ eq_n ] ) are not small , i.e. the small - amplitude approximation is no longer valid , the profiles of the three - component solitons should be found numerically . however , some of the lowest - order states can be calculated approximately using the approach of the `` almost exact '' solution ( [ eq_s]),([eq_p ] ) described above , which is presented in detail elsewhere  @xcite .
moreover , a number of the solutions and their families can be obtained in _
an explicit analytical form_. for example , for @xmath183 , there exist two _ families of three - component solitary waves _ for any @xmath152 , that describe soliton branches starting at the bifurcation points @xmath184 at : ( i )  the soliton with a zero - order guided mode for @xmath185 : @xmath186 , @xmath187 , @xmath188 , and ( ii )  the soliton with a first - order guided mode for @xmath103 : @xmath189 , @xmath190 , @xmath188 , where @xmath191 and @xmath192 .    for a practical realization of
the dpm processes and the soliton waveguiding effects described above , we can suggest two general methods .
the first method is based on the use of _ two commensurable periods _ of the quasi - phase - matched ( qpm ) periodic grating . indeed , to achieve dpm
, we can employ the first - order qpm for one parametric process , and the third - order qpm , for the other parametric process .
taking , as an example , the parameters for linbo@xmath193 and aa - s @xmath194 and bb - s @xmath195 processes @xcite , we find two points for dpm at about 0.89 @xmath196 m and 1.25 @xmath196 m .
this means that a single qpm grating can provide simultaneous phase - matching for two parametric processes .
for such a configuration , we obtain @xmath197 or , interchanging the polarization components , @xmath198 . the second method to achieve the conditions of dpm processes
is based on the idea of _ quasi - periodic qpm grating _
specifically , fibonacci optical superlattices provide an effective way to achieve phase - matching at _ several incommensurable periods _ allowing multi - frequency harmonic generation in a single structure .
we describe the properties of such structures in the next section .
for many years , solitary waves have been considered as _ coherent localized modes _ of nonlinear systems , with particle - like dynamics quite dissimilar to the irregular and stochastic behavior observed for chaotic systems  @xcite . however , about 20 years ago akira hasegawa , while developing a statistical description of the dynamics of an ensemble of plane waves in nonlinear strongly dispersive plasmas , suggested the concept of a localized envelope of random phase waves  @xcite . because of the relatively high powers required for generating self - localized random waves , this notion remained a theoretical curiosity until recently , when the possibility to generate spatial optical solitons by a partially incoherent source was discovered in a photorefractive medium  @xcite .
the concept of incoherent solitons can be compared with a different problem : the propagation of a soliton through a spatially disordered medium . indeed , due to random scattering on defects , the phases of the individual components forming a soliton experience random fluctuations , and the soliton itself becomes _ partially incoherent _ in space and time . for a low - amplitude wave ( linear regime )
spatial incoherence is known to lead to a fast decay . as a result
, the transmission coefficient vanishes exponentially with the length of the system , the phenomenon known as anderson localization  @xcite .
however , for large amplitudes ( nonlinear regime ) , when the nonlinearity length is much smaller than the anderson localization length , a soliton can propagate almost unchanged through a disordered medium as predicted theoretically in 1990 @xcite and recently verified experimentally  @xcite .    these two important physical concepts , spatial self - trapping of light generated by an incoherent source in a homogeneous medium , and suppression of anderson localization for large - amplitude waves in spatially disordered media , both result from the effect of strong nonlinearity .
when the nonlinearity is sufficiently strong it acts as _ an effective phase - locking mechanism _ by producing a large frequency shift of the different random - phase components , and thereby introducing _ an effective order _ into an incoherent wave packet , thus enabling the formation of localized structures . in other words ,
both phenomena correspond to the limit when the ratio of the nonlinearity length to the characteristic length of ( spatial or temporal ) fluctuations is small . in the opposite limit , when this ratio is large , the wave propagation is practically linear .    below we show that , at least for aperiodic inhomogeneous structures , solitary waves can exist in the intermediate regime in the form of _ quasiperiodic nonlinear localized modes_. as an example , we consider shg and nonlinear beam propagation in _ fibonacci optical superlattices _ , and demonstrate numerically the possibility of spatial self - trapping of quasiperiodic waves whose envelope amplitude varies quasiperiodically , while still maintaining a stable , well - defined spatially localized structure , _ a quasiperiodic envelope soliton_.      we consider the interaction of a fundamental wave with the frequency @xmath69 ( fh ) and its sh in a slab waveguide with quadratic ( or @xmath66 ) nonlinearity . assuming the @xmath66 susceptibility to be modulated and the nonlinearity to be of the same order as diffraction , we write the dynamical equations in the form @xmath199     { \displaystyle i\frac{\partial w}{\partial z } + \frac{1}{4 }     \frac{\partial^2 w}{\partial x^2 } + d(z ) u^2 e^{i\beta z } =   0 , }    \end{array}\ ] ] where @xmath200 and @xmath201 are the slowly varying envelopes of the fh and sh , respectively .
the parameter @xmath202 is proportional to the phase mismatch @xmath203 , @xmath204 and @xmath205 being the wave numbers at the two frequencies .
the transverse coordinate @xmath158 is measured in units of the input beam width @xmath206 , and the propagation distance @xmath0 in units of the diffraction length @xmath207 .
the spatial modulation of the @xmath66 susceptibility is described by the quasi - phase - matching ( qpm ) grating function @xmath208 . in the context of shg
, the qpm technique is an effective way to achieve phase matching , and it has been studied intensively  @xcite .    here
we consider a qpm grating produced by a quasiperiodic nonlinear optical superlattice .
quasiperiodic optical superlattices , one - dimensional analogs of quasicrystals  @xcite , are usually designed to study the effect of anderson localization in the linear regime of light propagation .
for example , gellermann _ et al .
_ measured the optical transmission properties of quasiperiodic dielectric multilayer stacks of sio@xmath209 and tio@xmath209 thin films and observed a strong suppression of the transmission  @xcite . for qpm gratings , a nonlinear quasiperiodic superlattice of litao@xmath193 , in which two antiparallel ferro - electric domains are arranged in a fibonacci sequence ,
was recently fabricated by zhu _
et al . _
@xcite , who measured multi - colour shg with energy conversion efficiencies of @xmath210 .
this quasiperiodic optical superlattice in litao@xmath193 can also be used for efficient direct third harmonic generation @xcite .
the quasiperiodic qpm gratings have two building blocks a and b of the length @xmath211 and @xmath212 , respectively , which are ordered in a fibonacci sequence [ fig .
[ fig : d_z](a ) ] .
each block has a domain of length @xmath213=l ( @xmath214=l ) with @xmath215=@xmath216 ( shaded ) and a domain of length @xmath217=@xmath218 [ @xmath219=@xmath220 with @xmath215=@xmath221 ( white ) . in the case of @xmath66
nonlinear qpm superlattices this corresponds to positive and negative ferro - electric domains , respectively .
the specific details of this type of fibonacci optical superlattices can be found elsewhere  @xcite . for our simulations presented below
we have chosen @xmath222= @xmath223= 0.34 , where @xmath224= @xmath225 is the so - called _
golden ratio_. this means that the ratio of length scales is also the golden ratio , @xmath226= @xmath224 .
furthermore , we have chosen @xmath227=0.1 .
the grating function @xmath208 , which varies between @xmath216 and @xmath221 according to the fibonacci sequence , can be expanded in a fourier series @xmath228 where @xmath229=@xmath230=0.52 for the chosen parameter values .
hence the spectrum is composed of sums and differences of the basic wavenumbers @xmath231=@xmath232 and @xmath233=@xmath234 .
these components fill the whole fourier space densely , since @xmath231 and @xmath233 are incommensurate .
figure  [ fig : d_z](b ) shows the numerically calculated fourier spectrum @xmath235 .
the lowest - order `` fibonacci modes '' are clearly the most intense .      to analyze the beam propagation and shg in a quasiperiodic qpm grating one
could simply average eqs .
( [ dynam ] ) .
to lowest order this approach always yields a system of equations with constant mean - value coefficients , which does not allow to describe oscillations of the beam amplitude and phase .
however , here we wish to go beyond the averaged equations and consider the rapid large - amplitude variations of the envelope functions .
this can be done analytically for periodic qpm gratings  @xcite . however , for the quasiperiodic gratings we have to resolve to numerical simulations .
thus we have solved eqs .
( [ dynam ] ) numerically with a second - order split - step routine . at the input of the crystal
we excite the fundamental beam ( corresponding to unseeded shg ) with a gaussian profile , @xmath236 we consider the quasiperiodic qpm grating with matching to the peak at @xmath237 , i.e. , @xmath238=@xmath237=82.25 .
first , we study the small - amplitude limit when a weak fh is injected with a low amplitude .
figures  [ fig : soliton](a , b ) show an example of the evolution of fh and sh in this effectively linear regime . as is clearly seem from fig .
[ fig : soliton](b ) the sh wave is excited , but both beams eventually diffract .
when the amplitude of the input beam exceeds a certain threshold , self - focusing and localization should be observed for both harmonics .
figures  [ fig : soliton](c , d ) show an example of the evolution of a strong input fh beam , and its corresponding sh . again
the sh is generated , but now the nonlinearity is so strong that it leads to self - focusing and mutual self - trapping of the two fields , resulting in a spatially localized two - component soliton , despite the continuous scattering of the quasiperiodic qpm grating .
it is important to notice that the two - component localized beam created due to the self - trapping effect is quasiperiodic by itself . as a matter of fact , after an initial transient its amplitude oscillates in phase with the quasiperiodic qpm modulation @xmath208 .
this is illustrated in fig .
[ fig : oscillations ] , where we show in more detail the peak intensities in the asymptotic regime of the evolution . the oscillations shown in fig .  [
fig : oscillations ] are in phase with the oscillations of the qpm grating @xmath208 , and we indeed found that their spectra are similar .
our numerical results show that the quasiperiodic envelope solitons can be generated for a broad range of the phase - mismatch @xmath238 .
the amplitude and width of the solitons depend on the effective mismatch , which is the separation between @xmath238 and the nearest strong peak @xmath235 in the fibonacci qpm grating spectrum [ see fig .  [
fig : d_z](b ) ] .
thus , low - amplitude broad solitons are excited for @xmath238-values in between peaks , whereas high - amplitude narrow solitons are excited when @xmath238 is close to a strong peak , as shown in fig .
[ fig : soliton](c , d ) .
to analyse in more detail the transition between the linear ( diffraction ) and nonlinear ( self - trapping ) regimes , we have made a series of careful numerical simulations  @xcite . in fig .
[ fig : transmission ] we show the transmission coefficients and the beam widths at the output of the crystal versus the intensity of the fh input beam , for a variety of @xmath238-values .
these dependencies clearly illustrate the universality of the generation of localised modes for varying strength of nonlinearity , i.e. a quasiperiodic soliton is generated only for sufficiently high amplitudes .
this is of course a general phenomenon also observed in many nonlinear isotropic media .
however , here the self - trapping occurs for quasiperiodic waves , with the quasiperiodicity being preserved in the variation of the amplitude of both components of the soliton .
we have overviewed several important physical examples of the multi - component solitary waves which appear due to multi - mode and/or multi - frequency coupling in nonlinear optical fibers and waveguides .
we have described several types of such multi - component solitary waves , including : ( i )  multi - wavelength solitary waves in multi - channel bit - parallel - wavelength fiber transmission systems , ( ii )  multi - colour parametric spatial solitary waves due to multistep cascading in quadratic materials , and ( iii )  quasiperiodic envelope solitons in fibonacci optical superlattices .
these examples reveal some general features and properties of multi - component solitary waves in nonintegrable nonlinear models , also serving as a stepping stone for approaching other problems of the multi - mode soliton coupling and interaction .
the work was supported by the australian photonics cooperative research centre and by a collaborative australia - denmark grant of the department of industry , science , and tourism ( australia ) .          for an overview of quadratic spatial solitons ,
see l.  torner , in : _
beam shaping and control with nonlinear optics _ , f. kajzer and r. reinisch , eds .
( plenum , new york , 1998 ) , p. 229 ; yu .  s.  kivshar , in : _ advanced photonics with second - order optically nonlinear processes _ , a.  d. boardman , l. pavlov , and s. tanev , eds .
( kluwer , dordretch , 1998 ) , p. 451