short clean one - dimensional ( 1d ) mesoscopic wires , often referred to as quantum point contacts , show conductance quantization@xcite as a function of the channel width .
the quantization is well described by the theory of adiabatic propagation of free electrons.@xcite for non - interacting particles , conductance quantization should occur in longer channels too , as long as there is no backscattering off inhomogeneities within the channel .
a lot is known about the role of electron - electron interaction of 1d channels .
electron - electron repulsion in a wire enhances dramatically the reflection coefficient , making it energy - dependent.@xcite however , interaction between electrons does not alter the quantization ( in units of @xmath1 ) of an ideal channel conductance in the limit of zero temperature.@xcite what is still an open question is whether there are other manifestations of interactions due to inelastic processes , which influence the transport properties .    in the absence of interactions , left- and right - moving particles in a wire are at equilibrium with the reservoirs they originate from .
if a bias is applied between the reservoirs , then these equilibria differ from each other , giving rise to a particular form of the non - equilibrium distribution inside the channel . on the other hand ,
in a long ideal channel and in the presence of interactions one may expect equilibration to occur between the left- and right - movers into a single distribution characterized by an equilibrium with respect to a reference frame moving with some drift velocity .
interestingly , in a model with momentum - independent electron velocity for left- and right - movers ( as it is the case in tomonaga - luttinger model ) there is no difference between the two distributions .
effects originating from the particle - hole asymmetry , however , may discriminate between the two .
thermopower and coulomb drag  @xcite are examples of such effects .    .
this fabrication method has the advantage of producing reflectionless contacts to the leads,@xcite so that the boundary conditions of the distribution function is given by the fermi function of the reservoirs .
we define the thermopower as @xmath2 , i.e.  the voltage @xmath3 required to counteract a current due to the temperature difference @xmath4 .
[ fig : qpc - picture],scaledwidth=40.0% ]    at present , little is known about equilibration in a 1d electron system . in higher dimensions
the electron - electron interaction provides the most effective relaxation mechanism at low temperatures and therefore we include this relaxation mechanism as the first approach .
however , in 1d pair collisions can not change the distribution function for quadratic dispersion , since the momentum and energy conservation@xcite laws result in either zero momentum exchange or an interchange of the two momenta .
@xcite in either case the distribution function remains the same .
thus , the leading equilibration mechanism is due to three - particle collisions , which we study in this paper.@xcite    we investigate here the effects of three - particle collisions in reasonably short wires ( see fig .
[ fig : qpc - picture ] ) , where electron - electron scattering can be considered perturbatively . as measurable quantities ,
we evaluate the temperature dependence of the thermopower and conductance .
note that for more than one mode pair collisions become important for certain fillings.@xcite    the paper is organized as follows : first we review the non - interacting limit of thermopower and give a qualitative explanation of the effects due to three - particle collisions .
then we describe how to include the electron interactions using the boltzmann equation .
next we calculate the main ingredient for our perturbation theory , namely the three - particle matrix element and scattering rate using a @xmath5-matrix expansion .
we note several interesting properties of this scattering rate . finally , we derive the conductance and thermopower corrections and discuss the deviation from the so - called mott formula .
furthermore , some technical details are put in two appendices , and in appendix [ appendix - a - meso - princip ] we show that the number of left and right movers have to change in a scattering event for the current to change .      for a wire without interactions the distribution function @xmath6 is determined solely by the electron reservoirs @xmath7 where @xmath8 is the dispersion relation for momentum @xmath9 and spin @xmath10 ( suppressed in the notation ) , and @xmath11^{-1}$ ] is the fermi function with @xmath12 and @xmath13 denoting the chemical potential and temperature of the left / right contact , respectively ( see fig .
[ fig : qpc - picture ] ) .
the electric current for low temperature @xmath14 and in linear response to the applied bias @xmath3 and temperature difference @xmath15 then follows as ( @xmath16)@xcite @xmath17 from this the well - known leading - order results for conductance @xmath18 and for thermopower @xmath19 for a fully open channel are obtained . here
@xmath20 is the fermi temperature .      ) and cold right moving electrons ( @xmath21 ) . due to the temperature difference of the initial distributions
the scattering process creating left movers dominate compared to the opposite scattering and therefore it gives a positive correction to the thermopower .
[ fig : scattering - picture],scaledwidth=45.0% ]    one of the main results of this paper is that the three - particle collisions give a _ positive _ contribution to the thermopower , i.e.  the current due to a temperature difference is _ increased _ by the three - particle scattering .
this can be explained in simple terms .
firstly , to change the current the number of left and right moving electrons need to change , since it is the number of electrons going through a mesoscopic structure that determines the current and not their velocity ( see the appendix [ appendix - a - meso - princip ] ) .
secondly , we find the dominant scattering process at low temperature to only involve a single electron changing direction .
this occurs near the bottom of the band as pictured on fig .
[ fig : scattering - picture](a ) . for the initial electronic distribution
the left moving electrons have a higher temperature then the right moving ones , which favors scattering into the warmer distribution as seen on fig .
[ fig : scattering - picture](b ) .
this thus creates more left moving electrons and thereby increases the particle current towards the colder reservoir , i.e.  increasing the thermopower .
another important point is that the thermopower and conductance corrections are exponential in temperature , i.e.  proportional to @xmath22 .
this is a direct consequence of the dominant three - particle scattering process requiring an empty state near the bottom of the band .
we find the form of the thermopower correction at low temperatures due to the three particle scattering to be given by @xmath23 where @xmath3 is the electron - electron interaction strength and @xmath24 the fermi temperature .
this is found perturbatively in the short wire limit .
the long - wire limit remains an open question , and we expect that the length dependence of thermopower saturates once @xmath0 exceeds some relaxation length ( which increases for decreasing temperature ) .
in contrast , the conductance correction is negative . to understand this , note that the chemical potential of the initial distribution is higher for the right moving electrons then the left moving ones .
this favors scattering into the left moving branch ( still with the process shown in fig .
[ fig : scattering - picture](a ) ) for non - zero temperature and thereby decreasing the current .
the form of the conductance correction is similar to the thermopower correction : @xmath25
to model the current through a short 1d quantum wire including perturbatively the three particle interactions , we use the boltzmann equation @xmath26 , \label{eq : boltzmann - eq}\end{aligned}\ ] ] where @xmath27 is the distribution function at a space point @xmath28 between zero and @xmath0 ( see fig .
[ fig : qpc - picture ] ) , @xmath29 is the velocity and @xmath30 $ ] is the three - body electronic collision integral , i.e.  no impurity or interface roughness effects are included here .
we include the voltage and temperature difference in the boundary conditions of the reflectionless contacts,@xcite i.e.    [ eq : boundary - conditions ] @xmath31    and therefore omit the term @xmath32 in the boltzmann equation , allowed in the linear response regime.@xcite a similar method has been used to investigate electron - phonon interactions in short quantum wires,@xcite quantum hall effect in quantum wires @xcite and ballistic coulomb drag.@xcite    the three - particle collision integral is assumed to be local in space and is given by    @xmath33= & -\!\!\!\!\!\!\sum_{\substack{\sigma_2\sigma_3\\ \sigma_{1^{\prime}}\sigma_{2^{\prime}}\sigma_{3^{\prime } } } } \sum_{\substack{k_2k_3\\k_{1^{\prime}}k_{2^{\prime}}k_{3^{\prime}}}}\!\!\!\ ! w_{123;1^{\prime}2^{\prime}3^{\prime } } \left[f_1f_2f_3(1-f_{1^{\prime}})(1-f_{2^{\prime}})(1-f_{3^{\prime } } ) -f_{1^{\prime}}f_{2^{\prime}}f_{3^{\prime}}(1-f_1)(1-f_2)(1-f_3)\right ] , \label{eq : collision - integral - full}\end{aligned}\ ] ]    where the quantum numbers are primed / unprimed after / before the scattering event , @xmath34 and the scattering rate @xmath35 is found in the next section . without interactions ( @xmath36 ) the solution of the boltzmann equation is simply given by @xmath6 in eq .
( [ eq : f-0 ] ) .
when interactions are included it becomes a very difficult task to solve the boltzmann equation to all orders in the interaction .
however , for a short wire the interactions only have short time to change the distribution function away from the initial distribution @xmath6 and therefore we expand the distribution function in orders of @xmath35 as @xmath37 to find @xmath38 to the first order in @xmath39 , we insert the expansion of @xmath40 in the boltzmann equation and realize that only @xmath6 is necessary in the collision integral . since @xmath41=\mathcal{i}_{k}[f^{(0)}]$ ] is independent of @xmath28 , we find that    [ eq : f-1kk ] @xmath42   \quad & \textrm{for}\quad   & k>0,\\ f_{k}^{(1)}(x)&=\frac{x - l}{v_{k}}\mathcal{i}_{k}[f^{(0 ) } ]   \quad & \textrm{for}\quad & k<0 , \label{eq : f-1}\end{aligned}\ ] ]    using the boundary conditions eq .
( [ eq : boundary - conditions ] )
. therefore the current to the first order in @xmath39 is @xmath43\equiv i^{(0)}+i^{\textrm{int } } , \label{eq : i - int - def}\end{aligned}\ ] ] where @xmath44 is the non - interacting ( landauer ) part of the current from eq .
( [ eq : i-0 ] ) and @xmath45 is the part due to interactions .      the form of the interacting part of the current is now known and the next step is therefore to evaluate it to linear response to @xmath3 and @xmath4 to obtain the thermopower and conductance corrections . to this end , we define @xmath46 via @xmath47 where @xmath48 is the fermi function with temperature @xmath5 and fermi level @xmath49 .
it turns out that @xmath46 is proportional to either @xmath3 or @xmath4 .
this is seen by using the identity @xmath50 so we can identify @xmath46 by expanding the non - interacting distribution function @xmath51 ( see eq .
( [ eq : f-0 ] ) and fig .
[ fig : qpc - picture ] )    @xmath52ev,\\ f^0_{\textrm{r}}({\varepsilon}_k)&\simeq f^0({\varepsilon}_k)+\left[-{\partial}_{{\varepsilon}}f^0({\varepsilon}_k)\right]({\varepsilon}-{\varepsilon_{\textrm{f}}^{}})\frac{\delta t}{t } , \ ] ]    i.e. @xmath53 therefore to get @xmath45 in linear response to @xmath3 and @xmath4 , we linearized the collision integral @xmath54 $ ] eq .  ( [ eq : collision - integral - full ] ) with respect to @xmath55 and insert it into @xmath45 eq.([eq : i - int - def ] ) to obtain @xmath56 , \label{eq : i-1-linear - response}\end{aligned}\ ] ] where we defined @xmath57 using the shorthand notation @xmath58 and @xmath59 . to linearize the collision integral and thereby the correction to the current due to interactions @xmath45
, we have used the relation @xmath60 valid at @xmath61 .
since @xmath62 is different for positive and negative @xmath63 , we need to divide the summation in @xmath45 eq .
( [ eq : i-1-linear - response ] ) into positive and negative @xmath9 sums , which gives @xmath64 terms . for this purpose
, we introduce the notation @xmath65 and similarly for other combinations of the summation intervals .
the 32 terms can be simplified to only three terms using energy conservation and symmetry properties of @xmath66 in eq .
( [ eq : delta - def ] ) under interchange of indices .
there are pairwise exchanges of indices @xmath67 etc . and
interchanges between primed and unprimed indices , @xmath68 , using eq .
( [ eq : h - teorem ] ) and the fact that @xmath35 contains a matrix element squared .
this leads to six terms .
furthermore , @xmath66 is invariant under @xmath69 for all @xmath70 simultaneously due to time reversal symmetry , also seen explicitly from the form of @xmath35 ( derived below ) .
an example of how the simplifications occurs can be seen in eq .
( [ eq : simplification - of - sums ] )
. thus we obtain the more compact result    @xmath71\nonumber\\ & + 4(-e)\sum_{\textrm{spin}}\sum_{\substack{--+\\+++}}\delta_{123;1^{\prime}2^{\prime}3^{\prime } } \left[\frac{\delta t}{{k_{\textrm{b}}^{}}t^2 _ { } } \left\{({\varepsilon}_1-{\varepsilon_{\textrm{f}}^{}})+({\varepsilon}_2-{\varepsilon_{\textrm{f}}^{}})\right\}-\frac{2ev}{{k_{\textrm{b}}^{}}t}\right]\nonumber\\ & + 3(-e)\sum_{\textrm{spin}}\sum_{\substack{++-\\+-- } } \delta_{123;1^{\prime}2^{\prime}3^{\prime } } \left[\frac{\delta t}{{k_{\textrm{b}}^{}}t^2 _ { } } \left\{-({\varepsilon}_3-{\varepsilon_{\textrm{f}}^{}})+({\varepsilon}_{2^{\prime}}-{\varepsilon_{\textrm{f}}^{}})+({\varepsilon}_{3^{\prime}}-{\varepsilon_{\textrm{f}}^{}})\right\}-\frac{ev}{{k_{\textrm{b}}^{}}t}\right ] , \label{eq : i-1-linear - simplificeret - exact}\end{aligned}\ ] ]    where the definition of @xmath62 in eq .
( [ eq : psi-0-value ] ) was inserted .
an important point is that the number of positive / negative wave - vector intervals is not the same before and after the scattering .
therefore , we note that only scattering events that change the number of left and right moving electrons contributes to the interaction correction to the current .
the origin of this is the cancellation of the velocity in the definition of the current and in the distribution functions .
this cancellation thus leads to an expression for the interaction correction to the current in eq .   where all the in - going and out - going momenta enter on equal footing . in the appendix [ appendix - a - meso - princip ] ,
we show that this is valid to all orders in perturbation theory . due to this property and momentum conservation
, there are no processes that alter the current possible near the fermi level .
consequently , states far away from the fermi level have to be involved in the scattering , which , as we will see , leads to a suppression of @xmath45 by a factor @xmath22 . the distribution function on the other hand , _
can _ be changed by scattering processes near the fermi level .    to identify the important processes we find in the next section the scattering rate @xmath35 .
the three particle scattering rate @xmath35 is calculated using the generalized fermi golden rule inserting the @xmath5-matrix , @xmath72 , iterated to second order in the interaction @xmath3 to get the three particle interaction amplitude , i.e. , @xmath73 where @xmath74 is the initial energy , @xmath75 the final energy , @xmath76 is the resolvent operator ( or free green function ) , @xmath77 is short hand for @xmath78 , and the subscript @xmath79 means connected in the sense that the scattering process can not be effectively a two particle process , where one of the incoming particles does not participate in the scattering .
explicitly @xmath76 and @xmath3 are given by @xmath80 here @xmath81 is the unperturbed hamiltonian ( i.e.  kinetic energy with some dispersion ) , @xmath82 the fourier - transformed interaction potential , and @xmath83 ( @xmath84 ) is the annihilation ( creation ) operator . to calculate the matrix element @xmath85
, we write the initial and final states as @xmath86 where @xmath87 is the empty state . using the anti - commutator algebra @xmath88
, we obtain @xmath89 where we introduced @xmath90 ( the last equality is only valid for a quadratic dispersion ) , and where the set of permutations is given by @xmath91 . here
the signs of the permutation , @xmath92 , are shown as superscripts .    in order to exclude the effectively two - particle processes when multiplying eq .
( [ eq : vg|i > ] ) by @xmath93 from the left , @xmath94 ( @xmath95 ) needs to be different from @xmath96 ( @xmath97 ) .
the result is    @xmath98    where @xmath99 is the symmetrized interaction .
the matrix element consist of 36 terms and the scattering rate thus has @xmath100 terms . to obtain this result
we did not use energy conservation . for a quadratic dispersion ,
the denominator is only zero if we have an effective pair collision or if the momentum transfer is zero , as seen from the expression @xmath101 . a picture of the matrix element
is found in fig .
[ fig : diagram](b ) , where the exchange processes ( including the sign ) are visualized as different ways to connect two interaction lines and an intermediate propagation ( @xmath76 ) seen on fig .
[ fig : diagram](a ) .
the inclusion of the fermi statistics makes a substantial difference for the properties of the scattering rate as compared to the case described in ref .
[ ] , which is obtained by setting all @xmath102 .
we can rewrite the matrix element eq .
( [ eq : matrix - element - vgv - genneral ] ) in a more transparent way in terms of quantum mechanical exchange symmetry .
first , we introduce the following combination of three - particle scattering amplitudes @xmath103 , \label{eq : v - def}\end{aligned}\ ] ]    and after some rewriting , we then obtain @xmath104.\nonumber\end{aligned}\ ] ] we interpret this result in a way similar to a two - particle matrix element , @xmath105 , \label{eq:2-particle - matrix - element}\end{aligned}\ ] ] which contains a direct ( first term ) and a exchange term ( where @xmath106 ) .    in the three - particle case
@xmath107 is the direct term and one can make five exchange processes ( instead of one ) by exchanging the three final states @xmath108 , @xmath109 and @xmath110 .
this gives eq .
( [ eq : matrix - element - vha - v ] ) .
the sign in front of each @xmath111 is determined by the number of exchanges made , e.g.  in @xmath112 a single exchange , @xmath113 , gives a minus @xmath114 whereas for @xmath115 two exchanges ( @xmath116 followed by @xmath117 ) gives a positive sign @xmath118 .
furthermore , the arguments in @xmath107 are ordered in three pairs such that the differences between the elements in each pair are the only arguments of the interaction potential , see eq .
( [ eq : v - def ] ) .
this is useful when constructing approximations having a specific scattering process in mind .    how the matrix element was rewritten into the form of eq .
( [ eq : matrix - element - vha - v ] ) can also be described in terms of the drawings of fig .
[ fig : diagram ] .
the direct term @xmath107 is the sum of the six terms having mirror - symmetric exchanges before and after the scattering .
the other terms in eq.([eq : matrix - element - vha - v ] ) then can be obtained by suitable changes of outgoing lines .
the expressions we obtain for the three - particle scattering rates eq .
( [ eq : w - fermi - rule ] ) are quite cumbersome .
nevertheless , the obtained results allow for some consistency checks
. remarkably , for some two - body potentials , scattering of the particles of an @xmath119-body system is exactly equivalent to a sequence of two - body collisions .
such `` special '' potentials were studied in the context of integrable quantum many - body problems.@xcite we recall now that for a quadratic band a pair collision does not change the momenta of the incoming particles or simply permutes the two momenta .
therefore , three - particle scattering for a the integrable potentials may result _ only _ in permutations within the group of three momenta of the colliding particles ; all other three - particle scattering amplitudes must be zero for such potentials . in the context of this work it means that even three - particle ( or higher - order ) collisions would not bring electron equilibration for such types of electron - electron interaction .    in this section
we check that the three - particle scattering amplitudes are indeed zero for two special potentials .      in the case of contact interaction , @xmath120 , and for any kind of electron dispersion relation ( i.e. , not necessarily quadratic ) , we find by using the energy conservation law that @xmath121 , \label{eq : v - const - e - general}\end{aligned}\ ] ] where @xmath122 .
this is a major simplification from @xmath100 to @xmath123 terms by performing the spin summation .
if furthermore , the dispersion is quadratic , @xmath124 , then we find the ( at first sight ) surprising cancellation @xmath125 this can be seen directly from eq.([eq : v - const - e - general ] ) or by noting that @xmath126 for a quadratic dispersion and constant interaction for @xmath127 , i.e.  each term of eq .
( [ eq : matrix - element - vha - v ] ) is zero .
@xmath107 cancels in such a way , that the three first terms of eq .
( [ eq : v - def ] ) cancel each other ( the even permutations of @xmath128 combined with the same primed permutation ) and the three last terms cancel each other ( the odd permutations of @xmath128 combined with the same primed permutation ) .
in fact , the described above cancellation is in agreement with the general factorization results for the @xmath129matrix of one - dimensional @xmath119body problem with @xmath130function interaction in real space.@xcite in this context , it is crucial that the particles have quadratic dispersion relation ; if we use e.g.  @xmath131 , then the cancellation does not occur .
notice also that the cancellation we demonstrate is not a trivial zero .
indeed , the underlying two - particle amplitudes eq .
( [ eq:2-particle - matrix - element ] ) are finite for a @xmath132independent potential , if one includes spins .
( for spinless fermions and contact interaction the matrix element would be zero because the direct and the exchange terms cancel in accordance with the pauli principle . )
we checked also that the energy conserving part of the matrix element @xmath85 in the case of _ spinless _ fermions , quadratic dispersion , and the fourier transformed interaction potential of the form @xmath134 becomes equal to zero .
this is also possible to expect because of the relation of the potential eq .
( [ liebliniger ] ) to the integrable 1d bosonic lieb - liniger model  @xcite .
indeed , the bosonic model with contact interaction potential @xmath135 may be exactly mapped  @xcite onto the spinless fermionic model with interaction @xmath136 .
the integrability of the bosonic model guarantees the integrability of the corresponding fermionic one . adding a contact interaction to @xmath137 does no harm , as we are considering spinless fermions . finally , fourier transformation takes us to eq .
( [ liebliniger ] ) .
we checked that including the spin degree of freedom , spoils the remarkable cancellation for a three - particle amplitude .    in the following sections , we assume a general case interaction potential for which the three - particle scattering amplitudes lead to a non - trivial re - distribution of the momenta between the particles .
in this section , we go though the main ideas and approximations in evaluating the current correction due to interactions @xmath45 eq .  ( [ eq : i-1-linear - simplificeret - exact ] ) to lowest order in the temperature , @xmath14 .
we give a more detailed calculation in the appendix [ appendix - b - detailed - calcualtion ] .    as noted previously , all three terms in @xmath45 eq .
( [ eq : i-1-linear - simplificeret - exact ] ) are exponentially suppressed , since momentum conservation @xmath138 forbids scattering processes near the fermi level for the given combinations of positive and negative @xmath9 intervals . to be more specific
, it is the phase space restrictions of the fermi functions that give the exponential suppression , i.e. @xmath139 we begin by identifying the most important three particle scattering process .
the three terms in @xmath45 eq .  ( [ eq : i-1-linear - simplificeret - exact ] ) are : ( i ) two right movers backscattering a left mover while remaining right movers ; ( ii ) one right mover keeping its direction , while backscattering two left movers , and ( iii ) a left and a right mover keeping their directions , while backscattering the third particle .
from now on , we will concentrate on the case of coulomb interaction @xmath140 , which is the largest for small @xmath132 , therefore we want to identify processes where the initial and final state are close in momentum space  @xcite .
further , the process(es ) should not require more than one electron in states suppressed exponentially by the fermi functions .
one can see that due to the constraints stemming from momentum and energy conservation , in fact only process ( iii ) allows both initial and final states to be close to each other in momentum space and at the same time having only a single exponentially suppressed factor .
the corresponding scattering process is of the type shown in fig .
[ fig : scattering - picture](a ) .
therefore to the first order in @xmath22 , we include only the third one in eq .
( [ eq : i-1-linear - simplificeret - exact ] ) .
this leads to @xmath141 .
\label{eq : i - int - vigtigt}\end{aligned}\ ] ] here @xmath66 expresses the available phase space in form of the fermi functions and the three - particle scattering rate , see eq .
( [ eq : delta - def ] ) .
one essential approximation is that for the scattering process depicted in fig .
[ fig : scattering - picture](a ) , we may replace the full fermi distribution functions by the exponential tales or the low - temperature limit expressions , i.e.    [ f0s ] @xmath142    note that @xmath143 , @xmath144 and @xmath145 are all positive .
we see that the product of the fermi functions is indeed exponentially suppressed , i.e.  @xmath146 .
the second essential approximation is that for the scattering process seen in fig .
[ fig : scattering - picture](a ) the initial and final states differ by a small momentum .
therefore the matrix element in the transition rate @xmath35 is dominated by the direct term @xmath107 in eq .
( [ eq : matrix - element - vha - v ] ) , since the five exchange terms are suppressed by the coulomb interaction @xmath147 , i.e. @xmath148    the direct term eq .
( [ eq : matrix - element - direct - only ] ) would be zero for @xmath149 . in the case of quadratic dispersion relation and general @xmath150
, it vanishes in the limit @xmath151 and @xmath152 for @xmath153 due to the pauli principle . for a quadratic dispersion and for a general symmetrized interaction @xmath154
, the direct term @xmath107 simplifies to the following expression @xmath155}{(k_1-k_3+q_1)q_1q_3(k_1-k_3-q_3 ) } \label{eq : v - direct - quadratic - sec4}\end{aligned}\ ] ] where we used energy conservation and introduced @xmath156 and @xmath157 .
next , we give a qualitative explanation for the power law in @xmath5 for the interacting current correction eq .
( [ eq : i - int - vigtigt ] ) using the quadratic dispersion .
first , we consider the phase space constraint .
to do the sum over all @xmath9 in eq .
( [ eq : i - int - vigtigt ] ) we use the momentum and energy conservation and introduce new variables @xmath156 and @xmath157 , i.e. , change the summation variables , @xmath158 the energy conservation for a quadratic dispersion gives a factor of @xmath159 ( see e.g.  eq .
( [ eq : energy - conservation ] ) ) . for the process
at hand , the @xmath143 and @xmath160 are close to the fermi level and each of their sums contribute with a factor of @xmath161 and @xmath162 , respectively .
the fermi functions give the exponential suppression and a contribution to the phase space in form of an exponential tail , i.e. @xmath163 see eq .
( [ eq : fermi - function-1 - 1m]-[eq : fermi - function-3 - 3 m ] ) . to get the low temperature result for @xmath45 , eq .
( [ eq : i - int - vigtigt ] ) , we use the method of steepest decent to calculate the integral . to this end
, we note that the exponent @xmath164 is a function of @xmath161 and @xmath162 and in the limit @xmath165 the most important part is around the origin @xmath166 . here
@xmath164 vanishes as @xmath167 ( see appendix [ appendix - b - detailed - calcualtion ] for details ) . therefore collecting the phase space factors the current correction due to three particle interactions eq .
( [ eq : i - int - vigtigt ] ) becomes @xmath168.\end{aligned}\ ] ] in the limit @xmath14 .
furthermore , it turns out that the constraints @xmath169 and @xmath170 in the sum eq .
( [ eq : i - int - vigtigt ] ) only leaves phase space close to @xmath171 for @xmath165 , so we can set @xmath172 in the integrand and do the integral over @xmath162 , which is @xmath173 due to the phase space limits .
to lowest order in temperature this yields @xmath174 .
\label{eq : i - int - uden - fase - rum}\end{aligned}\ ] ] from this we conclude that phase space alone ( i.e.  assuming @xmath175 to be a constant ) gives a temperature dependence of the form
@xmath176 \quad \textrm{(phase space only)}.\end{aligned}\ ] ] however , as we have seen the three - particle interaction rate has delicate momentum dependence that needs to be taken into account . therefore , to calculate the direct interaction term @xmath107 we expand the symmetrized potential @xmath140 for small @xmath132 , as @xmath177 , \label{eq : vq - taylor - series}\end{aligned}\ ] ] where the parameter @xmath178 describes the screening due to the metallic gates near the quantum wire and @xmath179 is ( twice ) the @xmath180 fourier transform of the coulomb potential cut - off by the screening . setting @xmath181 into the three - particle scattering rate eq .
( [ eq : v - direct - quadratic - sec4 ] ) , we obtain @xmath182 to lowest order in @xmath132 . inserting this into the eq .
( [ eq : i - int - uden - fase - rum ] ) the final result for the current correction , including both phase space factors and the momentum dependent scattering rate , becomes @xmath183.\end{aligned}\ ] ] ( here we noticed that the non - constant three particle scattering rate gave rise to four extra powers in temperature . ) the detailed calculation given in appendix [ appendix - b - detailed - calcualtion ] yields a prefactor , and the end result is @xmath184 + \mathcal{o}\left[\bigg(\frac{t}{{t_{\textrm{f}}}}\bigg)^8\right ] .
\label{eq : i - int - resultat}\end{aligned}\ ] ] combining this result with the zero - order in the interaction terms , see eqs .
( [ g0 ] ) and ( [ s0 ] ) , we find for the thermopower and conductance in the low temperature limit , @xmath185 , \label{eq : s - result}\\ g= & \frac{2e^2}{h }     -\frac{2e^2}{h}e^{-{t_{\textrm{f}}}/t}\left[1+\frac{l}{\ell_{\textrm{eee}}}\right ] .
\label{eq : g - result}\end{aligned}\ ] ] here we introduced the effective length @xmath186 by the relation@xcite @xmath187 which may be viewed as a mean free path with respect to backscattering for a hole near the bottom of the band .    to recapitulate , the temperature dependence @xmath188 in eq .
can be understood in the following way : the three - particle scattering of a single particle leaves five free momenta , and since two are taken by energy and momentum conservation this gives @xmath189 .
in addition , the interaction , @xmath82 , is proportional to @xmath190 , and when squared it gives rise to four more powers , which results in the @xmath191 dependence .    in the limit of a point - like interaction , @xmath192 ,
the corrections are zero in agreement with the result of section [ subsec : point - like - int ] .
it is known from the luttinger liquid theory that in the limit of linear spectrum , which corresponds to @xmath193 , the conductance remains finite even if the wire is infinitely long ( @xmath194 ) .
therefore it is tempting to speculate that the two terms in the square brackets of eq .
( [ eq : g - result ] ) are the first terms of an expansion in @xmath195 of some function @xmath196 which saturates at a constant value in the limit @xmath197 .
one may also have a similar speculation generalizing eq .
( [ eq : s - result ] ) for the thermopower , @xmath198\to f_s(l/\ell_{\textrm{eee}})$ ] .    as a final remark
, we note that the so - called mott formula  @xcite relating the thermopower to the low - temperature conductance , @xmath199 is clearly violated by eqs .
( [ eq : g - result ] ) and ( [ eq : s - result ] ) .
this violation could be expected , because the conventional derivation of the mott formula ( for the non - interacting case ) assumes that the main contribution to the conductance _ and _ thermopower comes from the states around the fermi level , in an energy interval of the order of temperature.@xcite however , in the considered case the main contribution to @xmath129 comes from the `` deep '' states , even in the zeroth order with respect to the interaction potential .
correspondingly there is no surprise that eqs .
( [ eq : s - result ] ) and ( [ eq : g - result ] ) being substituted , respectively , in the left- and right - hand sides of eq .
( [ mott ] ) produce a parametrically large mismatch @xmath200 .
we have calculated the leading interaction correction to the transport properties of a clean mesoscopic wire adiabatically connected to the leads , using perturbation theory in the length of the wire .    for a single - mode wire ,
the leading interaction corrections turns out to be given by three - particle scattering processes .
this is because two - particle processes can not change the current due to momentum and energy conservation . to calculate the effect of the three - body processes
, we have utilized the boltzmann equation formalism , with three - particle scattering events defining the collision integral .
we have identified the leading - order scattering processes and found that they involve at least one state near the bottom of the band , i.e. far from the fermi level .
the involvement of such `` deep '' states results in an exponentially small , @xmath201 , interaction - induced correction to thermopower and conductance at low temperatures .
the account for interaction in this paper is performed for relatively short wires , where perturbation theory in the interaction or equivalently in the wire length is valid . for longer wires one needs to find the distribution function by treating the collision integral in the boltzmann equation non - perturbatively .
it is not clear whether the relaxation of the distribution function would instead yield non - exponential corrections to the transport coefficients for longer wires .
however , since the scattering processes that contribute to the current must involve a particle that changes direction ( which is proven in appendix [ appendix - a - meso - princip ] ) , one might speculate that the exponential suppression is valid for all lengths , as long as electron - electron scattering is the only active relaxation mechanism .
the question of what the relaxed distribution function looks like for a mesoscopic wire is an interesting and unsolved problem .
here we have only given a partial answer for the leading contributions for a short wire , i.e. to lowest order in the interaction .
further studies should involve a self - consistent determination of the distribution function .    since thermopower is sensitive to the electron distribution function
, it might be a good experimental tool for answering the fundamental questions regarding the effect of electron - electron collisions .
indeed refined measurements of thermopower of short 1d quantum wires have been performed yielding a reasonably good agreement with the free - electron theory.@xcite it remains an open question whether the accuracy of thermopower measurements is high enough to see the interaction effects in longer wires .
we acknowledge illuminating discussions with m. garst , a. kamenev , m. khodas , t. l. larsen and m. pustilnik . a.m.l .
appreciates and enjoyed the hospitality of the william i. fine theoretical physics institute , university of minnesota .
this work is supported by nsf grants dmr 02 - 37296 and dmr 04 - 39026 .
in this appendix , we show that _ the particle current change due to electronic scattering if and only if the scattering changes the number of left and right moving electrons_. in the main text ( see eq.([eq : i-1-linear - simplificeret - exact ] ) ) , this was shown to first order in the transition rate , but here we show it to _ all orders in the interaction_.    we show it explicitly in the boltzmann equation framework , however , suspect it to be a general feature of mesoscopic systems .
intuitively , the statement means that it is the _ number _ of particles that passes though the mesoscopic system that matters and not their velocity .
in contrast to this is e.g.  a long 1d wire or a bulk metal , where a velocity change of the particles is enough to change the current .    to show the above statement explicitly , we formally rewrite the boltzmann equation ( [ eq : boltzmann - eq ] ) including the boundary conditions eq.([eq : boundary - conditions ] ) as @xmath202}{v_{k } } \ \ \
\textrm{for}\ k>0 , \label{eq : solution - to - be - genneral - positive - k}\\ f_{k}(x)&=f^0_{\textrm{r}}({\varepsilon}_{k})+\int_l^x\!\!\ ! { \textrm{d}}x^{\prime}\frac{\mathcal{i}_{kx^{\prime}}[f]}{v_{k } } \ \ \ \textrm{for}\ k<0 .
\label{eq : solution - to - be - genneral - negative - k}\end{aligned}\ ] ] note that this is not a closed solution of the boltzmann equation , since the distribution function is still contained inside the collision integral .
however , this rewriting enables us to find the current without finding the distribution function first , i.e.  by inserting eq .
( [ eq : solution - to - be - genneral - positive - k],[eq : solution - to - be - genneral - negative - k ] ) into the current definition @xmath203 and obtain ( after a few manipulations ) : @xmath204\nonumber\\ & -\frac{(-e)}{l } \int_0^l\!\!\ ! { \textrm{d}}x \sum_{\sigma k<0}\mathcal{i}_{kx}[f ] + \frac{(-e)}{l } \int_0^x\!\!\ ! { \textrm{d}}x^{\prime } \overbrace{\sum_{\sigma k}\mathcal{i}_{kx^{\prime}}[f]}^{=0}\nonumber\\ \equiv&\   i^{(0)}+i^{(\textrm{int})},\end{aligned}\ ] ] where the @xmath28-dependent part can be seen to be zero by changing variables .
we note the cancellation of the velocity in the distribution function eqs .
( [ eq : solution - to - be - genneral - positive - k ] ) , ( [ eq : solution - to - be - genneral - negative - k ] ) and the current definition eq .
( [ eq : current - def ] ) , which is the origin of the statement we are showing ( as in the first order calculation ) .
a similar cancellation occurs in the landauer formula thus relating the transmission to the conductance . by using the explicit form of the collision integral eq .
( [ eq : collision - integral - full ] ) the current from the interactions is    @xmath205.\end{aligned}\ ] ]    we can divide the summation over @xmath9 quantum number into positive and negative intervals as in the main text ( see section [ sec : linear - response - limit ] ) . the essential point is now , that all terms that have the same number of positive ( and negative ) intervals for the primed and unprimed wave numbers @xmath9 are zero .
in other words , if the number of left and right moving electrons does not change then the contribution is zero by symmetry of the transition rate .
we show this cancellation in practice by an example ( using the notation of eq.([eq : notation - def-++ ] ) ) : @xmath206\label{eq : simplification - of - sums}\\ = & \sum_{\textrm{spin } } \sum_{\substack{-+-\\-+- } } w_{123;1^{\prime}2^{\prime}3^{\prime } } \left[f_1f_2f_3(1-f_{1^{\prime}})(1-f_{2^{\prime}})(1-f_{3^{\prime } } ) -f_{1^{\prime}}f_{2^{\prime}}f_{3^{\prime}}(1-f_1)(1-f_2)(1-f_3)\right]\nonumber\\ = & \sum_{\textrm{spin } } \sum_{\substack{-+-\\-+- } } w_{123;1^{\prime}2^{\prime}3^{\prime } } f_1f_2f_3(1-f_{1^{\prime}})(1-f_{2^{\prime}})(1-f_{3^{\prime } } ) - \sum_{\textrm{spin } } \sum_{\substack{-+-\\-+- } } \underbrace{w_{123;1^{\prime}2^{\prime}3^{\prime } } f_{1^{\prime}}f_{2^{\prime}}f_{3^{\prime}}(1-f_1)(1-f_2)(1-f_3)}_{\textrm{interchange}\ ( 123)\leftrightarrows(1^{\prime}2^{\prime}3^{\prime})}=0,\nonumber \ ] ]    interchanging @xmath108 and @xmath109 at the first equality using @xmath207 and interchanging @xmath208 in the second term as indicated .
thereby we have shown to all orders that to change the current by electronic interactions the number of left and right movers have to change .
the statement is _ not limited _ to only three particle scattering and can be show equivalently for pair interaction including several bands , electron phonon coupling or any other interaction with the same kind of symmetry under particle interchange .
furthermore , the statement is still true if the collision is non - local in space , since that only introduce some spatial integrals in the collision integral that can be handled similarly .
note however that the distribution function can be changed by processes that does not change the number of left and right movers .
the purpose of this appendix is to calculate @xmath45 in eq .
( [ eq : i - int - vigtigt ] ) @xmath141 .
\label{eq : i - int - vigtigt - appendix - b}\end{aligned}\ ] ] in the low - temperature limit , @xmath14 , step by step to find the prefactor given in eq .
( [ eq : i - int - resultat ] ) . as already mentioned , we preform the calculation with the scattering process seen in fig .  [
fig : scattering - picture](a ) in mind .
therefore we use the fermi functions as given in eq .   and the matrix element entering in the scattering rate from eq .
( [ eq : matrix - element - direct - only ] ) and eq .
( [ eq : v - direct - quadratic - sec4 ] ) , i.e.  using a quadratic dispersion .
we preform the summation over all the @xmath9 in eq .
( [ eq : i - int - vigtigt - appendix - b ] ) in the following way : first of all , we note that due to the momentum and energy conservation in the interaction process described , the scattering of @xmath160 to @xmath209 has to be from above to below the fermi level , i.e. @xmath210 this is due to the signs of @xmath145 and @xmath211 and can be understood as a sign of the difference between the curvature of the dispersion near the bottom of the band and near the fermi level .
next , we introduce the momentum transfer around the fermi level @xmath212 for @xmath213 and using the momentum conservation to do the @xmath211 summation , we obtain @xmath214 remembering the constraint @xmath215 , @xmath216 and @xmath217
. the fermi factors @xmath218 and @xmath219 restrict the momentum transfer @xmath161 and @xmath162 to be much smaller then @xmath220 and the @xmath143 and @xmath160 to be near the fermi level for the process in mind .
therefore we can use the fermi functions @xmath218 to do the summation over @xmath143 . assuming slow variation of the scattering rate over a range of @xmath221 at the fermi level , the @xmath143 summation becomes @xmath222 similarly the @xmath160 summation is done using the phase space constraint in eq .
( [ eq : fase - rum - restriktion - k3 ] ) @xmath223 we see that since @xmath143 and @xmath160 are restricted to the fermi level , we can insert @xmath224 and @xmath225 in the rest of the integrand .
to do the @xmath169 summation , we use the energy conservation contained in the scattering rate .
it is rewritten as ( inserting @xmath226 and @xmath227 ) : @xmath228 .
\label{eq : energy - conservation}\end{aligned}\ ] ]    we have now done the summation over @xmath143 , @xmath145 and @xmath160 and are left with the summation over @xmath161 and @xmath162 of the scattering rate , some fermi functions and the phase factors described above . to this end , we introduce @xmath229 by inserting @xmath226 and @xmath230 in eq .
( [ eq : v - direct - quadratic - sec4 ] ) @xmath231 for a general symmetrized interaction @xmath154 .
furthermore , we collect the exponential tales of the fermi functions eq .
( [ eq : fermi - function-2 - 2m]-[eq : fermi - function-3 - 3 m ] ) , @xmath232 defining @xmath233 inserting @xmath234 , @xmath145 from the energy conservation eq .
( [ eq : energy - conservation ] ) and @xmath235 and @xmath230 .
therefore we finally get the interacting contribution to the current in eq .
( [ eq : i - int - vigtigt - appendix - b ] ) as : @xmath236.\nonumber\end{aligned}\ ] ] here only the step functions , that restricts the integral are included .
next we introduce the dimensionless integration variables @xmath237 for @xmath213 and the dimensionless functions @xmath238 in the integral : @xmath239 , \label{eq : three - particle - contribution - to - current - midway-2-dimless}\end{aligned}\ ] ] where @xmath240 is the integration area shown in fig .
[ fig : int - area](a ) . note that this expression is valid for a general interaction @xmath140 and that it is not possible to extract a power law in temperature times some integral by defining new integration variables .     for the integral eq .
( [ eq : three - particle - contribution - to - current - midway-2-dimless ] ) to calculate the current due to interactions .
the two boundaries for the integration area _ close to the origin _ stemming from the signs of @xmath145 and @xmath211 are indicated .
( b ) the @xmath241 function , a dimensionless version of @xmath242 eq .
( [ eq : xi - def ] ) , important in the calculation using the method of steepest decent .
[ fig : int - area],scaledwidth=35.0% ]    to proceed , we consider the low temperature limit @xmath243 by using the method of steepest decent : due to the exponential function @xmath244 the maximum of @xmath245 will dominate the integral for @xmath246 , since @xmath247 .
the maxima are @xmath248 and @xmath249 and @xmath245 is shown in fig .
[ fig : int - area](b ) . for a decreasing interaction the area of @xmath250 and @xmath251 dominates even though the integrand is zero for @xmath252 .
therefore we expand around the maximum @xmath253 to get the lowest order result in @xmath254 . in view of the integration region fig .
[ fig : int - area](a ) , we use @xmath255 in the integral eq.([eq : three - particle - contribution - to - current - midway-2-dimless ] ) and thereby do the @xmath256 integral using the approximate limits seen in fig .
[ fig : int - area](a ) , i.e. @xmath257 to model the symmetrized potential @xmath140 for small @xmath132 , we include the deviation from a constant as described in eq .
( [ eq : vq - taylor - series ] ) .
this gives @xmath258 to lowest order in @xmath259 . in the exponential we keep @xmath260 to lowest order in @xmath259 ,
i.e. @xmath261 so using the lowest order in @xmath259 in the integrand ( leading to lowest order in @xmath5 ) the interacting contribution to the current is : @xmath262\\ \simeq & \frac{8505}{2048\pi^4}e^{-\frac{{t_{\textrm{f}}}}{t } } \frac{e}{{\hslash}}\frac{(v_0{k_{\textrm{f}}^{}})^4}{{\varepsilon}_{\textrm{f}}^3}\nonumber \\ & \times ( l{k_{\textrm{f}}^ { } } ) \bigg(\frac{{k_{\textrm{f}}^{}}}{q_0}\bigg)^4 \bigg(\frac{t}{{t_{\textrm{f}}}}\bigg)^7 \bigg[\frac{\delta t}{t}+\frac{ev}{{\varepsilon_{\textrm{f}}^ { } } } \bigg]\end{aligned}\ ] ] to lowest order in temperature .
this is the result stated in the text in eq .
( [ eq : i - int - resultat ] ) .
this statement is still true for energy bands like @xmath263 or @xmath264 for @xmath265 .
a special case is linear dispersion where two - particle collisions _ can _ change the distribution .
however , even for this case the two - particle collisions _ can not _ change the number of left- and right movers and hence not the current , as shown in appendix a.                        see _ e.g. _ , m.a .
cazalilla , physical review a * 67 * , 053606 ( 2003 ) .
the interaction potential between fermions is written there in a slightly different form , which however can be reduced ( m. khodas , a. kamenev , l.i .
glazman , unpublished ) to the one we present in the text .
we assume that electrons in the channel produce some image charges in a distant gate .
as the result , the coulomb interaction is cut off at some distance @xmath266 large compared to the fermi wavelength .
note here that the sign of the conductance correction is negative even though there is a positive sign in eq.([eq : i - int - resultat ] ) .
this is due to the way the bias is defined the left contact having the higher chemical potential .      surprisingly , mott formula still yields a reasonably good approximation for thermopower even if the temperature is comparable to the energy scale over which the transmission amplitude changes value , see a.  m. lunde and k. flensberg , j. phys . : condens .
matter * 17 * , 3879 ( 2005 ) .


woven graph codes can be considered as a generalization of low - density parity - check ( ldpc ) block codes @xcite .
their structure as graph codes makes them suitable for iterative decoding .
moreover , the ldpc block codes are known as codes with low - complexity decoding and they can be considered as competitors to the turbo codes @xcite which are sometimes called parallel concatenated codes . as mentioned in @xcite ,
the underlying graph defines a permutation of the information symbols which resembles the interleaving in turbo coding schemes .    on the other hand , similarly to the ldpc codes , graph codes usually have minimum distances essentially smaller than those of the best known linear codes of the same parameters . at a first glance , the minimum distance of a graph code does not play an important role in iterative decoding since the error - correcting capability of this suboptimal procedure is often less than that guaranteed by the minimum distance .
however , in general , the belief - propagation decoding algorithms work better if the _ girth _ of the underlying graph is large , that is , if the minimum distance of the graph code is large @xcite .    in the sequel we distinguish between graph , graph - based , and woven graph codes .
we say that a _ graph code _ is a block code whose parity - check matrix coincides with the incidence matrix of the corresponding graph . _ graph - based codes
_ constitute a class of concatenated codes with constituent block codes concatenated with a graph code ( see , for example , @xcite ) .
each vertex in the underlying graph corresponds to a constituent block code .
the main feature of these codes is that the block length of their constituent block codes coincides with the degree of the underlying graph .
we introduce _ woven graph codes _ which are , in fact , graph - based codes with constituent block codes whose block length is a multiple of the graph degree @xmath7 , that is , their block length is @xmath8 , where @xmath9 is an integer .
in particular , when @xmath9 tends to infinity we obtain convolutional constituent codes .
distance properties of bipartite graph - based codes with constituent block codes were studied in @xcite .
it was shown that if the minimum distance of the constituent block codes is larger than or equal to @xmath10 , then there exist asymptotically good codes with fixed constituent codes among these graph - based codes .
also it was shown in @xcite that for some range of rates , random graph - based codes with block constituent codes satisfy the vg bound when the block length of the constituent codes tends to infinity .
one disadvantage of graph - based codes that becomes apparent in the asymptotic analysis is that good performances can only be achieved when the block length of the constituent block codes ( which in this case coincides with the graph degree @xmath7 ) tends to infinity . in practice
this leads to rather long graph - based codes with not only rather high decoding complexity of the iterative decoding procedures but also high encoding complexity .    in this paper
, we consider a class of the generalized graph - based codes which we call woven graph codes with constituent block and convolutional codes .
they are based on _
@xmath0-partite , @xmath0-uniform hypergraphs_. notice that graph - based codes with constituent block codes based on hypergraphs were considered in @xcite , @xcite .
it is mentioned in @xcite that gallager s ldpc codes are graph codes over hypergraphs .
we consider first woven graph codes with constituent @xmath11 block codes .
a product - type lower bound on the minimum distance of such codes is derived . in order to analyze their asymptotic performances we modify the approach used in @xcite to @xmath0-partite , @xmath0-uniform hypergraphs and constituent @xmath11 block codes .
it is shown that when @xmath9 grows to infinity in the random ensemble of woven graph codes with binary constituent block codes we can find @xmath12 such that there exist codes satisfying the vg lower bound on the minimum distance for any rate .    in order to generalize the asymptotic analysis to woven graph codes with constituent convolutional codes we assume that the binary constituent block code is chosen as a zero - tail ( zt ) terminated convolutional code and consider a sequence of zt convolutional codes of increasing block length @xmath9 .
it is shown that when the overall constraint length of the woven graph code tends to infinity in the random ensemble of such convolutional codes we can find @xmath13 such that there exist codes satisfying the costello lower bound on the free distance for any rate .
we also describe the constituent convolutional codes as block codes over the field of binary laurent series @xcite .
this description as well as the notion of _ block _ hamming distance @xcite of convolutional codes is used to derive a product - type lower bound on the free distance of woven graph codes with constituent convolutional codes and to construct examples of such woven codes with rate @xmath1 . for a given hypergraph the free distance of the woven graph code depends on the numbering of code symbols associating to the hypergraph vertices . by a search over all possible permutations of the constituent code we found an example of a rate @xmath1 woven graph code with overall constraint length @xmath14 and free distance @xmath2 .
the rate @xmath1 woven graph code is based on heawood s bipartite graph @xcite , @xcite and contains constituent convolutional codes with overall constraint length @xmath5 and free distance @xmath15 .
we consider also the encoding problem for graph and woven graph codes .
the traditional encoding technique for graph codes has complexity @xmath16 , where @xmath17 is the blocklength .
we show by examples that some regular block graph codes are quasi - cyclic and thereby can be interpreted as tailbiting ( tb ) codes ( see , for example , @xcite , @xcite ) .
it is known that the encoding complexity of such codes is proportional to the overall constraint length of the parent convolutional code .    by using a tb representation for the graph code we can construct an example of an encoder for a woven graph code that is also represented in the form of a tb code but with overall constraint length less than or equal to @xmath18 , where @xmath19 is the number of constituent convolutional codes with overall constraint length @xmath6 each .    in section
ii , we consider some properties of @xmath0-partite , @xmath0-uniform , @xmath7-regular hypergraphs . we define woven graph codes with constituent block codes as well as with constituent convolutional codes and
obtain product - type lower bounds on their minimum and free distances .
then , in section iii , we derive a lower bound on the free distance of the random ensemble of woven graph codes . in section iv , examples of woven graph codes are given .
we conclude the paper by considering encoding techniques for graph codes and woven graph codes in section v.
a _ hypergraph _ is a generalization of a graph in which the edges are subsets of vertices and may connect ( contain ) any number of vertices .
these edges are called hyperedges .
a hypergraph is called @xmath0-uniform if every hyperedge has cardinality @xmath0 or , in other words , connects @xmath0 vertices . if @xmath20 the hypergraph is simply a graph .
the _ degree of a vertex _ in a hypergraph is the number of hyperedges that are connected to ( contain ) it . if all vertices have the same degree we say that this is the _ degree of the hypergraph_. the hypergraph is _ @xmath7-regular _ if every vertex has the same degree @xmath7 .
let the set @xmath21 of vertices of an @xmath0-uniform hypergraph be partitioned into @xmath22 disjoint subsets @xmath23 , @xmath24 .
a hypergraph is said to be _
@xmath22-partite _ if no edge contains two vertices from the same set @xmath23 , @xmath24 .    in the sequel we consider @xmath0-partite , @xmath0-uniform , @xmath7-regular hypergraphs .
such a hypergraph is a union of @xmath0 disjoint subsets of vertices .
each vertex has no connections in its own set and is connected with @xmath25 vertices in the other subsets .    in fig .
[ figure : graph4 ] a @xmath10-partite , @xmath10-uniform , @xmath26-regular hypergraph is shown .
it contains three sets of vertices .
they are shown by triangles , rectangles , and ovals , respectively .
there are no edges connecting vertices inside any of these three sets .
the vertices are connected by hyperedges each of which connects three vertices .
a _ cycle _ of length @xmath27 in the hypergraph is an alternating sequence of @xmath28 vertices and @xmath27 hyperedges where all vertices are distinct except the initial and the final vertex , which coincide , and all edges are distinct .
the _ girth _ of a hypergraph is the length of its shortest cycle .    in fig .
[ figure : hypergirth ] we show a subgraph that contains the shortest cycle of the @xmath10-partite , @xmath10-uniform , @xmath26-regular hypergraph in fig .
[ figure : graph4 ] .
it consists of the vertices @xmath29 , @xmath30 , and @xmath29 and has girth equal to 2 .
we introduce the notion of a _
compact @xmath31-connected subgraph _ in the hypergraph .
it is a connected subgraph in which each vertex is incident with at least @xmath32 hyperedges .
we call the length ( number of hyperedges ) of the shortest compact subgraph its _ ( @xmath0,@xmath32)-girth_. in fig .
[ figure : hypergirth ] the hyperedges belonging to the shortest @xmath33-compact subgraph are marked by circles . it is easy to see that ( @xmath10,@xmath34)-girth is @xmath35 .
a @xmath34-partite , @xmath34-uniform hypergraph is a bipartite graph .
for such a hypergraph the ( @xmath34,@xmath34)-girth is equal to the girth and a compact subgraph is a cycle .
heawood s bipartite graph @xcite , @xcite with @xmath36 vertices and @xmath37 edges is shown in fig .
[ figure : graphh ] .
this graph contains a set of @xmath3 black and a set of @xmath3 white vertices .
each vertex has no connections within its own set and is connected with @xmath38 vertices from the other set .
the girth of the heawood graph is @xmath35 .      in order to illustrate the structure of a binary graph - based block code with constituent block codes we represent the heawood bipartite graph using a so - called tanner graph @xcite as shown in fig .
[ figure : graph ] .
we introduce a set of @xmath39 ( variable ) vertices which correspond to the code symbols .
each of the @xmath40 ( constraint ) vertices on the right- and left - hand sides corresponds to one of 14 parity checks .
the @xmath38 edges leaving one constraint vertex correspond to a codeword of the constituent @xmath41 block code of rate @xmath42 .
the parity - check matrix of the corresponding graph - based code with binary constituent block codes is @xmath43 where the parity - check matrix @xmath44 of size @xmath45 has the form @xmath46 where @xmath47 is a size @xmath48 parity - check matrix of the constituent block code , and @xmath49 is a size @xmath45 parity - check matrix which is the permutation of the columns of @xmath44 determined by the graph .
notice that in general by choosing @xmath50 and assigning constituent block codes of different rates @xmath42 to the same graph we can obtain graph - based codes of different rates . in general , since in an @xmath0-partite , @xmath0-uniform , @xmath7-regular hypergraph the total number of parity checks is equal to @xmath51 , the code rate @xmath52 of the graph - based code is @xmath53 with equality if and only if all parity - checks are linearly independent . if @xmath20 , then we get @xmath54 .    the simplest example of a heawood graph - based code can be obtained by choosing as constituent block codes a single - parity - check code of rate @xmath55 .
then the parity - check matrix @xmath47 has the form @xmath56 and the parity - check matrix of the graph - based code is @xmath57 in this case the graph - based code coincides with the graph code since ( [ incdmatr ] ) is the incidence matrix of the heawood graph . in @xcite
it is proved that the minimum distance of the bipartite graph - based code with single - parity - check constituent codes is @xmath58 , where @xmath59 is the girth of the corresponding graph .
notice that for the tanner graph we have @xmath60 .
the parity - check matrix ( [ incdmatr ] ) is a @xmath61 parity - check matrix . taking into account
that one check is linearly dependent on the other , we obtain a @xmath62 binary block code .
its minimum distance is @xmath63 .
consider the hypergraph shown in fig .
[ figure : graph4 ] .
its incidence matrix has the form @xmath64 and is a @xmath65 parity - check matrix of a hypergraph - based code which coincides with the parity - check matrix of the hypergraph code .
each column represents a hyperedge and each row represents a vertex of this hypergraph . for example , the first four rows represent the vertices 1 , 2 , 3 , and 4 ( triangles ) , the next four rows he vertices 5 , 6 , 7 , and 8 ( rectangles ) , and the last four rows the vertices 9 , 10 , 11 , and 12 ( ovals ) .
the first column represents the hyperedge which connects the vertices 1 , 5 , and 9 , the second column the hyperedge connecting vertices 1 , 6 , and 10 _ etc_. the rows of ( [ hyper ] ) are linearly dependent . by removing two parity checks
we obtain a @xmath66 linear block code with the minimum distance @xmath67 , where @xmath68 is the ( @xmath10,@xmath34)-girth of the hypergraph .
the rate of this hypergraph code is @xmath69 , which satisfies inequality ( [ rate ] ) ,    @xmath70    the tanner version of this hypergraph is shown in fig .
[ figure : tannerhypergraph ] .    for an @xmath0-partite , @xmath0-uniform ,
@xmath7-regular hypergraph - based code with constituent block codes we have the following theorem .
[ graph - based ] the minimum distance of a hypergraph - based code based on an @xmath0-partite , @xmath0-uniform , @xmath7-regular hypergraph with ( @xmath0,@xmath71)-girth @xmath72 and containing constituent block codes with minimum distance @xmath73 is    @xmath74    .
any nonzero codeword in an @xmath0-partite , @xmath0-uniform , @xmath7-regular hypergraph - based code always corresponds to a connected @xmath75-subgraph or a set of disjoint connected subgraphs .
these subgraphs are called _ active _ @xcite , @xcite .
all hyperedges and vertices in an active subgraph are also called _
the number of hyperedges in the shortest connected subgraph is equal to @xmath72 .
any nonzero symbol in a codeword corresponds to an active hyperedge in the graph . by using the arguments given above
, we conclude that for any codeword @xmath76 , @xmath77 where @xmath78 is the hamming weight of @xmath76 . minimizing over @xmath76 completes the proof .
now assume that the constituent code assigned to the hypergraph vertices is a binary @xmath11 linear block code determined by a parity - check matrix @xmath79 where @xmath80 is a size @xmath81 matrix , @xmath82 is the set of all possible binary matrices of size @xmath81 .
let @xmath83 denote such a binary @xmath11 constituent block code determined by the matrix ( [ nonbin ] ) .
we call the corresponding hypergraph - based code with @xmath83 as constituent codes a _ woven graph _ code with constituent block codes .
consider an example of a woven graph code based on the bipartite graph with girth @xmath84 shown in fig .
[ figure : bgraph ] . the tanner version of this the so - called `` utility '' bipartite graph
is shown in fig .
[ figure : tgraph ] .
the incidence matrix of this graph is @xmath85 we use a constituent ( @xmath86 ) linear block code with @xmath87 determined by the parity - check matrix @xmath88 by searching over all possible permutations of the matrices @xmath89 , @xmath90 , and @xmath91 we found the following parity - check matrix of the woven graph code with the best minimum distance    @xmath92    the matrix ( [ wovbl ] ) describes a @xmath93 linear block code with @xmath94 .
any codeword @xmath95 of the @xmath11 constituent block code can be represented as a sequence of @xmath7 blocks of length @xmath9 , that is , @xmath96 , where @xmath97 , @xmath98 .
we define the _ minimum hamming block distance _ between the codewords @xmath95 and @xmath99 of the constituent block code as @xmath100 where @xmath101 , @xmath98 .
next we will prove the following theorem .
[ block ] the minimum distance of woven graph codes based on @xmath0-partite , @xmath0-uniform , @xmath7-regular hypergraphs with ( @xmath0,@xmath102)-girth @xmath103 and containing constituent block codes with minimum distance @xmath71 and minimum block distance @xmath104 can be lower - bounded by @xmath105    _ proof_. any nonzero codeword corresponds to an active connected subgraph or a set of disjoint connected subgraphs and the number of hyperedges in the shortest subgraph is @xmath103 .
any nonzero symbol in a codeword activates a hyperedge in the graph , that is , not less than @xmath0 constituent subcodes correspond to a codeword . since at most @xmath7 hyperedges
are connected with any hypergraph vertex then the number of active constituent subcodes can be lower - bounded by @xmath106 taking into account that any codeword of block weight greater than or equal to @xmath102 in the constituent block code has a weight at least equal to @xmath71 we obtain the following inequality @xmath107 for any codeword @xmath76 and the proof is complete .    from theorem [ block ] for the woven graph code determined by ( [ wovbl ] )
we obtain that @xmath108      woven graph codes with constituent convolutional codes can be considered as a straightforward generalization of the woven graph code with constituent block codes .
assume that the @xmath83 code is chosen as a zero - tail terminated ( zt ) convolutional code and consider a sequence of zt convolutional codes with increasing @xmath9 .
it is evident that when @xmath9 tends to infinity the @xmath11 constituent code @xmath83 can be chosen as a rate @xmath42 binary convolutional code with constraint length @xmath6 .
then the corresponding woven graph code has rate @xmath109 and its constraint length is at most @xmath110 .
another description of woven graph codes with constituent convolutional codes follows from the representation of the constituent convolutional code in polynomial form .
let @xmath111 be a minimal encoding matrix @xcite of a rate @xmath42 , memory @xmath112 convolutional code , given in polynomial form , that is , @xmath113 where @xmath114 , @xmath115 , @xmath116 , are binary polynomials such that @xmath117 .
the overall constraint length is @xmath118 .
the binary information sequence @xmath119 is encoded as @xmath120 where @xmath121 is a binary code sequence .
let @xmath122 denote a parity - check matrix for the same code , @xmath123 where @xmath124 is the redundancy of the constituent code .
we denote by @xmath125 the field of binary laurent series and regard a rate @xmath42 constituent convolutional code as a rate @xmath42 block code @xmath126 over the field of binary laurent series encoded by @xmath111 .
then its codewords @xmath127 are elements of @xmath128 , which is the @xmath7-dimensional vector space over the field of binary laurent series @xcite .
the minimum hamming _ block _ distance between the codewords @xmath129 and @xmath130 is defined @xcite as    @xmath131    where @xmath132 is the hamming ( block ) weight of @xmath133 .    representing a convolutional code as a block code over the field of binary laurent series
we can obtain a woven graph code with constituent convolutional codes as a generalization of a graph - based code with binary constituent block codes .
for example , a parity - check matrix @xmath134 of the rate @xmath135 heawood s graph - based code with @xmath4 constituent convolutional codes has the form @xmath136 where @xmath137 and @xmath138 are short - hand for @xmath139 and @xmath140 , respectively , and @xmath141 is a parity - check matrix of the rate @xmath4 constituent convolutional code and @xmath142 is one of six possible permutations of @xmath143 .    exploiting the above definitions we can interpret this bipartite woven graph - based code with constituent convolutional codes as follows .
the left column of vertices in fig .
[ figure : graph ] represents @xmath19 parity checks each of which determines one of @xmath19 constituent fixed and identical convolutional codes and their @xmath144 branches represent the elements @xmath145 , @xmath146 even , @xmath147 , @xmath148 . similarly , the right column of vertices represents same convolutional codes and their @xmath144 branches represent the elements @xmath149 , @xmath146 odd , @xmath150 , @xmath148 , where the set @xmath151 is a random permutation of the set @xmath152 determined by the graph .
we can also regard the @xmath19 left constituent convolutional codes as a warp with @xmath144 threads .
each of the @xmath19 right constituent convolutional codes are tacked on @xmath7 of the threads in the warp such that each thread of the warp is tacked on exactly once .
thus , our construction is a special case of a woven code @xcite and we call this graph - based code a _ woven graph _ code .
[ wovbound ] the free distance of a woven graph code based on an @xmath0-partite , @xmath0-uniform , @xmath7-regular hypergraph with the ( @xmath0,@xmath102)-girth @xmath103 and containing constituent convolutional codes with free distance @xmath153 and minimum block distance @xmath154 can be lower - bounded by @xmath155    . since
woven graph codes with constituent convolutional codes can be considered as a generalization of woven graph codes with constituent block codes , the theorem follows from theorem [ block ] when @xmath9 tends to infinity .    for a woven graph code based on a bipartite graph with girth @xmath59 and containing constituent convolutional codes with minimum block distance @xmath156 and free distance @xmath153 by a straightforward generalization of the approach of @xcite we obtain the following tighter bound on the free distance @xmath157
we will show that the ensemble of random woven graph codes based on random @xmath0-partite , @xmath0-uniform , @xmath7-regular hypergraphs with a fixed degree @xmath7 and with a fixed number of vertices @xmath19 in each subgraph contains asymptotically good codes . in order to prove this
we will modify the approach in @xcite .
first we consider the ensemble of random woven graph codes with rate @xmath42 constituent block codes determined by the edges of a random @xmath0-partite , @xmath0-uniform , @xmath7-regular hypergraph corresponding to the time - varying random parity - check matrix @xmath158 where @xmath159 , @xmath160 , is a block matrix of size @xmath161 ( or a binary matrix of size @xmath162 ) and @xmath163 denotes a random permutation of the columns of @xmath164 , @xmath165 where @xmath166 , @xmath167 , denotes the random parity - check matrix ( [ nonbin ] ) which determines the @xmath11 constituent block code and @xmath19 is the number of constituent codes in each subgraph .
_ remark _ : in @xcite a more restricted ensemble of random codes is studied in which all matrices are identical random matrices . in the proof of theorem 1
we need that the syndrome components are independent random variables in the product probability space of random matrices and random permutations .
the following simple example shows that this is not always the case if all matrices are identical .
consider @xmath168 constituent block codes of block length @xmath169 with @xmath170 information symbols .
this example is rather artificial since the rate of the constituent block code @xmath171 and therefore the rate of the graph - based code with @xmath20 is @xmath172 . in this case
the parity - check matrix of the code has the form @xmath173 where @xmath174 is a random permutation of @xmath7 elements .
first assume that all matrices are identical , that is , @xmath175 .
there are only 8 equiprobable elements in the product space , namely , @xmath176 for any vector @xmath177 of weight 1 we have the following set of random equiprobable syndromes : @xmath178 therefore , @xmath179 if @xmath180 and @xmath181 are both random and independent this probability is equal to 1/4 .
although this remark contradicts the proof of theorem 3 in @xcite , there exists another ( combinatorial ) way to prove the same statement for identical @xmath182 @xcite .
next we prove the following theorem .
[ vgtheorem ] ( varshamov - gilbert lower bound ) for any @xmath183 , some @xmath184 , some integer @xmath185 and for all @xmath186 in the random ensemble of length @xmath187 woven graph codes with @xmath11 binary block constituent codes of rate @xmath42 there exist codes of rate @xmath188 such that their relative minimum distance @xmath189 satisfies the inequalities @xmath190    where @xmath191 is a root of the equation @xmath192 and @xmath193 is the solution of @xmath194 , and @xmath195 denotes the binary entropy function .
_ let @xmath196 be the hamming weight of the codeword @xmath76 of the random binary woven graph code @xmath197 .
we are going to find a parameter @xmath32 such that the probability @xmath198 tends to 0 for all @xmath199 .
we can rewrite @xmath200 as @xmath201 where @xmath202 and @xmath203 denotes the number of nonzero constituent codewords in the @xmath146th subgraph corresponding to the codeword of weight @xmath196 .    in the ensemble of random parity - check matrices
@xmath204 , @xmath205 , of size @xmath206 the probability that a nonzero vector @xmath207 is a codeword of the corresponding constituent random binary code @xmath83 is equal to @xmath208 since the syndromes of the constituent codes are equiprobable sequences of length @xmath209 . taking into account that in the @xmath146th subgraph we have @xmath203 nonzero constituent codewords the probability @xmath210 can be upper - bounded by @xmath211    in order to estimate the probability @xmath212 we prove the following lemma .    for the ensemble of binary woven graph codes with constituent block codes described in theorem [ vgbound ] ,
the probability @xmath212 that a codeword of weight @xmath196 contains @xmath202 nonzero constituent codewords in the @xmath0 subgraphs can be upper - bounded by @xmath213 [ lemma ]    _ proof .
_ taking into account that in the @xmath146th subgraph the number of nonzero component codewords is equal to @xmath203 and that the subgraphs are random and independent we can rewrite the probability @xmath214 as @xmath215 the probability @xmath216 can be upper - bounded as @xmath217 where @xmath218
. the cardinality of @xmath219 can be upper - bounded as @xmath220 @xmath221 where the sum is upper - bounded by the maximal term times the number of terms @xmath222 .
notice that in the above derivations we ignored the fact that @xmath223 can be noninteger since we consider the asymptotic behaviour of ( [ main ] ) .
it follows from lemma 1 that @xmath224 @xmath225 @xmath226 @xmath227    consider the asymptotic behaviour of ( [ main ] ) when @xmath228 tends to infinity .
introduce the notations @xmath229 and @xmath230 and the function @xmath231 after simple derivations we obtain @xmath232 } ( 1-s)h(\delta)-(1-r_{\rm wg})\gamma+s\gamma h\left(\frac{\delta}{\gamma}\right ) \right\ } \label{asympt}\ ] ] where @xmath188 is the rate of binary woven graph code . maximizing ( [ asympt ] ) over @xmath233
gives @xmath234 inserting @xmath235 and @xmath236 into ( [ asympt ] ) we obtain @xmath237 which coincides with ( 9 ) and ( 10 ) in @xcite for @xmath20 , that is , if the graph is bipartite .
for any @xmath238 and @xmath239 from @xmath240 , it follows that there exist codes of rate @xmath238 with relative minimum distance @xmath241 .
let @xmath191 denote the solution of the equation @xmath242 for @xmath243 and let @xmath244 be the solution of @xmath194 . solving ( [ fun ] ) for @xmath235 and @xmath236
we obtain that there exist woven graph codes of rate @xmath238 with the relative minimum distance @xmath245 satisfying the inequalities : @xmath246    in fig .
[ figure : rd ] the lower bound ( [ vgbound ] ) on the relative minimum distance for the ensemble of binary woven graph codes with block constituent codes as a function of the code rate is shown .
it is easy to see that when @xmath0 grows the ensemble of binary woven graph codes contains codes meeting the vg bound for almost all rates @xmath247 .
[ figure : rd2 ] demonstrates the gap @xmath248 between the vg bound and the code rate as a function of the relative minimum distance @xmath245 for different values of @xmath0 .
it follows from fig .
[ figure : rd2 ] that for @xmath249 the difference in code rate compared to the vg bound is negligible .
consider a zt convolutional woven graph code with constituent zt convolutional codes of rate @xmath42 .
the length of a zt woven graph codeword in @xmath144-tuples is equal to @xmath250 where @xmath9 is the number of @xmath144-tuples influenced by information symbols and @xmath251 is the memory of the woven graph code of rate @xmath188 .
denote by @xmath252 the free distance of the corresponding woven graph code .
now we can prove the following    ( costello lower bound ) for any @xmath183 , some @xmath253 , some integer @xmath13 , and for all @xmath254 in the random ensemble of rate @xmath188 woven graph codes over @xmath0-partite , @xmath0-uniform , @xmath7-regular hypergraphs with constituent convolutional codes of rate @xmath42 there exists a code with memory @xmath251 such that its relative free distance @xmath255 satisfies the costello lower bound @xcite , @xmath256    _ proof _ : analogously to the derivations in the proof of theorem [ vgtheorem ] let @xmath202 where @xmath203 denotes the number of nonzero constituent codewords in the @xmath146th subgraph corresponding to the codeword of weight @xmath196 , @xmath257 . in order to evaluate the number of nonzero constituent codewords among the @xmath19 constituent codewords , notice that the set of such codewords is a union of sets of nonzero constituent codewords belonging to each of the @xmath0 subgraphs .
the cardinality of the union is at least @xmath258 .
therefore the all - zero `` tail '' required to force the encoder into the zero state has length at least @xmath259 .
the total number of redundant symbols consists of two parts : the number @xmath260 of parity - check symbols for the nonzero constituent codewords in the @xmath0 subgraphs and at least @xmath259 redundant symbols required for zero - tail terminating of the woven graph code .
thus , formula ( [ redunds ] ) can be rewritten as    @xmath261    the statement of the lemma [ lemma ] is changed in a following way @xmath262    instead of ( [ mainf ] ) we now have @xmath263    by introducing the notations @xmath264 we obtain from ( [ redunds_conv])([mainf2 ] ) that @xmath265 } \left\ { ( 1-s)(1+\mu)h\left(\frac{\delta}{1+\mu}\right)-\gamma(1+\mu-\mu r_{\rm wg})+\gamma ( 1+\mu ) sh\left(\frac{\delta}{\gamma(1+\mu)}\right ) \right\}. \label{mainfconv}\end{aligned}\ ] ] maximizing ( [ mainfconv ] ) over @xmath233 , we obtain @xmath266 where @xmath267    if @xmath0 is large enough , then @xmath236 .
it follows from ( [ mainfconv ] ) that @xmath268 maximization of @xmath269 over @xmath270 gives @xmath271 where @xmath272 we can find a bound on @xmath273 by solving @xmath274 .
thus , we can conclude that for any @xmath275 we can find a woven graph code such that ( [ costellobound ] ) holds
we start with considering a graph code determined by the parity - check matrix ( [ incdmatr ] ) . as mentioned before
, the matrix ( [ incdmatr ] ) can be considered as a @xmath276 parity - check matrix . since the parity checks defined by the graph are linearly dependent ( the sum of the rows of ( [ incdmatr ] ) is equal to zero ) it turned out that by ignoring one parity check we obtain a parity - check matrix of a @xmath62 linear block code .
for simplicity we consider the rate @xmath277 code that is obtained by ignoring the eighth information symbol which yields a @xmath278 subcode of this code .
it is easy to see that renumbering the graph vertices by adding to each vertex number some fixed number modulo the total number of vertices preserves both the incidence and adjacency matrices of the graph .
for example , in fig .
[ figure : graphh ] , by adding 2 modulo 14 we will get exactly the same graph . when we have a similar property for linear codes we call such codes quasi - cyclic codes and these block codes can be described as tailbiting ( tb ) convolutional codes .
renumbering the vertices corresponds to permuting the rows of ( [ incdmatr ] ) . by row permutations , ( [ incdmatr ] )
can be reduced to the form @xmath279 it follows from ( [ modifmatr ] ) that the graph shown in fig .
[ figure : graphh ] corresponds to a @xmath278 tb code with a parent convolutional code determined by the parity - check matrix @xmath280 it means that this code is `` tail - bitten '' at the 7th level of the trellis diagram .
a corresponding polynomial generator matrix of the parent convolutional code has the form @xmath281    the minimum distance of the @xmath278 tb code is equal to the graph girth , that is , @xmath63 .
notice that many regular bipartite graphs look very similar to the heawood graph in the sense that by manipulating the incidence ( parity - check ) matrices and truncating lengths we can obtain infinite families of graphs .
some properties of these graphs can be easily predicted from the properties of the corresponding parent convolutional codes .
consider the parity - check matrix ( [ woven ] ) of the woven graph code based on the heawood bipartite graph with constituent convolutional codes of rate @xmath4 .
this woven graph code has the rate @xmath282
.    let the rate @xmath4 constituent convolutional code of memory @xmath283 and overall constraint length @xmath5 with @xmath15 be given by the generator matrix @xmath284 a corresponding parity - check matrix @xmath122 is @xmath285 notice that the constituent code @xmath286 considered as a block code over @xmath287 represents a @xmath288 block code with the minimum distance @xmath289 .    by using the product - type lower bound ( [ bipartbound ] )
we obtain @xmath290 on the other hand , it was verified by computer search that any codeword of the woven graph code determined by ( [ woven ] ) consists of at least three nonzero codewords of the component code @xmath126 described by ( [ 90 ] ) .
moreover , it was found by computer search that each of these nonzero codewords of @xmath126 has the minimum block weight @xmath291 .
note that the codewords of the block code over @xmath287 with block weight @xmath291 corresponds to the codewords of the convolutional code belonging to its subcodes of rate @xmath292 .
these three subcodes have generator matrices @xmath293 @xmath294 @xmath295 where @xmath296 , @xmath297 , and @xmath298 .
the minimum free distance over all these subcodes of rate @xmath292 is equal to 8 .
taking into account that all other codewords of the woven graph code contain at least four nonzero codewords of @xmath299 of block weight @xmath300 we obtain an improved lower bound on the free distance of the woven graph code as @xmath301 .    in order to obtain an upper bound on the free distance of the woven graph code we consider the parity - check matrix ( [ woven ] ) in more detail .
it also describes a quasi - cyclic code and can by row permutations be reduced to a parity - check matrix of a two - dimensional code , a tb ( block ) code in one dimension and a convolutional code in the other , @xmath302 a parity - check matrix of the parent convolutional code for the tb code ( [ modwoven ] ) given in symbolic form is @xmath303 where @xmath304 and @xmath305 are formal variables .
the matrix ( [ 97 ] ) can be considered as a parity - check matrix of a two - dimensional convolutional code .
the variable @xmath304 corresponds to the parent convolutional code of the heawood graph code ( [ convmatr ] ) , the variable @xmath305 is used for the constituent convolutional code ( [ 90 ] ) .    a generator matrix of the two - dimensional convolutional code with the parity - check matrix ( [ 97 ] )
has the form @xmath306 where @xmath307 , @xmath308 , @xmath309 , @xmath310 , @xmath311 , and @xmath312 .    the generator matrix ( [ 97gen ] ) tail - bitten over variable @xmath304 at length @xmath37 yields the generator matrix @xmath313 of the code ( [ modwoven ] ) , @xmath314 where @xmath315 is short - hand for @xmath316 .
notice that any of the six permutations of the columns @xmath139 , @xmath317 generates a woven graph code .
the permutation @xmath318 , @xmath319 , and @xmath320 describes the woven graph code with the largest free distance .
the overall constraint length of this generator matrix is equal to 70 but the matrix is not in minimal form .
a minimal - basic generator matrix @xcite has the overall constraint length equal to 64 and differs from ( [ tb13 ] ) by one row which can replace any of the rows of @xmath313 and has the form @xmath321 where @xmath322 where @xmath323 and @xmath324 .
the matrix ( [ tb13 ] ) is a generator matrix of a convolutional code of rate @xmath325 . by applying the beast algorithm @xcite to the minimal - basic generator matrices corresponding to the different permutations of the columns @xmath139 , @xmath317 we obtained the free distance and a few spectrum coefficients of the corresponding woven graph codes .
the parameters of the best obtained woven graph codes are presented in table 1 .
generally speaking , encoding of graph - based block codes has complexity @xmath16 , where @xmath17 is the blocklength .
this technique implies that we find a generator matrix corresponding to the given parity - check matrix and then multiply the information sequence by the obtained generator matrix .
however , we showed by examples that some regular graph codes as well as woven graph codes are quasi - cyclic codes and thereby they can be interpreted as tb codes . for this class of codes
the complexity of the encoding is proportional to the constraint length of the parent convolutional code .    in this section
we are going to illustrate by an example an encoder of a woven graph code with constituent convolutional codes having encoding complexity proportional to the overall constraint length of the corresponding woven graph code @xmath326 .
consider again the woven graph code in our example .
it is based on the heawood graph and uses constituent convolutional codes of rate @xmath4 and overall constraint length @xmath5 . taking into account the representation ( [ tb13 ] ) of the woven graph code as a rate @xmath325 two - dimensional code , a tb ( block ) code in one dimension and a convolutional code in the other , we can draw its encoder as shown in fig .
[ figure : encoder ] .    the input symbols @xmath327 enter the encoder once per each cycle of duration seven time instants . at each time moment
the contents of each register is rewritten into the next ( modulo seven ) register and the three output symbols @xmath328 are generated . in other words , each of the registers corresponding to the constituent code can be considered as an enlarged delay element of the encoder of the `` tb - dimension '' code determined by the graph .
the sequence @xmath327 determines a transition between the states of this encoder . after a cycle of seven time instants we return to the starting state of the enlarged encoder and a tb - codeword ( or a word from one of its cosets ) of length 21 has been generated .
then the following seven input symbols @xmath329 enter and after seven time instants another word of length 21 has been generated , etc .
the asymptotic behavior of the woven graph codes with block as well as with convolutional constituent codes has been studied .
it was shown that in the random ensemble of such codes based on @xmath0-partite , @xmath0-uniform , @xmath7-regular hypergraphs we can find a value @xmath13 such that for any code rate there exist codes meeting the vg and the costello lower bound on the minimum distance and free distance , respectively .
product - type lower bounds on the minimum distance of graph - based and woven graph codes have been derived .
example of a rate @xmath1 woven graph code with free distance above the product bound is presented .
it is shown , by an example , that woven graph codes can be encoded with a complexity proportional to the constraint length of the constituent convolutional code .
this work was supported in part by the royal swedish academy of sciences in cooperation with the russian academy of sciences and in part by the swedish research council under grant 621 - 2007 - 6281 .
r. g. gallager , _ low - density parity - check codes _ , cambridge , ma : mit press , 1963 . c. berrou , a. glavieux , and p. thitimajshima , `` near shannon limit error - correcting coding and decoding : turbo codes , '' _ in proc .
( icc ) _ , geneva , switzerland , pp . 10641070 , 1993 .
a.  barg and g.  zemor , `` distance properties of expander codes , '' _ ieee trans .
inf . theory _
52 , no.1 , pp .
7890 , jan . 2006 .
hu and m. f.  fossorier `` on the computation of the minimum distance of low - density parity - check codes , '' _ int .
conference on communications ( icc04 ) _ , paris , june 2004 .
g. schmidt , v.  v. zyablov , and m. bossert , `` on expander codes based on hypergraphs , '' _ isit 2003 _ , yokohama , japan , june 29july 4 , 2003 .
s. hoory and y. bilu , `` the codes from hypergraphs , '' _
european journal of combinatorics _ vol .
25 , issue 3 , pp . 339354 , april , 2004 .
a. barg , private communication , may 2007 .
r. johannesson and k. sh .
zigangirov , _ fundamentals of convolutional coding _ , ieee press , piscataway , nj , 1999 .
m. bossert , c. medina , and v.  sidorenko , `` encoding and distance estimation of product convolutional codes , '' _ in proc .
information theory ( isit05 ) _ , adelaide , australia , sept . 4 - 9 , pp .
10631067 , 2005 .        j. h. ma and j. k. wolf , `` on tail - biting convolutional codes , '' _ ieee trans .
com-34 , pp .
104111 , 1986 .
r. m.  tanner , `` minimum - distance bounds by graph analysis , '' _ ieee trans . inf . theory _ ,
2 , pp . 808820 , feb .
r. m.  tanner , `` a recursive approach to low complexity codes , '' _ ieee trans .
inf . theory _ ,
5 , pp . 533547 , sep . 1981 .
l. bassalygo , private communication , may 2007 .
s.  h@xmath330st , r. johannesson , and v.  v. zyablov , `` woven convolutional codes i : encoder properties , '' _ ieee trans .
inf . theory _
1 , pp . 149161 , jan .
r. m.  tanner , d. sridhara , a. sridharan , t. e. fuja , and d. j.  costello , jr .
, `` ldpc block and convolutional codes based on circulant matrices , '' _ ieee trans .
inf . theory _ ,
29662984 , dec . 2004 .
i. e. bocharova , m. handlery , r. johannesson , and b.d .
kudryashov , `` a beast for prowling in trees , '' _ ieee trans .
inf . theory _ ,
50 , no . 6 , pp .
12951302 , june 2004 .