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: the third author of this paper still vividly recollects his visit to victor klee at the university of washington ( seattle ) in august 2000 .    _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ we had a memorable drive in vic s car from the hotel to the math department . vic skipped all the small talk and immediately asked : `` do you think that the simplex method is polynomial ? '' i still remember how awkward i felt to be asked this question by the very person who had provided the first and seminal insights into it . my ( then as now quite irrelevant ) opinion shall remain between vic and me forever . _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _    the seminal work referred to above is of course the 1972 paper of klee and minty entitled `` how good is the simplex algorithm ? '' it deals with the number of iterations that may be required in the worst case to solve a linear program ( lp ) by dantzig s _ simplex method _ @xcite , which was at that time the only available practical method . klee and minty exhibited a family of lps ( now known as the _ klee - minty cubes _ ) for which the number of iterations is _ exponential _ in the number of variables and constraints . dantzig s specific _ pivot rule _ for selecting the next step is just one of many conceivable rules . the question left open by the work of klee and minty is whether there is some _ other _ rule that provably leads to a small number of iterations . this question is as open today as it was in 1972 . nowadays , the simplex method is no longer the only available method to solve lps . in particular , there _ are _ proven polynomial - time algorithms for solving lps  khachiyan s _ ellipsoid method _ @xcite and karmarkar s _ interior point method _ @xcite  that are based on techniques developed originally for nonlinear optimization . it is still unknown , though , whether there is a _ strongly _ polynomial - time algorithm for solving lps , that is , an algorithm for which the number of arithmetic operations does not depend on the bit lengths of the input numbers . [ [ the - p - matrix - linear - complementarity - problem . ] ] the p - matrix linear complementarity problem . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    in this paper we are concerned with pivoting methods for the _ p - matrix linear complementarity problem _ ( p - lcp ) , a prominent problem for which neither polynomial - time algorithms nor hardness results are available .    in general , an lcp is given by a matrix @xmath0 and a vector @xmath1 , and the problem is to find vectors @xmath2 such that @xmath3 it is @xmath4-complete to decide whether such vectors exist @xcite . however , in a p - lcp , @xmath5  is a p - matrix ( meaning that all principal minors are positive ) , and then there are unique solution vectors @xmath6 for every right - hand side  @xmath7  @xcite . the problem of efficiently finding them is unsolved , though . megiddo has shown that the p - lcp is unlikely to be np - hard . for this , he considers the following ( more general ) variant : given @xmath5 and @xmath7 , either find @xmath6 that satisfy ( [ eq : lcpdef ] ) , or exhibit a nonpositive principal minor of @xmath5 . np - hardness of the latter variant would imply that @xmath8  @xcite . there is a different notion of hardness that might apply to the p - lcp as a member of the complexity class  @xmath9  @xcite . this class has complete problems , and no polynomial - time algorithm is known for any of them . it has recently been shown that the famous problem of finding a nash equilibrium in a bimatrix game is @xmath9-complete  @xcite . some researchers consider this to be an explanation why ( despite many efforts ) no polynomial - time algorithm has been found so far . incidentally , this is the _ second _ problem for which megiddo showed in his technical report  @xcite that it is unlikely to be np - hard . however , it is currently not known whether the p - lcp is also ppad - complete .    there are various algorithms for solving p - lcps , among them the classic method by lemke @xcite as well as interior point approaches @xcite . still , there is no known polynomial - time algorithm for p - lcps . for example , the complexity of interior point algorithms depends linearly on a matrix parameter @xmath10 @xcite that is not bounded for p - matrices .    in this work , we focus on combinatorial methods for p - lcps and in particular on _ simple principal pivoting methods _ that share their essential idea with the simplex method . [ [ simple - principal - pivoting - methods . ] ] simple principal pivoting methods . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    let @xmath11 , and let @xmath12 be the @xmath13 matrix whose @xmath14th column is the @xmath14th column of  @xmath15 if @xmath16 , and the @xmath14th column of the @xmath13 identity matrix  @xmath17 otherwise . if @xmath5 is a p - matrix , then @xmath12 is invertible for every set @xmath18 . we call @xmath18 a _ basis_. if @xmath19 , we have discovered the solution : let @xmath20    if on the other hand @xmath21 , then @xmath22 and @xmath23 defined above satisfy @xmath24 and @xmath25 , but @xmath26 fails . in _ principal pivoting _ , one tries to improve the situation by replacing the basis @xmath18 with the symmetric difference @xmath27 , where @xmath28 is some nonempty subset of the `` bad indices '' @xmath29 . the greedy choice is to let @xmath28 be the set of all bad indices in every step . for some p - lcps , however , this algorithm cycles and never terminates  @xcite . , one such example is the p - lcp ( [ eq : morrismq ] ) whose digraph model ( see next paragraph ) appears in figure [ fig : uso ] . cycling occurs for any of the three start bases just below the top vertex in figure  [ fig : uso ] . ] in _ simple principal pivoting _ , @xmath28 is of size one , and a _ pivot rule _ is employed to select the bad index . simple principal pivoting methods are sometimes called _ bard - type _ methods and were first studied by zoutendijk  @xcite and bard  @xcite . for a detailed exposition on simple principal pivoting methods see , for instance , ( * ? ? ? * section  4.2 ) or  ( * ? ? ? * chapter  4 ) . [ [ the - digraph - model - and - unique - sink - orientations . ] ] the digraph model and unique sink orientations . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    there is a natural digraph model behind principal pivoting , first studied by stickney and watson @xcite . the graph s vertices are all bases  @xmath18 , with a directed edge going from @xmath18 to @xmath30 if @xmath31 and @xmath32 . to be nondegenerate , meaning that @xmath33 for all @xmath18 and @xmath14 this is no serious restriction since a nondegenerate problem can always be obtained from a degenerate one by a symbolic perturbation of @xmath7 . ] the underlying undirected graph is the graph of the @xmath34-dimensional cube . a principal pivot step considers the cube that is spanned by the outgoing edges at the current vertex , and it jumps to the antipodal vertex of some subcube of it . a simple principal pivot step can be interpreted as the traversal of an outgoing edge , and this is the analogy with the simplex method for linear programming . the above digraph has the following specific property : every nonempty subcube ( including the whole cube ) has a unique sink , and the global sink corresponds to the solution of the p - lcp @xcite . in the terminology of szab and welzl @xcite , we are dealing with a _ unique - sink orientation _ ( uso ) . the goal of this paper is to deepen the understanding of the uso description of simple principal pivoting methods for the p - lcp . on a general level , we want to understand whether and how algebraic properties of a p - lcp ( which are well studied ) translate to combinatorial properties of the corresponding uso ( which are much less studied ) , and what the algorithmic implications of such translations are . the principal motivation behind this research is the continuing quest for a polynomial - time p - lcp algorithm . in this paper , the concept of a uso serves as the main tool to obtain two more positive results for simple principal pivoting . we believe that this combinatorial approach has some untapped potential for the theory of ( simple ) principal pivoting methods for lcps .    [ [ the - randomized - method - of - murty - and - the - morris - orientations . ] ] the ( randomized ) method of murty and the morris orientations . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    the simple principal pivoting method of murty ( also known as the _ least - index _ rule ) @xcite works for all p - lcp instances , but it may take exponentially many iterations in the worst case @xcite . as a possible remedy , researchers have considered _ randomized _ methods . the two most prominent ones are randomizedmurty ( which is just the least - index rule , after randomly reshuffling the indices at the beginning ) , and randomedge ( which performs a purely random walk in the uso ) . however , morris found a family of p - lcp instances ( we call their underlying digraphs the _ morris orientations _ ) on which randomedge spends a long time running in cycles and thus performs much worse than even the exhaustive search algorithm  @xcite . for randomizedmurty there are as yet no such negative results on p - lcp instances . in particular , on murty s worst - case example , this algorithm takes an expected _ linear _ number of steps @xcite ; the expected number of steps becomes quadratic if we allow arbitrary start vertices  @xcite .    on general p - lcp instances , the expected number of iterations taken by randomizedmurty is _ subexponential _ if the underlying digraph is acyclic  @xcite , but it is unknown whether this also holds under the presence of directed cycles ( as they appear in the morris orientations , for example ) . we prove that not only the randomized variant , but actually _ any _ variant of the least - index rule obtained by some initial reshuffling of indices leads to a _ number of iterations on the morris orientations . in particular , this `` kills '' another family of potentially bad instances for the randomizedmurty rule .    [ [ k - matrix - linear - complementarity - problems . ] ] k - matrix linear complementarity problems . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    a _ k - matrix _ is a p - matrix so that all off - diagonal elements are non - positive . lcps with k - matrices ( k - lcps ) appear for example in free boundary problems ( * ? ? ? * section  5.1 ) and in american put option pricing @xcite . it is known that every k - lcp instance can be solved in polynomial time ( even by a simple principal pivoting method  @xcite , see also  ( * ? ? ? * section  4.7 ) ) , but we prove something stronger : _ every _ simple principal pivoting method takes only a linear number of iterations . we obtain this result by showing that in k - lcp - induced usos , _ all _ directed paths are short . our approach is to extract combinatorial structure from the algebraic properties of k - matrices . subsequently , we use the distilled combinatorial information to get our structural results .    [ [ lcps - with - sufficient - matrices - and - the - criss - cross - method . ] ] lcps with sufficient matrices and the criss - cross method . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    let us step back and take a broader view that brings together lp , the p - lcp , pivoting methods , and computational complexity : linear complementarity problems with _ sufficient _ matrices generalize both lp and the p - lcp while still being amenable to pivoting approaches . the _ criss - cross method _ @xcite may be considered an extension of murty s algorithm . in fact , the results of  @xcite make lcps with sufficient matrices the largest known class of lcps for which megiddo s techniques of  @xcite still apply , so that np - hardness is unlikely . determining the complexity of sufficient matrix lcp remains a tough challenge for the future . now we formally introduce the digraph model behind p - lcps and show some of its basic properties . the model was first described by stickney and watson in 1978 @xcite . we use the following notation . let @xmath35:=\{1,2,\ldots , n\}$ ] . for a bit vector @xmath36 and @xmath37 $ ] , let @xmath38 be the element of  @xmath39 defined by @xmath40 this means that @xmath38 is obtained from @xmath41 by flipping all entries with coordinates in  @xmath42 . instead of @xmath43 we simply write @xmath44 .    under this notation , the ( undirected ) _ @xmath34-cube _ is the graph @xmath45 with @xmath46\}.\ ] ]    a _ subcube _ of  @xmath47 is a subgraph @xmath48 of  @xmath47 where @xmath49 for some vertex @xmath41 and some set @xmath50 $ ] , and @xmath51 . [ def : uso ] let @xmath52 be an orientation of the @xmath34-cube ( a digraph with underlying undirected graph @xmath47 ) . we call @xmath52 a _ unique - sink orientation _ ( uso ) if every nonempty subcube has a unique sink in @xmath52 .    figure  [ fig : uso ] depicts a uso of the 3-cube . the conditions in definition [ def : uso ] require the orientation to have a unique global sink , but in addition , all proper nonempty subcubes must have unique sinks as well . in figure  [ fig : uso ] there are six 2-dimensional subcubes , twelve 1-dimensional subcubes ( edges ) and eight 0-dimensional subcubes ( vertices ) . for edges and vertices , the unique - sink condition is trivial . the figure also shows that usos may contain directed cycles . ( @xmath53  for an outgoing edge ; @xmath54  for an incoming edge ) . the unique global sink and a directed cycle are highlighted . [ fig : uso ] ]    if @xmath55 is the vertex set of a subcube ( as defined above ) , then the directed subgraph of  @xmath56 induced by  @xmath55 is denoted by  @xmath57 . it is convenient to view the orientation of the @xmath34-cube  @xmath47 as a mapping @xmath58 , where @xmath59 if the edge between @xmath41 and  @xmath44 is oriented towards @xmath44 , and @xmath60 if it is oriented towards  @xmath41 . to encode a vertex  @xmath41 along with the orientations of its incident edges , we can then use a _ configuration _ of @xmath34  bits and @xmath34  signs , where the @xmath14th sign is @xmath61 ( see figure  [ fig : uso ] ) . if an orientation @xmath52 of @xmath47 contains the directed edge @xmath62 , we write @xmath63 , or simply @xmath64 if @xmath52  is clear from the context . let @xmath52 be an orientation of the @xmath34-cube and let @xmath65 $ ] . then @xmath66  is the orientation of the @xmath34-cube obtained by reversing all edges in coordinates contained in  @xmath67 ; formally @xmath68    if @xmath52  is a uso , then @xmath66 is also a uso for an arbitrary subset  @xmath65 $ ] .    this proposition can be proved by induction on the size of  @xmath67 ; it essentially suffices to show that under @xmath69 , the whole cube still has a unique sink . now we can derive the _ unique - completion property_. its meaning is that we can prescribe the bits for some coordinates and the signs for the remaining coordinates , and in a uso there will always be a unique vertex that satisfies the prescription . in particular , it follows that any uso ( and any of its subcubes ) also has a _ unique source_.    in fact , the unique - completion property characterizes unique - sink orientations . [ lem : ucp ] an orientation  @xmath52 of the @xmath34-cube  @xmath47 is a unique - sink orientation if and only if for every partition @xmath35=a\cup b$ ] and every pair of maps @xmath70 , @xmath71 there exists a unique vertex  @xmath41 such that    [ eq : ucp ] @xmath72    first suppose that @xmath52  is a uso . let @xmath73 and consider the subcube  @xmath74 . let @xmath75 . as we have just noted , @xmath66  is a unique - sink orientation , and thus the subcube  @xmath74 has a unique sink  @xmath41 with respect to  @xmath66 . this vertex  @xmath41 is the only vertex that satisfies  ( [ eq : ucp]a ) and ( [ eq : ucp]b ) . conversely , let @xmath52 be an orientation of the @xmath34-cube which satisfies the unique - completion property and let @xmath76 . set @xmath77 and @xmath78\setminus c$ ] , and let @xmath79 for @xmath80 and @xmath81 for @xmath16 . then the unique vertex  @xmath41 satisfying  ( [ eq : ucp]a ) and ( [ eq : ucp]b ) is the unique sink of the subcube  @xmath74 .    [ [ p - lcp - induced - usos . ] ] p - lcp - induced usos . + + + + + + + + + + + + + + + + + + +    we now formally define the uso that is induced by a nondegenerate p - lcp instance @xmath82 @xcite . for @xmath36 , let @xmath83 : v_j=1\}$ ] be the canonical `` set representation '' of  @xmath41 . then the unique - sink orientation @xmath52 induced by @xmath82 is @xmath84 recall that @xmath12 is the ( invertible ) @xmath13 matrix whose @xmath14th column is the @xmath14th column of  @xmath15 if @xmath16 , and the @xmath14th column of the @xmath13 identity matrix  @xmath17 otherwise .    in this way , we have reduced the p - lcp to the problem of finding the sink in an implicitly given uso . access to the uso is gained through a _ vertex evaluation oracle _ that for a given vertex @xmath36 returns the orientations @xmath85 of all incident edges . in this setting , the aim is to find the sink of the uso so that the number of queries to the oracle would be bounded by a polynomial in the dimension  @xmath34 . in the case of a p - lcp - induced uso , a ( strongly ) polynomial - time implementation of the vertex evaluation oracle is immediate from  . from the sink of  @xmath52 , we can reconstruct solution vectors @xmath22 and  @xmath23 as in  . the fact that the @xmath34-cube orientation defined in this way is indeed a unique - sink orientation was proved by stickney and watson  @xcite .    not every uso is p - lcp - induced ; for @xmath86 , the p - lcp - induced usos are characterized in  @xcite , but for @xmath87 , we have only a necessary condition . this condition was originally proved by holt and klee for polytope graphs oriented by linear functions  @xcite . every p - lcp - induced uso of the @xmath34-cube satisfies the _ holt - klee condition _ , meaning that there are @xmath34 directed paths from the source to the sink with pairwise disjoint sets of interior vertices.[thm : hk ] in the uso setting resulting from p - lcps , simple principal pivoting can be interpreted as follows : start from any vertex , and then proceed along outgoing edges until the global sink ( and thus the solution to the p - lcp ) is reached . a _ pivot rule _ r determines which outgoing edge to choose if there is more than one option . here is the generic algorithm , parameterized with the rule r. it outputs the unique sink of the given uso  @xmath52 . @xmath88    choose @xmath36 choose @xmath89 according to the pivot rule r @xmath90 @xmath41    [ alg : pivot ]    note that when @xmath52 contains directed cycles as in figure [ fig : uso ] , this algorithm may enter an infinite loop for certain rules r , but we consider only rules for which this does not happen . in the uso setting we are restricted to _ combinatorial _ rules . these are rules that may access only the orientations of edges @xmath91 ( as given by the vertex evaluation oracle ) but not any numerical values such as @xmath92 that define these orientations . let us now introduce the combinatorial pivot rules that are of interest to us . we write them as functions of the set @xmath93 of candidate indices . the first one is murty s least - index rule @xmath94 that simply chooses the smallest candidate . in the cube , this may be interpreted as a `` greedy '' approach that always traverses the first outgoing edge that it finds ( in the order of the cube dimensions ) .    to    it is easy to prove by induction on @xmath34 that this rule leads to a finite number of iterations , even if the uso contains directed cycles . the rule has a straightforward generalization , using a fixed permutation @xmath95 : choose the index @xmath96 with the smallest @xmath14 such that @xmath97 . this just reshuffles the cube dimensions by applying @xmath98 , and with @xmath99 , we recover murty s least - index rule . the natural randomized variant of @xmath100 is randomizedmurty : at the beginning , choose the permutation  @xmath98 uniformly at random from the set  @xmath101 of all permutations , and then use @xmath100 throughout . finally , randomedge stands for the plain random walk : in each iteration , simply choose a random candidate .    to    unlike the previous rules , randomedge may lead to cycling in @xmath102 , but the algorithm still terminates with probability @xmath103 , and with an expected number of at most @xmath104 iterations . this is a simple consequence of the fact that there is always a short directed path to the sink . [ lem : hamming ] let @xmath105 be the global sink of an @xmath34-cube uso  @xmath52 , and let @xmath36 be any vertex . if @xmath106 is the hamming distance between @xmath41 and  @xmath107 , then @xmath52  contains a directed path of length  @xmath106 from @xmath41 to  @xmath107 . morris proved that under @xmath108 , algorithm  [ alg : pivot ] may be forced to perform an expected _ exponential _ number of iterations @xcite . more precisely , at least @xmath109 iterations are required on average to find the sink of the @xmath34-cube uso generated ( as described in section [ sec : uso ] ) by the @xmath110 , with @xmath111 here , @xmath34 must be an odd integer in order for @xmath5 to be a p - matrix .    to prove this result morris first extracted the relevant structure of the underlying uso , and then showed that on this uso , randomedge runs in cycles for a long time before it finally reaches the global sink . it was left as an open problem to determine the expected performance of randomizedmurty , the other natural randomized rule , on this particular uso . we solve this open problem by showing that for _ any _ permutation @xmath98 , @xmath100 incurs only @xmath112 iterations on the @xmath34-dimensional morris orientation . in particular , this bound then also holds for randomizedmurty . [ thm : general ] for any @xmath113 , algorithm [ alg : pivot ] with @xmath114 finds the sink of the @xmath34-dimensional morris orientation after at most @xmath115 iterations , from any start vertex . computational experiments suggest that the worst running time is attained for the identity permutation if the algorithm starts from @xmath116 ( the global source ) . therefore we analyze this case thoroughly and compute the precise number of iterations . [ thm : id ] algorithm [ alg : pivot ] with @xmath117 finds the sink of the @xmath34-dimensional morris orientation after @xmath118 iterations , starting from @xmath116 . experimental data suggest that this is an upper bound on the number of iterations for any permutation @xmath98 and start vertex @xmath116 . allowing start vertices different from @xmath116 , we get slightly higher iteration numbers , but @xmath99 continues to be the worst permutation . in view of this , we conjecture that @xmath119 iterations actually suffice for all @xmath98 and all start vertices . this would mean that the bound in theorem [ thm : general ] is still off by a factor of  @xmath120 . [ [ the - combinatorial - structure . ] ] the combinatorial structure . + + + + + + + + + + + + + + + + + + + + + + + + + + + +    the @xmath34-cube uso resulting from ( [ eq : morrismq ] ) can be described in purely combinatorial terms ( * ? ? ? * lemma 1 ) . here we make its structure even more transparent by exhibiting a _ finite state transducer _ ( finite automaton with output ) with just two states that can be used to describe the uso , see figure [ fig : morrisfst ] . for the remainder of this section , we fix @xmath34 to be an odd integer and we are considering the @xmath34-dimensional morris orientation  @xmath56 .    ] the orientation is then determined for each vertex @xmath121 as follows . if @xmath122 , then @xmath41  is the global sink . otherwise choose some  @xmath14 in  @xmath35 $ ] so that @xmath123 and consider the bit string @xmath124 as the input string to the transducer . in other words , the transducer is reading the input from right to left , starting immediately to the left of some @xmath116 . the output string @xmath125 then determines the orientation of each edge incident to  @xmath41 . the choice of  @xmath14 does not matter , because after reading any zero , the transducer is in the state  @xmath126 .    for instance , let @xmath127 and let @xmath128 . then the transducer takes @xmath129 , outputs @xmath54 and switches to state  @xmath130 ; reads @xmath131 , outputs @xmath53 and switches to  @xmath126 ; reads @xmath132 , outputs @xmath53 and stays in  @xmath126 ; reads @xmath133 , outputs @xmath54 and switches to  @xmath130 ; and finally reads @xmath134 , outputs  @xmath54 and switches to  @xmath126 . the resulting configuration for  @xmath41 is @xmath135 . figure [ fig : uso ] depicts the morris orientation for @xmath86 . it can easily be checked that the above procedural definition of the orientation is equivalent to the one given in ( * ? ? ? * lemma 1 ) . in particular , we obtain a uso . it has the remarkable symmetry that the set of configurations is closed under cyclic shifts . [ [ pivoting . ] ] pivoting . + + + + + + + + +    from the transducer we can easily derive the sign changes in the configuration  @xmath136 that are caused by a pivot step @xmath137 . we always have @xmath138 and @xmath139 . if @xmath140 , we get @xmath141 . in processing @xmath41 and @xmath142 , the transducer then performs exactly the same steps , except that at the @xmath14th coordinate different transitions are used to get from  @xmath130 to  @xmath126 . but this implies @xmath143 for @xmath144 . here is an example for such a pivot step ( with @xmath145 ; affected signs are indicated ) : @xmath146    if @xmath123 , we get @xmath147 , meaning that at the @xmath14th coordinate the transducer stays in  @xmath126 for  @xmath41 but switches to  @xmath130 for  @xmath142 ( assuming that @xmath142  is not the sink ) . this implies that signs change at all coordinates @xmath148 , @xmath149 down to ( and including ) the next coordinate  @xmath106 ( wrap around possible ) for which @xmath150 . in our example we have such a step for @xmath151 : @xmath152 in both @xmath41 and @xmath142 , the signs in the block of @xmath103 s at indices @xmath153 alternate .    [ [ levels . ] ] levels . + + + + + + +    it will be useful to define the following function on the vertices of the cube . let @xmath154 be the number of coordinates where @xmath155 appears . formally , @xmath156 the number @xmath154 is called the _ level _ of  @xmath41 . from the two types of pivots , it is easy to see that the value of  @xmath157 does not increase along any directed edge  ( * ? ? * lemma  3 ) . indeed , if @xmath142  is an out - neighbor of  @xmath41 , then either @xmath158 ( this happens in every pivot of type , and in pivots of type with an odd block of @xmath103 s at indices @xmath153 ) , or @xmath157 decreases to the next possible value : @xmath159 , or @xmath160 if @xmath161 . in particular , the values of  @xmath157 lie in the set @xmath162 . the sink is the only vertex at level  0 . let us make a small digression and briefly explain why randomedge is slow on  @xmath52 . let @xmath163 . so the outdegree of  @xmath41 is @xmath164 . now let @xmath41 be a vertex at level  @xmath103 . in pivots of type , @xmath165  decreases by one , whereas in pivots of type that do not reach the sink , @xmath165  increases by one . the latter occurs if and only if randomedge pivots at the unique @xmath155 , which happens with probability @xmath166 .    at level  @xmath103 , we can thus interpret randomedge as a random walk on the integers , where @xmath167 with probability @xmath168 , and @xmath169 with probability @xmath170 . therefore the walk is strongly biased towards  @xmath116 , but in order to reach a neighbor  @xmath142 of the sink , we need to get to @xmath171 . this takes exponentially long in expectation  @xcite . [ [ a - quadratic - bound - for - textscmurty_pi . ] ] a quadratic bound for @xmath100 . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    to prove theorem  [ thm : general ] we first derive a ( somewhat surprising ) more general result : for _ any _ pivot rule , starting from any vertex @xmath41 with @xmath172 for some @xmath106 , it takes @xmath112 iterations to reach a vertex @xmath142 with @xmath173 . this can be viewed as a statement about the uso induced by the @xmath34-dimensional morris orientation on the cube facet @xmath174 . the statement is that in this induced uso , all directed paths are short ( in particular , the induced uso is acyclic ) . [ lem : poten ] let @xmath36 be a vertex , and let @xmath106 be any index for which @xmath175 . starting from  @xmath41 , algorithm [ alg : pivot ] ( with an arbitrary pivot rule r ) reaches a vertex  @xmath142 satisfying @xmath176 after at most @xmath177 iterations . we employ a nonnegative potential that decreases with every pivot step . by a cyclic shift of coordinates we may assume that @xmath178 , and thus @xmath179 . let us define @xmath180 : v_j = 0 \text { and } \phi(v)_j={-}\},\\ n_1(v ) & : = \{j\in[n ] : v_j = 1 \text { and } \phi(v)_j={-}\}.    \end{aligned}\ ] ] the _ potential _ of @xmath41 is the number @xmath181    now we will observe how the potential changes during a pivot step @xmath182 . we always have @xmath138 , @xmath139 . there are two cases . if @xmath140 , as in , the configurations of @xmath41 and @xmath142 differ only at index @xmath14 . hence @xmath183 and @xmath184 . therefore @xmath185 .    in the other case we have @xmath123 , as in . if @xmath186 , we are done , since then @xmath187 . otherwise , @xmath41 is of the form @xmath188 with @xmath189 . if @xmath190 is odd , there is an even number of @xmath103 s between @xmath191 and @xmath192 , implying that @xmath193 and @xmath194 . we then have @xmath195 . if @xmath190 is even , there is an odd number of @xmath103 s between @xmath191 and @xmath192 , implying that @xmath196 and @xmath197 . since @xmath198 in this case , we get @xmath199 .    to summarize , the potential decreases by at least 1 in every pivot step @xmath200 , and as long as @xmath41 is not the sink , we have @xmath201 . the highest potential is attained by the source , @xmath202 . therefore after at most @xmath177 steps we have @xmath133 . if @xmath41 happens to be at level  1 already , we obtain a better bound , since @xmath203 in this case . in fact , the vertex with largest potential in level  1 such that @xmath179 is @xmath204 , and its potential is @xmath205 . [ cor : poten ] let @xmath36 be a vertex at level  1 , and let @xmath106 be any index for which @xmath175 . starting from  @xmath41 , algorithm [ alg : pivot ] ( with arbitrary pivot rule r ) reaches a vertex  @xmath142 satisfying @xmath176 after at most @xmath206 iterations . now @xmath100 comes in : we will apply either lemma  [ lem : poten ] or corollary  [ cor : poten ] to vertices @xmath41 with @xmath207 and @xmath208 , for some @xmath14 . the next lemma shows that the @xmath103 s at coordinates @xmath209 will stay . [ lem : visonly ] let @xmath210 $ ] . after algorithm [ alg : pivot ] with @xmath211 visits a vertex  @xmath41 with @xmath208 , it will visit only vertices  @xmath142 satisfying @xmath212 .    at every iteration , the algorithm replaces the current vertex @xmath41 with @xmath213 , where @xmath214 is the smallest index such that @xmath215 . it follows that coordinates @xmath209 will be touched only if @xmath216 for all @xmath217 . when this holds for the first time , we have @xmath218    the unique - completion property ( lemma [ lem : ucp ] ) then implies that we have already reached the sink . now we can put it all together . let @xmath41 be the start vertex , and let @xmath14 be the largest index such that @xmath207 ( if no such index exists , @xmath41 is the sink , and we are done ) . lemmas  [ lem : poten ] and [ lem : visonly ] show that after at most @xmath219 iterations , we additionally have @xmath220 ( since the next pivot step flips @xmath221 ) . at this point @xmath41 is at level  1 , since there is a unique  @xmath155 : at coordinate  @xmath96 . after the flip at coordinate  @xmath96 , we repeatedly apply corollary  [ cor : poten ] together with lemma [ lem : visonly ] to also produce @xmath103 s at coordinates @xmath222 , at which point we have reached the sink . this takes at most @xmath223 more iterations , summing up to a total of @xmath224    a better bound in the case of the identity permutation can be achieved by thoroughly examining the run of the algorithm , as we do next . [ [ the - identity - permutation . ] ] the identity permutation . + + + + + + + + + + + + + + + + + + + + + + + + +    let us now go on to prove theorem  [ thm : id ] . for this , we first identify certain `` milestone '' vertices that are visited by algorithm [ alg : pivot ] with @xmath117 , and then we count the number of iterations to get from one milestone to the next . let us define @xmath225    [ lem : l4 ] if algorithm [ alg : pivot ] with @xmath117 starts at the vertex @xmath226 , @xmath227 , it gets to the vertex @xmath228 in @xmath229 iterations . the algorithm starts at @xmath230 , and then proceeds by the rules of pivoting and in four phases ( recall that @xmath231 always pivots on the leftmost  @xmath53 ) :    * it pivots on the coordinates @xmath232 , @xmath233 ,   ,  @xmath103 and thus reaches the vertex + @xmath234 after @xmath235  iterations . * it pivots on the coordinates @xmath236 , @xmath120 , @xmath237 ,   ,  @xmath238 , @xmath239 and thus reaches the vertex + @xmath240 after @xmath235  iterations ; * it pivots on the coordinates @xmath238 , @xmath241 ,   ,  @xmath236 and thus reaches the vertex + @xmath242 after @xmath14  iterations ; * it pivots on the coordinates @xmath103 , @xmath243 , @xmath244 ,   ,  @xmath232 and thus reaches the vertex + @xmath245 after @xmath235  iterations . therefore it takes @xmath229 iterations to get from @xmath226 to @xmath228 . the previous lemma allows us to count the number of iterations from @xmath246 to @xmath247 . the following lemma takes care of the remaining iterations . [ lem : l6 ] if algorithm [ alg : pivot ] with @xmath117 starts at the vertex @xmath247 , it gets to the global sink in @xmath248 iterations . the proof is similar to the proof of lemma [ lem : l4 ] , except that here only the first phase takes place : if the algorithm starts at the vertex @xmath249 , it pivots on coordinates @xmath34 , @xmath250 ,   ,  @xmath103 and thus reaches the sink after @xmath248 iterations . by lemmas  [ lem : l4 ] and [ lem : l6 ] , the overall number of iterations from the source to the sink is @xmath251 in this section we examine _ k - lcps _ , the linear complementarity problems @xmath110 where @xmath5  is a k - matrix . we introduce two simple combinatorial conditions on unique - sink orientations and prove that one of them implies that all directed paths from  @xmath116 are short , and the other one implies that all directed paths in the orientation are short . we show that k - lcp orientations satisfy these conditions and conclude that a simple principal pivoting method with an arbitrary pivot rule solves a k - lcp in linearly many iterations , starting from any vertex . a unique - sink orientation  @xmath52 is _ 2-up - uniform _ if whenever @xmath253 is such that @xmath254 and @xmath142  is the source of  @xmath255 , then the orientation of the subcube  @xmath74 is uniform ( see figure  [ fig:2uu ] ) . the essential fact we use is that a subcube of a k - uso is also a k - uso , as was observed already by stickney and watson  ( * ? ? ? * lemma  1 ) . thus it suffices to prove 2-up - uniformity for k - lcps of dimension  2 . hence suppose that @xmath257 is a @xmath258 k - matrix , so @xmath259 and @xmath260 . then @xmath261 and nonsingular . now 0  is the source of the orientation induced by the  @xmath110 if and only if @xmath262 . hence @xmath263 , which proves uniformity . one particular nice property of a 2-up - uniform orientation is that all paths from the vertex  0 to the sink  @xmath107 have the shortest possible length , equal to the hamming distance @xmath264 of  @xmath107 from  @xmath116 . to prove this statement , we use the following lemma . suppose for the sake of contradiction that @xmath140 , @xmath266 and @xmath269 , and that @xmath267 , @xmath270 and @xmath271 . let @xmath272 and consider the subcube  @xmath74 ( see figure  [ fig : uni - proof1 ] , left ) : we observe that @xmath273 because otherwise the subcube would not contain a sink . but then it follows from 2-up - uniformity of  @xmath52 that @xmath274 , a contradiction . it follows from lemma  [ lem : crufa - up ] by induction that if @xmath142  is a vertex on any directed path starting from  @xmath116 , and if @xmath147 , then @xmath274 . hence all edges on a directed path from  @xmath116 to  @xmath107 are oriented from  0 to  1 , and therefore the length of a directed path from  @xmath116 to  @xmath107 is  @xmath264 .        what is the actual strength of this theorem ? we know that from any vertex of a unique - sink orientation there _ exists _ a path to the sink of length at most  @xmath34 ( lemma [ lem : hamming ] ) . it has also long been known how to find such a path from the vertex  0 in the case of a k - uso  @xcite . the novelty is that _ any _ directed path starting from  0 reaches the sink after the least possible number of iterations . hence a simple principal pivoting method with an _ arbitrary _ pivot rule finds the solution to a k - lcp in at most @xmath34  iterations .    in view of this result , it is natural to ask whether it also holds if the path starts in some vertex different from  0 . imposing 2-up - uniformity is insufficient , as we can observe by looking at the morris orientations . if we swap 0 s and 1 s in any morris orientation , we get a 2-up - uniform orientation . hence all directed paths from  0 to the sink are short ; indeed , 0  is the sink of the orientation after swapping , so all such paths have length  0 . nevertheless , long directed paths ( and even directed cycles ) do exist in this orientation . therefore we introduce the following stronger combinatorial property of usos , which turns out to be satisfied by k - usos as well . [ cor:2n ] algorithm [ alg : pivot ] with an arbitrary pivot rule  r , starting at an arbitrary vertex of the corresponding uso , finds the solution to an @xmath34-dimensional k - lcp in at most @xmath280  iterations .    as a consequence , we get a result for a larger class of lcps . if @xmath5  is a _ principal pivotal transform _ * section 2.3 ) of a k - matrix  @xmath281 , then @xmath5  is a p - matrix ( * ? ? ? * theorem 4.1.3 ) , and the uso @xmath52 induced by the @xmath110 is isomorphic to the uso  @xmath282 induced by the @xmath283 for suitable @xmath284 . it follows that corollary  [ cor:2n ] also applies to the @xmath110 , even though @xmath5  is not necessarily a k - matrix itself . our proof of theorem  [ thm:2lu ] is based on the following crucial fact , which extends lemma  [ lem : crufa - up ] . informally , it asserts that once we have a @xmath285 in some coordinate , we will always have a @xmath285 . if @xmath266 , the claim follows from lemma  [ lem : crufa - up ] . so let us suppose that we have @xmath140 , @xmath286 and @xmath269 , and that @xmath267 , @xmath270 and @xmath287 . let @xmath288 . we observe that @xmath273 so that the subcube  @xmath74 contains a sink . but then it follows from 2-up - uniformity of  @xmath278)}$ ] that @xmath64 , a contradiction ( see figure  [ fig : uni - proof1 ] , right ) . let @xmath14 be a fixed coordinate . suppose in the considered directed path there is an edge @xmath269 such that @xmath140 . as a consequence of lemma  [ lem : crufa ] , all vertices @xmath142 on a directed path starting at @xmath41 satisfy @xmath147 . it follows that in any directed path no more than two edges appear of the form @xmath289 for any fixed  @xmath14 : possibly one with @xmath123 and one with @xmath140 . this fact then implies that the length of any directed path is at most  @xmath280 . the strength of 2-uniformity can perhaps be explained by its being equivalent to what we call _ local uniformity_. a uso  @xmath52 is _ locally up - uniform _ if for every @xmath290 such that @xmath291 and @xmath292 for all @xmath293 , the orientation of the subcube  @xmath74 is uniform . a uso  @xmath52 is _ locally uniform _ if both @xmath52 and @xmath278)}$ ] are locally up - uniform . local up - uniformity allows for a slight improvement in solving a k - lcp . when following a directed path starting in  0 , instead of traversing one edge at a time we can perform a `` greedy '' principal pivot : jump straight to the sink of the subcube induced by outgoing edges . unlike the general p - matrix case , it can not lead to cycling . due to local up - uniformity , it will never increase the number of iterations and may sometimes reduce it . the answer to this question is negative . it can be shown that the number of @xmath34-dimensional locally uniform usos is asymptotically much larger than the number of k - usos . a lower bound of @xmath294 on the number of locally uniform usos can be derived by an adaptation of develin s construction  @xcite of many orientations satisfying the holt - klee condition . an upper bound of  @xmath295 on the number of p - usos ( and thus also k - usos ) follows from develin s proof of his upper bound on the number of lp - realizable cube orientations . a.  borici and h .- j . pricing american put options by fast solutions of the linear complementarity problem . in e. j. kontoghiorghes , b.  rustem , and s.  siokos , editors , _ computational methods in decision - making , economics and finance _ , volume  74 of _ applied optimization_. kluwer acad . dordrecht , 2002 . f.  holt and v.  klee . a proof of the strict monotone @xmath120-step conjecture . in _ advances in discrete and computational geometry _ , volume 223 of _ contemp . _ , pages 201216 .
Abstract:

Answer: we study the behavior of simple principal pivoting methods for the p - matrix linear complementarity problem ( p - lcp ) . we solve an open problem of morris by showing that murty s least - index pivot rule ( under any fixed index order ) leads to a quadratic number of iterations on morris s highly cyclic p - lcp examples . we then show that on k - matrix lcp instances , _ all _ pivot rules require only a linear number of iterations . as the main tool , we employ _ unique - sink orientations _ of cubes , a useful combinatorial abstraction of the p - lcp .


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: the primordial cosmic microwave background ( cmb ) was generated when photons first decoupled from the baryonic fluid when the universe was only 400,000 years old . the vast majority of these photons travel unperturbed to the present day , and features of their angular power spectrum such as acoustic peaks and the damping tail @xcite record valuable information about cosmological parameters @xcite . baryons and dark matter evolve from small inhomogeneities at decoupling into increasingly complicated large - scale structure which can subtly perturb the observed pattern of cmb anisotropies . assuming that the primordial cmb is gaussian , non - gaussian correlations in the observed map can be used to reconstruct the intervening large - scale structure @xcite . in addition to the importance of learning about the large - scale structure itself , reconstruction of the weak - lensing potential generated by structure is essential to constraining tensor perturbations . weak lensing converts a fraction of the e - mode polarization generated by scalar perturbations at the last - scattering surface into b - mode polarization in the observed map . only by subtracting this b - mode polarization can one conclusively detect the primordial b - modes which serve as a model - independent signal of tensor perturbations @xcite . understanding lensing reconstruction requires a more detailed discussion of how weak lensing affects the cmb . weak gravitational lensing deflects the paths of cmb photons as they travel from the last - scattering surface to the observer . this deflection is accomplished by a projected lensing potential which is a weighted line - of - sight integral of the gravitational potential between the observer and the surface of last - scattering . at each point on the sky , lensing remaps the temperature and polarization to that of a nearby point at the last - scattering surface , the deflection angle being the gradient of the aforementioned projected lensing potential . assuming that this deflection angle is small , the temperature at any point can be expanded in a taylor series in the gradient of the lensing potential . in fourier space , this expansion appears as a series of convolutions of individual temperature and projected potential modes . the observed temperature - squared map in fourier space also appears as a convolution of individual fourier modes . subject to an overall normalization dependent on the scale of the fourier mode , these convolutions cancel in such a manner that each fourier mode of the temperature - squared map acts as an estimator for the _ same _ fourier mode of the projected lensing potential . lensing reconstruction as outlined above has been considered previously @xcite . in these works , two sources of noise were identified , and a filter of the temperature - squared map in fourier space was chosen to minimize the variance associated with lensing reconstruction subject to these noise sources . the first source is intrinsic signal variance ; the observed large - scale structure is one arbitrary member of an ensemble of realizations allowed by theory . the second source of noise , endemic to this method of lensing reconstruction , is a consequence of the nature of the primordial cmb . like the large - scale structure itself , the pattern of cmb anisotropies at the last - scattering surface is only one of many possible realizations allowed by theory . we do not know _ a priori _ which of these realizations nature has provided us , and this uncertainty hinders our ability to deconvolve the effects of lensing from true anisotropies at the last - scattering surface . even if the true pattern of anisotropies at the last - scattering surface was known , the finite amount of power in the cmb at small scales would still constrain lensing reconstruction . silk damping at the last - scattering surface suppresses cmb power at small scales , while the finite resolution of any real experiment would limit the detection of any signal that is present at small scales . lensing reconstruction fails below scales at which there is sufficient power , for the same reason that any remapping is indistinguishable given a uniform background .    here , we consider a third source of noise neglected in previous studies . the filtered temperature - squared map is an unbiased estimator for the lensing potential in the approximation that a correlation between two given temperature modes is induced only by the single lensing mode whose wavevector is the sum of that of the two temperature modes . in actuality , any combination of two or more lensing modes whose wavevectors sum to this total induce correlations between the two temperature modes . there are many such combinations , but since these correlations add incoherently we do not expect a systematic bias . nonetheless , for estimators of each individual lensing mode we must use our knowledge of other lensing modes to subtract off this unwanted bias . this is an iterative process , and since our knowledge of the lensing map is imperfect it induces noise in lensing reconstruction . we calculate this additional variance for various estimators constructed from cmb temperature and polarization maps , and show how it compares to the dominant noise sources for an all - sky cmb experiment with a noise - equivalent temperature of 1 @xmath1k @xmath2 . since the lensing - potential power spectrum is a measure of the theoretical uncertainty with which we can predict the value of a given lensing mode , this noise associated with lensing reconstruction causes a systematic overestimation of the lensing - potential power spectrum . this systematic bias must be accounted for in order to compare observations with theoretical predictions . this paper is organized as follows . in  [ s : lensing ] , we define the formalism we will use to explore the effects of weak lensing on the cmb . the taylor expansion of the lensed cmb map in gradients of the lensing potential is given in both real and fourier space , and the power spectra and trispectra of various components of the cmb temperature map are listed for later use . in  [ s : estimators ] , we show that the fourier modes of the temperature - squared map when properly filtered can serve as estimators for the fourier modes of the displacement map with the same wave vector . using the power spectrum and trispectrum given in the preceding section , we calculate the variance associated with this estimator including a new component neglected in previous studies . this variance is evaluated numerically using the currently favored @xmath0cdm cosmological model with baryon density @xmath3 , matter density @xmath4 , cosmological constant density @xmath5 , hubble parameter @xmath6=0.65 , and power - spectrum amplitude @xmath7 . we then use the displacement estimator for individual fourier modes to construct an unbiased estimator for the lensing - potential power spectrum in   [ s : power ] , and calculate the variance and covariance associated with this estimator . a few concluding remarks about the implications of our work for future studies are given in  [ s : disc ] . the appendix contains useful formulae related to additional estimators of lensing based on polarization and a combination of temperature and polarization . we consider weak lensing under the flat - sky approximation following refs . @xcite . as discussed before @xcite , weak lensing deflects the path of cmb photons resulting in a remapping of the observed temperature pattern on the sky , @xmath8 \nonumber\\          & \approx & { \theta}({\hat{\bf n } } ) + \nabla_i { \phi}({\hat{\bf n } } ) \nabla^i { \theta}({\hat{\bf n } } ) + \frac{1}{2 } \nabla_i { \phi}({\hat{\bf n } } ) \nabla_j { \phi}({\hat{\bf n } } ) \nabla^{i}\nabla^{j } { \theta}({\hat{\bf n } } ) + \ldots\end{aligned}\ ] ] where @xmath9 is the unlensed primary component of the cmb in a direction @xmath10 at the last scattering surface . the observed , gravitationally - lensed temperature map @xmath11 in direction @xmath10 is that of the unlensed map in direction @xmath12 where @xmath13 represents the lensing deflection angle or displacement map . although a real cmb map will include secondary contributions such as the sz effect @xcite , we assume that such effects can be distinguished by their frequency dependence @xcite . they will not be further considered in this paper . a noise component denoted by @xmath14 due to finite experimental sensitivity must be included as well . thus the total observed cmb anisotropy will be @xmath15 . taking the fourier transform of the lensed map @xmath11 under the flat - sky approximation , we write @xmath16 where @xmath17 + \frac{1}{2 } { \int \frac{d^2 { \bf l''}}{(2\pi)^2}}{\phi}({{\mathbf{l^{\prime\prime } } } } ) \\ & & \quad \times { \phi}({{\mathbf{l}}}- { { \mathbf{l^{\prime}}}}- { { \mathbf{l^{\prime\prime } } } } ) \ , ( { { \mathbf{l^{\prime\prime}}}}\cdot { { \mathbf{l^{\prime } } } } )                  \left [ ( { { \mathbf{l^{\prime\prime}}}}+ { { \mathbf{l^{\prime}}}}- { { \mathbf{l}}})\cdot                               { { \mathbf{l^{\prime}}}}\right ] + \ldots \ , .   \nonumber\end{aligned}\ ] ]    cmb correlations in fourier space can be described in terms of a power spectrum and trispectrum as defined in the usual manner , @xmath18 where the angle brackets denote ensemble averages over possible realizations of the primordial cmb , large - scale structure between the observer and the surface of last - scattering , and instrumental noise . the subscript @xmath19 denotes the connected part of the four - point function and the superscript @xmath20 denotes the temperature map being considered ( @xmath21 or @xmath22 ) . the lensing - potential power spectrum can be defined analogously , @xmath23 where here the angle brackets denote an average over all realizations of the large - scale structure . we make the assumption that primordial fluctuations at the last - scattering surface are gaussian . gaussian statistics are fully described by a power spectrum ; the gaussian four - point correlator , @xmath24 is zero . the instrumental noise @xmath22 is also assumed to be gaussian , as is the lensing potential @xmath25 . this second assumption is justified because the dominant contributions to the lensing potential come from intermediate redshifts @xmath26 at which linear theory holds . using these definitions , we can calculate the anticipated power spectrum and trispectrum of the observed cmb map . because the instrumental noise is uncorrelated with the signal , the power spectrum of the observed map is the sum of signal and noise power spectra , @xmath27 the power spectrum of the noise component is given by : @xmath28 where @xmath29 is the fraction of the sky surveyed , @xmath30 is the variance per unit area on the sky , and @xmath31 is the effective beamwidth of the instrument expressed in terms of its full - width half - maximum resolution @xmath32 . a cmb experiment that spends a time @xmath33 examining each of @xmath34 pixels with detectors of sensitivity @xmath35 will have a variance per unit area @xmath36 @xcite . the power spectrum of the lensed cmb can be determined by inserting eq . ( [ e : thetal ] ) into eq . ( [ e : specdef ] ) as discussed in @xcite , @xmath37    \ ,                               c_l^{{\theta}{\theta } }          + { \int \frac{d^2 { \bf l}_1}{(2\pi)^2 } } c_{| { { \bf l}}- { { \bf l}}_1|}^{{\theta}{\theta } } c^{\phi\phi}_{l_1 }                  [ ( { { \bf l}}- { { \bf l}}_1)\cdot { { \bf l}}_1]^2   \ , . \label{e : lenpow}\end{aligned}\ ] ] this result is given to linear order in the lensing - potential power spectrum @xmath38 . lensing neither creates nor destroys power in the cmb , but merely shifts the scales on which it occurs as seen by the fact that @xmath39 the observed cmb trispectrum can be calculated in a similar manner ; under our assumptions of gaussian instrumental noise and no secondary anisotropies the trispectrum of the lensed component @xmath40 is the sole contribution to the total observed trispectrum , @xmath41 [ ( { { { \bf l}_1}}+ { { { \bf l}_3 } } ) \cdot { { { \bf l}_4 } } ]     + c^{{\phi\phi}}_{|{{\bf l}}_2+{{\bf l}}_3| } [ ( { { { \bf l}_2}}+{{{\bf l}_3 } } ) \cdot { { { \bf l}_3 } } ] [ ( { { { \bf l}_2}}+{{{\bf l}_3 } } ) \cdot { { { \bf l}_4 } } ] \big ] + \ , { \rm perm . } \ , . \nonumber \\ \label{e : trilens}\end{aligned}\ ] ] the term shown above is manifestly symmetric under the interchange @xmath42 , while the `` + perm . '' represents five additional terms identical in form but with the replacement of @xmath43 and @xmath44 with the other five combinations of pairs . the total trispectrum is symmetric under the interchange of any given pair as one would expect . having established a formalism within which to analyze weak lensing , we now consider the problem of reconstructing the lensing potential from an observed cmb temperature map . in this section , we examine lensing reconstruction following the approach of ref . @xcite , largely adopting their notation as well . the only important difference in notation is that we use @xmath40 to denote the lensed temperature field and @xmath45 for the unlensed field following ref . @xcite and most recent papers . ref . @xcite uses the opposite convention . for @xmath46 and to linear order in @xmath25 , @xmath47 where @xmath48 and @xmath49 . note that @xmath50 differs from the unmarked @xmath51 that first appeared in eq . ( [ e : specdef ] ) in that it denotes an ensemble average only over different gaussian realizations of the primordial cmb and instrument noise ; a fixed realization of the large - scale structure is assumed . for the purposes of estimating the large - scale structure actually realized in our observable universe , this is the appropriate average to take to ensure that our estimators are truly unbiased for a typical realization of the primordial cmb . when calculating the noise associated with lensing - potential estimators and again for power spectrum estimation in   [ s : power ] , we will return to the full unmarked ensemble average . ( [ e : ea ] ) , an immediate consquence of eq . ( [ e : thetal ] ) , suggests that a temperature - squared map appropriately filtered in fourier space can serve as an estimator for the deflection field @xmath52 . hu and okamoto define five different estimators for the deflection field constructed from various combinations of the temperature and polarization ; we discuss the temperature - squared estimator in this section and relegate the analogous formulae for polarization estimators to the appendix . the minimum - variance temperature - squared estimator derived in ref . @xcite is @xmath53 where @xmath54 @xmath55^{-1 } \ , , \ ] ] and @xmath56 . substitution of eqs . ( [ e : ea ] ) , ( [ e : bias ] ) , ( [ e : filt ] ) , and ( [ e : norm ] ) into eq . ( [ e : est ] ) shows the desired result , @xmath57 namely that @xmath58 is indeed an unbiased estimator for the deflection field in fourier space . we now proceed to calculate the variance of this estimator . at first , we assume a complete knowledge of all lensing modes not examined by this estimator . in that case , we find @xmath59 where @xmath60 . evaluating the four - point function in the integrand of eq . ( [ e : estvar , cmb ] ) to second order in the lensing field , we obtain @xmath61 \nonumber \\ & & \quad \times \bigl [ { { \mathbf{l^{\prime}}}}\cdot ( { { { \bf l}_2}}- { { \mathbf{l^{\prime } } } } ) \bigr ] - \frac{1}{2 } { \int \frac{d^2 { \bf l'}}{(2\pi)^2}}{\phi}({{\mathbf{l^{\prime } } } } ) { \phi}(-{{\mathbf{l}}}- { { \mathbf{l^{\prime } } } } )   \bigl\ { c_{l_1}^{{\theta}{\theta } } ( { { { \bf l}_1}}\cdot { { \mathbf{l^{\prime } } } } ) \bigl [ { { { \bf l}_1}}\cdot ( { { \mathbf{l}}}+ { { \mathbf{l^{\prime } } } } ) \bigr ] + c_{l_2}^{{\theta}{\theta } } ( { { { \bf l}_2}}\cdot { { \mathbf{l^{\prime } } } } ) \bigl [ { { { \bf l}_2}}\cdot ( { { \mathbf{l}}}+ { { \mathbf{l^{\prime } } } } ) \bigr ] \bigr\ } \biggr ] \nonumber \\ & & \quad \times \biggl [ \bigl ( c_{l_{1}^{\prime}}^{{\theta}{\theta } } + c_{l_{1}^{\prime}}^{{\theta}{\theta}{{\rm n } } } \bigr ) ( 2 \pi)^2 \delta_{{\rm d}}({{\mathbf{l^{\prime } } } } ) + { \phi}({{\mathbf{l^{\prime } } } } ) { f_{\theta\theta}}({{{\bf l}_1 } } ' , { { { \bf l}_2 } } ' ) - { \int \frac{d^2 { \bf l'}}{(2\pi)^2}}c_{l^{\prime}}^{{\theta}{\theta } } { \phi}({{{\bf l}_1 } } ' - { { \mathbf{l^{\prime } } } } ) { \phi}({{{\bf l}_2 } } ' + { { \mathbf{l^{\prime } } } } ) \bigl [ { { \mathbf{l^{\prime}}}}\cdot ( { { { \bf l}_1 } } ' - { { \mathbf{l^{\prime } } } } ) \bigr ] \nonumber \\ & & \quad \times \bigl [ { { \mathbf{l^{\prime}}}}\cdot ( { { { \bf l}_2 } } ' + { { \mathbf{l^{\prime } } } } ) \bigr ] + \frac{1}{2 } { \int \frac{d^2 { \bf l'}}{(2\pi)^2}}{\phi}({{\mathbf{l^{\prime } } } } ) { \phi}({{\mathbf{l^{\prime}}}}- { { \mathbf{l^{\prime } } } } ) \bigl\ { c_{l_{1}^{\prime}}^{{\theta}{\theta } } ( { { { \bf l}_1 } } ' \cdot { { \mathbf{l^{\prime } } } } ) \bigl [ { { { \bf l}_1 } } ' \cdot ( { { \mathbf{l^{\prime}}}}- { { \mathbf{l^{\prime } } } } ) \bigr ] + c_{l_{2}^{\prime}}^{{\theta}{\theta } } ( { { { \bf l}_2 } } ' \cdot { { \mathbf{l^{\prime } } } } ) \bigl [ { { { \bf l}_2 } } ' \cdot ( { { \mathbf{l^{\prime}}}}- { { \mathbf{l^{\prime } } } } ) \bigr ] \bigr\ } \biggr ] { \rm + perm . } \nonumber \\\end{aligned}\ ] ] the terms given explicitly in eq . ( [ e:4pt , cmb ] ) correspond to the correlations between @xmath62 and @xmath63 and those between @xmath64 and @xmath65 . the `` + perm . '' stands for two additional terms , identical in form , arising from the pairings @xmath66 and @xmath67 . this expression indicates how uncertainty in the cmb at the last - scattering surface propagates into uncertainty in lensing reconstruction for a particular realization @xmath68 of the large - scale structure . to linear order in eq . ( [ e:4pt , cmb ] ) , correlations between the modes @xmath69 and @xmath64 are induced by those lensing modes whose wavevectors are the sums of any pair of wavevectors of these four modes . these lensing modes are precisely those forming the diagonals of the quadrilaterals depicted in fig . [ f:2quad ] .    in practice , we do not know the large - scale structure between us and the last - scattering surface , so we assume a variance given by eq . ( [ e : lenpowdef ] ) with a model - dependent power spectrum @xmath70 . we must average eq . ( [ e : estvar , cmb ] ) over different realizations of the large - scale structure ( denoted by @xmath71 ) to obtain the total expected variance of our estimator , @xmath72 where @xmath73 the assumption that the lensing potential is gaussian imposes the constraint @xmath74 which closes the quadrilaterals of fig .  [ f:2quad ] . the average of the four - point correlation function in eq . ( [ e : estvar ] ) can be calculated by further averaging eq . ( [ e:4pt , cmb ] ) over the large - scale structure . terms linear in the lensing field vanish when averaged over different realizations of the large - scale structure . quadratic terms in the lensing field arise as products either of two linear terms or of a zeroth and second - order term . averages over the product of two linear terms produce the connected part of the four - point correlation function , the trispectrum defined in eq . ( [ e : specdef ] ) . averages over the product of a zeroth and second - order term have no connected portion , but instead furnish an implicit dependence on @xmath70 in the total observed power spectrum of eq . ( [ e : temptot ] ) . the final result of averaging over the large - scale structure can be expressed in terms of the observed power spectrum and trispectrum , @xmath75 \ , , \nonumber \\\end{aligned}\ ] ] where the trispectrum can written in terms of @xmath76 as @xmath77 this form of the trispectrum is consistent with that of eq . ( [ e : trilens ] ) given directly in terms of power spectra . since @xmath78 , @xmath79 and the first term of eq . ( [ e:4pt ] ) vanishes . the remaining two terms containing pairs of delta functions , inserted into eq . ( [ e : estvar ] ) , yield the dominant contribution to the variance , @xmath80 \ , , \ ] ] where @xmath81 . notice that @xmath82 is zeroth order in the lensing potential @xmath25 ; it depends on the lensing potential power spectrum @xmath83 only implicitly though the total observed power spectrum @xmath84 . the ellipsis represents terms of higher order in @xmath83 that we now proceed to calculate . these terms arise from the trispectrum term of eq . ( [ e:4pt ] ) after @xmath85 is removed . substituting these results into eq . ( [ e : estvar ] ) , we find that to first order in @xmath83 , @xmath86 \ , , \ ] ] where @xmath87 is given by , @xmath88 the first - order contribution to the noise @xmath87 involves integrals over the lensing - potential power spectrum , and thus probes lensing modes with wave vectors different from that of the estimator @xmath58 . it can be interpreted physically as interference from these other modes in the determination of the mode @xmath89 being estimated . the filter @xmath90 was chosen to optimize the signal - to - noise ratio in the absence of the first - order contribution @xmath87 ; it is no longer an optimal filter once this additional noise is taken into account . as long as @xmath91 , the noise reduction that can be attained by re - optimizing our filter will not be significant . formulas analagous to those presented here relevant to the construction of estimators using polarization data are given in the appendix . the significance of @xmath87 for two different experiments is shown in fig . [ f : plrec ] using the currently favored @xmath0cdm cosmological model with baryon density @xmath3 , matter density @xmath4 , cosmological constant density @xmath5 , hubble parameter @xmath6=0.65 , and power - spectrum amplitude @xmath7 . the planck experiment is equivalent to a one - year , full - sky survey with temperature and polarization sensitivities of 12.42 and 26.02 @xmath92 respectively and resolution @xmath93 arcminutes as described in   [ s : lensing ] . the reference experiment has the same resolution but superior sensitivities of 0.46 and 0.65 @xmath94 for temperature and polarization . these estimates of experimental parameters are identical to those given for the planck and reference experiments of ref . the @xmath95 and @xmath96 estimators have noise power spectra intermediate to those of the @xmath97 and @xmath98 estimators , while the @xmath99 estimator has substantially higher noise because the primordial cmb lacks true b - modes in the absence of inflationary gravitational waves . we see that for planck , with its comparatively inferior polarization sensitivity , the @xmath97 estimator will be best although it will be unable to detect individual fourier modes of the deflection field at the @xmath100 level . the reference experiment , and further experiments with similar sensitivity and even higher resolution , should be able to push @xmath100 detection of individual @xmath89 modes to @xmath101 by primarily relying on the @xmath98 estimator . in these cases the secondary noise @xmath102 is only smaller than the dominant noise @xmath103 by a factor of a few , whereas for higher sensitivity ( noisier ) experiments like planck it is smaller by at least an order - of - magnitude . this illustrates an interesting point , apparent from fig . [ f : plrec ] , that the zeroth - order noise @xmath82 declines dramatically with decreasing sensitivity until it becomes dominated by cosmic variance while the @xmath87 is largely unaffected by instrument sensitivity . the reasons for this trend are that instrument noise appears in @xmath81 through its contribution to the denominator of @xmath104 as shown by eqs . ( [ e : filt ] ) and ( [ e : norm ] ) . decreasing instrument noise raises the value of @xmath104 thereby lowering @xmath82 . by contrast instrument noise is reflected in @xmath87 through its effects on both @xmath105 and @xmath104 as shown in eq . ( [ e:1st ] ) . smaller instrument noise raises @xmath104 as before , driving @xmath87 up in this case , but this is compensated by a decrease in @xmath105 which appears as a prefactor outside the integrals . these two effects largely cancel each other out , rendering @xmath87 remarkably insensitive to instrument noise . although complete reconstruction of the deflection field @xmath89 can be an enormously powerful tool , such as for b - mode subtraction @xcite , for some purposes estimates of the lensing - potential power spectrum @xmath106 are sufficient . this power spectrum is a model - dependent prediction of theories of large - scale structure formation , and therefore estimates of the power spectrum from real data could be used to test these theories as well as the consistency of other determinations of cosmological parameters . furthermore , since estimates of all the modes @xmath89 with @xmath107 can be combined to estimate @xmath106 , @xmath100 detection of the power spectrum can be pushed to much higher @xmath108 than can that of individual modes . the deflection - field estimator @xmath58 derived in the preceding section can be used to construct an estimator for @xmath106 . our first guess for an appropriate lensing - potential power spectrum estimator is @xmath109 where @xmath110 is the area of the sky surveyed and @xmath111 is an annulus of radius @xmath108 and width @xmath112 . we ensemble average our estimator over different realizations of the cmb and large - scale structure using eq . ( [ e : prodave ] ) by bringing @xmath113 from the left to the right - hand side . this yields @xmath114 \ , .\ ] ] the definition of the dirac delta function , @xmath115 implies that @xmath116 . furthermore , in the limit that @xmath112 is small compared to the scales on which @xmath117 is varying , we can evaluate the integrand of eq . ( [ e : naive2 ] ) at its central value @xmath118 and extract it from the integral . the integral over the annulus @xmath111 cancels the factor @xmath119 in the denominator , reducing eq . ( [ e : naive2 ] ) to @xmath120 \ , .\ ] ] @xmath121 is indeed an estimator for the lensing - potential power spectrum @xmath106 , albeit a biaed one . note that the _ bias _ in the power - spectrum estimator @xmath121 is precisely the same as the _ variance _ shown in eq . ( [ e : prodave ] ) with which we were able to determine each individual lensing mode . this is no coincidence ; it reflects the fact that there are no grounds _ a priori _ on which to differentiate the variance with which we can reconstruct individual modes @xmath89 from the intrinsic variance @xmath122 of the underlying distribution from which they are drawn . to obtain an unbiased estimator to compare with theoretical predictions , we subtract off this unwanted reconstruction variance , @xmath123 \ , .\ ] ] since @xmath87 as defined in eq . ( [ e:1st ] ) itself depends on @xmath106 , this subtraction and evaluation must be performed iteratively until a self - consistent solution is obtained . the variance of our estimator @xmath124 can be calculated in the usual manner , @xmath125 evaluating this expression requires us to calculate @xmath126 \left [ { { \mathbf{d}}_{\theta\theta}}({{{\bf l}_2 } } ' ) \cdot { { \mathbf{d}}_{\theta\theta}}(-{{{\bf l}_2 } } ' ) \right ] \rangle \ , .\ ] ] since @xmath58 is a quadratic estimator in the temperature map , eq . ( [ e : dl2 ] ) includes the following integral over the eight - point correlation function in fourier space , @xmath127 \left [ { { \mathbf{d}}_{\theta\theta}}({{{\bf l}_2 } } ' ) \cdot { { \mathbf{d}}_{\theta\theta}}(-{{{\bf l}_2 } } ' ) \right ] \rangle = \frac{{a_{\theta\theta}}^2(l_{1}^{\prime } ) { a_{\theta\theta}}^2(l_{2}^{\prime})}{(l_{1}^{\prime})^2 ( l_{2}^{\prime})^2 } { \int \frac{d^2 { \bf k}_1}{(2\pi)^2 } } { \int \frac{d^2 { \bf k}_3}{(2\pi)^2 } } { \int \frac{d^2 { \bf k}_5}{(2\pi)^2 } } { \int \frac{d^2 { \bf k}_7}{(2\pi)^2 } } \nonumber \\ & & \quad \times \big\ { { f_{\theta\theta}}({{{\bf k}_1 } } , { { { \bf k}_2 } } ) { f_{\theta\theta}}({{{\bf k}_3 } } , { { { \bf k}_4 } } ) { f_{\theta\theta}}({{{\bf k}_5 } } , { { { \bf k}_6 } } ) { f_{\theta\theta}}({{{\bf k}_7 } } , { { { \bf k}_8 } } ) \langle { \theta}^{{\rm t}}({{{\bf k}_1 } } ) { \theta}^{{\rm t}}({{{\bf k}_2 } } ) { \theta}^{{\rm t}}({{{\bf k}_3 } } ) { \theta}^{{\rm t}}({{{\bf k}_4 } } ) { \theta}^{{\rm t}}({{{\bf k}_5 } } ) { \theta}^{{\rm t}}({{{\bf k}_6 } } ) { \theta}^{{\rm t}}({{{\bf k}_7 } } ) { \theta}^{{\rm t}}({{{\bf k}_8 } } ) \rangle \big\ } \ , , \nonumber \\\end{aligned}\ ] ] where @xmath128 , @xmath129 , @xmath130 , and @xmath131 . a fully general eight - point correlation function consists of a connected part , as well as terms proportional to the product of lower - order correlation functions . under the assumption that both the primordial cmb and the lensing potential are governed by gaussian statistics , all correlation functions higher than the four - point have vanishing connected parts @xcite . the temperature eight - point correlation function will therefore be composed of three groups of terms ; membership in a group being determined by whether the term contains zero , one , or two factors of the trispectrum . since the trispectrum given in eq . ( [ e : trilens ] ) is first order in the lensing - potential power spectrum @xmath106 , terms of these three groups are zeroth , first , and second order respectively in @xmath106 . combinatorics determines the number of terms in each group . there are : @xmath132 different ways of dividing @xmath133 into four pairs , and hence there will be 105 terms in the group containing no trispectra . similar calculations reveal that there are @xmath134 terms in the second group and @xmath135 terms in the third group . many terms in all three groups will vanish for the same reason that the first term vanished in the four - point correlation function of eq . ( [ e:4pt ] ) ; these terms are proportional to a dirac delta function evaluated at nonzero argument . consider now the first group of terms , those that are zeroth order in @xmath106 . the 60 nonvanishing terms in this group each contain four dirac delta functions ; they can be further segregated into the 12 terms that allow two of the integrals over @xmath136 appearing in eq . ( [ e:8pt ] ) to be immediately evaluated via dirac delta functions , and the 48 terms that allow evaluation of three @xmath136 integrals . the first 12 terms , inserted into eq . ( [ e:8pt ] ) and appropriately evaluated using the normalization of eq . ( [ e : norm ] ) , yield @xmath137 + \ldots \ , , \ ] ] while the remaining 48 terms give the final result to zeroth order in @xmath106 , @xmath138   \nonumber \\ & & \quad \quad \quad \quad \quad \quad + \frac{(2 \pi)^2}{a l^4 } \frac{2}{(2 \pi l \delta l)^2 } \int_{a_l } \frac{d^2 { { { \bf l}_1}}'}{(2 \pi)^2 } \ , \int_{a_l } \frac{d^2 { { { \bf l}_2}}'}{(2 \pi)^2 } \ , \frac{{a_{\theta\theta}}^2(l_{1}^{\prime } ) { a_{\theta\theta}}^2(l_{2}^{\prime})}{(l_{1}^{\prime})^2 ( l_{2}^{\prime})^2 } { \int \frac{d^2 { \bf k}_1}{(2\pi)^2 } } f_{{\theta}{\theta}}({{{\bf k}_1 } } , { { { \bf k}_2 } } ) p({{{\bf k}_1 } } , { { { \bf k}_2 } } , { { { \bf l}_1 } } ' , { { { \bf l}_2 } } ' ) \ , , \nonumber \\\end{aligned}\ ] ] where @xmath139 when @xmath140 is subtracted from @xmath141 in eq . ( [ e : sigpv ] ) , the first term in the square brackets of eq . ( [ e : zeropv ] ) will be eliminated . minimizing the variance associated with this estimator then consists of making an optimal choice of @xmath142 . the first noise term is proportional to @xmath143 . for a survey of area @xmath110 , @xmath144 is the specific area of an individual mode in @xmath145 space and @xmath146 is the area in @xmath145 space over which the power - spectrum estimator takes an average . this ratio is therefore the inverse of the number of individual @xmath58 modes whose inverse variances are added to determine the inverse variance of @xmath124 . it is obviously minimized by choosing @xmath147 . the second term , that involving @xmath148 , differs from the first noise term in that a dirac delta function has been used to evaluate an additional @xmath136 integral rather than an annulus integral . since the integrands are of the same order , we expect the second noise term to be suppressed relative to the first by a factor @xmath149 where @xmath150 is set by the resolution @xmath32 of the survey . under the conservative assumption @xmath151 , namely that we are probing scales well above our resolution , this term is assured to be small . we neglect such terms for the remainder of this paper . if we insert the portions of the eight - point correlation function that are first and second order in @xmath106 into eq . ( [ e:8pt ] ) and evaluate using eq . ( [ e:1st ] ) , we find @xmath152 \ , , \ ] ] and @xmath153 this result agrees with that given in ref . @xcite after subtracting our newly derived term @xmath87 . the term @xmath87 and corresponding terms for polarization - based estimators have two principal effects on power spectrum estimation . as shown in eq . ( [ e : finvar ] ) , they provide a fractional contribution to the variance of roughly @xmath154 when @xmath122 dominates the variance as in the reference experiment in the right - hand panel of fig . [ f : plrec ] . for the @xmath0cdm cosmological model considered here this represents an increase of @xmath155 in the variance of the @xmath98 estimator for @xmath156 . more importantly , @xmath87 acts as a bias for the naive estimator @xmath121 as shown by eq . ( [ e : naive3 ] ) . if this bias is not calculated and subtracted iteratively to form the unbiased estimator @xmath124 as in eq . ( [ e : pestsub ] ) , @xmath122 will be _ systematically _ overestimated by @xmath157 at low @xmath108 and by increasingly larger amounts at @xmath158 as the signal @xmath159 begins to plummet while @xmath87 remains comparatively flat . having evaluated the variance of our estimator @xmath124 , we consider whether this estimator has a substantial covariance @xmath160 the estimator @xmath121 as defined in eq . ( [ e : naive ] ) implies that : @xmath161 \left [ { { \mathbf{d}}_{\theta\theta}}({{{\bf l}_2 } } ' ) \cdot { { \mathbf{d}}_{\theta\theta}}(-{{{\bf l}_2 } } ' ) \right ] \rangle \ , , \ ] ] which can be evaluated using the same integral over the eight - point correlation function described in eq . ( [ e:8pt ] ) . whereas 60 of the 105 zeroth - order terms in @xmath106 coming from this equation were nonvanishing for the variance , for the covariance only 52 terms are nonzero provided that the widths @xmath112 and @xmath162 are chosen so that the annuli @xmath111 and @xmath163 do not overlap . this leads to a result analogous to eq . ( [ e : zeropv ] ) , @xmath164 note that the eight terms missing from the covariance when compared to the variance have altered the first term of eq . ( [ e : zerocv ] ) , and that it was precisely these terms that provided the dominant contribution to the variance of eq . ( [ e : finvar ] ) when @xmath142 was chosen appropriately . we therefore find that to zeroth order in @xmath106 , the covariance is given by @xmath165 where we have extracted the @xmath166 from the annular integrals since they are slowly varying over the widths @xmath112 and @xmath162 . for the same reasons that terms of this form were a subdominant contribution to the variance as discussed previously , we expect the covariance to be suppressed as well . if we define the ratio @xmath167 we can quantify this suppression . the triple integral of eq . ( [ e : fincov ] ) , appearing in the numerator of @xmath168 , involves integration over annuli with radii @xmath108 and @xmath169 and one integration over all fourier space . the triple integrals in the variances @xmath170 and @xmath171 appearing in the denominator of @xmath168 each consist of a single integration over an annulus of radius @xmath108 and @xmath169 respectively and two integrations over all fourier space . if we make the crude assumption that the integrand is constant , the ratio @xmath168 will simply be the ratio of these areas , @xmath172    the ratio @xmath168 is evaluated numerically for planck in fig . [ f : cov ] as a function of @xmath108 for various fixed values of @xmath169 . the estimators @xmath124 and @xmath173 were chosen such that @xmath174 , while integrals over fourier space were cut off at @xmath175 . substituting these values into eq . ( [ e : rcovest ] ) , we expect @xmath176 . this crude estimate is surprisingly close to the numerically obtained results of fig . [ f : cov ] ; in particular the slope of the curves is approximately 1/2 on this log - log plot . even at @xmath101 , @xmath177 suggesting that covariance in power - spectrum estimation can safely be neglected for planck . the estimate of eq . ( [ e : rcovest ] ) implicitly depends on the experimental resolution @xmath32 because the integrand appearing in expressions for the variance and covariance decreases rapidly for @xmath178 . for future experiments with better resolution than planck , @xmath179 will be higher implying by eq . ( [ e : rcovest ] ) that covariance will be even more negligible . weak gravitational lensing induces non - gaussian correlations between modes of the observed cmb temperature map as shown in eq . ( [ e : ea ] ) . these correlations , and assumptions about the gaussian nature of the primordial cmb , can be used to construct several temperature and polarization - based estimators of the fourier modes @xmath89 of the deflection field . this procedure was outlined in ref . @xcite , however in calculating the noise associated with this reconstruction , an assumption was made that the observed temperature map was gaussian . in the presence of lensing this assumption is invalid ; when calculating the variance of quadratic estimators all permutations of the observed trispectrum must be taken into account . one such permutation reflects the desired correlation making our estimator sensitive to @xmath89 , but the remaining two permutations induce additional variance proportional to the lensing - potential power spectrum @xmath106 . while subdominant , this variance will become increasing significant for future experiments as shown in the right - hand panel of fig . [ f : plrec ] . since the power spectrum @xmath106 is itself a measure of uncertainty in the deflection field , this additional variance in lensing reconstruction acts as a bias during power - spectrum estimation because there is no _ a priori _ way to distinguish it from the intrinsic variance of the underlying distribution . our calculation of the dependence of this variance on @xmath106 allows it in principle to be subtracted iteratively , which will prevent a systematic @xmath157 overestimate of @xmath106 at low @xmath108 . we close by considering several possible observational obstacles to the scheme for lensing reconstruction and power - spectrum estimation presented above . one hindrance is secondary contributions to the cmb such as the sz and isw effects . these effects increase the total temperature power spectrum appearing in the denominator of the optimum filter @xmath104 of eq . ( [ e : filt ] ) as would additional instrumental noise . they also correlate with the large - scale structure at low redshifts inducing further non - gaussian couplings and additional variance to lensing reconstruction . fortunately for our purposes the frequency dependence of the thermal sz effect differs from that of a blackbody . it can therefore be separated in principle from the lensed primordial cmb by an experiment with several frequency channels @xcite . the isw effect can not be removed in this manner , but is too small to significantly inhibit lensing reconstruction . polarization - dependent secondary effects are expected to appear at higher orders in the density contrast @xcite , and we therefore anticipate that they will not make a contribution at the levels considered here . a potentially more serious problem is that of galactic foregrounds , which though uncorrelated with the lensing signal may be substantial at certain frequencies . significant polarization has also been observed in some of these sources @xcite . we hope to understand and minimize the effects of galactic foregrounds in future work , and to pursue further refinements of lensing reconstruction . we thank wayne hu for useful discussions . this work was supported in part by nasa nag5 - 11985 and doe de - fg03 - 92-er40701 . kesden acknowledges the support of an nsf graduate fellowship and ac acknowledges support from the sherman fairchild foundation . here we provide the appropriate formulas for deriving the variance associated with polarization - based estimators of the deflection field @xmath89 . the cmb polarization can be decomposed into e and b - modes @xcite . these modes are mixed by weak lensing such that to linear order in @xmath68 , @xmath180 \phi({{\bf l}}-{{\bf l}}_1 ) \left [ ( { { \bf l}}-{{\bf l}}_1 ) \cdot { { \bf l}}_1 \right ] \ , , \nonumber \\ \tilde b({{\bf l } } ) & = & b({{\bf l } } ) - { \int \frac{d^2 { \bf l}_1}{(2\pi)^2 } } \left[e({{\bf l}}_1 ) \sin 2 ( \varphi_{{{\bf l}}_1 } -\varphi_{{\bf l } } ) + b({{\bf l}}_1 ) \cos 2 ( \varphi_{{{\bf l}}_1 } -\varphi_{{\bf l}})\right ] \phi({{\bf l}}-{{\bf l}}_1 ) \left [ ( { { \bf l}}-{{\bf l}}_1 ) \cdot { { \bf l}}_1 \right ] \ , .\end{aligned}\ ] ] we can exploit the sensitivity of the polarization modes to the lensing potential to construct lensing estimators from quadratic combinations of polarization modes . we generalize eq . ( [ e : ea ] ) to arbitrary combinations @xmath181 of @xmath45 , e , and b - modes as first derived in ref . @xcite , @xmath182 where @xmath183   \cos 2 ( \varphi_{{{\mathbf{l } } } } -\varphi_{{\mathbf{l^{\prime } } } } ) \ , , \nonumber \\ f_{eb}({{\mathbf{l } } } , { { \mathbf{l^{\prime } } } } ) & = & \left [ c_{l}^{ee } ( { { \mathbf{l}}}\cdot { { \mathbf{l } } } ) + c_{l^{\prime}}^{bb } ( { { \mathbf{l}}}\cdot { { \mathbf{l^{\prime } } } } ) \right ]   \sin 2 ( \varphi_{{{\mathbf{l } } } } -\varphi_{{\mathbf{l^{\prime } } } } ) = f_{be}({{\mathbf{l^{\prime } } } } , { { \mathbf{l } } } ) \ , , \nonumber \\ f_{bb}({{\mathbf{l } } } , { { \mathbf{l^{\prime } } } } ) & = & \left [ c_{l}^{bb } ( { { \mathbf{l}}}\cdot { { \mathbf{l } } } ) + c_{l^{\prime}}^{bb } ( { { \mathbf{l}}}\cdot { { \mathbf{l^{\prime } } } } ) \right ]   \cos 2 ( \varphi_{{{\mathbf{l } } } } -\varphi_{{\mathbf{l^{\prime } } } } ) \ , . \nonumber \\\end{aligned}\ ] ] in deriving these results we used parity considerations to demand @xmath184 . using these relations , we follow the approach of ref . @xcite to derive symmetric lensing estimators @xmath185 and @xmath186 { f_{xx'}}({{{\bf l}_1 } } , { { { \bf l}_2 } } ) \ , .\ ] ] we have explicitly symmetrized our estimators for @xmath187 to simplify the form of the optimal filter . the normalization bias of the estimators is removed by choosing @xmath188^{-1 } \ , , \ ] ] and @xmath189 { f_{xx'}}({{{\bf l}_1 } } , { { { \bf l}_2 } } ) \right]^{-1 } \ , .\ ] ] the minimum - variance filters @xmath190 for the various cases @xmath181 are given by @xmath191 using these optimal filters for the estimators defined in eqs . ( [ e : estxx ] ) and ( [ e : estxx ] ) , we can calculate the variances for these estimators in a fashion entirely analogoug to eq . ( [ e : estvar ] ) , @xmath192 @xmath193 \left [ x^{{\rm t}}({{{\bf l}_1 } } ' ) x'^{{\rm t}}({{{\bf l}_2 } } ' ) + x'^{{\rm t}}({{{\bf l}_1 } } ' ) x^{{\rm t}}({{{\bf l}_2 } } ' ) \right ] \rangle { f_{xx'}}({{{\bf l}_1 } } , { { { \bf l}_2 } } ) { f_{xx'}}({{{\bf l}_1 } } ' , { { { \bf l}_2 } } ' ) \nonumber \\ & & \quad \quad - ( 2 \pi)^2 \delta_{{\rm d}}({{\mathbf{l}}}- { { \mathbf{l^{\prime } } } } ) c_{l}^{dd } \ , . \nonumber \\\end{aligned}\ ] ] as for that of the temperature estimator , these variance will consist of zeroth - order terms in @xmath106 , @xmath194 and @xmath195 , and first - order terms , @xmath196 @xmath197 \nonumber \\ & & \quad \quad \quad + c_{|{{{\bf l}_1}}-{{{\bf l}_2}}'|}^{\phi\phi }   \bigl [ f_{xx'}(-{{{\bf l}_1 } } , { { { \bf l}_2 } } ' ) f_{x'x}(-{{{\bf l}_2 } } , { { { \bf l}_1 } } ' ) + f_{xx}(-{{{\bf l}_1 } } , { { { \bf l}_2 } } ' ) f_{x'x'}(-{{{\bf l}_2 } } , { { { \bf l}_1 } } ' ) \nonumber \\ & & \quad \quad \quad \quad + f_{x'x'}(-{{{\bf l}_1 } } , { { { \bf l}_2 } } ' ) f_{xx}(-{{{\bf l}_2 } } , { { { \bf l}_1 } } ' ) + f_{x'x}(-{{{\bf l}_1 } } , { { { \bf l}_2 } } ' ) f_{xx'}(-{{{\bf l}_2 } } , { { { \bf l}_1 } } ' ) \bigl ] \big\ } \ , , \nonumber \\\end{aligned}\ ] ] the filters given in eq .  ( [ e : filtxx ] ) are no longer optimal in the presence of this additional noise , but the difference between these filters and the optimal filters should be negligible provided that @xmath198 , @xmath199 . for the purposes of power - spectrum estimation , the terms @xmath200 and @xmath201 are not only an additional contribution to the variance , but are also a systematic bias if not subtracted iteratively following eq . ( [ e : naive3 ] ) . a final point to consider is that the six different estimators @xmath202 defined in this paper are not independent , as they are constructed from only three distinct maps . the covariance matrix for the six estimators will therefore not be diagonal , and this needs to be taken into account if the estimators are to be linearly combined to produce a single minimum - variance estimator . the off - diagonal elements of the covariance matrix can be evaluated in a straightforward manner involving pairs of double integrals similar to those of eqs . ( [ e : xxvar ] ) and ( [ e : xxvar ] ) .      p. j. e. peebles and j. t. yu , astrophys . j.  * 162 * , 815 ( 1970 ) ; r. a. sunyaev and ya . b. zeldovich , astrophys . space sci . * 7 * , 3 ( 1970 ) ; j. silk , astrophys . j.  * 151 * , 459 ( 1968 ) ; w. hu and s. dodelson , ann . astrophys .  * 40 * , 171 ( 2002 ) . g. jungman _ et al . _ , d  * 54 * , 1332 ( 1995 ) ; j. r. bond , g. efstathiou , and m. tegmark , mon . not . r. astron . soc .  * 291 * , l33 ( 1997 ) ; m. zaldarriaga and d. n. spergel , astrophys . j.  * 488 * , 1 ( 1997 ) ; d. j. eisenstein , w. hu , and m. tegmark , astrophys . j.  * 518 * , 2 ( 1999 ) . d. n. spergel and d. m. goldberg , phys . d  * 59 * , 103001 ( 1999 ) ; d. m. goldberg and d. n. spergel , phys . rev . d  * 59 * , 103002 ( 1999 ) ; a. cooray and w. hu , astrophys . j.  * 534 * , 533 ( 2000 ) ; m. zaldarriaga and u. seljak , phys . rev . d  * 59 * , 123507 ( 1999 ) ; h. v. peiris and d. n. spergel , astrophys . j.  * 540 * , 605 ( 2000 ) . g. giardino _ et al . _ , astrophys . * 387 * , 82 ( 2002 ) ; f. r. bouchet , s. prunet , and s. k. sethi , mon . not . r. astron . 302 * , 663 ( 1999 ) ; s. prunet , s. k. sethi , and f. r. bouchet , mon . not . r. astron . 314 * , 348 ( 2000 ) .
Abstract:

Answer: weak gravitational lensing by intervening large - scale structure induces a distinct signature in the cosmic microwave background ( cmb ) that can be used to reconstruct the weak - lensing displacement map . estimators for individual fourier modes of this map can be combined to produce an estimator for the lensing - potenial power spectrum . the naive estimator for this quantity will be biased upwards by the uncertainty associated with reconstructing individual modes ; we present an iterative scheme for removing this bias . the variance and covariance of the lensing - potenial power spectrum estimator are calculated and evaluated numerically in a @xmath0cdm universe for planck and future polarization - sensitive cmb experiments .