Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: 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 .
Abstract:

Answer: constructions of woven graph codes based on constituent block and convolutional codes are studied . it is shown that within the random ensemble of such codes based on @xmath0-partite , @xmath0-uniform hypergraphs , where @xmath0 depends only on the code rate , there exist codes satisfying the varshamov - gilbert ( vg ) and the costello lower bound on the minimum distance and the free distance , respectively . a connection between regular bipartite graphs and tailbiting codes is shown . some examples of woven graph codes are presented . among them an example of a rate @xmath1 woven graph code with @xmath2 based on heawood s bipartite graph and containing @xmath3 constituent rate @xmath4 convolutional codes with overall constraint lengths @xmath5 is given . an encoding procedure for woven graph codes with complexity proportional to the number of constituent codes and their overall constraint length @xmath6 is presented .    * index terms*convolutional codes , girth , graphs , graph codes , hypergraphs , ldpc codes , tailbiting codes , woven codes .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: vesicle dynamics has long been investigated both experimentally and theoretically  @xcite . the decay rate of the thermally excited shape fluctuation provides information of vesicle elasticity . the knowledge of the microscopic relaxation can also be used to predict the macroscopic rheological property of the vesicle solution @xcite . the vesicles are mostly non - equilibrium system . the external perturbation , be it thermal , electrical , sonication  @xcite , or flow  @xcite , transforms the lamellar structure into the vesicles . in such situation , one often obtains mixture of uni - lamellar vesicle and multi - lamellar vesicles ( mlv ) of various sizes . however , in sharp contrast with the very detailed calculations and scattering experiments on the unilamellar vesicles , there are relatively few studies of the mlv dynamics in the literature . in this paper , we consider the two - bilayer vesicle which is the simplest mlv . we calculate the dynamic relaxation rates , as well as the rheological response of the two - layer - vesicles . from the experimental point of view , it is hard to prepare vesicles with exactly two bilayers . nonetheless , such a concrete calculation provides a clear picture of how the interaction between the two membranes affect the relaxation rates , as well as the squeezing ( lubrication ) flow which arises exclusively in mlv . since our method can be extend to vesicles with more than two bilayers , a more specific calculation can be performed similarly . the squeezing dynamics of the sandwiched solvent between the two layers have been analyzed in two related problems . this relaxation process of the flat lamellar system was calculated by brochard and de gennes as the `` slip mode ''  @xcite , or later measured and called as `` baroclinic mode ''  @xcite . in the soap film system , the `` squeezing mode '' dispersion has been calculated and measured @xcite . in this work , we obtain the similar relaxation in which the solvent is squeezed between the two adjacent bilayers to relax the bilayer curvature energy and the mutual interaction energy between the two bilayers . as this mode arises only in mlv , we are interested in its dispersion relation and its coupling strength with the applied shear . due to the strong lubrication resistance , the squeezing mode is often the slowest relaxation mode . this makes the squeezing mode an important candidate for the vesicle rheology . however , not every relaxation mode is equally excited by the applied shear . therefore whether the shear can drive the squeezing mode with a large amplitude is an important problem to consider . in this work , we will try to build our understanding of the coupling strength through the concrete calculation . based on the analysis below , we find that the squeezing mode makes the dominant contribution to the complex viscosity for strongly interacting bilayers . interestingly , when the gap between the two bilayers is smaller than a characteristic crossover gap , the fast undulation mode becomes the dominant mode for the complex viscosity . we find that the crossover gap gets very small when the bilayer interaction is strong . on the opposite limit where the bilayer has extremely weak interaction , the crossover gap becomes comparable to the inner vesicle radius . below in sec . [ sec : model ] , we define our model . section [ sec : elasticforce ] discusses the elastic force . in sec . [ sec : normalforce ] and appendices , the flow resistance between the bilayers is solved . in sec . [ sec : spectrum ] , we shall analyze the relaxation rate . in sec . [ sec : viscoelasticity ] , we consider the viscoelasticity of the dilute vesicle suspension . finally in sec . [ sec : summary ] , we summarize and discuss our results . we consider a two - layer - vesicle with two bilayers located at the spheres with mean radius @xmath0 and @xmath1 as shown in fig .  [ fig1 ] . the deformation of the bilayer shapes are described by the radial layer displacements @xmath2 and @xmath3 relative to the two reference spheres respectively . both of the bilayer membranes are surrounded by a solvent of viscosity @xmath4 . let @xmath5 denote the area per molecule projected on the reference sphere , and @xmath6 the averaged projected area per molecule . for each membrane , we define the dimensionless surface density @xmath7 ( @xmath8 ) given by the ratio @xmath9 . notice that @xmath7 becomes unity at equilibrium . in a fluctuating vesicle , @xmath7 is not uniform in general . in this work we propose the free energy of a two - layer - vesicle as @xmath10 { \rm d } a_n ,   \label{f_cont}\end{aligned}\ ] ] where @xmath11 is the layer compression modulus , @xmath12 the surface tension , @xmath13 the bending modulus , and @xmath14 the area stretching / compression modulus . here both @xmath13 and @xmath14 are taken to be the same for the two membranes . moreover , @xmath15 is the differential solid angle , and the surface area element is approximately given by @xmath16 r_n^2 \ , { \rm d } \omega$ ] . the mean curvature @xmath17 is given by @xmath18 up to linear order in @xmath19 . two comments should be made about the first term of eq . ( [ f_cont ] ) . the compression modulus @xmath11 should be a function of @xmath0 and @xmath1 . in principle a microscopic statistical model for mlv should provide the functional form . here we focus on the layer dynamics , therefore to build such a model is beyond the scope of this work . as a rough estimate for discussion , below in subsection [ sec : interaction ] , we will use the @xmath11 , which is calculated from the flat layers , and replace the distance @xmath20 between two layers by @xmath21 . this approximation is justified for @xmath22 , while it may deviate considerably when @xmath20 is comparable to @xmath1 . the interaction term proposed here is proportional to @xmath23 , which arises naturally at small @xmath20 from the square of the strain . at large @xmath20 , the curvatures and the areas for the two bilayers are very different . presumably a microscopic model may derive a more suitable weighted interaction energy , which is proportional to @xmath24 ^ 2 $ ] with a weighting exponent @xmath25 . below we present the calculation without this extra weighting factor . nonetheless , the calculation procedure is exactly the same for the case @xmath26 .    at the dilute phase boundary of the lamellar phase , the surface tension @xmath12 should vanish . inside the lamellar phase , the tension depends on the applied osmotic pressure . to simplify the discussion , we shall only consider the case where the surface tension vanishes @xmath27 . in fact , the calculation with finite surface tension can be carried out in the same way . the layer displacement @xmath19 is related to the radial velocity @xmath28 at @xmath29 by @xmath30 where we neglect the solvent permeation . the surface density obeys @xmath31 , \label{phi}\ ] ] where @xmath32 is the two - dimensional ( 2d ) surface derivative and @xmath33 is the tangential velocity at @xmath29 . here we have neglected the amphiphile exchange flux between the neighboring bilayer ( the amphiphile permeation ) .     and @xmath1 . the surrounding solvent viscosity is @xmath4 everywhere . ] the flow field obeys the stokes equation which is presented here as the form @xmath34 where @xmath4 is the viscosity , and @xmath35 the pressure . here we write down the force balance conditions that are satisfied on the two layers . the normal force balance is given by @xmath36 where @xmath37 ( @xmath38 being the radial distance ) , and the superscripts @xmath39 and @xmath40 indicate that the stresses are evaluated at the exterior and interior of the bilayer , respectively . on the other hand , the tangential force balance is given by @xmath41 where @xmath42 is the 2d tangent displacement of the layers , and the 2d stress is defined as @xmath43 ( @xmath44 and @xmath45 being the polar and azimuthal angles , respectively ) . the viscoelastic response of a dilute two - layer - vesicle suspension can be calculated by considering the stress response to an external flow at large distances from the vesicles @xmath46 e^{i\omega t } , \label{v_shear}\ ] ] where @xmath47 is the strength of the elongational flow , @xmath48 are the spherical harmonics , and @xmath49 is the angular frequency . each suspended vesicle contributes to the averaged stress . in the dilute solution , the effective complex viscosity is calculated as  @xcite @xmath50 where @xmath51 is the volume fraction occupied by the vesicles ( @xmath52 is the number density of the vesicles ) , and @xmath53 is the frequency dependent complex coefficient of the spherical harmonic expansion of the pressure ( see eqs . ( [ p_expansion ] ) and ( [ pc ] ) later ) . the effective viscosity can also be expressed in terms of the complex modulus as @xmath54 . the elastic forces are calculated by evaluating the derivatives @xmath55 and @xmath56 . from eq . ( [ phi ] ) , we see that the bilayer tangential displacement and the radial displacement produce the first order surface density perturbation as @xmath57 the tangential elastic force is given by @xmath58 where the tension perturbation @xmath59 is induced by the surface density perturbation , i.e. , @xmath60 . notice that the 2d gradient has the usual component form ; @xmath61    up to linear order in @xmath19 and @xmath59 , the normal force on bilayer 1 is given by @xmath62 whereas that for bilayer @xmath63 is @xmath64 in the large stretching modulus limit ( @xmath65 ) , the tension perturbations @xmath59 become lagrange multipliers , so that they ensure that the right hand side of eq . ( [ delta_phi ] ) vanishes . by using eqs . ( [ tangent_force_balance ] ) and ( [ tangent_force ] ) , the values of @xmath59 are determined from the viscous stress @xmath66 ( see appendix c ) . there are several interactions which contribute to the layer compression modulus @xmath11 . in this subsection , we indicate their physical origins and give the simplest formulae to describe them . for charged bilayers , the electrostatic interaction , together with the counter - ion and co - ion entropy produce the free energy per area  @xcite @xmath67 where @xmath68 is the salt concentration , @xmath69 the potential energy of the counter - ion @xmath70 at the surface , @xmath71 the inverse of the debye screening length , @xmath72 the distance between the two surfaces , and @xmath73 the thermal energy . the van der waals attraction potential per unit area can be calculated by summing the dipoles to obtain @xmath74,\ ] ] where @xmath75 is the hamaker constant , and @xmath76 is the membrane thickness . a more complicated lifschitz theory calculation can provide a more accurate description  @xcite . the sum of these two interactions consists the standard dlvo theory . when the electrostatic repulsion and van der waals attraction stabilize mlv , @xmath11 can be evaluated from @xmath77    for a flexible bilayer where the undulation entropy depends strongly on the membrane separation @xmath20 , helfrich estimated the free energy per area and @xmath11 by a self - consistent argument to obtain  @xcite @xmath78 where @xmath79 is an order unity numerical constant . in the following examples below , we use @xmath80 . when the van der waals interaction becomes important , a more elaborate calculation is required to combine the undulation entropy and the van der waals attraction  @xcite . using the spherical harmonic expansion , the stokes equation can be solved in terms of the radial functions which are simple polynomial of the radius . the detailed calculation is presented in the appendix a. in general , the solution depends on the boundary velocities , which are the velocities of the membranes and the applied external shear . in the usual case where the bilayers have a large stretching modulus @xmath14 , we can simplify the discussion by considering the large stretching modulus limit @xmath65 . in this limit , the surface density approaches unity , and the combinations @xmath81 become the lagrange multipliers to ensure that surface densities are constant . hence the constant @xmath14 can be eliminated . the detailed calculation is presented in appendix b. at the boundaries , this limit also turns the tangential velocity into the function of the normal velocity . this means that the normal stresses @xmath82 and @xmath83 , as well as the tension perturbation @xmath59 are all proportional to the bilayer normal velocities and the external shear . the detailed calculation is presented in appendix c.    the bilayer displacements and velocities are expanded as the sum of the spherical harmonics @xmath48 times the time varying amplitudes , i.e. , @xmath84 and @xmath85 . without cluttering the notation with the angular quantum numbers @xmath86 and @xmath87 , hereafter we use @xmath19 and @xmath88 to express the time varying amplitudes for an arbitrary set of @xmath89 . after performing some calculations to express @xmath90 and @xmath59 by @xmath88 ( see eqs .  ( [ c1 ] ) , ( [ c5 ] ) , ( [ pa ] ) , ( [ pc ] ) and ( [ pb ] ) ) , we find that the normal force balance eq . ( [ normal_force ] ) at @xmath0 and @xmath1 can be written as @xmath91 where @xmath92 . here we prefer to use the variables @xmath93 and @xmath94 which make the matrices @xmath95 and @xmath96 symmetric . with this variable choice , the free energy per solid angle is then the quadratic form of @xmath95 . the components of the matrix @xmath95 are given by @xmath97 where have defined @xmath98 ( @xmath99 ) and @xmath100 whereas the components of the matrix @xmath96 are @xmath101 , \nonumber \\ d_{12}&=\frac{\eta r_2^{4l+1 } ( 4l+2)}{g_0   ( l^2+l ) } \times \nonumber \\ &   \left [ ( 2l^3 + 5l^2+l-2)(\rho^{3l+2 } -\rho^{l-1 } ) \right . \nonumber \\   & \left . + ( 2l^3+l^2 - 3l)(\rho^{l+1}-\rho^{3l})\right ] , \nonumber \\ d_{21}&=d_{12 } , \nonumber \\ d_{22}&=\frac{\eta r_2^{4l+1 } ( 2l+1 ) } { g_0   ( l^2+l ) } \times \nonumber \\ &   \left [ -(l+2)^2(4l^2 - 1)\rho^{2l+4 } \right . \nonumber \\   &   + 2l(l+2)(2l-1)(2l+3)\rho^{2l+2 } \nonumber \\ &    + l^2 ( 4l^2 + 8 l + 3)\rho^{2l } \nonumber \\ &   \left . + ( 8l^2 + 8l-4 ) \rho   \right],\end{aligned}\ ] ] where @xmath102.\end{aligned}\ ] ] notice that the components of @xmath96 are the product of @xmath103 and the functions of @xmath104 and @xmath86 . when the two bilayers are well separated , i.e. , @xmath105 , they are not hydrodynamically coupled . in this limit , @xmath106 and @xmath107 become small , and we recover the isolated vesicle damping given by  @xcite @xmath108    for the more general bilayer interaction which is proportional to @xmath109 ^ 2 $ ] , the similar calculation will give rise to an extra factor @xmath110 for the @xmath11 term in @xmath111 , and an extra factor @xmath112 for the @xmath11 terms in @xmath113 and @xmath114 . the form of @xmath115 is not affected . the relaxation rates , denoted as @xmath116 ( @xmath117 ) , are the eigenvalues of the matrix @xmath118 . the normal force balance condition eq . ( [ normal_force ] ) on the two bilayers gives the eigenvector equation @xmath119 where @xmath120 are the ( non - orthogonal ) right - eigenvectors . to work with the orthogonal eigenvectors , here we define the symmetric decay rate matrix @xmath121 which has the ( normalized ) eigenvectors @xmath122 @xmath123 here the eigenvalues are the same as in eq . ( [ ei1 ] ) , and @xmath122 is proportional to @xmath124 . we can decompose the decay rate matrix as @xmath125 which will be used in the later discussion . hereafter the faster and the slower rates are denoted by @xmath126 and @xmath127 , respectively . the full expression of the decay rates are straightforward , but complicated . hence we shall obtain some simplified expressions to gain some understanding on the nature of the relaxation . we first discuss the fast mode @xmath126 . when @xmath128 , the fast undulation mode has a well - known dispersion relation  @xcite @xmath129 where the bending modulus is doubled as there are two bilayers . for @xmath104 slightly less than unity , a series expansion can be made if desired . for @xmath130 , the decay rate can be approximated by @xmath131 this corresponds to a simple picture that the fast mode consists mainly the inner bilayer relaxation , whereas the outer bilayer does not move much . next we consider the slow mode @xmath127 . when @xmath132 , the slow mode can be approximately obtained by the series expansion of @xmath133 . based on the relation @xmath134 , in which @xmath127 only starts from the second order term , we get @xmath135 for the zeroth and the first order terms . because the product @xmath136 is given by @xmath137 , we use the leading two terms of @xmath138 to obtain @xmath139 the leading term @xmath140 is similar to the slip mode of the planar smectic , which has the decay rate @xmath141 , where @xmath142 is the wave vector projected on the bilayer plane  @xcite . this expression is also analogous to the squeezing mode dispersion obtained for soap film  @xcite . in the other extreme of @xmath143 , the slow mode has the dispersion such that @xmath144 l(l+1 ) } { 3 \eta ( 2l+1)(2l^2 + 2l-1)}. \label{b2a}\ ] ] if we set the shear perturbation as @xmath145 , the numerator contains a factor @xmath146 . hence the dimensionless parameter @xmath147 becomes important when its value is much larger than 72 . we now discuss another approximate expression for the slow mode valid for the whole range of @xmath104 . since the sum of the two decay rates @xmath134 is dominated by the fast mode for any @xmath104 , the slow mode can be approximated by the ratio @xmath148 when the inner bilayer 1 relaxes fast , one can apply the adiabatic approximation to the fast relaxing inner membrane . the approximated expression has even a simpler denominator given by @xmath149     as a function of the dimensionless size ratio @xmath150 between the two layers . here @xmath151 is the vesicle relaxation rate with the rigidity @xmath152 given by eq . ( [ omega0 ] ) . the two solid lines represent the two vesicle relaxation rates @xmath126 and @xmath127 obtained numerically . the dotted line is an approximation for the fast mode given by eq . ( [ omega1 ] ) . the dashed line represents the approximate slow relaxation rate @xmath127 given by eq . ( [ slow_mode ] ) . in this plot , we set @xmath153 j / m@xmath154 , @xmath155 j , @xmath156 pa@xmath157s , @xmath158 m , @xmath145 , so that @xmath159 . ]    in fig . [ fig2 ] , we plot the decay rates as a function of @xmath150 for @xmath145 , keeping a constant @xmath11 . the two solid lines represent the numerical calculated decay rates . the upper one is the fast mode @xmath126 which has the limit @xmath151 at @xmath128 . it coincides with the approximation eq . ( [ omega1 ] ) ( dotted line ) for @xmath160 . the lower solid line represents the slow mode @xmath127 . the approximate slow rate eq . ( [ slow_mode ] ) follows qualitatively the exact value of @xmath127 for the full range of @xmath104 . it also coincides with eq . ( [ b2a ] ) in the limit of @xmath161 . incidentally the inner layer relaxation eq . ( [ omega1 ] ) also fits the slow mode @xmath127 at @xmath132 . [ fig2 ] is useful to illustrate the nature of the relaxation for the fast and slow modes . in a real system , @xmath11 is a function of @xmath162 and @xmath1 ( or @xmath104 and @xmath1 ) , which is hard to be kept as a constant . one will make a plot with a function for @xmath11 which is suitable for the specific mlv system . if one uses the flat layer formulas described in subsection [ sec : interaction ] as approximations , caution should be kept for the accuracy at small @xmath104 . given an accurate function for @xmath11 , together with the proper weighting factor ( to combine @xmath110 and @xmath112 factors in @xmath163 ) , the above decay rate approximations should work better . defined by eq .  ( [ shear ] ) as a function of the dimensionless size ratio @xmath104 between the two layers . ] we now consider two - layer - vesicles under the external oscillatory shear flow . rearranging eq . ( [ effective_eta ] ) together with eq . ( [ pc ] ) , we obtain @xmath164 it is apparent that the dilute hard sphere limit is recovered when @xmath165 . by putting @xmath166 and solving eq . ( [ linear_equation ] ) for @xmath167 , we obtain @xmath168 where one should set @xmath145 for the matrixes @xmath95 and @xmath96 . here @xmath169 appears twice , because the shear directly affects the bilayer 2 through eq . ( [ linear_equation ] ) , and the velocity of the bilayer 2 carries the stress contribution of the vesicle according to eqs . ( [ effective_eta ] ) and ( [ pc ] ) .    considering the high - frequency limit , we now define the dimensionless shear coupling strength @xmath170 and the shear deformation unit vector @xmath171 by @xmath172 notice that @xmath170 and @xmath171 only depend on the ratio @xmath104 as long as we fix to @xmath145 . as shown in fig . [ fig3 ] , @xmath170 is weakly depended on @xmath104 because its value varies only between 2.1 and 2.19 . hence the right hand side of eq . ( [ high_shear_eta2 ] ) varies between 0.31 and 0.4 , whereas the hard sphere case gives 2.5 . this means that , at high frequencies , the coupling of a vesicle to the flow is weaker compared to that of a hard sphere .    ) as a function of the dimensionless size ratio @xmath104 between the two layers . the solid and the dashed lines corresponds to @xmath173 and @xmath174 , respectively . we chose the value @xmath175 . the crossing of the two modes occurs at @xmath176 . ] the high frequency viscosity asymptotically approaches @xmath177 = \eta^0-\eta s \phi_{\rm v},\ ] ] where @xmath178 $ ] . using the eigenmode decomposition @xmath179 , we can separate the viscosity contributions from each mode as @xmath180 so that we have @xmath181 then the full expression for the complex viscosity is given by @xmath182 the relative strength of the two modes are shown in fig . [ fig4 ] . alternatively ( but equivalently ) , the complex modulus can be expressed as @xmath183 where the two modes have the positive amplitudes @xmath184 . ( divided by @xmath185 ) of various 2d vectors as a function of the dimensionless size ratio @xmath104 between the two layers . the solid line represents @xmath171 , while the dotted and the dashed lines represent @xmath186 and @xmath187 , respectively . we chose the value @xmath188 . ]    even for different bilayer interaction strength @xmath11 , we always find that the slow mode has the larger viscosity amplitude at @xmath189 , and smaller one at @xmath132 . the latter is reasonable because the shear perturb the two bilayers similarly for @xmath132 , where @xmath190 is along the @xmath191-direction . in terms of the polar angle @xmath192 between @xmath190 and the first axes on the @xmath193-plane , the @xmath191-direction corresponds to @xmath194 . notice that this direction also corresponds to the undulation eigenmode direction . since the sum of the two viscosity amplitudes is roughly a constant as shown in eq . ( [ sum ] ) , the slow squeezing mode takes a small viscosity amplitude @xmath174 for @xmath132 . in the other limit of @xmath189 , the shear mainly perturbs the outer layer . in this limit , we have @xmath195 or @xmath196 , and the slow mode is due to the outer layer relaxation . as a result , @xmath174 becomes the dominant contribution for @xmath189 . in fig . [ fig5 ] , we plot the angle @xmath192 as a function of @xmath104 . the shear vector @xmath190 always coincides with the fast mode at @xmath128 and the slow mode at @xmath161 . therefore the mode switching always takes place between @xmath197 . it is also worth mentioning that figs . ( [ fig4 ] ) and ( [ fig5 ] ) , like fig . ( [ fig2 ] ) , are plotted by assuming a constant @xmath11 . these plots are used to analyze the behavior of the mode amplitudes . the more physical plots will be the ones using a suitable function for @xmath11 .     as a function of the dimensionless bilayer interaction @xmath198 . for @xmath199 , the fast mode is the dominant contribution to the viscosity , while the slower relaxation is dominant for @xmath200 . ]    here we define a crossover size ratio @xmath201 at which the two modes have the same viscosity amplitude , i.e. , @xmath202 . in fig . [ fig6 ] , we plot the calculated @xmath201 as a function of @xmath147 . when @xmath203 , the crossover happens at @xmath204 . this means that for non - interacting case @xmath205 , the crossover happens when @xmath206 is roughly the same as @xmath0 . when @xmath207 , on the other hand , @xmath201 approaches unity . as mentioned in eq . ( [ b2a ] ) , the parameter @xmath208 will have a noticeable effect only when it is greater than 72 . at the lower right corner of fig . [ fig6 ] , where the layer interaction is strong and the layer separation is not too small , the slow mode has the dominant viscosity contribution . at the upper left corner where the layer interaction is relatively weak , the fast mode is more excited as compared to the slow mode . in summary , we have calculated the slow relaxation rates and the viscoelasticity of a dilute two - layer - vesicle solution . we have found the following points : ( i ) at small gap @xmath132 , the slowest mode is the squeezing mode . the undulation mode appears to be faster . ( ii ) when the inner bilayer radius is small @xmath189 , the slow mode becomes the relaxation of the outer bilayer , and the faster mode is the relaxation of the inner bilayer . ( iii ) when one decreases @xmath104 , the relaxation spectrum changes from ( i ) to ( ii ) . ( iv ) for the complex viscosity , the low frequency viscosity approaches the hard sphere limit . the high frequency viscosity increment is between 1216 % of the hard sphere viscosity increment . the difference between the two limits comes from the contributions of the two modes . ( v ) at small gap @xmath199 , the slow ( squeezing ) mode has a small viscosity amplitude , while the fast ( undulation ) mode has a large viscosity amplitude . for large gap @xmath200 , on the other hand , the slow mode has the dominant viscosity amplitude . ( vi ) the crossover size ratio @xmath201 depends on the interaction between the two bilayers . as @xmath198 is increased , @xmath201 increases from 0.52 toward unity .     as a function of the size ratio @xmath104 . the crossover ratio @xmath201 in fig . [ fig6 ] is plotted by the solid black line . the blue lines are from the dlvo theory eq . ( [ dlvo ] ) with @xmath209 m , @xmath210 m , @xmath211 j , @xmath212 m , @xmath213 j. the surface potentials are @xmath214 mv and @xmath215 mv for the blue solid and dashed lines , respectively . the red lines are the helfrich repulsion eq . ( [ helfrich ] ) with @xmath80 , @xmath158 m. the red solid and dashed lines are for @xmath216 j and @xmath217 j , respectively . ] we have determined the crossover ration @xmath218 in fig . [ fig6 ] . the actual bilayer interaction strength @xmath11 depends on the separation @xmath219 . for the qualitative discussion at non - small @xmath104 , we use the flat layer results in sec .  [ sec : interaction ] . in fig . [ fig7 ] we present both information in the @xmath220-plot to compare a series of systems with the same outer bilayer size @xmath1 but with different size ratios @xmath104 . as one varies @xmath104 , the dimensionless bilayer interaction strength @xmath147 changes according to the dlvo theory eq . ( [ dlvo ] ) ( blue lines ) or the helfrich repulsion eq . ( [ helfrich ] ) ( red lines ) . if such lines happen to cross the line @xmath218 ( black line ) , one expects that the two viscosity amplitudes change their relative magnitudes . we have indicated such scenarios by the dashed lines .    as shown in the blue dashed line in fig . [ fig7 ] , an electrostatically stabilized system may get large @xmath11 by having a high surface charge density . at small gap , the van der waals attraction may also lower @xmath11 , causing an interesting multiple crossing . this means that the viscosity amplitudes switch their magnitudes more than once . for sterically stabilized system ( red lines ) , we find that the dimensionless interaction parameter depends strongly on the bending rigidity as @xmath221 ( see eq . ( [ helfrich ] ) ) . for a soft surfactant bilayer of @xmath222 , we find that the slow mode always dominates the viscosity . whereas for a lipid bilayer whose bending rigidity @xmath223 is one order of magnitude larger , the red dashed line indicates that the squeezing mode is not much excited by shear when @xmath224 .    in this paper , we have only presented the calculation for two - layer - vesicles . for mlv with more than two layers , the calculation can be performed in the same procedure , but with a greater algebraic complexity . for mlv with @xmath225 layers , we expect that there are @xmath225 relaxation modes . when the gap is small ( @xmath226 ) , we expect that the majority of the @xmath225 relaxations to bear some resemblance to the squeezing mode or the spherical version of the `` slip mode ''  @xcite . for larger gap , on the other hand , the relaxations of the @xmath225 layers may decouple from each other . as for the viscoelasticity , we speculate that the viscosity amplitudes of the squeezing modes are small for weakly interacting mlv . neutral or weakly charged lipid mlv should be an interesting system to investigate in this direction .    in polymer rheology calculation , one may consider a step strain for @xmath227 , where the time relaxation of the stress gives the relaxation modulus @xmath228 . this quantity can be further converted to the complex modulus @xmath229 by the fourier transform . right after the step strain , one often assumes that the polymer deforms in an affine way . do we implicitly use the affine deformation approximation for the mlv rheology ? a step strain for vesicle solution will induce a short but fast bilayer movement , where the @xmath96 terms dominate the left hand side of eq . ( [ linear_equation ] ) . therefore the initial layer displacements will be proportional to @xmath230 which is compatible with the incompressible constraints ( see eq . ( [ 2dincompressible ] ) ) at both layers . in the small gap limit , @xmath231 behaves like @xmath232 . compared with the affine deformation @xmath233 , or in terms of our chosen variables @xmath234 , it is clear that our theory does not use the affine deformation approximation . we thank m.  cates and s.  fujii for many useful discussions . we also acknowledge the support from the national science council of taiwan and the center of theoretical physics of national taiwan university . sk also acknowledges the supported by grant - in - aid for scientific research ( grant no . 24540439 ) from the mext of japan , and the jsps core - to - core program _ `` international research network for non - equilibrium dynamics of soft matter''_. for the incompressible solenoidal flow , the velocity can be expressed as @xmath235 where the scalar functions @xmath236 and @xmath237 are the defining function for the poloidal and toroidal flow fields , respectively . note that our definition of the defining function differs by a factor @xmath38 compared with the ones in the book by chandrasekhar  @xcite . for an incompressible fluid , the divergence of the pressure gradient vanishes , i.e. , @xmath238 therefore the pressure gradient is also a solenoidal field , and can be described by the above decomposition . since the toroidal part can be written as @xmath239 , the condition @xmath240 requires that @xmath241 or @xmath242 , where @xmath100 hence the pressure gradient is only poloidal , and written as @xmath243 where @xmath244 is the defining function of @xmath245 . comparing the @xmath246 and @xmath247 directions , we can set @xmath248 for these two directions . the laplace equation for the pressure eq . ( [ p ] ) then implies that @xmath249 therefore the radial component of the gradient pressure can be either expressed as @xmath250 or @xmath251 . the latter expression is consistent with the identification eq . ( [ p2 ] ) .    for the velocity field , we will drop the toroidal part by setting @xmath252 . this is justified because the tangential force field is a surface gradient and drives only the poloidal flow . to obtain the defining function @xmath236 , we take the curl of eq . ( [ stokes ] ) to get @xmath253 . since @xmath254 holds for incompressible flow , we have @xmath255=0,\ ] ] where in the square bracket , an equivalent form of @xmath256 is used . as detailed in ref . @xcite , each curl will switch poloidal and toroidal parts . the double curl will preserve the type and modify the defining function by @xmath257 . the curl of eq . ( [ stokes ] ) becomes @xmath258 or @xmath259    to find the coupling between the pressure @xmath35 and the velocity poloidal function @xmath236 in eq . ( [ stokes ] ) , we rewrite @xmath260 $ ] . then eq . ( [ stokes ] ) is satisfied when @xmath261 here we prefer to use @xmath35 instead of @xmath244 . operating the above equation by @xmath262 and substituting the combination @xmath250 by @xmath263 , we obtain @xmath264    using the standard method of separation of variables , the general solution of eqs . ( [ p ] ) , ( [ psi ] ) and ( [ ppsi ] ) are @xmath265 @xmath266 where @xmath48 are the spherical harmonics . the surface incompressibility condition also constrains the velocity field as @xmath267 at both @xmath268 and @xmath1 . because the fluid is incompressible , i.e. , @xmath269 , eq . ( [ 2dincompressible ] ) can also be expressed as the equivalent form @xmath270 at the bilayers . we prefer this condition because it simplifies the calculation . we now consider the relaxation rates of a small perturbation described by the spherical harmonics @xmath48 . to simplify the notation , we drop the subscript `` @xmath271 '' of the coefficients @xmath272 , @xmath273 , @xmath274 , and @xmath275 . for the region within the bilayer 1 ( @xmath276 ) , we set @xmath277 and @xmath278 to drop the functions which are singular at @xmath279 . we denote @xmath272 and @xmath274 as @xmath280 and @xmath281 , respectively . the 2d incompressibility condition eq .  ( [ 2dincompressible2 ] ) relates the two remaining coefficients as @xmath282 in terms of the radial velocity amplitude @xmath283 , we can further express the pressure coefficient as @xmath284    for the region exterior to the bilayer 2 ( @xmath285 ) , the coefficients @xmath272 are zero , except @xmath286 , which needs to be chosen to give the far flow eq . ( [ v_shear ] ) . comparing the radial component @xmath28 from eq . ( [ v_shear ] ) and @xmath287 , we set @xmath288 . we also set @xmath289 so that @xmath35 does not diverge at @xmath290 , and denote @xmath273 and @xmath275 as @xmath291 and @xmath292 , respectively . then the 2d incompressibility condition eq . ( [ 2dincompressible2 ] ) relates the two remaining coefficients as @xmath293 using the radial velocity amplitude @xmath167 , we can express the pressure coefficient as @xmath294    for the region between the two bilayers ( @xmath295 ) , the expressions are more complex . the incompressibility eq . ( [ 2dincompressible2 ] ) , evaluated at @xmath0 and @xmath1 , provides two conditions between the four coefficients @xmath296 where @xmath297 we prefer to use the variables @xmath283 and @xmath167 instead of @xmath298 and @xmath299 . then the pressure coefficients can be expressed as @xmath300 where @xmath301 \right\ } , \nonumber \\ g_{12}&=   \frac{\eta}{l g_0 } \left\ { ( -8l^2 - 4l+12)r_1 r_2^{3l+2 } \right . \nonumber \\   & - \left . ( 6l+4)r_1^{2l } r_2^{l+1}[l(2l+1)r_2 ^ 2-(l+2)(2l-1)r_1 ^ 2 ] \right\ } ,   \nonumber \\ g_{21}&=   \frac{(4l-2)\eta}{(l+1)g_0 } \left\ { -(2l^2 + 3l+1)r_1^{3l+4 } r_2^{2l } \right . \nonumber \\   & +   \left . ( 2l^2+l-3)r_1^{3l+2 } r_2^{2l+2}+ ( 2l+4)r_1^{l+1 } r_2^{4l+3 } \right\ } , \nonumber \\ g_{22}&= \frac{(4l-2)\eta}{(l+1)g_0 } \left\ { -(2l^2 + 3l+1)r_2^{3l+4 } r_1^{2l } \right . \nonumber \\   & +   \left . ( 2l^2+l-3)r_2^{3l+2 } r_1^{2l+2}+ ( 2l+4)r_2^{l+1 } r_1^{4l+3 } \right\},\end{aligned}\ ] ] with @xmath102.\end{aligned}\ ] ] when the two bilayer are well separated ( @xmath105 ) , both @xmath302 and @xmath303 become small . in this limit , @xmath304 becomes the coefficient of eq . ( [ pa ] ) with @xmath0 replaced by @xmath1 , whereas @xmath305 becomes the coefficient of eq . ( [ pc ] ) ( without the @xmath47 term ) with @xmath1 replaced by @xmath0 . the radial component of the normal stress appears in eq . ( [ normal_force ] ) . from its component form @xmath306 and the 2d incompressibility condition eq . ( [ 2dincompressible2 ] ) , it is just the negative pressure . therefore the stress differences at the two bilayers are @xmath307 y_{lm } , \nonumber \\ \sigma_{rr}(r_2^+)-\sigma_{rr}(r_2 ^ - ) & = \sum_{lm}\left [ p_{\rm b}^{\rm i } r_2^l   + ( p_{\rm b}^{\rm ii}-p_{\rm c}^{\rm ii})r_2^{-l-1 } \right ] y_{lm}. \label{c1}\end{aligned}\ ] ]    the large bilayer stretching modulus @xmath14 suppresses the surface density perturbation . for the slow relaxation , we will take the limit @xmath65 to eliminate the parameter @xmath14 . within this limit , both @xmath308 and @xmath309 approach unity , so that the tension perturbation @xmath310 and @xmath311 become the lagrange multipliers . then the values of the lagrange multipliers are determined by eqs . ( [ tangent_force_balance ] ) and ( [ tangent_force ] ) instead of their original definitions , as discussed below .    in eq . ( [ tangent_force_balance ] ) , we are interested in the difference between the tangential stress across the bilayer . we express the tangential stress as a function of @xmath236 : @xmath312.\ ] ] we now replace @xmath313 using eq . ( [ ppsi ] ) , and limit our discussion to @xmath314 mode . the velocity field @xmath315 has the components similar to eq . ( [ component ] ) . we can use the @xmath38 component @xmath316 to eliminate @xmath236 as @xmath317 so that @xmath318,\ ] ] where @xmath319 for spherical harmonics with nonzero @xmath86 . because @xmath28 is continuous across the bilayer and @xmath320 vanishes on the bilayer , only the first term can be different across the bilayer . this is the only important term for the tension perturbation in eqs . ( [ tangent_force_balance ] ) and ( [ tangent_force ] ) : @xmath321 we then obtain @xmath322
Abstract:

Answer: the dynamic shape relaxation of the two - layer - vesicle is calculated . in additional to the undulation relaxation where the two bilayers move in the same direction , the squeezing mode appears when the gap between the two bilayers is small . at large gap , the inner vesicle relaxes much faster , whereas the slow mode is mainly due to the outer layer relaxation . we have calculated the viscoelasticity of the dilute two - layer - vesicle suspension . it is found that for small gap , the applied shear drives the undulation mode strongly while the slow squeezing mode is not much excited . in this limit the complex viscosity is dominated by the fast mode contribution . on the other hand , the slow mode is strongly driven by shear for larger gap . we have determined the crossover gap which depends on the interaction between the two bilayers . for a series of samples where the gap is changed systematically , it is possible to observe the two amplitude switchings .