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: two sorts of brownian motion are defined by statistical physicists . translational brownian motion is the irregular motion exhibited by a massive particle as it collides with molecules in a fluid . collisions also affect the spin and orientation of particles ; theories of rotational brownian motion are concerned with the evolution of these quantities due to encounters . rotational brownian motion was first discussed by @xcite in the context of dielectric theory . the polarization of a dielectric material is a competition between torques due to the imposed field , which tend to align polar molecules , and collisions , which tend to destroy the alignment . the likely existence of black holes with masses much greater than those of stars has led to a renewed interest in brownian motion among astrophysicists @xcite . translational brownian motion of black holes in galactic nuclei may be directly observable ( e.g. @xcite ) , allowing constraints to be placed on black hole masses . translational brownian motion may also affect the rate of tidal disruption of stars by a nuclear black hole ( e.g. @xcite ) or the efficiency with which a binary black hole can interact with stars as it wanders through a nucleus @xcite . in paper i @xcite , the translational brownian motion of a massive binary was discussed ; it was shown that the amplitude of the wandering can be increased by a modest factor compared with that of a single point mass due to exchange of energy between field stars and the binary . here we consider the second kind of brownian motion experienced by a binary , the irregular variation of the orientation of the binary s spin axis due to encounters . a single field star that passes within a distance @xmath3 of the binary , with @xmath4 the binary s semi - major axis , will exchange orbital angular momentum with the binary , leading both to a change in the binary s orbital eccentricity as well as a change in the orientation of the binary s spin axis . repeated encounters will cause the binary s orbital eccentricity to evolve ( usually in the direction of increasing eccentricity ) and will also cause the orientation of the binary to undergo a random walk . the first process has been discussed extensively ( e.g. @xcite ) ; the second , discussed here for the first time , will be called rotational brownian motion , by analogy with the similar process that occurs in a fluid of polar molecules ( e.g. @xcite ) .    the orbital orientation of a black - hole binary influences the direction of the spin axis of the single black hole that forms after coalescence of the binary due to emission of gravitational radiation ( e.g. @xcite ) . the spin axis , in turn , is thought to determine the plane of the inner accretion disk which forms around the black hole @xcite , the direction of jets launched from the accretion disk @xcite , etc . thus changes in the orientation of a black - hole binary as it coalesces may have observable consequences .    in  2 the conservation equations for interaction of a field star with a massive binary are presented ; most of these have been seen before and were used in earlier studies to derive the rates of change of the binary s separation and eccentricity . the rotational diffusion equation is presented in  3 and the order - of - magnitude of the rotational diffusion coefficient @xmath5 is derived .  4 describes the numerical derivation of @xmath5 from scattering experiments . @xmath6 presents solutions of the diffusion equation and derives a relation between the degree of hardening of a binary and the change in its orientation . consider an encounter of a single particle ( the `` field star '' ) with a binary . assume that the binary remains bound during the encounter , appropriate if the mass of the field star is much less than the mass of each of the components of the binary or if the binary is sufficiently hard . at early times , the field star has velocity @xmath7 with respect to the center of mass of the field star - binary system and its impact parameter is @xmath8 . a long time after the encounter , the velocity @xmath9 of the field star attains a constant value . conservation of linear momentum implies that the change @xmath10 in the velocity of the binary s center of mass is given by = -12 where @xmath11 is the mass of the field star and @xmath12 is the mass of the binary . this relation was used in paper i in conjunction with scattering experiments to investigate the translational brownian motion of a binary due to encounters . the energy of the field star - binary system , expressed in terms of pre - encounter quantities , is    @xmath13    with @xmath14 the initial velocity of the binary s center of mass and @xmath15 the binary s initial semi - major axis . after the encounter , e= 12(1 + 12)v^2 - gm_1m_22a and @xmath16 , so that    @xmath17    where the latter relation assumes @xmath18 . this result was used by @xcite , @xcite , and @xcite to compute the hardening rate , @xmath19 , of a massive binary due to encounters . the angular momentum of the field star - binary system about its center of mass , expressed in terms of pre - encounter quantities , is _ 0=(1 + 12)**_0 + _ 12**_b0 where @xmath20 and @xmath21 with @xmath22 the binary s spin angular momentum . conservation of angular momentum during the encounter gives    @xmath23    changes in @xmath24 correspond to changes in the binary s orbital eccentricity @xmath25 via the relation @xmath26 @xcite . changes in the direction of @xmath27 correspond to changes in the orientation of the binary , leading to rotational diffusion . let @xmath28 be the probability that the spin axis of the binary is oriented within the solid angle @xmath29 at time @xmath30 . we seek an equation describing how @xmath31 evolves with time due to encounters of field stars with the binary , if each encounter is assumed to produce a negligibly small change in @xmath32 and @xmath33 . if we imagine that @xmath22 is initially oriented parallel to the @xmath34 axis , and that the encounters leading to changes in the orientation of the binary are isotropic in velocity and direction , then @xmath31 will evolve in such a way as to remain a function of @xmath32 alone , @xmath35 . these are the same assumptions made by @xcite in his theory of the rotational brownian motion of spherical molecules . debye showed that the evolution equation for @xmath31 , in the absence of an external torque , is = 1 ( ) [ eq_evol1 ] with @xmath5 the rotational diffusion coefficient , = ( d , d)^2d ; @xmath36 is the probability that , during a unit interval of time , a binary whose angular momentum @xmath22 is directed toward @xmath29 will reorient itself such that @xmath22 lies within @xmath37 , and @xmath38 is the angular separation between @xmath29 and @xmath37 .    from its definition , @xmath5 is the sum , over a unit interval of time , of @xmath39 due to encounters with field stars . it can be expressed in terms of mean changes in the components of the binary s angular momentum as follows . let the binary be initially oriented with its spin vector along the @xmath40 axis , @xmath41 . after an encounter , the angle @xmath38 between initial and final spins is = _ b , z_b = _ b , z . [ eq_cost ] expanding equation ( [ eq_cost ] ) to second order in @xmath38 and taking averages , we find = _ b , x^2 + _ b , y^2_b0 ^ 2 which can further be expressed in terms of changes in the angular momentum components of the field star via equation ( [ eq_dl ] ) as    @xmath42 . \label{eq_diff1}\end{aligned}\ ] ]    the last relations assume @xmath18 . specializing to the case of an equal - mass , circular - orbit binary , equation ( [ eq_diff1 ] ) becomes = 16(12)^2 _ x^2 + _ y^2g12a . [ eq_diff2 ] equation ( [ eq_diff2 ] ) is a prescription for computing the rotational diffusion coefficient from scattering experiments . this prescription will be applied in the next section . we can estimate the order of magnitude of the rotational diffusion coefficient by noting that @xmath43 , the relative velocity of @xmath44 and @xmath45 , and that a single close encounter produces a change in @xmath46 of the field star of order @xmath47 . the rate of close encounters is @xmath48 with @xmath49 and @xmath50 the number density and 1d velocity dispersion of the field stars ( paper i ) . thus @xmath51 and 32 where @xmath52 . the same encounters that induce rotational diffusion will also cause the binary to harden , according to equation ( [ eq_harden ] ) . the hardening rate may be written ( 1a ) = gh [ eq_defh ] with @xmath53 for a hard binary ( e.g. @xcite ) . thus [ eq_compare ] with @xmath54 . equation ( [ eq_compare ] ) suggests that rotational diffusion occurs on a time scale that is of order @xmath55 times the hardening time , or that the binary will rotate by an angle of order @xmath56 in the time it takes to harden significantly . this result will be assessed more quantitatively in  5 . the same set of scattering experiments described in paper i were used here to evaluate @xmath57 and @xmath58 . as described more fully in that paper , field stars were treated as massless particles moving in the potential of the binary whose center of mass and orbital parameters remained fixed . orbits were integrated using the routine d0p853 of @xcite , an 8(6)th order embedded runge - kutta integrator ; each field star was assumed to begin at a position @xmath59 and was advanced from @xmath60 to @xmath61 along a keplerian orbit about a point mass @xmath62 . the integrations were terminated when the star had moved a distance from the binary that was at least @xmath63 times its initial distance with positive energy . the binary s orbit was assumed to be circular ; the orientation of the binary s orbital plane with respect to the initial velocity vector of the field star , and the initial phase of the binary , were chosen randomly for each integration . unless otherwise indicated , distances and velocities are given below in program units of @xmath4 and @xmath64 respectively . the orbit integrations , and the expressions given above relating changes in the field - star s angular momentum to the orientation of the binary , were checked by carrying out a limited set of experiments using a fully general , three - body integrator and various values for @xmath65 . the expressions derived above for @xmath66 etc . were found to be accurately reproduced in the limit of small @xmath65 . consider first the rotational diffusion coefficient corresponding to interactions with stars of a single velocity @xmath67 at infinity . multiplying @xmath39 in a single encounter by @xmath68 , with @xmath8 the impact parameter and @xmath49 the density of field stars , and integrating over @xmath8 gives    @xmath69    ( cf . equation [ eq_diff2 ] ) , where @xmath70 are defined as mean square changes at fixed @xmath71 , averaged over many encounters in which the orientation of @xmath72 with respect to the binary s spin axis is taken to be random . figure 1 shows the distribution of field - star velocity changes as a function of @xmath73 for @xmath74 . we define a dimensionless , velocity - dependent diffusion coefficient @xmath75 as    @xmath76    with @xmath77 . figure 2 plots @xmath75 as a function of @xmath67 and @xmath78 . the most effective encounters are those with small @xmath67 and @xmath78 ; encounters with impact parameters @xmath79 have a closest - approach distance to the binary of @xmath80 and produce almost no change in the binary s orientation . hence @xmath75 reaches a maximum value at @xmath81 before levelling off at larger @xmath8 .    multiplying @xmath82 by @xmath83 , the distribution of field - star velocities , and integrating @xmath84 gives @xmath5 . we assume a maxwellian distribution of field star velocities , ( v_f ) = 1(2 ^ 2)^3/2 e^-v_f^2/2_f^2 , [ max ] so that = 4(2 ^ 2)^3/2_0^ dv  v^2 e^-v^2/2 ^ 2 ^ 2 . [ eq_defl1 ] this can be written in terms of a dimensionless coefficient @xmath85 ,    @xmath86    @xmath85 is defined in analogy with the dimensionless coefficients @xmath87 that describe respectively the rates of change of the binary s energy , the rate of mass ejection by the binary , and the binary s eccentricity growth rate ( cf . @xcite ) . note that @xmath85 is a function of two parameters , @xmath88 , where r = p_max^2g12 ,  s. @xmath89 is the maximum impact parameter in units of the radius of gravitational influence of the binary and @xmath90 is the dimensionless hardness of the binary . figure 3 plots @xmath91 . at a given hardness @xmath90 , @xmath85 reaches its limiting value @xmath92 by @xmath93 . for all but the softest binaries , this limiting value is reached before @xmath94 , which is roughly the value of @xmath89 we expect for a massive binary at the center of a steeply - falling density cusp ( paper i ) . hence we ignore the dependence of @xmath85 on @xmath78 in what follows and write @xmath95 . @xmath96 varies from @xmath97 for the softest binaries ( @xmath98 ) to @xmath99 for @xmath100 . equation ( [ eq_evol1 ] ) describes the evolution of @xmath101 , the probability distribution for the binary s orientation , given the rotational diffusion coefficient @xmath5 . the latter ( equation 20 ) is a function of the binary separation @xmath4 , directly via the factor @xmath102 , and indirectly via the weak dependence of @xmath85 on binary hardness ( figure 3 ) . furthermore the timescale for the binary to reorient itself is shorter than the hardening time @xmath103 ( cf . equation [ eq_compare ] ) , hence the change in @xmath4 can not be neglected when solving for the evolution of the orientation . we nevertheless begin by considering an idealized model in which @xmath5 , @xmath4 , and the parameters @xmath104 that describe the background are assumed constant in time . the evolution equation ( [ eq_evol1 ] ) becomes = 12 , where @xmath105 , @xmath106 , and t_0 2(12)ga l. consider the evolution of @xmath107 , the expectation value of @xmath108 . we have      or ( ) = _ 0 e^- . hence @xmath110 is the time constant for relaxation to a uniform distribution of orientations , @xmath111 , given a fixed binary separation and unchanging stellar background . it may be expressed in physical units as t_0 = 6.110 ^ 9 yr ( 1210 ^ 6 ) ( l50)^-1 ( 10 ^ 6 ^ -3)^-1 ( 200  km  s^-1 ) ( a1  pc)^-1 .    next consider the more interesting case where @xmath4 , as well as the parameters @xmath112 and @xmath50 that describe the field stars , may be changing with time . this time dependence is not easily specified but we can make progress by changing evolution variables from @xmath30 to @xmath113 . first combining equations ( [ eq_evol1 ] ) and ( [ eq_defdv ] ) , we find = l412ga . now changing variables , = fxdxdt = fx 1adadt = -fx gha ( cf . equation ( [ eq_defh ] ) ) and so = -14 lh 12 . this equation gives the evolution of the binary s orientation in terms of changes in its semi - major axis @xmath4 , with no explicit dependence on the parameters that describe the stellar background . ignoring the weak dependence of @xmath114 and @xmath85 on @xmath4 , appropriate for a hard binary , we find as before ( a ) = _ 0(aa_0)^12lh12 . the exponent in this expression is of order @xmath115 . hence we can write ( a)1 + 12lh12(aa_0 ) where @xmath116 has been set to unity , corresponding to an initial orientation parallel to the @xmath34 axis .    for the case of binary supermassive black holes , if we define @xmath15 as the separation when the binary first forms a bound pair , we expect gravitational radiation coalescence to occur when @xmath117 @xcite . since @xmath118 for a hard binary , we find for the expectation value of @xmath108 at coalescence : 1 - 1012 . [ eq_approx ] writing @xmath119 , the rms change in the angle defined by the binary s spin axis , this becomes ,  12 . for @xmath120 , e.g. a binary of mass @xmath121 surrounded by @xmath122 stars , this predicts @xmath123 . for a binary of so - called `` intermediate - mass '' black holes with @xmath124 , @xmath125 . milosavljevic & merritt ( 2001 ) give information about the rotational brownian motion of a massive binary in a set of @xmath126-body simulations . they show ( in their figure 11 ) the evolution of the orientation of an equal - mass , circular - orbit binary in three simulations with @xmath127 . the simulations all follow the evolution of the binary until @xmath128 ; hence equation ( [ eq_approx ] ) predicts @xmath129 . these values are in excellent agreement with the angular deflections shown in figure 11 of milosavljevic & merritt ( 2001 ) .    under what circumstances could the roeorientation angle of a supermassive black hole binary be much greater than the modest values predicted here ? since @xmath130 , one possibility is for the perturbing objects to be much more massive than @xmath122 ; for instance , they could be supergiant stars ( @xmath131 ) , star clusters ( @xmath132 ) , or giant molecular clouds ( @xmath133 ) . such physically large objects would presumably produce most of their perturbing torques when the binary separation was greater than or comparable to their own size . an even more extreme possibility is that the perturbing objects are themselves supermassive black holes . consider for example a binary supermassive black hole at the center of a giant elliptical galaxy . the stellar densities at the centers of such galaxies are low and one might expect the coalescence of the binary to stall due to a shortage of interacting stars ( e.g. @xcite ) . in this case , further decay of the binary might await the infall of a third and subsequent black holes from later galactic mergers ( e.g. @xcite ) . mass ratios @xmath65 would then be of order @xmath134 and signficant reorientation of the binary could occur before coalescence . orientations of radio jets in seyfert galaxies are known to be almost random with respect to the plane of the stellar disk ( e.g. @xcite ) . if jets are launched parallel to the spin axes of rotating black holes , this implies that the spins of nuclear black holes bear almost no relation to disk spins . while the modest reorientation rates derived here are probably insufficient to explain this phenomenon , a related mechanism may suffice . the spin orientation of a supermassive black hole that forms via _ repeated _ mergers with smaller black holes would undergo a random walk due to the random orbital inclinations of the infalling black holes . if nuclear black holes in seyfert galaxies grew through multiple mergers before the formation of the stellar disks , one might expect their spin axes to be uncorrelated with those of the stars and gas . backer , d. c. & sramek , r. a. 1999 , , 520 , 137 bardeen , j. m. , & petterson , j. a. 1975 , , 195 , l65 debye , p. 1913 , ber . phys . ges . 15 , 777 debye , p. 1929 , polar molecules ( new york : the chemical catalog company ) dorband , e. n. , hemsendorf , m. & merritt , d. 2001 , astro - ph/0112092 flanagan , e. e. & hughes , s. a. 1998 , phys . d. , 57 , 4535 hairer , e. , nrsett , s. p. & wanner , g. 1991 , solving ordinary differential equations i. nonstiff problems ( berlin : springer ) hills , j. g. 1983 , , 88 , 1269 hills , j. g. 1992 , , 103 , 1955 kinney , a. l. , schmitt , h. r. , clarke , c. j. , pringle , j. e. , ulvestad , j. s. & antonucci , r. r. j. 2000 , , 537 , 152 mcconnell , j. 1980 , rotational brownian motion and dielectric theory ( london : academic press ) merritt , d. 2000 , in dynamics of galaxies : from the early universe to the present , eds . f. combes , g. a. mamon , & v. charmandaris ( asp conference series , vol . 197 ) , p. 221 . merritt , d. 2001 , , 556 , 245 ( paper i ) milosavljevic , m. & merritt , d. 2001 , , 562 , 000 ( astro - ph/0103350 ) mikkola , s. & valtonen , m. j. 1992 , , 259 , 115 quinlan , g. d. , 1996 , newa , 1 , 35 quinlan , g. d. & hernquist , l. 1997 , new a , 2 , 533 rees , m. j. 1978 , nature , 275 , 516 reid , m. j. , readhead , a. c. s. , vermeulen , r. c. & treuhaft , r. n. 1999 , , 524 , 81 valtonen , m. j. 1996 , comm . , 18 , 191 young , p. j. 1977 , , 215 , 36
Abstract:

Answer: the orientation of a massive binary undergoes a random walk due to gravitational encounters with field stars . the rotational diffusion coefficient for a circular - orbit binary is derived via scattering experiments . the binary is shown to reorient itself by an angle of order @xmath0 during the time that its semi - major axis shrinks appreciably , where @xmath1 is the binary mass and @xmath2 the perturber mass . implications for the orientations of rotating black holes are discussed .


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: one of the most important constraints on supersymmetric models comes from considerations of the mass of the lightest neutral higgs , @xmath0 . in the minimal supersymmetric standard model ( mssm ) , the mass of the @xmath0 field at tree level @xcite must be less than the @xmath5 mass , but radiative corrections ( primarily due to top quark / squark loops ) can lift this to @xmath6 gev @xcite . while non - minimal supersymmetric models can evade this bound somewhat , it is very difficult @xcite to raise the upper bound by more than a few tens of gev . another constraint comes from the mass of the charged higgs boson , @xmath1 . at tree level , one has @xcite @xmath7 where @xmath3 is the pseudoscalar higgs mass . thus , the charged higgs mass can not be much heavier than the pseudoscalar , and the decay @xmath4 is kinematically forbidden@xcite . similar considerations involving the mass of the heavier neutral scalar , @xmath2 , show that at tree level @xmath8 which also shows that the charged higgs can not decay into @xmath9 . radiative corrections to these formulae have been calculated @xcite and , as in the neutral case , are typically tens of gev .    in this brief report , we note that one of the simplest extensions of the mssm substantially changes these bounds ; the addition of a sequential fourth generation . interest in a fourth chiral generation ( with a sufficiently heavy neutrino to avoid contributing to the z width ) has fluctuated over the years @xcite . in the early part of this decade , it was believed by many that electroweak radiative corrections ruled out a chiral fourth generation @xcite . however , this belief was based on the assumption that the quark masses were degenerate . it has been shown @xcite that relaxing this assumption allows a fourth generation , and for a mass splitting of @xmath10 , the s and t corrections fit inside the one - sigma error ellipse @xcite .    since corrections to the above mass relations vary as the fourth power of the top quark mass , it is clear that they could be substantially enhanced by a fourth generation . since the mass splitting between the @xmath11 and @xmath12 quarks is much smaller than their masses , we will assume that they are degenerate in our analysis . we will consider masses between @xmath13 and @xmath14 gev . the lower bound is somewhat lower than current tevatron limits @xcite , since it might be possible to weaken these limits if supersymmetric decay modes are allowed . the upper bound comes from the requirement that perturbation theory be valid . generically , the requirement that the yukawa couplings , @xmath15 in the mssm satisfy the condition @xmath16 implies @xcite that @xmath17 , where @xmath18 , @xmath19 is the ratio of vacuum expectation values , @xmath20 is the fermion mass , and @xmath21 gev . this window of @xmath19 closes for masses greater than @xmath22 gev . since the yukawa coupling increases with scale , and we want the model to be valid up to at least a tev , we will set the upper bound to @xmath14 gev .    in the next section , we look at the lightest neutral scalar , and then consider the mass relations involving the charged higgs scalar in section iii . section iv contains our conclusions . as noted above , @xmath19 must be fairly close to one in order for perturbation theory be valid . in this case , the tree level value of the lightest neutral higgs mass is negligible , and the entire mass must come from radiative corrections . in the three generation case , these corrections have been calculated , first in ref . @xcite , and then with increasing precision by many others @xcite . generally , the contributions from radiative corrections will not be sufficient to increase the mass above @xmath23 gev ( if the tree level mass is negligible ) .    in the case of four generations , this is not the case . some early works exploring four generation supersymmetry models@xcite did not really address the possibility of fourth generation quarks above 300 gev . more recently , an exploration of the electroweak phase transition in the four generation mssm appeared@xcite , but it used a value of @xmath24 which put the @xmath11 yukawa coupling far outside the region of perturbative validity . last year , murdock , et al.@xcite discussed ways of extending the validity of perturbation theory in the four generation mssm . the only recent discussion of the lightest higgs mass in the four generation mssm was the work of fok and kribs @xcite , who plotted the lightest higgs mass as a function of the heavy quark mass . they did restrict the model to the case in which the scalar @xmath12 and @xmath11 masses were within twenty percent of the @xmath12 and @xmath11 masses , since that was necessary to obtain sufficient electroweak baryogenesis .    the expression used by fok and kribs @xcite in the limit in which mixing between scalar quarks is ignored , all squarks have the same mass , the heavier higgs bosons are all very heavy and @xmath25 is @xmath26 where @xmath21 gev . a complete expression which includes two - loop leading logs , as well as squark mixing , can be found in the review article of carena and haber @xcite . we have used that complete expression , and have plotted the results in figure 1 for various values of the susy scale . the squark masses are presumed to be equal and with a value of @xmath27 the solid ( dashed ) lines are for a pseudoscalar mass of @xmath28 ( @xmath22 ) gev , and the curves , from the lower curve up , correspond to @xmath29 gev , respectively ]    in the graph shown , we have taken the mass of the pseudoscalar to be either @xmath22 or @xmath28 gev , and have fixed @xmath25 . the results are very insensitive to squark mixing and the precise values of @xmath19 and @xmath30 . however , they are sensitive to the pseudoscalar mass . for a pseudoscalar mass of several tev , the results agree with that given in eq . ( 3 ) , but for pseudoscalar masses of @xmath28 ( @xmath22 ) gev , the results ( for @xmath31 gev ) are smaller than that of eq . ( 3 ) by @xmath32 ( @xmath33 ) percent . we see that a large mass for the lightest higgs boson in the four generation mssm not only possible , but is likely for susy scales above several hundred gev . even a mass of close to @xmath14 gev can not be excluded . as was noted earlier , eq . ( 1 ) prevents the decay @xmath34 at tree level . this makes detection of the charged higgs quite difficult , since the @xmath35 decay has very large backgrounds . radiative corrections to eq . ( 1 ) can be found in ref . @xcite , and when these are modified to account for a fourth generation , the decay may no longer be kinematically forbidden . one obtains an expression@xcite for @xmath36 which depends primarily on @xmath37 and @xmath38 . the @xmath19 dependence is very small due to the above requirement of the validity of perturbation theory . the dependence on @xmath37 is illustrated in figure 2 for several values of @xmath38 and @xmath19 fixed to 1 . we see that the @xmath34 decay becomes kinematically allowed for much of parameter space . the squark masses are presumed to be equal and with a value of @xmath27 the curves , from the lower curve up , correspond to @xmath29 gev , respectively . @xmath39 is the average of the charged scalar and pseudoscalar masses . ] the situation is different with respect to the decay @xmath40 . radiative corrections to eq . ( 2 ) are given in ref . @xcite , and these can easily be generalized to a fourth generation . we find that @xmath41 is never less than @xmath42 , and in fact usually @xmath43 by tens of gev . therefore this decay process is still forbidden , even with higher order corrections . in addition to the possibility of the decay process @xmath44 , one has the conventional processes @xmath45 and , for sufficiently heavy charged higgs bosons , @xmath46 furthermore , one can also have @xmath47 . the latter decay will be small in this case . the reason is that in the limit of @xmath48 gev , one has@xcite @xmath49 , and thus the coupling of the light higgs to the charged higgs and a @xmath50 , which is proportional to @xmath51 , vanishes as @xmath52 . of course , one might also have various supersymmetric decay modes . we will only consider decays into @xmath53 , @xmath54 and @xmath55 .    , @xmath35 and @xmath46 versus @xmath37 , choosing @xmath56 gev and @xmath57 gev . ]    in figure 3 , we plot the branching ratios of @xmath4 , @xmath58 and @xmath59 versus @xmath37 , choosing @xmath56 gev and @xmath60 gev . note that if the pseudoscalar mass is @xmath28 gev , then one can see from figure 2 that @xmath4 is only allowed for very large values of @xmath37 , and the branching ratio for this mode is always less than a percent . for @xmath61 gev , the @xmath59 decay rapidly becomes kinematically inaccessible , and the branching ratio for @xmath4 can be as large as @xmath62 . note that in this case , the pseudoscalar will primarily decay into @xmath63 , and observation of the associated @xmath50 may be the only way to detect it@xcite . in the mssm , the higgs sector is very tightly restricted . there is an upper bound of approximately @xmath64 gev on the lightest higgs mass and the charged higgs mass can not differ by more than tens of gev from the pseudoscalar and heavier neutral higgs masses , thus ruling out two - body decays into these states . in this brief report , motivated by recent interest in four generation models , we have considered the effects of radiative correction in the four generation mssm . the upper bound on the lightest higgs mass is substantially increased , and can be as large as @xmath14 gev for reasonable parameter values . the mass splitting between the charged higgs and pseudoscalar is increased to the point where the @xmath44 decay is kinematically allowed , whereas the splitting between the charged higgs and heavy neutral scalar remains small . should a fourth generation exist , it will likely be discovered shortly at the lhc . our main point is that radiative corrections from these quarks can have major implications for the higgs sector , not only for the mssm , but other extensions of the standard model . we are grateful to chris carone , david reiss and graham c. kribs for useful discussions . this work was supported by the national science foundation grant nsf - phy-0757481 and the reu grant nsf - phy-0755262 . 99 k.  inoue , a.  kakuto , h.  komatsu and s.  takeshita , prog . phys .   * 67 * , 1889 ( 1982 ) ; r.  a.  flores and m.  sher , annals phys . * 148 * , 95 ( 1983 ) . h.  e.  haber and g.  l.  kane , phys . rept .   * 117 * , 75 ( 1985 ) . j.  f.  gunion , h.  e.  haber , g.  l.  kane and s.  dawson , _ the higgs hunter s guide _ ( addison - wesley , reading 1990 ) . h.  e.  haber and r.  hempfling , phys . lett  * 66 * , 1815 ( 1991 ) ; j.  r.  ellis , g.  ridolfi and f.  zwirner , phys . b * 257 * , 83 ( 1991 ) ; y.  okada , m.  yamaguchi and t.  yanagida , prog . theor .  phys .   * 85 * , 1 ( 1991 ) . g.  l.  kane , c.  f.  kolda and j.  d.  wells , phys . * 70 * , 2686 ( 1993 ) [ arxiv : hep - ph/9210242 ] . for a discussion of the symmetries giving these relations , see j.  m.  gerard and m.  herquet , phys . lett .   * 98 * , 251802 ( 2007 ) [ arxiv : hep - ph/0703051 ] . m.  a.  diaz and h.  e.  haber , phys .  rev . d * 45 * , 4246 ( 1992 ) . for a review and an extensive set of references of work before 2000 , see p.  h.  frampton , p.  q.  hung and m.  sher , phys . rept .   * 330 * , 263 ( 2000 ) [ arxiv : hep - ph/9903387 ] . for a review and extensive set of references of work after 2000 , see b.  holdom , w.  s.  hou , t.  hurth , m.  l.  mangano , s.  sultansoy and g.  unel , arxiv:0904.4698 [ hep - ph ] . c.  amsler _ et al . _ [ particle data group ] , phys . b * 667 * , 1 ( 2008 ) . m.  maltoni , v.  a.  novikov , l.  b.  okun , a.  n.  rozanov and m.  i.  vysotsky , phys .  lett . b * 476 * , 107 ( 2000 ) [ arxiv : hep - ph/9911535 ] ; h.  j.  he , n.  polonsky and s.  su , phys .  rev . d * 64 * , 053004 ( 2001 ) [ arxiv : hep - ph/0102144 ] ; b.  holdom , phys . d * 54 * , 721 ( 1996 ) [ arxiv : hep - ph/9602248 ] ; g.  d.  kribs , t.  plehn , m.  spannowsky and t.  m.  p.  tait , phys . d * 76 * , 075016 ( 2007 ) [ arxiv:0706.3718 [ hep - ph ] ] ; j.  alwall _ et al . _ , eur . j.   c * 49 * , 791 ( 2007 ) [ arxiv : hep - ph/0607115 ] . www-cdf.fnal.gov/physics/exotic/r2a/20090327.bprime/ ; www-cdf.fnal.gov/physics/new/top/2008/tprop/tprime 2.8/public.html e.  de pree , g.  marshall and m.  sher , arxiv:0906.4500 [ hep - ph ] . r.  barbieri and m.  frigeni , phys . b * 258 * , 395 ( 1991 ) ; m.  drees and m.  m.  nojiri , phys .  rev . d * 45 * , 2482 ( 1992 ) ; j.  a.  casas , j.  r.  espinosa , m.  quiros and a.  riotto , nucl . b * 436 * , 3 ( 1995 ) [ erratum - ibid . b * 439 * , 466 ( 1995 ) ] [ arxiv : hep - ph/9407389 ] ; j.  r.  espinosa and r.  j.  zhang , jhep * 0003 * , 026 ( 2000 ) [ arxiv : hep - ph/9912236 ] ; j.  r.  espinosa and r.  j.  zhang , nucl . b * 586 * , 3 ( 2000 ) [ arxiv : hep - ph/0003246 ] ; a.  brignole , g.  degrassi , p.  slavich and f.  zwirner , nucl . b * 631 * , 195 ( 2002 ) [ arxiv : hep - ph/0112177 ] ; d.  m.  pierce , j.  a.  bagger , k.  t.  matchev and r.  j.  zhang , nucl .  phys . b * 491 * , 3 ( 1997 ) [ arxiv : hep - ph/9606211 ] ; h.  e.  haber , r.  hempfling and a.  h.  hoang , z.  phys . c * 75 * , 539 ( 1997 ) [ arxiv : hep - ph/9609331];s . heinemeyer , w.  hollik and g.  weiglein , phys . b * 440 * , 296 ( 1998 ) [ arxiv : hep - ph/9807423 ] ; s.  heinemeyer , w.  hollik and g.  weiglein , eur . j.   c * 9 * , 343 ( 1999 ) [ arxiv : hep - ph/9812472 ] . j.  f.  gunion , d.  w.  mckay and h.  pois , phys . d * 53 * , 1616 ( 1996 ) [ arxiv : hep - ph/9507323 ] ; m.  s.  carena , h.  e.  haber and c.  e.  m.  wagner , nucl . b * 472 * , 55 ( 1996 ) [ arxiv : hep - ph/9512446 ] . s.  w.  ham , s.  k.  oh and d.  son , phys . d * 71 * , 015001 ( 2005 ) [ arxiv : hep - ph/0411012 ] . z.  murdock , s.  nandi and z.  tavartkiladze , phys . b * 668 * , 303 ( 2008 ) [ arxiv:0806.2064 [ hep - ph ] ] . r.  fok and g.  d.  kribs , phys . d * 78 * , 075023 ( 2008 ) [ arxiv:0803.4207 [ hep - ph ] ] . m.  s.  carena and h.  e.  haber , prog . phys .   * 50 * , 63 ( 2003 ) [ arxiv : hep - ph/0208209 ] . a.  g.  akeroyd and s.  baek , phys . b * 525 * , 315 ( 2002 ) [ arxiv : hep - ph/0105228 ] . for a discussion of the phenomenology at the lhc , see s.  de visscher , j.  m.  gerard , m.  herquet , v.  lemaitre and f.  maltoni , arxiv:0904.0705 [ hep - ph ] .
Abstract:

Answer: in the minimal supersymmetric standard model ( mssm ) with three generations of fermions , there is a stringent upper bound on the mass of the lightest neutral higgs @xmath0 , and the mass of the charged higgs @xmath1 , must be close ( within tens of gev ) to the heavier neutral higgs , @xmath2 , and the pseudoscalar higgs @xmath3 . in this brief report , we show that in the four generation mssm , the upper bound on the @xmath0 mass is much higher , as high as 400 gev , and @xmath1 is generally much heavier than the @xmath3 , allowing the @xmath4 decay , potentially changing search strategies for both the charged higgs and the pseudoscalar . the @xmath2 mass , on the other hand , remains within tens of gev of the charged higgs mass .


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 standard model ( sm ) has been very successful in explaining the data taken from former colliders such as cern lep or fermilab tevatron . testing sm at the cern lhc will either lead to additional confirmation of the sm or give some hints for new physics beyond the sm . because of the large mass of the top quark , its couplings are expected to be more sensitive to new physics than other particles @xcite . especially , wtb vertex deserves special attention since the top quark is expected to decay almost completely via this interaction . thus , studying top quark couplings will be substantial to test the sm and a deviation of the top couplings from the expected values would imply the existence of new physics .    in this work we have analyzed anomalous wtb couplings via single top quark production in @xmath2 collision at the lhc . this reaction is probable at the lhc via elastic photon emission from one of the incoming protons . the emitted photon can collide with the other proton and produce a final state of @xmath3 through deep inelastic scattering ( fig.[fig1 ] ) . we employ the equivalent photon approximation ( epa ) @xcite for elastic photon emission from the proton . in the epa , emitted photons have a low virtuality and it is a good approximation to assume that they are on - mass - shell . for this reason these photons are sometimes called quasi - real photons . when a proton emits a quasi - real photon it remains intact and scatters with a very small angle from the beam pipe . the atlas and cms collaborations at the lhc , have a program with very forward detectors . it is aimed to investigate soft and hard diffraction , low x dynamics with forward jet studies , high energy photon - induced interactions , large rapidity gaps between forward jets , and luminosity monitoring @xcite . these detectors will be located in a region nearly 220 - 420 m from the interaction point and they can detect protons in a continuous range of momentum fraction loss @xcite . momentum fraction loss of the proton is defined as @xmath4 . here @xmath5 is the momentum of the incoming proton and @xmath6 is the momentum of the intact scattered proton . therefore equipped with very forward detectors , lhc can to some extend be considered as a high - energy photon - photon or photon - proton collider . photon induced rections have been experimentally observed through @xmath7 processes in hadron - hadron collisions @xcite , @xmath8 in ep collisions @xcite , and pair production in aa collisions @xcite . these experiments raise interest on the potential of lhc as a photon - photon and photon - proton collider and motivate phenomenological works on photon - induced reactions at the lhc as a probe of new physics @xcite . anomalous @xmath9 couplings can be investigated in a model independent way by means of the effective lagrangian approach @xcite . we employ the following effective lagrangian describing anomalous @xmath9 couplings : @xmath10+h.c.\end{aligned}\ ] ] where @xmath11 it should be noted that lagrangian ( [ lagrangian ] ) also give rise to anomalous @xmath12 couplings . in the sm , the ( v -a ) coupling @xmath13 corresponds to the cabibbo - kobayashi - maskawa ( ckm ) matrix element @xmath14 , which is very close to unity and @xmath15 , @xmath16 and @xmath17 are equal to zero .    for off - shell top and/or bottom quarks the lagrangian in eq.([lagrangian ] ) is not the most general one , it should be extended with @xmath18 and @xmath19 terms where @xmath20 is the sum of the momenta of the t and b quarks . however , if wtb couplings arise from gauge invariant effective operators , single top production and decay can be described in full generality using the on - shell lagrangian in eq.([lagrangian ] ) for the wtb vertex , even in the process where the top and bottom quarks are far from their mass shell @xcite . if the w boson is off - shell , then there are additional terms containing @xmath21 @xcite . these terms are omitted in the lagrangian , they can be recovered by applying the equation of motion through operators of the original lagrangian @xcite . measurements at d0 detector at fermilab tevatron provide stringent direct constraints on these couplings @xcite . the most stringent bounds on the anomalous couplings @xmath22 and @xmath16 are given by @xmath23 and @xmath24 at 95% c.l . assuming that @xmath25 @xcite . recent results from early lhc data set comparable , but still weaker bounds with respect to tevatron @xcite . we see from the tevatron and lhc data that the bound on the coupling @xmath16 is weaker then the bound on @xmath17 and the bound on the coupling @xmath15 is weaker then the others @xcite . the ( v + a ) coupling @xmath15 is stringently bounded by the cleo @xmath26 data @xcite from an indirect analysis . limit from the cleo data is given by @xmath27 at @xmath28 level@xcite . it is explicit that indirect constraints are much more resrictive than direct constraints @xcite .    in the literature there has been a great amount of work on wtb couplings through single and pair top quark production . the single top quark production cross section was discussed below and above the @xmath29 threshold for the processes @xmath30 @xcite and @xmath31 at the cern lep @xcite . top quark single and pair production processes were studied for future linear @xmath32 collider and its @xmath33 and @xmath34 modes @xcite and also for @xmath2 collisions in tesla+herap and clic+lhc options @xcite . anomalous wtb couplings have also been probed at the lhc and tevatron @xcite . the equivalent photon spectrum of virtuality @xmath35 and energy @xmath36 is given by the following formula @xcite @xmath37\end{aligned}\ ] ] where @xmath38 in eq.(3 ) , e is the energy of the incoming proton beam and @xmath39 is mass of the proton . the magnetic moment of the proton is taken to be @xmath40 and @xmath41 @xcite . the photon spectrum which is integrated from a kinematic minimum @xmath42 up to @xmath43 is given by @xcite    @xmath44\end{aligned}\ ] ]    here the function @xmath45 is defined as follows @xmath46 + \frac{(1-b)y}{4x(1+x)^{3}}+c\left(1+\frac{y}{4}\right ) \nonumber \\\times\left[ln\frac{1+x - b}{1+x}+\sum^{3}_{k=1 } \frac{b^{k}}{k(1+x)^{k}}\right]\end{aligned}\ ] ] where @xmath47    @xmath48    in the epa emitted photons have a low virtuality and photon spectrum has a asymptotic behavior for large values of virtuality @xmath35 . in the epa that we have considered typical photon virtuality is @xmath49 @xcite . above this average virtuality value , spectrum function rapidly decreases and the contribution to the integral above @xmath50 is negligible . to be precise , the difference between sm cross sections for @xmath51 and @xmath52 is at the order of @xmath53pb . therefore during calculations we set @xmath54 . we consider the subprocesses @xmath55 and @xmath56 of our main process @xmath57 . in the sm single production of the top quark via the process @xmath58 is described by three tree level diagrams . each of the diagrams contains a wtb vertex . in the effective lagrangian approach , there are four tree level diagrams ; one of them contains an anomalous @xmath59 vertex , which is absent in the sm ( fig.2 ) .    the total cross section for the process @xmath0 can be obtained by integrating the cross section for the subprocesses over the photon and quark distributions : @xmath60\end{aligned}\ ] ] in this formula , @xmath61 and @xmath62 is the fraction which represents the ratio between @xmath63 @xmath64 quark and incoming proton s momentum . @xmath65 is the @xmath63 @xmath64 quark distribution function . we ignore interactions between different family quarks since the cross sections are suppressed due to small off diagonal elements of the cabibbo - kobayashi - maskawa matrix . in the total cross section calculations we have used martin , stirling , thorne and watt distribution functions @xcite .    in our calculations three different forward detector acceptance ranges have been discussed : @xmath66 , @xmath67 and @xmath68 . the former one was proposed by the atlas forward physics ( afp ) collaboration @xcite . the second acceptance range was proposed by the cms - totem forward detector scenario @xcite . since the forward detectors can detect protons in a continuous range of @xmath69 one can impose some cuts and choose to work in a subinterval of the whole acceptance region . hence , we also consider an acceptance of @xmath70 which is a subinterval of the cms - totem acceptance range .    in figs . 3 - 6 we show the integrated total cross section of the process @xmath71 as a function of anomalous couplings @xmath72 , @xmath73 , @xmath15 and @xmath74 for the acceptances of @xmath66 , @xmath75 and @xmath76 . here , the anomalous coupling @xmath74 is defined by @xmath77 . in figs . 3 and 4 we observe that cross section approximately has a same dependence to both @xmath78 and @xmath73 couplings . we see from fig . 5 that sensitivity of the cross section to anomalous coupling @xmath15 is comparably weak . we estimate the sensitivity of the process @xmath79 to anomalous couplings @xmath72 , @xmath73 , @xmath15 and @xmath74 using a simple one parameter @xmath80 criterion for integrated luminosities of @xmath81 and @xmath82 = 14 tev . the @xmath80 function is given by @xmath83 where @xmath84 is the statistical error . the expected number of events has been calculated considering the leptonic decay channel of the w boson and leptonic decay of the top quark as the signal @xmath85 , where @xmath86 or @xmath87 . atlas and cms have central detectors with a pseudorapidity coverage @xmath88 . therefore we consider an acceptance window of @xmath88 for final state electrons , muons and b quark . branching ratios appearing in the number of events are defined as @xmath89 where @xmath90 is the full width and @xmath91 is the decay rate for the corresponding channel with a cut of @xmath92 for final decay products . the limits for the anomalous coupling parameters are given in tables 1 - 4 for integrated luminosities of @xmath93 and forward detector acceptances of @xmath66 , @xmath94 and @xmath68 . we see from the tables that @xmath95 case provides more restrictive bounds on both @xmath96 and @xmath73 couplings compare to the @xmath66 and @xmath94 cases . on the other hand , limits on @xmath15 and @xmath74 couplings in @xmath97 case are weaker then the limits in @xmath66 and @xmath67 cases . the limits presented in tables are reasonable . in the effective lagrangian ( 1 ) anomalous couplings @xmath15 and @xmath13 are originated from dimension 4 effective operators but anomalous couplings @xmath17 and @xmath16 are originated from dimension 5 effective operators . therefore energy dependence of the terms in the cross section proportional to @xmath17 and @xmath16 are expected to be higher than the standard model and also the terms in the cross section proportional to @xmath15 and @xmath13 . hence , the main new physics contribution from couplings @xmath17 and @xmath16 comes from high energy region @xmath97 . on the contrary , the main standard model contribution comes from low energy region @xmath98 . the limits on @xmath17 and @xmath16 are expected to be better in the less forward region @xmath97 since the invariant mass of the incoming photon and b quark is large and the standard model background is low in that region . our limits on the couplings @xmath17 and @xmath16 are approximately a factor from 2 to 4 and 3 to 6 better then the limits from direct constraints at the tevatron respectively depending on the luminosity @xcite . on the other hand , our best limit on @xmath15 is a factor of 3.7 more restricted compared to tevatron direct constraint @xcite . physics studied at the lhc is significantly enhanced via the forward physics programs of atlas and cms collaborations . equipped with forward detectors lhc gives us new options to examine high energy photon - proton interactions . with respect to pure deep inelastic scattering processes , photon - proton interactions provide a quite clean channel due to absence of one of the incoming proton remnants . furthermore , detection of the intact scattered protons in the forward detectors allows us to reconstruct quasi - real photons momenta . this provides an advantage in reconstruction of the kinematics . the integrated total cross section of the process @xmath99 as a function of anomalous coupling @xmath17 for three different detector acceptances stated in the figure . the center of mass energy is taken to be @xmath100 tev . ] the integrated total cross section of the process @xmath99 as a function of anomalous coupling @xmath16 for three different detector acceptances stated in the figure . the center of mass energy is taken to be @xmath100 tev . ] the integrated total cross section of the process @xmath99 as a function of anomalous coupling @xmath15 for three different detector acceptances stated in the figure . the center of mass energy is taken to be @xmath100 tev . ] the integrated total cross section of the process @xmath99 as a function of anomalous coupling @xmath74 for three different detector acceptances stated in the figure . the center of mass energy is taken to be @xmath101 tev . ] k. piotrzkowski , phys . rev . d 63 , 071502 ( 2001 ) . v. m. budnev , i. f. ginzburg , g. v. meledin and v. g. serbo , phys . 15 , 181 ( 1975 ) . g. baur et al . 364 , 359 ( 2002 ) . c.royon , mod.phys.lett.a,18,2169(2003 ) . m.boonekamp , r.peschanski , and c.royon , phys.rev.lett.87 , 251806(2001 ) . m.boonekamp , r.peschanski , and c.royon , nucl.physb , 669,277 ( 2003 ) . m.boonekamp , a de roeck , r.peschanski , and c.royon , phys.lett.b 550,93 ( 2002 ) . v.a.khoze , a.d.martin , and m.g . ryskin , eur . c 19 , 477 ( 2001 ) . v.a.khoze , a.d.martin , and m.g . ryskin , eur . c 24 , 581 ( 2002 ) . v.a.khoze , a.d.martin , and m.g . ryskin , eur . c 55 , 363 ( 2008 ) . v.a.khoze , a.d.martin , and m.g . ryskin , phys . b 650 , 41 ( 2007 ) . a.b.kaidalov , v.a.khoze , a.d.martin , and m.g . ryskin , eur . c 33 , 261 ( 2004 ) . a.b.kaidalov , v.a.khoze , a.d.martin , and m.g . ryskin , eur . c 31 , 387 ( 2003 ) . o.kepka and c.royon , phys.rev.d 78 , 073005 ( 2008 ) . v.a.khoze , a.d.martin , and m.g . ryskin , eur . c 23 , 311 ( 2002 ) . n.schul and k.piotrzkowski , nucl.phys.b , proc.suppl . 179 - 180 , 289 ( 2008 ) . s.r.klein and j.nystrand phys.rev.lett . 92 , 142003 ( 2004 ) . gay ducati and g.g.silveira , phys.rev . d 82 , 073004 ( 2010 ) . v.p . goncalves and m.v.t.machado , phys.rev.d75 , 031502 ( 2007 ) . m.v.t.machado , phys.rev.d 78 , 034016 ( 2008 ) . x. rouby , phd thesis _ measurements of photon induced processes in cms and forward proton detection at the lhc _ , ucl - thesis 135 - 2008 , cms ts-2009/004 . a.abulencia _ et al . _ ( cdf collaboration ) , phys . 98 , 112001 ( 2007 ) . ( cdf collaboration ) , phys . 102 , 242001 ( 2009 ) . ( cdf collaboration ) , phys . 102 , 222002 ( 2009 ) . a.aktas _ et al . _ , phys.lett.b 583 , 428 ( 2004 ) . a.aktas _ et al . c 31,17 ( 2003 ) . f.d.aaron _ et al . _ , eur . c 71,1578 ( 2011 ) . zeus collaboration , nucl.phys.b 800,1 - 76 ( 2008 ) . zeus collaboration , nucl.phys.b 816,1 - 61 ( 2009 ) . h1 collaboration , eur.phys.j.c 48,749 - 766 ( 2006 ) . a.belkacem _ et al . _ , phys.rev.a 56 , 2806 ( 1997 ) . et al . _ , phys.rev.a 50 , 2313 ( 1994 ) . r.baur _ et al . _ , phys.lett.b 332 , 471 ( 1994 ) . j.adams _ et al . _ , phys.rev.c 70 , 031902 ( 2004 ) . w. buchmuller , d. wyler , nucl . b 268 , 621 ( 1986 ) k. hagiwara , s. ishihara , r. szalapski , d. zeppenfeld , phys . rev . d 48 , 2182 ( 1993 ) g.j . gounaris , f.m . renard , c. verzegnassi , phys . rev . d 52 , 451 ( 1995 ) g.j . gounaris , f.m . renard , n.d . vlachos , nucl . b 459 , 51 ( 1996 ) k.whisnant , j. yang , b. young , x. zhang , phys . rev . d 56 , 467 ( 1997 ) j.m . yang , b. young , phys . rev . d 56 , 5907 ( 1997 ) g. kane , g. ladinsky , c .- p . yuan , phys . d 45 , 124 ( 1992 ) j.  a.  aguilar - saavedra , nucl .  phys . b * 804 * , 160 ( 2008 ) [ arxiv:0803.3810 [ hep - ph ] ] e. boos , a. pukhov , m. sachwitz , h.j . schreiber , phys.lett . b404 , 119 - 123 ( 1997 ) v.m.abazov _ et al._(d0 collaboration ) phys.rev.lett . 101 , 221801 ( 2008 ) ; arxiv:0807.1692 [ hep - ex ] v.m.abazov _ et al._(d0 collaboration ) phys.rev.lett . 102 , 092002 ( 2009 ) . d0 collaboration , d0 conference note 5838 ( 2009 ) v.m.abazov _ et al._(d0 collaboration ) phys . b * 708 * , 21 ( 2012 ) ; arxiv:1110.4592 [ hep - ex ] . j.a.aguilar-saavedra , n.f.castro , and a.onofre , phys.rev.d 83 , 117301 ( 2011 ) cleo collaboration , m.s . alam et al . , phys . lett.74 , 2885 ( 1995 ) f. larios , m.a . perez , c .- p . yuan , phys . lett . b 457 , 334 ( 1999 ) j.drobnak , s.fajfer , and j.f.kamenik , nucl .  phys . b * 855 * , 82 ( 2012 ) [ arxiv:1109.2357 [ hep - ph ] ] k.kolodziej , phys . b * 710 * , 671 ( 2012 ) [ arxiv:1110.2103 [ hep - ph ] ] . s. ambrosanio , b. mele , z. phys . c 63 , 63 ( 1994 ) n.v . dokholian , g.v . jikia , phys . b 336 , 251 ( 1994 ) k. hagiwara , m. tanaka , t. stelzer , phys . lett . b 325 , 521 ( 1994 ) e. boos et al . , phys . lett . b 326 , 190 ( 1994 ) b. grzadkowski , z. hioki , phys . rev . d 61 , 014 013 ( 1999 ) b. grzadkowski , z. hioki , nucl . phys . b 585 , 3 ( 2000 )                        d. o. carlson and c .- p . yuan , phys . lett . b 306 , 386 ( 1993 ) . d. o. carlson , e. malkawi and c .- p . yuan , phys . * b 337 * , 145 ( 1994 ) . t. stelzer and s. willenbrock , b 357 , 125 ( 1995 ) . m. c. smith and s. willenbrock , phys . rev . d 54 , 6696 ( 1996 ) . z.hioki , k.ohkuma , phys . rev . d 83 , 114045 ( 2011 ) [ arxiv:1104.1221 [ hep - ph ] ] . s.k.gupta , a.s.mete , g.valencia , phys . d 80 , 034013 ( 2009 ) [ arxiv:0905.1074 [ hep - ph ] ] . sensitivity bounds of the coupling @xmath72 for various forward detector acceptances and integrated lhc luminosities . the center of mass energy of the proton - proton system is taken to be @xmath102 tev . [ cols="^,^,^,^,^,^,^",options="header " , ]
Abstract:

Answer: we study the possibility for the process @xmath0 with anomalous wtb couplings in a model independent effective lagrangian approach at the lhc . we find @xmath1 confidence level bounds on the anomalous coupling parameters for various values of the integrated luminosity . the improved constraints on the anomalous wtb couplings have been obtained compared to current limits .


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 interaction between light and matter is one of the fundamental topics in physics . modern techniques allow for fabricating artificial atoms with highly controllable parameters , which provide an ideal testing ground for quantum optics . in the emerging field of circuit quantum electrodynamics ( qed ) @xcite , where superconducting qubits , playing the role of artificial atoms , are strongly coupled to a superconducting transmission line resonator , many quantum optics effects have been demonstrated with unprecedented accuracy , and some novel phenomena have been observed . an example is single - qubit lasing @xcite , where quantum noise influences the linewidth of the emission spectrum in a characteristic way @xcite . recently proposals to replace in circuit qed setups the superconducting qubits by semiconductor quantum dots were put forward @xcite . quantum dots have the advantages that many of the parameters , e.g. , the level structure , coupling strength , and tunneling rates , can be varied _ in situ _ , and details of the transport properties can be resolved via adjacent quantum point contacts @xcite . fast experimental progress has been made towards coupling quantum dots to a superconducting resonator @xcite .    in this paper we consider such a circuit qed setup with a semiconductor double quantum dot coupled to a high - q transmission line . we find that by driving a current through a suitably biased double dot , a population inversion can be created between the dot levels , which then may lead to a lasing state in the resonator . this lasing state in the resonator correlates with the transport properties through the double dot system . this property not only allows probing the lasing state via a current measurement , which may be easier to perform in an experiment , but also , because of the narrow resonance window , opens perspectives for applications for high resolution measurements . the paper is organized as follows . in sec . [ sec : mm ] we present the model for the coupled quantum dot - resonator system and the pumping mechanism , provide estimates for the parameters , and we introduce the theoretical tools for the analysis of the system . in sec . [ sec : sp ] we investigate the stationary properties of the radiation field inside the resonator , while its emission spectrum and linewidth are discussed in sec .  [ sec : lw ] . in sec . [ sec : clt ] we analyze the correlation between the lasing and transport properties of this coupled system . we conclude with a brief summary . we consider a system shown schematically in fig . [ fig : dot ] . the double quantum dot is biased such that three charge states play a role , a reference neutral state and two states differing by a single extra electron in the left or right dot . the three states are denoted by @xmath0 , @xmath1 , and @xmath2 , respectively . the two singly charged states , with energy difference @xmath3 ( denoted as dot detuning ) , are coupled by coherent interdot tunneling with strength @xmath4 . the double dot system is therefore described by a two - level hamiltonian , @xmath5 with pauli matrices @xmath6 and @xmath7 . we model the superconducting transmission line as a harmonic oscillator with frequency @xmath8 , and @xmath9 and @xmath10 denoting the annihilation and creation operators for the radiation field . in typical circuit qed experiments , the resonator frequency is of the order of several ghz , and the vacuum - fluctuations - induced voltage drop @xmath11 between its central strip and ground plane is of the order of @xmath12v @xcite . for definiteness we choose @xmath13 throughout this paper . for the geometry shown in fig . [ fig : dot ] , the superconducting resonator is capacitively coupled to the charge states of the double dot which have different dipole moments . it leads to the interaction @xmath14 when the interdot tunneling is weak and the single - particle wavefunctions are strongly localized in either dot , the coupling strength can be estimated as @xmath15 , where @xmath16 is the distance between the centers of the two dots and @xmath17 the electric field at the antinode of the resonator mode . this estimate agrees with the experimental observations  @xcite .    in the eigenbasis of the two - level hamiltonian ( [ eq : dd ] ) , within the rotating wave approximation , the hamiltonian for the coupled system reduces to @xmath18 here @xmath19 ( similar for @xmath20 ) is defined in this subspace spanned by eigenstates @xmath21 with angle @xmath22 characterizing the mixture of the charge states . the energy splitting between the eigenstates , @xmath23 , can be tuned via gate voltages , allowing control of the detuning relative to the resonator frequency @xmath24 . the effective dot - resonator coupling strength , @xmath25 , reaches its maximum at @xmath26 . to achieve a lasing state , the detuning @xmath27 should be small and the coupling @xmath28 strong enough to overcome the dissipative processes . this requires both the energy difference @xmath3 and the interdot tunneling @xmath4 to be of the order of several @xmath29 .     and @xmath30 of the two dots are assumed to be arranged as indicated.,scaledwidth=35.0% ]    in order to create a population inversion , which is usually required for lasing behavior , we consider a pumping scheme involving the third state @xmath0 , similar to an optical laser . ( note that lasing without inversion is also possible in circuit qed setups where the dissipative environment enhances the photon emission as compared to absorption @xcite . population inversion induced by other pumping mechanisms , e.g. , driving the two - level system coherently @xcite , were also discussed . ) as indicated in fig . [ fig : pump ] , we assume that a bias voltage is applied across the double dot such that only the chemical potentials of the states @xmath1 and @xmath2 lie within the bias window . at low temperature compared to the charging energy , the only possibility for an electron to tunnel into the dot system is from the left lead to the left dot , leading to the transition from state @xmath0 to @xmath1 with tunneling rate @xmath31 . similarly , an electron can tunnel out into the right lead , creating a transition from @xmath2 to @xmath0 with tunneling rate @xmath32 . the pumping leads to a non - equilibrium state where the population in the state @xmath1 with higher energy is enhanced compared to that in @xmath2 . the population inversion persists when we go to the eigenbasis , where it is given by @xmath33 $ ] @xcite .    to analyze the dynamics of the coupled dot - resonator system , we employ a master equation for the reduced density matrix @xmath34 in the born - markovian approximation , @xmath35    + \mathcal l_{\rm r}\ , \rho + \mathcal l_{\rm \downarrow}\rho+\mathcal l_{\rm \varphi}\ , \rho    + \mathcal l_{\rm l}\ , \rho+\mathcal l_{\rm r}\ , \rho                              \nonumber \\   & \equiv & \mathcal l_{\rm tot } \ , \rho \ , .\end{aligned}\ ] ] we assume that the dissipative environment is characterized by smooth spectral functions @xcite and can be described by lindblad operators of the form @xmath36 for the oscillator we take the standard decay terms @xmath37 with rate @xmath38 . throughout this paper we consider low temperatures , @xmath39 , with vanishing thermal photon number and excitation rates . for the two - level system we account for the relaxation by @xmath40 with rate @xmath41 and for pure dephasing by @xmath42 with rate @xmath43 . for the discussion of lasing behavior , it is worthwhile to introduce a total decoherence rate for the dot system , namely , @xmath44/2 + \gamma_\varphi^*$ ] . the last two terms account for the incoherent tunneling between the electrodes and the left and right dots . in the basis of the pure charge states , the lindblad operators describing these processes are given by @xmath45 and @xmath46 . for simplicity , we set @xmath47 throughout the paper . we start with analyzing the stationary properties of the radiation field inside the resonator , namely the average photon number @xmath48 and the fano factor @xmath49 , which serve as good indicators of the state of the radiation field .     and fano factor @xmath50 as functions of dot detuning @xmath3 for vanishing decoherence of the dot levels ( @xmath51 ) . the analytical approximation @xmath52 ( dotted line ) in eq . ( [ eq_approximate_result_for_n ] ) is compared to the numerical result ( solid line ) . the decay rate of the resonator is chosen to be @xmath53 . throughout this paper we choose the incoherent tunneling rate @xmath54 , the coherent interdot tunneling strength @xmath55 , and the bare coupling strength @xmath56 . lower panel : photon number distribution at points a , b , c , and d marked in the upper panel . dotted line indicates a poissonian distribution for comparison with the sub - poissonian distribution at point d. , title="fig:",scaledwidth=53.0% ]   and fano factor @xmath50 as functions of dot detuning @xmath3 for vanishing decoherence of the dot levels ( @xmath51 ) . the analytical approximation @xmath52 ( dotted line ) in eq . ( [ eq_approximate_result_for_n ] ) is compared to the numerical result ( solid line ) . the decay rate of the resonator is chosen to be @xmath53 . throughout this paper we choose the incoherent tunneling rate @xmath54 , the coherent interdot tunneling strength @xmath55 , and the bare coupling strength @xmath56 . lower panel : photon number distribution at points a , b , c , and d marked in the upper panel . dotted line indicates a poissonian distribution for comparison with the sub - poissonian distribution at point d. , title="fig:",scaledwidth=40.0% ]    we plot both as functions of the detuning in fig . [ fig : sp ] . to begin with we neglect the relaxation and decoherence of the dot levels , @xmath51 . when the detuning is large ( e.g. , the point / curve marked by ` a ' in fig . [ fig : sp ] ) , the quantum dot effectively does not interact with the resonator , and the photon number is low with a thermal distribution . for weaker detuning ( e.g. , point ` b ' ) the photon number increases sharply and saturates at resonance , @xmath57 . approximately , the photon number can be estimated as @xcite , @xmath58 in the typical parameter regime of experiments the peak value of the photon number is mainly given by the first term , i.e. , the ratio between the incoherent pumping and the decay rate of the resonator . for large photon number good agreement is reached between the analytical expression ( dotted line ) and numerical results ( solid line ) . statistical properties of the radiation field can be obtained by investigating the fano factor . when the resonator is in a lasing state ( e.g. , point ` c ' ) , the radiation field is in the coherent state and the photon number distribution is poissonian . in this case the fano factor equals to 1 . in the strong - coupling regime , the fano factor can even become smaller than 1 ( e.g. , point ` d ' ) , which indicates a sub - poissonian distribution of the radiation field . in this non - classical regime , the photon number distribution is squeezed compared to the poissonian distribution .     and fano factor @xmath50 as functions of detuning . solid lines represent the results without decoherence , blue dotted lines are those for relaxation rate @xmath59 and vanishing pure dephasing , while red dashed lines are for pure dephasing rate @xmath60 but vanishing relaxation . other parameters are chosen to be the same in fig . [ fig : sp ] . , title="fig:",scaledwidth=44.0% ] +   and fano factor @xmath50 as functions of detuning . solid lines represent the results without decoherence , blue dotted lines are those for relaxation rate @xmath59 and vanishing pure dephasing , while red dashed lines are for pure dephasing rate @xmath60 but vanishing relaxation . other parameters are chosen to be the same in fig . [ fig : sp ] . , title="fig:",scaledwidth=44.0% ]    we now turn to the effects of decoherence of the dot levels on the lasing state , which are illustrated in figs . [ fig : deccharge ] . the relaxation process reduces the pumping efficiency and deteriorates the lasing state . hence the photon number decreases . at the same time the fano factor around resonance becomes larger than 1 , resulting from an increase of amplitude fluctuations . pure dephasing , on one hand , diminishes the efficiency of energy exchange between the dot levels and the resonator . on the other hand , it effectively broadens the window in which the dot levels can interact with the resonator . as a result , the peak of the photon number is broadened , while the peak value decreases . ( in arbitrary units ) for different detunings . , scaledwidth=43.0% ]    we now turn to the spectral properties of the radiation field . the emission spectrum of the radiation field , @xmath61 which can be measured directly in experiments , provides further insights into the state of the radiation field . to calculate the phase correlator , we make use of the quantum regression theorem @xcite , @xmath62,\end{aligned}\ ] ] with @xmath63 being the stationary solution of the master equation ( [ eq : me ] ) . the spectral properties of the radiation field in the quantum dot - resonator system are similar to those of a single - qubit maser discussed in refs . @xcite . in fig . [ fig : spectrum ] we plot the emission spectrum for different values of the detuning . at resonance , @xmath57 , the emission spectrum shows a sharp peak at the resonator frequency @xmath8 . for nonvanishing detuning the peak of the emission spectrum is shifted in frequency by @xmath64 . for weak detuning , this shift grows linearly with the detuning , @xmath65 , with @xmath66 characterizing the total dissipation rate of the quantum dot system . at the same time , the peak value of the emission spectrum decreases with growing detuning , related to a lowering of the average photon number .    , the blue dotted line for relaxation rate @xmath67 but vanishing pure dephasing , and the red dashed line for pure dephasing rate @xmath68 but @xmath69 . here the decay rate of the resonator is chosen to be @xmath70 . , scaledwidth=43.0% ]    the linewidth of the emission spectrum is plotted as a function of detuning in fig . [ fig : lw ] . in a typical situation where the coupling @xmath28 is not too strong compared to the decoherence of the coupled system @xmath71 , the linewidth grows monotonously with increasing detuning , as indicated , e.g. , by the blue dotted line . at resonance , deep in the lasing regime , the linewidth can be estimated as @xmath72 @xcite . it becomes very small for large photon numbers . here , due to the restrictions of our numerical calculation , we illustrate to a situation with low photon number around 15 and not very small values of @xmath73 . for stronger detuning the system leaves the lasing regime , and the linewidth is approximately given by @xmath74 , which approaches the bare linewidth of the resonator , @xmath75 , for strong detuning . when the coupling is much stronger than the decoherence rate , the linewidth depends nomonotonously on the detuning exhibiting an enhancement at the lasing transition ( green solid line ) @xcite . different from the conventional lasing situation where the emission spectrum is mainly broadened by phase fluctuations , in the strong coupling regime at the lasing transition , the contribution from amplitude fluctuations is magnified via its coupling to the phase fluctuations , leading to the peak shown in the figure . in this transition regime , the phase correlator @xmath76 exhibits a non - exponential decay in time . a further surprising phenomenon occurs for a strong pure dephasing rate @xmath77 . instead of broadening the emission spectrum , the pure dephasing actually reduces the linewidth ( red dashed line ) , leading to a much sharper spectrum peak . an intuitive interpretation can be obtained by examining the linewidth at resonance , where a semi - quantum approach is accessible @xcite . in this case , the linewidth at resonance in the deep lasing regime can be estimated as @xmath78 . with growing pure dephasing , the photon number changes little while the total decoherence rate @xmath79 increases linearly with the pure dephasing rate , which leads to the decrease in the linewidth . the parameters are chosen to be the same as in fig . [ fig : sp ] . the dot levels become resonant with the resonator at @xmath80 . , scaledwidth=40.0% ]    a current flows in the system if the tunneling cycle @xmath81 is completed . this can be achieved by either coherent interdot tunneling or a coherent process involving the excition of photons in the resonator . we calculate the current using @xmath82 where the index @xmath83 refers to the states @xmath84 , @xmath85 , and @xmath0 , and @xmath86 denotes the transition rate from state @xmath87 to @xmath88 . in fig . [ fig : curr ] we plot the current as a function of the dot detuning @xmath3 . as expected , two resonance peaks show up in the current . the coherent interdot tunneling leads to the broad peak around @xmath26 with width given by the tunneling rate @xmath4 ( here @xmath89 )  @xcite . interestingly , the second peak due to transitions involving the excitation of photons is  for realistic values of the parameters  much narrower . as shown in fig . [ fig : currdec ] , it correlates with the lasing state , since the excitation of a photon in the resonator is caused by the electron tunneling between the two dots . both the lasing state and the current peak exist only in a narrow  resonance window " @xmath90 . this window can be estimated from eq . ( [ eq_approximate_result_for_n ] ) with condition @xmath91 , namely , @xmath92\kappa\,\gamma}-1},\end{aligned}\ ] ] which reduces to @xmath93 for small @xmath94 . to estimate the height of the narrow current peak , we adopt an adiabatic approximation assuming the dynamics of the resonator to be much slower than that of the quantum dots @xcite . for @xmath95 and small @xmath94 , the peak value of the current is given by @xmath96.\end{aligned}\ ] ] here @xmath97^{-1}$ ] denotes the probability of having @xmath98 photons in the resonator . when the coupling to the resonator is strong compared to the incoherent tunneling @xmath99 , the peak current approaches @xmath100 . the correlation between the lasing state and the transport current is remarkable in two ways . on one hand , the current peak , which may be easier to measure than the photon state of the resonator , can be used as a probe of the lasing state . on the other hand , the rather sharp resonance condition needed for the lasing makes the current peak narrow , while at the same time the value of the current is reasonably high . this allows resolving in an experiment small details of the dot properties .     with the green solid line for the result with vanishing decoherence @xmath101 , the blue dotted line with relaxation rate @xmath67 but vanishing pure dephasing , and the red dashed line with pure dephasing rate @xmath68 but @xmath69 . the photon number ( gray dot - dashed line ) is also shown for comparison . , scaledwidth=42.0% ]    the effects of decoherence on the transport current are illustrated in fig . [ fig : currdec ] . relaxation , on one hand , reduces the efficiency of the coherent transition between the excited and ground states . hence within the resonance window the current induced by the coherent transition decreases . on the other hand , relaxation also opens up an incoherent channel which increases the current . pure dephasing leads to a broadening of the resonance peak in the current , as well as of the photon number . we have investigated a double quantum dot coherently coupled to a transmission line resonator . a population inversion between the dot levels is created by a pumping process due to a voltage applied across the double dot . a lasing state of the radiation field develops within a sharp resonance window . it correlates with a peak in the transport current . the sharp resonance condition allows for resolving small differences in the dot properties . this opens perspectives for applications of the setup and operation principle for high resolution measurements . we analyzed the effects of dissipation on the dot levels . both relaxation and pure dephasing deteriorate the lasing state , accompanied by a decrease in the photon number and an increase of amplitude fluctuations . relaxation processes shorten the life time of the two - level system and reduce the population inversion and hence the lasing effect . pure dephasing , which also shortens the life time , broadens the effective resonance window between the quantum dots and resonator . surprisingly , the emission spectrum can become even sharper with moderate pure dephasing rate . for the transport current , the relaxation opens up an extra incoherent channel which increases the current outside the resonance window with the resonator . we acknowledge fruitful discussions with s. andr , a. romito and j. weis , as well as the support from the baden - wrttemberg stiftung via the kompetenznetz funktionelle nanostrukturen and the dfg via the priority program semiconductor spintronics.
Abstract:

Answer: we investigate a double quantum dot coupled to a transmission line resonator . by driving a current through the double dot , a population inversion between the dot levels can be created , and a lasing state of the radiation field is generated within a sharp resonance window . the transport current correlates with the lasing state . the sharp resonance condition allows for resolving small differences in the dot properties . dissipative processes in the quantum dots and their effect on the lasing and transport behavior are investigated .


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: grs 1915 + 105 is a well - known and somewhat unusual example of a galactic x - ray binary system . it was first observed in outburst on 1992 august 15 @xcite . it is believed to contain a black hole that is both maximally rotating and rather massive at @xmath3 solar masses @xcite . grs 1915 + 105 was also the first known galactic example of a  superluminal source " @xcite , so called because the relativistic jets that are ejected from it move with an apparent velocity greater than @xmath4 . unlike most transient low - mass x - ray binaries , grs 1915 + 105 has never been observed in quiescence since it first became active . high absorption along the line of sight to grs 1915 + 105 makes observation in optical and ir wavelengths difficult . however , a @xmath0-band counterpart has been discovered @xcite . here we report daily monitoring of the @xmath0-band counterpart of grs 1915 + 105 . in  [ secobs ] we review our observation and data reduction methods , and present @xmath0-band and x - ray light curves . section [ seccor ] features analysis of the correlation between ir flux , x - ray flux , and x - ray hardness ratio . in  [ secper ] we discuss extraction of the orbital period , its phase and amplitude , and physical implications . finally ,  [ secsup ] concerns the second ,  superhump " period , leading to limits on the mass ratio of the system . section [ seccon ] reiterates our conclusions . our observations of grs 1915 + 105 span seven seasons , from mjd 51652 to 53969 ( 2000 april 18 to 2006 august 21 ) , with a total of 604 individual data points ( see table [ datatable ] . ) data from 2000 - 2002 were taken using the andicam instrument and the yale 1.0 m telescope at the cerro tololo inter - american observatory ( ctio ) , which was operated at the time by the yalo consortium @xcite . subsequent data were obtained using the same instrument mounted on the 1.3 m telescope originally constructed for the two micron all sky survey ( 2mass ) , also at ctio and operated by the small and moderate aperture research telescope system ( smarts ) consortium . all observations are internally dithered by a tiltable mirror to remove artificial effects such as bad pixels and to allow sky subtraction . due to the large extinction , our observations are in the @xmath0-band only . @xmath0-band images of grs 1915 + 105 were taken at five spatially offset positions , with five internally dithered images at each position for a total of 25 images per observation . when combining the images the spatial offset was removed first , then the internal dithering ; this process served to reduce artifacts introduced by irregularities in the filters during internal dithering . combination was done using the iraf software package , along with sky and flat subtraction and bias removal . exposure time for each individual image was 60 s for data taken from 2000 - 2002 , and 40 s for that from 2003 - 2006 , for a total exposure time of 1500 s and 1000 s , respectively . due to the larger aperture of the 1.3 m telescope , the signal - to - noise ratio generated by these exposure times is similar . differential photometry was performed on the reduced images in iraf , by comparing the ir counterpart of grs 1915 + 105 to a set of non - variable reference stars . three different sets of reference stars were used in all : one for data from 2000 , one from 2001 - 2003 , and one from 2004 - 2006 . reference star sets were changed in the first instance to account for a bad quadrant on the detector , and in the second instance due to the change in the field of view as the instrument moved to a different telescope . calibrated magnitudes for the most recent data were obtained by comparison with 2mass catalog values of several reference stars . earlier differential photometry done with different sets of reference stars was then standardized by subtracting the difference between the average instrumental magnitudes of each reference set . x - ray counts are taken from the all - sky monitor ( asm ) instrument aboard the _ rossi x - ray timing explorer _ ( _ rxte _ ) satellite , over an energy range of 1.5 - 12 kev . for the purpose of considering x - ray hardness we define  hard " and  soft " x - ray bands as those above and below 5 kev , respectively , corresponding to the asm s own @xmath5 and @xmath6 bands . the combined @xmath0-band and x - ray light curve is presented in figure [ lc ] . the large gaps seen in the @xmath0-band data correspond to seasons when the line of sight to grs 1915 + 105 was near the sun , making observation impossible . the apparent rough correlation between the two bands , particularly near peak values , is characteristic of an x - ray binary system . note that the x - ray counts per second are always significant , so that the system never truly reaches a quiescent state over the course of our observations . the dotted line at @xmath7 indicates a magnitude cutoff used in the frequency analysis . in performing a frequency analysis , we are looking in particular for luminosity variations related to the orbital period of the secondary star . regardless of the exact nature of such magnitude variations , we expect generically that they will be drowned out by emission from the accretion disc when the system is brightest . including data in which the period is not detectable will only reduce the signal - to - noise ratio , so we excluded all points above a fixed magnitude cutoff from the calculation of the frequency spectrum . interpretation of this method in light of our particular model of the system , as well as a justification for the choice of the value @xmath8 mag , are given in  [ secper ] below .    to search for periodic structure in the data , the raw @xmath0-band results were transformed to the frequency domain . for time - domain data sets with large , irregular gaps such as ours , attempting a naive fourier analysis would introduce a large number of artifacts into the frequency domain , obscuring the true signature from the data . to avoid this effect , we used a one - dimensional implementation of the well - known clean algorithm @xcite that reduces the number of artifacts introduced into the power spectrum . the effect of  cleaning " on the raw , or  dirty , " spectrum is shown in figure [ clean ] , including a demonstration of the effect of the magnitude cut . conversion of @xmath0-band apparent magnitudes into fluxes was accomplished with reference to @xcite , according to which the two quantities are related by @xmath9 . here zp@xmath10 is a zero - point , which is given to be @xmath11 for the @xmath0-band . in the above formula , @xmath12 is taken to be in units of j s@xmath13 @xmath14 hz@xmath13 , which is equivalent to @xmath15jy . least - squares linear fitting of our results demonstrates a correlation between the x - ray and ir data . in particular both the x - ray counts per second and hardness ratio are positively correlated with the @xmath0-band flux ; this result is consistent with previous analysis of grs 1915 + 105 @xcite . the correlation coefficients obtained are @xmath16 and @xmath17 , respectively , so the ir / hardness correlation is stronger . both correlation coefficients are found to be different from zero at a @xmath18-level below @xmath19 @xcite , indicating that our result is statistically significant . the relation between @xmath0-band flux and x - ray hardness ratio is depicted in figure [ scatter ] . it is clear from figure [ scatter ] that there is no significant correlation at large @xmath0-band flux . since the correlation is established entirely by the points with low @xmath0 flux , our linear fit to the full data set is not physically relevant . this suggests that there are several sources of ir flux , not all of which are correlated with the x - ray properties of the source . thermal ir flux from the outer parts of the accretion disk is expected to correlate with the x - rays , since both depend on the accretion flow through the disk . on the other hand , large non - thermal ir flares without an easily observed x - ray response have been reported in several other transients @xcite . we therefore suggest that the ir flux observed in grs 1915 + 105 includes both a thermal component associated with the disk , and occasionally larger , possibly non - thermal emission uncorrelated with the x - rays . we also expect a contribution from the secondary star , which will be discussed in more detail in the next section . unfortunately , our ability to interpret the ir flux in greater detail is limited by the fact that we have ir observations in only one band , and so have no ir color information . conversion of our ir observations into frequency space , cleaned and cut at 12.5 mag , reveals a significant peak corresponding to a period of @xmath20 days . the error is estimated by fitting a gaussian to the frequency spectrum at the peak frequency and taking the standard deviation of the fit to be the associated error @xmath21 . monte carlo statistical testing @xcite over @xmath22 trials yields a significance level of the period that is indistinguishable from @xmath23 . folding and binning the data according to phase at this period makes the periodic structure manifest in the time domain ( see figure [ foldphase ] . ) furthermore , our period agrees to @xmath24 with previous spectroscopic observations of grs 1915 + 105 , which determined an observed period of @xmath25 days @xcite .    for @xmath26 days , we find that our phase corresponds to a point of minimum brightness at mjd @xmath27 . early in our light curve , minimum brightness occurred at mjd @xmath28 , in extremely good agreement with the blue - to - red radial velocity crossing observed by @xcite at mjd @xmath29 . this agreement is consistent with interpretation of the orbital modulation as being due to the heated face of the secondary star .    before proceeding with our analysis of the frequency spectrum , we introduce a simple physical model of the system , assuming that the heated face of the secondary is responsible for the orbital variation . in addition to lending physical significance to our measurements , the model allows us to validate our method of omitting data below a particular cutoff from our analysis . in particular we seek to describe the size of the variation of the observed magnitude due to the orbital period of the donor star . if the side of the donor star facing the compact object is heated by emission from the accretion flow , then brightness variations will occur at the orbital period due to the rotation of the star , which is tidally locked to the orbital period . although the detailed behavior is complex , this model features a simple and robust prediction of the phase dependence . by purely geometric reasoning , we can see that the phase corresponding to minimum measured brightness also corresponds to the point of blue - to - red crossing in the companion star s radial velocity , when the heated face is completely occluded . as noted above , the observed phase of our photometric observations agrees with that of the spectroscopy of @xcite , confirming this prediction . based on this model of grs 1915 + 105 , our use of a magnitude cutoff seems to be reasonable . when the accretion disc is emitting strongly we would expect the brightness variations due to the heated face of the star to be diluted . the exact choice of @xmath7 for the cut is arbitrary , but not unreasonable . through experimentation with various cuts , we determined that the detected period of interest is quite robust against the choice of cut ( see figure [ cuts ] ) . we adopted the cut at @xmath8 mag because it yields the strongest orbital signal . presumably , cuts below this value reduce the number of relevant data points included , while cuts above this value include points at which the orbital variation is suppressed by strong emission from the compact object ; both result in a reduced signal - to - noise ratio . the cut is depicted against the light curve in figure [ lc ] .    in order to interpret the amplitude of the fourier component at @xmath30 , we must adopt some additional assumptions . we model the heated face of the star as a sphere which has one hemisphere brighter than the other . adding a coordinate system and doing some simple vector analysis then yields an expression for @xmath31 , the variable part of the flux from the system . the resulting formula is    @xmath32    where @xmath33 is the flux of the total system in the absence of heating , @xmath34 is a constant of proportionality relating the variation @xmath31 to the total emitted flux @xmath33 , and @xmath35 is the angle between the line of sight and the line from the compact object to the star . note that the range of @xmath35 is determined by the inclination angle @xmath36 of the orbital system : @xmath37 $ ] . in particular if @xmath38 there is no variation in the flux , since at all times half of each of the bright and dark hemispheres will be visible .    converting to magnitudes and taking the total range of variation , we find @xmath39 where @xmath36 should be expressed in radians . at @xmath26 days , we find the amplitude of the magnitude variation ( i.e. the magnitude of the cleaned frequency spectrum component at the period of interest ) to be @xmath40 mag ; the total range of variation @xmath41 is equal to twice this . using an inclination angle of @xmath42 @xcite , we find @xmath43 . in other words , our model predicts that the flux emitted by grs 1915 + 105 is roughly @xmath44 greater when the heated face of the secondary star is maximally visible . the magnitude variation can also be related to the temperature difference @xmath45 between the two hemispheres of the star . we assume a blackbody spectrum , so that the distribution of intensity with respect to wavelength follows planck s law . integrating over the range of wavelengths corresponding to our @xmath0-band filter then yields the @xmath0-band flux .    by determination of the spectral type , the temperature of the secondary star has been estimated to be @xmath46 k @xcite . our @xmath0-band filter includes wavelengths on the order of @xmath47 m , which is well into the tail of the planck distribution at the given temperature . therefore we take the @xmath48 approximation , which gives a flux directly proportional to temperature . in terms of @xmath41 , we find    @xmath49\ ] ]    where @xmath50 is the temperature of the unheated side of the star , and @xmath51 is the ratio between the flux emitted by the compact object and accretion flow ( @xmath52 ) and the unheated side of the star ( @xmath33 ) . taking @xmath53 as above and @xmath46 k as noted , we find a lower bound of @xmath54 k on the temperature difference if @xmath55 , which is expected since grs 1915 + 105 has never reached a quiescent state . if the quiescent magnitude of the system is measured , then @xmath33 will be known , and @xmath45 and @xmath52 will be able to be determined for typical outburst values . also apparent in our frequency analysis was the existence of a period slightly longer than the orbital period . exact determination of this period is not possible , as it wanders over a range of values depending on which subset of the data is examined , as seen in figure [ cuts ] . specifically , it appears to range from @xmath56 to @xmath57 days . separating our data into  rising " and  falling " subsets , corresponding to how the magnitude is varying over sets of three adjacent data points , produces a striking difference in the frequency spectrum , depicted in figure [ superhump ] . the secondary period is strongly visible in the falling data , but vanishes completely in the rising data . the peak observed in the falling plot corresponds to a period of @xmath58 days . a period such as this one , differing by a few percent from the orbital period of the system , is indicative of a superhump , a resonance in the accretion disc itself frequently observed in cataclysmic variables ( cvs ; @xcite ) , as well as low - mass x - ray binaries @xcite . the mass ratio of grs 1915 + 105 is certainly well below the maximum value of @xmath59 for the existence of a superhump resonance @xcite . the difference between the  rising " and  falling " data is also consistent with previous observations of superhumps , which are typically only seen in systems declining from outburst , although this may be somewhat of an empirical rule since observation of the rising magnitudes leading up to outburst is difficult in cvs . since @xmath60 is strictly larger than @xmath30 , we can conclude that the superhump precession is prograde with respect to the rotation of the accretion disc @xcite . finally , we determine the quantity @xmath61 , the fractional difference between the superhump and orbital periods , which is defined as @xmath62 . given our range of values for @xmath60 , we find that @xmath61 lies between @xmath63 and @xmath64 . there is a simple empirical formula relating @xmath61 to the mass ratio of the system @xmath65 , derived from observation of a number of low - mass x - ray binary and cv systems @xcite :    @xmath66    by this rule , the mass ratio @xmath65 of grs 1915 + 105 should lie between @xmath67 and @xmath68 . this is in agreement with @xcite , who find a mass ratio of @xmath69 . their result depends on assigning a mass to the observed spectral type of the secondary star , rather than the more reliable line - broadening technique ( e.g. @xcite ) . our result confirms that the masses deduced by greiner et al . for the components of this system are approximately correct . conversely , if we assume @xmath70 then @xmath71 , corresponding very closely to the superhump peak seen in the  falling " data subset ( figure [ superhump ] . ) the most significant periodicity seen in the frequency spectrum of our data is a feature at a period of approximately @xmath72 days . this is consistent with the spectroscopically determined orbital period of the secondary star , determined by @xcite , of @xmath25 days . therefore we claim @xmath20 days to be the orbital period , with a phase corresponding to minimum light at @xmath73 mjd @xmath27 . in addition , we find in our data another period ranging from @xmath56 to @xmath57 days , which we attribute to a  superhump " resonance in the accretion disc @xcite . the superhump period is longer than the orbital period , indicating that the precession is prograde with respect to the rotation of the disc @xcite . the fractional difference @xmath61 between the orbital and superhump periods is determined to lie between @xmath63 and @xmath64 . our measurement of @xmath61 allows us to make a prediction of the mass ratio @xmath65 of the system , based on an empirical formula @xcite . we find that @xmath2 , providing an independent verification of the result @xmath74 due to @xcite . the verification of the mass ratio @xmath65 via the superhump period is of particular importance . the previous determination of the mass ratio was based on spectroscopic identification of the secondary star @xcite , which requires a number of assumptions ; our method is based on an unrelated set of assumptions . conversely , our result can be seen as a verification of the empirical relation between @xmath61 and @xmath65 , in a system with an extremely long orbital period compared to other known superhumpers @xcite . we thank the asm/_rxte _ teams at mit and the _ rxte _ science operations facility and guest observer facility at nasa s goddard space flight center for providing the x - ray results . we are grateful for the efforts of yalo / smarts service observers d. gonzalez and j. espinoza , who collected the data reported here . this work was supported by nsf grant ast 04 - 07063 to c.d.b . cc 51652.3878 & 13.15 + 51653.3471 & 13.16 + 51654.3760 & 13.16 + 51655.4107 & 13.13 + 51657.3792 & 13.25 + 51658.4115 & 13.12 + 51659.3678 & 13.31 + 51660.3831 & 13.25 + 51661.3867 & 13.34 + 51663.3795 & 13.20 +
Abstract:

Answer: we present the results of 7 years of @xmath0-band monitoring of the low - mass x - ray binary grs 1915 + 105 . positive correlations between the infrared flux and the x - ray flux and x - ray hardness are demonstrated . analysis of the frequency spectrum shows that the orbital period of the system is @xmath1 days . the phase and amplitude of the orbital modulation suggests that the modulation is due to the heating of the face of the secondary star . we also report another periodic signature between 31.2 and 31.6 days , most likely due to a superhump resonance . from the superhump period we then obtain a range for the mass ratio of the system , @xmath2 .