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 idea that quantum mechanical systems can efficiently simulate physical systems @xcite is at the heart of the theory of quantum information and computation . deutsch @xcite constructed the first example of an algorithm which would require two queries to solve on a classical computer but that can be solved with only one quantum query . subsequently , deutsch and jozsa @xcite , bernstein and vazirani @xcite and simon @xcite demonstrated the striking difference between the classical and quantum query complexity . this culminated with grover s search algorithm @xcite achieving a quadratic speedup over any classical algorithm for the unstructured search problem along with shor s factoring algorithm @xcite that could factor integers efficiently and also calculate discrete logarithms . grover s algorithm gives an optimal solution to the unstructured search problem , the problem of deciding whether a black - box boolean function has any input that evaluates to 1 . it provides a quadratic temporal speedup over the best classical search algorithms , even when they both require the same spatial resources to perform the search task . zalka @xcite proved the optimality of the grover search algorithm . grover s algorithm was subsequently generalized to the framework of amplitude amplification and to counting the number of solutions by brassard _ grover s search algorithm has since been applied to a wide variety of similar problems @xcite@xcite@xcite . the usual paradigm of computation ( quantum or classical ) is defined in a discrete setting . however , adiabatic quantum computation @xcite provides a continuous time model for quantum computing by using the quantum adiabatic theorem . here , one begins by finding an appropriate hamiltonian whose ground state describes the solution to the problem of interest . next , a system with a simple hamiltonian is prepared and initialized to the ground state . finally , the simple hamiltonian is adiabatically evolved to the desired hamiltonian . by the quantum adiabatic theorem , the system remains in the ground state at all times during this adiabatic evolution , depending on the spectral properties of the hamiltonian @xcite@xcite . at the end , the final state of the system describes the solution to the problem . the time complexity for an adiabatic algorithm is simply defined as the time taken to complete the adiabatic evolution , which is known to depend on the spectral gap of the hamiltonian @xcite . the adiabatic model of quantum computing is universal and equivalent to standard quantum computation @xcite . farhi and gutmann were the first to construct an analog algorithm for solving the grover search problem @xcite referred to hereafter as the _ analog analogue of grover search algorithm _ or simply the _ analog grover search algorithm_.    it is understood that the speedup in quantum algorithms comes from quantum mechanical features like quantum coherence and quantum correlations like quantum entanglement , that exist amongst the qubits . quantum correlations , especially quantum entanglement , is one of the crucial resources in quantum information theory and has been studied extensively @xcite . the indispensable role of entanglement in quantum information processing tasks such as quantum teleportation @xcite , superdense coding @xcite and remote state preparation @xcite @xcite , etc . has been quantified and understood deeply . however , despite showing that quantum entanglement is _ necessary _ for the pure state implementation of grover search @xcite @xcite and shor s factoring algorithm @xcite , the innate role of entanglement is not clear for general quantum computing tasks . one of the fundamental attributes of quantum systems is their ability to exist in linear superpositions of different physical states . this physical phenomenon is called quantum superposition . quantum coherence arises from superposition and is at the heart of several quantum features like quantum interference , multiparticle entanglement @xcite , quantum biology @xcite @xcite @xcite , quantum thermodynamics @xcite @xcite @xcite @xcite , quantum game theory @xcite @xcite @xcite etc . , which in turn are some of the most important applications of quantum physics and quantum information science . inspired from the resource theory of entanglement @xcite @xcite , there has been a lot of effort recently to quantify coherence as a resource theory @xcite . the role of coherence in the deutsch - jozsa algorithm has also been explored recently @xcite . in a standard resource - theoretic treatment of coherence , the `` incoherent '' states are those that are diagonal in some fixed reference basis . the amount of coherence is then defined as the distance from these reference states . a ` coherence measure ' is a real - valued function over the quantum state - space , such that it vanishes for all the states that are deemed to be incoherent and can not increase under some class of operations that preserve incoherence .    in this paper , we quantitatively analyze the role of quantum coherence and monogamy of entanglement @xcite @xcite in the discrete analogue of the analog grover search algorithm . to quantify the dynamics of coherence in this analog grover search , we use two coherence monotones namely the @xmath1-norm and the relative entropy of coherence @xcite . we then discretize the analog grover search algorithm to study how the monogamy of quantum entanglement evolves parallel to the system s adiabatic evolution . we use an entanglement monotone called the concurrence  @xcite to study the same . we find that the amount of coherence ( as quantified by the coherence measures ) is non - zero at all times during the search . the only time this coherence goes to zero is when the search algorithm attains the maximum success probability equal to one . since we begin with a maximally coherent state , in a resource theoretic sense , @xmath0-maximal coherence is actually consumed during the quantum search . using coherence monotones , we relate the success probability of the search algorithm with the amount of coherence . further , during the analog search , the final state is a product state and so all entanglement monotones must go to zero at the end of the search algorithm . the bipartite entanglement entropy , concurrence and the monogamy score , all behave in accord with this , since for both qubit partitions , these measures of quantum correlations go to zero and mark the completion of the search algorithm . furthermore , the two qubit concurrence and the monogamy score peak simultaneously with the increasing rate of success probability , hence implying that entanglement monogamy is indeed satisfied for the discrete analogue of the analog analogue of grover search algorithm . the paper is organized as follows . in sec . [ sec : review ] , we review the analog grover search algorithm . in sec . [ sec : coherence ] , we use the two coherence monotones to calculate the dynamics of coherence and elucidate their behavior . we also define the mapping through which we discretize the analog grover search algorithm . in sec . [ sec : entanglement ] , we explore the entanglement entropy and the concurrence , and demonstrate the synonymous behavior between the rate of change of concurrence in the continuous and the discrete grover search . in the next section , we calculate the two - qubit concurrence and show its connection to the success probability of the search algorithm . we then use the monogamy inequality to bound the bipartite entanglement at all times during the search . we conclude our paper with the discussions in sec . [ sec : conclusion ] . adiabatic quantum computation can be described as a controlled hamiltonian evolution of a system obeying the schrdinger equation @xmath2 we briefly recapitulate the analog analogue of grover search algorithm  @xcite , where the problem is to use quantum evolution to find a marked state among @xmath3 orthonormal states . + imagine that we are given a hamiltonian in an @xmath3 dimensional hilbert space such that it has only one non - zero eigenvalue , @xmath4 . the task then is to find the corresponding eigenvector @xmath5 which has eigenvalue @xmath6 . the hamiltonian can be represented as @xmath7 with @xmath5 unknown and normalized , i.e. , @xmath8 . now we choose some normalized vector @xmath9 , which is independent of @xmath5 ( since we do not know what @xmath5 is yet ) . we then add to @xmath10 , the  driving " hamiltonian @xcite @xmath11 so that the full hamiltonian is given by @xmath12 now , starting from the initial state @xmath9 at @xmath13 , we calculate the time evolution of the state @xmath14 . since the total hamiltonian is time - independent , we can simply write the time evolved state as @xmath15 we work in the units @xmath16 . it suffices to confine our attention to the two dimensional subspace spanned by @xmath9 and @xmath17 . the vectors @xmath9 and @xmath17 are not orthogonal ( in general ) and let us denote their inner product as @xmath18 , i.e. , @xmath19 , where @xmath18 can be taken to be real and positive since any phase in the inner product @xmath20 can be absorbed in @xmath9 . now the vectors @xmath21 and @xmath17 are orthonormal . in the \{@xmath17 , @xmath22 } basis , the hamiltonian is given by @xmath23\ ] ] and @xmath24.\ ] ] now the state of the system at time @xmath25 is found to be @xmath26 .\ ] ] thus , we can see that at time @xmath25 , the probability of finding the state @xmath17 is given by @xmath27 and the probability is one at time @xmath28 given by @xmath29 the inner product between the vectors @xmath17 and @xmath9 ( defined as _ x _ above ) , is assumed to be non - zero , else it will take infinite time for the quantum search to complete . this analog model has been generalized further . for example , the same algorithm was recast in the form of a spatial lattice search problem @xcite , where there is an @xmath3-dimensional lattice and the basis state @xmath30 is localized at the @xmath31th lattice site . the on - site potential energy , @xmath6 is zero everywhere except at @xmath32 , where it takes the value 1 . the objective is same as before , to reach the marked state @xmath32 starting from an equal superposition of all the @xmath33 s . the kinetic term is formulated through the laplacian of the lattice , which effectively introduces uniform hopping to all the nearest neighbors from any given lattice site , and is kept constant . our results will also hold in these generalized models . quantum coherence is one of the salient features of the quantum world . as this crucial feature drives several quantum technologies , it is very desirable to quantify the usefulness of coherence as a resource . this is done using the mathematical framework called a ` resource theory ' . to characterize something as a resource , we must first impose certain constraints on what specific physical operations are allowed ( e.g. the local operations and classical communication(locc ) framework for quantum entanglement @xcite restricts one from performing joint quantum operations between distant laboratories ) , which define the freely accessible operations . to be able to execute general quantum operations under such a constraint then requires some `` special '' quantum states that contain a relevant resource ( e.g. entangled states ) and can be consumed in the process . in fact , a quantitative resource theory for entanglement is already in place @xcite @xcite and has later expanded to encompass a wider range of operational phenomenon @xcite @xcite . however , for a long time there did not exist a resource theoretic framework to quantify the physical aspects of coherence . coherence was often defined in terms of functions of the off - diagonal entries of a density matrix , and the definition justified on the basis of physical intuition , which in turn also lead to the idea of decoherence developing along similar lines @xcite . eventually , in 2014 , baumgratz _ et al_. @xcite defined a resource theory of coherence and rigorously quantified the role of coherence in close analogy to the resource theory of entanglement . let us formally introduce the measures of coherence in the framework of a resource theory that is based on the set of incoherent operations and incoherent states @xcite . coherence is a basis dependent property and hence first , we must fix a _ reference basis_. the choice of the reference basis may be dictated by the underlying physics of the problem ( say the energy eigenbasis ) or the task to be performed ( wherever we wish to use quantum coherence ) . let \{@xmath34 } be a basis for @xmath35 , a @xmath36-dimensional hilbert space . the set of all density matrices that are diaogonal in this chosen basis form the so called `` incoherent states '' . the set of all the states of the form @xmath37 , where @xmath38 and @xmath39 , forms the set @xmath40 of all incoherent states . any quantum state that does not belong to the set @xmath40 will be called a coherent state and will act as a resource . a maximally coherent state in this basis is then given by @xmath41 , as any other state can be created from @xmath42 using only the set of ` incoherent operations ' . et al_. @xcite define incoherent operations as the set of trace preserving completely positive maps @xmath43 admitting a set of kraus operators @xmath44 @xcite such that @xmath45 ( trace preservation ) and , for all @xmath0 and @xmath46 , @xmath47 } \in \mathcal{i}.\ ] ] the definition of the set of incoherent operations is not unique and so we redirect the reader to a recent review for further discussions @xcite . we shall use two coherence monotones introduced in @xcite , namely , the @xmath1-norm of coherence and the relative entropy of coherence .    given a density matrix @xmath48 for the system , the @xmath1-norm of coherence is defined as @xmath49    for the analog grover search , the @xmath1-norm of coherence for the state given in eq . ( [ eqn : wavefunction ] ) is found to be    @xmath50    similarly , the relative entropy of coherence is given by @xmath51 which is found to be    & _ r ( ) = - ( 1 - x^2 ) ^2(e x t ) _ 2 ( ( 1- x^2 ) ^2(e x t ) ) & + & -(x^2 ^2(e x t)+^2(e x t ) ) _ 2 ( x^2 ^2(e x t)+^2(e x t ) ) . & -norm , the relative entropy of coherence and the probability of success as a function of time at parameter values @xmath52 and @xmath18 = 0.707 . ] the @xmath1-norm of coherence has a closed form expression in terms of the probability of success given by , @xmath53 where the form of @xmath54 is given in eq . ( [ eqn : successprobability ] ) . similarly , the relative entropy of coherence expressed in terms of the probability of success takes the form of a binary entropy function , i.e. , @xmath55    it can be observed analytically from eq . ( [ eqn : l1intermsofp ] ) and eq . ( [ eqn : relintermsofp ] ) and numerically from the fig . [ fig : relcoherence ] that both coherence measures go to zero _ iff _ the probability of success peaks to one , hence acting as a clear signature for the success of the search algorithm . even for other general initial states ( i.e. , with a different value of @xmath18 ) , after evolving under the analog grover search hamiltonian , the coherence attains a minima as the success probability peaks to 1 . note that although @xmath56 would also imply both coherence measures going to zero as can be observed from eq . ( [ eqn : l1intermsofp ] ) and eq . ( [ eqn : relintermsofp ] ) , the probability of success can not be zero since that would require @xmath18 to be zero ( which is excluded else the search would take infinite time to complete , see eq . ( [ eqn : time ] ) ) .      to quantify coherence in the analog grover search algorithm , we consider a mapping through which the @xmath3 eigenstates of the system are mapped to the logical states of @xmath0 qubits ( where @xmath57 ) . the hamiltonian initially acting on the state @xmath58 of the system now acts on the @xmath0 qubit system . throughout this paper , we refer to this mapping between the @xmath3 eigenstates and the @xmath0 qubits as the _ discrete analogue of analog grover search algorithm_. + now , consider the initial state to be an equal superposition of all the basis states , i.e. , + @xmath59 and the final state @xmath17 is one of the orthonormal basis states . the @xmath60 can be expressed in the \{@xmath61 } basis as @xmath62 where @xmath63 , @xmath64 and @xmath65 . + therefore , the density matrix in the \{@xmath61 } basis is @xmath66    after discretizing the analog grover search , the @xmath1-norm of coherence for the density matrix in eq . ( [ eqn : densitymatrix ] ) is analytically calculated to be    @xmath67    similarly the relative entropy of coherence is calculated to be    @xmath68    the @xmath1-norm and the relative entropy of coherence display a similar behavior as before , and it is important to note that the initial state for the analog grover search algorithm is an equal superposition of all orthonormal basis states , i.e. , @xmath69 which corresponds to a maximally coherent state in this eigenbasis . however , at time @xmath70 ( which is the time at which the search algorithm succeeds ) , the coherence is reduced to zero . this implies that maximally coherent states are actually consumed during the search algorithm . entanglement entropy is a measure of entanglement for many body quantum states . bipartite entanglement entropy is defined with respect to a bipartition of a state into two parts say @xmath71 and @xmath72 . for this , consider a quantum system whose hilbert space @xmath73 . the bipartite entanglement entropy @xmath74 of a state @xmath75 is defined as the von neumann entropy of either of its reduced states . that is , for a pure state @xmath76 , the entanglement entropy is given by + @xmath77 +   + where @xmath78 and @xmath79 are the reduced density matrices for each partition . the von neumann entanglement entropy in the eigenbasis is defined as @xmath80    it is easier to quantify the role of entanglement in the original grover search due to its discrete nature , however it is not so in the analog version . therefore , we use the discrete mapping to calculate the entanglement entropy for the analog grover search . we partition the @xmath0 qubit state into one qubit ( say @xmath81 ) vs the other @xmath82 qubits . the reduced density matrix for the @xmath81th qubit , @xmath83 is obtained by tracing out the other @xmath84 qubits and is given as    @xmath85    the eigenvalues of the matrix , @xmath86 are given by @xmath87 from which one can find the @xmath88 . concurrence is an entanglement measure , which for the special case of a pair of qubits is closely related to the entanglement of formation @xcite . for an arbitrary two qubit density matrix @xmath48 , concurrence is defined as follows : we first define a `` spin - flipped '' density matrix , @xmath89 as @xmath90 where @xmath91 is the pauli matrix @xmath92 . then we calculate the square root of the eigenvalues of the matrix ( @xmath93 ) and arrange them in decreasing order ( say @xmath94,@xmath95,@xmath96,@xmath97 ) . the concurrence is then defined as @xmath98 for pure states , the concurrence is also defined via another quantity known as the `` tangle '' , where @xmath99 and the @xmath100 = @xmath101 . concurrence for the pure state ( in eq.([eqn : wavefunction ] ) ) across the bipartition @xmath102 is found to be @xmath103    ) and the probability of success as a function of time at parameter values @xmath52 and @xmath104 . ] the entanglement entropy and concurrence across the system partition of @xmath102 are calculated analytically and their variation with time along with the probability of success is shown in fig . [ fig : entanglement ] . note that both entanglement measures go to zero as the success probability peaks to 1 . rate of change of the concurrence with time is found to be @xmath105    one can see that @xmath106 by neglecting o(@xmath107 ) terms for n @xmath108 1 . therefore , for a large database size ( n @xmath108 1 ) , the rate of concurrence goes to zero , hence implying that a very small amount of entanglement is generated during the search process . [ fig : dcbydt ] shows the analogous rate of change between the discrete and the continuous grover search and suggests that the mapping chosen in subsection [ sec : mapping ] preserves the properties of the original grover search algorithm . +      the amount of quantum correlations that can be shared amongst the subsystems of a multipartite quantum state is captured by the idea of monogamy . the central idea of monogamy is that entanglement can not be freely shared . precisely put , _ if two qubits a and b are maximally entangled then they can not be entangled at all with a third qubit c _ @xcite .    for the tripartite case , if @xmath109 is a bipartite quantum correlation measure , then this measure is said to be monogamous ( or satisfy monogamy ) for a tripartite quantum state @xmath110 , if the following condition holds @xmath111 here @xmath112 is the quantum correlation ( with respect to the correlation measure @xmath113 ) between subsystems @xmath71 and @xmath72 , @xmath114 is the quantum correlation between subsystems @xmath71 and @xmath115 , and @xmath116 is quantum correlation between subsystem @xmath71 and subsystems @xmath72 and @xmath115 taken together . for example , if we have a three - qubit pure state @xmath117 and a quantum correlation measure say concurrence then it is known that @xmath118 even though @xmath119 is a two qubit subsystem with a four dimensional subspace , it can be shown that the support of @xmath120 is spanned by the eigenvectors corresponding to _ at most _ two non - zero eigenvalues of the reduced density matrix @xmath120 , and hence it effectively becomes a two - dimensional space . this allows one to treat the bipartition of @xmath71 and @xmath119 as an effective two - qubit system whose concurrence , @xmath121 , is simply given by @xmath122 . substituting the value of @xmath121 above , we get the coffman , kundu , wootters ( ckw ) inequality @xcite @xmath123 which suggests that the concurrence is a monogamous entanglement measure . + this also leads to the concept of quantum monogamy score , which , for a given bipartite quantum correlation measure , is defined as @xmath124 if a tripartite state satisfies entanglement monogamy then its monogamy score is positive or else it is negative . + another measure of bipartite entanglement is the entanglement of formation ( eof)@xcite , which is closely related to two - qubit concurrence . consider a bipartite quantum state @xmath125 , and the ensemble @xmath126 denoting a possible pure state decomposition of @xmath125 , satisfying @xmath127 . the eof is defined as @xmath128 where @xmath129 is the von neumann entropy of the reduced density matrix corresponding to the @xmath71 subsystem of @xmath125 . for a two - qubit mixed state @xmath125 , @xmath130 , where @xmath131 is the binary entropy function . the eof , being a concave function of squared concurrence , does not obey the ckw inequality . however , the square of the eof does obey the same relation as the squared concurrence for tripartite systems @xcite . the discrete analogue of analog grover search algorithm satisfies the squared concurrence monogamy calculated as @xmath132    the reduced density matrix for two qubits is calculated by tracing out the @xmath133 qubits    @xmath134    the concurrence of an arbitrary two - qubit state is calculated according to the formula in subsection [ subsec : concurrence ] , using the spin flipped qubit . the eigenvalues for the density matrix @xmath93 ( i.e , the density matrix obtained after multiplying with the spin - flipped qubit for @xmath125 ) are +   + @xmath135 , @xmath136 , @xmath137 and @xmath138 . therefore , the two qubit concurrence is given by    @xmath139    this is the pairwise entanglement in the analog grover search . the evolution of the pairwise entanglement is calculated numerically and the result is shown in fig . [ fig : concurrenceandrate ] along with the rate of success probability in the search algorithm . the two peak simultaneously suggesting that entanglement is indeed necessary for the discrete analogue of analog grover search algorithm . for the multipartite system , in particular , the pairwise entanglement sharing and other pairwise correlations are monogamous ; when @xmath0 tends to infinity all of the pairwise entanglement vanishes as seen from eq . ( [ eqn : twoqubitconc ] ) .    ) and the rate of success probability as a function of time at parameter values @xmath52 and @xmath104 . ]    since in the analog grover search case , the concurrence between any two pair of qubits is the same , i.e. , @xmath140 ; as a result , the monogamy score between the @xmath0 qubits reduces to the following @xmath141    the monogamy score is , thus , given by @xmath142     and @xmath3 = 32 . ] similarly , the squared entanglement of formation satisfies a monogamy inequality . we do not provide the expression for this as it is simply too long . monogamy inequality bounds the amount of pairwise entanglement that can be shared between multiple qubits , and we can see that the discrete analog of grover search algorithm satisfies two monogamy inequalities for arbitrary times . to summarize , in this paper we have explored the role of quantum coherence and monogamy of entanglement in the discrete analogue of the analog analogue of grover search algorithm . using the @xmath1-norm and the relative entropy of coherence , it was shown that coherence acts as a signature for the success of the analog grover search algorithm . it was also shown that a maximally coherent state ends up into an incoherent state as the search algorithm evolves and hence @xmath0-maximal coherence is actually consumed during the search process . the variation of entanglement was also quantified and the analogous rate of change of concurrence between the discrete and analog grover search algorithms suggests that our mapping preserves the original behavior of the algorithm . the pairwise entanglement was shown to peak simultaneously with the rate of success probability as evidence that entanglement is indeed necessary , for the pure state implementation of analog grover search algorithm . the pairwise entanglement also suggested a monogamous behavior of quantum correlations in the analog grover search and it is then shown that the discrete analogue of analog grover search satisfies the entanglement monogamy inequality for both entanglement measures namely the concurrence and the squared entanglement of formation , for all times during the search algorithm . + note : after the completion of this work , the authors noticed the paper @xcite , where similar results about coherence have been obtained independently by hai - long shi , si - yuan liu , xiao - hui wang , wen - li yang , zhan - ying yang and heng fan . namit anand would like to acknowledge the hospitality of the harish - chandra research institute for allowing him to use their facilities during several visits made as a summer student over the last year during the preparation of this manuscript . 100 r. p. feynman , int . * 21 * , 467 ( 1982 ) . d. deutsch , proc . london a * 400 * , 97 ( 1985 ) . d. deutsch and r. jozsa , proc . london a * 439 * , 553 ( 1992 ) e. bernstein and u. vazirani , siam j. comput . * 26*(5 ) , 1411 ( 1997 ) . d. simon , siam j. comput . * 26*(5 ) , 1474 ( 1997 ) . l. grover , phys . lett . * 79 * , 325 ( 1997 ) . p. shor , siam j. comput . * 26*(5 ) , 1484 ( 1997 ) . c. zalka , phys . rev . a * 60 * , 2746 ( 1999 ) . g. brassard , p. hyer , m. mosca , and a. tapp , in _ quantum computation and information _ , ams contemporary mathematics series vol . * 305 * , 53 ( 2002 ) . g. brassard , p. hyer , and a. tapp , sigact news * 28 * , 14 ( 1997 ) . c. drr , m. heiligman , p. hoyer , and m. mhalla , siam j. comput . * 35 * , 1310 ( 2006 ) . a. ambainis and r. palek , symposium on theoretical aspects of computer science , lecture notes in computer science vol . * 3884 * ( springer , berlin ) , 172 ( 2006 ) . e. farhi , j. goldstone , s. gutmann , and m. sipser , https://arxiv.org/abs/quant-ph/0001106[arxiv:quant-ph/0001106 ] . m. born and v. fock , zeitschrift fr physik * 51 * , 165 ( 1928 ) . s. jansen , m. ruskai , and r. seiler , j. math . phys . * 48 * , 102111 ( 2007 ) . d. aharonov , w. dam , j. kempe , z. landau , s. lloyd , and o. regev , siam j. comput . * 37*(1 ) , 166 ( 2007 ) . e. farhi and s. gutmann , phys . a * 57 * , 2403 ( 1998 ) . r. horodecki , p. horodecki , m. horodecki , and k. horodecki , rev . * 81 * , 865 ( 2009 ) . c. bennett , g. brassard , c. crpeau , r. jozsa , a. peres , and w. wootters , phys . * 70 * , 1895 ( 1993 ) . c. bennett and s. wiesner , phys . lett . * 69 * , 2881 ( 1992 ) . a. k. pati , phys . a * 63 * , 014302 ( 2000 ) . c. bennett , d. divincenzo , p. shor , j. smolin , b. terhal , and w. wootters , phys . * 87 * , 077902 ( 2001 ) . s. l. braunstein and a. k. pati , quant . info . comput . * 2 * , 399 ( 2002 ) . a. k. pati and s. l. braunstein , j. of indian inst . of sci . , * 89 * , 295 ( 2009 ) . a. ekert and r. jozsa , rev . phys . * 68 * , 733 ( 1996 ) . d. meyer and n. wallach , j. math . phys . * 43 * , 4273 ( 2002 ) . n. lambert , y. chen , y. cheng , c. li , g. chen , and f. nori , nat . phys . * 9 * , 10 ( 2013 ) . d. abbott , p. davies , and a. k. pati , quantum aspects of life , imperial college press , london , ( 2008 ) . s. huelga and m. plenio , contemp . * 54 * , 181 ( 2013 ) . j. goold , m. huber , a. riera , l. rio , and p. skrzypczyk , j. phys . a : math . theor . * 49 * , 143001 ( 2016 ) . g. gour , m. mller , v. narasimhachar , r. spekkens , and n. yunger halpern , phys . , * 583 * , 1 ( 2015 ) . d. janzing , int . . phys . * 39 * , 2717 ( 2000 ) . e. ruch and a. mead , theor . acta * 41 * , 95 ( 1976 ) . d. meyer , phys . 82 * , 1052 ( 1999 ) . j. eisert , m. wilkens , and m. lewenstein , phys . lett . * 83 * , 3077 ( 1999 ) . n. anand and c. benjamin , quant . * 14 * , 4027 ( 2015 ) . plenio and s. virmani , quant . * 7 * , 1 ( 2007 ) . t. baumgratz , m. cramer , and m. plenio , phys . lett . * 113 * , 140401 ( 2014 ) . m. hillery , phys . a * 93 * , 012111 ( 2016 ) . b. terhal , ibm j. res * 48 * , 71 ( 2004 ) . h. s. dhar , a. k. pal , d. rakshit , a. sen de , and u. sen , https://arxiv.org/abs/1610.01069[arxiv:1610.01069 ] . w. wootters , phys . lett . * 80 * , 2245 ( 1998 ) . a. childs and j. goldstone , phys . a * 70 * , 022314 ( 2004 ) . m. nielsen and i. chuang , quantum computation and quantum information , cambridge university press , cambridge , ( 2004 ) . e. chitambar , d. leung , l. maninska , m. ozols , and a. winter , commun . * 328 * , 303 ( 2014 ) . b. coecke , t. fritz , and r. spekkens , inf . & comp . * 250 * , 59 ( 2016 ) . m. horodecki and j. oppenheim , int . b * 27 * , 1345019 ( 2013 ) . m. schlosshauer , rev . phys . * 76 * , 1267 ( 2005 ) . a. streltsov , g. adesso , and m. b. plenio , https://arxiv.org/abs/1609.02439[arxiv:1609.02439 ] . v. coffman , j. kundu , and w. k. wootters , phys . a * 61 * , 052306 ( 2000 ) . k . bai , y .- f . xu , and z. d.wang , phys . * 113 * , 100503 ( 2014 ) . shi , s .- y . liu , x .- h . wang , w .- l . yang , z .- y . yang , and h. fan ,  coherence depletion in the grover quantum search algorithm " , https://arxiv.org/abs/1610.08656[arxiv:1610.08656 ] .
Abstract:

Answer: grover s search algorithm is the optimal quantum algorithm that can search an unstructured database quadratically faster than any known classical algorithm . the role of entanglement and correlations in the search algorithm have been studied in great detail and it is known that entanglement between the qubits is necessary to gain a quadratic speedup , for pure state implementation of the grover search algorithm . here , we systematically investigate the behavior of quantum coherence and monogamy of entanglement in the discrete analogue of the _ analog analogue of grover search algorithm_. the analog analogue of grover search is a continuous time quantum algorithm based on the adiabatic hamiltonian evolution that gives a quadratic speedup , similar to the original grover search algorithm . we show that the decrease of quantum coherence , quantified using various coherence monotones , is a clear signature of attaining the maximum success probability in the analog grover search . we also show that for any two qubit reduced density matrix of the system , the concurrence evolves in close vicinity to the increasing rate of success probability . furthermore , we show that the system satisfies a @xmath0-party monogamy inequality for arbitrary times , hence bounding the amount of @xmath0-qubit entanglement during the quantum search .


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: in this work we study the ( real ) eigenfunctions of the laplace - beltrami operator @xmath3 on a riemannian surface @xmath4 with dirichlet boundary conditions ( if @xmath4 has boundaries ) . consider a real eigenfunction which satisfies @xmath5 the nodal domains are the maximally connected domains in @xmath4 where @xmath6 has a constant sign . the nodal set ( or the set of nodal lines ) is the zero set : @xmath7 , which also forms the boundaries of the nodal domains . we shall denote the nodal count ( i.e. the number of nodal domains ) of @xmath8 by @xmath9 . + the investigation of quantum signatures of classical chaos and integrability has been a hot topic in quantum chaos for a long time @xcite . in the past few years , the interest in nodal domains , their counting and their morphology increased after blum _ et al _ @xcite proposed a quantitative method which distinguishes between the distributions of nodal counts in domains where the underlying classical dynamics is integrable ( separable ) or chaotic . this added a new approach , the statistical investigation of nodal patterns , to the more common investigation methods of spectral or wavefunction statistics , which are often connected to random - matrix theory @xcite . et al _ showed that if @xmath9 is the nodal count of the @xmath10th energy eigenstate of a domain , then @xmath11 has a limiting distribution @xmath12 , where the characteristics of the distribution depend on the classical properties of the domain . for separable domains , @xmath13 has a square root singularity at a ( system - dependent ) maximum value , while for chaotic systems @xmath13 is ( approximately ) normally distributed . comparison between the numerical results for chaotic billiards and the random - wave ensemble supports berry s conjecture @xcite - wave functions in a chaotic system behave in the limit of high energy like a random superposition of plane waves . by that , the qualitative observation of miller _ et al _ @xcite , that nodal sets can be used to distinguish between wave functions in chaotic and integrable domains , could be tested in a quantitative way . other studies of various quantities - which pertain to the morphology and complexity of the nodal network - were published in the mathematical and physical literature , building upon the older results regarding the bounds on the total lengths of the nodal lines and their curvature @xcite . e.g. , in @xcite the distribution of the curvature is calculated , in addition to the mean and the variance of the total length of the nodal set . the distribution of the avoidance distances between nodal lines was also computed @xcite , to mention few examples . + an important breakthrough has been achieved by bogomolny and schmit @xcite who implemented a critical percolation model that explains the large scale structure of nodal domains in chaotic wave functions . this model is supported by a variety of numerical calculations . for example : the expectation value and variance of the nodal count for chaotic billiards , as well as the distribution of areas of nodal domains , follow the predictions of the model @xcite ; the nodal lines in the high energy limit seem ( on large scales ) to be sle@xmath14 curves @xcite as it is proved for the boundaries of percolation clusters @xcite . despite the good agreement , the percolation description is a priori insensitive to the structure of the nodal set in scales of the order of a wavelength . in addition , it was demonstrated by foltin _ et al _ @xcite that there are some special measures with a scaling behavior which is different for percolation and the nodal set of the random - wave ensemble . the latter special measures , in general , probe subwavelength scales at two points at a ( large ) distance .    in this work we suggest a new ( quantum mechanical ) method for the classification of billiards according to their classical properties . we will discuss below in what sense the signatures in nodal patterns differ from the scenario known for the more common spectral and wavefunction analysis . + our method provides yet another test to the conjectures by berry and bogomolny . the parameter which we use in order to interrogate the morphology of nodal lines is defined as follows - we consider the @xmath10th eigenfunction of ( [ eq : helmholtz ] ) , and its nodal domains sequence @xmath15 . the indices @xmath16 specify a nodal domain ; for this domain we define the area - to - perimeter ratio @xmath17 by : @xmath18 where @xmath19 and @xmath20 are the area and perimeter of @xmath21 and the ratio is measured in units of the wavelength @xmath1 . we shall define for different ensembles two different probability measures on the parameter @xmath0 . + for wave functions which satisfy ( [ eq : helmholtz ] ) on a compact domain @xmath4 , we consider a spectral interval @xmath22 , \ g>0 $ ] with @xmath23 , and define : @xmath24 note that in the above , the weights of nodal domains which belong to the same eigenfunction are equal , but not necessarily the same as the weight of domains which belong to another eigenfunction . + the second probability measure pertains to an ensemble of wave functions on unbounded domains , in our case - the gaussian random - wave ensemble ( which will be described in section [ sec : random and chaotic ] ) . since the wave functions do not satisfy any boundary condition , we consider them over an arbitrarily large and fixed domain @xmath25 , and include only the nodal domains which are strictly inside @xmath26 . we denote their number for a given member of the ensemble by @xmath27 , and define : @xmath28 the reason for using two different measures is the different nature of the problems at hand . however in the limit the two measures coalesce . we shall investigate the existence and the features of a high energy limiting distribution @xmath29 the choice of the area - to - perimeter ratio @xmath0 as a parameter to characterize the geometry of nodal domains is is inspired by the following considerations : the nodal pattern for separable surfaces is a checker - board , where a nodal domain @xmath30 is asymptotically a rectangle with sides of the order of a wavelength . therefore @xmath31 and @xmath32 will be of the order of one . similarly , according to the percolation model , a nodal domain of a chaotic surface is asymptotically shaped as a chain @xcite with @xmath33 cells , where for each cell @xmath34 where @xmath35 is the cell s contribution to the nodal domain s perimeter ( see [ eq : a_to_l_perco ] ) . we get that in both cases the parameter @xmath0 for a typical nodal domain will be of the order of unity , yielding localized distributions for the two types of surfaces . however , as will be shown below , the distributions differ substantially for systems with different classical properties . + in addition , the area - to - perimeter ratio @xmath0 is relevant not only to the study of the high energy limit , but arises as a natural parameter in the study of isoperimetric inequalities ( see e.g @xcite ) . the restriction of a wave function @xmath6 to one of its nodal domains @xmath36 , is an eigenfunction of the laplace - beltrami operator on the domain @xmath36 , with dirichlet boundary conditions . since it consists of a single nodal domain , courant theorem @xcite implies that it is the ground - state of @xmath36 . therefore , knowing @xmath17 we can express the ground - state energy in terms of the area and perimeter of @xmath36 . in the mathematical literature there are known bounds for such expressions - a relevant example is the bound for convex domains , derived by makai @xcite ( lower bound ) and plya @xcite ( upper bound , which was generalised to all simply or doubly connected domains by osserman @xcite ) : @xmath37 in order to derive a distribution function @xmath2 between the extreme values , some measure on domains should be defined . here we confine ourselves to well defined families of domains - those obtained as nodal domains of a given ensemble , and due to this restriction we are able to define the measures ( [ eq : gendis1]),([eq : gendis2 ] ) and study the limiting distribution for different classes of systems . + in this paper we shall examine the distribution function @xmath2 for separable and chaotic domains and for the gaussian random - wave ensemble . we will show that    * the limit distributions we obtain have strong  universal  features . that is , they depend crucially on the type of classical dynamics the manifold supports , and only to a lesser extent on the idiosyncratic details of the actual system . * the limiting distribution for the random - wave ensemble is similar to the one for chaotic domains , as predicted by berry s conjecture . * the limiting distribution for random - waves ( chaotic billiards ) is consistent with the percolation model but contains ( universal ) information beyond percolation as short length scales on the order of a wavelength are probed for small nodal domains . the limitation to ( quantum mechanically ) separable systems is due to the checkerboard structure in the nodal patterns of their wavefunctions . a generalization to all integrable or pseudo - integrable domains , where , in general , the checkerboard structure is lost , would be desirable . this highlights a general difference in the scenario known e.g.  from spectral statistics where all integrable systems ( separable or not ) share the same ( poissonian ) statistics . quite contrary statistical properties of nodal domains in integrable systems are very different for seperable systems with a checkerboard structure and non - separable systems with generically no nodal crossings @xcite . note , that when a separable system is slightly perturbed , all nodal crossings will open ( in an often highly correlated way which makes the introduction of a percolation model at this point quite difficult ) and statistical properties will change singularly ( again in contrast to what is known from spectral statistics where a small perturbation smoothly changes the statistics ) . as was mentioned above , the nodal network of eigenfunctions of separable surfaces has a checkerboard structure . this follows from the fact that one can always choose a basis in which all the eigenfunctions can be brought into a product form . therefore , one might expect that the main features of the limiting distributions of different surfaces will be similar . we will show that this is the case , and therefore we begin this section by explicit calculation of @xmath38 for a rectangular billiard . the discussion of this simple example will pave the way to computing @xmath2 for other systems ( the disc billiard and a family of surfaces of revolution ) and to the identification of some common features which we assume to be universal for all the separable systems . the detailed computations are presented in [ app : a ] . + the dirichlet eigenfunctions for a rectangular billiard with side lengthes @xmath39 are : @xmath40 the corresponding eigenvalues are : @xmath41 \equiv e_m+e_n\end{aligned}\ ] ] where @xmath42 can be interpreted classically as the energy stored in each degree of freedom ( a formal definition can be found in [ app : a ] ) . the nodal domains are rectangles of size : @xmath43 , therefore : @xmath44 and @xmath45 where : @xmath46 . the mere form of ( [ eq : rhorec ] ) implies that @xmath0 is bounded by : @xmath47 thus , the support of the distribution function @xmath2 is an interval which is narrower than ( [ eq : mpbound ] ) . + in the limit @xmath48 , the distribution ( [ eq : gendis1 ] ) can be approximated ( neglecting corrections of order @xmath49 ) by an integral . performing the integration over the variables : @xmath50 the distribution function is : @xmath51     the limiting distribution for a rectangular billiard ( [ eq : prec ] ) compared to the calculated distribution for eigenfunction with @xmath52,width=384 ]    the explicit form of @xmath38 suggests the following qualitative and quantitative conclusions :    1 .   the existence of a limiting distribution function ( which is independent of the aspect ratio @xmath53 of the billiard ) is demonstrated . 2 . it is supported by the compact interval @xmath54 $ ] . @xmath38 is an analytic and monotonic decreasing function in the interval where it is supported . 4 .   @xmath55 . @xmath56 , hence , @xmath38 is discontinuous at both boundaries of the support . the fact that @xmath57 depends solely on the partition of the energy between the modes was a key element in the construction above , and plays a similar role in computing @xmath2 for the other separable systems . when @xmath58 , we get from equation ( [ eq : rhorec ] ) that @xmath59 , while @xmath60 means that the energy is concentrated completely in one degree of freedom . the concentration of probability near @xmath58 shows that equal partition of energy is prevalent among the nodal domains .    an explicit derivation of the limiting distribution @xmath2 for the family of simple surfaces of revolution ( following bleher @xcite ) can be found in [ app : a ] , in addition to a separate derivation for the disc billiard . in both cases it is proven that : @xmath61 where @xmath62 is given by ( [ eq : prec ] ) , and @xmath63 is a finite , positive and smooth function of @xmath0 . therefore the features which characterize @xmath38 dominate @xmath64 for all the systems considered ; thus @xmath64 is supported on the same interval and demonstrates the same type of discontinuities at its boundaries . + following the striking similarity of the distributions for all of the investigated manifolds , we suggest that properties ( i)-(v ) of @xmath38 which were derived for the rectangular billiard , are universal features of @xmath64 for all two - dimensional separable surfaces . we support this assumption by a heuristic model which is presented in [ app : a ] .    , @xmath65 and @xmath66 where @xmath67,width=384 ]    numerical simulations for the rectangle and the disc billiards for several energy intervals show good agreement with the analytic derivation ( see figures [ fig : prec],[fig : pdisc ] ) . for numerically obtained results at finite energies , two kinds of deviations from the limiting distributions can be observed :    1 . fluctuations along the entire range of @xmath0 ( for the disc ) or discrete jumps ( for the rectangle ) in the value of @xmath2 , which vanish in the limiting distribution due to the convergence of the corrections to the semi - classical approximations ( i.e. turning sums over quantum numbers into integrals and the neglect of terms of order @xmath68 ) . cusps near @xmath69 and @xmath70 for the disc : the origin for the appearance of these features is due to nodal domains with exceptional geometry - the inner domains of all wavefunctions ( which are asymptotically triangles ) for the former , the domains of @xmath71 ( which are ring shaped ) for the latter . it is verified ( analytically and numerically ) that these differences converge to zero as @xmath72 , or faster . while for chaotic wavefunctions there is no known analytic expression for the nodal lines , we will use known results about the morphology of the nodal set in order to propose some physical arguments for the expected distribution . the explanations we propose are all in agreement with numerical simulations - a detailed information about the numerical techniques and the reliability of the results can be found in appendix b. + a frequently used model for eigenfunctions in a chaotic billiard is that of the gaussian random - wave ensemble . this is based on a conjecture by berry @xcite that eigenfunctions of a chaotic billiard in the limit of high energies have the same statistical properties as the gaussian random - wave ensemble . + a solution @xmath73 for the helmholtz equation ( [ eq : helmholtz ] ) with a given energy @xmath74 on a given domain , can be written as a superposition of functions @xmath75 which span a complete basis , for example : @xmath76 since the solutions of ( [ eq : helmholtz ] ) are real , we are restricted ( for this choice of basis ) by @xmath77 . according to berry s conjecture , expanding the eigenfunctions of chaotic billiards ( in the high energy limit ) in terms of ( [ eq : randomwave ] ) , the coefficients @xmath78 distribute for @xmath79 as independent gaussian random variables with @xmath80 and therefore can be modeled statistically by this ensemble of independently distributed gaussian random - waves . +        as suggested by bogomolny and schmit @xcite , the nodal domains of a random wave are shaped as critical percolation clusters ( see figure [ fig : bond_percolation ] ) , where each site is of an average area @xmath81 where , as before , @xmath82 . the area ( or alternatively , the number of sites ) of the nodal domains ( see e.g. @xcite ) distributes ( asymptotically ) as a power law : @xmath83 where ( for 2d percolation ) : @xmath84 . for bond - percolation model over a lattice ( as illustrated in fig . [ fig : bond_percolation ] ) , the area of a cluster which spreads on @xmath33 sites is @xmath85where @xmath86 is the area of a single site and @xmath87 is the area of the connection between two sites ; the average perimeter is @xmath88where @xmath89 are the average contributions to the perimeter of a site and a connection ( since a cluster may contain loops which affect its perimeter , we must speak about average ) . the average area - to - perimeter ratio can be written as @xmath90 where @xmath91 . this relation can be used as a guideline to the desired distribution of @xmath0 . + indeed , as was confirmed numerically , the distribution @xmath2 for nodal domains follows ( [ eq : a_to_l_perco ] ) in several aspects . we have examined the restricted distribution for nodal domains with a given number of sites - we define @xmath92 to be the distribution for nodal domains with area @xmath93 we found out that @xmath92 is roughly symmetric about a mean value : @xmath94 . as in ( [ eq : a_to_l_perco ] ) , @xmath94 is increasing with @xmath33 , and converging to a limiting value : @xmath95 . since the percolation model is assumed to provide an exact description of the system in the high - energy limit , we expect ( [ eq : a_to_l_perco ] ) to serve as a good approximation for large domains ( see fig . [ fig : rho_of_n ] ) . + the value of @xmath96 is a direct result of a theorem by cauchy for the average chord length of a domain : @xmath97 where @xmath98 is the chord length . the original theorem ( which was stated for convex domains ) is extended in @xcite , to include nonconvex and multiply connected domains . + for nodal domains of infinite size , the statistics of the average chord length should follow that of the entire nodal set . that in turn is known to be @xcite : @xmath99 , therefore : @xmath100 the value of @xmath101 can also be estimated : as shown in @xcite , the single cell nodal domains are mild deformations of a circle of radius @xmath102 ( where @xmath103 is the first zero of @xmath104 ) . therefore : @xmath105     is plotted for several values of @xmath33 . the mean value is increasing with @xmath33 , while the variance decreases ( except for @xmath106 ) as suggested by the heuristic model.,width=432 ]     compared to the best numerical fitting ( for high values of @xmath33 ) to ( [ eq : a_to_l_perco ] ) : @xmath107.,width=336 ]    in order to study the impact of the deformations on the value of @xmath101 , we have calculated @xmath0 for a variety of domains , like ellipses , rounded shapes with corners ( a quarter of a circle or a stadium etc . ) and others . the results show that stretching of the nodal domain ( e.g.  increasing the eccentricity of an ellipse ) increases @xmath0 , while turning it  polygonal  ( i.e.  having points on the nodal line with very high curvature ) reduces @xmath0 . + derivation of  @xmath94   for other values of @xmath33 seems to be more complicated . however , fitting between the numerical results for high @xmath33 values and ( [ eq : a_to_l_perco ] ) equips us with the empirical result ( which is valid for @xmath108 ) : @xmath109 another interesting feature is the width of the distribution around @xmath94 . equation ( [ eq : chord ] ) implies that the variance is proportional to the variance in the average chord length between different nodal domains of the ( approximately ) same area . therefore , the variance is expected to be smaller for larger domains , which follows the statistics of the entire nodal network to a larger extent . the only exception is the variance for single site domains , which as was mentioned @xcite , have strong limitations on their shape , and therefore a relatively small variation in the average chord length . + the bounds ( [ eq : mpbound ] ) on @xmath0 should not hold in general for the nodal domains of ( [ eq : randomwave ] ) ; however , the numerical bounds seem to agree with ( [ eq : mpbound ] ) for all of the measured nodal domains , including multiply connected domains , suggesting that ( [ eq : mpbound ] ) is valid for the nodal domains of the ensemble with probability 1 .     calculated for the random - wave ensemble and for inner domains of a sinai and stadium billiards.,width=384 ] the distribution of @xmath0 for all of the nodal domains is given by : @xmath110 where @xmath111 is given asymptotically by ( [ eq : fisherexp ] ) . [ fig : rw_vs_inner ] shows the calculated distribution for 3 different systems - a random - wave ensemble , the inner domains of a sinai billiard and those of a stadium billiard . comparing the functions , we find additional strengthening to berry s conjecture . there are boundary effects of chaotic billiards - e.g. a peak in the distribution @xmath2 near @xmath112 and @xmath113 , and a lower probability for @xmath114 , however they vanish in the semi - classical limit ( see fig . [ fig : inner_vs_outer ] ) . in our study it was easy to put those effects aside - if we consider only inner nodal domains for the chaotic billiards ( as in fig . [ fig : rw_vs_inner ] ) , we observe no prominent differences between the distributions .    ) derived for the rectangle billiard and the random - wave ensemble.,width=384 ] the main results of this work can be summarized as follows :    1 . the distribution function @xmath2 of the area - to - perimeter ratio @xmath0 , distinguishes between billiards with separable or chaotic classical limit ( see fig . [ fig : rw_vs_separable ] ) . the distribution ( [ eq : pofrdsor ] ) for the examined separable billiards has some universal features , such as a common support and a square - root divergence at the lower support . in all studied cases the distribution is a mild deformation of the distribution that we found for a rectangle . 3 .   in accordance with the random - wave conjecture , we find numerically that chaotic billiards ( stadium and sinai ) have a universal limiting distribution @xmath2 , and it converges to the distribution found for the random - wave ensemble . by considering only the inner nodal domains for billiards , the agreement can be shown also for finite energies . 4 .   the numerical results suggest that for nodal domains of a random wave or of eigenfunctions of chaotic billiards , the area - to - perimeter ratio is bounded by ( [ eq : mpbound ] ) , i.e. @xmath115 , including nonconvex and multiply connected nodal domains ( for which these bounds have not been proven ) , with probability one . we examined the percolation model for the nodal set of random waves from the perspective of the area - to - perimeter ratio . it is shown that on the wavelength scale the geometry of the nodal domains can only be poorly characterized by percolation arguments . however , for large domains the geometry can be described by heuristic expressions like ( [ eq : a_to_l_perco ] ) , which are consistent with percolation theory . in this appendix we suggest a heuristic model for the universal features of the limiting distributions of the area - to - perimeter ratio for two - dimensional separable domains . the model is supported by an explicit derivation of the limiting distribution for the disc billiard and for simple surfaces of revolution . we begin by considering the classical geodesic flow in a two - dimensional compact domain ( e.g. a billiard ) . for a separable domain , a trajectory can be specified by its action - variables : @xmath116 where @xmath117 are the ( separable ) coordinates , @xmath118 are the conjugated momenta and the integration is over one period of the specified coordinate . at every point along the trajectory , the energy @xmath119 can be expressed as @xmath120 , where : @xmath121 where @xmath122 are the local unit vectors - if we consider circular domain for example , then @xmath123 . in general @xmath42 are not constants of motion . the only exceptions are the trajectories in a rectangular billiard .    from a quantum point of view , the eigenstates of schrdinger equation ( [ eq : helmholtz ] ) for a separable domain , can be written as @xmath124 , while the laplace - beltrami operator can be written as : @xmath125 , where @xmath126 . this allow us to define the quantum analogue to ( [ eq : parte ] ) : @xmath127 in the semi - classical limit , the spectrum of a separable domain is given by @xmath128 , where @xmath129 is the energy of the classical trajectory specified by the action - variables @xmath130 ( as emerging from bohr - sommerfeld quantization ) . in addition the semi - classical value of ( [ eq : parteq ] ) converges to the classical value ( [ eq : parte ] ) for every point in the domain . +     as a function of the partition of energy between the two degrees of freedom , for the rectangular billiard.,width=288 ]        in section [ sec : separable ] , the value of the area - to - perimeter ratio for a given realization @xmath131 for the rectangular billiard , was derived to be @xmath132 therefore , for a rectangular billiard , the value of the ( quantum ) parameter @xmath0 , has also an immediate classical interpretation ( see figures [ fig : rho_of_ei],[fig : p_of_em_rec ] ) . + in order to generalize the limiting distribution which was derived for a rectangular domain to other separable domains , we suggest the following heuristic model :    * in the high energy limit , ( almost all of ) the nodal domains of a separable domain are converging to rectangles , and the wave function in the close neighborhood of a nodal domain is converging to ( [ eq : psirec ] ) . therefore ( in the limit ) equation ( [ eq : rhorec2 ] ) should hold . however , since for general separable domains @xmath133 is not an eigenfunction of the operators @xmath134 , the value of @xmath135 for a given realization @xmath131 will depend on the interrogated nodal domain @xmath136 . * in the continuum limit , the limiting distribution is of the form : @xmath137 where the first integral is over the energy interval , and the second is over the domain . the function @xmath138 is the quotient of the appropriate jacobian and @xmath139 . since in the vicinity of @xmath140 two solutions for ( [ eq : rhorec ] ) coalesce , we expect the square root singularity at @xmath141 to be a universal feature . * in addition , since @xmath142 is convex ( see equation [ eq : rhorec2 ] and figure [ fig : rho_of_ei ] ) , the distribution function should be monotonically decreasing . this supports the assumption that the properties ( i)-(iv ) for @xmath143 and @xmath38 , which were derived in section [ sec : separable ] for the rectangle , are universal features of @xmath2 for all two dimensional separable surfaces . moreover , this model suggest - at least for separable domains - that a geometric feature of the nodal pattern , i.e. the area - to - perimeter ratio of a given domain , can be deduced directly from the underlying classical dynamics . + the suggested model is supported by an explicit derivation of the limiting area - to - perimeter distribution for several separable domains . in these calculations we approximate eigenfunctions and eigenvalues by the wkb method . in addition , we approximate sums over quantum numbers ( see eq . [ eq : gendis1 ] ) by integrals and neglect terms of order @xmath68 . the error resulting from these approximations is of the order of @xmath68 and therefore converge to zero in the limit . + the theme of the derivations is similar . the hamiltonian @xmath144 for these systems is homogeneous i.e. : @xmath145 this implies that the energy of the state @xmath131 can be expressed as : @xmath146 therefore , integration over the quantum number @xmath147 becomes trivial . we will also use the first term in the weyl series : @xmath148 in order to estimate @xmath149 . equation ( [ eq : helmholtz ] ) can be written in polar coordinates as @xmath150 for the disc billiard the boundary condition are : @xmath151 . the eigenfunction and eigenvalues of ( [ eq : phidisc ] ) are : @xmath152 where @xmath153 is an arbitrary phase and @xmath154 is the @xmath33th zero of @xmath155 . the nodal domains of @xmath131 will be @xmath156 replicas ( or one for @xmath157 ) of a slice containing @xmath33 domains ; we will enumerate them as @xmath158 , where @xmath159 is the most inner domain . the area and perimeter of @xmath160 are : @xmath161 an implicit semi - classical expression for @xmath154 can be deduced by applying the wkb approximation to ( [ eq : phidisc ] ) : @xmath162 where : @xmath163 setting @xmath164 , and substituting ( [ eq : bsdisc ] ) into ( [ eq : discfactors ] ) we get : @xmath165 which implies that @xmath166 depends on @xmath167 solely . since @xmath168 varies with @xmath33 as @xmath68 , we can approximate : @xmath169 . expressing @xmath170 in terms of ( [ eq : bsdisc],[eq : discfactors ] ) yields : @xmath171 bearing in mind that for a point @xmath172 : @xmath173 it can be shown that equation ( [ eq : rho(c ) ] ) is equivalent to ( [ eq : rhorec2 ] ) . + integrating ( [ eq : gendis1 ] ) over the variables @xmath174 , we get : @xmath175 performing the integration we get the limiting distribution : @xmath176 where : @xmath177    for the disc . @xmath178 is monotonic , finite and positive.,width=384 ]      the investigated surface of revolution @xmath4 is generated by a rotation of the analytic profile curve @xmath179 ( where @xmath180 $ ] ) around the @xmath181 axis . we restrict @xmath179 by : @xmath182 , which ensures smoothness at the poles of the surface . in addition , we request that @xmath183 , so @xmath179 has a single maximum : @xmath184 . + the lagrangian of the surface is given by : @xmath185 from which the action variables can be deduced : @xmath186 where @xmath187 , @xmath188 are the classical turning points which satisfy : @xmath189 . + as in [ subsec : disc ] the nodal domains will be @xmath190 copies of a slice with @xmath33 domains and will be denoted by @xmath191 . the homogeneity of the hamiltonian follows directly from ( [ eq : actionvar ] ) . we will follow the notations @xmath192 to get : @xmath193 the wkb approximation to the eigenfunctions is @xmath194 where @xmath195 in the limit of large @xmath33 , the nodal points density on the curve is high . therefore , applying the wkb approximation , successive nodal points @xmath196 , @xmath197 should satisfy : @xmath198 in addition , due to the homogeneity of @xmath144 , @xmath199 depend on @xmath167 solely , therefore : @xmath200 therefore , the location of the @xmath201 zero will be an ( implicit ) function @xmath202 and will not be depended on @xmath147 . the area and perimeter of the nodal domains are : @xmath203 therefore @xmath204 substituting ( [ eq : sorquantization ] ) in ( [ eq : rhosor ] ) we get : @xmath205 where @xmath206 at the turning points . since for @xmath172 : @xmath207 equation ( [ eq : rhosor ] ) is equivalent to ( [ eq : rhorec2 ] ) as well . + integrating ( [ eq : gendis1 ] ) over @xmath208 we get : @xmath209 setting @xmath210 , we get that for @xmath211 $ ] and fixed @xmath167 , there are two allowed values of @xmath212 : @xmath213 since @xmath214 is doubly valued , there are four allowed values of @xmath215 : @xmath216 for @xmath217 . therefore @xmath218 where @xmath219 @xmath220 is a monotonic increasing function , therefore the integrals in ( [ eq : psor3 ] ) diverges at @xmath221 . for @xmath222 , @xmath223 , therefore the integrals converge at infinity . + for @xmath224 , the coefficient of the integrals in ( [ eq : psor3 ] ) has a square root singularity , while the integral is converging to a finite positive value , therefore @xmath225 . + for @xmath226 , the first coefficient in ( [ eq : psor3 ] ) is of the order of @xmath227 . the lower limit of integration is defined by : @xmath228 consequently , the first term in ( [ eq : psor3 ] ) will be of order : @xmath229 the value of @xmath215 depends on the profile curve @xmath179 , however for @xmath230 , @xmath231 , while for @xmath232 , we get that @xmath233 , therefore @xmath234 is bounded for all possible values of @xmath215 . + the value of the second term is due to contributions of nodal domains which satisfy @xmath235 therefore , only eigenfunctions for which @xmath236 contribute to the second term , as a result , it is bounded by : @xmath237 + therefore @xmath238 is finite , and the universal features specified in section [ sec : separable ] are all fulfilled by ( [ eq : psor3 ] ) . in order to evaluate the area - to - perimeter ratios and their distribution for the random - wave ensemble and chaotic billiards , we have simulated the appropriate wavefunctions on a grid . + we have calculated the statistics for 5000 realizations of random waves , where in each realization we summed over 70 terms in ( [ eq : randomwave ] ) ; for chaotic billiards we have reproduced the first 2430 eigenfunctions of a sinai billiard and the first 2725 eigenfunctions of a stadium billiard . + the ( seemingly simple ) task of perimeter evaluation must be carried out carefully ; it can be shown that naive methods , like perimeter s pixels counting , produce an error which is independent of the sampling resolution . in order to avoid this error , we have approximated the nodal line as a polygon , where the vortices are calculated using a linear approximation . we have set the sampling resolution to contain 85 pixels along the average distance between two nodal lines ( @xmath239 ) . this resolution was proved to produce an error which is less then a percent . the measured perimeter is expected to be shorter then the real one , as we are approximating a curve by a polygon . + the accuracy of the method was tested by calculating the ratio between @xmath240-the total length of the nodal set , and the area . it is known that for the random - wave ensemble @xcite : @xmath241 the numerical values for this ratio were between : @xmath242 it should be noted that when we calculate the total nodal length , we have to calculate the perimeter of nodal domains at the edge of the grid . in many cases ( i.e. where the edge domains are very small ) the perimeter calculated for them is larger then the real value . it seems likely that the error due to this effect is of the order of the error due to polygonal approximation , therefore the two compensate each other , to yield a total error which is relatively small , of order @xmath243 . + * bibliography *
Abstract:

Answer: in an attempt to characterize the distribution of forms and shapes of nodal domains in wave functions , we define a geometric parameter - the ratio @xmath0 between the area of a domain and its perimeter , measured in units of the wavelength @xmath1 . we show that the distribution function @xmath2 can distinguish between domains in which the classical dynamics is regular or chaotic . for separable surfaces , we compute the limiting distribution , and show that it is supported by an interval , which is independent of the properties of the surface . in systems which are chaotic , or in random - waves , the area - to - perimeter distribution has substantially different features which we study numerically . we compare the features of the distribution for chaotic wave functions with the predictions of the percolation model to find agreement , but only for nodal domains which are big with respect to the wavelength scale . this work is also closely related to , and provides a new point of view on isoperimetric inequalities .