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: quantum tunneling is one of the major aspects of quantum mechanics that sets it apart from classical mechanics . the ability of a quantum system to enter a classically disallowed region and , through it , reach an energetically favorable region has inspired many applications . in this article , we are concerned with tunneling as it pertains to quantum annealing , which is a quantum algorithm that uses quantum effects , such as tunneling , to relax a system into a desired ground state . quantum adiabatic optimization ( qao ) @xcite is a form of quantum annealing that solves optimization problems by adiabatically changing the hamiltonian for a system from one with a known ground state to one with the optimization solution encoded in its ground state . many studies @xcite have sought to characterize what , if any , speed - up this algorithm has over classical alternatives such as simulated annealing or simulated quantum annealing . a popular belief among the qaocommunity is that at least part of quantum annealing s power comes from its ability to tunnel through barriers more easily than classical algorithms can go over them , but several recent results @xcite have shown that tunneling alone is not enough to ensure a quantum speed - up . nevertheless , generalized qaois universal for quantum computation @xcite , and an understanding of tunneling in a quantum annealing setting can shed light on the power and limitations of qao , even if tunneling alone is not enough for a quantum speed - up . to that end , this article studies quantum tunneling as part of qaoin the simplest possible situation : tunneling in position - space through a one - dimensional barrier . in a previous article @xcite , we showed that this simple problem is equivalent to barrier tunneling of @xmath0 qubits in a bit symmetric setting in the large @xmath0 limit .    despite the simplicity of one - dimensional position - space tunneling , we find that an asymptotic analysis of qaospectral gaps for @xmath2 , which give bounds on the adiabatic run time , fail to accurately describe the behavior of the tunneling for moderately large @xmath0 . we find that @xmath4 is needed to ensure that asymptotic expressions accurately describe behavior , to the point where polynomial and exponential run time scaling can be confused for lower @xmath0 . given the current and near - term sizes of quantum computer implementations @xcite , this failure of asymptotics casts doubt on the efficacy of using asymptotics to analyze the output of these systems or to predict their capabilities .    in section  [ sec : background ] we provide a brief review of qaoand introduce the details of the barrier tunneling problem to be studied . we mention both the original hypercube problem and the simpler continuous one - dimensional problem . section  [ sec : trough ] gives the first evidence that the asymptotic limit may not be accurate even for moderately large @xmath0 . we show that the failure of the asymptotics is the same in both the hypercube and one - dimensional cases , proving that the failure is in the tunneling problem itself and not in the approximation scheme from @xcite . we compare the exact spectral gap to its asymptotic approximation in section  [ sec : exact_scaling ] to see where the issues discussed in the previous section arise . this section also provides estimates for the @xmath0 necessary to see asymptotic behavior and explores what effect next - order corrections to the asymptotics have on the analysis . section  [ sec : scaling ] uses properties of the derivatives of log - log data to analyze how well a curve approximates a power law versus and exponential fit . furthermore , we show how easily finite data for tunneling could lead to a mislabeling of the scaling behavior of the algorithm . finally , section  [ sec : conc ] concludes and provides an outlook for the future of such asymptotic analysis . quantum annealing derives its name by comparison with classical simulated annealing , which starts a system at a high temperature and then slowly lowers the temperature , trying to relax the system into its ground state . quantum annealing works on the same principle except it uses quantum effects instead of temperature . as it is usually formulated @xcite , the quantum annealing hamiltonian is a linear time interpolation between two static hamiltonians : @xmath5 where @xmath6 is the time parameter , @xmath7 , divided by the total predetermined runtime of the algorithm , @xmath8 . @xmath9 is the initial hamiltonian and is usually chosen to have an easily prepared ground state . @xmath10 encodes the solution to an optimization problem in its ground state configuration . the initial formulation of qao@xcite set this in a hilbert space of @xmath0 qubits where @xmath9 is a uniform field in the negative @xmath11 direction : @xmath12 and @xmath10 is diagonal the in @xmath13 basis with a cost function encoded along its diagonal : @xmath14    these canonical choices for @xmath9 and @xmath10 make the the algorithm easier to create and analyze , but they also ensure that @xmath15 is stoquastic , meaning its off - diagonal terms are all non - positive . the universality proof @xcite for qaorequires non - stoquastic hamiltonians , so it is currently unclear how much power stoquastic hamiltonians have in qao . key to qaois that the algorithm must run adiabatically , where adiabatically means that the total time , @xmath8 , scales according to a combination of matrix norms of time derivatives of @xmath15 and the spectral gap , @xmath16 , between the ground state and first excited state . typically , the adiabatic condition is stated as adiabaticity is ensure if @xmath17}||{\frac{d \hat{h}(s)}{d s}}||/(g(s))^2 $ ] . the derivation of this version of the adiabatic theorem is subtly incorrect , which has led to more robust statements and derivations of an adiabatic condition @xcite , but all of these more accurate conditions claim that the worst - case bound on @xmath8 depends on matrix norms , which in our case are linear in @xmath0 , and the inverse of @xmath18 } g(s)$ ] . therefore , the majority of this paper focuses on the calculation of @xmath19 since this is the primary changing quantity that determines the worst - case bound on the run time of qao . our barrier tunneling problem is inspired by tunneling using @xmath0 qubits in a bit - symmetric setting . a thin spike version of this bit - symmetric problem was originally proposed by farhi et al . @xcite , and the general bit - symmetric barrier tunneling problem has been analyzed extensively by other groups @xcite . the final cost function @xmath20 is a function only of the hamming weight @xmath21 as @xmath22 where @xmath23 is a function centered around @xmath24 with width and height proportional to @xmath1 .    in a previous article @xcite , we showed that this hamiltonian can be mapped onto the symmetric @xmath25-dimensional subspace where it describes a single spin @xmath26 particle . furthermore , using a modified version of the villain transformation @xcite , we showed that this system can be mapped , in the @xmath2 limit , onto a semiclassical wave equation in one continuous variable @xmath11 . this mapping only works for low - lying energy states , but because we only care about the spectral gap between the ground state and first excited state , this restriction is acceptable . in @xcite , we restrict ourselves to a square barrier so that our final , semiclassical schrdinger equation becomes @xmath27\psi(x),\ ] ] where @xmath28 the wave function @xmath29 is the continuous version of the eigenvector components in the spin-@xmath30 problem , and @xmath31 is the eigenenergy . the variable @xmath11 is a rescaled and shifted version of the hamming weight , taken in the continuum limit , and the barrier has also been rescaled to have width @xmath32 and height @xmath33 . the constants can either be taken as arbitrary or set to @xmath34 and @xmath35 if we want a direct correlation to the qubit problem . while eq . [ eq : main_de ] was derived from the @xmath0 qubit problem , it can also stand on its own , describing the critical tunneling moment in a one - dimensional continuous system . in this context , we just have a particle moving in a one - dimensional quadratic well with a barrier in the middle . for most of this article , we treat the schrdinger equation in eq . [ eq : main_de ] as its own independent problem .    the continuous well differential equation , eq . [ eq : main_de ] , can be solved using exponentials and parabolic cylinder functions , @xmath36 , up to a transcendental equation for the eigenenergies : @xmath37 where @xmath38 , @xmath39 , and @xmath40 refers to even ( @xmath41 ) or odd ( @xmath42 ) energy solutions . in the limit of @xmath2 , eq . [ eq : trans_con ] can be solved for the lowest two energy levels , resulting in a gap expression @xmath43 for @xmath44 , and @xmath45 for @xmath46 . these asymptotic gap expressions match results for the original symmetric hypercube tunneling problem from other groups using various methods @xcite . in subsection  [ sub : problem ] we outline three steps in the approximation procedure from which the gap can be extracted . in the rest of the paper we will refer back to these three different models , so we will take a moment here to label these three models and define exactly what we mean by each term . * * discrete analysis * : the original @xmath0 qubit hypercube problem consists of a discrete hilbert space that can be interpreted in terms of hamming weight or spin states . the gap can be extract from this problem through exact diagonalization of the @xmath25-dimensional symmeterized problem . while the discrete problem was the original inspiration for this work , much of the subsequent text will focus exclusively on the continuous and asymptotic gaps . * * continuous analysis * : after the villain transformation , the discrete problem transforms into a parabolic well potential for a continuous , position - like variable . the schrdinger equation governing this continuous problem is shown in eq . [ eq : main_de ] , and the exact spectral gap can be extracted by numerically solving the transcendental equation for the energies , eq . [ eq : trans_con ] . for most of this article , the continuous problem will be treated as an independent problem , and most of our efforts will focus on comparing its exact spectral gap to the asymptotic limit . * * asymptotic analysis * : in the large @xmath0 limit , the asymptotic expression for the spectral gap of the discrete and continuous problems can be derived . we did this derivation in a previous paper @xcite , and the gap expressions are reproduced in eqs . [ eq : gap_poly ] and [ eq : gap_exp ] . this version of the gap is the true asymptotic limit of the discrete and continuous gaps , but much of this article will be focused on how large @xmath0 must be for the asymptotic limit to be a good approximation of the true gap . when analyzing an algorithm numerically or experimentally , a typical question to ask is how its run time , or in our case the gap , scales with the input size @xmath0 . therefore , in this section we analyze the gap numerically mostly ignoring our asymptotic expressions from eqs . [ eq : gap_poly ] and [ eq : gap_exp ] and try to determine what asymptotic information we would extract from numerical data of @xmath0 vs. @xmath19 . we then compare these finite @xmath0 extrapolations to see how close these results are to our analytic expressions for the asymptotic behavior of @xmath19 .     for the symmetric qubit tunneling problem with a binomial shaped barrier . each data point shows what scaling power , @xmath47 , is extracted from a fit to @xmath48 , assuming the fit is done on data from a finite range of @xmath0 versus the scaling power of the barrier size , @xmath49 . each type of point represents data binned from a different range of @xmath0 , and the dashed line represents the power @xmath47 predicted by the asymptotic expression . based on this data alone , the scaling of @xmath19 does not seem to approach the asymptotic expression . , scaledwidth=48.0% ]    we use symmetry to reduce the discrete qubit hamiltonian to a tridiagonal @xmath25-dimensional matrix and then numerically diagonalize it for a finite range of @xmath0 before performing power law fits of the form @xmath50 in order to extract the exponent @xmath47 . [ fig : trough_qubit ] outlines the resulting @xmath47 values for these fits done on simulations with @xmath44 , leading to polynomial scaling asymptotically .    also , note that the data in fig  [ fig : trough_qubit ] comes from a problem with a binomial shaped barrier rather than the normal square barrier we consider . to collect these data points , we resitricted ourselves to @xmath0 such that the discrete width of the barrier has just increased , so the binomial barrier s full width , which grows with @xmath51 , allows us to collect more data points than the square barrier , whose width grows with @xmath52 . therefore , we used the binomial barrier in order to collect more data points and create a cleaner data set .    in fig . [ fig : trough_qubit ] , we show @xmath47 values for fits done on different bins of @xmath0 data , so the red filled triangle should most closely approximate the asymptotic value of @xmath47 . unfortunately , the scaling powers do not seem to approach their asymptotic values , shown as a dashed black line in the figure , even for the relatively large problem size of one million qubits . the lack of agreement with the asymptotics , even for relatively large numbers of qubits casts doubt on the validity of the asymptotic analysis techniques used in @xcite .     for the continuous well de , eq . [ eq : main_de ] , with a square barrier . each data point shows what scaling power , @xmath47 , is extracted from a fit to @xmath48 , assuming the fit is done on data from a finite range of @xmath0 versus the scaling power of the barrier size , @xmath49 . each type of point represents data binned from a different range of @xmath0 , and the dashed line represents the power @xmath47 predicted by the asymptotic expression . the trough shows the same behavior as seen in fig . [ fig : trough_qubit ] , indicating that the problem in that plot was not in the approximations made for the villain transformation . , scaledwidth=48.0% ]    however , the only possible source of error in the asymptotic analysis in @xcite is in the villain transformation and the other assumptions that lead to the derivation of the continuous well de , eq . [ eq : main_de ] . therefore , if the spectral gap for the differential equation shows the same difficulty in reaching its asymptotic value , we know that extremely lare @xmath0 are necessary to observe asymptotic behavior because this is a simple tunneling problem , not because the villain transformation was used in the derivation of these asymptotics . we do in fact see that if we perform the same procedure that gave us fig . [ fig : trough_qubit ] for numerical @xmath19 values found from solving eq . [ eq : trans_con ] , we obtain a very similar plot , fig .  [ fig : trough_de ] .    the failure of @xmath19 to quickly reach its asymptotic scaling even in fig . [ fig : trough_de ] indicates that there is something special about simple tunneling problems that defies the conventional wisdom that asymptotic behavior sets in quickly . based on these figures alone , an experimental or numerical approach to this problem would mislabel the polynomial scaling of this algorithm . from this point , we focus exclusively on the continuous problem , ignoring the discrete problem that motivated it . in the next section , we examine what @xmath0 are necessary to bring the scaling close to the asymptotic expressions . given the large discrepancy seen in section  [ sec : trough ] between finite @xmath0 runtime scaling and asymptotic scaling , it makes sense to explore what values of @xmath0 are necessary to get close to the true asymptotic scaling behavior . the discrepancy seen in the last section is largest in the polynomial scaling region near the middle of the region far from the boundaries with the constant and exponential regions . much of the analysis in this section is done for @xmath53 which is a simple value in the middle of the polynomial region where this asymptotic discrepancy effect is large .      since our asymptotic analysis of @xmath19 in eqs . [ eq : gap_poly ] and [ eq : gap_exp ] is exact enough to get the constant prefactor , we can compare the numerical gap directly to the asymptotic expression to visually and numerically see where agreement occurs .     vs.  @xmath0 for @xmath53 for our three different gap extraction methods with a square barrier . the discrete gap and continuous gap well approximate each other at these relatively modes @xmath0 , but neither of them are well approximated by the asymptotic expression for the gap . this difference seems to confirm the failure of the asymptotics observed in figs . [ fig : trough_qubit ] and [ fig : trough_de ] . , scaledwidth=48.0% ]    in fig . [ fig : comparison ] we show the discrete , continuous , and asymptotic gaps for @xmath53 up to @xmath54 . since this is a log - log plot , a straight line corresponds to a power law with the slope equal to the scaling power . the continuous and discrete gaps approach each other for the reasonable values of @xmath0 shown in this plot and have very similar slopes . however , neither the continuous nor the discrete gaps are close to the asymptotic expression , nor do their slopes appear similar to the asymptotic slope . fig . [ fig : comparison ] shows how far away the actual gap is from the asymptotic expression , leading to the trough in fig . [ fig : trough_de ] .     vs.  @xmath0 for @xmath53 for our three different gap extraction methods with a square barrier . this plot covers a much larger range than fig . [ fig : comparison ] , allowing us to see that the continuous gap does eventually approach the asymptotic gap . however , this approach occurs at extremely large @xmath0 that are impractical for most current applications . in addition , we have included the second order expansion of the asymptotic large @xmath0 expression from eq . [ eq : next_order ] that more closely approximates the continuous gap . , scaledwidth=48.0% ]    with our current computational limits , we can not push exact diagonalization of the symmeterized @xmath0 qubit hamiltonian to much higher @xmath0 , but the condition in eq . [ eq : trans_con ] can be solved for extremely large @xmath0 with rounding error being the only block . therefore , we can extend fig .  [ fig : comparison ] to higher @xmath0 to get fig . [ fig : comparison_big ] from which we can see that the continuous gap does eventually approach its asymptotic value , but only for extremely large @xmath4 . [ fig : comparison_big ] also shows another curve , labeled `` second order gap . '' this curve shows the asymptotic expression for the gap in the large-@xmath0 limit , but it keeps both the leading order term and the second highest order term so that this second order gap expression is given by @xmath55 for @xmath44 where @xmath56 is the polygamma function .    the second order term does improve agreement with the continuous gap , but this fact also shows that the second order ( and higher ) terms are still relevant and important for determining run time behavior even at the very large end of what quantum computers are expected to be able to accomplish in the next few decades . since even a simple tunneling problem needs more than just its asymptotic limit to accurately describe its behavior at any reasonable @xmath0 , this data casts doubt on the efficacy of asymptotic analysis . we have exclusively used parabolic wells up to this point so the argument could be made that this large-@xmath0 behavior is an artifact of our choice of wells . to counter this argument , we have done this same analysis using an appropriately sized square well instead of a parabolic well . the asymptotic scaling behavior of @xmath19 has the same @xmath0 and @xmath49 dependence in the square well case , and extremely large @xmath0 are still needed for the asymptotics to accurately describe the exact gap for the square well .    as a general rule , we have noted that the smoother the potential is the faster @xmath19 begins to approximate its asymptotic value . we have chosen the quadratic potential since it makes the differential equation solvable , but we have run a few other simulations with different shaped barriers . while a smoother barrier ( such as a gaussian ) does level off faster , it still takes significantly large @xmath0 to have approximately the correct slope . exact analysis is more difficult for different barrier and potential shapes since we can not get exact asymptotic scaling expressions for them up to the constant prefactors . necessary for the ratio of the continuous gap and the asymptotic gap to reach the indicated percentages as a function of @xmath49 . especially close to the boundaries of the polynomial region , the system requires extremely large @xmath0 to reach its asymptotic value . when coupled with fig . [ fig : trough_de ] , this plot indicates that for the middle values of @xmath49 the continuous gap reachs its asymptotic value faster by deviating more from its asymptotic scaling . , scaledwidth=48.0% ]    we may also want to know what size of @xmath0 is necessary before the continuous gap is a certain fraction of the asymptotic gap . in fig .  [ fig : ntoreach ] , we plot the @xmath0 necessary for the ratio of the continuous gap to the asymptotic gap to reach the specified value . when we want 50% , the asymptotic gap is a good approximation relatively quickly for higher @xmath49 , but even for lower @xmath49 , extremely large @xmath0 are required . for ratios around 80% or 90% the necessary @xmath0 becomes prohibitively large , even for the best @xmath49 values . additionally , for these higher ratios , high @xmath49 also start requiring astronomical @xmath0 in addition to the low @xmath49 values . the behavior noted in fig . [ fig : ntoreach ] at first seems to contradict the information in fig . [ fig : trough_de ] since these two plots claim that opposite regions are worse . this is not a contradiction because the bad cases in fig . [ fig : trough_de ] have scaling behavior very different from the asymptotic behavior that paradoxically allows them to approach the asymptotic value of the gap faster . thus , the middle range of @xmath49 approachs the true asymptotic values of the gap faster at the cost of having vastly different scaling behavior until they reach the appropriate value . however , this faster approach to the asymptotic value is still extremely slow in comparison to the concerns of practical quantum computers , meaning that the benefit of a faster approach can not be realized .    to explore this phenomena more , let us analytically consider the ratio between the leading order asymptotic gap and the second order asymptotic gap @xmath57 where @xmath58 .    in our case @xmath59 , so the prefactor , @xmath60 , is close to @xmath61 , meaning any galactic issues with this expression are a direct result of the @xmath62 dependence . if we want this ratio to reach some threshhold value @xmath63 , then the necessary @xmath64 becomes @xmath65 since @xmath66 , the expression @xmath67 is greater than one for most possible values of @xmath63 , so the size of @xmath64 depends mostly on the exponent . when @xmath68 is close to zero , meaning the problem is close to the boundary between the constant and polynomial scaling regions , the exponent is extremely large , leading to large @xmath0 necessary to reach a specific @xmath69 . however , the expression in eq . [ eq : nofv ] would seem to indicate that the necessary @xmath0 should be relatively small for @xmath68 farther from zero as we are close to the boundary between the polynomial and exponential scaling regions . therefore the second order expansion can not explain the large spike in @xmath64 seen for high @xmath68 in fig . [ fig : ntoreach ] for large @xmath69 . given the discrepancies between the actual gap and the asymptotic gap , an important question becomes how easily scaling behavior can be mislabeled . given a data set , coming from an experiment or numerical simulation , how difficult is it to tell apart a polynomial scaling algorithm from an exponential or super - polynomial algorithm . [ fig : trough_de ] already explored how easily the asymptotic scaling power could be misidentified , but that analysis assumed a power law fit from the beginning . in this section , we explore ideas for analyzing how much a data set follows a power law model or an exponential model to judge if finite data indicates polynomial or super - polynomial scaling with @xmath0 .    whereas fig . [ fig : trough_de ] looked at a broad sweep of data in an attempt to extract scaling information , this section will be concerned with small scale features of curves , asking how well a curve can locally be approximated by a certain functional form . our ultimate goal is to be able to take data like that presented in fig .  [ fig : comparison ] and judge how likely it is that the data is obeying a power law or exponential relationship . we will look at both a power law and exponential relationship , and for ease of writing we define @xmath70 and @xmath71 : @xmath72 our analysis technique relies on properties of @xmath73 and its derivatives with repspect to @xmath11 . in fig . [ fig : trough_de ] , when we did linear fits to log - log data , we were already looking at the properties of @xmath74 . two natural other quantities to look at would be @xmath75 and the curvature of the function @xmath76 . the curvature is not a good metric since the curvatures of both @xmath77 and @xmath78 go to zero for @xmath79 . the second derivative @xmath75 is a better metric since @xmath80 goes is zero ; whereas , @xmath81 goes to @xmath82 as @xmath79 . the problem with @xmath75 is that we can not immediately extract any additional information from it about the exponential . our proposed analysis technique relies on another important feature of the derivatives of these log - log functions . specifically , we can look at the ratio , @xmath83 , of the second derivative to the derivative : @xmath84 note that only exponential scaling gives a non - zero , constant derivative ratio which gives us another metric to judge if a function has exponential scaling and allows us to discover @xmath85 . in general , sub - exponential scaling takes this ratio to zero for large @xmath0 . true polynomial scaling has a ratio of exactly zero , but given that we do not immediately reach asymptotic scaling , asymptotic polynomial and sub - exponential scaling are not immediately distinguishable through this method .     where @xmath70 and @xmath71 for the continuous gap . for true exponential scaling , this ratio should approach a constant equal to the power of @xmath0 in the exponential . the dashed lines are what this ratio should be asymptotically . none of these curves do a good job approximating their asymptotic constant value , and even the very large barrier cases require large @xmath0 to even approach the correct constant . this plot indicates that even the exponential scaling region does not approach its asymptotic expression at reasonable @xmath0 . , scaledwidth=48.0% ]    in fig . [ fig : exp_ratio ] , we have plotted the ratio @xmath86 for the continuous gap for @xmath49 values that should have asymptotically exponential scaling . the thin dashed lines represent what the curve of the same color should approach in its asymptotic limit . for lower @xmath49 , close to the boundary with the polynomial region , the @xmath83 vs.  @xmath0 curves do not come close to their asymptotic values , even for the extremely large @xmath0 shown , up to @xmath87 . some of the more strongly exponential systems do approach their asymptotic values for the @xmath0 displayed , but these systems still require extremely large @xmath0 . based on the nature of the ratio @xmath83 as a function of @xmath88 , an asymptotically increasing or non - zero constant ratio would imply a super - exponential or exponential scaling of @xmath19 with @xmath0 , so it could be argued that these plots imply non - polynomial behavior much earlier . even with this argument , many of the borderline @xmath49 cases continue decreasing for a large range of @xmath0 and do not have an upturn until much larger @xmath0 . additionally , if our goal is to extract information on the type of exponential scaling , fig . [ fig : exp_ratio ] implies that we need very large @xmath0 .     where @xmath70 and @xmath71 for the continuous gap . these @xmath49 values were chosen in the polynomial scaling region , and for true power law scaling , we would expect the ratio @xmath89 . we do see the curves approaching zero , but all of them require extremely large @xmath0 to do so , with several of the smaller barrier cases being far away from zero even at the extreme end of our simulation range . the solid , red line is for @xmath90 in the exponential region ; notice that the ratio appears to approach zero faster than many of the ratios for @xmath49s in the polynomial scaling region . , scaledwidth=48.0% ]    asymptotically polynomial scaling should give us an @xmath83 ratio that is zero , so in fig . [ fig : poly_ratio ] , we have plotted @xmath83 as a function of @xmath0 for several @xmath49 values that asymptotically have polynomial scaling . all of these curves , except the solid line for @xmath90 , should approach zero , but as we see , many of them take extremely large @xmath0 to actually approach this asymptotic value . notice especially that @xmath49 closer to @xmath91 , the boundary between constant and polynomial scaling , require even larger @xmath0 to reach zero . this behavior of small @xmath49 corresponds well with what is seen in fig . [ fig : ntoreach ] in that lower @xmath49 take a long time to reach the correct asymptotic value . also notice that the solid , red line for @xmath90 in fig . [ fig : poly_ratio ] appears to approach zero faster than many of the ratios for @xmath49 actually in the polynomial scaling region . therefore , for a large range of @xmath0 , it would be easier to interpret @xmath90 as having a power law scaling behavior than many of the @xmath49 close to @xmath91 which take very large @xmath0 to display the appropriate ratio behavior . the slowness of lower @xmath49 systems in reaching their asymptotic value is at least partly due to the size of their scaling powers . [ fig : trough_de ] indicates that these @xmath49 values close to the constant boundary reach their asymptotic scaling behavior more quickly . however , their asymptotic scaling behavior is a very small power law , meaning that any change in these functions will be small . by reaching their asymptotic scaling earlier , it will take these low @xmath49 systems even longer to reach the correct asymptotic value for their scaling , as seen in figs . [ fig : ntoreach ] and [ fig : poly_ratio ] . the problem of tunneling through a barrier in a quadratic well is a simple quantum mechanical system with elements that are not unfamiliar to an introductory quantum mechanics class , so it is surprising that asymptotics have difficulty describing such an elementary system at reasonable system sizes . asymptotic analysis is deemed useful because it is assumed that the run time of an algorithm is well approximated by its asymptotic behavior starting at relatively small system sizes , @xmath0 . our study shows that barrier tunneling is a case where asymptotic behavior is not a good descriptor of spectral gap and therefore run time scaling until extremely large @xmath0 . this failure of asymptotic analysis is worsened by the direction of the failure : asymptotically polynomial algorithms appear to locally scale in a super - polynomial manner as seen in fig . [ fig : poly_ratio ] , making it more difficult to justify the use of asymptotic analysis for small system sizes . quantum tunneling is a key part of many near - term quantum algorithms , including quantum adiabatic optimization . these algorithms are gaining popularity because of their applicability on smaller , possibly non - error - corrected , quantum computers . our results on the difficulty of describing barrier tunneling in small systems with asymptotic analysis suggest that care should be taken when analyzing the potential power of these algorithms for near - term experiments and systems .
Abstract:

Answer: we study the asymptotic behavior of the spectral gap of simple barrier tunneling problems , which are related to using quantum annealing to find the global optimum of cost functions defined over @xmath0 bits . specifically we look at the problem of having an @xmath0 qubit system tunnel through a barrier of width and height proportional to @xmath1 . we show that with these quantum annealing problems , the asymptotic , @xmath2 , behavior of the spectral gap does not accurately describe the behavior of the gap at finite @xmath0 until extremely large values of @xmath0 ( @xmath3 ) . we prove that this deficiency of the asymptotic expression is a feature of simple one - dimensional tunneling problems themselves , casting doubt on whether asymptotic analysis is an appropriate tool for studying tunneling problems in quantum annealing for reasonably sized systems .


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: detecting planetary transits is highly valuable not only because it makes an independent confirmation of a planet from radial - velocity observations but also because the photometric transits provide unambiguous information on planet masses and radii , and thus mean density and interior structure ( e.g. * ? ? ? * and references therein ) . various high - precision follow - up studies for transiting planets can also uncover planetary atmospheres and the ( mis-)alignment between the stellar spin and planetary orbital axes via the rossiter - mclaughlin effect ( e.g. * ? ? ? * ; * ? ? ? these properties of transiting planets provide valuable hints for planet formation and evolution processes including dynamical interaction between planets . although transiting planets are as such important , it is difficult to detect them around evolved stars , especially giant stars , due to the large sizes of the host stars . targets for on - going radial - velocity surveys for planets around giants have typical radii of @xmath110 @xmath6 ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) , and then the relative flux variation of such a giant host star caused by a transit of a jupiter - sized planet is only @xmath7 , which is comparable to that for transit of an earth - sized planet across a solar - type star . it is impossible to detect such a transit from ground - based photometry and thus no transiting planets had been found around giant stars . our understanding of properties of planets around such stars are thus far behind from those around solar - like stars , although tens of planets have been found around giants by precise radial - velocity surveys ( e.g. * ? ? ? * and references therein ) . the _ kepler _ mission was successful in detecting transiting planets with very high photometric precision ( @xmath8 ) since 2009 . to date it has discovered more than 4000 transiting planet candidates from sub - earth - size to super - jupiter - size in @xmath10.32000 d orbits ( e.g. * ? ? ? although most of them are orbiting solar - type stars , _ kepler _ s photometric precision is high enough to detect a jupiter - sized planet transiting a giant star . as expected , _ kepler _ has identified several planet candidates around giant stars with radii larger than @xmath9 @xcite . the planet candidates have radii comparable to or larger than jupiter s , and it is interesting that many of them have short - period orbits . since such short - period planets have rarely been found by radial - velocity surveys around evolved stars ( e.g. * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) , they will provide us unique opportunities to investigate planet formation and evolution processes around such stars .    here we report the results of radial - velocity follow - up observations using subaru 8.2 m telescope for two of the _ kepler _ evolved stars , kepler-91 ( koi-2133 ) and koi-1894 . thanks to the high precision in their radial - velocity measurements , we independently confirmed a jovian planet ( kepler-91b , koi-2133.01 ) previously reported around kepler-91 @xcite and newly found a hint for the existence of additional companion in the system . we also detected a possible sub - saturn - mass planet around koi-1894 ( koi-1894.01 ) with 23@xmath4 confidence level . the rest of the paper is organized as follows . the adopted stellar parameters for the two stars are presented in section [ targets ] . the observations are described in section [ obs ] and the results of global analysis for the radial velocities and light curves are presented in section [ analysis ] . section [ summary ] is devoted to discussion and summary . kepler-91 ( kic 8219268 , koi-2133 ; @xmath10=12.495 ) was identified as a candidate star hosting a transiting short - period jupiter - sized planet ( koi-2133.01 ; @xmath11 , @xmath12 d ) by @xcite . after that , @xcite reported the confirmation of its planetary nature based on the detailed analysis of orbital brightness modulation seen in the light curve caused by ellipsoidal variation , doppler boosting , and reflected / emitted light from planet ( e.g. * ? ? ? * ; * ? ? ? * ) , and then the planet was named kepler-91b with the radius @xmath13 and mass @xmath14 . although @xcite and @xcite claimed a possible non - planertary nature for the system , @xcite and @xcite very recently reconfirmed the planetary nature via radial - velocity measurements with a precision of @xmath1100 m s@xmath2 and @xmath120 m s@xmath2 , respectively , and obtained the planetary mass of 1.09@xmath150.20@xmath16 and 0.73@xmath150.13@xmath16 , respectively .    the stellar parameters ( effective temperature @xmath17 , surface gravity @xmath18 , mass @xmath19 , and radius @xmath20 ) of kepler-91 were reported by @xcite to be @xmath21 k , @xmath22 cgs , @xmath23 , and @xmath24 . after that , the parameters have been updated by the spectroscopic and asteroseismic analyses in @xcite and @xcite , which are consistent with each other . thus , we here adopted the values listed in @xcite : @xmath25 k , @xmath26 cgs , @xmath27 , and @xmath28 . koi-1894 ( kic 11673802 ; @xmath10=13.427 ) was also reported to be a planet - host candidate having a transiting short - period jupiter - sized planet ( koi-1894.01 ; @xmath29 , @xmath30 d ) by @xcite , though radial - velocity follow - up observations for the star have not been reported yet . the stellar parameters of koi-1894 were derived by @xcite to be @xmath31 k , @xmath32 cgs , @xmath33 , and @xmath34 , and they have been updated by the spectroscopic and asteroseismic analyses in @xcite as @xmath35 k , @xmath36 , and @xmath37 . here we adopted the updated values in this paper . recently @xcite reported a non - detection of blended stars for koi-1894 , which could have been physically associated companions and/or responsible for transit false positives if they were within @xmath10.15@xmath382.5@xmath38 separation and with magnitude difference up to @xmath39 , by high - angular - resolution ao imaging . the parameters for the two stars we adopt here are summarized in table [ tbl : stars ] . lrr [ tbl : stars ] @xmath17 ( k ) & 4550@xmath1575 & 4992@xmath1575 + @xmath18 ( cgs ) & 2.953@xmath150.007 & 2.87 + @xmath19 ( @xmath40 ) & 1.31@xmath150.10 & 1.410@xmath150.214 + @xmath20 ( @xmath41 ) & 6.30@xmath150.16 & 3.790@xmath150.190 we obtained high - precision radial - velocity data for the stars with the 8.2 m subaru telescope and the high dispersion spectrograph ( hds ; * ? ? ? * ) in 2013 and 2014 . we used the setups of stdi2b ( 2013 june 28 - 30 and 2014 july 8 , 10 , 14 - 16 ) and stdi2a ( 2013 december 6 and 2014 july 12 - 13 ) , which simultaneously cover a wavelength region of 35006200@xmath42 and 49007600@xmath42 respectively , the image slicer no.2 ( is#2 ; * ? ? ? * ) yielding a spectral resolution ( @xmath43 ) of 80000 , and an iodine absorption cell ( i@xmath44 cell ; * ? ? ? * ) for precise radial - velocity measurements . we obtained a total of 29 and 18 data points for kepler-91 and koi-1894 with typical signal - to - noise ratio of s / n@xmath453080 pix@xmath2 and 3045 pix@xmath2 , respectively , by an exposure time of 1200 sec depending on weather condition . the reduction of echelle data ( i.e. bias subtraction , flat - fielding , scattered - light subtraction , and spectrum extraction ) was performed using the iraf software package in the standard manner . we performed radial - velocity analysis for i@xmath44-superposed stellar spectra ( star+i@xmath44 ) by the method described in @xcite and @xcite , which is based on the method by @xcite and @xcite . a star+i@xmath44 spectrum is modeled as a product of a high resolution i@xmath44 and a stellar template spectrum convolved with a modeled instrumental profile ( ip ) of the spectrograph . we obtain the stellar template spectrum by deconvolving a pure stellar spectrum with the ip estimated from an i@xmath44-superposed flat spectrum . we achieved a radial - velocity precision of 517 m s@xmath2 for kepler-91 and 915 m s@xmath2 for koi-1894 . the derived radial velocities are listed in table [ tbl : rv - koi2133 ] and table [ tbl : rv - koi1894 ] together with the estimated uncertainties , and are plotted in figure [ fig : rv_koi2133_notrend ] and figure [ fig : rv_koi1894 ] . lrr [ tbl : rv - koi2133 ] 6473.00247 & 43.20 & 4.61 + 6473.01719 & 33.85 & 4.58 + 6473.03192 & 45.26 & 5.03 + 6473.96075 & 17.35 & 7.36 + 6473.97549 & 8.93 & 7.34 + 6473.99022 & 30.04 & 6.39 + 6474.95764 & @xmath4624.28 & 5.43 + 6474.97235 & @xmath4622.62 & 5.03 + 6474.98708 & @xmath4639.46 & 5.50 + 6633.69753 & @xmath4642.12 & 17.38 + 6847.93547 & 42.30 & 10.56 + 6847.95020 & 37.46 & 13.63 + 6847.96513 & 11.49 & 8.76 + 6847.98006 & 29.77 & 7.96 + 6847.99479 & 17.07 & 6.12 + 6848.00952 & 19.94 & 8.61 + 6848.08121 & 23.25 & 6.59 + 6849.75412 & @xmath4652.12 & 5.23 + 6849.76886 & @xmath4636.10 & 5.06 + 6849.78360 & @xmath4634.47 & 5.48 + 6853.04854 & @xmath4640.23 & 13.48 + 6853.97273 & 5.73 & 5.72 + 6853.98746 & 9.02 & 4.73 + 6854.00219 & 7.66 & 5.17 + 6854.84234 & 2.14 & 4.99 + 6854.85708 & 0.36 & 5.25 + 6854.87181 & @xmath466.69 & 6.03 + 6855.80523 & @xmath4643.37 & 6.63 + 6855.81828 & @xmath4641.91 & 5.28    lrr [ tbl : rv - koi1894 ] 6473.04710 & @xmath469.95 & 9.45 + 6473.06182 & @xmath460.82 & 11.01 + 6473.07656 & @xmath4611.81 & 9.94 + 6474.00523 & @xmath465.31 & 10.97 + 6474.01995 & 1.95 & 10.25 + 6474.03468 & 17.33 & 10.36 + 6474.04941 & 5.73 & 10.49 + 6475.00210 & @xmath463.24 & 10.03 + 6475.01683 & @xmath4621.07 & 9.68 + 6475.03159 & @xmath4613.24 & 8.79 + 6475.04632 & @xmath466.92 & 9.62 + 6848.09625 & @xmath464.14 & 15.02 + 6849.94269 & 0.52 & 10.34 + 6849.95743 & 13.69 & 9.86 + 6849.97216 & 10.09 & 9.04 + 6855.88939 & 0.65 & 9.63 + 6855.90411 & @xmath463.06 & 8.95 + 6855.91885 & 6.48 & 9.09 the variations in the observed radial velocities for both targets show a sign of planetary companions to kepler-91 and koi-1894 . in order to obtain accurate and precise estimates for system parameters of those systems , we here present a global analysis that makes use of all the available information from the _ kepler _ photometry and our spectroscopy . as is well known , a very precise light curve of a star orbited by planet(s ) shows a periodic modulation due to several astrophysical effects : the ellipsoidal variation , doppler boosting , and reflected / emission light from the planet ( e.g. , * ? ? ? * ; * ? ? ? the last two effects ( doppler boosting and planetary reflection / emission ) are synchronous with the planet s orbital period @xmath47 , with different peak locations along the planet phase ( @xmath48 for boosting and @xmath49 for the planetary reflection / emission , respectively ) . on the other hand , the ellipsoidal variation , which is caused by a tidal distortion by the planet s gravity , have two flux peaks ( @xmath48 and 0.75 ) within one orbit . thanks to the different phase dependence of these three effects , a very precise light curve and its modeling enable us to distinguish these three effects , and we can extract physical parameters such as the planet - to - star mass ratio @xmath50 and scaled semi - major axis @xmath51 , where @xmath52 is semi - major axis and @xmath20 is stellar radius . moreover , for transiting systems as in the present cases , incorporating the transit and/or secondary eclipse model into the above phase - curve variation lets us learn more about the planet properties ( e.g. , planet - to - star radius ratio @xmath53 and orbital inclination @xmath54 ) .    to obtain the phase - folded light curves for kepler-91 and koi-1894 , and estimate system parameters , we downloaded all available public light curves ( q0q17 ) from _ kepler _ mast archive . while only long - cadence data were available for kepler-91 , koi-1894 s light curves involve some short - cadence data . adopting the pdc - sap flux data , for which unphysical artifacts are detrended , we reduced the light curves by the following procedure . first , after removing planetary transits , we further detrended and normalized the light curve for each quarter by fitting it with a fifth - order polynomial so as to remove the long - term trends that were not removed in the pdc - sap flux . this process was repeated implementing a @xmath55 clipping to remove outliers . we then combined the light curves for all quarters and phase - folded them with the ephemerides derived by the official _ kepler _ team using the q1q16 data ; orbital period @xmath56 d and transit center @xmath57 bjd for kepler-91 , and @xmath58 d and @xmath59 bjd for koi-1894 ( from nasa exoplanet archive ) . we later consider the impact of incorrect ephemeris . finally , the folded light curves were binned into 300 and 250 phase bins for kepler-91 and koi-1894 , respectively . these bin numbers were adopted so that each bin approximately covers a cycle span of long - cadence data ( @xmath60 minutes ) . the flux error for each bin @xmath61 , which is the standard deviation of the mean flux , was computed based on the dispersion of the flux values within the bin . the long - cadence and short - cadence flux data were separately folded and binned for koi-1894 , but the binned light curve of the short - cadence data was much noisier than that for long - cadence . we thus decided to ignore the short - cadence data in the following analysis . the analysis below is based on the method by @xcite , who employed the evil - mc model @xcite for the phase - curve variation ( i.e. , ellipsoidal variations , doppler boosting , and planetary reflection ) with some revisions ( e.g. , application to an eccentric orbit ) . this phase - curve model is multiplied by the analytic transit model by @xcite , and the relevant parameters ( e.g. , @xmath51 and @xmath54 ) are simultaneously determined . in addition , we also model and fit the observed radial velocities @xmath62 . thus , the @xmath63 statistics in the present case is @xmath64 where @xmath65 and @xmath66 are the @xmath67-th observed radial velocity and its error , and @xmath68 and @xmath69 are the @xmath67-th observed light curve flux and its error , which are the mean flux and its standard deviation in the @xmath67-th bin , respectively ( e.g. , * ? ? ? * ) . the radial velocity is modeled as @xmath70 where @xmath71 , @xmath72 , @xmath73 , @xmath74 , and @xmath75 are the radial velocity semi - amplitude , true anomaly , orbital eccentricity , argument of periastron , and radial velocity offset of our dataset , respectively . for the light curve model @xmath76 , we refer the readers to @xcite , for details . in the next subsections , we describe the fitting procedure for each of the targets . the binned light curve of kepler-91 in figure [ fig : phase_koi2133 ] shows a clear pattern of phase - curve variation ; the double peaks at @xmath48 and 0.75 are representative of the ellipsoidal variation . in addition , figure [ fig : phase_koi2133 ] suggests a possible detection of secondary eclipse at @xmath77 , which means the reflected / emitted light from the planet is visible in the folded light curve . therefore , following @xcite and @xcite , we model the planet s light @xmath78 by the following expression : @xmath79 where @xmath80 is the flux variation amplitude of planet s reflected / emitted light and @xmath81 is the planetary flux offset arising from the homogeneous surface emission . following @xcite , @xcite , and @xcite , we here introduce the  phase - delay " @xmath82 of the brightest part on the planet surface from the substellar point . the planetary light @xmath78 is added to the flux model for the beaming and ellipsoidal variation , integrated over the stellar visible hemisphere and wavelength through the _ kepler _ band . in integrating the local flux calculated from the evil - mc model , we assume that the stellar effective temperature is 4550 k , and adopt the gravity darkening exponent of @xmath83 for equation ( 10 ) of @xcite based on the theoretical calculation by @xcite . the free parameters relevant to the light curve model are @xmath51 , the transit impact parameter @xmath84 , limb - darkening coefficients @xmath85 and @xmath86 for the quadratic limb - darkening law , @xmath50 ( @xmath87 ) , @xmath81 , @xmath80 , @xmath82 , the overall normalization factor @xmath88 for the folded light curve , @xmath53 , @xmath89 , @xmath90 , and @xmath91 , which represents the small time deviation of the transit center from the ephemeris reported by the _ kepler _ team . in our updated ephemeris , the initial transit center becomes @xmath92 . we fix the orbital period to be @xmath93 d that are derived by the _ kepler _ team based on the q1q16 data ( see section [ method ] ) .    assuming that the likelihood is proportional to @xmath94 in equation ( [ eq : chisq ] ) , we simultaneously model the observed radial velocity and light curve , and compute the posterior distribution for each fitting parameter by implementing markov chain monte carlo ( mcmc ) simulation . in addition to the above mentioned twelve parameters , we add @xmath75 in equation ( [ eq : rvmodel ] ) to the fitting parameters . note that the radial velocity semi - amplitude @xmath71 in equation ( [ eq : rvmodel ] ) is related to the planet - to - star mass ratio @xmath50 by @xmath95 in computing the posteriors , we do not impose priors on the fitting parameters except for @xmath96 and @xmath97 ; due to the sparse sampling of the long - cadence data ( @xmath60 minutes ) and quality of the binned light curve , the limb - darkening coefficients are poorly constrained in the absence of priors , and thus we decide to put gaussian priors on the limb - darkening coefficients based on the theoretical table by @xcite as @xmath98 and @xmath99 . in our mcmc algorithm , originally developed in @xcite , the step size of each fitting parameter is iteratively scaled so that the overall acceptance ratio falls between 15% and 35% . after running 1,000,000 chains , the best - fit value and uncertainty for each fitting parameter are estimated from the median , and 15.87 and 84.13 percentiles of the marginalized posterior distribution of that parameter . shortly after we performed the first mcmc trial , we noticed a possible trend or drift in the observed radial velocities . figure [ fig : rv_koi2133_notrend ] plots the best - fit radial velocity curve as a function of the orbital phase @xmath100 . for clarity , the data taken in 2013 are shown in green and those by the 2014 campaign are plotted in red . while we have detected a clear modulation with an amplitude of @xmath101 m s@xmath2 , radial velocities take lower values for the 2014 data , which is evident in the bottom panel indicating the residual of the observed radial velocities from the best - fit model .    in order to confirm the presence of the radial velocity drift , we try to fit the data with an additional parameter : the radial velocity drift @xmath102 . adding the drift term @xmath103 to the right side of equation ( [ eq : rvmodel ] ) , we performed again the mcmc simulation to fit both radial velocity and light curve . the derived best - fit parameters and their uncertainties are summarized in table [ tbl : best - fit - kepler91 ] . to compare between the radial velocity models with and without a drift term , we compute bayesian information criteria ( bic ) for the two cases , which is computed as @xmath104 , where @xmath105 is the number of free parameters and @xmath106 is the number of data points . in the absence of a trend , we obtain @xmath107 , and for the case of the radial velocity model with @xmath102 , bic becomes 383 ; @xmath108 is much larger than 10 , meaning that the model with @xmath102 is strongly favored . therefore , we conclude that a radial velocity drift ( or trend ) is present in our dataset , and report the best - fit parameters for the case with @xmath102 as the final result . we also tried several periods around @xmath47 and found that the @xmath109 gave slightly better results in the global fitting than @xmath47 did ( @xmath110 ) . however , the resultant parameters for the two cases are well consistent within 0.2@xmath4 level . differs by @xmath111 level . ]    figures [ fig : phase_koi2133 ] and [ fig : rv_koi2133 ] plot the folded light curve and radial velocities ( red points ) along with their best - fit models ( blue line ) . the constant radial velocity offset and trend are both removed from the radial velocity data in figure [ fig : rv_koi2133 ] . the bottom panel in each figure shows the residuals from the best - fit model . from the posterior distribution of the fitting parameters , we also estimate the orbital and planetary parameters ( e.g. , the orbital inclination @xmath54 , planet mass @xmath112 and radius @xmath113 ) assuming the stellar properties reported by asteroseismology ( table [ tbl : stars ] ) to be @xmath114 deg , @xmath115 , and @xmath116 . these estimates are also shown in table [ tbl : best - fit - kepler91 ] . lrr [ tbl : best - fit - kepler91 ] @xmath51 & 2.253@xmath117 & 2.238@xmath118 + @xmath84 & 0.9041@xmath119 & 0.9062@xmath120 + @xmath85 & 0.737@xmath121 & 0.7370@xmath122 + @xmath86 & 0.493@xmath123 & 0.493@xmath124 + @xmath125 ( 10@xmath126 ) & 4.82@xmath127 & 4.72@xmath128 + @xmath81 ( 10@xmath129 ) & 5.65@xmath130 & 5.63@xmath131 + @xmath80 ( 10@xmath129 ) & 1.41@xmath132 & 1.42@xmath133 + @xmath82 & 0.394@xmath134 & 0.395@xmath135 + @xmath53 & 0.02286@xmath136 & 0.02294@xmath137 + @xmath89 & 0.0280@xmath138 & 0.0277@xmath139 + @xmath90 & @xmath460.043@xmath140 & @xmath460.046@xmath140 + @xmath91 ( 10@xmath141 day ) & @xmath461.09@xmath142 & @xmath461.11@xmath142 + @xmath75 ( m s@xmath2 ) & @xmath4642.5@xmath143 & @xmath4633.2@xmath144 + @xmath102 ( m s@xmath2 day@xmath2 ) & @xmath460.0612@xmath145 &  + @xmath63 & 296 & 396 + bic & 383 & 477 + @xmath54 ( deg ) & 67.37@xmath146 & 67.20@xmath147 + @xmath73 & 0.0519@xmath148 & 0.0535@xmath149 + @xmath74 ( deg ) & @xmath4657.2@xmath150 & @xmath4658.8@xmath151 + @xmath112 ( @xmath16 ) & 0.66@xmath150.06 & 0.65@xmath150.06 + @xmath113 ( @xmath152 ) & 1.40@xmath150.04 & 1.41@xmath150.04      the radial - velocity variations are less visible for koi-1894 owing to the large radial - velocity error compared to the small semi - amplitude . the phase - folded light curve also shows a very tiny variation , if any , along the orbital phase . to extract possible planetary signals , we simultaneously model the radial velocities and light curve as in the case of kepler-91 . since the observed transit depth and koi-1894 s stellar radius suggest that the radius of the transiting companion ( koi-1894.01 ) is no more than @xmath153 , the reflected / emitted light from the planet is expected to be very small . according to equations ( 1 )  ( 3 ) in @xcite , the planetary reflection is estimated to be @xmath154 ppm at the location of koi-1894.01 , which is smaller than the expected amplitude of the ellipsoidal variation ( @xmath155 ppm ) . visual inspection of the binned light curve also suggests the absence of secondary eclipse . thus for simplicity , we here neglect the planet flux @xmath78 for koi-1894 , and model the folded light curve with the ellipsoidal variation ( including doppler boosting ) and transit only . neglecting @xmath78 also helps to avoid @xmath63 being stuck at a local minimum during the optimization . fixing the system parameters as @xmath156 k , @xmath157 , the gravity darkening exponent of @xmath158 @xcite , and @xmath159 derived with the q1q16 data ( see section [ method ] ) , we perform mcmc simulations for koi-1894 to infer the posterior distributions for the fitting parameters . in the present case , we have the eleven free parameters : @xmath51 , @xmath84 , @xmath85 , @xmath86 , @xmath50 , @xmath88 , @xmath53 , @xmath89 , @xmath90 , @xmath91 , and @xmath75 . again , we assume gaussian priors for the limb - darkening parameters as @xmath160 , and @xmath161 from the table by @xcite . due to the weak radial - velocity signal of koi-1894.01 , we found that the global fit to the radial velocities and light curve exhibits a degeneracy between the fitting parameters ( e.g. , @xmath162 and @xmath50 ) . thus , for koi-1894 , we also impose an additional prior on the host star s density from table [ tbl : stars ] ( @xmath163 ) , assuming a gaussian distribution @xcite . otherwise , the fitting algorithm is exactly the same as for kepler-91 . the result of the fit is summarized in table [ tbl : best - fit - koi1894 ] . the best - fit value for the mass ratio is @xmath164 , implying @xmath165 detection of the planet . the mass and radius ratios are translated as @xmath166 and @xmath167 assuming the stellar mass and radius in table [ tbl : stars ] . table [ tbl : best - fit - koi1894 ] also shows the result of our fit in the absence of the gaussian prior on @xmath168 . as expected , the mass constraint is slightly weaker for this case . the orbital eccentricity is poorly constrained due to the small planetary signal and small number of radial velocity data , but we can rule out a high eccentricity ( @xmath169 ) with our current datasets . we also try to fit the data with a circular orbit ( @xmath170 ) , and compare between the eccentric and circular cases . as shown in table [ tbl : best - fit - koi1894 ] , the circular orbital fit yields a slightly better constraint on the mass ratio ; @xmath171 with the stellar - density prior , leading to a @xmath172 detection . the mass and radius ratios are translated as @xmath173 and @xmath174 assuming the stellar mass and radius . based on the @xmath63 values for the best - fit parameters , we compare bic values for the two cases ; bic values of @xmath175 and @xmath176 are obtained for the eccentric and circular cases ( with the stellar - density prior ) , respectively , so that fitting with a circular orbit is favored . the observed radial velocities and whole light curve are plotted in figures [ fig : rv_koi1894 ] and [ fig : phase_koi1894 ] , with their best - fit models for @xmath170 . we also show zoomed - in versions of the transit and phase - curve variation in figures [ fig : transit_koi1894 ] and [ fig : ev_koi1894 ] , respectively . in figure [ fig : ev_koi1894 ] , the binned flux data used for the fit ( 250 bins ) are plotted by the red crosses , and black points with errorbars indicate the flux data binned into 30 bins . we also tested several periods around @xmath47 and found that about @xmath177 gave the minimum @xmath63 value in the global fitting ( @xmath178 ) , suggesting that the true period might exist around it . nonetheless , the resultant parameters are well consistent with those for the case of @xmath47 within 0.2@xmath4 level . differs by @xmath179 level . ] our global analysis indicates the detection of koi-1894.01 is still marginal , with only @xmath180 level , and the planet mass seems to be mainly constrained by the radial velocity data . in order to see if an independent estimate for the planet mass from the light curve alone gives a comparable result , we fit the folded light curve and estimate system parameters without radial velocity data . as a result , we obtain @xmath181 , consistent with the above result by the global fit . we reported the results of high - precision radial - velocity measurements with subaru / hds for the two _ kepler _ evolved stars kepler-91 and koi-1894 , which are candidates to host transiting planets . based on the simultaneous modeling of the radial - velocity data and _ kepler _ light curves , we independently confirmed the planetary nature of kepler-91b . the radial velocity semi - amplitude of @xmath182 m s@xmath2 we derived is consistent with what @xcite found ( @xmath183 m s@xmath2 ) with @xmath184 , but incompatible with the result by @xcite who reported @xmath185 m s@xmath2 with @xmath186 , although the reason for this discrepancy is unknown . we newly detected a drift of @xmath187 m s@xmath2 between the radial velocities taken at @xmath1 1-yr interval thanks to our better measurement precision and longer period of time of observations compared to the previous ones . the persistent radial velocity drift is suggestive of a possible presence of another companion ( planet ) outside of kepler-91b . due to the lack of data , however , we are not able to pin down the period nor mass of the companion . intensive high - precision radial - velocity monitoring of the star will uncover the unseen companion and provide hints for formation and evolution of the planetary system . the estimated parameters for kepler-91b in table [ tbl : best - fit - kepler91 ] are reasonably in good agreement with the previous study by @xcite except for the mass ratio @xmath50 and scaled semi - major axis @xmath51 . this is likely because our radial velocity data show a smaller amplitude and planet mass is estimated to be small . since the flux amplitude due to ellipsoidal variation is approximately proportional to @xmath188 @xcite , a small @xmath50 resulted in the smaller @xmath51 in order to explain the observed flux amplitude . the westward phase - shift ( @xmath189 ) of the flux maximum in planetary light ( from the substellar point ) is indicative of the inhomogeneous cloud coverage on the planetary surface . this is consistent with the trend that @xcite found , claiming that westward phase - shifts are preferentially seen for relatively cool close - in planets with equilibrium temperatures of @xmath190 k ; @xcite estimated the equilibrium temperature of kepler-91b to be @xmath191 k , depending on the assumed heat redistribution parameter . the magnitude of the phase - shift ( @xmath192 ) is , however , unusually large in our best - fit model compared with the previously reported values for other close - in planets @xcite . the reason is unknown , but the huge coverage of the planet surface illuminated by the host star of the planet surface is always illuminated by the host star due to its huge radius . ] might be responsible . note that we also detected a dip in the folded light curve ( figure [ fig : phase_koi2133 ] ) around @xmath193 reported in @xcite , and this dip should have more or less affected the fitting result . lrrrr [ tbl : best - fit - koi1894 ] @xmath51 & 3.79@xmath194 & 3.68@xmath195 & 4.25@xmath196 & 3.24@xmath197 + @xmath84 & 0.634@xmath198 & 0.56@xmath199 & 0.48@xmath200 & 0.68@xmath201 + @xmath85 & 0.724@xmath124 & 0.723@xmath202 & 0.724@xmath123 & 0.7271@xmath203 + @xmath86 & 0.391@xmath204 & 0.3908@xmath205 & 0.392@xmath206 & 0.391@xmath207 + @xmath125 ( 10@xmath126 ) & 1.19@xmath208 & 1.02@xmath209 & 1.29@xmath210 & 0.70@xmath211 + @xmath81 ( 10@xmath129 ) & 0 ( fixed ) & 0 ( fixed ) & 0 ( fixed ) & 0 ( fixed ) + @xmath80 ( 10@xmath129 ) & 0 ( fixed ) & 0 ( fixed ) & 0 ( fixed ) & 0 ( fixed ) + @xmath82 & 0 ( fixed ) & 0 ( fixed ) & 0 ( fixed ) & 0 ( fixed ) + @xmath53 & 0.01777@xmath212 & 0.01739@xmath213 & 0.01717@xmath214 & 0.01790@xmath215 + @xmath89 & 0 ( fixed ) & 0.052@xmath216 & 0 ( fixed ) & 0.047@xmath217 + @xmath90 & 0 ( fixed ) & 0.085@xmath218 & 0 ( fixed ) & 0.10@xmath219 + @xmath91 ( 10@xmath141 day ) & @xmath460.25@xmath220 & @xmath460.53@xmath221 & @xmath460.24@xmath222 & @xmath460.57@xmath223 + @xmath75 ( m s@xmath2 ) & @xmath4610.6@xmath224 & @xmath4610.0@xmath225 & @xmath4611.3@xmath226 & @xmath466.8@xmath227 + @xmath63 & 234 & 233 & 235 & 233 + bic & 285 & 295 & 286 & 294 + @xmath54 ( deg ) & 80.4@xmath228 & 80.2@xmath229 & 83.5@xmath230 & 76.2@xmath231 + @xmath73 & 0 ( fixed ) & 0.149@xmath232 & 0 ( fixed ) & 0.186@xmath233 + @xmath74 ( deg ) & 0 ( fixed ) & 61@xmath234 & 0 ( fixed ) & 65@xmath235 + @xmath112 ( @xmath16 ) & 0.18@xmath150.07 & 0.15@xmath236 & 0.19@xmath150.08 & 0.10@xmath237 + @xmath113 ( @xmath152 ) & 0.66@xmath150.04 & 0.64@xmath238 & 0.63@xmath239 & 0.66@xmath240    as for koi-1894 , we did not detect any statistically significant radial - velocity variations with our measurement precision of  915 m s@xmath2 and the number of data points . we excluded the possibility of a grazing transit by a binary companion for the star and set an upper limit on the mass of koi-1894.01 to be @xmath241 by our non - detection of radial - velocity variations with 3@xmath4 level . however , we detected possible radial - velocity variations with a semi - amplitude of @xmath115 m s@xmath2 in phase with ellipsoidal variations of the star with 23@xmath4 level . although we can not say the detection is statistically significant at this stage , it suggests that the koi-1894.01 could be a sub - saturn - mass planet . figure [ fig : plm - strad ] shows distribution of mass of exoplanets currently known plotted against their host star s radius . as seen in the figure , koi-1894.01 could be one of the lowest mass planets ever discovered around evolved stars together with kepler-56b , a super - neptune - mass planet ( @xmath242 ) detected via ttv ( transit timing variation ) method @xcite . actually the stellar parameters for koi-1891 are similar to those of kepler-56 ( * ? ? ? * @xmath243 , @xmath244 ) . confirmation of koi-1894.01 is highly encouraged in order to uncover such a new population of sub - saturn and super - neptune planets around relatively massive evolved stars , which have rarely been found so far either by radial - velocity surveys or transit ones . this research is based on data collected at subaru telescope , which is operated by national astronomical observatory of japan ( naoj ) . we are grateful to all the staff members of subaru for their support during the observations . this research has made use of the simbad database , operated at cds , strasbourg , france . is supported by japan society for promotion of science ( jsps ) fellowship for research ( pd:25 - 3183 ) . is supported by jsps research fellowships for young scientists ( no . 26 - 7182 ) and by the leading graduate course for frontiers of mathematical sciences and physics . barclay , t. et al . 2014 , , submitted ( arxiv:1408.3149 ) batalha , n.m . , et al . 2013 , , 204 , 24 butler , r. p. , marcy , g. w. , williams , e. , mccarthy , c. , dosanjh , p. , & vogt , s. s. 1996 , , 108 , 500 claret , a. 1998 , a&as , 131 , 395 claret , a. , & bloemen , s. 2011 , , 529 , a75 da silva , l. , et al . 2006 , , 458 , 609 demory , b .- o . , et al . 2013 , , 776 , l25 esteves , l.j . , de mooij , e.j.w . , & jayawardhana , r. 2013 , , 772 , 51 esteves , l. j. , de mooij , e. j. w. , & jayawardhana , r. 2014 , , submitted ( arxiv:1407.2245 ) faigler , s. , & mazeh , t. 2011 , mnras , 415 , 3921 faigler , s. , & mazeh , t. 2014 , , submitted ( arxiv:1407.2361 ) hirano , t. , narita , n. , shporer , a. , sato , b. , aoki , w. , & tamura , m. 2011 , , 63 , 531 hirano , t. , et al . 2012 , , 759 , l36 hirano , t. , et al . 2015 , , 799 , 9 huber , d. et al . 2013a , , 767 , 127 huber , d. et al . 2013b , science , 342 , 331 jackson , b. k. , lewis , n. k. , barnes , j. w. , drake deming , l. , showman , a. p. , & fortney , j. j. 2012 , , 751 , 112 johnson , j.a . , et al . 2007 , , 665 , 785 jones , m.i . , jenkins , j.s . , bluhm , p. , rojo , p. , & melo , c.h.f . 2014 , , 566 , 113 kambe , e. , et al . 2002 , , 54 , 865 law , n.m . 2014 , , 791 , 35 lillo - box , j. et al . 2014a , , 562 , 109 lillo - box , j. et al . 2014b , , 568 , 1 liu , y. , sato , b. , takeda , y. , ando , h. , & zhao , g. 2010 , , 62 , 1071 mazeh , t. , nachmani , g. , sokol , g. , faigler , s. , & zucker , s. 2012 , , 541 , a56 noguchi , k. et al . 2002 , , 54 , 855 ohta , y. , taruya , a. , & suto , y. 2009 , , 690 , 1 sato , b. , kambe , e. , takeda , y. , izumiura , h. , & ando , h. 2002 , , 54 , 873 sato , b. , et al . 2008 , , 60 , 539 sato , b. , et al . 2012 , , 64 , 97 sato , b. , et al . 2013 , , 65 , 85 seager , s. & malln - ornelas , g. 2003 , , 585 , 1038 seager , s. & deming , d. 2010 , , 48 , 631 shporer , a. , et al . 2011 , , 142 , 195 sliski , d. h. & kipping , d.m . 2014 , , 788 , 148 tajitsu , a. , aoki , w , & yamamuro , t , 2012 , , 64 , 77 takeda , y. , sato , b. , & murata , d. , 2008 , , 60 , 781 valenti , j.  a. , butler , r.  p. & marcy , g.  w. 1995 , , 107 , 966 . wang , l. , liu , y. , zhao , g. , & sato , b. 2011 , , 63 , 1035 winn , j. n. 2011 , in exoplanets , ed . s. seager ( tucson , az : university of arizona press ) , p.5577 zielinski , p. , niedzielski , a. , wolszczan , a. , admow , m. , & nowak , g. 2012 , , 547 , 91
Abstract:

Answer: we present results of radial - velocity follow - up observations for the two _ kepler _ evolved stars kepler-91 ( koi-2133 ) and koi-1894 , which had been announced as candidates to host transiting giant planets , with the subaru 8.2 m telescope and the high dispersion spectrograph ( hds ) . by global modeling of the high - precision radial - velocity data taken with subaru / hds and photometric ones taken by _ kepler _ mission taking account of orbital brightness modulations ( ellipsoidal variations , reflected / emitted light , etc . ) of the host stars , we independently confirmed that kepler-91 hosts a transiting planet with a mass of @xmath0 ( kepler-91b ) , and newly detected an offset of @xmath120 m s@xmath2 between the radial velocities taken at @xmath3-yr interval , suggesting the existence of additional companion in the system . as for koi-1894 , we detected possible phased variations in the radial velocities and light curves with 23@xmath4 confidence level which could be explained as a reflex motion and ellipsoidal variation of the star caused by the transiting sub - saturn - mass ( @xmath5 ) planet .