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 .


in recent years , systems where strong correlations between particles are induced by pair - wise scattering , have attracted much attention . in the canonical example , which is the parametric down - conversion ,
photon pairs are generated during the propagation of a laser beam through a non - linear medium .
the outcoming pairs of photons are entangled , and can serve as a probe of fundamental properties of quantum mechanics @xcite , such as the einstein - podolsky - rosen paradox , or violation of the bell inequalities @xcite . on the other hand , entanglement can be exploited in practical applications , such as teleportation @xcite or metrology beyond the shot - noise limit ( snl ) @xcite .    in this latter context , recent experiments with entangled states of atoms were a major breakthrough @xcite . in @xcite ,
two - body interactions were utilized to prepare non - classical squeezed states of atoms trapped in a double - well potential , which implies presence of many - body entanglement @xcite .
a similar idea was exploited to generate squeezing in the internal @xcite degrees of freedom . in @xcite , squeezing of a large spin of a collection of two - level atoms was achieved , using an intense laser field interacting with particles trapped in an optical cavity .
simultaneously , a substantial experimental effort was put in order to generate entangled pairs of atoms scattered out of a bose - einstein condensate ( bec ) . in @xcite , a collision of two becs lead to weak scattering of correlated atomic pairs onto a three - dimensional sphere of initially unoccupied modes .
although moderate number - squeezing between the opposite regions of the halo , and the related violation of the cauchy - schwarz inequality were experimentally demonstrated @xcite , entanglement was never directly observed .
alternatively , pair - production schemes were developed , where only few modes are strongly populated in a stimulated process , making the system somewhat easier to handle .
stimulated four - wave - mixing processes have been implemented using different spin states of atoms @xcite or bragg scattering @xcite . also , dynamic instabilities in moving optical lattices , populating modes with opposite quasi - momenta , have been used @xcite . in @xcite ,
a bec was transferred into the first excited state of a trapping potential and subsequent two - body collisions created a `` twin - beam '' system , where stronger - than - classical correlations could directly be observed .
analogous schemes have been implemented in internal atomic states , building upon spin - changing collisions @xcite .
furthermore , in @xcite it was shown that particles scattered in this process into a pair of @xmath0 zeeman sub - levels are usefully entangled from the metrological point of view .    in this work
we develop a theoretical model for the generic type of experiments , where particles scatter in pairs into two well - separated regions .
if these regions are separated in the momentum space , they could also be set apart by a sufficiently long expansion of the cloud . on the other hand , in cases
when the regions are defined as two zeeman sub - levels , they can be separated in space using a stern - gerlach scheme .
our model applies to any such possible configuration .
therefore , the general conclusions of this study , concerning the correlations , number - squeezing and entanglement , are valid for recent experiments @xcite in the regime , where the depletion of the source bec is low .
we derive the bogoliubov equations governing the dynamics of pair formation , and applying the bloch - messiah reduction @xcite , we write the state in terms of pairs of independently squeezed modes .
we calculate the density and the number of scattered atoms , and the two body correlation between them .
we demonstrate the presence of ideal number - squeezing between the opposite regions , violation of the cauchy - schwarz inequality and , using the fisher information criterion known from quantum metrology , we show that the atoms from the twin - beam system are entangled .
we also provide a simple yet useful lower bound for the fisher information in terms of the average number of scattered atoms , and the number of modes they occupy .
finally , we apply the above formalism to the twin - beam experiment of @xcite .
this paper is organized as follows . in sec .  [ universal ]
we discuss the general properties of the solutions of the bogoliubov equations .
these observation allow to easily calculate the density and the second order coherence of the system in sec .  [ density ] and the fluctuations of the population imbalance between the opposite regions in sec .  [ squeezing ] . in sec .
[ sec_csi ] we take the first step towards the demonstration of particle entanglement present in the system , by showing that the second order correlation function violates the cauchy - schwarz inequality . in sec .
[ entanglement ] we demonstrate that the scattered atoms are usefully entangled from the metrological point of view . in sec .
[ outline ] we briefly describe the experimental setup of @xcite and in [ theory ] derive the corresponding effective bogoliubov equations . finally , in sec .
[ results ] we review the most relevant properties of the twin - beam system by showing the results of the numerical simulation .
we conclude in sec .
[ conclusions ] .
we first present the general properties of the solution of the bogoliubov equation , in cases where particles scatter pair - wise into well - separated regions .
our theoretical description of the pair - production process starts with a many - body hamiltonian with contact two - body interactions    [ ham ] @xmath1    here @xmath2 is an external trapping potential and @xmath3 is the strength of the two - body interactions , @xmath4 is the scattering length , @xmath5 is the atomic mass and the field operator @xmath6 satisfies the bosonic commutation relations . to derive the bogoliubov equation
, we first find the c - number ( mean field ) wave function of the bec using the gross - pitaevskii equation ( gpe ) @xmath7 we then write the field operator as a sum of the c - number part and the bogoliubov correction , @xmath8 and insert this expression into ( [ ham ] ) . by keeping only the terms up to quadratic in @xmath9
we obtain the bogoliubov hamiltonian @xmath10 the resulting bogoliubov equation of motion is linear @xmath11 usually , a numerical solution of this equation is found in a following way .
the field operator is expanded in a basis of wave - functions @xmath12 which match the geometry of the scattering problem @xmath13 this expression is inserted into eq .
( [ bog ] ) , the resulting equation is multiplied by @xmath14 and the outcome is integrated by sides over the whole space . in effect , what we obtain is an equation of motion , which , through the matrices @xmath15 and @xmath16 , couples the evolution of the @xmath17-th operator @xmath18 , with ( in general ) all others operators @xmath19 this equation is linear  a consequence of the linearity of the bogoliubov equation ( [ bog ] )  so the general solution of ( [ eq_mot ] ) reads @xmath20 where the matrices @xmath21 and @xmath22 satisfy @xmath23 and @xmath24 . later , we will apply this method to solve the bogoliubov dynamics of the twin - beam production .
however , we will show in the following , that in cases where the detailed form of the hamiltonian ( [ bog_ham ] ) drives the scattering of atomic pairs into opposite regions ( as indeed happens in twin - beam experiments ) , the basic properties of the system can be deduced analytically if an appropriate set of mode functions @xmath25 is chosen .
let us denote the two separate regions into which the particles are scattered by l ( left ) and r ( right ) .
particles populate l and r in a process of elastic scattering , so the regions are usually separated in momentum space . from this point of view , it is convenient to switch to the space of wave - vectors @xmath26 and decompose the field operator as follows @xmath27 the operators @xmath28 annihilate a particle in a mode characterized by the time - dependent wave function @xmath29 , which is localized in the right / left region in momentum space .
we underline , that this kind of separation is also present in position space after expansion of the cloud , or after application of a stern - gerlach pulse in internal - state experiments , respectively .
moreover , the vector @xmath26 might denote the quasi - momentum , if the scattering takes place in an optical lattice .
formally , the only difference between the formulation ( [ op_nondiag ] ) and ( [ op_diag ] ) is the splitting of the field operator into the r and l modes .
however , for a linear equation of motion such as ( [ bog ] ) , there exists a unique basis of mode functions for which the evolution equations of the mode pairs decouple from each other :    [ pair_evolution ] @xmath30    where @xmath31 .
this form of the bogoliubov equation has a simple physical interpretation : atoms scatter pair - wise into opposite regions , and the total field operator ( [ op_diag ] ) is a sum of independent mode pairs , which are squeezed in their relative population fluctuations , as will be explained in detail below .
although the diagonal form ( [ pair_evolution ] ) is much clearer than ( [ nondiag ] ) , it is not obvious at the moment how this particular basis ( [ op_diag ] ) can be found .
this is done in two steps , applying the procedure of the bloch - messiah reduction @xcite .
first , using equations ( [ op_diag ] ) and ( [ pair_evolution ] ) , we evaluate the one - body density matrix ( first - order correlation function ) and obtain @xmath32 where @xmath33 .
note that @xmath34 is the average occupation of the @xmath35-th eigen - mode , and that a pair of modes @xmath36 and @xmath37 is degenerate ( has the same eigen - value @xmath34 ) due to the assumed symmetry between the left and the right region . since we are using the heisenberg picture , the average value in eq .
( [ den_mat ] ) and all equations that follow are calculated in the initial vacuum state of scattered atoms .    in any practical approach ,
if the basis ( [ op_diag ] ) is not known a priori , this step first requires a numerical evaluation of the density matrix ( [ den_mat ] ) in any convenient basis ( [ op_nondiag ] ) , and subsequent diagonalization .
once this is done , then according to eq .
( [ den_mat ] ) , the basis functions @xmath29 are the momentary eigen - functions of the one - body density matrix ( natural orbitals ) .
however , a second step is necessary to fully determine the functions @xmath29 , because the density matrix  contrary to the field operator ( [ op_diag ] )  is insensitive to global phases of the mode functions . to retrieve this additional information
, we calculate the anomalous density @xmath38 multiply it by sides with the eigen - functions of the density matrix and integrate over space . as a result , we retrive the information about the phases and obtain the full form of the mode functions @xmath29 of the diagonal basis .    to summarize , we have outlined the structure of the solution of the bogoliubov equation for cases where atoms are scattered into two opposite regions .
we will now show that the extra step , which is the transition from the `` numerical approach '' ( [ op_nondiag ] ) to the diagonal basis ( [ op_diag ] ) , allows to easily determine the basic properties of the system of scattered atoms , like its density or higher correlation functions .
the simplest observable characterizing the pair - production process is the density @xmath39 which is , consistently with our derivation , localized in the two opposite regions . by integrating the above function over space ,
we obtain the information about the expected number of scattered atoms as a function of time @xmath40 additional information about the system is carried by the correlations between the scattered particles .
the probability of simultaneous detection of two atoms at momenta @xmath41 and @xmath42 can be obtained from the normalized second - order correlation function @xmath43 according to the wick s theorem , this function can be written in terms of the one - body density matrix ( [ den_mat ] ) and the anomalous density ( [ anom ] ) as follows @xmath44 the transition from eq .
( [ g2 ] ) to ( [ g2_w ] ) might seem an unnecessary complication , however we will argue that it allows for a simple and intuitive interpretation of the second - order correlation function . according to eq .
( [ den_mat ] ) , the density matrix is non - vanishing only when @xmath41 and @xmath42 are both either in the right or left region , so @xmath45 governs the hanbury - brown and twiss ( hbt ) type of local correlations . on the other hand , as can be seen from eq .
( [ anom ] ) , the anomalous density is non - zero only when @xmath41 and @xmath42 are in the opposite regions , so it describes the cross - correlations between the two members of the scattered pair . clearly , this simple interpretation of the second order correlation function as a sum of local- and opposite - momentum correlations would have been much more difficult if we had not applied the diagonalization procedure and the wick s theorem .
another property characterizing the scattering process are the fluctuations of the population imbalance between the two regions .
if these fluctuations are suppressed below the properly defined shot - noise level , the system is _ number squeezed _ , which proves that the atoms scatter in pairs rather then independently to the left and to the right region . a quantitative description of the number squeezing involves the left and right atom number operators defined as the integrals of the density operators over the corresponding volumes , i.e. @xmath46
the population imbalance operator is then simply defined as @xmath47 and using eq .
( [ op_diag ] ) we obtain @xmath48 the number squeezing factor is defined as @xmath49 where @xmath50 is the variance of the population imbalance operator . if the fluctuations between the two regions are suppressed below the shot - noise level defined as @xmath51 , the system is called `` number - squeezed '' . in our case , since @xmath52 does not depend on time and the initial state is a vacuum , we obtain that @xmath53 . therefore , the two - region bogoliubov system is perfectly number - squeezed , as anticipated in the previous section .
the ideal number - squeezing is a result of clear separation of the two scattering regions . in such case ,
it is natural to define the local atom - number operators ( [ at_num ] ) and the population imbalance operator ( [ pop_imb ] ) .
it is important to note that not all systems , where particles are scattered in pairs are perfectly number squeezed .
for instance , when two bose - einstein condensates collide , they produce a halo of atoms due to two - body elastic scattering into the initially unoccupied modes @xcite . in this system
however , there is no simple way to define two separate regions .
one can instead measure the number of atoms in two bins lying on the opposite sides of the halo .
moderate number - squeezing of the atom number difference between these bins has been observed experimentally @xcite , but it is impossible to reach the limit @xmath54 @xcite .
in contrast , the twin matter wave configurations @xcite , are ideal sources of correlated atomic pairs occupying two well - defined areas .      apart from the number squeezing
, the twin - region system can be characterized by another expression , which is called the cauchy - schwarz inequality .
it relates the strength of the local and opposite correlations to witness the pair - scattering process .
following @xcite , we define averaged second - order correlations as @xmath55 where @xmath56 . in the symmetric case ,
the cauchy - schwarz inequality @xmath57 can now be re - written as @xmath58 using expressions ( [ den_mat ] ) and ( [ anom ] ) we obtain    @xmath59    thus the cauchy - schwarz inequality reads @xmath60 which is true only for all @xmath61 .
as soon as particles start to scatter into the two regions , the cauchy - schwarz inequality is clearly violated . to quantify the degree of violation ,
a coefficient @xmath62 was introduced in @xcite , which reads @xmath63 when @xmath64 the system is in the `` classical '' regime , while @xmath65 signify correlations which are stronger than allowed by the classical physics . in our case
this coefficient reads @xmath66 clearly , always @xmath65 , because it is a sum of unity and a non - negative part . for high mode populations @xmath34 , the second term , which is inversely proportional to the number of scattered particles tends to zero , restoring the classical limit
nevertheless , as demonstrated with photons in @xcite , the confidence by which the cauchy - schwarz inequality can be violated in the presence of classical noise still _ increases _ with more strongly populated modes .
however , it is the fisher information , which is the quantity highly sensitive to particle entanglement in the high - gain regime , as we show in the following section .
this measure quantifies the potential for sub - shot - noise interferometry , and increases with rising mode population , in spite of the decreasing `` granularity '' of the matter wave @xcite that leads to all second - order correlation functions approaching equal values .
we now show that atoms occupying the two regions are entangled , and could be used as an input of a quantum interferometer operating below the shot - noise level .
we first recall how the precision of the phase estimation is related to the entanglement of input states using as an example the standard two - mode mach - zehnder interferometer ( mzi ) .
then , we extend these concepts to the case , where the interferometer operates between two regions , each having a multi - mode structure determined by the bogoliubov equations .
when speaking about two - mode interferometers , it is convenient to introduce a set of three operators    [ ang ] @xmath67    which obey the same commutation relations as the angular momentum operators .
the mzi , which is an interferometric device , where the imprint of the phase @xmath68 onto the input state is preceded and followed by a pair of symmetric beam - splitters , can be represented by a unitary evolution operator @xmath69 .
if the phase is estimated in a series of @xmath70 measurements performed on the output state , the precision of the phase estimation is limited by the cramer - rao lower bound ( crlb ) @xcite , @xmath71 here , @xmath72 is the quantum fisher information ( qfi ) , which is related to the unitary transformation @xmath73 . for pure states transformed by the mzi it is equal to @xmath74 , where the variance is calculated in the input state of the interferometer @xcite .
the crlb states , that if @xmath68 is determined using any possible type of measurement and estimator , then the precision @xmath75 is bounded as in eq .
( [ crlb ] ) .    apart from providing a lower bound for the error of the phase estimation
, the @xmath76 is an entanglement measure .
namely , when the input state has an average number of @xmath77 particles , then if @xmath78 , the state is particle - entangled @xcite .
we now show , that a natural extension of the two - mode picture allows to employ the concept of the qfi as an entanglement measure also in our multi - mode system of interest . to this end , we introduce the following analog of the two - mode angular momentum operators ( [ ang ] ) ,    [ ang_twin ] @xmath79    where we dropped the explicit time - dependence of the @xmath80 to simplify the notation . also , for simplicity , we choose the well - separated regions r and l to be localized symmetrically on the opposite sites of @xmath81 .
the construction of these operators , which satisfy the same commutation relations as ( [ ang ] ) , is based on the analogy between the two - mode systems and the twin - beam configuration . in the former case ,
the operators connect the right and left modes , while in the latter the left and right sub - spaces .
such a definition ( [ ang_twin ] ) is meaningful only in situations , where the system consists of two well - separated sub - systems . using the decomposition of the field operator into the set of independent modes , equations ( [ op_diag ] ) and ( [ pair_evolution ] ) , the above integrals yield , that each angular momentum operator is a sum of operators acting on each mode independently , that is    [ ang2 ] @xmath82    these expressions show again that it is natural to describe the two - region system using the diagonal basis ( [ op_diag ] ) . in this language ,
the angular momentum operators are simply a sum of operators acting on each pair of modes independently , which vastly simplifies the further analysis .    to establish a direct relation between the two - mode and two - region case
, we now assume that the system is transformed in the multi - mode analog of the mach - zehnder interferometer . as outlined above , to demonstrate the presence of useful entanglement between the atoms in the left and in the right , it is necessary to calculate the qfi . using eq .
( [ ang2 ] ) we obtain that @xmath83 since the construction of the basis ( [ op_diag ] ) explicitly assumes that each mode is independent from all others , the second term in the last equality is @xmath84 , because the symmetry between the @xmath85 and @xmath86 regions implies that @xmath87 for all @xmath35 .
therefore we obtain that the qfi is equal to @xmath88 where the last equality comes directly from the substitution of ( [ pair_evolution ] ) into the definition of the @xmath89 operator .
also , we used @xmath90 , according to eq .
( [ nscat ] ) .
clearly @xmath78 , so the system is entangled .
moreover , one can refer the qfi to the ultimate bound for the precision of the parameter estimation , which is the heisenberg limit . for a system with fluctuating number of particles ,
this upper bound is equal to @xmath91 . using ( [ pair_evolution ] ) again
, we obtain , that @xcite @xmath92 for a large number of scattered particles , when @xmath93 , we obtain @xmath94 . the value of the qfi , which is only one - half smaller than the heisenberg limit is a clear indication of very strong entanglement present in the system in the high - gain regime . at intermediate times , @xmath95 due to mode competition , which has a negative impact on the entanglement as witnessed by the qfi @xcite . to picture this ,
consider a `` frustrated case '' , where all atoms scatter uniformly into @xmath96 pairs of modes , so that all @xmath97 are equal . in this case , the number of scattered atoms is simply @xmath98 , and the qfi is @xmath99 .
the qfi normalized to the snl is @xmath100 when , on average , there is less than a particle per a set of modes , i.e. @xmath101 , the qfi surpasses the snl only by a factor of 2 , a natural reminiscence of atoms being scattered in pairs .
equation ( [ snl ] ) is a simple yet intuitive estimation of the lower bound of useful entanglement in terms of the number of scattered atoms and occupied modes .
we now apply the above formalism to the twin - beam system of @xcite .
first , we describe the physical mechanism which leads to the creation of the two correlated beams .
as shown below , some basic information about the dynamics of the pair production allow to construct a simple one - dimensional bogoliubov model , which can be easily solved numerically .
the experimental sequence applied in @xcite to produce correlated atom pairs was following .
first , an almost pure bose - einstein condensate ( bec ) of @xmath102 @xmath103rb atoms with scattering length equal to @xmath104 nm was created at temperature @xmath105nk .
the cloud was trapped in an approximately harmonic potential @xmath106 where atomic mass is equal to @xmath107 , and the frequency @xmath108 is much smaller than @xmath109 and @xmath110 , so the bec is strongly elongated along the @xmath111-axis .
after the bec was created , the trapping potential was shaken in a controlled way , so the atoms were transferred to the first excited state along the @xmath112-direction . in order to achieve the maximal transfer efficiency ,
the shaking was optimized using quantum optimal control theory @xcite . afterwards , binary collisions transfered particle pairs to the ground state of the potential , and the excess energy @xmath113 was converted into back - to - back movement of the two atoms along @xmath111 .
momentum conservation ensured , that their wave vectors had equal lengths @xmath114 and point in opposite directions .
small corrections to the value of @xmath115 may arise from an effective mean - field potential , as will be discussed below .      neglecting thermal phase fluctuations along the elongated direction @xmath111 , which is valid at very low temperatures only  @xcite , the condensate wave function acting as a source for the pair - production
can be found by solving the stationary gpe @xmath116 where @xmath117 is the chemical potential .
this function can be evaluated numerically , by referring to the description of the experiment from the previous section , and noting that after the shaking of the trap , the bec is in the first excited state @xmath118 along the @xmath112 axis and in the ground state @xmath119 along @xmath120 . however ,
this can be approximated by an analytical expression , as argued below .
first note , that since the characteristic energies @xmath121 and @xmath122 are large , and the number of atoms in the bec is small , the non - linear term can be safely neglected in evaluation of the eigenstates along @xmath112 and @xmath120 . as a result , assuming that the total wave - function @xmath123 separates in three directions ( which has been confirmed numerically ) , we obtain @xmath124 where the functions @xmath125 and @xmath126 are the eigen - states of the one - dimensional harmonic potential in @xmath112 and @xmath120 correspondingly .
the function @xmath127 is found by inserting the above expression into eq .
( [ stat_gpe ] ) and integrating out the orthogonal directions . as a result
, we obtain an effective equation @xmath128 where zero - point energy equals @xmath129 and the non - linearity reads @xmath130\times \left[\int\!\ !
\mathrm{d}z\ , \left(\psi^{\rm ( ho)}_{n_z=0}(z)\right)^4\right]\nonumber\\    & = & \frac{3}{8\pi}\frac g{a_{\mathrm{ho},y}\ , a_{\mathrm{ho},z}}.\end{aligned}\ ] ] here @xmath131 are the harmonic oscillator lengths for @xmath132 .
since the trap is shallow in the @xmath111 direction , the solution of the stationary gpe ( [ gpex ] ) can be well approximated by the thomas - fermi ( tf ) formula @xcite @xmath133 where the effective chemical potential is @xmath134 leading to a tf radius of @xmath135 .
within the approximation of neglecting thermal phase fluctuations , we have fully determined the wave - function of the bec , which we insert into the bogoliubov hamiltonian ( [ bog_ham ] ) .
next , we expand the field operator @xmath136 in an orthonormal basis . along the @xmath112 and @xmath120 directions ,
it is natural to use the eigen - states of the harmonic oscillator as the basis functions , since it matches the geometry of the source bec . along the @xmath111 direction
, we use a plane - wave basis , and get @xmath137 since the atom pairs are emitted into the ground state along @xmath112 only ( which is ensured by the anisotropy and anharmonicity of the potential ) , the sum over the eigen - states can be safely truncated at @xmath138 and @xmath119 .
this reduces the dynamics of the pair - production to one - dimensional problem along the @xmath111 axis , with the orthogonal directions frozen out , i.e. @xmath139 where @xmath86 is the quantization volume .
we insert this field operator into eq .
( [ bog_ham ] ) , evaluate the spatial integrals and upon the change of variables @xmath140 obtain @xmath141 where @xmath142 .
we solve the resulting bogoliubov equation numerically , using a fourth order runge - kutta method , and find the matrices @xmath21 and @xmath22 as defined in eq .
( [ nondiag ] ) .    using the above hamiltonian
, one can also analytically determine @xmath115 , i.e. the position of the central peak . to this end
, we employ a two - mode approximation by replacing the function @xmath143 with a dirac delta , and obtain the bogoliubov equation @xmath144 where @xmath115 is shifted with respect to the harmonic excitation energy due to the mean - field repulsion and reads @xmath145 this result is in good agreement with the experimentally measured position of the peak density , i.e. @xmath146 .      in this section , we display the most important characteristics of the twin - beam system , starting from the solution of the eigen - problem of the density matrix ( [ den_mat ] ) . in fig .  [ fig_vals ] we plot the first four eigen - values of the density matrix , as a function of time .
the inset shows the total number of scattered atoms @xmath77 normalized to the occupation of the bec , as a function of time .
the bogoliubov approximation is valid for as long as @xmath147 , so we interrupt the simulation at @xmath148ms , when @xmath149 . for longer times ,
when the depletion of the bec can not be neglected , a atom - number conserving method , such as the one introduced in @xcite must be used .
ms , when @xmath150 . ]    in fig .
[ fig_vecs ] we plot the first four eigen - vectors of @xmath151 localized in the right half - space , i.e. @xmath152 with @xmath153 , calculated at an early time @xmath154ms and at @xmath148ms . due to the time - energy uncertainty relation ,
the eigen - modes localize around @xmath155 at later times .     with @xmath153 .
the solid black lines are results of diagonalization of the density - matrix at @xmath154ms while the dashed red lines at @xmath148ms .
the figure shows how due to the time - energy uncertainty relation , the eigen - vectors narrow in the course of time around the central wave - vector @xmath115 . ]
this can be seen even more clearly , by plotting the density @xmath156 at these two instants , as shown in fig .
[ fig_den ] ( dashed lines ) . at @xmath154ms , two broad beams start to form on top of the uniform density .
later , at @xmath148ms , strongly localized peaks clearly dominate over the flat background . on top of these curves
, we plot the normalized second - order correlation function as defined in eq .
( [ g2 ] ) , with one of the arguments set equal to the resonant wave - vector @xmath115 , i.e. @xmath157 . at @xmath154ms , the cross - correlation , which is governed by the anomalous density , is very large , i.e. @xmath158 .
this is a characteristic property of the bogoliubov system in the low - occupation regime @xcite , and indicates strong violation of the cauchy - schwarz inequality  ( [ csi_one_body ] ) . also , for this early time , the width of both @xmath159 peaks are much more narrow than the beam size .
this is consistent with the results shown in fig .
[ fig_vals ] , where at early times many eigen - mode pairs of the density matrix are almost equally occupied . at later times , when a single pair of modes start to become dominant , the width of the peak in @xmath159 and the system size approach each other . while this corresponds to beams that are single - mode with respect to their local one - body properties , the local averaged correlation function as introduced in eq .
( [ eq : summed_g2 ] ) reaches the limit of @xmath160 , exceeding the number fluctuations of a coherent state by a factor of two .     for fixed @xmath161 ( solid lines , left @xmath112-axis ) , and density profiles @xmath162 ( dashed lines , right axis ) in momentum space .
the results are calculated at @xmath154ms ( a ) and @xmath148ms ( b ) . at early times , many momentum modes are occupied and the width of @xmath159 is much smaller than the beam size .
later , two distinct peaks emerge , which are almost single - mode . ]    .
contributions of @xmath151 and @xmath96 in eq .
( [ g2_w ] ) are shown in blue and red hues , respectively .
the dashed line indicates the position of the cut shown in fig .
[ fig_den](b ) .
arrows indicate the axes of fig .
[ fig_avgd ] . ]    , as obtained in experiments .
solid line : momentum - space peak near @xmath163 , shown along the difference coordinate @xmath164 as indicated by the arrow in fig .
[ fig_k1k2 ] .
dotted line : peak near @xmath165 , along the sum coordinate @xmath166 . dashed and
dash - dotted lines : respective functions , taking into account the finite expansion time in time - of - flight momentum measurements .
the grey shaded area is proportional to the normalization @xmath167 . ]
the typical tool to capture local second - order correlations in experiments are collinearly integrated functions of the type @xmath168 where the integrals run over an appropriately chosen momentum region  @xcite . for symmetric ,
non - local correlations , back - to - back integration of the type @xmath169 is used .
the corresponding normalized functions for our system @xmath170 are shown in fig .
[ fig_avgd ] at @xmath171 . for both local ( solid ) and
non - local ( dotted ) functions , correlations peaks , which do not span the entire populated range ( grey area ) and follow a gaussian shape , are clearly present .    in the next step
we take towards future comparison with experiments , we present the results not in momentum space , but rather using real - space data calculated after some finite time @xmath172 of ballistic expansion . only in the limit of @xmath173 ( far field ) ,
the real - space data is equivalent to the initial momentum space distribution ( if the expanding clouds are sufficiently dilute , so that the mean - field repulsion can be safely neglected ) . in  @xcite ,
the expansion time was @xmath174 , which was sufficient to resolve the twin - beam peaks . nevertheless , the system was not fully in the far - field regime yet , which may have some impact on the correlation functions . as shown in fig .
[ fig_avgd ] , the finite expansion time affects the back - to - back peak at @xmath175 much more strongly than the collinear hbt peak , leading to smearing of the measured @xmath176 ( dash - dotted line ) over the entire size of the twin - beam packets .
this observation is consistent with some previous results @xcite .
intuitively , the broadening effect is related to the random position of scattering events along @xmath111 within the size of the initial cloud , which is non - vanishing with respect to the expanded size of the twin - beam peaks . on the other hand ,
the local correlation function @xmath177 ( dashed line ) remains largely unaffected .
note that although at every instant of the evolution , the field operator @xmath80 can be written as a sum of independently squeezed modes , at very early times the division between the right and left modes is unjustified , because the two peaks are not yet fully separated .
however , at @xmath154ms when the density distribution is broad , the number of scattered atoms is @xmath178 .
therefore , the system at such early time is hardly accessible experimentally so the quantum state of much less then a single particle is not of interest .
as soon as the two peaks are well - formed , at @xmath179ms , with @xmath180 scattered atoms , all the general considerations from sec .
[ general ] apply .    finally , in fig .
[ fig_qfi ] we plot the qfi from eq .
( [ qfi ] ) as a function of time and normalized to the heisenberg limit , i.e. @xmath181 . instead of interrupting the simulation at @xmath182ms , where the scattered fraction of atoms becomes non - negligible and particle number
conservation is strongly violated , we extend the calculation up to @xmath183ms , when the number of scattered atoms significantly exceeds 15% of @xmath184 .
this is done solely to illustrate that , once the population of one of the modes dominates , @xmath185 , as argued in sec .
[ entanglement ] .
note that the dominance of a single mode pair at long times is also predicted by the number - conserving theory @xcite , justifying this proceeding .
indeed , in the inset , we show the number of pairs of right / left modes which have an occupation bigger or equal to 10% of the largest mode .
this approximately tells , how many modes are significantly occupied in the system . at early times , there are over 100 pairs of modes . at @xmath182ms
there are still 5 significantly occupied pairs , and only around 4.2ms a single pair of modes starts to dominate . at the same time the qfi approaches its upper bound .    .
the horizontal grey dashed line denotes the best possible value for the bogoliubov system , which is achieved in a regime , where only a single pair of left / right modes is relevant .
the vertical grey dashed line denotes the time @xmath148ms , when the bogoliubov simulation should be interrupted . in the inset ,
we show the number of pairs of left / right modes that have at least 10% of occupation of the largest modes .
we see that only around 4.2ms , the two - mode approximation is valid , as denoted by the horizontal dashed line . ]
we have developed a simple bogoliubov model describing twin - atom beam experiments similar to refs .  @xcite . due to the elongated geometry of the trapping potential ,
the dynamics is one - dimensional . as a consequence
, the final step of our method can be easily solved numerically without the need for any stochastic method .
furthermore , basic information about the scattered particles can directly be drawn from the general properties of the solution of the bogoliubov equations . in this way
, we can quantitatively characterize the mode structure and correlation functions of the scattered atoms .
also , quite generally , we can show that the population imbalance between the two beams is ideally squeezed and that the system is strongly entangled .
these general observations can be applied to most recent experiments , where the atomic pairs scatter into two well - separated regions . finally , using the notion of the quantum fisher information
, we have derived a simple lower bound for the useful entanglement of the system .
this expression employs only the average number of scattered particles and the number of occupied modes .
having understood the fundamental properties of few - mode twin beams , further steps can be made to take into account more specific issues of experimental implementations .
a general feature of strongly elongated bose gases at realistic temperatures , such as the source cloud in @xcite , is the quasi - condensation  @xcite , where the coherence along the @xmath111-axis is limited due to thermal phase fluctuations .
although these fluctuations do not alter the general considerations of sec .
[ general ] , they affect the emission dynamics @xcite , and also might have influence on the spatial properties of both the density and the correlation functions @xcite . in future
, our method could be applied to the analysis of some more complicated schemes , basing on the twin - beam setup , such as the rarity - tapster - type experiments @xcite .
also , according to our results , the two - region state could be used as an input to the mach - zehnder - like interferometer , similarly to @xcite .
finally , note that the bogoliubov approximation neglects the secondary collisions between the scattered particles and the atoms from the source cloud . when a scattered atom propagates through the bec , the number of secondary collisions is bounded from above by @xmath186 , where @xmath187 is the peak density of the bec . by plugging in the experimental numbers , we obtain @xmath188 , which well justifies the use of the bogoliubov approximation
j. ch . acknowledges the foundation for polish science international team programme co - financed by the eu european regional development fund .
t.w . and p. sz .
acknowledge the foundation for polish science international ph.d .
projects programme co - financed by the eu european regional development fund . r.b . acknowledges support from the austrian science fund ( fwf ) projects cap ( 1607-n16 ) , and atom chip ( z118-n16 ) , and the fwf doctoral programme coqus ( w1210 ) .
this research was supported by the national science center grants no .
dec-2011/03/d / st2/00200 and n202 167840 .
paul g. kwiat , klaus mattle , harald weinfurter , anton zeilinger , alexander v. sergienko , and yanhua shih , phys .
lett . * 75 * , 4337 ( 1995 ) r. horodecki , p. horodecki , m. horodecki , and k. horodecki , rev .
phys . * 81 * , 865 ( 2009 )      charles h. bennett , gilles brassard , claude crpeau , richard jozsa , asher peres , and william k. wootters , phys .
70 * , 1895 ( 1993 ) charles h. bennett , gilles brassard , sandu popescu , benjamin schumacher , john a. smolin , and william k. wootters , phys .
lett . * 76 * , 722 ( 1996 )    a.  perrin , h.  chang , v.  krachmalnicoff , m.  schellekens , d.  boiron , a.  aspect , c.i .
westbrook , phys .
lett .  * 99 * , 150405 ( 2007 ) .
v.  krachmalnicoff , j .- c .
jaskula , m.  bonneau , v.  leung , g.  b.  partridge , d.  boiron , c.  i.  westbrook , p.  deuar , p.  zi , m.  trippenbach , k.  v.  kheruntsyan , phys .
lett .  * 104 * , 150402 ( 2010 ) .
jaskula , m.  bonneau , g.  b.  partridge , v.  krachmalnicoff , p.  deuar , k.  v.  kheruntsyan , a.  aspect , d.  boiron , c.  i.  westbrook , phys .
105 * , 190402 ( 2010 ) k.v .
kheruntsyan , j .- c .
jaskula , p.  deuar , m.  bonneau , g.b .
partridge , j.  ruaudel , r.  lopes , d.  boiron , c.i .
westbrook , phys .
lett .  * 108 * , 260401 ( 2012 ) .
r. g. dall , l. j. byron , a. g. truscott , g. r. dennis , m. t. johnsson , and j. j. hope , phys .
a * 79 * , 011601 ( 2009 ) d. pertot , b. gadway , and d. schneble , phys .
* 104 * , 200402 ( 2010 ) wu rugway , s. s. hodgman , r. g. dall , m. t. johnsson , and a. g. truscott , phys .
lett . * 107 * , 075301 ( 2011 )        r. bcker , j. grond , s. manz , t. berrada , t. betz , c. koller , u. hohenester , t. schumm , a. perrin and j. schmiedmayer , nat . phys . * 7 * , 608 ( 2011 ) b. lcke , m. scherer , j. kruse , l. pezz , f. deuretzbacher , p. hyllus , o. topic , j. peise , w. ertmer , j. arlt , l. santos , a. smerzi and c. klempt , science * 11 * , 773 ( 2011 )      j. m. vogels , k. xu , and w. ketterle , phys .
lett . * 89 * , 020401 ( 2002 ) a. perrin , h. chang , v. krachmalnicoff , m. schellekens , d. boiron , a. aspect , and c. i. westbrook , phys .
lett . * 99 * , 150405 ( 2007 )                  r. bcker , t. berrada , s. van frank , j .- f .
schaff , t. schumm , j. schmiedmayer , g. jger , j. grond and u. hohenester , j. phys .
b * 46 * , 104012 ( 2013 ) d.s . petrov , g.v .
shlyapnikov and j.t.m .
walraven , phys .
lett . * 87 * , 050404 ( 2001 ) f. dalfovo , s. giorgini , l. p. pitaevskii and s. stringari , rev .
71 * , 463 ( 1999 )