Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: in the standard model of low - mass star formation , molecular clouds have an initially subcritical mass - to - flux ratio , and spend a relatively long time ( @xmath2 years ) undergoing quasi - static fragmentation by ambipolar diffusion until cores come to be supercritical . after a core becomes supercritical , it collapses on a dynamical timescale ( @xmath3 years ) to form either one or a small multiple system of stars ( e.g. , * ? ? ? * ; * ? ? ? models of nonaxisymmetric fragmentation by this mechanism have been presented recently in several papers ( * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * hereafter kboy ) . this scenario is supported by the observed relatively low efficiency of star formation in molecular clouds , i.e. , only @xmath4 for clouds as a whole but several times greater within cluster - forming cores ( see * ? ? ? an alternative model is that star formation begins in clouds with supercritical mass - to - flux ratio that have additional support due to turbulence ( e.g. , * ? ? ? * ; * ? ? ? * ) . in this model , star formation occurs relatively rapidly , as turbulence dissipates over a dynamical timescale , and the large - scale magnetic field plays a minor role . the rapid star formation model is supported by the observational results that the age spreads of young stars in nearby molecular clouds are often a few @xmath5 years @xcite , and the low fraction of clouds that are observed to be in a pre - star - formation state . recently , @xcite and @xcite have modeled mildly subcritical clouds and shown that the timescale of cloud fragmentation is reduced by supersonic turbulence . they showed this by performing two - dimensional simulations in the thin - disk approximation . such a model can explain both relatively rapid star formation and the relatively low star formation efficiency in molecular clouds that is not well explained if star formation starts from a supercritical cloud . in this letter , we study the three - dimensional extension of the @xcite model , by including a self - consistent calculation of the vertical structure and dynamics of the cloud . we confirm that supersonic flows can significantly shorten the timescale of core formation in three - dimensional subcritical clouds , and clarify how the scaling of various physical quantities allows this to occur . we solve the three - dimensional magnetohydrodynamic ( mhd ) equations including self - gravity and ambipolar diffusion , assuming that neutrals are much more numerous than ions . instead of solving a detailed energy equation , we assume isothermality for each lagrangian fluid particle @xcite . for the neutral - ion collision time and associated quantities , we follow @xcite . the basic equations are summarized in kboy .    as an initial condition , we assume hydrostatic equilibrium of a self - gravitating one - dimensional cloud along the @xmath6-direction in a cartesian coordinate system @xmath7 . though nearly isothermal , a molecular cloud is usually surrounded by warm material , such as neutral hydrogen gas . hence , we assume that the initial sound speed makes a transition from a low value @xmath8 to a high value @xmath9 at a distance @xmath10 , with a transition length @xmath11 ( see eq . 16 of kboy ) . we take @xmath12 , @xmath13 , and @xmath14 , where @xmath8 is the initial sound speed at @xmath15 , @xmath16 , and @xmath17 is the initial density at @xmath15 . a numerical solution for the initial density distribution shows that it is almost the same as the @xcite solution for an equilibrium isothermal layer in the region @xmath18 . we also assume that the initial magnetic field is uniform along the @xmath6-direction . a set of fundamental units for this problem are @xmath8 , @xmath19 , and @xmath17 . these yield a time unit @xmath20 . the initial magnetic field strength ( @xmath21 ) introduces a dimensionless free parameter @xmath22 the ratio of gas to magnetic pressure at @xmath15 . in the above relation , we have also used @xmath23 , the initial alfvn speed at @xmath15 . in the sheet - like equilibrium cloud with a vertical magnetic field , @xmath24 is related to the mass - to - flux ratio for spitzer s self - gravitating layer . the mass - to - flux ratio normalized to the critical value is @xmath25 , where @xmath26 is the column density of the spitzer cloud . therefore , @xmath27 . although the initial cloud we used is not exactly the same as the spitzer cloud , @xmath24 is a good indicator of whether or not the magnetic field can prevent gravitational instability . for the model presented in this letter , we choose @xmath28 so that the cloud is slightly subcritical . for this choice , @xmath29 . dimensional values of all quantities can be found through a choice of @xmath17 and @xmath8 . for example , for @xmath30 km s@xmath31 and @xmath32 @xmath33 , we get @xmath34 pc , @xmath35 yr , and @xmath36 g if @xmath28 . the level of magnetic coupling in the partially ionized gas is characterized by numerical values of the ion number density @xmath37 and neutral - ion collision timescale @xmath38 . from eqs . ( 8) and ( 9 ) of kboy , and using standard values of parameters in that paper as well as the values of units used above , we find an initial midplane ionization fraction @xmath39 and a corresponding neutral - ion collision time @xmath40 . the ionization fraction @xmath41 and timescale @xmath38 at other densities can be found from the initial midplane values given that they both scale @xmath42 @xcite .    in this equilibrium sheet - like gas layer , we input a nonlinear perturbation to @xmath43 and @xmath44 at each grid point . independent realizations of @xmath43 and @xmath44 are generated in fourier space with amplitudes drawn from a gaussian distribution and consistent with power spectrum @xmath45 . here , @xmath46 is the total wavenumber . our adopted power spectrum means that the large scale has more energy than the small scale . this perturbation is consistent with that implemented in @xcite . the rms value of the initial velocity perturbation in physical space , @xmath47 , is about @xmath48 , so that @xmath49 as well . the method of solution and boundary conditions are described by kboy ( see also * ? ? ? * ; * ? ? ? the computational region is @xmath50 and @xmath51 , with a number of grid points for each direction @xmath52 . the top panel of figure 1 shows the time evolution of the maximum density @xmath53 at @xmath15 . the simulation is stopped when @xmath54 . the solid line shows the result when we input an initially nonlinear supersonic perturbation to the @xmath28 model . the dashed line ( @xmath28 ) and the dash - dotted line ( @xmath55 ) show the result in the case of the small ( linear ) initial perturbation that we have applied in kboy . this figure shows that the timescale of collapsing core formation for the nonlinear perturbation case is much shorter than that for the linear perturbation case , when @xmath24 is the same . even when the initial cloud is subcritical ( @xmath28 ) , the core formation occurs on almost the same timescale as that of the supercritical ( @xmath55 ) linear perturbation case . , from the initial flow . this may be _ too _ rapid to agree with observations . ] the bottom panel of figure 1 shows the time evolution of the maximum value of density at @xmath15 ( @xmath53 ) and @xmath56 at the location of maximum density ( @xmath57 ) . at first , @xmath57 increases rapidly up to @xmath58 due to rapid ambipolar diffusion in the highly compressed regions caused by the initial supersonic perturbation . however , there is enough stored magnetic energy in the compressed region that it rebounds and starts oscillations , with @xmath57 around @xmath59 and increasing gradually . eventually , @xmath57 becomes @xmath60 and the dense region collapses to form a core . this figure implies that @xmath56 is a good indicator to see whether a subregion of the cloud is supercritical or not . the evolution of @xmath53 confirms that there is an initial compression followed by a rebound to a lower density ( still greater than the initial background value ) and subsequent oscillations until a runaway collapse starts when continuing ambipolar diffusion has created a region with @xmath61 . figure 2 shows an image of the logarithmic density at the last snapshot ( @xmath62 , when @xmath63 ) for the nonlinear perturbation case of the subcritical cloud ( @xmath28 ) . the top panel shows the cross section at @xmath15 , and the bottom panel shows the cross section at @xmath64 . the value of @xmath65 for the bottom panel is chosen so that the vertical cut passes through the maximum density point . a collapsing core is located in the vicinity of @xmath66 . the size of the core is similar to that created by linear initial perturbations ( see fig . 2 in kboy ) , although the shape is notably less circular . the top panel of figure 3 shows the density and @xmath67-velocity of neutrals and ions along an @xmath67-axis cut at @xmath68 , taken from the snapshot illustrated in figure 2 . the @xmath67-velocities show infall motion toward the center of the core , although the core itself is moving with nonzero negative @xmath67-velocity . the relative infall speed to the core is subsonic and about @xmath69 . it is comparable to the case of initial linear perturbation ( see fig . 9 in kboy ) . however , there are systematic motions throughout the simulation region that are still supersonic at this time , even though the initial turbulent energy has decayed somewhat . this is qualitatively distinct from the corresponding case with linear perturbations , and should be observationally testable . for the single core that is formed in this simulation , the systematic core motion is about @xmath70 in the @xmath67-direction . however , it can be even larger in other realizations . the ion velocity @xmath71 can be related to the neutral velocity @xmath72 and the magnetic acceleration @xmath73 in the assumed limit of low ion inertia by the relation @xmath74 figure 3 reveals that the ion infall motions are smaller in magnitude than the neutral motions due to the retarding magnetic acceleration ; however the relative drift in the @xmath67-direction between ions and neutrals is typically within @xmath75 . in order to see the force balance in the core , we plot the @xmath67-accelerations in the bottom panel of figure 3 along the same @xmath67-axis as the top panel . these values are normalized by the computational units . the inward gravitational acceleration ( @xmath76 ) dominates around the core . the acceleration from the magnetic force ( @xmath77 ) and thermal pressure force ( @xmath78 ) resist the contraction . the net acceleration ( @xmath79 ) is working toward the center of the core . the top panel of figure 4 shows the density and @xmath6-velocities along a @xmath6-axis cut at @xmath80 from the snapshot illustrated in figure 2 . the @xmath6-velocities also show infall motion toward the center of the core . the relative neutral infall speed to the core is also subsonic and about @xmath81 . the ion velocity in the vertical direction is almost identical to the neutral velocity but actually very slightly _ greater _ in magnitude at the peak . the bottom panel of figure 4 shows the @xmath6-accelerations along the same @xmath6-axis as in the top panel . in the @xmath6-direction , @xmath76 is nearly balanced by @xmath78 , and @xmath77 is almost negligible , _ but working in the same direction as gravity_. the net acceleration ( @xmath79 ) is also working toward the center of the core . the inward pointing vertical magnetic acceleration associated with an hourglass - type magnetic field morphology explains why the ion velocity in the @xmath6-direction is greater in magnitude than the neutral velocity . our three - dimensional mhd simulations have shown that the supersonic nonlinear flows significantly reduce the timescale of collapsing core formation in subcritical clouds . it is of order several @xmath5 years for typical parameters , or @xmath82 times less than found in the linear initial perturbation studies of @xcite , @xcite , and kboy . our supersonic perturbations are also trans - alfvnic , in agreement with analysis of observed magnetic field strengths @xcite . it is also in line with the theoretical result of @xcite that turbulent clouds will settle into a trans - alfvnic state when global motion / expansion is allowed , even if driven with turbulence that is initially super - alfvnic .    to see how accelerated ambipolar diffusion can occur , we note that the magnetic induction equation ( see , e.g. , eq . 3 of kboy ) and our assumptions of ionization balance can be used to estimate the diffusion time @xmath83 , where @xmath84 is the gradient length scale introduced by the initial turbulent compression and @xmath85 is the magnetic field strength . because the compression by the nonlinear flow is nearly one - dimensional , the magnetic field scales roughly as @xmath86 within the flux freezing approximation . if the compression is rapid enough that vertical hydrostatic equilibrium can not be established ( unlike in previous calculations using the thin - disk approximation ) , then @xmath87 as well ( i.e. , one - dimensional contraction without vertical settling ) , and @xmath88 . this means that diffusion can occur quickly ( and lead to a rapidly rising value of @xmath56 ) if the turbulent compression creates small values of @xmath84 . if diffusion is so effective during the first turbulent compression that a dense region becomes magnetically supercritical , then it will evolve directly into collapse . alternately , the stored magnetic energy of the compressed ( and still subcritical ) region may lead to a reexpansion of the dense region . the timescale for this , in the flux - freezing limit , is the alfvn time @xmath89 , which scales @xmath90 for the above conditions . thus , @xmath91 decreases more rapidly than @xmath92 , and sufficiently small turbulent - generated values of @xmath84 may lead to enough magnetic diffusion that collapse occurs before any reexpansion can occur . see @xcite for some similar discussion along these lines . ultimately , whether or not reexpansion of the first compression can occur depends on the strength of the turbulent compression , mass - to - flux ratio of the initial cloud , and neutral - ion collision time . in a preliminary study of such effects , we found that keeping @xmath47 and @xmath24 fixed but allowing significantly poorer or stronger neutral - ion coupling yielded differing results . for @xmath93 ( poorer coupling ) , a collapsing core formed due to the first compression , at @xmath94 , while for @xmath95 it happened after more oscillations than in our standard model , at @xmath96 . a full parameter study will be required to elucidate . if reexpansion of the initial compression does occur , as in the standard model presented in this paper , then there is enough time for the vertical structure to settle back to near - hydrostatic equilibrium , in which case @xmath97 . since the compressed and reexpanded region executes oscillations about a new mean density , it is convenient to analyze the scalings in terms of the density @xmath98 . the diffusion time now scales as @xmath99 . this yields a scaling of @xmath91 that is the traditionally used one ( and is satisfied by design in the thin - disk approximation ) . however , the diffusion occurs more rapidly than it would in the initial state due to the elevated value of @xmath98 in the compressed but oscillating region ( see figure 1 ) .
Abstract:

Answer: we demonstrate that the formation of collapsing cores in subcritical clouds is accelerated by nonlinear flows , by performing three - dimensional non - ideal mhd simulations . an initial random supersonic ( and trans - alfvnic ) turbulent - like flow is input into a self - gravitating gas layer that is threaded by a uniform magnetic field ( perpendicular to the layer ) such that the initial mass - to - flux ratio is subcritical . magnetic ambipolar diffusion occurs very rapidly initially due to the sharp gradients introduced by the turbulent flow . it subsequently occurs more slowly in the traditional near - quasistatic manner , but in regions of greater mean density than present in the initial state . the overall timescale for runaway growth of the first core(s ) is several @xmath0 yr , even though previous studies have found a timescale of several @xmath1 yr when starting with linear perturbations and similar physical parameters . large - scale supersonic flows exist in the cloud and provide an observationally testable distinguishing characteristic from core formation due to linear initial perturbations . however , the nonlinear flows have decayed sufficiently that the relative infall motions onto the first core are subsonic , as in the case of starting from linear initial perturbations . the ion infall motions are very similar to those of neutrals ; however , they lag the neutral infall in directions perpendicular to the mean magnetic field direction and lead the neutral infall in the direction parallel to the mean magnetic field .


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: nonuniform growth is inherently present in a broad class of phenomena including the development of biological populations , communication networks or economic systems like incomes of persons or companies @xcite . in many cases it is obvious to assume that in a system consisting of groups or clusters of units the attachment of a new entity to one of the groups depends on the already achieved strength or size of that particular group . simon  @xcite analysed a simple model of this kind where the growth probability was proportional to the cluster size and he gave exact results for the time dependent size distribution . referring to the examples of words in a book or personal incomes simon derived a power law distribution of cluster sizes . recently , in the search for an explanation of the widely observed scale invariance of large networks like the www @xcite , the internet or power networks @xcite , scientific citation @xcite the idea of preferential growth has been applied to evolving graphs @xcite . it turned out that such graphs behave remarkably : they have `` small world '' properties @xcite and the distribution of the strength of vertices ( number of edges from or to a vertex ) is scale free , provided that the probability of linking a vertex with a new one is proportional to its strength @xcite this class of models represent a new mechanism for `` self - organized criticality''@xcite . the idea of preferential growth seems to be essential in economic systems too where clustering of companies , e.g. , according to their market seem to follow such a pattern @xcite . these models have been treated by different tools including simulations , continuum or mean field theories @xcite and exact calculations @xcite by which information has been accumulated about the asymptotic behavior and the time dependence of the global distribution functions . however , much less attention has been paid to the full time - dependent solution of the problem . the aim of the present work is to give such a solution of a particular model . the paper is organized as follows . in section ii we define the model and the quantities of interest as well as we present the basic master equation . in section iii the main steps of the full time dependent analytic solution is given and the consequences for the steady state and the integrated distributions are drawn . section contains the analysis about the asymptotics and scaling . in section v we present a discussion of our results . the paper terminates with two appendices containing some details of the calculations . we model a growing system which consists of groups of different sizes . at the beginning ( @xmath5 ) we have one group with one element in it . at each time step we add a new element to the system . with probability @xmath6 it will belong to one of the existing groups . the probability that it joins the @xmath2-th group is proportional to the size of the group ( @xmath7 ) , see fig  [ fig_explain ] . ( the number of elements is equal to the time , @xmath8 , because the system size is rising by one in each time step . ) with probability @xmath9 the new element will belong to a new group .    the process can be described by the following master equation : [ master ] _ i(k , t ) & = & p(k-1 ) t-1 _ i(k-1,t-1 ) + p(1-k t-1 ) _ i(k , t-1 ) + + & & +  ( 1-p ) _ i(k , t-1 ) + ( 1-p ) _ i-1(t-1 ) _ k,1(1-_i,1 ) , where @xmath1 is the probability that at time @xmath3 there are @xmath4 elements in the group @xmath2 , and @xmath10 is the probability that at time @xmath3 there are @xmath2 groups in the system : [ boxprob ] _ i(t ) = t-1 i-1  p^t-1-(i-1 )  ( 1-p)^i-1 . in the following we introduce some important quantities and their definitions .    given the size distribution of the individual groups , @xmath1 , the size distribution of the total system can be calculated as their average : [ def_p_k , t ] ( k , t ) = 1 t _ i=1^t  _ i(k , t ) . in the long time limit this quantity approximates to a stationary value : @xmath11 . the mean of the @xmath2-th group size : [ def_mean ] ( t ) = _ k=1^t - i+1  k  _ i(k , t ) . the reason that the upper limit of the above sum is not infinity is that @xmath12 if @xmath13 . in the first step we calculate the group size distribution in the asymptotic case , @xmath14 . the exact analytic formula for @xmath14 was already calculated in @xcite , we present it here to see the dependence of the exponent on the parameter @xmath6 .    if we sum up eq . ( [ master ] ) for @xmath15 , we get : [ recurs_p ] t  ( k , t ) = ( t-1-pk ) ( k , t-1 ) + p(k-1 ) ( k-1,t-1 ) + ( 1-p)_k,1  , since : _ i-1(t-1 ) ( 1-_i,1 ) & = & 1 , + _ i(k , t-1 ) = _ i=1^t-1 _ i(k , t-1 ) & = & ( t-1 ) ( k , t-1 ) . the stationary behavior of @xmath16 , mentioned in the previous section , can be checked from eq . ( [ recurs_p ] ) . replacing the stationary quantity @xmath14 into eq . ( [ recurs_p ] ) one gets : ( k ) = -pk(k ) + p(k-1)(k-1 ) + ( 1-p ) _ k,1 , which can be solved for @xmath14 : [ p_k , inf ] ( k ) = ( k)(2 + 1 p ) ( k+1 + 1p ) 1-p 1+p_k k^-1 - 1/p . [ fig_p_k , inf ]      in the model the first group has an accentuated role since it always has at least one element because of the initial conditions . therefore the master equation ( [ master ] ) for the first group ( @xmath17 ) has the following simpler form : [ i1_master ] _ 1(k , t ) = _ 1(k , t-1 )  -  p t-1  k  _ 1(k , t-1 )  +  ( k-1 )  _ 1(k-1,t-1 ) . for @xmath18 in the above equation on the r.h.s the last term vanishes so the probability @xmath19 can be calculated easily : _ 1(1,t ) = ( t - p ) ( t ) ( 1-p ) . for @xmath20 one can prove ( see appendix [ sec_sum ] ) that the following equality holds : [ k>1_i1_biz ] _ k=1^l  ( -1)^k-1  l-1 k-1  _ 1(k , t ) = ( t - lp ) ( t ) ( 1-lp ) .    the analytic form of @xmath21 can be received from eq . ( [ k>1_i1_biz ] ) by multiplying both sides with @xmath22 and summing up for @xmath23 , [ pm_1,k ] _ 1(k , l=1^k  ( -1)^l-1  k-1 l-1 . in the case of @xmath24 we have to look at the hole master equation ( [ master ] ) . in this case the equality ( [ k>1_i1_biz ] ) does nt hold because of the last factor in ( [ master ] ) . our assumption is that the probability @xmath1 will have a modified form : [ pm_i , k , t ] _ i(k , t )  = ^_1(k , t ) . the validity of the above form can be checked by replacing it back in eq . ( [ master ] ) , see appendix [ sec_i>1 ] .      replacing the analytic formula ( [ pm_i , k , t ] ) into ( [ def_mean ] ) one gets : ( t ) = _ k=1^t - i+1 k _ l=1^k  ( -1)^l-1 k-1 l-1 ( t - lp ) ( t)(1-lp ) , the two sums can be transposed @xmath25 , and @xmath26 so the mean value will have the following form : [ mean - k , t ] ( t ) = _ l=1^t - i+1 ( -1)^l-1 l t - i+2 l+1 ( t - lp ) ( t)(1-lp ) _ b = i^t ( b)(1-lp ) ( b - lp ) b-2 i-2 p^b - i ( 1-p)^i-1 .      in sec . [ sec_p_k , inf ] we calculated the stationary group size distribution directly from the master equation . now we are interested in its dynamic . in order to compute that , we start from the definition ( [ def_p_k , t ] ) of @xmath16 , and replace the solution we got for @xmath1 in the previous sections ( [ pm_i , k , t ] ) . [ p_k , t _ _ ] ( k , t ) & = & 1 t _ l=1^k ( -1)^l-1 k-1 l-1 ( t - lp ) ( t)(1-lp ) transposing the two sums : @xmath27 , and taking into account that : @xmath28 one finally arrives at the time dependent distribution : [ p_k , t ] ( k , t )  =   _ ( k , ) +  p+lp 1+lp ] in the long time limit we will get back our result ( [ p_k , inf ] ) since the second term in @xmath29 $ ] decays for large @xmath3 values with @xmath30 , and the sum transforms into : ( k , )  =  1-p 1+p   ( k)(2 + 1/p ) ( k+1 + 1/p ) . we study the @xmath31 limes of @xmath1 and @xmath32 . in the analytic formula for @xmath1 , see . ( [ pm_i , k , t ] ) , there are two components , @xmath33 and @xmath34 , that depend on the time . [ beta_pm_i , k , t ] _ i(k , t ) = ( 1-p)^i-1 _ l=1^k ( -1)^l-1 k-1 l-1 _ a(l , p , t ) _ b(i , l , p , t ) the limes of the first term , @xmath33 , can be easily calculated _ t a(l , p , t ) = 1 ( 1-lp ) t^-lp the second term in the long time limit @xmath35 will converge to a hypergeometric sum : _ t b(i , l , p , t ) = ( i , l , p ) = ( i ) ( 1-lp ) ( i - lp )  _ 2f_1(i , i-1;i - lp;p ) .    for large time values the only time dependent term in ( [ beta_pm_i , k , t ] ) will be @xmath36 which in case of large @xmath3 is a fast decaying function of @xmath37 . so in the case of @xmath38 we can assume that only the first term of the sum gives non - negligible component for @xmath1 . [ pm_i , k , inf ] _ t _ i(k , t ) = t^-p  ( 1-p)^i-1 ( i ) ( i - p )  _ 2f_1 ( i , i-1;i - p;p ) + o(t^-2p ) for large @xmath2 values the above formula simplifies further , because in that case @xmath39 , and @xmath40 : [ pm_inf , k , inf ] _ i(k , t ) = ( i t ) ^p    to study the asymptotic behavior of @xmath41 we start from the fact , that for small @xmath4 values , @xmath42 , the individual group size distribution , @xmath1 , can be described by the first term of the sum , see eq . ( [ pm_i , k , inf ] ) , and for larger values @xmath43 it has a fast decay , fig . [ fig_pm - k , inf ] . a cut - off parameter , @xmath44 , can be defined and we can assume that ( [ def_mean ] ) transforms into [ mean_inf ] ( t ) _ k=1^k^ * k  _ i(k , t ) = _ i(1,t ) k^ * ( k^ * + 1 ) 2 . the definition of @xmath44 can be done in many ways . we defined @xmath44 as the inflection point of @xmath1 , hence : [ def_k_star ] k^ * = t^p  ( i,3,p ) ( i,4,p )  ( 1 - 4p ) ( 1 - 3p ) + 2 + o(t^-p)= t^p  ( i-4p ) ( i-3p ) + 2 + o(t^-p ) . replacing @xmath44 into ( [ mean_inf ] ) ( t ) t^p ( 1-p)^i-1 ( i ) ( i - p )  _ 2f_1 ( i , i-1;i - p;p ) ^2  . for large @xmath2 values the above formula gets a simpler form , because in this case @xmath40 , @xmath45 , @xmath46 , @xmath39 . [ mki_inf ] ( t )  _ i  ( t i ) ^p in this paper we presented a simple preferential growth model consisting of a system of clusters with different sizes . we gave exact solutions for the main characteristic quantities as the distribution , @xmath1 , and the mean value , @xmath32 , of the individual group size as well as for the distribution of the average group size , @xmath16 . the question rises why are such time dependent quantities of interest since most of the asymptotic scaling behavior can be obtained with much less labor . in fact the growth models and network usually provide only a background for some dynamic process  an aspect which has not yet paid enough attention to . if there is a strong separation of time scales , i.e. , the growth is much smaller than the process itself then it is satisfactory to concentrate on the asymptotics only . this is probably the case with the internet or the www . however , in some cases such a separation of scales could be approximate only or even missing and then the importance of the full time dependence becomes apparent . we expect that in certain economic processes this will be the case . an important aspect in the asymptotic scaling is universality . similarly to other preferential growth models , our system exhibits nonuniversal parameter dependent scaling : the exponents depend on the parameter @xmath0 ( the probability of creating a new group ) . it is worth mentioning that the examples quoted in the introduction also show a wide variety of scaling exponents . further interesting study would be to analyse a model where this parameter @xmath0 depends on the time of the growth . the presented system is not a network , the different groups are not linked to each other . however , for a specific value of the parameter , @xmath47 , it can be interpreted as a kind of mean field network model . the clusters then denote the different nodes , and the particles are the links . the value @xmath47 means that in average in every second time step one new group and two elements are created ( in the odd time steps the new element joins to an old group and in even time steps it will create a new group . ) the new group is the new node while the two new elements are the two ends of the new link , one is pointing to the old node , the other is to the new one . this case corresponds to the barabasi s network model with parameter @xmath48 which means that the new node connects to _ one _ old sites . for this particular parameter choice our results agree with them got for the barabasi s network model : @xmath49 , see eq . ( [ p_k , inf ] ) , and @xmath50 , see eq . ( [ mki_inf ] ) . 0.5truecm      99 r. albert , h. jeong , a. l. barabsi , nature * 401 * , 130 ( 1999 ) a. l. barabsi , r. albert science * 286 * , 509 ( 1999 ) b. a. huberman , l. a. adamic , nature * 401 * , 131 ( 1999 ) h. a. simon , biometrica * 42 * , 425 ( 1955 ) , s. n. dorogovtsev , j. f. f. mendes , a. n. samukhin , phys . rev . * 85 * , 4633 ( 2000 ) m. faloutsos , p. faloutsos , c. faloutsos , comput . comm . rev . * 29 * 251 ( 1999 ) s. redner , eur . j. b. * 4 * , 131 ( 1998 ) a. l. barabsi , nature - progress report , _ in press _ m. kochen ( ed . ) , _ the small world _ ( ablex , norwood , nj , 1989 ) , s. milgram , psychology today * 1 * , 61 ( 1967 ) , d. j. watts , s. h. strogatz , nature * 393 * , 440 ( 1998 ) p. l. krapivsky , s. redner , f. leyvraz , phys . * 85 * , 4629 ( 2000 ) p. bak : how nature works ( copernicus , new york , 1996 ) h.j.christensen : self - organized criticality ( cambridge university press , 1998 ) since the absence of characteristic sizes occurs in preferential growth at a natural choice of the parameters it can be considered as an example of soc . of course , the mechanism is totally different from that of the `` sand pile models '' . a similarity with the usual soc models is that change in some of the model parameters leads out of criticality , e.g. , the preferential growth with a nonlinear probability introduces a characteristic size @xcite . n. kiyotaki , r. wright , american economic review * 83 * , 63 ( 1993 ) a. l. barab ' asi , r. albert , and h. jeong , physica a * 272 * , 173 ( 1999 ) s. n. dorogovtsev , j. f. f. mendes , a. n. samukhin , cond - mat/0011115 , s. n. dorogovtsev , j. f. f. mendes , cond - mat/0012009 a. erdlyi , higher transcendental functions , mcgrow - hill book co. , inc . ( 1953 )        if one multiplies ( [ i1_master ] ) by @xmath51 and sums it up for @xmath52 one gets : [ k>1_i1_rekurz ] & & _ k=1^l  ( -1)^k-1  l-1 k-1  _ 1(k , t ) = _ k=1^l  ( -1)^k-1  l-1 k-1  _ 1(k , t-1 )  -   + & & -  p t-1 , where in the first term ( @xmath53 ) we detach the last term of the sum : x = _ k=1^l-1 ( -1)^k-1  k l-1 k-1 _ 1(k , t-1 ) + ( -1)^l-1  l _ 1(l , t-1 ) . taking into account that @xmath54 the second term ( @xmath55 ) can be rewritten as : y = _ k=2^l ( -1)^k-1 l-1 k-1 ( k-1 ) _ 1(k-1,t-1 ) = -  _ k=1^l-1 ( -1)^k-1   ( l - k ) l-1 k-1 _ 1(k , t-1 ) . [ x_min_y ] x - y = l  _ k=1^l ( -1)^k-1 l-1 k-1 _ 1(k , t-1 ) . replacing the difference , @xmath56 , back to eq . ( [ k>1_i1_rekurz ] ) one gets the time evolution of the sum : _ k=1^l  ( -1)^k-1 l-1 k-1 _ 1(k , t )  = t-1-lp t-1 _ k=1^l  ( -1)^k-1 l-1 k-1 _ 1(k , t-1 ) , which leads us back to our assumption ( [ k>1_i1_biz ] ) . the l.h.s of the equation after detaching the last term ( @xmath57 ) of the sum : [ 1term ] _ i(k , t ) = _ l=1^k ( -1)^l-1 k-1 l-1 p^t - i + _ l=1^k ( -1)^l-1 k-1 l-1 ( t - lp ) ( t)(1-lp ) _ b = i^t-1 ( b)(1-lp ) ( b - lp ) b-2 i-2 p^b - i  . the first term of the r.h.s : [ 2term ] t-1-kp t-1  _ i(k , t-1 ) & = & ( -1)^k-1 ( t - kp ) ( t ) _ b = i^t-1  ( b ) ( b - kp ) b-2 i-2 p^b - i  + + & + & t-1-kp t-1   _ l=1^k-1  ( -1)^l-1  k-1 l-1 ( t-1-lp ) ( t-1 ) _ b = i^t-1  ( b ) ( b - lp ) b-2 i-2 p^b - i , taking into account that @xmath58 , the second term will be : [ 3term ] ( k-1)p t-1 _ i(k-1,t-1 ) = p t-1 _ l=1^k-1 ( -1)^l-1 ( k - l ) k-1 l-1 ( t-1-lp ) ( t-1 ) _ b = i^t-1 ( b ) ( b - lp ) b-2 i-2 p^b - i . the sum of ( [ 2term ] ) and ( [ 3term ] ) will be : [ 23term ] & & t-1-kp t-1  _ i(k , t-1 ) + ( k-1)p t-1  _ i(k-1,t-1 ) = + & & =  ( -1)^k-1 ( t - kp ) ( t ) _ b = i^t-1  ( b ) ( b - kp ) b-2 i-2 p^b - i  + _ l=1^k-1  ( -1)^l-1  k-1 l-1   ( t - lp ) ( t)_b = i^t-1  ( b ) ( b - lp )  p^b - i = + & & =  _ l=1^k  ( -1)^l-1  k-1 l-1   ( t - lp ) ( t)_b = i^t-1  ( b ) ( b - lp )  p^b - i , which will be equal to the second term of ( [ 1term ] ) . simplifying with this term the remaining equation : p^t - i   _ l=1^k  ( -1)^l-1  k-1 l-1  = t-2 i-2  p^t - i  _ k,1 . which is true , because the sum equals with @xmath59 .
Abstract:

Answer: we consider a preferential growth model where particles are added one by one to the system consisting of clusters of particles . a new particle can either form a new cluster ( with probability @xmath0 ) or join an already existing cluster with a probability proportional to the size thereof . we calculate exactly the probability @xmath1 that the size of the @xmath2-th cluster at time @xmath3 is @xmath4 . we analyze the asymptotics , the scaling properties of the size distribution and of the mean size as well as the relation of our system to recent network models . pacs numbers : 05.10.-a , 05.40.-a , 02.50.cw    = 0truecm


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 will not be long before the current set of interferometric gravitational wave detectors such as ligo , virgo , geo and tama achieve sensitivities at which detections are expected to be possible with a reasonable event rate . the ground - based interferometric detectors are now starting full - blown  science runs " and the plans for more sensitive second and even third generation detectors are in the drawing boards . even lisa , the space - based observatory , seems to be a possibility in a reasonably near future . this work tries to help in setting the stage for realistic calculations of what now seems to be a very likely occurrence in active galactic nuclei : the capture of compact stellar - sized objects by the supermassive black holes lurking in the center of most galaxies . this type of events will generate gravitational waves within lisa s planned sensitivity band . since the captured objects are small compared to the capturing black hole , they can be treated in a first approximation by black hole perturbation theory . but one thorny problem still remains that has stalled progress in this area . we can model the waves coming from many kinds of stationary particle orbits around the black hole , but real orbits will not be stationary . the emission of gravitational waves will back - react on the object making the orbit decay in time . one of the challenges in this area is to find workable prescriptions for these radiation reaction forces that can be used in numerical simulations that can produce accurate waveforms and energy fluxes for observers at large distances from the hole to be used as templates for the interferometric detectors . the prescriptions that have been worked out so far involve calculating quite complex integrals along the past light cone of the particle @xcite . since almost all perturbative treatments of particles around rotating holes are done in the frequency domain , applying these radiation reaction schemes has proved challenging . frequency domain computations become computationally expensive when one is dealing with realistic non - circular orbits . we try here to work out a simple case of equatorial circular particle orbits around a rapidly rotating black hole to show how to treat this problem perturbatively in the time domain . we hope this might help to open the door to the use of the proposed fully relativistic radiation reaction schemes in the calculation of the more complex but interesting problem of modeling the orbital decay and gravitational radiation emission of particles in more general elliptic and/or inclined orbits . up to now , almost all treatments have been based on the separability of the equation in the frequency domain and calculate the energy and waveforms for the first few @xmath0 multipoles of the spheroidal harmonics expansion once the radial part has been dealt with .    using the frequency domain approach , quite a few detailed simulations of the gravitational waves emitted by a particle in a bound orbit around a black hole have been carried out in recent years . the first results for radiation emitted by a particle in orbit around a schwarzschild black hole were carried out by detweiler , who pioneered these techniques @xcite . in a series of six papers , poisson and various collaborators studied in detail the gravitational wave emission to infinity and into the horizon of a particle in circular orbit around a non - rotating hole @xcite . other groups have been doing extensive work on post - newtonian expansions for various types of particle orbits around kerr holes ( see @xcite and references therein ) . work has been done on eccentric orbits around schwarszchild @xcite , circular orbits slightly inclined away from the equatorial plane of kerr holes @xcite , and particles with spin orbiting kerr holes on equatorial circular orbits @xcite . ( 1993,1995 ) & circular orbits around kerr bh & infall trajectories with kojima & nakamura ( 1984 ) & ang mom @xmath1 0 & infall along symmetry sasaki & nakamura ( 1982 ) & axis of kerr bh & spinning particles in tanaka _ et al _ ( 1996 ) & circ equatorial orbits & slightly eccentric and hughes ( 2000 ) & non - equatorial orbits    for single static orbits ( ignoring the decay by radiation reaction ) to get good accuracies in the frequency domain approach involves summing for values of @xmath0 up to 12 for quite a few of the above - harmonic frequencies , requiring close to 3000 iterations using a numerical algorithm that involve complex green s function calculations @xcite . further discussion of the number of @xmath0 modes needed , and also a discussion of the self - force in the case of a scalar field can be seen in the paper by burko @xcite . we will perform the numerical evolution with the latest version of a code originally developed by laguna and collaborators @xcite to evolve the teukolsky equation in the time domain from a generic initial data slice . the method analyzes the radiation at `` infinity '' by dealing with the @xmath2 version of the equation . then the equation is rewritten as a first order matrix equation and numerically integrated . the code has been tested treating scalar fields , scattering of gravitational waves and analysis of quasi - normal ringing and power law tails of the outgoing radiation . we start with the original teukolsky equation @xcite before separation is performed , @xmath3 \frac{\partial^2 \psi}{\partial t^2 } + \frac{4 m a r}{\delta } \frac{\partial^2 \psi}{\partial t \partial \phi } + \left [ \frac{a^2}{\delta } - \frac{1}{\sin^2 \theta } \right ] \frac{\partial^2 \psi}{\partial \phi^2 } - \delta^{-s } \frac{\partial}{\partial r } \left ( \delta^{s+1}\frac{\partial \psi}{\partial r}\right ) \nonumber \\ & & - \frac{1}{\sin \theta } \frac{\partial}{\partial \theta } \ , \left ( \sin \theta \frac{\partial \psi}{\partial \theta } \right ) -2 s \left [ \frac{a ( r - m)}{\delta } + \frac{i \ , \cos \theta}{\sin^2 \theta } \right ] \frac{\partial \psi}{\partial \varphi } - 2 s \left [ \frac{m(a^2 - r^2)}{\delta } - r - i\ , a\ , \cos \theta \right ] \frac{\partial \psi}{\partial t } \nonumber \\ & & + \left [ s^2 \cot^2 \theta - s \right ] \psi = 4 \pi ( r^2+a^2\cos^2\theta ) t\end{aligned}\ ] ] where @xmath4 postulating an ansatz for the @xmath5 function of the form , @xmath6 to match this @xmath7 dependence , we write the matter source term as a similar fourier decomposition @xmath8    with a modicum of algebra one can rewrite the teukolsky equation as , @xmath9    this is possible with the following definitions @xmath10 @xmath11 + i \left [ 4 m a r m + 2 s a \delta \cos \theta \right ] \right ) , \ ] ]    @xmath12    @xmath13   + f \frac{(r^2 + a^2)}{\sigma^2 \delta^s } \frac{\partial^2 f}{\partial r^{*2}}\end{aligned}\ ] ]    @xmath14    @xmath15    where @xmath16 . one now introduces an auxiliary field @xmath17 that converts the teukolsky equation to a set of 2 coupled pdes which are first order in space and time @xmath18 @xmath19 the code we are using produces stable evolutions of these equations using a two - step lax - wendroff method .    as in other applications of this teukolsky code , we impose boundary conditions at the edges of the polar computational domain , i.e. at the black hole horizon , at the rotation axis and at the far end of the radial grid . the condition near the horizon is @xmath20 , due to the known asymptotic behavior of the fields @xcite . at the outer boundary we impose usual outgoing boundary conditions . since the wave equation has a potential , errors due to reflection from this boundary bounce back into the computational domain . we deal with this by making the computational domain in the radial direction large enough that any numerical reflections will not make it back to the point where we compute the waveforms and energy flux in the time allotted for the simulation . at the axis of rotation of the black hole one imposes either the condition @xmath21 or @xmath22 depending on the parity of the field specified by the azimuthal integer @xmath23 . the initial data we take for the fields is zero . since we are considering the teukolsky equation with a source we are therefore violating the constraints of general relativity on the initial data this would correspond to the particle `` appearing suddenly '' and therefore generates an artificial burst of radiation . to get more realistic initial data slices , one could try constructing initial data sets that correctly satisfy the hamiltonian and momentum constraints with prescriptions like that of bowen and york @xcite in the limit in which one black hole is very small , for instance , or modifications of particle limit data sets such as those of lousto and price @xcite . these prescriptions specify the @xmath24 sets for use in standard adm evolutions . to convert these quantities to the set @xmath25 needed for the teukolsky evolution , one can use the formulae developed in reference @xcite . we plan to explore the extent of the benefits that this more accurate intial data specification may bring to enhance the orbital simulations in the particle limit case future work . the next step is to construct the source term for a general geodesic trajectory of the particle moving in the background kerr geometry . to include the effect of radiation reaction one treats it as a slight deviation of the background geodesic orbital motion .    in the time domain teukolsky equation the main source term is given by @xmath26 where @xmath27 + ( \delta^ * - \tau^ * + \beta^ * + 3 \alpha + 4 \pi ) \nonumber \\ & & \times [ ( \delta + 2 \gamma + 2 \mu^*)t_{nm^ * } - ( \delta^ * - \tau^ * + 2 \beta^ * + 2 \alpha ) t_{nn}].\end{aligned}\ ] ]    all the differential operators and newman  penrose coefficients , as well as the choice of kinnersley tetrad vectors follow the notation of @xcite . by choosing their values in boyer - lindquist coordinates and expanding everything we can end up with an explicit but quite complicated form of the source term that is amenable for numerical calculations@xcite . we obtain the @xmath28 , @xmath29 , and @xmath30 by discretizing and integrating the geodesic equations in the kerr background , @xmath31 ^ 2- \delta\left[r^2 + ( l_z - a e)^2 + q\right]\;,\label{eq : rdot } \\ % \sigma^2\left({d\theta\over d\tau}\right)^2 & & = q - \cot^2\theta l_z^2 -a^2\cos^2\theta(1 - e^2)\;,\label{eq : thetadot } \\ % \sigma\left({d\phi\over d\tau}\right ) & & = \csc^2\theta l_z + ae\left({r^2+a^2\over\delta } - 1\right ) - { a^2l_z\over\delta}\;,\label{eq : phidot } \\ % \sigma\left({dt\over d\tau}\right ) & & = e\left[{(r^2+a^2)^2\over\delta } - a^2\sin^2\theta\right ] + al_z\left(1 - { r^2+a^2\over\delta}\right)\;.\label{eq : tdot } \label{eq : geodesiceqns}\end{aligned}\ ] ]    the source term is singular in those points where the particle is located at each timestep . we treat this by using the approximation in which we substitute the delta function by a very narrow gaussian function ( whose width is only of a few gridpoints ) : @xmath32    one must show that this substitution is adequate by showing convergence of the obtained waveform as @xmath33 , and one must take special care in normalizing the mass correctly as the particle moves in the background kerr spacetime so that integration over the whole gaussian density profile gives a constant mass throughout the whole trajectory and to satisfy @xmath34 . this translates in that the gaussian must approximately integrates to one on the spatial manifold . to achieve this with the correct three - metric of the background black hole spacetime , we normalize the mass - density of the particle by the factor , @xmath35 where @xmath36 is the 3-metric of flat space and @xmath37 is the 3-metric of a slice of constant boyer - lindquist time of the kerr spacetime . we then normalize the particle mass as it moves thru the space by multiplying the quantity @xmath38 in the source term expression by this factor @xmath39 . we handle the dirac deltas in the radial and angular direction through the gaussian approximation , whereas the @xmath40 function we handled analytically . the normalizations presented are common in the sph literature @xcite . we now want to calculate the energy flux emitted in the form of gravitational waves from the binary that an orbiting observatory such as lisa would be able to detect . in this paper we focus on calculating fluxes from circular equatorial orbits without the effects of radiation reaction mainly to test the viability of this approach and to show stability , convergence , and agreement with previous results of high accuracy frequency domain calculations@xcite . in the next section we will explain the advantages of working in the time domain for future calculations of eccentric and inclined particle orbits . the code monitors the value of @xmath5 at a set point in the grid far from the horizon and computes the energy and angular momentum fluxes in gravitational radiation according to the formulas @xcite : @xmath41\right\}\ ; \;\;\;\;\ , { \rm d}\omega=\sin\theta{\rm d}\theta{\rm d}\varphi . \label{jfl}\end{aligned}\ ] ]    the code produces output files listing ( a ) the real and imaginary parts of the @xmath5 field at representative locations of the 2-d grid , and ( b ) the computed energy and angular momentum fluxes for different values of @xmath23 . and all this is done periodically at preselected times along the length of the simulation for a circular orbital geodesic and from here the fluxes are postprocessed .    [ cfigure ]    we made standard convergence tests by assuming that the numerical values included errors that decrease as a power of the grid spacing @xmath42 . specifically we assume that the value of the teukolsky function @xmath43 at any point in the computational domain goes like @xmath44 grids typically used in this work are @xmath45 and @xmath46 and @xmath47 and @xmath48 . therefore , if we run at various resolutions ( l : 2000x20 , m : 4000x40 , h : 8000x80 ) we can get a quantity @xmath49 that measures the order of convergence _ a _ of the code . this is defined as , @xmath50 and results for it are depicted in figure 1 where we see second order convergence of the evolution code . we also verified that the results are independent of the gaussian width @xmath51 for both the radial and angular gaussian distributions used to approximate the location of the particle in the discrete spatial grid ( see figure 5 ) .    since equatorial orbits have been treated with great precision in the frequency domain method , it is important to show that our method gives comparable results if one looks at quantities like the average gravitational energy flux at large distances from the horizon . finn and thorne @xcite published a comprehensive work that tabulates and includes all previous results for high - precision gravitational wave flux and wave amplitudes coming from many kinds of representative compact objects in equatorial circular orbits around massive central black holes . in table [ compft ] we present representative comparisons with the results of formulas ( 3.8 ) and ( 3.10 ) of that paper using the precise relativistic corrections included in tables iii - vi that they also provide . @xmath52 & @xmath53 & @xmath54 & & @xmath55 & @xmath56 & @xmath57 & & @xmath58 & @xmath59 & @xmath60 & & @xmath61 & @xmath62 & @xmath63 we have shown that it is possible to integrate the teukolsky equation in the time domain for the problem of a particle orbiting around a black hole with the particle was represented as a gaussian and obtain results reasonably close to those of the frequency - domain calculations . the time domain approach is expected to become more useful for the case of more complicated orbits . we plan in future papers to discuss how to how to use better approximations to the gaussian like those of reference @xcite , and how to incorporate radiation - reaction effects . also , the small discrepancies with the thorne  finn results need to be further investigated to define well a `` window of reliability '' for the results of the code . another element to be incorporated is to take into account the radiation flux falling into the black hole , since this will also influence trajectories corrected by radiation - reaction . the first approach will be to correct the orbit trajectory using energy - balance arguments , as for instance in @xcite . further possibilities include incorporating proposed radiation reaction forces directly into the integration of the orbit . we are grateful to pablo laguna who continually provided help with the use of his teukolsky code during this work and for comments on the manuscript . we also thank eric poisson for providing results for comparison and comments . we are indebted to lior burko , scott hughes , carlos lousto and karl martel for comments on the manuscript . rla is grateful for support in form of a fipi ( institutional fund for starting research ) grant of the office of the dean of research and graduate studies of the rio piedras campus at the university of puerto rico . gk acknowledges research support from long island university . this work was supported by nsf grants phy-0140236 , phy-9800973 and by funds of the horace hearne institute for theoretical physics at lsu .
Abstract:

Answer: we consider the problem of the gravitational waves produced by a particle of negligible mass orbiting a kerr black hole . we treat the teukolsky perturbation equation in the time domain numerically as a 2 + 1 partial differential equation . we model the particle by smearing the singularities in the source term by the use of narrow gaussian distributions . we have been able to reproduce earlier results for equatorial circular orbits that were computed using the frequency domain formalism . the time domain approach is however geared for a more general evolution , for instance of nearly geodesic orbits under the effects of radiation reaction .


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 prototypical hollow nanoparticle is the buckyball c@xmath4 , which crystallizes into _ fullerite_. many other morphologies for carbon sheets have been found since the early 1990s , including fullerenes c@xmath5 of varying size , multi - walled fullerenes ( _ carbon onions _ ) , and single- and multi - walled carbon nanotubes . closed structures of carbon are formed in order to avoid dangling bonds at the edges of finite sized carbon sheets . this mechanism is generic for anisotropic layered material of finite size , and up to now more than 30 other materials have been prepared as hollow nanoparticles ( see the reviews @xcite and references therein ) . this includes the metal dichalcogenides mex@xmath0 ( me = w , se , x = s , se ) , bn , and gan . usually _ inorganic fullerenes _ like ws@xmath0 and mos@xmath0 are multi - walled . in fig . [ tem ] we show transmission electron micrographs of mos@xmath0-particles which were synthesized by solid - gas reaction . control of size and shape is rather difficult , as evidenced by the irregularly faceted shapes . however , the methods used for the synthesis of hollow nanoparticles are developing very rapidly , and it is to be expected that control of size , thickness and shape will become much better in the future .    [ cols="^,^ " , ]     hollow nanoparticles combine covalent in - plane strength with flexible out - of - plane bending of thin films , which results in high mechanical stability . they interact by van der waals ( vdw ) forces , which for example determine the material properties of fullerite @xcite and the phase behavior of buckyballs @xcite . vdw interactions become even more important for larger fullerenes , carbon onions and inorganic fullerenes ; these materials , therefore , feature a large degree of non - specific adhesion . in particular , vdw adhesion can lead to deformation of the hollow nanoparticles . for example , it has been found that carbon nanotubes adhering to each other or to flat surfaces show observable deformations for radii as small as 1 nm @xcite . in the presence of local surface features like crystal steps , adhering nanoparticles will have spacially varying electronic properties due to spacially varying deformation . hollow nanoparticles may lead to new applications in nanoelectronics and -optics . in contrast to filled nanoparticles , they can also be used for storage and delivery systems . they show superior mechanical properties , like high flexibility , high tensile strength and light weight , which will lead to applications as ultra - strong fibers . larger fibers can be produced using the cohesiveness provided by vdw interactions . inorganic fullerenes have also been shown to offer favorable tribological properties @xcite . detailed investigations with the surface force apparatus demonstrated that low friction and wear is caused by material transfer onto the sliding surfaces @xcite .    for these applications , a good theoretical understanding of the physical properties of hollow nanoparticles is important . here we discuss structural properties of spherical hollow nanoparticles in a continuum approach that allows us to focus on unusual generic properties which result mainly from geometrical effects . in particular , we address mechanical properties , van der waals interaction , thermodynamic behavior and deformation by adhesion . in this paper , we focus on the main results of our work ; more details can be found elsewhere @xcite . during recent years , mechanical properties of fullerenes and carbon nanotubes have been extensively studied at the level of molecular calculations , that is by first principles , tight binding and force field techniques . the advantage of molecular models is that they can provide detailed quantitative predictions . for larger systems like multi - walled nanoparticles , their implementation becomes difficult due to computer time requirements . here we use classical continuum elasticity theory , which is asymptotically correct for large systems and allows us to treat all these systems using the same framework . the large advantage of this approach is that it provides insight into the generic properties of these systems . molecular calculations for carbon nanotubes have shown that the predictions of continuum elasticity theory persist well into the limit of radii smaller than 1 nm @xcite . recently this observation could be verified also for bn and mos@xmath0 nanotubes @xcite . the layered material discussed here usually has a hexagonal lattice structure , for which elasticity theory @xcite predicts that there are only two in - plane elastic constants and the sheet is elastically isotropic . its deformation energy consists of stretching and bending terms . for single - walled carbon sheets , the in - plane stretching modulus @xmath6 erg/@xmath7 and the out - of - plane bending rigidity @xmath8 erg . for both mos@xmath0 and ws@xmath0 , @xmath9 is smaller by a factor @xmath10 and @xmath11 is larger by a factor @xmath12 . for multi - walled nanoparticles , the elastic constants scale with thickness @xmath2 as @xmath13 and @xmath14 , where @xmath15 is the largest in - plane elastic constant of the corresponding layered material @xcite . the values for @xmath15 are @xmath16 , @xmath17 and @xmath18 erg/@xmath19 for c , mos@xmath0 and ws@xmath0 , respectively . stretching can be avoided completely by bending the sheet into a cylindrical nanotube . however , this is not possible for spherical nanoparticles , which we discuss here . the same holds true for saddle - like structures , and in general for all surfaces with non - zero gaussian curvature . on a microscopic level , the requirement for non - zero gaussian curvature also means that a flat sheet can not be bent into a sphere without introducing defects . continuum elasticity approaches to fullerenes therefore explicitly treated the 12 pentagonal defects necessary to close a sheet of carbon hexagons @xcite . in order to avoid this complication , here we assume that the curvature generating defects are distributed in a homogeneous manner and result in a vanishing internal strain ( like for a ping - pong ball ) . in this case , the hollow nanospheres can be modeled as elastic shells with thickness @xmath2 and prefered radius @xmath1 . mechanical collapse in tribological applications is likely to be caused by either direct forces or high pressure . in the latter case , an inward buckling instability occurs , since the pressure energy scales with a higher power of indentation @xmath20 than the restoring elastic energy ( @xmath21 versus @xmath22 ) @xcite . mechanical collapse is expected when indentation @xmath20 becomes of the order of radius @xmath1 . the critical pressure for collapse can then be estimated to be @xmath23 . for buckyballs and typical inorganic fullerenes , this yields values of 15 and 2 gpa , respectively . for nanoparticles , the corresponding forces are in the nn - range . since vdw interactions are very difficult to treat on a molecular level , a continuum approach is even more useful in this case . it is well known from colloid science that vdw interactions strongly depend on the geometry of the system under investigation . geometrical aspects of vdw interactions are well treated in the hamaker approach , in which an atomic lennard - jones interaction is summed in a pairwise manner over the relevant geometry and the molecular details of the interaction are lumped into the effective lennard - jones parameters @xmath24 and @xmath25 . for carbon , @xmath26 and @xmath27 erg @xcite . the corresponding hamaker constant @xmath28 is in the order of @xmath29 erg . for two buckyballs , the hamaker approach for thin sheets leads to the _ girifalco potential _ @xcite , which is considerably shorter ranged than a lennard - jones potential . since a small range of attraction leads to the disappearance of the gas - liquid coexistence from the phase diagram @xcite , it has been suggested that buckyballs constitute the first non - colloidal system for which the gas - liquid critical point does not appear in the phase diagram . however , large scale simulations have shown that buckyballs indeed do show a gas - liquid coexistence , which extends roughly from 1900 to 2000@xmath30c @xcite .    for hollow nanoparticles in general , the effect of a finite thickness @xmath2 becomes important . however , for nanoparticles in a crystal or adhering to a flat substrate or to each other , the distance of closest approach @xmath31 between interacting particles is on an atomic scale . under such conditions , the vdw interaction saturates on the atomic scale @xmath31 with increasing thickness @xmath2 . for example , the surface energy for the interaction with a halfspace scales as @xmath32 and @xmath33 for thin films and halfspaces , respectively , with the crossover occuring for @xmath34 . since @xmath35 in our case , the surface energy is @xmath33 independently of @xmath2 and no significant scaling with thickness is expected .    for a more detailed investigation , one has to consider the general case of two nanospheres of finite thickness interacting in a pairwise manner by a lennard - jones potential @xcite . the interaction energy can be written as @xmath36 , where @xmath37 is the interaction energy between two filled balls of radii @xmath38 and @xmath39 , respectively , which can be calculated analytically for a lennard - jones potential . in order to predict phase behavior , three quantities characterizing the resulting potential are especially important : the hard core diameter @xmath40 sets the density scale , the potential depth sets the temperature scale , and the interaction range determines if a gas - liquid coexistence occurs ( compare fig . [ tem ] ) . our full analysis has shown that for hollow nanoparticles of radius @xmath1 and thickness @xmath2 , the potential depth and interaction range scale linearly and inversely with @xmath1 . both quantities depend only weakly on @xmath2 . in fig . [ interactionrange ] we show our numerical results for the interaction range as a function of @xmath1 and @xmath2 . recent work for the double - yukawa potential suggests that the gas - liquid coexistence disappears for an interaction range smaller than @xmath41 @xcite . our analysis predicts that this happens for @xmath42 and @xmath43 for single- and multi - walled nanoparticles , respectively . for larger @xmath1 , the critical point will disappear from the phase diagram . we consider the case where a spherical hollow nanoparticle is indented a distance @xmath20 due to adhesion to a flat substrate . the vdw energy gained on adhesion can be estimated to be @xmath44 where @xmath28 is the hamaker constant and @xmath31 an atomic cutoff @xcite . the energy of deformation can be estimated to be @xmath45 . with increasing @xmath1 , deformation becomes more likely , since both the vdw energy and deformability increase . however , deformability decreases rapidly with increasing thickness @xmath2 , since @xmath46 . solving for the indentation @xmath20 and using the values given above , it is found that for hollow carbon nanoparticles with few walls , @xmath20 can be in the nm - range @xcite . however , for typical inorganic fullerenes deformations are suppressed by the larger shell thickness @xmath2 . finally we briefly discuss a nanotube adhering to a flat substrate . as explained above , nanotubes can be bent without being stretched . on deformation , one has to pay a bending energy per length of @xmath47 . the gain in vdw energy scales as @xmath48 . it follows that for carbon nanotubes with small thickness @xmath2 , the indentation @xmath20 can be in the nm - range for even smaller values of radius @xmath1 than for hollow carbon nanoparticles . it is a pleasure to thank r. tenne for many stimulating discussions . uss was supported by the minerva foundation and sk by the ministry of education , science and culture of japan . this estimate is valid for weak adhesion . in @xcite , we used an estimate for strong adhesion , @xmath49 . then the predicted deformation is larger by a factor of @xmath50 , yet the general conclusions are the same .
Abstract:

Answer: when sheets of layered material like c , ws@xmath0 or bn are restricted to finite sizes , they generally form single- and multi - walled hollow nanoparticles in order to avoid dangling bonds . using continuum approaches to model elastic deformation and van der waals interactions of spherical nanoparticles , we predict the variation of mechanical stability , adhesive properties and phase behavior with radius @xmath1 and thickness @xmath2 . we find that mechanical stability is limited by forces in the nn range and pressures in the gpa range . adhesion energies scale linearly with @xmath1 , but depend only weakly on @xmath2 . deformation due to van der waals adhesion occurs for single - walled particles for radii of few nm , but is quickly suppressed for increasing thickness . as @xmath1 is increased , the gas - liquid coexistence disappears from the phase diagram for particle radii in the range of 1 - 3 nm ( depending on wall thickness ) since the interaction range decreases like @xmath3 .


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: increasing interest is attracted by modified versions of general relativity @xcite . they have been proposed as serious alternatives to einstein s theory of gravitation , and could be used to describe more accurately the observed accelerated expansion of our universe @xcite . in addition , it has been shown @xcite that it is actually possible to reconstruct the explicit form of the postulated curvature function @xmath0 , from the universe expansion history . it is quite well known that modified versions of general relativity are mathematically equivalent to scalar fields models ( see e.g. @xcite ) , meaning that a solution in a modified gravity model can always be mapped into a solution of the corresponding scalar field theory . in spite of this mathematical correspondence , physical equivalence does not always follow . in fact , two corresponding solutions of two equivalent theories can actually exhibit rather different physical behaviors . furthermore , it is not necessary , in order to justify modified gravity , to do it by always using this relation with scalar field theories . because of the new situation , in the following we will disregard this mathematical equivalence and do consider in our analysis modified gravity as an independent theory aiming directly at some measurable physical properties . what is more , in our treatment modified gravity will in fact be viewed just as a different classical theory of gravitation . although other models have been considered @xcite with the gauss - bonnet scalar in the action , here we shall restrict our attention to pure @xmath0 models . as often discussed , there are limitations on the function @xmath0 when trying to construct a theory which is in agreement with the very precise solar system tests carried out so far , as well as with all the known cosmological bounds @xcite . recently , different models of that kind have been studied @xcite , the last three of them having been reported to pass all solar system tests ; in addition , they exhibit a number of very interesting features . in @xcite possible newton law corrections to such models have been considered . we present here an analysis of the models above based on the study of pure schwarzschild and also schwarzschild - de sitter black holes ( sbh , sdsbh ) , calculated through use of the noether charge method . we start with a discussion of two examples considered in @xcite and go over to study more recent ones @xcite . the direct confrontation of basic quantities , as the black hole entropy , with the well established , classical ( einsteinian ) result can offer a further insight into the construction of a general @xmath0 theory . we start with a short review of modified @xmath0 gravity and the noether charge method to compute the bh entropy and then calculate the bh entropy for the models in refs .  @xcite . after that , we extend our analysis to the models @xcite , taking due care of the sign of the bh entropy and also discussing the stability conditions as well as the existence of a schwarzschild bh solution . we conclude by providing a brief comparison of the different models considered . the action for @xmath0 gravity ( see e.g. @xcite for a review ) is @xmath1 with @xmath0 a generic function . as discussed in @xcite , in order to give rise to a realistic cosmology this function needs to fulfill some limiting conditions . however , for the moment we can ignore them , because they do not affect the considerations which will follow . the equations of motion for this theory , in the presence of matter , are @xmath2 where @xmath3 is the matter stress energy tensor . contracting the indices in the last equation , we obtain the relation @xmath4 since we are interested in the study of the schwarzschild solution ( in the case when there is a @xmath5 solution ) or either in schwarzschild - de sitter black holes , we restrict our reasoning to the case of metric tensors with constant scalar curvature in the vacuum . in that case , we simply have @xcite @xmath6 for completeness , let us recall that , in order to build a realistic modified gravity , @xmath0 needs satisfy the two conditions : @xmath7 the first condition corresponds to the existence of an effective cosmological constant at high curvature . the second one allows for vacuum solutions , as for example minkowski or schwarzschild space - times . then , although an effective cosmological constant exists , vacuum solution are preserved and it is legitimate to study them also in a large scale universe with nonzero @xmath8 . moreover , in order to give rise to stable solutions @xcite , the following additional condition needs to be fulfilled : @xmath9/f''(r ) > r.\ ] ] we can now consider the schwarzschild - de sitter metric , a spherically symmetric solution of ( [ eq2 ] ) with constant curvature @xmath10 ( see @xcite ) @xmath11 where @xmath12 and @xmath13 . we will use the noether charge method , as discussed in @xcite , in order to calculate the entropy for the schwarzschild - de sitter bh . the entropy formula reads @xmath14 and the integration is made on the external horizon of events surface . in the case of constant curvature , for a generic modified theory , the result is @xmath15 a_h/4g,\ ] ] where @xmath16 is the area of the bh horizon . this enables us to calculate the bh entropy for a generic @xmath0 theory . we must here stress the fact that the requirement of positive black hole entropy simply avoids the appearance of ghost or tachyon fields in the corresponding scalar field theory . then a negative entropy is simply a footprint of some instabilities in the einstein frame . what is new in this picture is that we do not need to involve the ( mathematical ) equivalence of these models in order to give a physically meaningful interpretation of such constrain . in order to illustrate the method with explicit examples of entropy calculation , we analyze here two modified gravity models that appeared some time ago and which have been quite successful up to now . in those models no @xmath5 solution occurs . the first one , introduced in @xcite , is given by @xmath17 with @xmath18 . the condition to obtain a sdsbh ( namely , @xmath19 ) leads to @xmath20 \nonumber \\ & & \hspace*{10 mm } = 2 \left [ - a\left(r_0 - \lambda_1 \right)^{-n } + b \left ( r_0 - \lambda_2 \right)^m \right],\end{aligned}\ ] ] so that we have for the entropy @xmath21.\ ] ] thus the sdsbh entropy is positive for all @xmath22 . the second model , studied in @xcite , is defined by @xmath23 this modified gravity model does not admit vacuum solutions , thus we can calculate the entropy for the sdsbh . the ricci scalar is such that it satisfies the relation @xmath24 the sdsbh entropy is given by @xmath25 and turns out to be positive for all values of @xmath26 . we now analyze three recent models @xcite which have been proven to comply with the solar - system as well as with other cosmological parameter constraints . their respective authors have given a complete discussion of each model , taking care to provide a range for the free parameters contained in the @xmath0 function , and have also produced stable solutions . here we just want to stress , with the help of these examples , how the corresponding bh entropy calculation offers a further tool in order to confront each of those modified gravity theories with einstein s general relativity , given the fact that the presence of spherically symmetric bh solutions is a necessary element of all local tests . in this model @xcite ( with @xmath27 ) we have @xmath28 in @xcite , @xmath29 is chosen such that , at cosmological scale , @xmath30 at the present epoch , and @xmath0 satisfies the condition @xmath31 for @xmath30 . this also ensure that solutions with @xmath30 are stable . moreover the requirement that @xmath32 at fixed @xmath33 gives a cosmological constant , in both cosmological and local tests of gravity . in spite of this fact , since @xmath34 , this theory admits the schwarzschild solution ( i.e. @xmath5 ) . by the way we note that the stability condition for the vacuum solution is not satisfied unless @xmath35 and @xmath36 . therefore , except of this case , vacuum solutions ( than also sbh ) are unstable . note also that @xmath35 corresponds to a lagrangian @xmath37 for small @xmath8 , so it is associated with a correction to the gravitational coupling constant for small @xmath8 , giving an effective @xmath38 ( see @xcite ) . the entropy formula gives for the sds metric @xmath39 the entropy for the schwarzschild solution is @xmath40 then , in the only stable case , with @xmath35 and @xmath36 , a correction to the classical einstenian bh entropy is found . from ( [ entropyhu ] ) it also follows that , for sds bh with @xmath41 , the entropy is positive and corrections to its einstenian value are of order @xmath42 .      in this model @xmath0 is @xcite @xmath43^{-n } -1 \right\},\ ] ] from where @xmath44^{-n-1}.\ ] ] note that in this case @xmath45 an thus , although this model admits a sbh solution , it is unstable together with all its vacuum solutions . in @xcite the author limits his analysis to solutions that satisfy the following stability conditions @xmath46 we can therefore consider the sdsbh solutions , with curvature @xmath10 given by @xmath47^{-n-1 } + 1 \right\}\nonumber \\ & & \hspace*{10 mm }   = - 2 \lambda c \left [ \left ( 1+(r_0/c)^2 \right)^{-n } -1 \right],\end{aligned}\ ] ] that satisfies also ( [ stabilitistarobinski ] ) . in this case , the entropy is just @xmath48^{-n-1 }      \right\}.\ ] ] therefore , in this case a _ non trivial _ correction of the sdsbh entropy is found . we just stress the fact that the sbh solutions have classical entropy but are _ unstable _ , and that the sdsbh ones have a modified entropy which , under the limitations stated in @xcite , is strictly positive .      here @xmath0 is given by @xcite @xmath49/2a\ ] ] and @xmath50/2.\ ] ] this model admits a sbh solution . the entropy for the sdsbh is simply @xmath51 a_h/8g.\ ] ] we can use the stability condition given in @xcite , @xmath52 , to obtain @xmath53 for the sbh , the stability condition is just @xmath54 . moreover , being @xmath31 for all @xmath8 , in this model the vacuum solutions are always stable and there are no substantial corrections to the classical result .      it is interesting to put together all three models and explicitly compare the behavior of the function @xmath0 , in particular , the stability of the euclidean limit and the asymptotic behavior at large curvature . to simplify the comparison , we do not play with the values of the different parameters and set all coefficients equal to 1 and the curvature powers equal to 2 or 4 . for the case of the hu - sawicki model , with @xmath55 being @xmath56 , the corresponding plot is given in fig .  1 .    for the case of the starobinsky model , with @xmath57 being @xmath58 , we obtain fig .  2 .    and for the case of the appleby - battye model , with @xmath59 being @xmath60 and @xmath61 , fig .  3 . in a first comparison of these different models , we note that the one of starobinsky , as remarked by the author himself @xcite , has unstable vacuum solutions . this seems true also for the hu - sawicki model , except for the case when @xmath62 , that leads to a negative and infinite bh entropy . the appleby - battye model has the important property to possess stable vacuum solutions for a suitable range of the free parameters . it also yields an unmodified expression of the bh entropy . very recently , a modified gravity model has been published @xcite , that unifies inflation and cosmic acceleration under the same picture and also complies with the solar system tests . it is @xmath63 dx,\ ] ] where @xmath64 and @xmath65 is a constant given by @xmath66 @xmath67 and @xmath68 is the ricci scalar at present . the effective cosmological constant in the early universe is simply @xmath69 and the present cosmological constant is @xmath70 . because of the fact that @xmath34 , this model allows for sbh . to be general , we first calculate the sdsbh entropy , which is given by @xmath71 \right\},\ ] ] where @xmath10 is the sdsbh curvature that fulfills condition ( 4 ) . thus , black holes are less entropic than in einstein s theory . for the sbh , we have @xmath72 note that the stability condition is here @xmath31 , thus @xmath73 to have stability for sbh , we need that @xmath74 this can be considered , together with the condition @xmath75 stated in @xcite , to avoid newton law corrections in our solar system and on the earth surroundings . comparison of bh entropy in modified gravity theories and in the usual einstenian gravity case have been carried out . we have here analyzed different suitable models , recently considered in the literature , and have shown explicitly how corrections to the ` classical ' bh entropy can in fact appear . we have also argued that the condition of positive entropy can be used as an extra condition in order to constrain the viability of modified gravity theories . of course this conditions is equivalent to the requirement that neither ghost nor tachyon fields appear in the equivalent scalar field models . anyhow , if referred to the bh entropy , this condition has a direct interpretation in the framework of modified gravity , without needing to pass through the ( mathematically but not physically equivalent ) scalar field theories . we hope that this quite simple considerations may be useful for future analysis . we thank sergei odintsov for helpful discussions . this investigation is partly based on work done by ee while on leave at the department of physics and astronomy , dartmouth college , 6127 wilder laboratory , hanover , nh 03755 , usa , and was completed during a stay of fb in barcelona . it has been supported by mec ( spain ) , projects fis2006 - 02842 , and by agaur 2007be-1003 and contract 2005sgr-00790 . s. nojiri and s.d . odintsov , int . mod . phys . * 4 * , 115 ( 2007 ) . s. nojiri and s.d . phys . rev . * d74 * , 086005 ( 2006 ) ; s. capozziello , s. nojiri , s.d . odintsov and a. troisi , phys . lett . * b639 * , 135 ( 2006 ) . f. faraoni , arxiv : gr - qc/0607116 ; arxiv : gr - qc/0511094 ; a. cruz - dombriz and a. dobado , arxiv : gr - qc/0607118 ; n. poplawski , gr - qc/0610133 ; a. brookfield , c. van de bruck and l. hall , arxiv : hep - th/0608015 ; t. sotiriou and s. liberati , arxiv : gr - qc/0604006 ; a. bustelo and d. barraco , arxiv : gr - qc/0611149 ; g. olmo , arxiv : gr - qc/0612047 ; f. briscese , e. elizalde , s. nojiri and s.d . odintsov , phys . lett . * b646 * , 105 ( 2007 ) ; b. li and j. barrow , arxiv : gr - qc/0701111 ; o. bertolami , c. boehmer , t. harko and f. lobo , arxiv:0704.1733 ; s. carloni , a. troisi and p. dunsby , arxiv:0706.0452 . s. nojiri and s. odintsov , phys . lett . * b599 * , 137 ( 2004 ) ; m. abdalla , s. nojiri and s.d . odintsov , class . * 22 * , l35 ( 2005 ) ; g. cognola , e. elizalde , s. nojiri , s. d. odintsov and s. zerbini , jcap * 0502 * , 010 ( 2005 ) ; phys . rev . * d73 * , 084007 ( 2006 ) ; d.a . easson , int . a19 * , 5343 ( 2004 ) ; s. capozziello , v. cardone and a. troisi , phys . rev . * d71 * , 043503 ( 2005 ) ; g. allemandi , a. borowiec , m. francaviglia and s.d . odintsov , phys . rev . * d72 * , 063505 ( 2005 ) ; j. kainulainen , m. juvela and j. alves , arxiv : astro - ph/0703607 ; j.a.r . cembranos , phys . rev . * d73 * , 064029 ( 2006 ) ; t. koivisto and h. kurki - suonio , arxiv : astro - ph/0509422 ; t. clifton and j. barrow , arxiv : gr - qc/0509059 ; o. mena , j. santiago and j. weller , arxiv : astro - ph/0510453 ; i. brevik , arxiv : gr - qc/0601100 ; r.woodard , arxiv : astro - ph/0601672 ; t. koivisto , arxiv : astro - ph/0602031 ; g. cognola , m. castaldi and s. zerbini , arxiv : gr - qc/0701138 ; s. nojiri , s.d . odintsov and p. tretyakov , arxiv:0704.2520[hep - th ] ; l. amendola and s. tsujikawa , arxiv:0705.0396[astro - ph ] ; k. uddin , j. lidsey and r. tavakol , arxiv:0705.0232[gr - qc ] . s. capozziello , int . phys . * d11 * , 483 ( 2002 ) ; s. capozziello , s. carloni and a. troisi , rec . res . astroph .- rsp / aa/21 - 2003 ; s.m . carroll , v. duvvuri , m. trodden and s. turner , phys . rev . * d70 * , 043528 ( 2004 ) .
Abstract:

Answer: * abstract *    an analysis of some modified gravity models , based on the study of pure schwarzschild and of schwarzschild - de sitter black holes , and involving the use of the noether charge method , is carried out . corrections to the classical einsteinian black hole entropy appear . it is shown explicitly how the condition of positive entropy can be used in order to constrain the viability of modified gravity theories .