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: tidal theories have gained great attention in the last decade due to the large number of discovered close - in planets which orbit their parent stars with periods of a few days . according to classical theories ( e.g. darwin 1880 ; kaula 1964 ; mignard 1979 ) , the long - term tidal evolution of these planets results in orbital decay and circularization , on time - scales which depend on the physical properties of the interacting systems ( hut 1981 ; dobbs - dixon , lin & mardling 2004 ; ferraz - mello , rodrguez & hussmann 2008 ; jackson , barnes & greenberg 2009 ; rodrguez & ferraz - mello 2010 ) .    the rotation of short - period planets is also affected by tidal interactions . it is nowadays accepted that the spin rate of planets soon after their formation was much higher than their orbital mean motion ( e.g. agnor , canup & levison 1999 ; kokubo & ida 2007 ) . however , the cumulative tidal effects causes the planet to spin down , reaching its stationary value in time - scales much shorter than those needed to circularize the planet orbit . this evolution ultimately results in the synchronous rotation ( the spin state in which the orbital and rotation periods coincide ) , which occurs when the planet orbit is circularized ( hut 1981 ; ferraz - mello et al . 2008 ) . for non - circular orbits , the synchronous rotation only occurs if an additional torque acts to counteract the tidal torque . the additional gravitational torque may be generated by planets with a rocky - like composition having a permanent equatorial deformation or frozen - in figure ( see goldreich & peale 1966 ; greenberg & weidenschilling 1984 ) . these torques , in contrast to the tidal ones , are conservative . in the context of the two - body problem , the gravitational interaction of the star with an asymmetric planet drive the planet rotation into different regimes of motion , including oscillations around exact spin - orbit resonances ( goldreich & peale 1966 ) . when dissipative tidal effects are taken into account , the oscillations are damped and the planet rotation is trapped into a resonance . indeed , the rotational motion of some solar system bodies can be explained by the combined effects of the tidal and asymmetric torques ( goldreich & peale 1966 ; greenberg & weidenschilling 1984 ; correia & laskar 2009 ) . mercury is an example of a rocky planet with a rotation state currently evolving in the 3/2 spin - orbit resonance ( pettengill & dyce 1965 ) . several models of the secular spin evolution of mercury have been developed , including i ) the chaotic evolution of its orbit ( correia & laskar 2004 ) ; ii ) effects of the core - mantle friction ( correia 2006 ; correia & laskar 2009 ) ; iii ) the evolution starting with retrograde and subsequent synchronous motion combined with impact events ( correia & laskar , 2012 ; wieczorek et al . 2012 ) ; iv ) the consideration of non - linear tidal theories ( makarov 2012 ) . however , despite the large amount of previous works , the question on how the current resonant state of mercury was achieved still remains under discussion . concerning the rotation of earth - like close - in extrasolar planets , a near synchronous motion is usually assumed in view of tidal evolution . however , correia , levrard & laskar ( 2008 ) have shown that when an additional torque is considered , the equilibrium rotation of tidally evolved planets can assume rotational states which depart from synchronism . moreover , as the capture probability in a non - synchronous spin - orbit resonance requires a non - circular orbit ( goldreich & peale 1966 ; celletti & chierchia 2008 ) , the rotation of several recently discovered super - earth planets in eccentric orbits can be currently evolving under a resonant motion , provided the planet figure can support a permanent equatorial deformation . as the eccentricity is tidally damped , several resonant trappings would become unstable and the planet rotation ultimately reaches synchronization for almost circular orbits . however , the locking in a non - synchronous resonance with the very short orbital period can protect certain planet s longitudes from solar radiation , providing more favorable conditions with respect to habitability as compared to non - resonant rotation ( dobrovolskis 2007 ) . moreover , heller , leconte & barnes ( 2011 ) studied the equilibrium rotation of super - earths and the implications on their habitability . they found that , under assumption of an evolution driven by tidal torques , the equilibrium value of the rotation period of gj 581 * d * is probably about two times smaller than the orbital period .    in a recent work , callegari & rodrguez ( 2012 ) have analysed the rotation of extrasolar planets with orbital periods less than 33 days and masses below the mass of neptune . depending on one parameter related to the equatorial ellipticity and also on the orbital eccentricity , they identified regions of chaotic rotation in the phase space of some particular planets like 55 cnc**e * * , gj 1214**b * * , gliese 876**d * * , among others . the present investigation aims at two main aspects : i ) how the orbital evolution of the planet , due to tidal interaction with the central star , is altered when the asymmetric torque is considered ; and ii ) to analyse the occurrence of captures in spin - orbit resonant motions under combined effects of both torques for the planets kepler-10**b * * , 55 cnc**e * * and gj 3634**b**. these planets are examples of hot super - earths , which are close - in planets having a few earth - masses with a rocky - like composition . we are interested in the past behavior in the history of the rotation and the subsequent evolution to synchronization . we also discuss the possibility of a current evolution of the rotation of these planets into a super - synchronous resonance state . moreover , we investigate the implications on the orbital evolution due to such temporary resonant trappings .    despite the continuously growing number of discovered planets , their rotation state and equilibrium figures remain unknown . in this work , we perform several numerical simulations of the evolution of the systems , assuming several values of unknown parameters , such as equatorial ellipticity and initial rotation . the results allow us to understand the past history of these systems , their current orbital configuration and final rotational states . it should be noted that the majority of previous works concerning spin - orbit interactions with tidal dissipation have been addressed in the context of a fixed keplerian orbit , with constant values of semi - major axis and eccentricity ( e.g. , celletti , froeschl & lega 2007 ; celletti & chierchia 2008 ) . this assumption is adequate for solar system bodies , where tidal interactions are weak and , consequently , variations of orbital elements are small .    on the other hand , the existence of extrasolar planets with orbital periods of a few days indicates that their current orbital configurations would be the consequence of strong tidal interaction with the host star . in this case , the keplerian approximation is no longer valid . thus , we analyse spin - orbit interactions as the orbit evolves due to tidal effects . solving the exact equations of motion of the problem , we show that the forces arising from the potential of a asymmetric planet also contributes to the orbital variation , accelerating the orbital decay and the eccentricity damping , but only in the cases in which the rotation is locked in resonance .    the spin - orbit problem can also be modeled by the use of averaged equations of motion , where only terms with secular and resonant arguments are retained ( e.g. celletti & chierchia 2008 ; correia & laskar 2009 ) . through this model , we compute the averaged variations of semi - major axis and eccentricity due to the combination of asymmetric and tidal forces . we also give an explicit expression for the critical value of the equatorial prolateness , for which a given spin - orbit resonance would become unstable . the study of the particular cases of 2/1 , 3/2 and 1/1 spin - orbit resonances indicates a good agreement with the results of the numerical simulations of the exact equations . this paper is structured as follows . in section [ model ] , we present the exact equations of motion of the spin - orbit problem . section [ num ] shows the results of the numerical simulations , illustrating the evolutions of kepler-10**b * * , gj 3634**b * * and cnc 55**e**. we show that , assuming a nonzero value of the equatorial ellipticity , different temporary resonant trappings of the rotation can occur , with implications in the orbital decay and circularization . an analytical approach in the context of the averaged equations is shown in section [ mean ] , explaining the results obtained in the numerical simulations . finally , section [ discussion ] is devoted to discussion and conclusions . with the reference axis @xmath0 , which is assumed coincident with the orbital pericenter direction . the angle between the long axis and the star - planet direction is given by @xmath1 , where @xmath2 is the true anomaly of the planet . ] we consider a planet orbiting a central star with a period of a few days . hereafter , we refer to @xmath3 and @xmath4 as the planet and star masses , respectively ; @xmath5 and @xmath6 as the planet radius and the instantaneous star - planet distance , respectively . the elements @xmath7 stand for the semi - major axis , eccentricity and mean - motion of the astrocentric orbit of the planet , respectively . we suppose that the rotating planet body can be approximated by a homogeneous triaxial ellipsoid , with the equatorial axes @xmath8 and @xmath9 and the rotation axis @xmath10 . we assume that @xmath11 and @xmath10 is normal to the orbital plane ( i.e. , no obliquity is considered ; planar problem ) . in addition , we assume that the equatorial asymmetry ( @xmath12 ) corresponds to a frozen - in body or permanent deformation of the planet . the assumption above can be supported under consideration of a rocky composition of the planet body . this is the case of super - earth planets , which are the subject of the present investigation . the gravitational potential of the planet at the distance @xmath6 from its centre is given by ( e.g. beutler 2005 , p. 102 ) @xmath13 where @xmath14 is the gravitational constant and the angles , which has not been included in the present paper ( see ferraz - mello et al . @xmath2 and @xmath15 are defined in figure [ elips ] . note that @xmath16 is the angular rotation velocity of the planet . the parameter @xmath17 describes the equatorial ellipticity of the planet s gravity field and is related to the principal momenta of inertia @xmath18 through @xmath19 , where @xmath20 is the structure constant given by @xmath21 and @xmath22 is the polar moment of inertia . we neglect the contribution of terms arising from the oblateness of the rotating body , that is , @xmath23 and higher order terms . in fact , in the planar problem , the gravitational potential corresponding to the @xmath23term depends only on the distance @xmath6 . in this case , the potential generates a radial force with zero torque that does not affect the rotation evolution .    at the distance @xmath6 from the centre of the planet , the gravitational field @xmath24 gives rise to a force on a point mass @xmath4 which is given by adopted in beutler ( 2005 ) , the corresponding force is computed as the positive gradient of the potential . ] @xmath25 where @xmath26 and @xmath27 are unit vectors with origin in the star ( see fig . [ elips ] ) . the force @xmath28 alters the orbital component of the angular momentum , whereas the reaction force governs the evolution of the spin through the creation of the torque on the planet given by @xmath29 where @xmath30 and @xmath31 . thus , @xmath32 note that , for @xmath33 and @xmath19 , the above equation agrees with the classical expression for the torque ( e.g. danby 1962 ; goldreich & peale 1966 ) . the effects of the dynamical tides are raised by the gravitational action of the star on the inelastically deformed rotating planet . due to the internal viscosity of the body , there exists a delay @xmath34 ( time lag ) between the star perturbation and the maximum deformation . the relationship between @xmath34 and tidal frequencies ( which are linear combinations of @xmath35 and @xmath36 ) depends on the internal rheology , which is poorly known even for solar system bodies . for sake of simplicity , in this work , we adopt a linear model with constant @xmath34 ( darwin 1879 ; mignard 1979 ; hut 1981 ; rodrguez et al . 2011a ) ; for a review of other tidal models , the reader is referred to efroimsky & williams ( 2009 ) and ferraz - mello ( 2012 ) .    according to mignard ( 1979 ) , the tidal force and torque are given by @xmath37,\ ] ] and @xmath38,\ ] ] where @xmath39 and @xmath40 is the second degree love number . in this paper , we only consider the tides raised by the star on the planet , neglecting the tides raised on the star by the planet . indeed , the strength of stellar tide is proportional to the planet mass and can be safely neglected in the case of super - earths ( see rodrguez & ferraz - mello 2010 ) . note however that , after circularization of the planet orbit , stellar tides may be the only source of orbital decay .    in the reference frame centered in the star ( see figure [ elips ] ) , the exact equations of motion of the planet can be now written as @xmath41 the perturbing forces @xmath28 and @xmath42 alter the keplerian orbit of the planet which is defined by the first term of the first equation , while the corresponding torques in the second equation produce variations of the planet rotation . note that , in the planar case , this last equation can be written in scalar form , because both torques are normal to the orbital plane . in this section we apply the analytical model described in the previous section , in order to investigate the motion of three known super - earths planets , namely kepler-10**b * * , 55 cnc**e * * and gj 3634**b**. for this task , the exact equations of motion ( [ mov1 ] ) were integrated numerically using the ra15 code ( everhart 1985 ) . we stress that no scaling parameter , often used to accelerate the tidal evolution , was used in our simulations ( see , for instance , tittemore & wisdom ( 1988 ) , for the validity of the usage of the scaling parameter ) . table 1 shows the physical parameters and current orbital elements of the systems considered in this paper . [ tab1 ]    .physical parameters and current orbital elements of the analysed systems ( @xmath43batalha et al . 2011 ; @xmath44gillon et al . 2012 ; @xmath45bonfils et al . 2011 ) . the first two planets were discovered through the transit technique , allowing us a good estimation of their radii . @xmath46 and @xmath47 stand for mass and radius of the earth , respectively . [ cols="^,^,^,^,^,^",options="header " , ]      numerical integrations require a setting of the initial configurations of the systems . in the case of problems including dissipative processes , the choice is complicated , mainly , due to the inadequacy of the simulations of the past evolution of dissipative systems through integrations backward in time ( see details in michtchenko & rodrguez 2011 ) . however , as shown in that paper , the final state of a migrating system is independent of its initial position in the phase space , defined by three values of parameters : the mass and semi - major axes ratios and the total angular momentum . as a consequence , to assess the orbit configuration in the past , we have to set solely the starting semi - major axis value . the initial eccentricity should satisfy the conservation of the total angular momentum of the system . the total angular momentum of the system is composed of the orbital component written as @xmath48 , and the rotational component given by @xmath49 . using kepler s third law , it is easy to show that the ratio of rotational and orbital components is of the order of @xmath50 ( see rodrguez et al . 2011a ) ; since @xmath51 , the contribution of rotation to the total angular momentum can be neglected . therefore , considering that the orbital angular momentum is nearly conserved during tidal migration , we obtain that @xmath52 where the subscripts  ini " and  current " stand for initial and current orbital elements , respectively . in this way , choosing an initial value for the semi - major axis and knowing the current values of @xmath53 and @xmath54 , equation ( [ am ] ) allows us to constrain the initial value of the eccentricity ( see next section and specially figure [ am - fig ] ) . it should be emphasized that the value obtained in such a way must be considered as an approximation of the initial value of the eccentricity , since the planets may belong to systems with two or more planets . the value of @xmath34 is determined by the relationship between the tidal frequency @xmath55 and the quality factor @xmath56 . in the case of a linear approximation , we have @xmath57 , where @xmath58const . it is known that , in the fourier spectrum of the tidal potential , the most important frequencies ( that is , those involved in terms of order zero and order one in the eccentricity ) are @xmath59 , @xmath60 and @xmath61 . the first one is associated to the term which is independent of the eccentricity , while the two last frequencies correspond to terms which are proportional to @xmath54 ( darwin , 1879 ; ferraz - mello et al . 2008 ) . using the frequency @xmath62 and a modified version of the quality factor written in terms of the love number @xmath40 as @xmath63 , we obtain @xmath64    actually , the spin rotation of exoplanets is not available from observations , thus we choose its initial value arbitrarily but considering a super - synchronous motion ( @xmath65 ) .    finally , according to callegari & rodrguez ( 2012 ) , the equilibrium shape attained by the planet due to the gravitational interaction with the star , results in a prolateness , or equatorial ellipticity , which can be computed as @xmath66 . knowing that @xmath19 we can obtain the magnitude of @xmath17 for the planets with estimated radius . for the planets kepler-10**b * * and 55 cnc**e * * ( see section [ k10 ] and [ 55cnc ] ) , we obtain @xmath67 and @xmath68 , respectively , assuming @xmath69 . for sake of comparison with some rocky solar system tidally evolved bodies , we note that @xmath70 in the case of mercury ( smith et al . 2012 ) , @xmath71 for titan ( iess et al . 2010 ) and @xmath72 for io ( schubert et al . 2004 ) . the integration of the exact equations of motion usually requires a long computation time . for this reason , we decided to work with a small set of initial spin - orbit configurations and physical parameters . a statistical analysis taking into account many different initial conditions could give us a more accurate global picture of the past spin - orbit evolution of the planets . as pointed out by correia & laskar ( 2009 ) , such statistical study is more adequate in the context of the averaged equations of the spin - orbit problem ( see section [ mean ] ) . ( solid black curve ) . the dotted line corresponds to the case @xmath73 , where the trappings 2/1 and 3/2 occur . after escaping from the 3/2 resonant motion , the rotation follows the sts towards a synchronization between the rotation and orbital motions . ] kepler-10**b * * is a very close - in planet with a period of @xmath74 0.84 days , orbiting an old ( @xmath75 12 gyr ) sun - like star ( batalha et al . the physical parameters and the current orbital elements of the planet are summarized in table 1 . the estimated mass and radius values ( @xmath76 , @xmath77 ) suggest a rocky composition of the planet , which can be considered as a super - earth ( valencia 2011 ) . the transit technique measurements do not allow a reliable estimation of the current orbital eccentricity , which was arbitrarily fixed at zero in this work . this choice is supported by our numerical simulations of the whole system , which includes also the exterior companion , kepler-10 * c*. the results have shown that the long - term evolution , due to the coupled effect of tides and mutual interaction , could result in the complete circularization of both planetary orbits ( rodrguez et al . , in preparation ) . we adopt the initial semi - major axis at @xmath78 au , which provides @xmath79 through equation ( [ am ] ) . in addition , we set @xmath80 and @xmath81 . for the tidal parameter , we assume @xmath82 ( a value usually adopted for rocky planets ) and , solving equation ( [ k2dt ] ) , we obtain @xmath83 1.1 min .    in section [ sett ] we have estimated that , for kepler-10**b * * on its current orbit , @xmath67 . it should be emphasized that @xmath17 would be smaller if the frozen - in shape of the planet was attained at larger distances than @xmath84 . hence , in order to investigate different possible scenarios , we consider three different values of @xmath17 , namely , @xmath85 , @xmath86 and @xmath87 .     for kepler-10**b * * depending on the value of @xmath17 . different evolutions are observed , including captures in 2/1 , 3/2 and 1/1 spin - orbit resonances ( sors ) . the solid black curve corresponds to the evolution for @xmath88 ( i.e. , when only the tidal torque exist ) and coincides with the stationary tidal solution ( sts , see text for discussion ) of the rotation . the captures in the synchronous motion are detailed on the small panel . ] figure [ stat - k10 ] shows the variation of @xmath89 as a function of the eccentricity of kepler-10**b**. the curve with segments line corresponds to the stationary tidal solution ( hereafter sts ) of the rotation ( @xmath90 ) , obtained when only the tidal torque is considered . the sts is a function which depends only on the planet eccentricity , and is given by    @xmath91    ( see equations ( 42 ) or ( 43 ) of hut ( 1981 ) ; see also section [ mean ] for more details ) . the stationary rotation evolves with the tidally damped eccentricity and the system reaches synchronization ( @xmath92 ) only for circular orbits . note that for small values of @xmath54 , there follows that @xmath93 ; for this reason , the sts is also called pseudo - synchronous solution ( but this is valid only for @xmath94 ) .    the solid black curve in figure [ stat - k10 ] corresponds to the result of the numerical simulation performed with @xmath88 ( only tidal forces are considered ) . in this case , the stationary rotation is attained almost instantaneously ( when compared to the time of circularization ) . note that both curves are in good agreement , at least for @xmath95 ; hence , in this eccentricity interval , we can refer to the solution obtained with @xmath88 as sts . we also note that the curve for @xmath88 reaches the sts for @xmath89 close to two ; however , this is incidental and unrelated with the 2/1 trapping .    figure . [ rot - k10 ] shows the time variation of @xmath89 . for the three considered values of @xmath96 , the rotation rapidly encounters the value @xmath97 and thus the planet rotation is captured in the 2/1 spin - orbit resonance ( hereafter sor ) . on one hand , for @xmath98 and @xmath99 , the rotation escapes from the 2/1 trapping approximately at 80 kyr , when the rotation drops to the value @xmath100 , implying capture into 1/1 sor and synchronous motion . however , it is worth noting that , for @xmath98 , the synchronism is attained in a form of jump , whereas for @xmath99 the rotation evolves without capture for a very short time before the 1/1 sor ( see small panel in fig . [ rot - k10 ] ) . we have also plotted the result for @xmath88 , showing that synchronization is attained around 400 kyr . it is interesting to note that the introduction of the equatorial asymmetry of the planet reduces the time to reach synchronization ; for instance , for @xmath98 , it is reduced by a factor of 5 .    on the other hand , when @xmath73 , the 2/1 sor is destabilized at around 48 kyr giving rise to the trapping in the 3/2 sor ( @xmath101 ) . thus , the permanence of the rotation in the 2/1 sor is temporary and its duration clearly depends on the @xmath17 value . the rotation evolves in a mercury - like state for some time , when a new instability occurs , after which an evolution free of trapping starts . note that during the period of non - capture , the rotation returns to the stationary value , as can be seen in the dotted curve of fig . [ stat - k10 ] when @xmath102 . finally , the synchronization is reached at around 270 kyr .     and @xmath54 depend on the specific capture in which the rotation is trapped in . as different values of @xmath17 lead to several captures , different scenarios are expected for the orbital evolution , despite the final values of the elements being the same for all @xmath17 ( see next figure ) . ]     shows the final evolution . the orbital decay continues until circularization is achieved . the current semi - major axis of kepler-10**b * * coincides with the one obtained in the simulation . ] figures [ ae - zoom - k10][ae - k10 ] show the time evolution of semi - major axis and eccentricity of kepler-10**b**. we note that , when the rotation evolves under capture in resonance ( that is , for @xmath96 ) , the rates of variations of the elements are larger than the one corresponding to the pure tidal case ( @xmath88 ) , with the difference increasing as the eccentricity is damped ( when the comparison is done for a given sor ) . moreover , for the same value of @xmath17 , the rate of orbital decay and circularization seems to depend on the specific sor in which the planet rotation is trapped in . for example , in the case of @xmath73 , the evolutions of @xmath53 and @xmath54 show two different regimes as the rotation evolves first in the 2/1 and then in 3/2 trappings ( see section [ mean ] ) . we note that the time necessary to get an orbit with @xmath103 is much greater for @xmath104 . in fact , when the rotation evolves without capture , the time - scale for circularization increases , according to figures [ ae - zoom - k10 ] and [ ae - k10 ] . it is worth noting that , in the case of @xmath96 , the synchronous motion can be attained even for @xmath105 . for instance , in the case of @xmath99 , the system is captured in 1/1 sor at roughly 120 kyr , when @xmath106 . this fact confirms the classical result that , for eccentric orbits , the synchronous motion is only possible if there exists an additional torque which counteracts the tidal torque ( see goldreich & peale 1966 ; ferraz - mello et al . 2008 and references therein ) . we also note that after circularization , the 1/1 sor is the only possible resonant motion . figure [ am - fig ] shows the evolution of the system in the plane of orbital elements @xmath53 and @xmath54 . the solid line shows the result of numerical simulation in the case of @xmath98 . the dashed line corresponds to equation ( [ am ] ) , which is independent on the value of @xmath17 . note as both curves are in very good agreement , in such a way confirming our assumption about the conservation of the orbital angular momentum ( section [ sett ] ) .    . the solution of equation ( [ am ] ) is shown by the dashed line and is compared to the result of the numerical simulation for the case @xmath98 , giving a good agreement ( the curves for the other two @xmath17 values also agree with the solution of equation ( [ am ] ) , but are not shown here for sake of simplicity ) . the starting position of the planet is indicated by the star symbol . ] gj 3634 b is a super - earth planet with mass @xmath107 orbiting a m2.5 dwarf each 2.65 days with an eccentricity of @xmath108 ( bonfils et al . since there is no transit detection for this planet , its radius is unknown . however , valencia , sasselov & oconnell ( 2007 ) have shown that there is a maximum radius that a rocky terrestrial planet may achieve for a given mass . a larger radius would imply a large amount of water and even the possibility of a massive volatile envelope . indeed , for a planet with mass @xmath109 the maximum terrestrial radius would be 11,600 km ( @xmath110 ) , which result in a mean density of @xmath111 . hence , we assume for gj 3634 b that @xmath112 to obtain @xmath113 . the adopted initial orbital elements are @xmath114 au and @xmath115 , while @xmath116 , @xmath81 and @xmath82 ( @xmath117 min ) . at the current value of the semi - major axis the estimated value of @xmath17 is @xmath118 ( section [ sett ] ) ; in this work , we perform the analysis using two values of @xmath17 : @xmath86 and @xmath87 .     as a function of @xmath54 ( left panel ) and the time evolution of @xmath53 ( right panel ) for the planet gj 3634**b**. depending on @xmath17 , the planet rotation may be currently evolving in either 2/1 or 3/2 sor . note that the time - scale of orbital decay is drastically reduced provided the planet can support a permanent deformation resulting in a nonzero @xmath17 value . the vertical and horizontal dotted lines indicates the current values of @xmath54 and @xmath53 , respectively , title="fig : " ]   as a function of @xmath54 ( left panel ) and the time evolution of @xmath53 ( right panel ) for the planet gj 3634**b**. depending on @xmath17 , the planet rotation may be currently evolving in either 2/1 or 3/2 sor . note that the time - scale of orbital decay is drastically reduced provided the planet can support a permanent deformation resulting in a nonzero @xmath17 value . the vertical and horizontal dotted lines indicates the current values of @xmath54 and @xmath53 , respectively , title="fig : " ]    figure [ rot - exce - gj3634 ] shows the variation of @xmath89 as a function of @xmath54 ( left panel ) and the time variation of @xmath53 ( right panel ) . for @xmath99 , the rotation is rapidly trapped in the 2/1 sor , remaining in such resonant configuration during the most time of the tidal eccentricity damping . only at low eccentricities , close to 0.032 , the system leaves the 2/1 sor dropping directly into the synchronous motion .    at variance , for @xmath73 , the rotation quickly reaches the stationary ( non - resonant ) state given by the sts curve ( @xmath88 ) in figure [ rot - exce - gj3634 ] ( left ) . the system evolves along the sts curve up to close approach to the 3/2 sor . then , the system is captured and remains in this resonance until dropping in the synchronized state , at @xmath119 .    in both cases , the final configuration of the planet s orbit is characterized by circularization and synchronization , as in the case of kepler-10**b**. however , the current eccentricity of the gj 3634**b * * is different from zero ( its location at 0.08 is shown by the dotted vertical line in the left panel of figure [ rot - exce - gj3634 ] ) . note that the current eccentricity is reached when the rotation still evolves in the 2/1 or 3/2 sor ( depending on @xmath17 ) . this fact suggests that gj 3634**b * * could be presently involved in one of such resonant trappings , in particular , it could be found in a mercury - like rotational state ( @xmath120 ) . however , it should be kept in mind , that this conclusion depends on the choice of specific values of @xmath17 and @xmath121 . indeed , the results of an additional simulation performed with @xmath122 , have shown that the trapping of the system in the 3/2 sor is destabilized at @xmath123 . since @xmath124 , the planet rotation would be currently evolving without capture ( along the sts ) for such a small @xmath17 value . the semi - major axis evolution ( figure [ rot - exce - gj3634 ] , right panel ) shows that the tidal decay is faster when the rotation of the planet is trapped in the 2/1 and 3/2 spin - orbit resonant motions . the fact that the current eccentricity of the planet orbit is different from zero suggests that , in the absence of other perturbations , the orbital decay when @xmath92 or @xmath125 ( see equation ( [ a-11-tot2 ] ) ) . ] of the planet is presently in progress and , consequently , the final semi - major axis will be smaller than the current one . 55 cnc**e * * is a hot super - earth planet , with the mass @xmath126 and the radius @xmath127 , orbiting a sun - like star each 0.74 days with @xmath128 ( winn et al . 2011 ; gillon et al . the averaged density of the planet is comparable to that of earth , suggesting a rock - iron composition supplemented by a significant mass of water , gas , or other light elements ( @xmath129 ) indicate that the possibility of 55 cnc**e * * to be a mini - neptune instead of a super - earth should not be ruled out . ] ( winn et al . 2011 ) .    for numerical simulations , we adopt @xmath130 au , @xmath131 , @xmath80 and @xmath81 . the tidal history of the system indicates that the magnitude of @xmath121 must be higher than the commonly adopted values for rocky planets ( ferraz - mello 2012 ) . in this work , we use @xmath132 ; this value allows us to obtain @xmath133 s. we perform numerical simulations for two different values of @xmath17 : @xmath85 and @xmath86 . for the purpose of further comparison with the results reported in callegari & rodrguez ( 2012 ) , we also consider the values @xmath134 and @xmath135 . figure [ rot - exce-55cnc ] shows the evolution of rotation as a function of the planet s eccentricity . the captures inside the 5/2 , 2/1 , 3/2 , 4/3 and 1/1 resonances can be observed . for @xmath134 , the planet rotation evolves under capture in the 3/2 sor almost up to the orbital circularization , with @xmath89 roughly varying between 1.46 and 1.54 .     as a function of @xmath54 ( in logarithmic scale ) , for the planet 55 cnc**e**. we have used @xmath98 in addition to the value adopted in callegari & rodrguez ( 2012 ) @xmath134 . also values ten times smaller were taken into account for sake of completeness . the vertical dotted line indicates the current value of @xmath54 , namely , 0.057 . the trappings in the rsos 5/2 , 2/1 , 3/2 , 4/3 and 1/1 appear during the evolutions . note the large amplitude oscillation of @xmath89 inside the 3/2 and 4/3 cases with @xmath17 of the order of @xmath85 . depending on @xmath17 , 55 cnc**e * * can be an example of a planet currently evolving in the 5/2 , 2/1 or 3/2 sor . ]     together with the sts curve . as for previously analysed planets , the orbital decay is faster when the planet rotation is trapped in a resonant motion . the current value of @xmath53 is obtained for all values of @xmath17 . ]    in a recent work , callegari & rodrguez ( 2012 ) have shown that , for @xmath128 and @xmath136 ( or @xmath134 ) , the phase space of the rotation ( @xmath137 ) of 55 cnc**e * * , between the main spin - orbit resonances , is dominated by a large region of chaotic motion . this region increases as larger values of the eccentricity are considered . capture in spin - orbit resonance can work as a mechanism to prevent a chaotic rotation , provided that the dissipation is large enough in order to ensure the corresponding trapping . according to our simulations ( see figure [ rot - exce-55cnc ] ) , destabilizations of different resonances occur as the eccentricity is tidally damped ( see also next section ) , giving rise to new temporary trappings with a subsequent evolution to synchronism . in this way , a chaotic behavior of the rotation can be avoided by a continuous evolution under resonant motion . it is important to note that , at least for the cases shown in figure [ rot - exce-55cnc ] , there is no interval of non - resonant motion between two successive trappings , as in the case of kepler-10**b * * ( see figure [ rot - k10 ] for @xmath99 and @xmath87 ) .    an additional simulation done with @xmath138 ( not shown in figures ) , has resulted in a 2/1 trapping until @xmath139 , indicating that the planet rotation may be currently evolving in resonant motion , for such a small value of @xmath17 . figure [ elem-55cnc ] shows the time variation of semi - major axis . it is worth noting that the orbital decay results in a more rapid evolution towards @xmath84 when the planet rotation is trapped in a sor . also , the rate of variation depends on the specific resonance . from figure [ elem-55cnc ] we see that the evolution when @xmath140 , ( capture in 2/1 ) is almost two times faster than the one for @xmath134 ( 3/2 capture ) . we will return to this discussion in the next section . in order to better understand the results previously obtained through numerical simulations of the exact equations of motion , in this section , we develop an analytical approach in the context of the averaged equations of the spin - orbit problem . the equation of motion for the spin in equation ( [ mov1 ] ) can not be solved analytically . however , because we are interested in the evolution inside spin - orbit resonances , we will consider planetary spin rates having values such that @xmath141 , where @xmath142 is a half - integer . our investigation is restricted to the cases @xmath143 , @xmath144 and @xmath145 but can be easily extended to any other sor . we introduce the `` resonant '' angle @xmath146 , where @xmath147 is the mean anomaly of the planet orbit and @xmath15 was defined in section [ model ] . note that the angle @xmath148 gives the orientation of the long axis whenever the planet passes through the pericenter of its orbit ( i.e. , @xmath149 when @xmath150 , see figure [ elips ] ) . it is easy to show that , inside the spin - orbit resonance given by the parameter @xmath142 , the resonant angle is a slowly varying variable when compared to the orbital period of the planet . indeed , we can write that @xmath151 , where @xmath152 and , consequently , @xmath153 . now , replacing @xmath15 by @xmath154 in equation ( [ t22 ] ) and averaging over the orbital period , we obtain the classical result for the averaged torque on the planet produced by the permanent equatorial deformation ( see goldreich 1966 ; goldreich & peale 1966 ) : @xmath155 where @xmath156 are power series of @xmath54 with the leading term of order @xmath157 ( see appendix a for a literal expansion up to seventh order ) . the averaged tidal torque on the planet is given by @xmath158 where @xmath159 and @xmath160 is a suitable function , whose expansion in powers of @xmath54 is shown in appendix a ( hut 1981 ) .    at resonance , the equilibrium solution of the rotation is given by the condition    @xmath161    and , resolving for @xmath162 , we obtain @xmath163 if we know the maximum value that @xmath148 can reach without disturbing the trapping in the @xmath142sor , we may obtain the minimum value of @xmath17 necessary for capture in the same resonant motion .      replacing @xmath143 in equation ( [ j22-gral ] ) and expanding the ratio @xmath164 in powers of @xmath54 , up to order 6 , we obtain    @xmath165    we note that , for @xmath166 , the above expression is the same obtained in ferraz - mello et al . ( 2008 ) , for synchronous asymmetric planets . should be done , where @xmath167 is the phase lag of the tidal wave whose frequency is @xmath168 . ] the averaged variation of the orbital elements can be obtained using the euler - gauss equations ( see beutler 2005 , p. 230 ) , which are written in terms of the radial and transverse components of the perturbing force . the result for the semi - major axis is @xmath169 replacing @xmath162 by the corresponding expression given in equation ( [ j22 - 11 ] ) , we obtain @xmath170    in the context of the linear tidal model , the averaged variation of the semi - major axis due to tides is given by equation ( 9 ) of hut ( 1981 ) , valid for any values of @xmath142 and @xmath54 . the expansion up to @xmath171 reads    @xmath172}+{\cal o}(e^8).\ ] ]    replacing @xmath143 in the above equation we obtain    @xmath173    the total averaged variation of the semi - major axis is then given as    @xmath174    and thus @xmath175 at the same order of approximation , the above equation coincides with the result of previous works ( see ferraz - mello et al . ( 2008 ) and references therein ) . when only the tidal torque on the planet is taken into account , the non - resonant sts is given by equation ( [ sync ] ) . in this case , the synchronism is possible only for circular orbits . if instead of @xmath143 , we replace @xmath176 in equation ( [ a - tide - gral ] ) , we obtain    @xmath177    we define @xmath178 as the ratio between averaged semi - major axis variation in the @xmath142sor ( @xmath179 ) and in the tidal stationary solution of the rotation ( @xmath180 ) , that is    @xmath181    comparing equations ( [ a-11-tot2 ] ) and ( [ a-11-tide - stat ] ) , we note that @xmath182 , indicating that the rate of semi - major axis decay is stronger when the rotation is locked in the 1/1 sor , compared to the case when it evolves following the sts . this result is in agreement with wisdom ( 2008 ) , who has shown that the rate of energy dissipation in a body with synchronous rotation is larger than in a asymptotic non - synchronous rotation . we repeat the procedure described in the previous section to obtain the averaged variation of the planet eccentricity . the resulting expressions are :    @xmath183    @xmath184    @xmath185    @xmath186    similarly to the semi - major axis variation case , @xmath187 , where @xmath188 is the analogous for the eccentricity of equation ( [ ra ] ) . we analyse in this section the spin - orbit behavior when the planet evolves in the 3/2 and 2/1 sors , employing the same method described in the previous section has no term independent of @xmath54 and thus equation ( [ j22-gral ] ) can not be expanded around @xmath189 . ] . we have    @xmath190    and    @xmath191     as a function of orbital eccentricity , for @xmath145 , 3/2 and 1 . in the cases of @xmath192 the ratio of rates of variation of @xmath53 in the sor and in the sts increases for small @xmath54 , showing that the orbital decay is enhanced ( when compared with the evolution in the sts ) whenever the planet rotation evolves under resonance with the orbital motion . ] figure [ razao ] shows a plot of @xmath178 for the three analysed resonances as a function of the eccentricity . on one hand , we note that @xmath193 and @xmath194 increase for small eccentricity , indicating that the major effect on the orbital evolution due to the resonance captures occurs as @xmath54 decreases , in agreement with the results of the numerical simulations ( see figures [ ae - zoom - k10 ] , [ rot - exce - gj3634 ] and [ elem-55cnc ] ) . on the other hand , for large @xmath54 , @xmath195 is small ( mainly in the case @xmath144 ) as also shown in figures [ ae - zoom - k10 ] and [ rot - exce - gj3634 ] , where there is a small difference between the evolutions with and without @xmath17 ( i.e , @xmath196 ) .    when the synchronous motion is analysed , we note in figure [ razao ] that @xmath197 is an increasing function of @xmath54 , indicating that the orbital decay when the rotation is trapped in the 1/1 sor is large for large @xmath54 . however , as simulations show , synchronism is attained for very low values of @xmath54 , making the orbital decay almost negligible . the results above are also valid for @xmath198 . it is interesting to note that @xmath178 and @xmath198 are only functions of the eccentricity and do not depend on any physical parameter . it is worthwhile mentioning that @xmath199 averages to zero unless the condition @xmath141 is assumed ( whit @xmath142 integer ) , that is , the averaged orbital evolution due to the triaxial torque is expected only for resonant cases . as it has been noted , equation ( [ j22-gral ] ) can be used to obtain the minimum value of @xmath17 necessary for the stability of the @xmath142sor . calling @xmath200 the one value given by the right hand of equation ( [ j22-gral ] ) , the stability condition for a resonant spin state requires that @xmath201 . the above condition is frequently found in literature as strength criterion of the resonance ( see murray & dermott 1999 , p. 205 ) . figure [ c22 ] shows the functional dependence of the absolute value of @xmath200 with the eccentricity , obtained for kepler-10**b**. the dashed horizontal lines represent three values of @xmath17 used in the numerical simulations shown in section [ k10 ] , namely , @xmath202 and @xmath87 . the case @xmath145 , shown by the blue curve , indicates that @xmath203 for @xmath204 , @xmath205 for @xmath206 and @xmath207 for @xmath208 . a good agreement is obtained doing the comparison with the results of the numerical simulations for kepler-10**b * * , shown in figure [ rot - exce - k10 ] . note as the 2/1 sor destabilizes when @xmath209 , @xmath206 and @xmath208 . the same can be observed for the 3/2 sor , in which the rotation escapes from the trapping when @xmath123 if @xmath73 , in concordance we use @xmath210 , taking also into account that truncated series in @xmath54 were used for the same figure . ] with the result of figure [ c22 ] ( orange curve ) .    since for roughly @xmath211 the condition @xmath201 is satisfied more easily for @xmath143 as compared to @xmath192 , there is no problem in understanding the occurrence of synchronous motion for very small eccentricity , as can be also verified in figure [ rot - exce - k10 ] . it is worth noting that instead to compute a critical @xmath17 value , the condition for stationary solutions @xmath212 also enables us to obtain a critical value of @xmath54 for a given fixed @xmath17 . the resonant motion would destabilize whenever the eccentricity becomes smaller than the critical value for each sor . this approach is more convenient in those cases for which the observations allow an estimation of @xmath17 ( see correia & laskar ( 2004 ) for the study of the rotation of mercury ) .     given by equation ( [ j22-gral ] ) and calculated with parameters of the planet kepler -10**b**. the calculation of @xmath200 can be used as a guide to investigate the stability of a given sor as a function of eccentricity ( see text for discussion ) . note that in the case of @xmath143 we have @xmath213 for @xmath214 ( see equation ( [ j22 - 11 ] ) ) , indicating that the synchronism is the only possible resonant motion for circular orbits ( see goldreich & peale 1966 ; celletti & chierchia 2008 ; callegari & rodrguez 2012 ) . ]     as a function of the orbital eccentricity for kepler-10**b * * , resulting from the analysis of numerical simulations of the exact equations of motion ( see section [ k10-res ] ) . if a comparison with the results of figure [ c22 ] is done , we obtain a good agreement for the values of @xmath54 for which each sor destabilizes ( see text for details ) . ] the numerical exploration of the planet rotation which is subject to the conservative triaxial gravitational torque of the star shows different behaviors , including the oscillation around the equilibrium points of spin - orbit resonances ( see callegari & rodrguez 2012 and references therein ) . when dissipative effects produced by tidal torques are included , the rotation of close - in bodies can be captured in a resonant motion ( goldreich 1966 ; goldreich & peale 1966 ; celletti & chierchia 2008 ) . the capture in a specific resonance depends on the eccentricity @xmath54 and on parameters related to tides and equatorial ellipticity ( @xmath56 and @xmath17 , respectively ) . hence , since @xmath54 is tidally damped during the inward migration , each specific capture should become unstable and the rotation can achieve another resonant state , which , at the same time , should result in a temporary trapping . when the orbit is circularized , the final evolution results in the synchronization between the orbital and rotational periods ( i.e , the 1/1 spin - orbit resonance ) .    in this paper , we analysed the coupled spin - orbit evolution of kepler-10**b * * , gj 3634**b * * and 55 cnc**e * * short - period super - earth planets . the numerical simulation of the exact equations of motion have shown that , depending on @xmath17 , the rotation displays a number of resonance trappings for which @xmath141 , including the cases @xmath215 , 2 , 3/2 , 4/3 and 1 . moreover , the results indicate that gj 3634**b * * can be an example of an exoplanet with a rotation state similar to mercury s rotation ( @xmath144 ) , provided the body figure can support a permanent equatorial asymmetry as small as @xmath73 . 55 cnc**e * * can currently evolve under 5/2 or 3/2 trappings if @xmath216 and under 2/1 if @xmath140 or even if @xmath17 is as small as @xmath217 . the rotation of 55 cnc**e * * may have a chaotic behavior for @xmath216 ( callegari & rodrguez 2012 ) . however , we have shown that a scenario is possible in which the rotation evolves to synchronism through successive trappings in sor . indeed , the eccentricity damping due to tides provides the destabilization mechanism required to escape from one given capture and subsequently evolve to lower order resonances ( see correia & laskar 2004 , 2009 ) . we have shown that different spin - states corresponding to various resonance captures can occur in extrasolar systems . although the 1/1 spin - orbit coupling is the final state of evolution , it is not reached in every case as present configuration . for current non - circular orbits , trappings in the 3/2 sor or other resonance configurations can be sustained over long times . the latter would depend on the orbital eccentricity , the planet s @xmath17 value and its dissipation factor . given that these parameters are not well - determined for terrestrial extrasolar planets , we note that the 1/1 coupling may not be the current rotation state of some planets . this would have consequences on the stellar insulation and therefore on the evolution of the atmosphere and surface of the planet . when only the tidal torque is taken into account ( @xmath88 ) , the rotation evolves to the stationary tidal solution ( sts ) . we have shown that under assumption of @xmath96 , the evolution along spin - orbit resonance ( sor ) leads to higher rates of orbital decay and circularization than the corresponding ones when the rotation evolves in the sts . moreover , for a given @xmath17 value , the variation rates of the elements depend on the specific trapping in which the rotation is trapped in , and increases with @xmath142 ( see figure [ ae - zoom - k10 ] ) . we have also calculated the averaged rates of variations of semi - major axis and eccentricity including the gravitational and tidal torques . moreover , through the analytical calculation of the critical value of @xmath17 it is possible to predict the values of the eccentricity for which a given sor would become unstable . the comparison with the direct numerical simulations have shown a good concordance for the analysed planets . the importance of the @xmath17 value on the rotational history of the planet is illustrated in the example of kepler-10**b**. it is worth noting that the planet rotation of kepler-10**b * * can evolve under the influence of the 2/1 sor for longer time - scales ( i.e. , small @xmath54 ) provided the planet figure can support a permanent equatorial deformation characterized by @xmath216 . for smaller @xmath17 , the 2/1 resonant state would become unstable for relatively large eccentricity ( figures [ c22 ] and [ rot - exce - k10 ] ) . we stress that the results of the numerical simulations provide different possible spin - orbit evolutions depending on the considered initial values of elements and rotation . the @xmath56 value and the specific tidal model would also play an important role in determining the rotation history of the planet , since the temporary trapping in sor depends on the rate of eccentricity damping , which is @xmath56-dependent . in fact , in a recent paper , makarov ( 2012 ) has applied the model by efroimsky ( 2012 ) and successfully explained the current resonant state of mercury within the considered model . our results concerning the association of the critical @xmath17 value with the capture probability are not new in the context of averaged equations of the spin - orbit problem ( see for instance correia & laskar 2004 , 2009 ) . however , in this work , we are able to compare the results of the averaged equations with those resulting from the numerical simulations of the exact equations of motion .    a model of planet interior including the effect of viscous friction at the core - mantle boundary would provide quantitatively different results ( see peale 2005 , and correia & laskar 2009 , for the case of mercury ) . according to previous works the chances of resonant capture increase if the latter effect is taken into account ( goldreich & peale 1967 ; peale & boss 1977 ) . however , in order to highlight the contribution of @xmath17 to the orbital evolution , in this paper we adopt a simple model neglecting a possible layered structure of the planet .    in the present study we have assumed constant @xmath17 values . however , depending on the adjustment to hydrostatic equilibrium @xmath17 may be itself time - dependent . the time - scale for adjustment to the completely relaxed state depends on material parameters , e.g. , rigidity and viscosity which are highly temperature dependent . furthermore , it will depend on the spin - orbit coupling the planet is trapped in . for the application to terrestrial planets given here the constant , i.e. frozen , @xmath17 configuration is justified because adjustment of solid rocky planets to equilibrium usually occurs on geological time - scales . therefore , the states described here do not necessarily represent hydrostatic configurations .    in this work we have considered a model of one interacting planet , despite the fact that kepler-10**b * * and 55 cnc**e * * belongs to multiple - planet systems with exterior ( more massive ) companions . previous works have shown that the mutual secular interaction with the companion may affect the eccentricity evolution of the tidally evolved planet , increasing the circularization time - scale ( mardling 2007 ; rodrguez et al . 2011a ; rodrguez et al . moreover , the planetary perturbations can produce large excursions of the innermost planet eccentricity , providing a mechanism for increasing or decreasing the capture probability of the spin - orbit resonances ( correia & laskar 2004 , 2009 ) . however , a numerical simulation of 55 cnc system including the exterior companion 55 cnc**b * * has shown little differences with the single planet case . indeed , the evolution of the inner orbit eccentricity is almost not affected by the mutual interaction between the planets and thus the rotation evolution of 55 cnc**e * * can be safely analysed within the one planet model . an investigation about the rotation of asymmetric close - in super - earths including an exterior companion will be addressed in a future investigation . the authors acknowledge the support of this project by fapesp ( 2009/16900 - 5 ( a.r ) ; 2006/58000 - 2 ( n.c ) ) and cnpq . the authors also gratefully acknowledge the support of the computation centre of the university of so paulo ( lcca - usp ) . we are grateful to the reviewer for his / her valuable comments . we show here the explicit expressions for the functions @xmath156 and @xmath218 . the exact form of @xmath156 is given by
Abstract:

Answer: we investigate the spin behavior of close - in rocky planets and the implications for their orbital evolution . considering that the planet rotation evolves under simultaneous actions of the torque due to the equatorial deformation and the tidal torque , both raised by the central star , we analyze the possibility of temporary captures in spin - orbit resonances .    the results of the numerical simulations of the exact equations of motions indicate that , whenever the planet rotation is trapped in a resonant motion , the orbital decay and the eccentricity damping are faster than the ones in which the rotation follows the so - called pseudo - synchronization . analytical results obtained through the averaged equations of the spin - orbit problem show a good agreement with the numerical simulations . we apply the analysis to the cases of the recently discovered hot super - earths kepler-10**b * * , gj 3634**b * * and 55 cnc**e**. the simulated dynamical history of these systems indicates the possibility of capture in several spin - orbit resonances ; particularly , gj 3634**b * * and 55 cnc**e * * can currently evolve under a non - synchronous resonant motion for suitable values of the parameters . moreover , 55 cnc**e * * may avoid a chaotic rotation behavior by evolving towards synchronization through successive temporary resonant trappings . celestial mechanics  planets and satellites : general .