wave scattering in surface - corrugated guiding systems remains a quite interesting problem , both from the theoretical viewpoint and due to many applications ranging from the sound and electromagnetic propagation , to the conduction of electrons in mesoscopic wires . in spite of many important results obtained to date ( see , e.g. , refs .
@xcite ) , there is a number of open questions to be resolved .
one of such problems is to understand the interplay between different mechanisms of the scattering that are due to principally different terms in an effective scattering potential constructed with the use of various perturbative approaches .    in the recent papers
@xcite it was shown that apart from the well known terms in the potentials appearing in the first order perturbation theory , one should keep some specific terms that were neglected due to their seeming smallness .
it was demonstrated that these terms , although formally belonging to the second order terms , can play a decisive role in the scattering in some region of parameters .
since these terms are proportional to the square of the derivative of surface profiles , they can be associated with the _ square - gradient scattering _ ( sgs ) mechanism .
in contrast with ref.@xcite where a particular case of one rough surface was considered , below we develop a general approach by concentrating our attention on quasi - one dimensional waveguides with two corrugated surfaces that are either statistically uncorrelated or correlated .
our main interest is in the expression for the _ attenuation length _ @xmath0 ( also known as the scattering length , or the total mean free path ) , corresponding to specific @xmath1-th propagating mode .
we assume that the total number @xmath2 of these modes ( or conducting channels ) is large , @xmath3 , and that the lower and upper surface profiles ( boundaries ) of such quasi-1d waveguides are described by two random functions @xmath4 and @xmath5 , respectively .
in what follows , we consider three configurations of a particular interest :    1 .
the waveguide with the _ uncorrelated _ boundaries , @xmath4 and @xmath5 .
2 .   the waveguide with the _ antisymmetric _ boundaries , @xmath6 3 .   the waveguide with the _ symmetric _ boundaries , @xmath7    in our approach we assume that due to the multiple scattering of a travelling wave from the rough boundaries , the longitudinal wave number of an @xmath1-th propagating mode can be written as @xmath8 , where @xmath9 is its unperturbed value , @xmath10 the real part @xmath11 is responsible for a roughness - induced correction to the phase velocity of a given mode . as is known , the shift @xmath11 does not change the transport properties of a disordered system . therefore
, our further analysis shall be focused on the attenuation length @xmath0 only .    in order to properly take into account the surface irregularities we use the approach developed in ref .
specifically , by means of the canonical transformation we formally reduce the problem of the _ surface _ scattering to the _ bulk _ one described by the hamiltonian @xmath12 , @xmath13 with the effective potential @xmath14 .
it should be stressed that such a transformation does not mean an equivalence between the surface and bulk scattering .
the very point is that the obtained `` bulk '' hamiltonian @xmath12 is of very specific form that does not allow to use standard random matrix approaches based on a complete randomness in effective potentials .
the paper is organized as follows . in sec .
[ sec - problem ] we start with the general expressions describing the waveguide with any form of rough boundaries . in sec .
[ sec - dysoneq ] we write the dyson - type equation and derive approximate expressions for the surface scattering potential in the form , most useful for our perturbative approach . in sec .
[ sec - xi - tilde ] we obtain the fluctuating part of the potential for three configurations of a waveguide . in sec .
[ sec - averagegf ] we present a brief description of the formalism for obtaining the average green s function developed in ref
.  @xcite . here
we also specify the conditions of applicability of our method . in sec .
[ sec - ln - an ] we perform a detailed analysis of the attenuation length for three specific configurations of the waveguide . in sec .
[ sec - num - analysis ] we report the data that illustrates our analytical results and make clear the role of various terms in the expression for the scattering potential . in this section
we also discuss the so - called `` repulsion '' phenomenon discovered in ref .
, in view of its relation to the square - gradient scattering mechanism .
a brief comparison of the three configurations , as well as our concluding remarks , are presented in sec .
[ sec - conclusion ] .
in this section we derive the general expressions describing quasi-1d waveguides with any type of the boundaries .
then , we present them in the form , most suitable for our analysis of three particular waveguide configurations .    specifically , we consider the standard model of an open waveguide of the average width @xmath15 , stretched along the @xmath16 axis . the lower and upper surfaces of the waveguide are assumed to be described , respectively , by the profiles @xmath17 , and @xmath18 .
here @xmath19 is the root - mean - square roughness height that is assumed to be identical for both boundaries . in other words ,
the waveguide occupies the region @xmath20 of the @xmath21-plane . the fluctuating width @xmath22 is defined by @xmath23 , \qquad \langle w(x)\rangle = d\,,\ ] ] where the random functions @xmath24 ( with @xmath25 ) describe the roughness of the corresponding boundary .    in what follows
we consider the model for which both of the boundaries are statistically identical .
they are assumed to be described by a statistically homogeneous and isotropic gaussian random processes with the zero mean and unit variance ,    [ gen - xi ] @xmath26    here the angular brackets stand for the statistical averaging over different realizations of the surface profiles @xmath27 .
we also assume that the binary correlator @xmath28 decreases on the scale @xmath29 , with the normalization , @xmath30 .    for three configurations we analyze below
, the cross - correlations are as follows . in the first case of the uncorrelated boundaries we obviously have , @xmath31 the waveguide with the antisymmetric boundaries implies : @xmath32 as for the symmetric boundaries case , the opposite boundaries hold the correlation equality : @xmath33    the _ roughness - height power _ ( rhp ) spectrum @xmath34 is defined by : @xmath35 since @xmath28 is an even function of @xmath16 , its fourier transform ( [ ft - w ] ) is an even , real and non - negative function of @xmath36 .
the rhp spectrum has the maximum at @xmath37 with @xmath38 , and decreases on the scale @xmath39 .    in order to analyze the surface scattering problem
we shall employ the method of retarded green s function @xmath40 . specifically , we start with the dirichlet boundary - value problem    [ gen - gfp1 ] @xmath41 & & { \cal g}(x , x';z=\sigma \xi_{\downarrow}(x),z')=0 \nonumber \\ & & { \cal g}(x , x';z = d+\sigma \xi_{\uparrow}(x),z')=0 .
\label{gen - bc1}\end{aligned}\ ] ]    here @xmath42 and @xmath43 are the dirac delta functions .
the wave number @xmath44 is equal to @xmath45 for an electromagnetic wave of the frequency @xmath46 and te polarization , propagating through a waveguide with perfectly conducting walls . for an electron quantum wire
, @xmath44 is the fermi wave number within the isotropic fermi - liquid model .
the equation ( [ gen - geq1 ] ) does not contain any scattering potential , therefore , the perturbation is hidden in the boundary conditions ( [ gen - bc1 ] ) . in this case it is convenient to perform the canonical transformation to new coordinates ,    [ gen - coortr ] @xmath47 z_{new}&=&\frac{d\ , [ z_{old}-\sigma\xi_{\downarrow}(x_{old } ) ] } { w(x_{old})}\notag \\ & = & \frac{d\ , [ z_{old}-\sigma\xi_{\downarrow}(x_{old } ) ] } { d+\sigma[\xi_{\uparrow}(x_{old})-\xi_{\downarrow}(x_{old})]},\end{aligned}\ ] ]    in which both boundaries are flat . in accordance with eqs .
( [ gen - coortr ] ) , the transformations for the first and second derivatives with respect to the old coordinates are given by ( below for convenience we drop the subscript @xmath48 for @xmath16 and @xmath49 ) :    [ gen - d1tr ] @xmath50 \frac{\partial}{\partial z}\ , , \nonumber\\[6pt ]
\\ \frac{\partial}{\partial z_{old}}&=&\frac{\partial x}{\partial z_{old}}\frac{\partial}{\partial x}+\frac{\partial z}{\partial z_{old}}\frac{\partial}{\partial z } \nonumber \\ & = & \frac{d}{w(x)}\frac{\partial}{\partial z}\ , , \end{aligned}\ ] ]    [ gen - d2tr ] @xmath51 \frac{\partial}{\partial z } \right\ } \nonumber \\ & \times & \left\ { \frac{\partial}{\partial x}- \left[\frac{\sigma\xi'_{\uparrow}(x)}{w(x)}z + \frac{\sigma\xi'_{\downarrow}(x)}{w(x)}(d - z)\right ] \frac{\partial}{\partial z}\right\ } \nonumber \\[6pt ] \\
\label{gen - d2ztr } \frac{\partial^2}{\partial z_{old}^2}&=&\frac{d^2}{w^2(x ) } \frac{\partial^2}{\partial z^2}\,.\end{aligned}\ ] ]    correspondingly , we introduce the canonically conjugate green s function , @xmath52 here the prefactor before @xmath53 is due to the jacobian of the transformation to new variables . as a result
, we arrive at the equivalent boundary - value problem governed by the equation ,    [ gen - gfp4 ] @xmath54 & & -\widehat{u}(x , z)\,{\cal g}(x , x';z , z ' ) = \delta(x - x')\delta(z - z'),\qquad\label{gen - geq4}\\[6pt ] & & { \cal g}(x , x';z=0,z')=0,\nonumber\\ & & { \cal g}(x , x';z = d , z')=0 .
\label{gen - bc4}\end{aligned}\ ] ]    here we also omitted the subscript @xmath48 for @xmath55 . as a result of the transformation
, we have formally reduced the problem of surface scattering to that of the bulk one specified by the effective potential @xmath56 , @xmath57 \frac{\partial^2}{\partial z^2 } \nonumber \cr \\ & & + \frac{\sigma}{w(x ) } \bigg\{\bigg[\xi_{\uparrow}''(x ) + 2\xi_{\uparrow}'(x)\frac{\partial}{\partial x}\bigg ] \bigg[\frac{1}{2}+z\frac{\partial}{\partial z}\bigg ] \nonumber \cr \\ & & -\bigg[\xi_{\downarrow}''(x)+ 2\xi_{\downarrow}'(x)\frac{\partial}{\partial x}\bigg ] \bigg[\frac{1}{2}+(z - d ) \frac{\partial}{\partial z}\bigg ]   \bigg\}\nonumber \cr \\ & & -\frac{\sigma^2}{w^2(x)}\bigg\{{\xi_{\uparrow}'}^2(x ) \bigg[\frac{3}{4}+3z\frac{\partial}{\partial z } + z^2\frac{\partial^2}{{\partial z}^2}\bigg]\nonumber \cr \\ & & + { \xi_{\downarrow}'}^2(x ) \bigg[\frac{3}{4}+3(z - d)\frac{\partial}{\partial z } + ( z - d)^2\frac{\partial^2}{{\partial z}^2}\bigg ] \nonumber \cr \\ & & -\xi_{\downarrow}'(x)\,\xi_{\uparrow}'(x)\bigg[\frac{3}{2}+3(2z - d ) \frac{\partial}{\partial z}+2z(z - d)\frac{\partial^2}{{\partial z}^2 } \bigg ]   \bigg\}.\nonumber \cr \\\end{aligned}\ ] ] note that the prime over the function @xmath27 denotes a derivative with respect to @xmath16 .
it is important to stress that eqs .
( [ gen - gfp4 ] ) and ( [ gen - u ] ) are _ exact _ ones , and valid for any form of the surface profiles @xmath4 and @xmath5 .
as is expected , the general expression ( [ gen - u ] ) reduces to that obtained previously in ref.@xcite for a particular case of one rough boundary , if one substitutes @xmath58 and @xmath59 .
here we start with the non - local dyson equation in @xmath36-representation , @xmath60 with @xmath61 being the unperturbed pole - factor . in this representation
the effective surface scattering potential @xmath62 is given by the expression , @xmath63 here @xmath36 and @xmath64 are the lengthwise and transverse wave numbers , @xmath65 their unperturbed eigenvalues are @xmath9 , @xmath66 and @xmath67 , respectively .
the total number @xmath68 of propagating waveguide modes ( conducting channels that have real values of @xmath9 ) is determined by the integer part @xmath69 $ ] of the ratio @xmath70 , @xmath71.\ ] ]    after the substitution of the expression ( [ gen - u ] ) for @xmath72 in eq .
, we realize that @xmath62 consists of three groups of terms of a different nature .
the first one has the factor @xmath73 $ ] , the second and third groups have the factors @xmath74 and @xmath75 , respectively .
one needs to note that while itself the kernel @xmath62 is hermitian , the latter two parts are non - hermitian ones . in order to present each part of @xmath62 in an hermitian form
, we perform the integration by parts for the term containing @xmath76 and @xmath77 in the second group . as a result
, we arrive at the final _ exact _ form for the perturbation potential , @xmath78 \frac{\partial^2}{\partial z_1 ^ 2 } + \frac{\sigma}{w(x_1)}\ , [ i(k_x+k_x ' ) ] \nonumber \cr \\ & & \times \bigg\{\xi_{\uparrow}'(x_1)\bigg[\frac{1}{2 } + z_1\frac{\partial}{\partial z_1}\bigg ] -\xi_{\downarrow}'(x_1)\bigg[\frac{1}{2}+(z_1-d ) \frac{\partial}{\partial z_1}\bigg]\bigg\}\nonumber \cr \\ & & - \frac{\sigma^2}{w^2(x_1)}\bigg\{{\xi_{\uparrow}'}^2(x_1 ) \bigg[\frac{1}{4}+2z_1\frac{\partial}{\partial z_1 } + z_1 ^ 2\frac{\partial^2}{{\partial z_1}^2}\bigg]\nonumber \cr \\ & & + { \xi_{\downarrow}'}^2(x_1 ) \bigg[\frac{1}{4}+2(z_1-d)\frac{\partial}{\partial z_1 } + ( z_1-d)^2\frac{\partial^2}{{\partial z_1}^2}\bigg ] \nonumber\cr\\ & & - \xi_{\downarrow}'(x_1)\,\xi_{\uparrow}'(x_1)\bigg[\frac{1}{2 } + 2(2z_1-d)\frac{\partial}{\partial z_1 } \nonumber\cr\\ & & + 2z_1(z_1-d)\frac{\partial^2}{{\partial z_1}^2 } \bigg ] \bigg\ } \bigg ) \exp\left(ik'_xx_1\right)\frac{\sin(k'_zz_1)}{k'_z } \label{gen - xi - exactf}\end{aligned}\ ] ] this equation has a peculiar structure , very useful for a further analysis .
the kernel written in this form consists also of three groups of terms , however , now they can be associated with different scattering mechanisms .
since we are interested in the averaged green s function , we have to calculate the binary correlator of @xmath62 .
therefore , in order to avoid very cumbersome calculations , it is reasonable to make a simplification of @xmath62 that does not destroy the hermitian structure of each group of terms .
taking also into account that our main interest is in the case of small surface corrugations @xmath79 , we can do the following : expand the factor @xmath80/d$ ] in the first term of eq .
( [ gen - xi - exactf ] ) and to put @xmath81 in all the others .
in such a way , we derive the following _ approximate _ expression for the surface scattering potential , @xmath82\frac{\partial^2}{\partial z_1 ^ 2 } + \frac{\sigma}{d}\ , [ i(k_x+k_x ' ) ] \nonumber \cr \\ & & \times \bigg\{\xi_{\uparrow}'(x_1)\bigg[\frac{1}{2 } + z_1\frac{\partial}{\partial z_1}\bigg ] -\xi_{\downarrow}'(x_1)\bigg[\frac{1}{2}+(z_1-d ) \frac{\partial}{\partial z_1}\bigg]\bigg\}\nonumber \cr \\ & & - \frac{\sigma^2}{d^2}\bigg\{{\xi_{\uparrow}'}^2(x_1 ) \bigg[\frac{1}{4}+2z_1\frac{\partial}{\partial z_1 } + z_1 ^ 2\frac{\partial^2}{{\partial z_1}^2}\bigg]\nonumber \cr \\ & & + { \xi_{\downarrow}'}^2(x_1 ) \bigg[\frac{1}{4}+2(z_1-d)\frac{\partial}{\partial z_1 } + ( z_1-d)^2\frac{\partial^2}{{\partial z_1}^2}\bigg ] \nonumber\cr\\ & & - \xi_{\downarrow}'(x_1)\,\xi_{\uparrow}'(x_1)\bigg[\frac{1}{2 } + 2(2z_1-d)\frac{\partial}{\partial z_1 } \nonumber\cr\\ & & + 2z_1(z_1-d)\frac{\partial^2}{{\partial z_1}^2 } \bigg ] \bigg\ } \bigg ) \exp\left(ik'_xx_1\right)\frac{\sin(k'_zz_1)}{k'_z}\end{aligned}\ ] ]    the expression ( [ gen - xi - approx ] ) contains three types of scattering terms of a different physical meaning : the terms depending on the roughness height @xmath83 , derivative of the roughness height @xmath84 , and square of the derivative of the roughness height @xmath85 . in this connection
one can speak about the _ amplitude _ , _ gradient _ and _ square - gradient scattering _ ( sgs ) mechanisms . note that these terms are separately associated with each of two surfaces . apart from this
, there is a specific group of terms that are due to the inter - correlations between two surfaces . in what follows
we refer to these terms depending on the product @xmath86 , as to the _ gradient inter - product scattering _
( gis ) mechanism .
it should be stressed that the gis - terms , due to their proportionality to @xmath87 , were neglected in previous studies of the transport properties of surface - disordered waveguides .
a similar situation occurs for the sgs - terms that have been analyzed only recently in ref .
as we show below , all these terms can give the main contribution in the scattering , although formally they may be treated as negligible ones .
we start our analysis with the presentation of the kernel as a sum of its average and fluctuating parts , i.e. , @xmath88 it can be tested that the average , @xmath89 , contributes only to the real part @xmath11 of the complex renormalization @xmath90 of the lengthwise wave number @xmath9 ( see eq .
( [ k+dk ] ) ) , and therefore , does not change static transport properties of the surface disordered waveguide .
thus , we will omit it because our interest is in the attenuation length @xmath0 , not in @xmath11 .
therefore , it is necessary to identify the terms that have a non - zero mean - value in the general expression .
one can see that the mean value of the kernel in eq .   is defined by the sgs - terms that are proportional to @xmath91 . in order to construct the fluctuating part of @xmath62
, we introduce the _ square - gradient functions _ @xmath92 for the lower ( @xmath93 ) and upper ( @xmath94 ) boundary , by subtracting the mean values , @xmath95 , from @xmath96 , @xmath97 the functions @xmath98 play an important role in our further consideration . in accordance with the gaussian nature of the surface - profile functions
@xmath24 , there are no correlations between @xmath98 and both @xmath27 and @xmath99 , @xmath100 as well as between the product @xmath101 and @xmath102 , @xmath27 or @xmath99 .
the binary correlators for @xmath98 are given by @xmath103 correspondingly , the fourier transform of @xmath104 reads , @xmath105 due to the standard relation , @xmath106 here we have introduced the _ roughness - square - gradient power _ ( rsgp ) spectrum @xmath107 , @xmath108 we also introduce the fourier transform of the product @xmath109 and its binary correlator , @xmath110 \label{corr - ft - xi1-xi2 } & & \langle v_{\updownarrow}(k_x)v_{\updownarrow}(k_x')\rangle= 2\pi\delta(k_x+k_x')\,t(k_x),\end{aligned}\ ] ] as well as the fourier transform of @xmath27 , @xmath111    it should be pointed out that by the integration by parts the power spectrum of the roughness gradients @xmath99 can be reduced to the rhp spectrum @xmath34 .
however it is not possible to do for the rsgp spectrum @xmath107 .
this very fact reflects a highly non - trivial role of the sgs- and gis - terms , since they can compete with the other terms , although they are proportional to @xmath87 .
we can write now the fluctuating part of the total scattering potential in the form , @xmath112 here the first term is associated with the terms depending in eq .
on @xmath113 and @xmath114 , therefore , on the amplitude and gradient scattering mechanisms , respectively ,    [ ub - xi1u ] @xmath115\right\}\exp(ik_x'x_1)\,\frac{\sin(k_z'z_1)}{k_z ' } \qquad \label{ub - xi1u - def}\\[6pt ] & & \xi_{1\uparrow}(k_n , k_{n'})= \nonumber \\ & & -\sigma\{\delta_{nn'}+\cos[\pi(n - n')](1-\delta_{nn'})\ } \tilde\xi_{\uparrow}(k_n - k_{n ' } ) .
\label{ub - xi1u - exp}\end{aligned}\ ] ]    the second term in eq .   corresponding to the lower boundary , is of the same form ,    [ ub - xi1d ] @xmath116\right\ } \exp(ik_x'x_1)\,\frac{\sin(k_z'z_1)}{k_z ' } \qquad\quad \label{ub - xi1d - def}\\[6pt ] & & \xi_{1\downarrow}(k_n , k_{n'})=\sigma\ , \tilde\xi_{\downarrow}(k_n - k_{n ' } ) .
\label{ub - xi1d - exp}\end{aligned}\ ] ]    next two terms in eq .
are the sgs - terms , with @xmath117 being replaced by @xmath118 ,    [ ub - xi2u ] @xmath119 \nonumber\cr\\ & & \times\exp(ik_x'x_1)\,\frac{\sin(k_z'z_1)}{k_z ' } \label{ub - xi2u - def}\\[6pt ] & & \xi_{2\uparrow}(k_n , k_{n'})= \nonumber \\ & & \frac{\sigma^2d}{2 } \bigg\{\left[\frac{1}{3}+ \frac{1}{(2\pi n)^2}\right]\delta_{nn ' } + \frac{4}{\pi^2}\,\frac{n^2+n'^2}{(n^2-n'^2)^2 } \nonumber\cr \\ & & \times \cos[\pi(n - n')]\ , ( 1-\delta_{nn ' } ) \bigg\}\ , v_{\uparrow}(k_n - k_{n ' } ) .
\qquad \label{ub - xi2u - exp}\end{aligned}\ ] ]    [ ub - xi2d ] @xmath120 \nonumber\cr\\ & & \times\exp(ik_x'x_1)\,\frac{\sin(k_z'z_1)}{k_z ' } \label{ub - xi2d - def}\\[6pt ] & & \xi_{2\downarrow}(k_n , k_{n'})= \nonumber \\ & & \frac{\sigma^2d}{2}\bigg\{\left[\frac{1}{3}+ \frac{1}{(2\pi n)^2}\right]\delta_{nn ' } \nonumber\cr\\ & & + \frac{4}{\pi^2}\,\frac{n^2+n'^2}{(n^2-n'^2)^2}\ , ( 1-\delta_{nn ' } ) \bigg\}\,v_{\downarrow}(k_n - k_{n'}).\qquad\quad \label{ub - xi2d - exp}\end{aligned}\ ] ]    finally , the last term in eq .
is due to the gis terms in eq .  ,    [ ub - xi2lu ] @xmath121 \nonumber\cr\\ & & \times\exp(ik_x'x_1)\,\frac{\sin(k_z'z_1)}{k_z ' } \label{ub - xi2lu - def}\\[6pt ] & & \xi_{2\updownarrow}(k_n , k_{n'})= \nonumber \\ & & \frac{\sigma^2d}{2}\bigg\{\left[\frac{1}{3}- \frac{1}{2(\pi n)^2}\right]\delta_{nn'}-\frac{8}{\pi^2}\ , \frac{n^2+n'^2}{(n^2-n'^2)^2 } \nonumber\cr\\ & & \times\cos^2[\pi(n - n')/2]\,(1-\delta_{nn ' } ) \bigg\}\ , v_{\updownarrow}(k_n - k_{n ' } ) .
\label{ub - xi2lu - exp}\end{aligned}\ ] ]    the kernel @xmath122 written as the sum of specially designed terms , has a very convenient form .
first , all the terms are chosen to have zero average .
second , since the functions @xmath27 , @xmath123 and @xmath124 do not correlate with each other , their fourier transforms are also non - correlated :    [ ft - corr ] @xmath125    therefore , all the terms in eq .   determining the potential , are also non - correlated .
their autocorrelators are ,    [ ub - xi1-corr - q1 ] @xmath126 \label{ub - q1 } & q_1(k_n , k_{n'})=\sigma^2w(k_n - k_{n'}),\end{aligned}\ ] ]    [ ub - xi2-corr - q2 ] @xmath127 \label{ub - q2 } & q_2(k_n , k_{n'})=\frac{\sigma^4d^2}{2}t(k_n - k_{n ' } ) \bigg\{\left[\frac{1}{3}+\frac{1}{(2\pi n)^2}\right]^2\delta_{nn ' } \notag \\ & + \frac{16}{\pi^4}\,\frac{(n^2+n'^2)^2}{(n^2-n'^2)^4}\ , ( 1-\delta_{nn'})\bigg\},\end{aligned}\ ] ]    [ ub - xi2lu - corr - q2lu ] @xmath128 \label{ub - q2lu } & q_{2\updownarrow}(k_n , k_{n'})=\frac{\sigma^4d^2}{4}t(k_n - k_{n ' } ) \left\{\left[\frac{1}{3}- \frac{1}{2(\pi n)^2}\right]^2\delta_{nn'}\right .
\nonumber \\ & \left.+\frac{64}{\pi^4}\,\frac{(n^2+n'^2)^2}{(n^2-n'^2)^4}\ , \cos^4[\pi(n - n')/2]\,(1-\delta_{nn'})\right\}. \raisetag{-5pt}\end{aligned}\ ] ]    note that in spite of the absence of correlations between the upper and lower surface profiles , there are the correlations due to their derivatives since @xmath129 ( see eq .  ) .    as a result
, the correlator of the fluctuating scattering potential in the @xmath36-representation reads ,    [ ub - tilde - xi - corr ] @xmath130 \label{ub - q } & q(k_n , k_{n'})=2q_1(k_n , k_{n'})+2q_2(k_n , k_{n ' } ) + q_{2\updownarrow}(k_n , k_{n'})\ , .
\raisetag{-4pt}\end{aligned}\ ] ]    the factor @xmath131 in the first and second terms arises due to independent contributions , @xmath132 and @xmath133 , from the lower and upper boundaries .      in this configuration
both boundaries of the waveguide are identical , see eq .
, therefore , antisymmetric with respect to the central line @xmath134 . for this specific case
the width of the waveguide ( the distance between upper and lower boundaries ) remains constant along the @xmath135direction . instead of eq .
we have to start now with the exact expression ( ) by substituting the relation ( ) . as a result , one gets the relatively simple expression for the scattering potential , @xmath136\ , \xi'(x_1)\,\frac{\partial}{\partial z_1 } \\ & -\sigma^2\,{\xi'}^2(x_1)\frac{\partial^2}{{\partial z_1}^2}\bigg ) \exp\left(ik'_xx_1\right)\frac{\sin(k'_zz_1)}{k'_z } \end{split}\ ] ] one can realize that all terms proportional to @xmath87 can be written as one sgs - term with a non - zero mean value . by making use of the functions defined by eqs .  , and
with @xmath94 , in what follows we omit the subscript @xmath137 .
now the fluctuating part of the scattering potential is , @xmath138 the first summand related to the gradient scattering terms in eq .   has the form ,    [ ccb - xi1 ] @xmath139\,\xi'(x_1)\,\frac{\partial}{\partial z_1 } \exp(ik_x'x_1)\,\frac{\sin(k_z'z_1)}{k_z ' } \label{ccb - xi1-def}\\[6pt ] & & \xi_1(k_n , k_{n'})=2\sigma\sin^2[\pi(n - n')/2]\ , \tilde\xi(k_n - k_{n ' } ) .
\label{ccb - xi1-exp}\end{aligned}\ ] ]    the second summand is due to the sgs - term , however , with the function @xmath140 instead of @xmath141 . with the use of the fourier
transform ( [ ft - v ] ) for the operator @xmath142 , one gets ,    [ ccb - xi2 ] @xmath143 & & \xi_2(k_n , k_{n'})=\frac{\sigma^2d}{2}\,\delta_{nn'}\ , v(k_n - k_{n ' } ) .
\label{ccb - xi2-exp}\end{aligned}\ ] ]    note that the expressions ( [ ccb - xi1 ] ) and ( [ ccb - xi2 ] ) for the scattering potential are still _ exact_. in contrast with other configurations discussed below , the expansion in @xmath144-smallness is not required because of the constant width of the waveguide .
the correlation properties of the scattering potential are described by the following formulas .
due to the condition ( [ ft - corr - xi - v ] ) , the terms @xmath145 and @xmath146 are non - correlated , @xmath147 on the other hand , the autocorrelators are non - zero ,    [ ccb - xi1-corr - q1 ] @xmath148 .
\nonumber\\\label{ccb - q1}\end{aligned}\ ] ]    [ ccb - xi2-corr - q2 ] @xmath149    finally , the correlator of the fluctuating part of the total scattering potential @xmath150 has the following form ,    [ ccb - xi - tilde - corr - q ] @xmath151      let us now consider the waveguide with symmetric boundaries , see eq .
( [ sbcase ] ) , therefore , with the width defined by @xmath152 by substituting the symmetric profile into eq .
, we obtain for the surface scattering potential , @xmath153 \nonumber\cr\\ & & -\left.\frac{\sigma^2{\xi'}^2(x_1)}{d^2}\left[1 + 4\left(2z_1-d\right)\ , \frac{\partial}{\partial z_1}+ \left(2z_1-d\right)^2\frac{\partial^2}{\partial z_1 ^ 2}\right]\right\ } \nonumber\cr\\ & & \times\exp\left(ik'_xx_1\right)\frac{\sin(k'_zz_1)}{k'_z}\,.\end{aligned}\ ] ] as for the previous configuration , see eq .  , the sgs - terms corresponding to the lower boundary and the gis - terms are written in terms of @xmath154 .
the potential in eq .
( [ sb - xi - approx ] ) can be written as the sum of two terms related to substantially different scattering mechanisms , @xmath155 here the first summand is related to the terms proportional to @xmath19 in eq .
( the amplitude and gradient scattering terms ) ,    [ sb - xi1 ] @xmath156\right\ } \exp(ik_x'x_1)\,\frac{\sin(k_z'z_1)}{k_z ' } \label{sb - xi1-def}\\[6pt ] & & \xi_1(k_n , k_{n'})=-2\sigma\{\delta_{nn'}+\cos^2[\pi(n - n')/2 ] \nonumber\\ & & \times(1-\delta_{nn'})\}\,\tilde\xi(k_n - k_{n'}).\nonumber\\ \label{sb - xi1-exp}\end{aligned}\ ] ]    the second summand is associated with the fluctuating part of the sgs - terms in eq .
( [ sb - xi - approx ] ) , and consequently , with the function @xmath157 being replaced by @xmath158 ,    [ sb - xi2 ] @xmath159 \nonumber\cr\\ & & \times\exp(ik_x'x_1)\,\frac{\sin(k_z'z_1)}{k_z ' } \label{sb - xi2-def}\\[6pt ] & & \xi_2(k_n , k_{n'})=\frac{\sigma^2d}{2 } \bigg\{\left[\frac{1}{3}+ \frac{1}{(\pi n)^2}\right]\delta_{nn'}+\frac{16}{\pi^2}\,\nonumber\cr\\ & & \times\frac{n^2+n'^2}{(n^2-n'^2)^2 } \cos^2[\pi(n - n')/2](1-\delta_{nn'})\bigg\}v(k_n - k_{n'}).\nonumber\\ \label{sb - xi2-exp}\end{aligned}\ ] ]    the correlation properties of the introduced surface scattering potentials @xmath160 and @xmath161 are described by the correlator ,    [ sb - xi1-corr - q1 ] @xmath162(1-\delta_{nn'})\}. \label{sb - q1}\end{aligned}\ ] ]    [ sb - xi2-corr - q2 ] @xmath163 ^ 2 \delta_{nn'}\nonumber\\ & & + \frac{256}{\pi^4}\,\frac{(n^2+n'^2)^2}{(n^2-n'^2)^4 } \cos^4[\pi(n - n')/2](1-\delta_{nn'})\bigg\}.\nonumber\\ \label{sb - q2}\end{aligned}\ ] ]    thus , we arrive at the expression ,    [ sb - xi - tilde - corr - q ] @xmath164
now we replace the problem for the random green s function @xmath165 with the problem for the green s function @xmath166 averaged over the surface disorder . to perform the averaging of eq .
( [ g - de ] ) with @xmath167 , given for the uncorrelated , antisymmetric and symmetric boundary configurations by eq .  , and , respectively , and obtain @xmath168 , we can apply one of the standard and well known perturbative methods .
for example , it can be the diagrammatic approach developed for surface disordered systems  @xcite , as well as the technique developed in ref .  .
both of the methods take adequately into account the multiple scattering from the corrugated boundaries , and allow one to develop the consistent perturbative approach with respect to the scattering potential . after quite cumbersome calculations ( see ref .
@xcite for details ) , we obtain the average green s function in the normal - mode representation , @xmath169 the quantity @xmath0 is the wave attenuation length or electron total mean - free - path of the @xmath1th conducting mode .
this quantity is determined by the imaginary part of the self - energy and describes the scattering from the @xmath1th mode into other propagating modes , @xmath170 & & \times\left[q(k_n ,- k_{n'})+q(k_n,+k_{n'})\right ] .
\label{ln - q}\end{aligned}\ ] ] here the explicit expressions for @xmath171 and @xmath172 are given by eqs .  , and
for the waveguide with uncorrelated , antisymmetric and symmetric boundaries .    in deriving eqs .
( [ < gf > ] ) and ( [ ln - q ] ) we have made the following simplifications .
first , the self - energy in the dyson - type equation for the average green s function was obtained within the second - order approximation in the perturbation potential . in terms of the diagrammatic technique
this is similar to the  simple vortex " or , the same , bourret approximation  @xcite which contains the binary correlator @xmath173 of the surface scattering potential and the unperturbed pole factor @xmath174 .
second , in order to extract the inverse attenuation length from the self - energy , we substituted the lengthwise wave number @xmath36 by its unperturbed value @xmath9 .
the justification of this substitution , as well as the conditions of applicability resulting from the above simplifications , were discussed in ref .
specifically , it was shown that the domain of applicability is restricted by two independent criteria of the _ weak surface scattering _
, @xmath175 here @xmath176 is the distance between two successive reflections of the @xmath1th mode from one boundary to the opposite one . the conditions in imply that the wave is weakly attenuated on both the correlation length @xmath29 and the cycle length @xmath176 .
they implicitly include the requirement that @xmath79 for the surface corrugations be small in height , used in the sec .
[ sec - xi - tilde ] when deriving the explicit form ( [ gen - xi - approx ] ) for the surface scattering potential @xmath167 .
in this section we discuss the attenuation length @xmath0 for three configurations of the waveguide . in the first part of the subsections we give explicit expressions for the @xmath177 for any form of the binary correlator @xmath28 . in the second part
, we assume that the random surface profile @xmath178 has the gaussian binary correlator , @xmath179 then the rhp spectrum ( [ ft - w ] ) is given by @xmath180 the rsgp spectrum @xmath107 defined by eq .
( [ t - def ] ) ) , can be written as @xmath181 \nonumber\\[6pt ] & \times&\exp\left[-(k_xr_c)^2/4\right]\ , .
\label{t - gaus - expl}\end{aligned}\ ] ]    in what follows , we express the inverse attenuation length in the form of the dimensionless quantity @xmath182 .
since @xmath0 depends on as many as four dimensionless parameters @xmath183 , @xmath184 , @xmath70 , and @xmath1 , a complete analysis appears to be quite complicated .
for this reason , below we restrict ourselves to an analysis of the interplay between different scattering mechanisms , as a function of the dimensionless correlation length @xmath184 for different values of @xmath183 .
we concentrate our attention on a multimode waveguide with a large number of propagating modes , @xmath185 .
our main interest is in two situations , namely , in a small - scale roughness of the `` white noise '' kind for which @xmath186 , and in large - scale random corrugations with @xmath187 .      in view of eq .
( [ ub - q ] ) , the general expression ( [ ln - q ] ) for the inverse attenuation length shows that for a waveguide with uncorrelated boundaries this quantity consists of three terms , @xmath188 these terms are originated from different mechanisms of the surface scattering .
the first attenuation length @xmath189 is related to the rhp spectrum @xmath34 through the expression for @xmath190 . in accordance with eq .
( [ ub - q1 ] ) it is given by @xmath191 & & \times\left [ w(k_n+k_{n'})+w(k_n - k_{n'})\right].\qquad\qquad\end{aligned}\ ] ] here the term @xmath192 corresponding to @xmath193 , is related to the amplitude scattering while the terms @xmath194 result from the gradient scattering .
these two mechanisms of surface scattering are due to the corresponding terms in the expression for @xmath195 , see eqs .   and .
the first one is related to the term depending on the amplitude of the roughness profile @xmath196 , and the second one is due to the terms depending on the roughness gradient @xmath197 .
the expression ( [ ub - ln1-def ] ) exactly coincides with that previously obtained by various methods ( see , e.g. , ref .  ) .
the second attenuation length @xmath198 related to the rsgp spectrum through @xmath199 , is associated solely with the sgs mechanism due to the operator @xmath98 , see eqs .   and .
in accordance with eq .
( [ ub - q2 ] ) , it is described by @xmath200 here the terms in the sum are @xmath201 ^ 2 \left[t(2k_n)+t(0)\right]\ ] ] and @xmath202 & \times\left[t(k_n+k_{n'})+t(k_n - k_{n'})\right ] .
\end{split}\ ] ]    the third term @xmath203 is also related to the rsgp spectrum through @xmath204 due to the product @xmath101 , see eq .
( [ ub - xi2lu - def ] ) . in accordance with eqs .
( [ ub - q2lu ] ) , it is described by @xmath205 where @xmath206 ^ 2 \\[6pt ] & \times\left[t(2k_n)+t(0)\right ] .
\end{split}\ ] ] and @xmath207 & \times\left[t(k_n+k_{n'})+t(k_n - k_{n'})\right ] \cos^4[\pi(n - n')/2 ] .
\end{split}\ ] ]    analogously , the terms in eq .
can be classified as partial terms describing the intramode ( @xmath192 , @xmath208 and @xmath209 ) , and intermode ( @xmath194 , @xmath210 and @xmath211 ) scattering . according to the discussion of eq .
, the amplitude and gradient terms result in the intramode and intermode scattering , respectively . as for the square - gradient terms , they lead to both types of scattering .      to continue our analysis , let us start with the widely used case of _ a small - scale boundary perturbation _
, @xmath186 . in this case
the surface roughness can be regarded as a delta - correlated random process with the correlator @xmath212 and constant power spectrum @xmath213 .
taking into account the evident relationship @xmath214 , one can get the following inequalities specifying this case , @xmath215 it is necessary to underline that in the regime of small - scale roughness ( [ ssr - def ] ) the second of the weak - scattering conditions in eq.([ws - ln ] ) is not so restrictive as the first one , and directly stems from it , @xmath216 .    in eqs .
( [ ub - ln1-def ] ) , ( [ ub - lnn2 ] ) , ( [ ub - lnn2 ] ) , ( [ ub - lnn2lu ] ) and ( [ ub - lnn2lu ] ) for the attenuation lengths the argument of the correlators @xmath34 and @xmath107 turns out to be much less than the scale of their decrease @xmath39 under the conditions ( [ ssr - def ] ) . thus , for any term in the summation over @xmath217 the argument
can be taken as zero .
therefore , the first attenuation length is determined as follows ,    [ ub - ln1-ssr ] @xmath218 & & \approx ( k\sigma)^2\frac{n}{kd/\pi}\,\frac{kw(0)}{2 } \label{ub - ln1-ssr - w}\end{aligned}\ ] ]    the approximation is made here for a large number of the conducting modes , with the change of the summation over @xmath217 by integration . for a further estimate , one can take into account the formula @xmath219 that directly follows from the definition ( [ ft - w ] ) for the fourier transform @xmath34 of the binary correlator @xmath28 .
the function @xmath220 is the dimensionless correlator of the dimensionless variable @xmath221 , with the scale of decrease of the order of @xmath222 . as a result ,
the function @xmath220 does not depend on @xmath29 .
therefore , the integral over @xmath221 entering eq .
( [ w(0 ) ] ) is a positive constant of the order of unity .
for example , in the case of the gaussian correlations ( see eq .  ( [ wft - gaus ] ) ) , we have , @xmath223 , therefore , the integral is @xmath224 .
for the second length in eq .
one gets ,    [ ub - ln2-ssr ] @xmath225 & & \times\bigg\{\left[\frac{1}{3}+\frac{1}{(2\pi n)^2}\right]^2 + \frac{16k_nd}{\pi^4(\pi
n / d)^2 } \nonumber\\[6pt ] & & \times\left(\sum_{n'=1}^{n-1}+\sum_{n'=n+1}^{n_d}\right ) \frac{(\pi n'/d)^2}{k_{n'}d}\frac{(n^2+n'^2)^2}{(n^2-n'^2)^4}\bigg\ } \qquad\label{ub - ln2-ssr - sum}\\[6pt ] & & \approx\frac{\pi^2}{4}\frac{(k\sigma)^4}{k_nd } \frac{n^3}{(kd/\pi)}\frac{t(0)}{k^3}. \label{ub - ln2-ssr - nt}\end{aligned}\ ] ]    in eq .
( [ ub - ln2-ssr - sum ] ) every term in the sum rapidly decreases with an increase of the absolute value of @xmath226 .
this can be seen by making use of the following estimate , @xmath227 therefore , the sum in eq .
( [ ub - ln2-ssr - sum ] ) can be correctly estimated by three terms with @xmath228 . for simplicity , in eq .
( [ ub - ln2-ssr - nt ] ) we assume @xmath229 , and replace the curly braces by factor @xmath230 .    the explicit form for @xmath231 directly follows from the definition ( [ t - def ] ) for the correlator @xmath107 , @xmath232 & = & r_c^{-3}\int_{-\infty}^\infty d\rho \left[\frac{d^2{\cal w}(r_c\rho)}{d\rho^2}\right]^2.\end{aligned}\ ] ] if the roughness correlations are of the gaussian form , then according to eq .
( [ t - gaus - expl ] ) , we have @xmath233 and the integral over @xmath221 entering eq .
( [ t(0 ) ] ) is equal to @xmath234 .    the second term in eq .
is bigger than the third one .
therefore , a good aproximation for eq .
is @xmath235 one can make sure of this if we consider in eq .   the product of the fast decreasing factor and @xmath236 $ ] ( the latter cancels the terms with @xmath237 ) .
thus , we arrive at the following estimate , @xmath238 by comparing eq .   with eq
.   one can confirm eq .  .
the expression works well in the whole domain of applicability of our theory , not only in the region of small - scale roughness . to support this , in fig .
[ fig : l2_vs_l2lu ] we compare the plot of @xmath239 with the plot of @xmath240 ; they are practically the same .
plot of @xmath239 ( solid line ) and @xmath240 ( dashed line ) versus @xmath184 , for @xmath241 , @xmath242 and @xmath243 . ]    therefore , in what follows for this case of uncorrelated boundaries we will neglect the third term in eq .  .    due to
, the relationship between the first and second attenuation lengths is expressed by @xmath244 according to a substantially different behavior of the quantities @xmath245 and @xmath246 with respect to @xmath184 , it becomes clear that they have to intersect at the crossing point @xmath247 .
if the crossing point falls onto the region of a small - scale roughness ( @xmath248 ) , its dependence on the model parameters is obtained by equating the expression ( [ ub - l1/l2-ssr ] ) to unity , @xmath249 to the left from this point @xmath247 the sgs length prevails , @xmath250 . to its right
the main contribution is due to the first attenuation length , @xmath251 .
the expression ( [ ub - cp - ssr ] ) shows that the crossing point is smaller for smaller values of @xmath252 , as well as for smaller mode indices @xmath1 , or for larger values of the mode parameter @xmath70 .
the intermediate situation arises when the correlation length @xmath29 becomes much larger than the wave length @xmath253 , but still remains much less than the cycle length @xmath176 , @xmath254 as before , the first of the weak - scattering conditions ( [ ws - ln ] ) is the most restrictive , @xmath216 .
since the distance @xmath176 is larger than the correlation length @xmath29 , the successive reflections of the waves from the rough surface are weakly correlated .
note that the distance between neighboring wave numbers @xmath9 and @xmath255 is much smaller than the variation scale @xmath39 of the correlators @xmath34 and @xmath107 , @xmath256 this implies that the correlators @xmath34 and @xmath107 are smooth functions of the summation index @xmath217 .
therefore , the sum in the expression ( [ ub - ln1-def ] ) for the first attenuation length @xmath189 can be substituted by the integral ,    [ ub - ln1-wc ] @xmath257 & & \times\frac{1}{k}\left[w(k_n+k_{n'})+w(k_n - k_{n'})\right ] \\[6pt ] \label{ub - ln1-wc - fin } & & = \frac{(k\sigma)^2}{\pi}\frac{n}{kd/\pi}\,\int_{-k}^{k } dk_x\sqrt{k^2-k_x^2}\,\frac{w(k_n - k_x)}{k}.\qquad\quad\end{aligned}\ ] ]    eq .
( [ ub - ln1-wc ] ) shows that the first attenuation length is contributed by the scattering of a given @xmath1-th mode into other propagating modes . in order to obtain this asymptotic result ,
we have used only the condition of week correlations , @xmath258 .
therefore , eq .  ( [ ub - ln1-wc ] ) provides a reduction to eq .
( [ ub - ln1-ssr ] ) for small - scale corrugations , @xmath186 . in the case of a large - scale roughness , @xmath259 ,
the formula ( [ ub - ln1-wc ] ) obtained for @xmath189 allows one for further simplifications as was done in ref .  .
in contrast with @xmath189 , due to a rapidly decaying factor ( [ factordeltan ] ) the sgs length @xmath198 can be still described by keeping tree terms only , @xmath260 , in the sum in eq .
( [ ub - ln2 ] ) .
taking into account the estimate ( [ kn - wc ] ) , for the case @xmath229 one can obtain ,    [ ub - ln2-wc ] @xmath261 & & \times\bigg\{\left[\frac{1}{3}+\frac{1}{(2\pi n)^2}\right]^2 + \frac{8}{\pi^4}\bigg\ } \label{ub - ln2-wc - nt}\\[6pt ] & & \approx\frac{\pi^2}{8}\frac{(k\sigma)^4}{k_nd } \frac{n^3}{(kd/\pi)}\frac{t(0)+t(2k_n)}{k^3}. \qquad\quad\label{ub - ln2-wc - fin}\end{aligned}\ ] ]    to get the final expression ( [ ub - ln2-wc - fin ] ) , we replaced the curly braces in eq .
( [ ub - ln2-wc - nt ] ) by the factor @xmath230 .
naturally , at small - scale corrugations , @xmath186 , the obtained result ( [ ub - ln2-wc ] ) reduces to eq .
( [ ub - ln2-ssr ] ) . for a large - scale roughness , @xmath259
, one should use eq .
( [ ub - ln2-wc ] ) because of an arbitrary value of the parameter @xmath262 .
we do not consider here this case in detail due to its intermediate character .
another limit case refers to the correlation length @xmath29 to be very large both in comparison with the wave length @xmath253 and with the cycle length @xmath176 , @xmath263 in this case the number of wave reflections on the scale of the correlation length @xmath29 is large .
therefore , the successive reflections are strongly correlated to each other .    under the relations ( [ lsr - sc - def ] ) the second of the weak - scattering conditions ( [ ws - ln ] ) is the most restrictive , therefore , the condition of the applicability reads as @xmath264 the latter requirement determines the upper limit for the value of the correlation length @xmath29 .    due to eq .
( [ lsr - sc - def ] ) , the distance between neighboring wave numbers @xmath9 and @xmath255 turns out to be much larger than the variation scale @xmath39 of the correlators @xmath34 and @xmath107 , @xmath265 this indicates that the probability of the intermode ( @xmath266 ) transitions is exponentially small and the attenuation lengths are mainly formed by the incoherent intramode ( @xmath193 ) scattering . formally , for strong correlations ( [ lsr - sc - def ] ) the correlators @xmath34 and @xmath107 are sharpest functions of the summation index @xmath217 . in the sums of eqs .
( [ ub - ln1-def ] ) and ( [ ub - ln2 ] ) for the attenuation lengths the main contribution is due to the diagonal terms with @xmath193 , for which @xmath267 and @xmath268 can be neglected in comparison with @xmath269 and @xmath231 .
thus , the first attenuation length reads @xmath270 correspondingly , for the sgs length one gets , @xmath271 ^ 2.\ ] ]    the ratio of the first attenuation length to the second one can be presented as @xmath272 according to the first inequality in eq .
( [ ws - ln ] ) , applied to the condition ( [ lsr - sc - def ] ) of the strong - correlations and to the evident relationship @xmath273 , one can see that the amplitude scattering length always prevails over the sgs length within the interval of strong correlations .
for this reason in the condition of applicability ( [ domain - sc ] ) one should substitute @xmath0 by @xmath189 .
this allows one to arrive at the inequalities in the explicit form , @xmath274    before we pass to the next configuration of the waveguide , it is interesting to compare the case of two uncorrelated boundaries with the waveguide with one rough boundary , studied in ref .
one can found that under the approximation , the expression for @xmath177 turns out to be as twice as bigger than that one for the one - rough - boundary configuration .
this fact is quite natural since in the case of uncorrelated boundaries the two rough surfaces are of the same statistical nature ( see the factor @xmath131 in the first two terms of eq .  ) .
the remarkable similarity between the behavior of the inverse attenuation length for both configurations can be observed in fig .
[ fig : q1d-2rb_l1l2_n_krc ] .
note that the normalization factor @xmath176 for one boundary is twice as much of the quantity for two boundaries , and both curves are practically superposed .
plot of @xmath182 versus @xmath184 for @xmath241 , @xmath242 and different values of @xmath183 .
different sets of curves are shown for three waveguide configurations , as well as for the waveguide with one rough boundary analyzed in ref .
@xcite ( dashed line ) .
curves without symbols depict eq .   for the waveguide with uncorrelated boundaries ; curves labelled by triangles and diamonds stand for eqs .
and for waveguides with antisymmetric and symmetric boundaries , respectively ]    .
according to eqs .
( [ ccb - q ] ) and ( [ ln - q ] ) , the inverse attenuation length can be presented as a sum of the partial inverse lengths @xmath275 that are responsible for the scattering from @xmath1th into @xmath217th modes , @xmath276 however , it is noteworthy that in this case the partial terms describing the intramode @xmath277 and intermode @xmath278 scattering are rather different since they are associated with different scattering mechanisms . indeed , the intermode scattering results from the rhp spectrum only , through the expression in eq .
( [ ccb - q1 ] ) for the kernel @xmath279 .
these partial lengths are formed by the gradient scattering mechanism , due to the term containing in the scattering potential the derivative @xmath280 , see eq .  .
they read as @xmath281 \nonumber\cr\\ & & \times[w(k_n+k_{n'})+w(k_n - k_{n'})].\end{aligned}\ ] ]    from the above expression one can draw an interesting conclusion .
as one can see , in the case of symmetric boundaries the transitions are forbidden between the modes for which the difference of the indexes @xmath282 is an even number .
there are transitions only between odd and even channels .
this fact can be explained by an existence of an additional integral of motion , namely , the parity .
specifically , if one solves the stationary problem for the effective hamiltonian , there are two sets of eigenstates of different parity .
correspondingly , there are two independent sets of eigenvalues , see refs .  @xcite .    the intramode scattering , related to the rsgp spectrum through the kernel @xmath283 in the eq .  ,
is related to the sgs mechanism via the function @xmath158 in eq .  .
therefore , we have @xmath284.\ ] ] note that in a single - mode waveguide with @xmath285 the sum over @xmath217 in eq .
( [ ccb - ln ] ) contains only one term with @xmath286 . in this case
the backscattering length @xmath287 entering in eq .
( [ ccb - lnn ] ) is in accordance with that obtained in ref .
@xcite .
bearing in mind very different dependencies of the @xmath34 and @xmath107 power spectra on the roughness correlation length @xmath29 , one can write the inverse attenuation length as the sum of two terms .
one is related to @xmath34 , @xmath288 and the other is just the inverse of the intramode attenuation length determined by eq .  .
therefore , @xmath289 , and one can write , @xmath290    in fig .
[ fig : q1d-2rb_l1l2_n_krc ] we depict the length defined by eq .   as a function of @xmath184 for different values of @xmath252 ; the curve is labelled with triangles .    as above
, it is convenient to consider separately three regions of the correlation length .
we start with the standard case of a small - scale boundary perturbation specified by the conditions . in this case for any term in the summation over @xmath217 the argument of the correlators @xmath34 and @xmath107 is much less than the scale of their decrease @xmath39 , therefore , it can be taken as zero .
thus , for the attenuation length describing the intermode scattering one gets ,    [ ccb - ln1-ssr ] @xmath291 & & \times\sum_{n'=1}^{n_d}\frac{(\pi n'/d)^2}{k_{n'}d } \sin^4[\pi(n - n')/2]\qquad \label{ccb - ln1-ssr - sum1}\\[6pt ] & \approx & ( k\sigma)^2\frac{n}{kd/\pi}\,\frac{kw(0)}{2 } \label{ccb - ln1-ssr - w}\end{aligned}\ ] ]    the factor @xmath292 $ ] gives alternating values of one and zero for @xmath293 and @xmath294 , respectively , therefore , for a large number of conducting modes @xmath295 it can be replaced by @xmath296 .
the estimate ( [ ccb - ln1-ssr - w ] ) is obtained with the change of the summation by integration . for a further analysis one
has to take into account the formula in eq .   for @xmath269 .    the second term in eq
reads as @xmath297 where the explicit form for @xmath231 can be taken as in eq .  .
it is easy to realize that the relationship between lengths defined by eqs .   and
is of the same order as in the case of uncorrelated boundaries , see eq .  .
in this way one obtains the expression for the crossing point @xmath247 at which the curves of @xmath245 and @xmath246 intercept each other .
this intermediate situation is specified by eq .   .
we replace @xmath292 $ ] by @xmath296 , and change the sum by the integral since the correlator @xmath34 is a smooth function of @xmath217 ( see the discussion before eq .  ) , @xmath298 note that eq .
coincides with eq .   for the case of two uncorrelated boundaries .    for the length @xmath198 in the large - scale - roughness region , @xmath259
, one should use the general expression ( [ ccb - lnn ] ) because of arbitrary value of the parameter @xmath262 .
its normalized version reads @xmath299      in this limit case specified by eq .
, the distance between wave numbers @xmath9 and @xmath255 is much larger than the variation scale @xmath39 of the correlator @xmath34 , see eq .  .
this means that the probability of the intermode ( @xmath266 ) transitions is exponentially small and the total attenuation length is mainly formed by the incoherent intramode ( @xmath193 ) scattering , @xmath300 for strong correlations the correlator @xmath107 is a sharp function of its argument @xmath36 .
thus , the term with @xmath268 can be neglected in comparison with @xmath231 , and one gets , @xmath301    in this region of @xmath184 the sgs mechanism is recovered , determining the total attenuation length . the second inequality in eq .
is automatically satisfied for any value of @xmath184      as above , for this configuration it is convenient to write down the inverse attenuation lengths as a sum of two terms , @xmath302 here each term is associated with the corresponding term in eq .  , resulting from different mechanisms of the surface scattering .    the first term in eq .
is related to the rhp spectrum @xmath34 through the expression for @xmath303 . according to eq .
( [ sb - q1 ] ) , we have @xmath304\nonumber\\[6pt ] & & \times\left [ w(k_n+k_{n'})+w(k_n - k_{n'})\right].\qquad\qquad\end{aligned}\ ] ] the term @xmath192 corresponding to @xmath193 , is related to the amplitude scattering while the terms @xmath194 result from the gradient scattering
. these two mechanisms of surface scattering are due to the corresponding terms in the expression for @xmath305 , see eq .  .
namely , the former follows from the term depending on the amplitude of the roughness profile @xmath196 , and the latter from the terms depending on the roughness gradient @xmath197 .    the second term in eq .
related to the rsgp spectrum through @xmath306 , is associated solely with the sgs mechanism due to the square - gradient function @xmath307 , see eq .
( [ sb - xi2 ] ) . in accordance with eq .
( [ sb - q2 ] ) , it is described by @xmath308 here @xmath309 ^ 2 \left[t(2k_n)+t(0)\right]\ ] ] and @xmath310 & \times\cos^4[\pi(n - n')/2]\\[6pt ] & \times\left[t(k_n+k_{n'})+t(k_n - k_{n'})\right ] .
\end{split}\ ] ] note that in a single - mode waveguide with @xmath285 ( therefore with @xmath286 in eqs .   and )
the backscattering length @xmath311 takes the expression derived in ref .
@xcite .
the diagonal part of the inverse attenuation length , @xmath312 , characterizes the intramode scattering .
the off - diagonal partial attenuation length @xmath313 describes the intermode scattering .    to compare with the above case of the antisymmetric boundaries , here there are no transitions between the modes with odd values for the difference @xmath282 .
therefore , only the transitions between odd or even modes are allowed only .
this fact is related to the underlying symmetry , and can be associated with specific integrals of motion of the total hamiltonian .      in the case of a small - scale roughness
the first attenuation length reads    [ sb - ln1-ssr ] @xmath314 & & \approx ( k\sigma)^2\frac{n}{kd/\pi}\,\frac{kw(0)}{2}. \label{sb - ln1-ssr - w}\end{aligned}\ ] ]    for the second length we have ,    [ sb - ln2-ssr ] @xmath315 & & \times\bigg\{\left[\frac{1}{3}+\frac{1}{(\pi n)^2}\right]^2 + \frac{256k_nd}{\pi^4(\pi
n / d)^2 } \nonumber\\[6pt ] & & \times\left(\sum_{n'=1}^{n-1}+\sum_{n'=n+1}^{n_d}\right ) \frac{(\pi n'/d)^2}{k_{n'}d}\frac{(n^2+n'^2)^2}{(n^2-n'^2)^4 } \nonumber\\[6pt ] & & \times\cos^4[\pi(n - n')/2]\bigg\ } \qquad\label{sb - ln2-ssr - sum}\\[6pt ] & & \approx\frac{\pi^2}{8}\frac{(k\sigma)^4}{k_nd } \frac{n^3}{(kd/\pi)}\frac{t(0)}{k^3}.\label{sb - ln2-ssr - fin}\end{aligned}\ ] ]    as well as in the configuration with uncorrelated boundaries , due to the estimate every term in the sum of eq .
( [ sb - ln2-ssr - sum ] ) rapidly decreases with an increase of the absolute value of @xmath316 .
in addition , the factor @xmath317 is zero for @xmath293 .
therefore , the sum in eq .
( [ sb - ln2-ssr - sum ] ) can be evaluated by three terms with @xmath318 , and one can replace the quantity in the curly braces by @xmath230 . here
we also assume , @xmath229 .
it is evident that in the small - roughness case the crossing point is determined by eq .  .
in the intermediate situation the first attenuation length @xmath189 coincides with that obtained for the previous configurations , @xmath319 as for the second length @xmath198 , for @xmath229 it can be written , @xmath320      in this region of @xmath184 the first attenuation length reads , @xmath321 for the sgs length one gets , @xmath322 ^ 2.\ ] ]    the ratio of the first attenuation length to the second one can be obtained due to eq .  .
one can see that the amplitude scattering length always prevails over the sgs length within the interval of strong correlations , @xmath323 . for this reason ,
the condition of applicability coincides with eq .  .    in fig .
[ fig : q1d-2rb_l1l2_n_krc ] by diamonds we show the behavior of the total dimensionless inverse scattering length as a function of the correlation parameter @xmath184 .
in this section we present numerical data for periodic waveguides ( billiards ) with the potential @xmath56 .
our main attention is payed to the structure of eigenstates for the discussed types of boundaries .
although this information is not directly related to the scattering properties of open waveguides , it is instructive to reveal a role of the sgs - terms in the potential , by comparing the structure of eigenstates for different boundaries . the numerical method is described in details in ref.@xcite ; below we discuss the most interesting results and relate them to the previous studies ( see , ref.@xcite and references therein ) .
our particular interest is in the `` repulsion '' effect discovered in ref .
@xcite for quasi-1d billiards with rough surfaces .
this phenomena can be observed for a subset of eigenfunctions with the smallest values of the transversal quantum number of the unperturbed billiard with perfectly flat boundaries . for these eigenstates it was found that the maximum of the eigenfunction intensity is strongly shifted towards from the rough boundary .
this  repulsion " was shown to increase with a decrease of the roughness correlation length @xmath29 . based on our numerical data ,
one can make a conclusion that this effect of  repulsion " is directly related to the sgs - terms discussed above .    to start with , we show in fig .
[ fig : psi(x , z)-sgs - nointersgs ] a typical eigenstate of this type for the billiard with one rough surface ( upper panel ) . to compare with
, the low panel presents the corresponding eigenstate for the same potential , however , without sgs - terms .
as one can see , by removing these terms , the repulsion disappears .
it is also instructive to compare the structure of such eigenstates with those computed for the case of both symmetric and antisymmetric rough boundaries , see fig .
[ fig : psi(x , z)-antisym - sym ] . in both cases
the potentials have the sgs - terms .
however , the repulsion phenomena is absent in the antisymmetric billiard .
a close inspection reveals that the difference between these two cases is entirely due to the fact that in the antisymmetric billiard there is no square - gradient scattering between different channels , in contrast to the symmetric billiard .
thus one can make a conclusion that the main contribution to the repulsion is due to the intermode sgs - terms .
this conclusion supports an observation found in ref .
@xcite according to which the effect of the repulsion can be explained due to a strong localization in the channel space .
this localization occurs due to a relatively strong interaction between different conducting channels ( waveguide modes ) .
( color online ) upper panel : intensity of a typical eigenstate manifesting a  repulsion " from the rough surface .
low panel : the corresponding eigenstate for the potential without the sgs - terms.,title="fig : " ]   ( color online ) upper panel : intensity of a typical eigenstate manifesting a  repulsion " from the rough surface .
low panel : the corresponding eigenstate for the potential without the sgs - terms.,title="fig : " ]     ( color online ) intensity of an eigenfunction revealing strong repulsion from two symmetric rough boundaries ( low panel ) .
the same for the antisymmetric billiard for which there is no interaction between different channels ( upper panel).,title="fig : " ]   ( color online ) intensity of an eigenfunction revealing strong repulsion from two symmetric rough boundaries ( low panel ) .
the same for the antisymmetric billiard for which there is no interaction between different channels ( upper panel).,title="fig : " ]
we have developed a perturbative approach to the surface scattering in multimode quasi-1d waveguides with rough surfaces .
the approach is based on the construction of an effective hamiltonian that formally corresponds to the  bulk " scattering described by the potential of a very complicated structure .
the detailed analysis of this potential has led to the conclusion that there are specific square - gradient scattering ( sgs ) terms that were not taken into account in the previous studies .
we show that in spite of their seeming smallness in comparison with the standard terms , they can give a major contribution to the value of the scattering length in the case of a sufficiently strong roughness characterized by a corresponding value of the derivative of boundary profiles .    in order to demonstrate a strong influence of this sgs - mechanism ,
we have applied the approach to the waveguides with the boundaries of three different types . specifically , we considered the waveguides with two uncorrelated boundaries , and compared the results with those obtained for other two cases of symmetric and antisymmetric rough boundaries . in this way one can see the influence of various terms very clearly .
according to our analysis , all terms in the expression for the effective potential can be regrouped in three parts .
the first one is due to the roughness mean amplitude @xmath19 .
another group is characterized by the derivative of profiles , reflecting the fact that the scattering depends on two characteristic parameters , namely , on the amplitude @xmath19 ( or , the variance @xmath87 ) , and on the correlation length @xmath29 inversely proportional to the mean derivative of the profiles . in the expression for the scattering potential
the above two types of terms have linear dependence on the root - mean - square roughness height @xmath19 .
there is one more group of terms that can not be reduced to the two previous ones , and that gives rise to the square dependence on the profile derivative @xmath324 .
for this reason , in the previous studies it was assumed that these sgs - terms can be safely neglected .    however , we show that apart from the @xmath19-dependence due to the amplitude and gradient terms , as well as the @xmath87-dependence due to the sgs - terms , there is also dependence on the correlation parameter @xmath29 .
as was found , the @xmath29-dependence for the sgs - terms turns out to be very different from those corresponding to the amplitude and gradient terms .
our analysis reveals that there is a region of parameters where the sgs - contribution to the inverse scattering length is compared ( or even larger ) with the standard contributions due to the amplitude and gradient scattering .    in accordance with the general theory of surface scattering @xcite the wave propagation through any @xmath1-th conducting channel
is determined by its total ( or outgoing ) attenuation length @xmath0 .
the inverse attenuation length @xmath177 can be presented as the sum of partial inverse lengths , see eqs .  , and .
the first term @xmath325 is associated with both the amplitude and gradient scattering , it depends on the roughness - height power spectrum @xmath34 . to the contrary ,
the partial length @xmath326 being connected with the sgs - terms , depends on the power spectrum @xmath107 that is principally different from @xmath34 .
for the convenience of comparison between three different types of boundaries , in eq .
we have separately written the term @xmath327 that appears in the specific case of two uncorrelated boundaries .
this term also depends on @xmath107 .
such a representation allows us to analyze the total dependence of the inverse attenuation length @xmath177 on two control parameters , the roughness hight @xmath19 and the correlation length @xmath29 .
our analysis shows that one can speak about three characteristic regions of @xmath29 , namely , the small - scale , large - scale weak - correlation and large - scale strong - correlation roughness .    as we noticed above , the @xmath34- and @xmath107-power spectra have very different dependencies on the roughness correlation length @xmath29 .
specifically , the inverse value of the @xmath328-length decreases with an increase of the correlation parameter @xmath184 ; this occurs in the whole interval of the roughness values . as for the inverse value of the @xmath329-length , it shows a different behavior in each of the three regions of roughness and strongly depends on the waveguide configuration .
in the region of a small - scale roughness the inverse length @xmath330 increases with an increase of @xmath184 .
therefore , the curves corresponding to the terms @xmath325 and @xmath331 ( @xmath332 for the case of uncorrelated boundaries ) intersect each other , see fig
if the crossing point @xmath247 falls into this region of small - scale roughness , it obeys the law given by eq .
( [ ub - cp - ssr ] ) . to the left from the crossing point @xmath247 ,
the sgs - length @xmath328 prevails over the standard length , @xmath250 .
contrary , to the right from the point @xmath247 , the @xmath329-length gives the main contribution to the scattering process , @xmath251 .
the detailed analysis of three types of boundaries shows that for the small - scale roughness region one can specify that the crossover ( with an increase of @xmath184 ) occurs due to the following interplays :    * uncorrelated boundaries : the two sgs - terms that depend on the @xmath107-power spectrum , versus the amplitude and gradient terms depending on the @xmath34-power spectrum ; * antisymmetric boundaries : the sgs - term versus the gradient term ; * symmetric boundaries : the sgs - term versus the amplitude and gradient terms .    thus , we show that in the region of small - scale roughness and for any fixed value of the roughness height @xmath19 , one can indicate the region of small values of the correlation length @xmath29 where the sgs - length @xmath333 predominates over the the length @xmath334 .
this predominance occurs in spite of the fact that @xmath335 is proportional to @xmath87 while @xmath336 is proportional to @xmath337 .
it is interesting to note that in the region to the left from the point @xmath247 we have observed that the length @xmath338 for symmetric boundaries is twice as much as the length @xmath339 for the uncorrelated boundaries . on the other hand ,
the latter length is twice as much as the length @xmath340 for the antisymmetric boundaries , i.e. , @xmath341 .
there is another border that is characterized by the transition to the large - scale roughness with strong correlations , where a very different behavior of the total inverse length @xmath177 was found for the antisymmetric boundaries , in comparison with the uncorrelated and symmetric ones .
for the two last cases , the major mechanism of the scattering is due to the amplitude terms in the potential , and a similar behavior with respect to @xmath184 ( @xmath342 ) occurs for both of them .
however , the attenuation length for the uncorrelated boundaries turns out to be twice as much as the length for the symmetric boundaries , i.e. , @xmath343 ( note that this relation is inverted in comparison with the region of small - scale roughness ) .    as for the waveguide with the antisymmetric boundaries for which the amplitude terms are absent in the potential , the sgs - mechanism prevails in this region , leading to the dependence @xmath344 .
note , that here there is the interplay between the sgs - term and gradient one .
the transition point was found to be @xmath345 , with the crossover specified by the following interplays :      finally , there is an intermediate region of a large - scale roughness and weak correlations where the first attenuation length mainly contributes , @xmath346 . here
one can observe a strong similarity of the inverse attenuation length behavior for all configurations considered .
it should be noted that this region is quite difficult for obtaining simple analytical estimates due to a strong dependence of the number of terms in the corresponding sums on the parameter @xmath184 .
numerical data manifest that the attenuation length roughly remains the same inside this intermediate region .