Instruction: Write an abstract for the following acadamic paper.

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

Paper: it is now well accepted that the giant planets of the solar system did not have the current orbits when the gas was removed from the protoplanetary disk . the first work showing that the planets had to migrate while interacting with the remaining planetesimal disk was @xcite . this work demonstrated that , in the planetesimal scattering process , on average saturn , uranus , and neptune gain angular momentum and thus migrate outwards , while jupiter loses angular momentum and moves inwards . since then , the process of planetesimal - driven migration was considered to be a fundamental one in shaping the current structure of the solar system . @xcite showed that planetesimal - driven migration can explain the origin of the orbit of pluto and of the resonant populations of kuiper belt objects . however , it appeared soon that planetesimal - driven migration alone is a process too smooth to explain several aspects of the kuiper belt s orbital distribution , e.g. the orbital excitation of the non - resonant objects . @xcite was the first paper showing that a dynamical instability of the giant planets could lead to a global orbital excitation of the kuiper belt , while delivering the giant planets on final orbits roughly similar to the current ones . building on these results , the so - called _ nice model _ @xcite combined the processes of planetesimal - driven migration and dynamical instability into a new global scenario . it suggested that the planets jupiter , saturn , uranus , and neptune formed in a compact configuration with jupiter around 5.45 au , saturn slightly interior to the 2:1 mean motion resonance ( mmr ) with jupiter ( @xmath3 @xmath4 8.65 au ) , and the initial semimajor axes of the ice giants ( uranus and neptune ) in the ranges 1113 au and 13.517 au . the planets were surrounded by a disk of leftover planetesimals , accounting for a total of about 35 @xmath5 , from @xmath016 to 30 au ( where @xmath6 is the mass of the earth ) . in this model the planets interacted with the planetesimal disk and , by planetesimal - driven migration , their mutual separations slowly increased . when jupiter and saturn crossed their mutual 2:1 mmr the full system of the giant planets became unstable and spread further through mutual close encounters . in particular , uranus and neptune were scattered onto orbits with relatively large eccentricities and semimajor axes , which crossed the original trans - neptunian disk . then , following a phase of dynamical friction caused by their interaction with the disk , the orbits of uranus and neptune damped in eccentricity and inclination . a residual planetesimal - driven migration placed the planets on orbits fairly similar to the current ones @xcite . the nice model had great success in explaining many features of our solar system , as : the capture of jupiter s trojans @xcite , the capture of the irregular satellites of the giant planets @xcite , the absence of regular moons beyond oberon s orbit at uranus @xcite , to cite a few . moreover , it was shown that , provided the inner edge of the trans - neptunian disk was at an appropriate distance from the original orbit of neptune , the giant planet instability could have occurred after hundreds of my of evolution , potentially explaining the origin of the so - called late heavy bombardment ( lhb ) of the terrestrial planets @xcite . however , the original version of the nice model presented a major problem ; the initial configuration of the planets was chosen in an ad hoc manner . @xcite performed hydrodynamical simulations of the evolution of the giant planets when they were still embedded in a disk of gas and found that the planets should have acquired a fully resonant configuration , with jupiter and saturn preferentially locked in their mutual 3:2 mmr with jupiter at @xmath05.4 au , saturn and uranus ( ice i ) in the 3:2 mmr , and uranus ( ice i ) and neptune ( ice ii ) in the 4:3 mmr ( although other resonant configurations are also possible ) . the authors also concluded that such a multi - resonant system is compatible with a subsequent evolution similar to that described in @xcite , although the instability in this case happens when a pair of planets leave the original resonant configuration , rather than when jupiter and saturn cross their 2:1 mmr . @xcite investigated whether , under some conditions , this initial multi - resonant configuration could still lead to a _ late _ giant planet instability , in order to explain the putative lhb as in the original @xcite paper . they found that late instabilities are possible , provided that the inner edge of the disk is several au s beyond the outermost planet and about 1,000 pluto - sized planetesimals are embedded in the disk , so to produce an important dynamical self - stirring of the planetesimals population .    during the dynamical instability of the giant planets , several dynamical paths are possible , even when restricting the evolutions to those leading to final planetary orbits similar to the current ones . the actual dynamical evolution of the giant planets at the instability time can be constrained by investigating the consequences on the other components of the solar system : the terrestrial planets and the small body populations . the terrestrial planets and the asteroid belt give similar constraints , as investigated in @xcite and @xcite . if the orbits of jupiter and saturn migrate away from each other smoothly and slowly , as in planetesimal - driven migration , the @xmath7 secular frequency ( dominant in the precession of jupiter s longitude of perihelion ) decreases slowly in value and becomes very similar for some time to the @xmath8 secular frequency , dominant in the precession of mercury s perihelion . as a consequence , the orbit of mercury is destabilized @xcite . similarly , the @xmath9 secular frequency , dominant in the precession of saturn s perihelion , as it slows down becomes equal to the frequencies characterizing the precession of the asteroids at low inclination in the inner asteroid belt , so that this region becomes severely depleted @xcite . both papers concluded that the ratio of the orbital periods of jupiter and saturn @xmath10 had to evolve almost discontinuously from @xmath11 to @xmath12 ( the current value is @xmath13 ) . this happens if jupiter and saturn encounter another planet ( uranus or neptune or a third planet of comparable mass ) . the sequence of encounters causes a divergent jump in the semimajor axes of the orbits of jupiter and saturn , and therefore a jump in the @xmath14 ratio . this scenario is commonly known as the _ jumping - jupiter _ model . notice that if the giant planet instability occurred early , the terrestrial planets were possibly not yet fully formed and therefore the constraint on the evolution of @xmath7 relative to @xmath8 does not apply @xcite . however , the constraint set by the asteroid belt is still valid , supporting the need for a _ jumping - jupiter _ evolution regardless of the time of the instability @xcite .    given that in a _ jumping - jupiter _ evolution the planet that encounters jupiter is often ejected onto a hyperbolic orbit , @xcite and @xcite suggested that the solar system had originally a third ice giant planet with a mass comparable to uranus or neptune . this scenario increases the probability that a _ jumping - jupiter _ evolution ends with four planets near their current orbits , once the putative fifth planet is ejected after its encounter with jupiter . then , @xcite considered many different initial multi - resonant configurations for jupiter , saturn , uranus , neptune , and the rogue fifth planet , and determined for each of them the probability that the @xmath14 ratio had a jump across the 2.12.3 range and the final planets reached orbits compatible with those of our current solar system . the successful simulations of @xcite have then been used to investigate again the capture of the trojans of jupiter @xcite , the capture of irregular satellites at jupiter @xcite , the excitation of the inclination of the moon iapetus at saturn @xcite , the evolution of the asteroid main belt and of its collisional families @xcite . they have also been shown to be compatible with the survival of the galileans satellites @xcite , despite of the encounters of jupiter with the ejected planet .    however , when investigating the dynamical sculpting of the kuiper belt , @xcite and @xcite realized that , of all possible evolutions of neptune consistent with the five - planet _ jumping - jupiter _ model , only some are consistent with the current orbital structure of the trans - neptunian population . in order to get the correct inclinations of the hot population of the kuiper belt objects ( kbo ) , @xcite found that neptune should have migrated more that 5 au from its original resonance by planetesimal - driven migration before that the planets became unstable . this migration should have occurred on a timescale @xmath15 @xmath16 10 my , with the planet on a quasi - circular planar orbit ( @xmath17 0.1 and @xmath18 ) to avoid excessive orbital excitation of the cold classical belt , in the @xmath04245 au region . the planetary instability should have happened when neptune was already at 28 au @xcite in order to explain the so - called _ kernel _ of the cold kuiper - belt ( a clump of objects around 44 au ; @xcite ) . in this case the _ kernel _ would have formed from objects transported outwards in the 2:1 mmr with neptune and released from the resonance when neptune s orbit had a jump of about 0.5 au in semimajor axis , due to a close encounter with another planet . finally , in order to reproduce the observed ratio between the resonant and non - resonant kuiper belt populations , @xcite concluded that the planetesimal - driven migration of neptune should have been characterized by several small amplitude jumps , as those caused by a population of 1,0004,000 pluto - size objects in the disk . the brief review reported above , sets the context that motivated the study reported in this paper . it is clear that the work by @xcite needs to be revisited to test which initial resonant configurations are the most compatible with the new constraints on the dynamical evolution of neptune . this is done in the first part of the paper ( section  [ early ] ) . in the second part of the paper ( section  [ late ] ) we investigate whether the initial configurations selected in the first part of the work can lead to a _ late _ planetary instability , while still showing an evolution of neptune consistent with kuiper belt constraints . the possibility of a _ late _ instability is far from trivial . in fact , in order to exhibit a long - range planetesimal - driven migration before the planet instability , neptune needs to be embedded in the planetesimal disk . but in this case , neptune would start migrating away from the other planets from the very beginning of the simulation ( corresponding to the gas - removal time ) . this is likely to lead to an instability after just a few 10s of my . if instead the planetesimal disk is relatively far from neptune , the planet can be extracted from its initial resonance at a _ late _ time as in @xcite , but the instability is likely to happen as soon as the initial resonant configuration is broken , i.e. before that neptune can migrate far away . for this problem we will introduce a new concept : that of dust - driven migration , i.e. , migration due to the gravitational interaction of a planet with the dust produced in a distant planetesimal disk , spiraling inward by poyning - robertson drag . we will study the evolution of the system as a function of the dust properties and the numerical resolution used in modeling the planetesimal disk . finally , our results are summarized and discussed in section [ conclusions ] . @xcite characterized the evolution of the giant planets from several possible initial resonant configurations . they assumed that the solar system could have formed with 4 , 5 , or 6 giant planets . the additional planets had a total mass comparable to that of uranus or neptune . they analyzed the rate of success of each configuration by measuring the fraction of the planetary orbital evolutions that fulfilled a set of constraints , namely : the final planetary system is made of four giant planets ( criterion a ) , the orbits of the 4 remnant planets resemble the present ones ( criterion b ) , the proper mode of jupiter s eccentricity ( @xmath19 ) is greater than 0.022 in the final system , i.e. , at least half of its current value ( criterion c ) , and the ratio @xmath14 evolves from @xmath11 to @xmath12 in less than 1 my ( criterion d ) . criterion c is important because it is found that , although encounters between jupiter and an ice giant can excite @xmath19 , not all of them can do it as much as needed . criterion d is important because it is linked with the evolution of the terrestrial planets and of the asteroid belt @xcite . if jupiter and saturn migrate too slowly through the period ratio 2.12.3 , the @xmath20 and @xmath21 resonances sweep through the terrestrial planet region and the inner asteroid belt , resulting in orbital distributions inconsistent with those observed . then , from the success rates , @xcite concluded that a planetary system initially having 5 planets should be preferred .    here , we revisit the @xcite work , adding a new constraint on the orbital evolution of neptune from the work of @xcite , namely that neptune should have migrated to @xmath028 au before the onset of the planetary instability . the works by @xcite and @xcite also give other constraints on the migration timescale of neptune and its eccentricity and inclination before the instability . the eccentricity and inclination constraints are met fairly easily because in this scenario neptune evolves in a planetesimal - driven migration without planetary close encounters until the instability phase . as for the migration timescale , @xcite claims that before the planetary instability , once out of resonance , neptune can migrate according to @xmath22 10 my for disks with 20 @xmath5 and according to @xmath22 50 my for disks with 15 @xmath5 . considering that those disks have the same extension , this exemplifies how sensitive the migration timescale can be with respect to the density of the disk . therefore , we decide not to consider the migration timescale as an important parameter at this point , because getting a correct timescale of the migration is probably only a matter of tuning the density and profile of the planetesimal disk . thus we add just one new criterion ( criterion e ) : at the time of the instability @xmath23 has to be in the interval 2729 au . while criteria a to d indicate that the initial solar system most probably hosted five giant planets @xcite , criterion e is a very important criterion that can be used to constrain in which resonant chain these five planets likely were at the end of the gas - disk phase .    for simplicity , we consider only the cases of @xcite that involved 5 planets . we also always considered jupiter initially @xmath05.4 au , a distance between the outermost ice giant and the inner edge of the planetesimal disk @xmath24 1 and 2 au and extending to 30 au @xcite , composed by 1,000 planetesimals ( which is the resolution considered by @xcite ) , and two values of surface density ( @xmath25 ) equal to 0.02 and 0.03 @xmath26 . such configuration implies that the mass of the planetesimal disk ranges from @xmath013 @xmath5 ( in the more relaxed configuration  3:2 , 3:2 , 2:1 , 2:1 5:4 au , and the three equal - mass ice - giant planets placed beyond them ] , with @xmath27 2 au and @xmath28 0.02 @xmath26 ) to @xmath065 @xmath5 ( in the more compact configuration  3:2 , 3:2 , 4:3 , 5:4 , with @xmath27 1 au and @xmath28 0.03 @xmath26 ) . to perform our simulations , we used the hybrid integrator in the _ mercury - package _ @xcite , with a time step of half an year . table [ table1 ] gives in column 6 the rate of success for our new criterion e , after 20 runs of each initial configuration from ( * ? ? ? * table 6 in their paper ) . those 20 simulations were composed by 5 simulations for each pair @xmath29 [ 1,2 ] au , and @xmath28 [ 0.02 - 0.03 ] @xmath26 . although we have used two different values of surface density and two values of @xmath30 for the planetesimal disk , we decided to compress the results of all 20 simulations , for each configuration , in table [ table1 ] , column 6 . this is because , concerning our criterion e , only minor differences were observed with these changes . thus , when important , those minor differences will be reported . for completeness and to increase statistics , table [ table1 ] also reports the average of the rate of success for the other criteria taken from ( * ? ? ? * table 6 in their paper ) . because the simulations performed in this work are essentially the same of @xcite , the success rate for criteria a to d are consistent ( within the small number statistics ) with those reported in @xcite thus , we did not report our own results for the criteria a , b , c , and d. rather than that , we decided to check only the success rate relative to the new criterion e , and then , combine our findings with those of @xcite , which have much higher statistics , in order to determine which of the initial planetary configurations are the most compatible with the constraint on the dynamical evolution of neptune , unveiled in @xcite . lccccc 3:2 , 3:2 , 4:3 , 5:4 & 24 & 12 & 0 & 2 & 5 + 3:2 , 3:2 , 3:2 , 3:2 & 30 & 15 & 2 & 5 & 15 + 3:2 , 3:2 , 4:3 , 4:3 & 47 & 23 & 3 & 3 & 5 + * 3:2 , 3:2 , 2:1 , 3:2 * & * 32 * & * 13 * & * 4 * & * 7 * & * 65 * + 3:2 , 3:2 , 2:1 , 2:1 & 38 & 18 & 4 & 4 & 25 + 2:1 , 3:2 , 3:2 , 3:2 & 53 & 40 & 2 & 18 & 0 + 2:1 , 4:3 , 3:2 , 3:2 & 73 & 33 & 0 & 18 & 40 +    our results can be interpreted as follows :    : in these configurations , as it also appears in ( * ? ? ? 16 ) , it is possible that neptune leaves the resonant configuration and migrates ahead the other planets , but not always in the way required in @xcite . indeed , a configuration with initial resonances 2:1 , 3:2 , 3:2 , 3:2 never lead to an instability that fulfilled criterion e in our 20 simulations . the evolutions from the 2:1 , 4:3 , 3:2 , 3:2 configuration , on the other hand , had one of the largest rate of success for criterion e. however , such configurations present other major problem , mostly concerning the final orbital separation between jupiter and saturn , which is typically too big . the point is that , because the ratio @xmath14 has to change less than 0.5 if jupiter and saturn are initially in the 2:1 mmr , the mass of the fifth planet should be @xmath01/2 of the mass of neptune . instead , we assumed for this planet a mass equal to that of uranus - neptune , as in @xcite . a smaller mass planet , however , would not excite jupiter s eccentricity enough , so that the criterion c would not be fulfilled . moreover , hydrodynamical simulations @xcite show that in a massive gas - disk , jupiter and saturn would be locked in their mutual 3:2 mmr , rather than 2:1 . the grand tack model @xcite also implies that jupiter and saturn were in their 3:2 mmr at the disappearance of the disk of gas . so , we do nt consider that configurations with jupiter and saturn in the 2:1 mmr are very realistic . : these two configurations are too compact . only one out of 20 simulations for each configuration fulfilled criterion e. usually , the evolution of the planetary system is as in ( * ? ? ? 7 ) , where the instability takes place as soon as one planet goes out of resonance . as a result , no large - range planetesimal - driven migration occurs before the planetary instability . instead , the instability is violent and the evolution of neptune is largely different from that constrained in @xcite ( figure [ fig1 ] top left , in this work ) . nevertheless , although they are low , none of the success rates reported in table [ table1 ] for the configuration 3:2 , 3:2 , 4:3 , 4:3 are null . thus , even being very unlikely , due to the reasons presented , this configuration can not be ruled out .    : this configuration is too loose . @xcite already pointed out that neptune stops farther than 30 au for whatever reasonable disk s parameters starting for this configuration . here , we found that even with the planetesimal disk of lowest mass ( @xmath31 [ 1324 ] @xmath5 ) , ranging from @xmath32 [ 25(26)30 ] au , neptune goes out of resonance and usually migrates all the way up to 30 au or more , before that the planetary instability happens ( figure [ fig1 ] top right ) . however , sometimes the instability can also occur before neptune reaches 30 au , say at @xmath028 au . therefore , although this configuration may also present some problems with the extension of the planetesimal disk , as it will be discussed in the following sections , it does have a good rate of success in fitting criterion e ( 25% success rate ) . : this configuration allows neptune to migrate away from the other planets before the instability , but usually not enough . the instability commonly happen when neptune is around 20 au ( see * ? ? ? * figs 8 and 9 , and figure [ fig1 ] bottom left in this work ) . only for high density disks neptune can migrate fast enough and reach @xmath028 au before the planetary instability . if the disk s density is low , neptune migrates slowly and ca nt go too far into the planetesimal disk . so , this configuration has a low , but non - negligible , success rate of 15% and can not be ruled out . : this is the configuration that presents the highest rate of success , considering all criteria together , so that it is highlighted in bold in table [ table1 ] . this is also the most successful configuration relative to criterion e , with a success rate of 65% independently of the disk setup . the evolution of neptune and of the rest of the planetary system is usually similar to that shown in ( * ? ? ? * figs 12 , 13 , and 14 ) , and satisfies the constraint identified in @xcite ( figure [ fig1 ] bottom right in this work ) . figure [ fig1 ] shows examples of the evolution of the systems corresponding to the descriptions given for each of the resonant chains discussed above ( with the exception of the 2:1 jupiter - saturn resonance , due to the considerations previously made ) . almost all evolutions depicted in figure [ fig1 ] fulfill the criteria a , b , c , and d. criterion c was not fulfilled only for the configuration 3:2 , 3:2 , 4:3 , 4:3 ( top left ) , where , although the eccentricity of jupiter was excited during the instability , it was quickly damped afterwards ( @xmath33 0.0075 ) due to the strong dynamical friction provided by the planetesimal disk . only the configuration 3:2 , 3:2 , 2:1 , 3:2 ( bottom right ) fulfilled all constraint together , including our new criterion e.    summarizing , when the new constraints on the evolution of neptune from @xcite are considered , the most probable initial configuration of the giant planets was 3:2 , 3:2 , 2:1 , 3:2 . the configurations 3:2 , 3:2 , 2:1 , 2:1 , and 3:2 , 3:2 , 3:2 , 3:2 , or 3:2 , 3:2 , 4:3 , 4;3 also appear possible , but less probable . this seems to be a consequence of the more relaxed orbital separation between the first and second ice giants . therefore , one would expect that any configuration having the first and the second ice giant in distance from the sun in a 2:1 mmr or similarly relaxed could perhaps also work . however , for simplicity , we decide to restrict ourselves only to the cases of @xcite , which we know the rate of success in filling criteria a , b , c , and d.    the analysis performed in this section does not take into account the timing of the planetary instability . with the location of the inner edge of the disk that we adopted , all evolutions are relatively fast , and lead to instabilities after a few 10s my . in the next section we analyze whether it is possible to get a late instability out of the 3:2 , 3:2 , 2:1 , 3:2 configuration , under some assumption on the geometry of the disk , still satisfying criterion e on the evolution of neptune . since @xcite presented the _ nice model _ , the existence of a late heavy bombardment of the terrestrial planets has been associated to the possibility that the giant planet instability occurred after hundreds of my of evolution @xcite . actually , the nice model , both in its initial version and from generic initial resonant configurations @xcite , is compatible with both early and late instabilities , depending on the distance of the inner edge of the original trans - neptunian disk and other disk properties ( mass , excitation , etc . ) . thus , the nice model does not prove the existence of a late instability . the timing of the instability can only be determined by looking at other constraints , such as the lunar impact chronology . unfortunately , the latter is still under debate . however , it is not obvious that the nice model can still be compatible with a late instability , if we fix a specific initial resonant configuration and impose a specific evolution of neptune . as explained in the introduction , if the disk is sufficiently close to neptune at the beginning , it is relatively easy that neptune migrates several au away before that the instability happens , but in this case the full dynamical evolution is rapid , as in the previous section . if the initial location of the disk is sufficiently far from neptune , the instability may happen late , but neptune is unlikely to migrate significantly before that the instability happens . moreover , in the case of the most relaxed resonant configurations ( e.g. 3:2 , 3:2 , 2:1 , 2:1 , which has neptune initially around @xmath024.5 au ) , the requirement that the inner edge of the disk is several au beyond neptune while the outer edge is at @xmath030  au reduces the disk to a ring . this seems implausible , particularly when considering the constraint that the disk should have hosted 1,000 - 4,000 pluto - size bodies @xcite , which would probably be problematic for the late instability as well , because the spreading of the narrow disk would not permit a delay .    if one could show that the most developed and constrained version of the nice model is incompatible with a late instability , this would be very important to clarify the chronology of events in the solar system . thus , in this section we try hard to find whether the instability could have happened late and under which conditions . we let the reader judge whether these conditions are realistic or plausible , with our considerations summarized in sect . [ conclusions ] .    for this goal , we proceed in steps . first , in sect . [ dr ] , we determine two critical values for the distance of the inner edge of the disk from neptune . one , @xmath34 , is the maximal distance that allows neptune to have a long - range planetesimal - driven migration away from the other planets , as in the previous section . the other one , @xmath35 is the minimal distance of the disk for which the five - planet system remains stable in resonance for all times . clearly @xmath36 and @xmath34 is larger than or equal to the neptune - disk distances considered in the previous section . the problem is that both @xmath35 and @xmath34 depend on disk resolution ( number of planetesimals and individual masses that are used to model the disk ) and that realistic values can not be used in practical simulations . so , we need to estimate trends to evaluate the values of the critical distances .    in a second step , in sect  [ de ] , we study how a flux of dust from the planetesimal disk , due to collisional planetesimal grinding and poynting - robertson drag , influences planet migration . this depends on total mass flux and individual dust size . the idea is that , under an appropriate dust - flux , neptune may migrate away from its initial location , so that its distance from the planetesimal disk decreases from @xmath37 to @xmath38 . at this point planetesimal - driven migration starts , as in the previous section . if the dust - driven migration is slow enough , the planet instability can occur late , even though it occurs only a few millions of years after the beginning of planetesimal - driven migration . this concept is then proven in section  [ dde ] , where we combine in a single simulation dust - driven migration and planetesimal - driven migration , and we obtain an instability after 300  my , with neptune s evolution satisfying criterion e. this simulation is of course performed with an unrealistically low disk - resolution , but is a proof of concept that , for whatever disk resolution , an appropriate dust flux can be found . we finally discuss what a realistic dust flux could be in sect . [ dgl ]      the interaction of a planetary system with a planetesimal disk can be very resolution - dependent , i.e. , sensitive on the number of planetesimals used to described the disk in the simulation and to their individual mass .    in some cases , as in those of section [ early ] , this problem of sensitivity on disk - resolution is not strong because the distance between the outermost planet and the inner edge of the disk is small , so that there are from the very beginning many planetesimals on unstable orbits that encounter the planet , force its migration and , consequently , bring many new planetesimals to planet - crossing orbit . this is the reason for which @xcite found similar results for test simulations conducted with 100 , 300 , 1,000 , 3,000 , and 10,000 planetesimals in the disk ( same total disk masses ) , so that they decided to use only 1,000 planetesimals in their final simulations . of course , even in these cases , the disk resolution has an influence on how grainy or smooth the planet migration will be , but does nt have a major influence on the amplitude and the timescale of the migration .    on the other hand , the resolution of the planetesimal disk becomes a real issue when the inner edge of the disk is far enough from the outermost planet . in this situation , the fraction of the disk that is unstable is small , mostly limited to the volume of the strongest mean motion resonances with the planet . the resonant planetesimals develop planet - crossing eccentricities and , with close encounters , force the planet to move a bit outwards . this moves outwards the resonances as well , destabilizing new planetesimals ( those which were not in resonance at the beginning but are swept by the resonant motion ) . depending on whether the mass of the new unstable planetesimals is larger or smaller than the mass of the original unstable planetesimals , the planet migration will accelerate or damp @xcite . in the damped migration mode , the planet eventually stops migrating while the resonances become empty .    if the number of planetesimals is small , the planetesimals are individually more massive to represent a disk with the same total mass . thus , the migration of the planet is more grainy because each encounter between the planet and a planetesimal causes a larger jump in the planet s semimajor axis . of course , there are fewer jumps , so that the planet migrates , on average , by the same amount . but the more jumpy migration of the planet causes the resonances to jump around as well , thus destabilizing more mass from the disk than in the case of a smoother planet migration with the same averaged range . because more mass is transferred from the disk to planet - crossing orbit , the migration of the planet can change , in some cases , from the damped to the accelerated mode .    with this issue in mind , we proceed in our study looking for two critical disk distances , for a series of disk resolutions . first , we look for the minimum distance ( @xmath1 ) that the inner edge of the planetesimal disk has to have from neptune ( the outermost planet for our case ) in order to maintain the planets stable for at least 400 my . second , we search for the maximal distance ( @xmath2 ) that neptune can have from the inner edge of the planetesimal disk to evolve as in section [ early ] . from these results we will extrapolate the values of @xmath1 and @xmath2 to realistic almost infinite resolution . once again we consider disks with fixed surface densities ( scaling as @xmath40 ) , because in this case we guarantee always the same amount of mass in the planets mmrs , independently of the resolution of the disk ( apart from stochastic sampling ) . the initial distribution of eccentricity and inclination of the disk was such that @xmath41 and @xmath42 ( the same initial distribution for all disks considered in this work ) . this is the simplest configuration that represents a cold planetesimal disk emerging from the end of the dissipation of the gas nebula . even if the disk was initially more dynamical excited not much differences in the overall evolution are expected because the mmrs inside the disk are very efficient in exciting the disk s eccentricities to planet - crosser values . so , in a very short timescale an initially cold disk would appear like an initially excited one . finally , we consider disks with resolution of 1,000 , 5,000 , 10,000 , 15,000 , and 20,000 planetesimals . the outer edge of the disk is fixed at 30 au . the distance @xmath30 between the inner edge of the disk and the planet ranges from 1 to 7 au , with increments of 1 au .    given the great demand of cpu - time required to simulate high - resolution disks , we only considered the most successful initial configuration in section [ early ] , with planets in the 3:2 , 3:2 , 2:1 and 3:2 mmrs . for this configuration , neptune is initially @xmath020 au . for each value of @xmath30 , we considered six different values of @xmath43 . .@xmath1 as a function of @xmath43 for the less resolved planetesimal disk ( 1,000 planetesimals ) . @xmath1 : minimum distance that the inner edge of the planetesimal disk has to have from neptune ( the outermost planet for our case ) in order to maintain the planets stable for at least 400 my . @xmath43 : density of the planetesimal disk , scaled as @xmath40 . the mass of the planetesimal disk , @xmath44 , in earth masses , for each pair @xmath43@xmath1 for the initial configuration 3:2 , 3:2 , 2:1 , 3:2 , with neptune at @xmath020 au and outer disk edge at 30 au is also shown for reference . [ cols="^,^,^,^,^,^,^,^,^",options="header " , ]     the results of the  proof of concept " simulation are shown in figures [ fig4 ] and [ fig5 ] . figure [ fig4 ] shows the evolution of the semimajor axis ( @xmath45 ) , the perihelion distance ( @xmath46 ) and the aphelion distance ( @xmath47 ) for all planets as a function of the time . figure [ fig5 ] brings snapshots of the ( @xmath45 @xmath39 @xmath48 ) and ( @xmath45 @xmath39 @xmath49 ) of the whole system for different times of the evolution ] . the evolution shown in these figures has an instability after about 300  my and the evolution of the planets is in a considerably good agreement with that deduced in @xcite from kuiper belt constraints . particularly , the instability happens when neptune is at @xmath50  au as required . moreover , as we can see from these figures , the evolution of neptune before the onset of the instability , follows closely the assumptions / predictions of tables [ table2 ] and [ table3 ] , and the dust mechanism of section [ de ] . the instability did not happen precisely at the time we targeted , due to the approximations inherent in our linear combination of different processes , which in reality have a non - linear interplay . nevertheless , the instability happens late enough to suggest that the association between a planetary instability and the trigger of the late heavy bombardment is indeed possible .      to assess whether the dust fluxes that are needed to migrate neptune into the disk ( see table [ table3 ] ) are reasonable , we turned to the simulation presented in the supplementary material ( s4 ) of @xcite . that simulation considered a disk with initially 60 earth masses , extended from 21 to 34 au the disk was quite excited in eccentricity and inclination within 27 au and dynamically cold beyond it . the planetesimals were assumed to have a size - frequency distribution ( sfd ) similar to that of the hot kuiper belt population in the inner part of the disk ( r @xmath4 27 au ) and similar to that of the cold kuiper belt population beyond it . both sfds were scaled up in number so that the masses in the inner disk and outer disk were respectively 1/3 and 2/3 of the total . the collision probabilities and impact velocities within particles in the inner part of the disk , the outer part or in different parts were evaluated to be @xmath51 km@xmath52y@xmath53 and 0.83 km / s , @xmath54 km@xmath52y@xmath53 and 0.25 km / s , @xmath55 km@xmath52y@xmath53 and 0.68 km / s , respectively . the collisional evolution simulations were conducted with the code boulder , developed in @xcite . the specific impact energy for collisional disruption was assumed to be 1/3 of that of strong ice @xcite . according to these simulations , the disk produced 32 earth masses of dust in 400 my . a total mass of 32 earth masses is big compared to the numbers in table [ table3 ] . on the other hand , not all this dust would survive all the way to neptune - crossing orbit , because the collisional cascade in the dust population comminute grains until they reach the ejection limit @xmath56 . this part of the collisional cascade is not modeled in the boulder code . also , the total amount of dust produced could be reduced by truncating the initial sfd to objects larger than some limit size ; in a new simulation where we retained only bodies with @xmath57 km in the initial sfd the total dust production is reduced to 15 earth masses . another possibility is to increase the specific energy of disruption for the bodies . in fact , @xcite revised the previous estimates of @xcite and concluded that the specific energy for disruption should be similar to that for strong ice in @xcite . a simulation with an initial truncated sfd but with a specific disruption energy 3 times larger produced only 7 earth masses of dust . even more simply , reducing the initial population in the disk would reduce the amount of dust produced , given that the dust production scales as the square of the number of planetesimals in the disk . in conclusion , we think that the dust fluxes listed in table [ table3 ] , required to cause a significant migration of neptune , are reasonable and in the good ballpark of what a massive planetesimal disk can produce .    a more relevant problem is the dust production timeline . in the boulder experiments described above , the dust production decays as @xmath58 . in the boulder simulations , the collisional evolution of the dust is not tracked . but the dust collisional lifetime should scale as the square root of the dust production . consequently , the total amount of dust of a given size ( or @xmath59 ) should evolve with time as @xmath60 . because the migration of the planet is linear in the amount of dust that it encounters , the evolution of its semimajor axis should be proportional to @xmath61 . this means that , of the total migration range spanned in -say- 400 my , 1/2 of it is covered in the first 100 my and 13% in the last 100 my . in other words , the time to reach @xmath2 is proportional to @xmath62 , instead of being proportional to @xmath63 as in the case of a linear evolution ( where @xmath64 is @xmath65 au ) . thus , small changes in the flux can make the curve @xmath66 to reach @xmath2 distance early or very late . obtaining a late instability around 400 my requires a more delicate tuning of the dust production than that done in sect . [ dde ] , where we considered that the dust production is constant and the collisional lifetime of dust infinite . it is very likely that , at the disappearance of the protoplanetary disk of gas , the giant planets of the solar system were all in mean motion resonances with each other , because this is the typical outcome of gas - driven planet migration @xcite . @xcite investigated which resonant configurations are more likely to reproduce the current orbital configuration of the giant planets via a phase of dynamical instability , as that described in the nice model . they considered four criteria of success and found a number of initial configurations that meet all of them in a non - negligible fraction of the simulations . all these configurations require the existence of a fifth giant planet , with a mass comparable to those of uranus or neptune , eventually ejected during the instability . more recently , @xcite showed that explaining and reproducing the current orbital structure of the kuiper belt requires that neptune had a specific evolution , reaching a heliocentric distance of @xmath028 au by planetesimal - driven migration before that the instability happened . thus , in section [ early ] of this paper , we have revisited the work of @xcite , adding a new criterion of success from the results of @xcite . of the configurations selected by @xcite , the one denoted 3:2 , 3:2 , 2:1 , 3:2 ( each integer ratio referring to the mean motion resonance between a pair of adjacent planets ) stands out as the more likely to fulfill the new criterion of success . other configurations ( the 3:2 , 3:2 , 2:1 , 2:1 , the 3:2 , 3:2 , 3:2 , 3:2 , and the 3:2 , 3:2 , 4:3 , 4:3 ) can not be excluded , but their success rate relative to the new criterion is significantly lower . an open issue with the nice model is the timing of the giant planet instability . if one wants to explain the late heavy bombardment of the moon as a spike in the impact crater rate @xcite , the instability has to have happened late @xcite , unlike in the simulations of section [ early ] , in which the instability occurred after a few 10s of my . it is not obvious that a 5-planet resonant configuration can remain stable for long time and become unstable then , particularly with a specific evolution of neptune as that required to sculpt the kuiper belt . thus , in section [ late ] we investigated whether a late instability from the 3:2 , 3:2 , 2:1 , 3:2 is possible and under which conditions . first , we looked for characteristics of the planetesimal disk that can maintain the planets in stable orbits for 400 my . we found that the evolution and stability of the planets are very sensitive to the resolution of the planetesimal disk . thus , we determined the minimum distance that the inner edge of the planetesimal disk has be away from neptune to keep the planets in stable orbit for 400 my as a function of the surface density of the disk and its numerical resolution . we named this distance @xmath1 ( table [ table2 ] ) . in a second step , we found a characteristic neptune - disk distance , which we named @xmath2 , from where neptune starts a long - range planetesimal - driven migration . thus , for the system to exhibit a late instability together with an evolution of neptune like the one that is needed to sculpt the kuiper belt , neptune needs to be extracted from the 3:2 , 3:2 , 2:1 , 3:2 resonance chain , and to migrate outwards over hundreds of mys until it reaches the distance @xmath2 . to understand how this could be possible , we introduced the concept of dust - driven migration , i.e. , migration due to the gravitational interaction of a planet with the dust produced in a distant planetesimal disk , spiraling inward by poyning - robertson drag . we have shown that a late instability can be achieved when considering a flux of dust ( @xmath67 ) ranging from 7 to 22 @xmath5 per 400 my , depending on the value of @xmath1 ( table [ table3 ] ) . continuing , in section [ dgl ] we evaluated how much dust could be produced from the disk - like configurations that we assumed and found good agreement with the dust flux described above . indeed , we found that the total amount of dust produced in 400 my may vary from @xmath07 @xmath5 to a maximum of about 32 @xmath5 depending on the size frequency distribution ( sfd ) and catastrophic disruption energies assumed for the planetesimal disk and its constituents . finally , we also found that the dust production rate expected for a self - grinding disk is not constant , but it decays as @xmath58 . the abundance of dust , which is also limited by the dust collisional lifetime , should decay as @xmath60 . thus , because the migration of the planet is linear in the amount of dust that it encounters , the evolution of its semimajor axis should be proportional to @xmath61 , which implies that small changes in the flux of dust can result in neptune reaching @xmath2 distance early or very late . therefore , obtaining a late instability around 400 my requires a more delicate tuning of the dust production .    summarizing , based on the work of @xcite , and considering the evolution of neptune required to explain the structure of the kuiper belt @xcite , our results indicate that the preferred initial configuration of the giant planets at the time of the dissipation of the gas nebula is the 3:2 , 3:2 , 2:1 , 3:2 , and that , although an early instability reproduces more easily the evolution of neptune required to explain the structure of the kuiper belt , such evolution is also compatible with a late instability . we thank an anonymous reviewer for useful comments and suggestions on the submitted manuscript . r.d . acknowledges support provided by grants @xmath682015/18682 - 6 and @xmath682014/02013 - 5 , so paulo research foundation ( fapesp ) and capes . a.m. acknowledges support by the french anr , project number anr-1313-bs05 - 0003 - 01 projet mojo ( modeling the origin of jovian planets ) . r.s.g . acknowledges his grant no . 307009/2014 - 9 from cnpq , conselho nacional de desenvolvimento cientfico e tecnolgico , brazil . the work of d.n . was supported by the nasa emerging worlds program . 99 batygin , k. , brown , m.  e. , & betts , h.  2012 , , 744 , l3 benz , w. , & asphaug , e.  1999 , , 142 , 5 bottke , w. f. , vokrouhlick , d. , minton , d. , et al . 2012 , , 485 , 78 brasil , p.  i.  o. , roig , f. , nesvorn , d. , et al . 2016 , , 266 , 142 brasser , r. , morbidelli , a. , gomes , r. , tsiganis , k. , & levison , f. h.  2009 , 134 , 1790 burns , j.  a. , lamy , p.  l. , & soter , s.  1979 , , 40 , 1 chambers , j. e.  1999 , , 304 , 793 deienno , r. , yokoyama , t. , nogueira , e. c. , callegari , n. jr . , & santos , m. t.  2011 , , 536 , a57 deienno , r. , nesvorn , d. , vokrouhlick , d. , & yokoyama , t. 2014 , , 148 , 25 deienno , r. , gomes , r.  s. , walsh , k.  j. , morbidelli , a. , & nesvorn , d.  2016 , , 272 , 114 fernandez , j.  a. , & ip , w .- h . 1984 , , 58 , 109 fernandez , j. a. , & ip , w. h. ,  1996 , , 44 , 431 gomes , r.  s. , morbidelli , a. , & levison , h.  f.  2004 , , 170 , 492 gomes , r. s. , tsiganis , k. , morbidelli , a. , & levison , h. f.  2005 , , 435 , 466 izidoro , a. , morbidelli , a. , raymond , s.  n. , hersant , f. , & pierens , a.  2015 , , 582 , a99 kaib , n.  a. , & chambers , j.  e.  2016 , , 455 , 3561 leinhardt , z.  m. , & stewart , s.  t.  2009 , , 199 , 542 leinhardt , z.  m. , & stewart , s.  t.  2012 , , 745 , 79 levison , h.  f. , morbidelli , a. , gomes , r. , & backman , d.  2007 , protostars and planets v , 669 levison , h. f. , morbidelli , a. , vanlaerhoven , c. , gomes , r. , & tsiganis , k.  2008 , , 196 , 258 levison , h.  f. , bottke , w.  f. , gounelle , m. , et al . 2009 , , 460 , 364 levison , h.  f. , morbidelli , a. , tsiganis , k. , nesvorn , d. , & gomes , r.  2011 , , 142 , 152 nesvorn , d. , vokrouhlick , d. , & morbidelli , a.  2007 , , 133 , 1962 nesvorn , d.  2011 , , 742 , 22 nesvorn , n. & morbidelli , a.  2012 , , 144 , 117 nesvorn , d. , vokrouhlick , d. , & morbidelli , a.  2013 , , 768 , 45 nesvorn , d. , vokrouhlick , d. , & deienno , r.  2014 , , 784 , 22 nesvorn , d. , vokrouhlick , d. , deienno , r. , & walsh , k.  j.  2014 , , 148 , 52 nesvorn , d.  2015 , , 150 , 73 nesvorn , d.  2015 , , 150 , 68 nesvorn , d. , & vokrouhlick , d.  2016 , , 825 , 94 malhotra , r.  1993 , , 365 , 819 malhotra , r.  1995 , , 110 , 420 morbidelli , a. , levison , h. f. , tsiganis , k. , & gomes , r.  2005 , , 435 , 462 morbidelli , a. , tsiganis , k. , crida , a. , levison , f. h. , & gomes , r.  2007 , , 134 , 1790 morbidelli , a. , bottke , w.  f. , nesvorn , d. , & levison , h.  f.  2009 , , 204 , 558 morbidelli , a. , brasser , r. , gomes , r. , levison , h. f. , & tsiganis , k.  2010 , , 140 , 1391 moro - martin , a.  2013 , planets , stars and stellar systems . volume 3 : solar and stellar planetary systems , 431 petit , j .- m . , kavelaars , j. j. , gladman , b. j. , et al . 2011 , , 142 , 131 . pierens , a. , & nelson , r.  p.  2008 , , 482 , 333 roig , f. , & nesvorn , d.  2015 , , 150 , 186 tera , f. , papanastassiou , d.  a. , & wasserburg , g.  j.  1974 , earth and planetary science letters , 22 , 1 thommes , e.  w. , duncan , m.  j. , & levison , h.  f.  1999 , , 402 , 635 toliou , a. , morbidelli , a. , & tsiganis , k.  2016 , , 592 , a72 tsiganis , k. , gomes , r. s. , morbidelli , a. , & levison , h. f.  2005 , , 435 , 459 walsh , k.  j. , & morbidelli , a.  2011 , , 526 , a126 walsh , k. j. , morbidelli , a. , raymond , s. n. , obrien , d. p. & mandell , a. m.  2011 , , 475 , 206
Abstract:

Answer: recent works on planetary migration show that the orbital structure of the kuiper belt can be very well reproduced if before the onset of the planetary instability neptune underwent a long - range planetesimal - driven migration up to @xmath028 au . however , considering that all giant planets should have been captured in mean motion resonances among themselves during the gas - disk phase , it is not clear whether such a very specific evolution for neptune is possible , nor whether the instability could have happened at late times . here , we first investigate which initial resonant configuration of the giant planets can be compatible with neptune being extracted from the resonant chain and migrating to @xmath028 au before that the planetary instability happened . we address the late instability issue by investigating the conditions where the planets can stay in resonance for about 400 my . our results indicate that this can happen only in the case where the planetesimal disk is beyond a specific minimum distance @xmath1 from neptune . then , if there is a sufficient amount of dust produced in the planetesimal disk , that drifts inwards , neptune can enter in a slow dust - driven migration phase for hundreds of mys until it reaches a critical distance @xmath2 from the disk . from that point , faster planetesimal - driven migration takes over and neptune continues migrating outward until the instability happens . we conclude that , although an early instability reproduces more easily the evolution of neptune required to explain the structure of the kuiper belt , such evolution is also compatible with a late instability .