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 emergence of cooperation among selfish individuals is an intensively studied problem @xcite . traditionally , the problem of cooperation is investigated by means of the game theoretical models of the prisoner s dilemma for pairwise interactions , and more generally public goods game for groups of interacting individuals . in particular , the public goods game is abundant in human society , e.g. , protecting the global climate and avoiding overfishing of the oceans @xcite . in the classical public goods game ( pgg ) , individuals engage in multiplayer interactions and decide simultaneously whether to contribute ( cooperate ) or not ( defect ) to a common pool . then the accumulated contributions by cooperators are multiplied by a factor large than one , i.e. , the so - called multiplication factor , and finally the resulting assets are shared equally among all group members irrespective of their initial decision . from the perspective of each individual , defection is clearly the rational decision to make as it yields the highest income compared to other members . thus , selfish individuals should decline to contribute and attempt to free ride on the other players contributions . however , if nobody decides to invest , the group fails to harvest the benefits of a collective investment , which drives the population into the tragedy of the commons @xcite . actually , the group is most successful if everybody cooperates , and hence the dilemma is caused by the selfishness of individual players .    to study the social dilemma in realistic situations , in the last decade the risk pgg @xcite , the optional pgg @xcite , the threshold pgg @xcite , the continuous pgg @xcite , and the ecological pgg @xcite have been developed based on the classical pgg from the viewpoint of realistic societies . on the other hand , several mechanisms for the evolution of cooperation in the pgg , such as punishment @xcite , reward @xcite , reputation @xcite , network reciprocity @xcite have been justified . in particular , complex interaction networks provide a natural and reasonable framework for studying the pgg in structured populations . within this framework , some aforementioned mechanisms , such as punishment and reward have been further studied @xcite . also , some other factors have been incorporated , such as noise @xcite , social diversity @xcite , and success - driven distribution @xcite . it is found that social diversity associated with the number and the size of the public goods game as well as the individual contribution to each game can greatly promote the emergence of cooperation @xcite . indeed , social diversity by means of the system s other feature information , e.g. , game payoffs @xcite , teaching activity @xcite and preferential selection @xcite in strategy updating , have been also demonstrated to facilitate cooperation in the pgg . it is worth mentioning that the inhomogeneities and social diversity about features of the system are widely existent in human society and animal world , which can characterize the asymmetric and different influence of individuals or interacting environments . however , they are introduced artificially in some previous studies mentioned above . indeed the inhomogeneities or social diversity can emerge spontaneously via the coevolutionary rules , since the values of property should be not invariable , but evolve based on the state of the system . in the context of evolutionary game theory , the adaptive features are often coupled with the strategy evolution . the coevolution of strategies and features of the model , e.g. , individual social ties ( for example see refs . @xcite ) , noise level @xcite , payoff matrices @xcite , capability of strategy transfer @xcite and individual learning rules @xcite , have been investigated in different evolutionary games , especially the prisoner s dilemma game ( see @xcite for a review ) . remarkably , recently lee et al . further proposed a multiadaptive prisoner s dilemma game where both the payoff matrices and the interaction structure are shaped by the behavior of the agents , and found that such multiadaptive mode can result in the coemergence of hierarchical structure and cooperation @xcite .    at present , we propose a coevolutionary rule in the pgg . we consider that each interacting group has its own multiplication factor , which evolves based on the local strategy distribution in the group and the global strategy distribution in the whole system . different from the setting in some previous works @xcite , this adaptive multiplication factor is used to measure the local interacting or cooperative environment in each group , rather than individual s feature or the whole system s interaction conditions . correspondingly , the multiplication factor represents the feedback return of the local investment to the public good , and a larger value of the multiplication factor enables a better investment return @xcite . structured populations provide a competent framework to describe this local feature , which is updated based on the local and the global level of cooperation . in the present study , in order to focus solely on the effects of the adaptive investment returns , we employ a square lattice where the number of group members is fixed and always the same . due to the diversity of local strategy distribution , the investing cooperators can get together in some groups , whereas in other groups they are sparse . the inhomogeneous distribution can induce different investment returns in different interacting groups , which may correspond to the phenomenon of uneven regional exploitation of the common resources in a society . for example , in some pastures the herdsmen may over - exploit the pasture resource by adding more and more animals to their herd , which may lead to the gradual desertification of the grassland . correspondingly , the socioeconomic returns from herding in these grazing areas decrease gradually and finally the economic losses are unavoidable . on the contrary , in other pastures the herdsmen may still use the grassland resources while at the same time considering the conservation of the ecosystem , and in such cases a higher socioeconomic return is likely . we adopt the state of the system , i.e. , the global cooperation level , as the criterion to measure whether the local cooperative environment is favorable or not . since the state of the system is evolving simultaneously , here we prefer the dynamical global cooperation level , instead of the static criterion in refs . moreover , in real situations the investment return is variable , but should be somewhat limited by external adjustment . in general , it should be limited in a certain range @xcite .    in this study , we assume that the enhancement factor in each group is updated based on the differences between the dynamical local and global cooperation levels , and limited between the lower @xmath0 and the upper @xmath1 limit . we study how the adaptive and bounded investment returns influence the evolution of cooperation in the spatial pgg . we find that this pgg model with the adaptive and bounded investment returns can effectively enhance cooperation in spatially structured populations , and that appropriately bounded limitations of the multiplication factor can result in the best cooperation level . we find further that , in comparison to the traditionally spatial pgg where only one invariable multiplication factor which is larger than one exists in the whole population , our proposed pgg model can produce a higher cooperation level when each multiplication factor is limited to change only between one and the group size . we start by presenting the results as obtained when the lower limit of the multiplication factor @xmath0 is equal to the opposite number of the upper limit @xmath1 , i.e. , @xmath2 . here @xmath3 is the limit value ( for the detailed definitions see the methods section ) . figure  [ fig1](a ) shows the cooperation level @xmath4 at equilibrium in the population in dependence on the feedback strength @xmath5 for different values of @xmath3 . we find that when @xmath6 cooperation can be promoted for larger @xmath3 . to be specific , when @xmath3 is small , e.g. , @xmath7 , full defection is achieved , irrespective of the values of feedback strength . this is because cooperators can not survive in structured populations if the maximum enhancement factor is not sufficiently large @xcite . when @xmath3 becomes larger , e.g. , @xmath8 , the cooperation level first increases and then slightly decreases . subsequently , it holds at about @xmath9 with increasing @xmath5 . while there is no limitation for the multiplication factor , although the cooperation level is very high , but a small amount of defectors can survive in the population even for high feedback strength , e.g. , @xmath10 . when moderate values of @xmath3 is set , e.g. , @xmath11 , full cooperation can be achieved for high feedback strength , and the cooperation level is very similar for other moderate values of @xmath3 . to further qualify the effects of @xmath3 on the evolution of cooperation , we present @xmath4 as a function of @xmath3 for different @xmath5 in fig . [ fig1](b ) . clearly , we see that for different values of @xmath5 , the cooperation level first increases dramatically from zero until reaching the maximum value at a moderate @xmath3 , then decreases slowly with increasing @xmath3 . here , we do not show the cooperation level in dependence on @xmath3 for large @xmath5 . in fact , we can still observe the nonmonotonous dependence of @xmath4 on @xmath3 even for large @xmath5 . these results suggest that the pgg model with the adaptive and bounded multiplication factor can effectively enhance cooperation in spatial structures , and higher cooperation level can emerges if an appropriate limitation is considered for the dynamical multiplication factor .    in order to intuitively understand the evolution of cooperation , we show some typical snapshots of the distribution of strategy and multiplication factor in the whole population in fig .  [ fig2 ] . we find that at the beginning of evolution , cooperators can form many small and isolated patches . but subsequently , some small compact cooperator clusters are embedded in the sea of defectors . as time increases , the cooperator clusters increase gradually , and finally cooperators may expand as a single ever growing cluster [ upper row in fig . [ fig2 ] ] . correspondingly , the multiplication factor in the full cooperation group can reach the upper limit , whereas the multiplication factor in the full defection group can reach the lower limit [ bottom row in fig . [ fig2 ] ] . however , the multiplication factor within the groups along the boundary of cooperators and defectors reaches a value between the lower and upper limit .    in combination with the above investigations , let us now explain the emergent results . indeed , a feedback mechanism is at work between the strategy distribution and the distribution of multiplication factors in all the groups . therefore , cooperators form compact clusters , and these clusters can become larger and larger , especially when the global cooperation level is not very high . meanwhile , the multiplication factor in the cooperators clusters becomes larger and larger , which provides cooperators , especially the ones on the boundary , with a higher payoff . whereas defectors also gather together , and the multiplication factor in the defectors clusters becomes smaller and smaller [ bottom row in fig . [ fig2 ] ] . in those interacting groups where the multiplication factor is negative , it is better for the players to choose the defective strategy such that they can have relatively higher payoffs . in a sense , this adaptive mode can induce a double - edged sword effect on the evolution of cooperation . however , under the social learning defectors are inclined to learn from their neighboring cooperators . as a consequence , the evolution of cooperation can be favored by this locally adaptive investment return .    if there is no limitation for the adaptive multiplication factor , due to the continuing negative feedback effects cooperators can not invade the defectors clusters , even if the feedback strength is high . in this situation , cooperators can thrive , but can not dominate the whole population . when the limitation is considered for the dynamical multiplication factor , the feedback mechanism , especially the negative feedback effects , can be effectively weakened . if there is too much restriction for the adaptive multiplication factor , that is , @xmath3 is not very large , e.g. , @xmath7 , the multiplication factor in some groups can reach the upper limit from a negative value due to the social learning , but this upper limit can not warrant the promotion of cooperation in spatial pgg even if the feedback strength is enough high . whereas for a larger @xmath3 , e.g. , @xmath11 , the multiplication factor in the group along the boundary reaching the upper limit can warrant a better promotion of cooperation . hence , this adaptive and bounded mode for the multiplication factor can provide a better environment for the evolution of cooperation . however , under this adaptive and bounded mode , the multiplication factor in some groups along the boundary can not rapidly become a positive and large value from a negative one when the feedback strength is not very high [ bottom row in fig . [ fig2 ] ] . although the average multiplication factor @xmath12 in the whole population can reach an enough high value which can make the cooperation level reach one in the traditionally spatial pgg @xcite , the average multiplication factor along the boundary of cooperators and defectors @xmath13 becomes negative as time increases [ fig . [ fig3](a ) ] . this does not provide a favorable environment for players interactions . correspondingly , the average payoffs of cooperators and defectors along the boundary are both negative . moreover , as time increases the average payoff of cooperators along the boundary is just a litter higher than the one of defectors along the boundary , but have larger fluctuations [ fig .  [ fig3](b ) ] . under the stochastic strategy updating , defectors do not always successfully imitate their neighboring cooperators , but sometimes may spread their strategy to the cooperators . as a result , cooperators can not defeat those defectors along the boundary , and they can only coexist with defectors for an exceedingly long time . on the contrary , when the feedback strength is high , by means of social learning the multiplication factor in the group along the boundary can have the opportunity to suddenly reach the upper limit , which can warrant the invading of cooperative behavior into the defectors clusters . finally , full cooperation can be achieved .    in what follows , we study how cooperation evolves if the lower limit is not the opposite number of the upper limit . figure  [ fig4](a ) shows the typical time evolution of cooperation for fixed @xmath14 and three different values of @xmath1 . we find that increasing the value of @xmath1 can make the system reach a higher cooperation level , but the cooperation level at equilibrium for @xmath15 is just slightly larger than the one for @xmath16 . in fig . [ fig4](b ) we show the fraction of cooperators as a function of time for fixed @xmath16 and three different values of @xmath0 . we find that increasing the value of @xmath0 can make the system reach a higher cooperation level . moreover , as time increases the fraction of cooperators first drops and then rapidly increases , but the larger values of @xmath1 or @xmath0 make the cooperation level increases faster . we also find that to have a favorable cooperation level , it is better to set the lower limit higher and it is not necessary to set the upper limit too high . finally , we study whether cooperation can be better promoted if the adaptive multiplication factor is constrained between @xmath17 and @xmath18 by means of a comparative investigation . previous work has reported that in the traditional pgg where only one invariable multiplication factor exists in the whole population , defectors outperform cooperators in any given mixed group for @xmath19 @xcite . we further find that in the traditionally spatial pgg , for noise value @xmath20 cooperators can dominate the whole population only if @xmath21 , and they can survive in the system only if @xmath22 , as shown in fig . [ fig5](a ) . it is worth pointing out that the traditionally spatial pgg corresponds to the situation of @xmath23 in this present model , where the multiplication factor in each interacting group is fixed at @xmath24 . in fig . [ fig5](b ) , we set @xmath25 , and show the cooperation level as a function of feedback strength @xmath26 for @xmath27 and different values of @xmath28 . we see that the cooperator density varying with @xmath5 displays two different behaviors : for smaller values of @xmath1 , e.g. , @xmath29 , the cooperation level first decreases and then increases until reaching the maximum value . subsequently , it decreases very slowly with increasing @xmath5 and its value approaches @xmath9 ; for larger values of @xmath1 , e.g. , @xmath30 , the cooperation level does not change too much for small values of @xmath5 , then monotonously increases to one with increasing @xmath5 . in addition , for smaller values of @xmath1 just a small amount of @xmath5 ( @xmath31 ) is needed to warrant a better promotion of cooperation in comparison to the traditionally spatial pgg , and the critical amount of @xmath5 becomes smaller if the @xmath1 is increased . for larger values of @xmath1 , the cooperation level for any value of @xmath26 is not less than the one for @xmath23 . in fact , the average multiplication factor in the population for @xmath26 is not larger than the one for @xmath23 , but these results suggest that cooperation can be better promoted in comparison to the traditionally spatial model . in addition , [ fig5](b ) shows that the cooperation level increases with increasing the upper limit of the multiplication factor and the initial value of the multiplication factor . we have also verified that increasing the initial fraction of cooperators is beneficial for the evolution of cooperation in this model . in summary , we have presented a coevolutionary rule where the multiplication factor in each interacting group is updated based on the local strategy distribution in the group and the global strategy distribution in the whole population , and studied its impact on the evolution of cooperation in the spatial public goods game . we found that this adaptive rule for the multiplication factor can effectively enhance the evolution of cooperation . when the appropriate bounded limitation for the dynamical multiplication factor is further considered , cooperation can be better promoted . in particular , full cooperation can be achieved in the system when the feedback strength is high enough . also , increasing the lower and upper limit values of the multiplication factor is favorable for the evolution of cooperation , but high cooperation level can be reached even if the upper limit is not very large . we further found that even if the multiplication factor is constrained to change between one and the group size , cooperation can be better promoted in the adaptive mode , in comparison to the classically spatial public goods game where the payoff parameter in each group is fixed and identical . the adaptive mode for the investment returns results in that a feedback mechanism is at work , that is , the matthew effect is introduced . from the viewpoint of this emergent feature , our model is related to the one proposed by perc @xcite , who considered that the reproductive success of each individual is updated by means of the enforcement of strategy and the distribution of public goods is driven by the reproductive success of individuals . under this success - driven mechanism , cooperation can be promoted . however , defectors can have a much higher payoff even in the sea of cooperators , and easily enforce their strategy choice to their neighbors . correspondingly , the superpersistent defector emergences spontaneously , and cooperators can not dominate the whole population . the complete dominance of cooperators is elusive even if the limitation factor about the value of reproductive success is considered . whereas in our model , the multiplication factor in each interacting group is updated based on the local and global strategy distribution , which characterizes the local investment environment for collective interactions , rather than individual s personality . also , we incorporate the limitation factor for the dynamical multiplication factor . in this framework , cooperators and defectors can form their own compact clusters respectively , and correspondingly cooperators along the boundary can have a higher payoff than the neighboring defectors . under the social learning , cooperators can easily spread their strategy even if the noise level for strategy updating is large . in particular , when the feedback strength is high , the interacting environment including some defectors can rapidly become favorable . thus , cooperators can gradually invade defector s clusters , and finally dominate the whole population . it could be concluded that our work further enriches the knowledgeless of coevolutionary rules in pggs , and importantly our spatial pgg model with adaptive and bounded investment returns not only can promote cooperation , but also make cooperators completely dominate the population . it is worth emphasizing the bounded values of the multiplication factor play a different role in the evolution of cooperation in our model in comparison to the one in ref . it is found that the main results remain qualitatively if the limit value of the payoff parameter in the prisoner s dilemma is large enough . moreover , it is demonstrated that the final cooperation level strongly depends on the values of initial payoff parameter and feedback strength in ref . in particular , with increasing the initial payoff parameter , the probability that the system ends in full cooperation state decreases . whereas in this work , we find that the introduced limitation factor can weaken the matthew effect , particularly the negative feedback effect on the unfavorable interacting groups , which makes the limit value a crucial model parameter . the limited negative effects can be overcome via social learning . hence , appropriate limitation can warrant the best promotion of cooperation . in addition , we show that the finial cooperation level not only depends on the the values of initial payoff parameter and feedback strength , but also depends on the limit values . with increasing the initial payoff parameter and the initial fraction of cooperators , cooperation can be better enhanced when the dynamical multiplication factor is limited . in a sense , this work further explores the effects of adaptive and bounded game payoffs on the evolution of cooperation .    during the coevolutionary process , the investment return in most of interacting groups reaches the upper or lower limit due to the feedback effects . this segregation and polarization of investment returns occurs spontaneously over time , which is different from the distribution in ref . @xcite . in the latter case , the distribution of the multiplication factor is artificially introduced by the authors and does not change during the evolutionary process . although the emergent values of the multiplication factor in all the interacting groups do not display too much diversity , we find that cooperation can be promoted in this adaptive and bounded mode . compared with the results in the traditionally spatial pgg model , cooperation can be better promoted even if the multiplication factor can only change between one and the group size .    in the present model , we consider the adaptive mode for the multiplication factor in a group based on the local cooperation level in the classical pgg where players just have two discrete strategy choices @xmath32 or @xmath33 , and correspondingly the local cooperation level only has several finite values . to make the local cooperation level change continuously between zero and one , we also introduce the adaptive and bounded investment returns into the spatial continuous pgg @xcite , and still find that this pgg with adaptive and bounded investment returns promotes cooperation . we also test our model in well - mixed populations as well as on other types of interactions networks , and still find that cooperation can be enhanced by the proposed coevolutionary rule . moreover , we would like to point out that in this work we fix the value of noise to one . in general , the qualitative behavior of the system remains unchanged for other values of noise , although for pairwise interactions there may exist an optimal value of noise at which the evolution of cooperation is most successful @xcite . it could be inferred that if we can choose the optimal noise value for strategy updating , the positive effects from social learning can be amplified . we also believe that , cooperation can be better promoted if we further incorporate the selection of noise level in strategy adoption @xcite into this adaptive and bounded mode for the multiplication factor . we consider the pgg on a square lattice of size @xmath34 with periodic boundary conditions . each individual who is a pure strategist can only follow two simple strategies : cooperate ( @xmath32 ) and defect ( @xmath33 ) . cooperators contribute a fixed amount ( here considered to be equal to @xmath17 without loss of generality ) to the public good while defectors contribute nothing . the sum of all contributions in each group @xmath35 is multiplied by the factor @xmath36 , and the resulting public goods are distributed among all the group members . correspondingly , the payoff of player @xmath37 from the group @xmath35 is    @xmath38    where @xmath39 denotes the strategy of player @xmath37 , @xmath40 denotes the number of cooperators in the group @xmath35 , and @xmath18 denotes the group size . here , we consider the square lattice with von neumann neighborhood . accordingly , the interacting group size is fixed at @xmath41 , and each individual belongs to five different groups . the payoff of each player is accumulated from the fixed five interacting groups , and thus player @xmath37 total payoff @xmath42 .    after playing the games , the multiplication factor in each group needs to be updated . specifically , we assume that the multiplication factor of the group centered on player @xmath37 at time @xmath43 is @xmath44,\ ] ] where @xmath5 controls the strength of feedback from the comparison between the local and global cooperative environments , @xmath45 is the multiplication factor of the focal individual @xmath37 s group at time @xmath46 , @xmath47 is the fraction of cooperators in the whole population at time @xmath46 , and @xmath48 is the local cooperation level in the group centered on @xmath37 at time @xmath46 . here , @xmath49 , where @xmath50 denotes the number of cooperators in the group where player @xmath37 is the focal individual at time @xmath46 . moreover , we consider the limitation for the adaptive multiplication factor following previous work @xcite , that is , letting @xmath51 if @xmath52 and letting @xmath53 if @xmath54 . here , @xmath1 and @xmath0 respectively represent the upper and lower limits of the multiplication factor in each group , and they satisfy the following inequalities : @xmath55 and @xmath56 . in particular , when @xmath57 , the multiplication factor is constrained in a symmetric interval , which is the same to the setting in ref.@xcite . moreover , when @xmath58 ( no bounded limitation for the multiplication factor ) or @xmath23 ( no updating for the multiplication factor ) , the average multiplication factor in all the interacting groups @xmath59 , where @xmath60 is the initial value of the multiplication factor in each group . subsequently , each player is allowed to learn from one of its neighbors and update its strategy . player @xmath37 adopts the randomly chosen neighbor @xmath61 strategy with a probability depending on the payoff difference as @xmath62},\ ] ] where @xmath63 denotes the amplitude of noise @xcite , accounting for imperfect information and errors in decision making . following previous work @xcite , we simply set @xmath20 representing that it is very likely that the better performing players will pass their strategy to other players , yet it is possible that players will occasionally learn also from the less successful neighbors . simulations of this spatial pgg model are performed by means of a synchronous updating rule , using @xmath64 to @xmath65 system size . initially , the two strategies of @xmath32 and @xmath33 are randomly distributed among the population with an equal probability , and the multiplication factor in each interacting group has the same value @xmath60 . the key quantity for characterizing the cooperative behavior of the system is the density of cooperators , which is defined as the fraction of cooperators in the whole population . the system can reach a dynamical equilibrium after a suitable transient time @xcite . then the density of cooperators reaches its asymptotic value @xmath4 and remains there within small fluctuations ( less than @xmath66 ) . this asymptotic value is taken to describe the cooperation level in the whole population , and all the simulation results are averaged over @xmath67 different realizations of initial conditions . yl acknowledges the financial support from the special fund for the basic scientific research of central colleges , changan university , the special fund for the basic research support planning of changan university , and the open fund of road traffic intelligent detection and equipment engineering research center , shaanxi province ( grant no . chd2010jc134 ) . yz acknowledges the financial support from the national natural science foundation of china ( nsfc ) ( grand no . 11161011 ) . lw acknowledges the financial support from 973 program ( 2012cb821203 ) and the national natural science foundation of china ( nsfc ) ( grant nos . 61020106005 and 10972002 ) . mp acknowledges the financial support from the slovenian research agency ( arrs ) ( grant j1 - 4055 ) . the funders had no role in study design , data collection and analysis , decision to publish , or preparation of the manuscript .                                                                                in dependence on the feedback strength @xmath5 for different values of @xmath3 . panel ( b ) depicts the fraction of cooperators in dependence on the boundary value @xmath3 for different values of @xmath5 . it can be observed that cooperation can be promoted for large values of feedback strength , and there exist moderate boundary values warranting the best promotion of cooperation . here , @xmath69.,width=480 ]         and different values of upper limit . panel ( b ) depicts the fraction of cooperators in the whole population as a function of time for fixed upper limit @xmath16 and different values of lower limit . increasing the values of lower and upper limit can provide more positive effects on the evolution of cooperation . here , @xmath70 and @xmath69.,width=480 ]    $ ] . * panel ( a ) shows the fraction of cooperators as a function of @xmath60 for @xmath23 . in this situation , the model recovers to the traditionally spatial pgg , where the multiplication factor in each group is fixed at @xmath60 and @xmath71 . for @xmath20 , cooperators can survive only if @xmath72 , and they can dominate the whole population only if @xmath73 . panel ( b ) shows the fraction of cooperators as a function of @xmath5 for fixed @xmath27 and different values of @xmath1 . initially , the multiplication factor in each interacting group is @xmath25 . dash lines are used to indicate the critical value of @xmath26 for a better promotion of cooperation in this adaptive and bounded mode for the enhancement factor.,width=480 ]
Abstract:

Answer: the public goods game is one of the most famous models for studying the evolution of cooperation in sizable groups . the multiplication factor in this game can characterize the investment return from the public good , which may be variable depending on the interactive environment in realistic situations . instead of using the same universal value , here we consider that the multiplication factor in each group is updated based on the differences between the local and global interactive environments in the spatial public goods game , but meanwhile limited to within a certain range . we find that the adaptive and bounded investment returns can significantly promote cooperation . in particular , full cooperation can be achieved for high feedback strength when appropriate limitation is set for the investment return . also , we show that the fraction of cooperators in the whole population can become larger if the lower and upper limits of the multiplication factor are increased . furthermore , in comparison to the traditionally spatial public goods game where the multiplication factor in each group is identical and fixed , we find that cooperation can be better promoted if the multiplication factor is constrained to adjust between one and the group size in our model . our results highlight the importance of the locally adaptive and bounded investment returns for the emergence and dominance of cooperative behavior in structured populations .


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: it is indicated by seven year wilkinson microwave anisotropy probe data combined with measurements of baryon acoustic oscillations and hubble constant that @xmath4@xmath5 of the matter content in the universe is non - baryonic dark matter @xcite . weakly interacting massive particles ( wimps ) @xcite , predicted by extensions of the standard model of particle physics , are a well - motivated class of candidates for dark matter . they are distributed in the halo surrounding the milky way . wimps may be directly detected through measuring nuclear recoils in terrestrial detectors produced by their scattering off target nuclei @xcite . the nuclear reoils is expected to have a roughly exponential energy distribution with a mean energy in a few tens of kev @xcite .    in direct searches for wimps , there are three different methods used to detect the nuclear recoils , including collecting ionization , scintillation and heat signatures induced by them . the background of this detection is made up of electron recoils produced by @xmath6 and @xmath7 scattering off electrons , and nuclear recoils produced by neutrons scattering elastically off target nucleus . it is very efficient discriminating nuclear recoils from electron recoils with pulse shape discrimination , hybrid measurements and so on . the rejection powers of these techniques can even reach @xmath8@xmath9 @xcite . for example , the cdms - ii @xcite and edelweiss - ii @xcite experiments measure both ionization and heat signatures using cryogenic germanium detectors in order to discriminate between nuclear and electron recoils , and the xenon100 @xcite and zeplin - iii @xcite experiments measure both ionization and scintillation signatures using two - phase xenon detectors . however , it is very difficult to discriminate between nuclear recoils induced by wimps and by neutrons . in some dark matter experiments and researches , tagging neutron is applied to reject neutron background . in the zeplin - iii experiment , the @xmath10 gadolinium ( gd ) doped polypropylene is used as the neutron veto device , and its maximum tagging efficiency for neutrons reaches about @xmath11 @xcite . in ref . @xcite , the @xmath12 gd - doped water is used as the neutron veto , and its neutron background can be reduced to 2.2 ( 1 ) events per year per tonne of liquid xenon ( liquid argon ) . in our past work @xcite , the reactor neutrino detector with @xmath13 gd - doped liquid scintillator ( gd - ls ) is used as the neutron veto system , and its neutron background can be reduced to about 0.3 per year per tonne of liquid xenon . cesium iodide(csi ) crystals as a kind of dark matter target have been applied to dark matter experiments , such as the kims experiment@xcite . in a dark matter experiment with csi target , csi crystals need nt to be cooled with the liquid nitrogen or xenon , so this detector can be of the simpler structure . in the sun , lu et al.s work@xcite , the rejection power against electron recoil can reach o(@xmath14 ) with na - doped csi crystals . the feasibility of direct wimps detection with a neutron veto based on a neutrino detector had been validated in our past work@xcite . so , in the present paper , a neutrino detector with gd - ls ( @xmath13 gd - doped ) is still used as a neutron - tagged device and wimp detectors with csi target are placed inside the gd - ls . here we designed an experimental configuration : four wimp detectors with csi(na ) target are individually placed inside four reactor neutrino detector modules which are used as a neutron veto system . the experimental hall of the configuration is assumed to be located in an underground laboratory with a depth of 910 meter water equivalent ( m.w.e . ) , which is similar to the far hall in the daya bay reactor neutrino experiment@xcite . collecting scintillation signals is considered as the only method of wimps detection in our work . the neutron background for this design are estimated using the geant4 @xcite simulation . four identical wimp detectors with csi(na ) target are individually placed inside four identical neutrino detector modules . the concentration of na is 0.02% in csi crystals . the experimental hall of this experimental configuration is assumed to be located in an underground laboratory with a depth of 910 m.w.e . , which is similar to the far hall in the daya bay reactor neutrino experiment . the detector is located in a cavern of 20@xmath020@xmath020 @xmath15 . the four identical cylindrical neutrino modules ( each 413.6 cm high and 393.6 cm in diameter ) are immersed into a 13@xmath013@xmath08 @xmath15 water pool at a depth of 2.5 meters from the top of the pool and at a distance of 2.5 meters from each vertical surface of the pool . the detector configuration is shown in fig.[fig : detector ] . each neutrino detector module is partitioned into three enclosed zones . the innermost zone is filled with gd - ls , which is surrounded by a zone filled with unload liquid scintillator ( ls ) . the outermost zone is filled with transparent mineral oil . 366 8-inch pmts are mounted in the mineral oil . each wimp detector ( 69.6 cm height , 59.9 cm in diameter ) consists of three components : a csi(na ) crystal array , pmts and copper vessels . the csi(na ) array is made up of 31 csi(na ) crystals whose sections are regular hexagons ( 34 cm height , 4.5 cm side length ) and is placed inside a 1 cm thick copper vessel ( 69.2 cm height , 59.5 cm in diameter ) filled with dry nitrogen gas . fig.[fig : detector ] shows the arrangement of these crystals . 62 3-inch pmts are individually mounted on two ends of these crystals . the copper vessel of each wimp detector is surrounded by an 0.2 cm thick aluminum reflector for photons produced in the gd - ls . the geant4 ( version 8.2 ) package@xcite has been used in our simulations . the physics list in the simulations includes transportation processes , decay processes , low energy processes , electromagnetic interactions ( multiple scattering processes , ionization processes , scintillation processes , optical processes , cherenkov processes , bremsstrahlung processes , etc . ) and hadronic interactions ( lepton nuclear processes , fission processes , elastic scattering processes , inelastic scattering processes , capture processes , etc . ) . the hadronic processes include the low energy(<20 gev ) , high energy(>20gev ) and neutron high - precision(<20 mev ) models . the cuts for the productions of gammas , electrons and positrons are 1 mm , 100 @xmath16 m and 100@xmath16 m , respectively . the quenching factor is defined as the ratio of the detector response to nuclear and electron recoils . the birks factor for protons in the gd - ls is set to 0.01 g/@xmath17/mev , corresponding to the quenching factor 0.17 at 1 mev , in our simulation . besides , we utilize 10 intel core i5 cpus ( four cores , each core offer a base of speed of 2.8ghz ) in the neutron background evaluation . the recoil energies for wimp interactions with csi(na ) nuclei were set to a range from 20 kev to 100 kev in this work . proton recoils induced by neutrons and neutron - captured signals are used to tag neutrons which reach the gd - ls . the energy deposition produced by proton recoils is close to a uniform distribution . neutrons captured on gd and h lead to a release of about 8 mev and 2.2 mev of @xmath6 particles , respectively . due to the instrumental limitations of the gd - ls , we assume neutrons will be tagged if their energy deposition in the gd - ls is more than 1 mev , corresponding to 0.17 mevee ( electron equivalent energy ) . in the gd - ls , it is difficult to distinguish signals induced by neutrons from electron recoils , which are caused by the radioactivities in the detector components and the surrounding rocks . but these radioactivities can be controlled to less than @xmath4 50 hz according to the daya bay experiment@xcite . if we assume a 100 @xmath16s for neutron tagging time window , the indistinguishable signals due to the radioactivities will result in a total dead time of less than 44 hours per year . neutrons are produced from the detector components and their surrounding rock . for the neutrons from the surrounding rock there are two origins : first by spontaneous fission and ( @xmath18 , n ) reactions due to @xmath19 and @xmath20 in the rock ( these neutrons can be omitted because they are efficiently shielded , see sec.4.2 ) , and secondly by cosmic muon interactions with the surrounding rock . we estimated the number of neutron background in the csi(na ) target of one tonne . this number has been normalized to one year of data taking and are summarized in tab.[tab : bg ] . neutrons from the detector components are induced by ( @xmath18 , n ) reactions due to u and th . according to mei et al.@xcite , the differential spectra of neutron yield can be expressed as    @xmath21@xmath22    where @xmath23 is the total number of atoms for the @xmath24 element in the host material , @xmath25@xmath26 refers to the @xmath18-particle production rate for the decay with the energy @xmath27 from @xmath28 or @xmath29 decay chain , @xmath30 refers to the @xmath18 energy , @xmath31 refers to the neutron energy , and @xmath32 is the mass stopping power of the @xmath24 element . the u and th contaminations in the @xmath33 material are considered as the only neutron source in the pmts in our work . neutrons from @xmath33 are emitted with their average energy of 2.68 mev@xcite . the pmts in the copper vessels of the four wimp detectors amounts to 248 . the u and th concentrations in the pmt components can reach ten or even less ppb@xcite , so a rate of one neutron emitted per pmt per year is conservatively estimated@xcite . consequently , there are 248 neutrons produced by all the pmts in the copper vessels per year . a simulated sample of 2.48@xmath0@xmath9 events is used to study this neutron background . these events generated isotropically are uniformly distributed in the @xmath33 material of the pmts . the simulation result is summarized in tab . 1 . 2.9 neutron events/(ton@xmath34yr ) reach the csi(na ) target and their energy deposition falls in the same range as that of the wimp interactions , as seen in tab . 1 . because 0.04 of them are not tagged in the gd - ls , these background events can not be eliminated . the uncertainties of the neutron background from the pmts in tab.[tab : bg ] are from the binned neutron spectra in ref.@xcite . but the neutron background errors from the statistical fluctuation ( their relative errors are less than 1@xmath35 ) are too small to be taken into account .      in the copper vessels , neutrons are produced by the u and th contaminations and emitted with their average energy of 0.81 mev@xcite . their total volume is about @xmath36 @xmath37 . the radioactive impurities th can be reduced to @xmath38 ppb in some copper samples@xcite . if we conservatively assume a 0.001 ppb u / th concentrations in the copper material@xcite , a rate of one neutron emitted per @xmath39 @xmath37 per year is estimated@xcite . consequently , there are 1.8 neutrons produced by the all copper vessels per year . a simulated sample of 1.8@xmath0@xmath40 events is used to study this neutron background . these events generated isotropically are uniformly distributed in the copper material of the copper vessels . the simulation result is summarized in tab.[tab : bg ] . 0.03 neutron events/(ton@xmath34yr ) reach the csi(na ) target and their energy deposition falls in the same range as that of the wimp interactions ( see tab.[tab : bg ] ) . as 0.001 of them are not tagged in the gd - ls , these background events can not be eliminated . the uncertainty of the neutron background from the copper vessels in tab.[tab : bg ] are from the binned neutron spectra in the ref.@xcite . but the neutron background errors from the statistical fluctuation ( their relative errors are less than 1@xmath35 ) are too small to be taken into account .      the u and th contaminations in other detector components also contribute to the neutron background in our experiment setup . neutrons from the aluminum reflectors are emitted with the average energy of 1.96 mev@xcite . the u and th contaminations in the carbon material are considered as the only neutron source in the gd - ls / ls . neutrons from the gd - ls / ls are emitted with the average energy of 5.23 mev@xcite . the u and th contaminations in the @xmath33 material are considered as the only neutron source in the pmts in the oil . neutrons from pmts are emitted with the average energy of 2.68 mev@xcite . the u and th contaminations in the iron material are considered as the only neutron source in the stainless steel tanks . neutrons from the stainless steel tanks are emitted with the average energy of 1.55 mev@xcite . we evaluated the neutron background from the above components using the geant4 simulation . all the nuclear recoils in the csi(na ) target , which fall in the same range as that of the wimp interactions , are tagged . the neutron background from these components can be ignored . the simulated samples that amount to 1000 years of data taking are used to evaluate these neutron backgrounds .      in the surrounding rock , almost all the neutrons due to natural radioactivity are below 10 mev @xcite . water can be used for shielding neutrons effectively , especially in the low energy range of less than 10 mev @xcite . the wimp detectors are surrounded by about 2.5 meters of water and more than 1 meter of gd - ls / ls , so these shields can reduce the neutron contamination from the radioactivities to a negligible level . neutrons produced by cosmic muon interactions constitute an important background component for dark matter searches . these neutrons with a hard energy spectrum extending to several gev energies , are able to travel far from produced vertices .    the total cosmogenic neutron flux at a depth of 910 m.w.e . is evaluated by a function of the depth for a site with a flat rock overburden @xcite , and it is 1.31@xmath0@xmath41 @xmath42 . the energy spectrum ( see fig.[fig : energy ] ) and angular distribution of these neutrons are evaluated at the depth of 910 m.w.e . by the method in @xcite . the neutrons with the specified energy and angular distributions are sampled on the surface of the cavern , and the neutron interactions with the detector are simulated with the geant4 package . a simulated sample of 1.25@xmath0@xmath43 events is used to study in this neutron background . tab.[tab : bg ] shows that 2.5 neutron events/(ton@xmath34yr ) reach the csi(na ) target and their energy deposition falls in the same range as that of the wimp interactions . 0.28 of them are not tagged by the gd - ls / ls . muon veto systems can tag muons very effectively , thereby most cosmogenic neutrons can be rejected . in the daya bay experiment , the contamination level can even be reduced by a factor of more than 30@xcite . we assume the neutron contamination level from cosmic muons decreases by a factor of 30 using a muon veto system . this could lead to the decrease of cosmogenic neutron contamination to 0.01 events/(ton@xmath34yr ) . the uncertainties of the cosmogenic neutron background in tab.[tab : bg ] are from the statistical fluctuation . besides neutron background , other background events are mainly from reactor neutrino events and electron recoils in the experimental design in the present paper . the contamination caused by electron recoils consists of bulk electron recoil events and surface events . since neutrino detectors are fairly close to nuclear reactors ( about 2 kilometers away ) in reactor neutrino experiments , a large number of reactor neutrinos will pass through the detectors , and nuclear recoils will be produced by neutrino elastic scattering off target nucleus in the wimp detectors . although neutrinos may be a source of background for dark matter searches , they can be reduced to a negligible level by setting the recoil energy threshold of 10 kev@xcite . besides , nuclear recoils may also be produced by low energy neutrons produced by the inverse @xmath7-decay reaction @xmath44 . but their kinetic energies are almost below 100 kev@xcite , and their maximum energy deposition in the wimp detectors is as large as a few kev . thus the neutron contamination can be reduced to a negligible level by the energy threshold of 10 kev . the intenal sources of electron recoils are mainly caused by the radioisotopes of @xmath45 and @xmath46 in the csi crystals , and the external ones are mainly from the radioisotopes of @xmath19 , @xmath20 and @xmath47 in the 3 " pmts inside the copper vessels . here we assume that the concentrations of the @xmath45 and @xmath46 in the csi crystals are about 2 mbq / kg and 1 ppb@xcite , respectively . then their background rates will be about 1.7 counts / kev / kg / day ( cpd ) in the region 10kev@xcite . while the concentrations of u / th / k are about 78 , 25 and 504 mbq / pmt , respectively@xcite , and their background rates will be about 0.5 cpd in the region 10kev . considering the fact that the rejection power against electron recoils can reach o(@xmath14)@xcite , we roughly estimate that the bulk electron recoil contamination due to the internal and external sources of the csi(na ) crystals is about 0.9 events/(ton@xmath34yr ) .      the surface events for csi(na ) crystals are caused by the deliquescence which reduces the @xmath48 concentration on the crystal surface@xcite . signals produced by @xmath18 particles from the progenies of @xmath49 in the air on the `` old '' surface of the crystal mimic nuclear recoils seriously . the surface events can be prevented from by avoiding the deliquescence of the crystal surface . for example , the copper vessels are filled with dry nitrogen gas@xcite , and thus there is no air between the copper vessels and csi(na ) crystals . hence the surface event contamination can be reduced to a negligible level by the above methods . neutron background can be effectively suppressed neutrino detectors used as a neutron veto system in direct dark matter searches . tab.[tab : bg ] shows the total neutron contamination are 0.05 events/(ton@xmath34yr ) . and compared to ref.@xcite , it is reduced by a factor of about 6 . if the electron recoils contaminations are considered , the total background amounts to about 1 events/(ton@xmath34yr ) .    according to our work , the neutron background is mainly from the pmts in the copper vessels in this configuration with the csi(na ) targets . after finishing a precision measurement of the neutrino mixing angle @xmath50 , we can utilize the existing experiment hall and neutrino detectors . this will not only save substantial cost and time for direct dark matter searches , but the neutron background could also decrease to o(0.01 ) events per year per tonne of csi(na ) in the case of the daya bay experiment . according to ref.@xcite , the neutron fluxes in the reno ( in an underground laboratory with a depth of 450 m.w.e . ) , double chooz ( in an underground laboratory with a depth of 300 m.w.e . ) experiments@xcite are respectively about 5 and 3 times more than that of the daya bay experiment . their neutron backgrounds are roughly estimated to be about 0.1 events/(ton@xmath34yr ) , if their detector configurations are the same as the one described above . the neutron contamination is one order of magnitude smaller than the electron recoil contamination , so neutron contamination can be ignored in this detector configuration .    to evaluate the detector capability of directly detecting dark matter , we assume a standard dark matter galactic halo@xcite , an energy resolution that amounts to 25@xmath35 for the energy range of interest and 40@xmath35 nuclear recoil acceptance@xcite . if no signals are significantly observed , a sensitivity to wimp - nucleon spin - independent elastic scattering can be calculated via the same method as ref.@xcite . our calculation shows that an exposure of one tonne @xmath0 year could reach a cross - section of about 3@xmath0@xmath51 at the 90@xmath35 confidence level ( see fig.[fig : limit ] ) . this work is supported in part by the national natural science foundation of china ( nsfc ) under contract no . 11235006 and the science fund of fujian university of technology under contract no . gy - z13114 . n. jarosik , et.al . , apjs , 192 , 14 ( 2011 ) w. j. percival , et al . , mnras , 401 , 2148 ( 2010 ) a. g. riess , et al . , apj , 699 , 539 ( 2009 ) g. steigman and m. s. turner , nucl . b253 , 375 ( 1985 ) m. w. goodman and e. witten , phys . d31 , 3059 ( 1985 ) g. jungman , m. kamionkowski , and k. griest , phys . 267 , 195 ( 1996 ) r. j. gaitskell , ann . nucl . part . 54 , 315 ( 2004 ) g. bertone , d. hooper and j. silk , phys . 405 , 279 ( 2005 ) j. d. lewin and p. f. smith , astropart . phys . 6 , 87 ( 1996 ) z. ahmed , et al . , science 327 , 1619 ( 2010 ) w.h.lippincott , et al . , c78 , 035801 ( 2008 ) e. armengaud , et al . b702 , 329 ( 2011 ) e. aprile , et al . , astropart . 35 , 573 ( 2012 ) d. akimov , et al . b709 , 14 ( 2012 ) d. akimov , et al . , astropart , phys . 34 , 151 ( 2010 ) a. bueno , m. c. carmona and a. j. melgarejo , jcap . 08 , 019 ( 2008 ) . ye xu , et al . , jcap06 , 009 ( 2011 ) s.c . kim , et al . 108 , 181301 ( 2012 ) xilei sun , junguang lu , et al . instrum . meth . a642 , 52 - 58 , 2011 f.p.an , et al . letts , 108 , 171803 ( 2012 ) ; daya bay collaboration , arxiv : hep - ex/0701029v1 s. agostinelli , et al . , instru . meth . a506 , 250 ( 2003 ) d. m. mei , c. zhang and a. hime , nucl . instrum . meth . a606 , 651 ( 2009 ) m.j.carson et al.,astropart . 21 667 ( 2004 ) a. bueno , m.c . carmona and a.j . melgarejo , jcap08 , 019 ( 2008 ) , arxiv : 0805.0694 marin e. keillor , et al . , j. radioanal . 282 , 703 - 708 ( 2009 ) d. s. leonard , et al . , nucl . instrum . and meth . a591 , 490 ( 2008 ) m. j. carson , et al . , astropart . 21 , 667 ( 2004 ) j. m. carmona , et al . , astropart . 21 , 523 ( 2004 ) d. m. mei and a. hime , phys . d73 , 053004 ( 2006 ) y.f.wang , et al . rev.d64 , 013012 ( 2001 ) jocelyn monroe and peter fisher , phys . d76 , 033007 ( 2007 ) m. apollonio , et al . , chooz collaboration , eur . j. c27 , 331 - 374(2003 ) , arxiv : hep - ex/0301017 t. y. kim , et al . instrum . and meth . a500 , 337 - 344 ( 2003 ) h. s. lee , et al . a571 , 644 - 655 ( 2007 ) j.k . ahn , et al . 108 , 191802 ( 2012 ) f. ardellier , et al . , double chooz : a search for the neutrino mixing angle @xmath50 , arxiv : hep - ex/0606025 g.j.feldman and r.d.cousins , phys . d57 , 3873 ( 1998 ) http://pisrv0.pit.physik.uni-tuebingen.de/darkmatter/limits/index.php    .estimation of neutron background from different sources for an underground laboratory at a depth of 910 m.w.e . the column labeled `` 20kev<@xmath52<100kev '' identifies the number of neutrons whose energy deposition in the csi(na ) target is in the same range as wimp interactions . the column labeled `` not tagged '' identifies the number of neutrons which are misidentified as wimp signatures ( their energy deposition in the csi(na ) target is in the same range as wimp interactions while their recoil energies in the gd - ls / ls are less than the energy threshold of 1 mev ) . the row labeled `` copper vessel '' identifies the number of neutrons from the copper vessels . the row labeled `` pmts in copper vessel '' identifies the number of neutrons from the pmts in the copper vessels . the row labeled `` cosmic muons '' identifies the number of cosmogenic neutrons in the case of not using the muon veto system . the row labeled `` muon veto '' identifies the number of cosmogenic neutrons in the case of using the muon veto system . we assume that neutron contamination level from cosmic muons decreases by a factor of 30 using a muon veto system . only the total background in the case of using the muon veto system is listed in this table . the terms after @xmath53 are errors . [ cols="^,^,^",options="header " , ]
Abstract:

Answer: a direct wimp ( weakly interacting massive particle ) detector with a neutron veto system is designed to better reject neutrons . an experimental configuration is studied in the present paper : a wimp detectors with csi(na ) target is placed inside a reactor neutrino detector . the neutrino detector is used as a neutron veto device . the neutron background for the experimental design has been estimated using the geant4 simulation . the results show that the neutron background can decrease to o(0.01 ) events per year per tonne of csi(na ) . we calculate the sensitivity to spin - independent wimp - nucleon elastic scattering . an exposure of one tonne @xmath0 year could reach a cross - section of about 3@xmath0@xmath1 pb . @xmath2department of mathematics and physics , fujian university of technology , fuzhou 350118 , china    @xmath3school of physics , nankai university , tianjin 300071 , china    dark matter , neutron background , neutrino detector , csi(na )    pacs numbers : 95.35.+d , 95.55.vj , 29.40.mc


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: high - mass stars are much less frequent than intermediate- or low - mass stars due to both the star formation process , which gives rise to an initial mass function declining with the mass ( e.g. * ? ? ? * ; * ? ? ? * ) , and to their shorter evolutionary times . however , high - mass stars are very important because they can affect their surroundings with their winds , their strong radiation fields , and their catastrophic death as supernovae , chemically enriching their environment and triggering star formation . they usually form within the dense cores of stellar clusters and/or associations where dynamical interactions play an important role . it is widely believed that massive stars may be the product of collisions between two or more intermediate - mass stars . this idea is supported by the fact that a large fraction of massive stars harbour close companions , as failed mergers . it has been recently estimated that at least 50% of massive stars are member of binary or multiple star system @xcite . this lucky occurrence allows to directly measure the masses by means of their radial velocity ( rv ) curves . in many cases , spectral lines of both components are visible ( sb2 systems ) , allowing to derive the orbital parameters like the period , @xmath5 , the projected semi - major axes , @xmath6 , and the masses , @xmath7 , apart from the factor @xmath8 . if our line - of - sight is close to the orbital plane and fractional radii of the components are not too small the stars display mutually eclipses . the orbital inclination and fractional radii of the component stars can be determined by the analysis of photometric light curves . therefore , eclipsing binaries are unique targets for determining the masses and radii from their combined light curves and radial velocities analyses . nevertheless , absolute radii were measured only for a rather small number of early - type b - stars which are members of eclipsing binary systems @xcite . thus , we started a systematic observing program devoted to the spectroscopic study of close eclipsing binary systems with at least one hot component . the light variability of v446cephei ( hd 210478 ; bd+60@xmath92348 ; hip 109311 ; v=7@xmath10.32 , b - v=0@xmath10.08 ) has been discovered by the @xmath0 photometry @xcite . it was classified as an algol - type eclipsing binary ( ea ) from visual inspection of light variation . the orbital period of v446cep has been determined by the @xmath0 satellite as 3.8084 days and the ephemeris is given as follows , @xmath11 the @xmath0 light curve ( lc ) displays a primary minimum with a depth of about 0.14 mag @xcite . @xcite designated it as v446cep and classified it as an ea according to the criteria of the general catalogue of variable stars . @xcite included the star in the list of probable members of the cepob2 association . in this paper we present new spectroscopic observations of v446cep . we evaluate effective temperatures and surface gravities of the component stars from an ad - hoc analysis of the spectra taken near the quadratures . combining the results obtained by the analysis of rv and light curves we determine directly absolute masses and radii of the components . we infer distance to the system , their fundamental parameters , and briefly discuss the main outcomes of this study . optical spectroscopic observations of v446cep were obtained at the tubitak national observatory using the turkish faint object spectrograph camera ( tfosc ) attached to the 1.5 m telescope . the observations were made from july 22 , 2012 to august 3 , 2013 , under good seeing conditions . further details on the telescope and the spectrograph can be found at http://www.tug.tubitak.gov.tr . the wavelength coverage of each spectrum was 4000 - 9000   in 11 orders , with a resolving power of @xmath12/@xmath13 @xmath147000 at 6500  . the average signal - to - noise ratio ( s / n ) was @xmath14120 . we also obtained high s / n spectra of early type standard stars 1cas ( b0.5iv ) , hr153 ( b2iv ) , @xmath15her ( b5iv ) , 21peg ( b9.5v ) and @xmath16lyr ( a0v ) which were used as templates in derivation of the radial velocities . the echelle spectra were extracted from the raw images following standard reduction steps involving electronic bias subtraction , flat field division , cosmic rays removal , optimal extraction of the echelle orders , and wavelength calibration thanks to the emission lines of a th - ar lamp . the reduction was performed using tasks of the iraf package our spectroscopic dataset contains 15 observations for v446cep . we have measured radial velocities ( rvs ) from the spectra , focusing on spectral segments containing the hei @xmath125876 ( order4 ) and @xmath126678 ( order3 ) lines which are the most prominent un - blended features in our spectra , apart from the balmer lines . we have employed the standard cross - correlation method for measuring the velocities of the component stars of the system . the cross - correlation technique @xcite is widely used for measuring rvs from the spectra of close binary systems . cross - correlation analyses were made using the spectra of @xmath15her and 21peg as templates . the principle spectral features showing splitting due to binarity were the hei lines at @xmath175876 and 6678 . we used also order 9 , containing the hei @xmath124471 line , for a few measurements of the radial velocities . the spectra taken close to the conjunctions , which display no double - lined feature , were disregarded . the balmer lines were not used in the measurements of radial velocities due to their extremely broad profiles . we obtained 15 radial velocities for each component of v446cep . the average radial velocities and their associated standard errors derived from the spectral segments containing hei @xmath174471 , 5876 , and 6678 lines are presented in table1 , along with the observation date and orbital phase . the mean error of radial velocities is 2.5kms@xmath18 for the primary , and 3.6 kms@xmath18 for the secondary star of v446cep . the rvs are plotted against the orbital phase in fig.1 , where the filled squares represent the primary and the empty squares the secondary stars . examination of the @xmath0 light curve show no evidence for any eccentricity in the orbit of the system . therefore , we have assumed circular orbit and analysed the rvs using the rvsim software programme @xcite . final orbital parameters are presented in table2 . intermediate - resolution optical spectroscopy permits us to derive most of the fundamental stellar parameters , such as the projected rotational velocity ( @xmath19 ) , spectral type ( s@xmath1 ) , luminosity class , effective temperature ( @xmath20 ) , surface gravity ( @xmath21 ) , and metallicity ( [ fe / h ] ) . the width of the cross - correlation function ( ccf ) is a good tool for the measurement of projected rotational velocity ( @xmath22 ) of a star . we use a method developed by @xcite to estimate the @xmath22 of each star composing the investigated sb2 system from its ccf peak by a proper calibration based on a spectrum of a narrow - lined star with a similar spectral type . the rotational velocities of the components were obtained by measuring the fwhm of the ccf peak related to each component in five high - s / n spectra acquired near the quadratures , where the spectral lines have the largest doppler - shift . the ccfs were used for the determination of @xmath19 through a calibration of the full - width at half maximum ( fwhm ) of the ccf peak as a function of the @xmath19 of artificially broadened spectra of slowly rotating standard star ( 21peg , @xmath2314 kms@xmath18 , e.g. , @xcite ) acquired with the same setup and in the same observing night as the target system . the limb darkening coefficient was fixed at the theoretically predicted values , 0.42 for the system @xcite . we calibrated the relationship between the ccf gaussian width and @xmath19 using the @xcite data sample . this analysis yielded projected rotational velocities for the components of v446cep as @xmath24=120  kms@xmath18 , and @xmath25  kms@xmath18 . the mean deviations were 3 and 9 kms@xmath18 , for the primary and secondary , respectively , between the measured velocities for different lines . we also performed a spectral classification for the components of the system using compo2 , an idl code for the analysis of high - resolution spectra of sb2 systems written by one of us ( see , e.g. , * ? ? ? * ) and adapted to the tfosc spectra . this code searches for the best combination of two reference spectra able to reproduce the observed spectrum of the system . we give , as input parameters , the radial velocities and projected rotational velocities @xmath22 of the two components of the system , which were already derived . the code then finds , for the selected spectral region , the spectral types and fractional flux contributions that better reproduce the observed spectrum , i.e. which minimize the residuals in the collection of difference ( observed@xmath26composite ) spectra . for this task we used reference spectra taken from the @xcite @xmath27 ( with a a resolution of @xmath281 ) that are representative of stars with spectral types from late - o type to early - a , and luminosity classes v , iv , and iii . the atmospheric parameters of these reference stars were recently revised by @xcite . we selected 198 reference spectra spanning the ranges of expected atmospheric parameters , which means that we have searched for the best combination of spectra among 39204 possibilities per each spectrum . the observed spectra of v446cep in the @xmath1765256720 spectral region were best represented by the combination of hd187459 ( b0.5 ii ) and hd178125 ( b8 iii ) . however , we have adopted , for each component , the spectral type and luminosity class with the highest score in the collection of the best combinations of templates , where the score takes into account the goodness of the fit expressed by the minimum of the residuals . we have thus derived a spectral types for the primary and secondary component of v446cep as b1 and b9 main - sequence stars , with an uncertainty of about 1 spectral subclass . the effective temperature and surface gravity of the two components of the system are obtained as the weighted average of the values of the best spectra at phases near to the quadratures combinations of templates adopting a weight @xmath29 , where @xmath30 is the average of residuals for the @xmath31-th combination . the standard error of the weighted mean was adopted for the atmospheric parameters . both stars appear to have a solar metallicity , within the errors . the atmospheric parameters obtained by the code and their standard errors are reported in table3 . the observed spectra of v446cep at phases near to the quadratures are shown in fig.2 together with the combination of two reference spectra which gives the best match . the light curve of v446cep was obtained by the @xmath0 spacecraft and is composed of 124 photometric points . the light curve is very similar to that of a detached algol - type binary . the brightness of the system at the maximum and depth of the primary minimum were estimated by @xcite to be 7.31 and 0.14 mag , respectively . however , we estimate the brightness of the system at the maximum light as hp=7.348 mag , with a mean error of about 0.013 mag , and the depth the primary minimum as 0.10 mag . the hp magnitudes were transformed to the johnson s v - passband using the coefficients given by @xcite . the @xmath0 light curve is plotted against the orbital phase in fig.3 . there is no indications of any asymmetry in the lc . as the secondary minimum occurs at phase 0.5 , we have adopted circular orbits for our analysis . the effective temperatures for the primary and secondary star were estimated from the spectra as 26580@xmath2880 and 12000@xmath21050 k , respectively . the apparent visual magnitude and colour indices were given by @xcite as @xmath32=7@xmath33.32 , ( @xmath34)=@xmath260@xmath33.71 , and ( @xmath35)=0@xmath33.08 . we obtained the reddening - free quantity @xmath36=@xmath260.768 , which corresponds to a b1 main - sequence star @xcite with an intrinsic colour of ( @xmath35)=@xmath260@xmath33.273 . this colour index corresponds to an effective temperature of 26800@xmath2800 k @xcite , which is consistent with that obtained directly from the spectra . a preliminary analysis of the light curve gives a light ratio of @xmath37/@xmath38=0.014 . using the intrinsic colour of the primary star , the light ratio and the observed composite colour of the system an interstellar reddening of @xmath39=0@xmath33.338 was determined for the system . we started to analyze the light curve using the wilson - devinney code ( hereafter wd ; e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) as implemented in the software phoebe @xcite . the wd code is widely used for determination of the orbital parameters of the eclipsing binaries . to run the code we need some initial parameters . the initial logarithmic limb - darkening coefficients were taken from the tables given by @xcite as @xmath40=0.41 and @xmath41=0.62 , @xmath42=0.24 @xmath42=0.30 and are automatically interpolated at each iteration by phoebe . the effective temperature of the primary star is taken as 26600 k and the ratio of the masses of two components @xmath43=0.147 . we have started with @xmath44 meant for detached binary systems , keeping the temperature of the primary and the mass - ratio as fixed parameters . according to the wd code we adjusted the following parameters : @xmath31 ( the orbital inclination ) , @xmath45 ( the potential for the primary ) , @xmath46 ( the potential for the secondary ) , @xmath47 ( the effective temperature of the cool star ) , l@xmath48 ( the luminosity of the primary ) , and the zero - epoch offset . the luminosity of the secondary star , l@xmath49 , was constrained by the model . after a few numbers of runs of the dc program in @xmath44 the sum of residuals squared showed a minimum and the corrections to the adjustable parameters became smaller than their probable errors . the results are given in the table4 . the corresponding computed light curve is shown in fig.3 as a continuous line . based on the results of radial velocities and light curves analyses we have calculated the physical properties of the v446cep . for this purpose , we used the @xmath50 code developed by @xcite . this code is now widely used for derivation of the absolute parameters of the eclipsing binary stars components . it calculates complete error budgets using a perturbation algorithm . the fundamental stellar parameters for the components such as masses , radii , luminosities and their standard deviations have been derived using this code . the astrophysical parameters of the components , and other properties for the stars of v446cep is presented in table5 . the separation between the components was found to be 28.10@xmath20.63 r@xmath3 for v446cep . the masses were measured to precision of about 67% , apart from the mass of the secondary star of v446cep , which has an uncertainty of about 11% . on the other hand the radii of the stars have been derived with a precision of better than 5% . the accuracy of any parameter of an eclipsing binary system depends mainly on the coverage of the both spectroscopic and photometric observations and their precision . in addition , the lc solutions are more accurate for totally eclipsing systems . the light curve of v446cep shows total eclipses , but the precision of the photometric measurements is not sufficiently high . despite these drawbacks , the physical parameters of the components of system could be determined with sufficient precision . we note that the effective temperature of the secondary star derived from the spectra is in good agreement with that obtained from the light curve analysis . the luminosities and absolute bolometric magnitudes are calculated directly from the radii and and effective temperatures of the components . the effective temperature of 5777  k and the absolute bolometric magnitude of 4.74 mag were adopted for the sun ( e.g. , * ? ? ? the bolometric corrections were interpolated from the tables of @xcite . the v - band magnitude of the system at out - of - eclipse phases is taken as 7.32 for v446cep . we have calculated the absolute visual magnitudes for the components using the fractional luminosities and bolometric corrections given in table4 and 5 . combining these values with the interstellar absorption of 1.09mag for v446cep we have estimated the distance to the system as 1100@xmath262pc .    in the log @xmath20@xmath51 ( left panel ) and log @xmath20-logg planes ( right panel ) of fig.4 we have plotted the positions of the components , with 1-@xmath52 error bars . the filled and empty circles represent the primary and secondary star of v446cep . the evolutionary tracks and isochrones for the non - rotating single stars with solar composition are taken from @xcite . v446cep was included in the list of probable members of the cepob2 association ( e.g. , * ? ? ? * and reference therein ) . according to @xcite the distance to the cepheus ob2 association is @xmath14900pc and it is embedded in the cepheus bubble , a giant shell structure of atomic and molecular gas extending in a radius of about 120pc , which is believed to be generated by an earlier generation of hot and massive stars in ngc7160 . @xcite derived a reddening in the range 0.291.12 for the stars in the cepob2 field , with a mean value of 0.59 for the members of the association . the interstellar reddening that we estimated for v446cep , @xmath39=0.34 , is lower than the value estimated for the association , but it is consistent with the value of 0.32 derived by @xcite for the binary system lzcep in this association . this result suggests that the binary is located near the ngc7160 in the border of cepheus bubble @xcite . the open clusters tr37 and ngc7160 with the ages of 4 and 12 myr locate in the association @xcite . our distance of v446cep agrees quite well with those distance values of cepob2 association . both components of v446cep appear to have an age of about 10 myr which is consistent with the age of stars in the open cluster ngc7160 . as seen from fig.4 , components of the system are located on the main - sequence , i.e. in the central hydrogen burning phase . they are still inside their lobes , corresponding to the detached algols . @xcite collected the physical parameters of 74 detached algols , mainly composed of hot ( baf - type ) stars . their mass - ratios , m@xmath53/m@xmath54 , are generally larger than 0.5 , with a mean value of 0.88 . the binary system arcas has the smallest mass - ratio in their table1 , with a value of 0.315 . we find v446cep as a detached algol with a lower mass - ratio . @xmath55 , which is among lowest values found for this class of binaries . the less massive stars in the semi - detached binaries fill their corresponding roche lobes and are oversized and over - luminous relative to a zero - age main - sequence star of the same mass @xcite . although the less - massive secondary star of v446cep does not fill its respective roche volume it is seen as oversized and over - luminous with respect to its mass , similar to the secondaries of the semi - detached systems . v446cep is close eclipsing binary containing high - mass star . we carried out spectroscopic observations of the system . the atmospheric parameters of the stars in the eclipsing pair have been determined from its spectra . the spectra were analyzed using cross - correlation for measuring the radial velocities of both components and with an ad - hoc code for deriving their atmospheric parameters . moreover , @xmath56 light curve was modeled using the wd code . the physical parameters for the system are measured to accuracies of 6 - 7% in mass , apart from the secondary of v446cep , and 5% in radius . the distance to the system v446cep was estimated as 1100@xmath262pc . a comparison of physical parameters of the components with the theoretical models of single stellar evolution models has been made and an age of about 10myr for v446cep has been derived . both the estimated age and distance of v446cep confirm its membership of cepheus ob2 association , locating close to the open cluster ngc7160 . we thank to tbitak national observatory ( tug ) for a partial support in using rtt150 telescope with project number 11brtt150 - 198 . we thank to ebiltem ege university research center for a partial support with project number 2013/bil/018 . we also thank to the staff of the bakirlitepe observing station for their warm hospitality . this study is supported by turkish scientific and technology council under project number 112t263 . this research was also partly supported by the scientific research projects coordination unit of istanbul university . project number 3685 . we thank canakkale onsekiz mart university astrophysics research center and ulupinar observatory together with istanbul university observatory research and application center for their support and allowing use of ist60 telescope . this work was partially supported by the italian _ ministero dellistruzione , universit e ricerca _ the following internet - based resources were used in research for this paper : the nasa astrophysics data system ; the simbad database operated at cds , strasbourg , france ; and the ar@xmath57iv scientific paper preprint service operated by cornell university . conti , p.  s. , & ebbets , d. , 1977 , apj , 213 , 438 drilling j. s. , landolt a. u. , 2000 , allen s astrophysical quantities , 4th ed . edited by arthur n. cox . isbn : 0 - 387 - 98746 - 0 . publisher : new york : aip press ; springer , 2000 , p.381 ekstrm s. , georgy c. , eggenberger p. , et al . , 2012 , a&a , 537 , a146 flower , p.  j. , 1996 , apj , 469 , 355 frasca , a. , guillout , p. , marilli , e. 2006 , a&a , 454 , 301 garmany , c. d. , & stencel , r. e. , 1992 , a&as , 94 , 211 harmanec , p.  1998 , a&a , 335 , 173 harries , t.  j.,hilditch , r.  w. , & hill , g. , 1998 , mnras , 295 , 386 hilditch , r.  w. , 2004 , aspc , 318 , 198 hovhannessian r.  k. , 2004 , ap , 47 , 499 ibanoglu c. , cakirli , . , & sipahi , e. , 2013a , mnras , 436 , 750 ibanoglu , c. , akirli , . , & sipahi , e. , 2013b , newa , 25 , 68 ibanoglu , c. , soydugan , f. , soydugan , e. , & derviolu , a. , 2006 , mnras , 373 , 435    kane , s.  r. , schneider , d.  p. , & ge , j. , 2007 , mnras , 377 , 1610 kazarovets e.  v. , samus n.  n. , durlevich o.  v. , et al . , 1999 , ibvs , 4659 , 1 kun , m. , kiss , z.  t. , & balog , z.  2008 , handbook of star forming regions , volume i , 136 kroupa , p. 2001 , mnras , 322 , 231 malkov , o.  y. , oblak , e. , snegireva , e.  a. , & torra , j. , 2006 , a&a , 446 , 785 patel , n.  a. , goldsmith , p.  f. , heyer , m.  h. , snell , r.  l. , & pratap , p. , 1998 , apj , 507 , 241 penny , l.  r. , 1996 , apj , 463 , 737 perryman m. a. c. , lindegren l. , kovalevsky j. et al . , 1997 , a&a , 323 , l49 pra a. , zwitter t. , 2005 , apj , 628 , 426p royer f. , gerbaldi m. , faraggiana r. , & gomez a. e. 2002 , a&a , 381 , 105 reed , b.  c. , 2003 , aj , 125 , 2531 sana , h. , dunstall , p.  r.,hnault - brunet , v. , et al . 2012 , proceedings of a scientific meeting in honor of anthony f.  j.  moffat , aspc , 465 , 284 salpeter , e.  e. 1955 , apj , 121 , 161 sicilia - aguilar , a. , hartmann , l.  w. , briceo , c. , muzerolle , j. , & calvet , n. , 2004 , aj , 128 , 805 simkin , s.  m. , 1974 , a&a , 31 , 129 simonson , s.  c. , iii , 1968 , apj , 154 , 923 southworth , j. , maxted , p.  f.  l. , & smalley , b. , 2005 , a&a , 429 , 645 tonry , j. , & davis , m.  1979 , aj , 84 , 1511 torres g. , andersen j. , & gimnez a. , 2010 , a&arv , 18 , 67 valdes , f. , gupta , r. , rose , j.  a. , singh , h.  p. , & bell , d.  j. , 2004 , apjs , 152 , 251 van hamme , w. , 1993 , aj , 106 , 2096 wilson , r.  e. , 1979 , apj , 234 , 1054 wilson , r.  e. , 2006 , astrophysics of variable stars , aspc , 349 , 71 wilson , r.  e. , & devinney , e.  j. , 1971 , apj , 166 , 605 wu , y. , singh , h. p. , prugniel , p. , gupta , r. , & koleva , m. 2011 , a&a , 525 , 71                  .heliocentric radial velocities of v446cep . the columns give the heliocentric julian date , orbital phase and the radial velocities of the two components with the corresponding standard deviations . [ cols="^,^,>,^,>,^,^ " , ]
Abstract:

Answer: we present new spectroscopic observations of the early type , double - lined eclipsing binary v446cep . the radial velocities and the photometric data obtained by @xmath0 were analysed for deriving the astrophysical parameters of the components . masses and radii were determined as m@xmath1=17.94@xmath21.16 m@xmath3 and r@xmath1=8.33@xmath20.29 r@xmath3 , m@xmath4=2.64@xmath20.30 m@xmath3 and r@xmath4=2.13@xmath20.10 r@xmath3 for the components of v446cep . our analyses show that v446cep is a detached algol - type system . based on the position of the components plotted on the theoretical hertzsprung - russell diagram , we estimate that the age of v446cep is about 10 myr , neglecting the effects of mass - loss and mass exchange between the components . using the ubvjhk magnitudes and interstellar absorption we estimated the mean distance to the system v446cep as 1100@xmath262pc . stars : binaries : eclipsing  stars : fundamental parameters  stars : binaries : spectroscopic  stars : v446cep