Instruction: Write an abstract for the following acadamic paper.

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

Paper: the discovery of the planet around 51 pegasi ( mayor & queloz 1995 ) has spurred an avalanche of observations over the last 15 years , leading to a census of almost 700 confirmed planets to date . within this data we can find a variety of exotic and unexpected systems , such as planets with short orbital periods , extreme orbital eccentricities , or systems with strong resonant interactions . these observations have also inspired a wealth of theoretical work , on both the origin and structure of such planets . furthermore , planet searches have continued to push the observational thresholds towards lower masses , with current observations now probing below @xmath5 ( mayor et al . 2009a , b ; leger et al . 2009 ; howard et al . 2011 ; batalha et al . 2011 ) , with the hope of eventually finding planets that may convincingly be argued to be hospitable to life .    despite the wealth of data and interpretation now available , it is still difficult to draw together the different elements into an overall synthesis , in large part because of the heterogeneous nature of the observations . however , in the next few years a variety of observations offer the possibility to provide new constraints , which may allow us to frame a rigorous context for planet formation and evolution . recently , howard et al . ( 2010 ) and mayor et al . ( 2011 ) published the results of volume limited radial velocity surveys of nearby sun - like stars . the importance of these studies is that one can now begin to compare the relative frequencies of different types of planets , illuminating which are the generic product of the planet formation process , and which are the exotica . the results suggest that lower mass planets ( @xmath6 ) are even more common than gas giants at short periods , and that they appear to form a distinct population for orbital periods @xmath7  days . particularly striking is the fact that the low - mass planet distribution appears to show an upper edge in a plot of mass versus orbital period ; planets are found with masses between the completeness limit of the survey and an upper limit ranging between @xmath8 at @xmath9  day and @xmath10 at @xmath11  days . this is in contrast to the results of the predominant planet population synthesis models ( e.g. ida & lin 2008 ; mordasini , alibert & benz 2009 ) , which predict that this mass range ( @xmath12 ) should be liberally populated with partially - accreted giant planet cores ( `` hot neptunes '' ) that have migrated in from larger radii . in addition to the radial velocity surveys , the kepler satellite has begun to release the contents of its large - scale transit catalogue ( borucki et al . 2011 , lissauer et al . the majority of the present sample are still unconfirmed ` candidates ' , but a variety of tests suggest that the true planetary yield from this sample is expected to be high ( borucki et al . 2011 ) . the kepler sample once again has the advantage of being drawn from a relatively homogeneous program , and the distribution of planetary radii supports the aforementioned claim that the dominant species of short period planets corresponds to radii and masses below that expected from gas giant planets ( howard et al . 2011 ; wolfgang & laughlin 2011 ; youdin 2011 ) . these new surveys highlight the importance of short - period low - mass planets , as the low - mass planets now appear to outnumber their more massive brethren . the origins of such hot neptunes and ` super - earths ' is still poorly understood . most models of their formation ( alibert et al . 2006 ; terquem & papaloizou 2007 ; raymond , barnes & mandell 2008 ; mordasini et al . 2009 ; ida & lin 2010 ; mcneil & nelson 2010 ) adopt the position that they form at larger radii , in a manner related ( if not similar ) to the planets uranus and neptune in our own solar system , and then migrate inwards , to be trapped closer to the star . however , quantitative comparisons with the properties of observed systems remain unsatisfactory ( mcneil & nelson 2010 ) .    ultimately , many of the problems with these calculations stem from the properties of the migration models assumed . although the notion of gaseous disk migration has a long , physically motivated provenance ( e.g. goldreich & tremaine 1980 ; ward 1997 ) , the nature of the migration and when it occurs is still a subject of significant uncertainty ( see kley & nelson 2012 for a recent review ) . the principal reason for the revival of interest in migration models was the notion that hot jupiters were emplaced in this fashion ( e.g. lin , bodenheimer & richardson 1996 ) . however , it is becoming clear that these systems represent a small minority of all planetary systems ( e.g. howard et al . 2010 ; mayor et al . 2011 ) and may therefore reflect rather special circumstances . indeed , there is a growing accumulation of recent results demonstrating the misalignment of stellar spin and planetary orbits in many hot jupiter systems ( e.g. winn et al . 2011a and references therein ) . the simplest interpretation is that at least some hot jupiters did not migrate through a gas disk . furthermore , it is no longer clear that migration implies a monotonic inward motion ( e.g. paardekooper & mellema 2006 ; oishi , mac low & menou 2007 ; adams & bloch 2009 ) . with the large numbers of @xmath13 planets now seen at short periods , where the simplest migration models suggest they should not be , indicates that large radial movement need not be ubiquitous in planetary systems .    motivated by this renewed skepticism , we have chosen to revisit the question of in situ accumulation of rocky planets on small scales , and to explore the consequences of planet assembly without significant late - time gaseous migration . however , we will consider the possibility that significant radial migration is possible for smaller bodies , before they accumulate into planets , and so we will allow for disks with a wide variety of masses , discussed in   [ inventory ] . in   [ accum ] we then present the results of a series of traditional in situ planetary accumulation models on scales of 0.051  au , from the resulting disks . in  [ accretion ] we investigate the consequences of gas accretion onto the massive rocky cores that form . we address the role of jupiter mass planets in  [ hotties ] and evaluate the possible role of photoevaporation in   [ evap ] . the prototypical initial condition for assembly of a planetary system is the ` minimum - mass solar nebula ' ( weidenschilling 1977a ; hayashi 1981 ) , based on an estimate of the material required to form our own planetary system . the required surface density of solids ( based here on the hayashi estimate ) is @xmath14 with the gas surface density @xmath15 times larger . this applies interior to the ` snow line ' ( taken to be at @xmath16  au ) . exterior to the snow line the available material in solids is a factor of four higher , reflecting the fact that more volatile elements are able to condense and form solids as the nebular temperature drops . this enhancement brings the solid / gas ratio closer to the estimated metallicity of the sun , which presumably reflects the original composition of the nebular gas . if we integrate this surface density we infer a mass in solids of 3.3@xmath17 interior to 1  au , @xmath18 inside the snow line , and 25.5 @xmath17 within 10  au . this inventory is the minimum amount of mass required to match the solar system planets , but factors of two or three more are not unreasonable , if we assume a certain inefficiency in the process of converting nebular material into planets . in addition , many planet hosting stars are known to show super - solar atmospheric abundances ( gonzalez 1998 ; santos , israelian & mayor 2001 ; fischer & valenti 2005 ; udry & santos 2007 ) , often up to a factor of two . as such , we can reasonably argue for an increase in the above inventory of solids by a factor of five . thus , the simplest extension of the mmsn mass accounting could plausibly provide 17@xmath17 inside 1  au , 29@xmath17 inside the snow line , and 128@xmath17 inside 10  au . the above estimates are still very much in the spirit of the original mmsn model , in which solid material condenses out of a gaseous nebula and accumulates in situ to form planets from the local reservoir of material . however , the modern view of nebular and planetary evolution is considerably more dynamic . recent models have focused largely on the migration of planet - sized bodies , but it has been known for a long time that the dust and small rocks from which planets form are also potentially subject to significant radial migration ( weidenschilling 1977b ; nakagawa , sekiya & hayashi 1986 ; youdin & shu 2002 ; bai & stone 2010 ) . if a significant fraction of the solid material can migrate interior to 1  au while it is still in the form of small bodies , then the potential reservoir of material for the accumulation of terrestrial - class planets could be considerably larger . we will consider mass reservoirs up to 100 @xmath17 , corresponding to a radial concentration of mass by factors of at least 10 . the nature of the processes that produce such radial concentrations is still a subject of active research , and it is not guaranteed that the solids will preserve the same radial profile as the original gas disk . thus we will consider solid surface density profiles with a range of slopes . there is a long history of studies of terrestrial planet accumulation in our own solar system . we adopt the oligarchic model of kokubo & ida ( 1998 ) as our starting point , which provides an analytic prescription for a protoplanetary system consisting of a sequence of annuli , each dominated by a single planetary embryo . simulations of the last stages of gravitational interaction and accumulation of embryos into planets in this model , using n - body techniques ( chambers & wetherill 1998 , agnor , canup & levison 1999 ; chambers 2001 ; raymond , quinn & lunine 2006 ) , shows that configurations of three or four planets emerge naturally from a minimum mass solar nebula disk , although the low masses of mercury and mars in the solar system may require a more localised initial distribution of mass ( hansen 2009 ) . the production of multi - planet systems is fairly robust with respect to varying the amount of mass and the disk density profile ( raymond , quinn & lunine 2005 ; kokubo , kominami & ida 2006 ) , although the masses of the planets tend to scale with the reservoir mass . prior simulations of this process have generally restricted themselves to distances @xmath19  au because of the lack of any closer planet in the solar system ( and shorter timesteps add unecessarily to the computational burden ) . this is also often motivated on the basis of arguments regarding the chemistry of volatiles and a presumed reduction in the mass fraction that condenses out of the nebula . however , our knowledge of nebula dynamics and chemistry is sufficiently poorly constrained that it is prudent to examine the dynamical consequences of different disk masses regardless , and argue about the underpinnings later . thus , in the following , we will examine the assembly of terrestrial planet - class systems from power law surface density distributions of solid material extending down to 0.05  au and up to 1  au . we use the n - body code mercury ( chambers 1999 ) to simulate the final stages of accumulation of planetary embryos , born of the oligarchic model , into planets around a @xmath20 star . we have performed a series of simulations of disks that extend from 0.05 to 1  au , while varying the total mass and surface density profile , to assess whether such an initial distribution can produce the planetary systems observed in a homogenously selected sample , such as that by howard et al . ( 2010 ) . in order to resolve the orbital dynamics at the inner edge of the disk , we use timesteps of 12  hours ! fortunately , most of the interesting dynamics in these collisional accumulation simulations occurs within a fixed number of orbits , so that the compact systems simulated here evolve very rapidly in absolute terms . we evolve our systems to 10  myr , which is sufficient to produce dynamically stable final results . we have verified this by evolution of selected systems to 100  myr and 1  gyr , and have observed little significant change . results of these and subsequent simulations are collected in the table  [ simparam ] . we consider three surface density profiles , all of the power law form @xmath21 , and examine models with @xmath22 and @xmath23 . this encompasses a wide range of behaviour and is taken to represent the range of uncertainty in the initial conditions , given the possibility that this might be sculpted by radial migration of solids from large scale to small . with each choice of @xmath24 , we normalise the surface density to have a particular mass of solids interior to 1  au . we choose three different representative normalisations , of 25 , 50 and 100  @xmath17 . in each case , we start with the surface density profile consisting of a series of protoplanetary embryos realised according to the oligarchic distribution described in equations  17 & 18 of kokubo & ida ( 2002 ) , spaced out between 0.05  au to 1  au . for the disk masses considered here , which are larger than usually assumed , the mass of the initial oligarchs can sometimes reach several @xmath17 , and one might be concerned that such bodies are massive enough that their mutual perturbations would invalidate the assumption of mutually independent radial annuli . in appendix [ democratic ] , we have also investigated models in which we restrict the mass of oligarchs to be less than @xmath25 , and we found no meaningful difference in the eventual outcomes . figures  [ snap0][snap2.5 ] show the characteristic evolution of three different disks . all have the same total mass ( @xmath26 ) within 1  au , but distributed according to surface densities with different slopes . the disk with the flattest profile ( @xmath27 ) is shown in figure  [ snap0 ] . in this case , the masses of the initial oligarchs increase rapidly with radius , scaling with the area of individual annuli . the consequence of this is that the most massive planets at the end are located from 0.51  au . this is somewhat too large a radius to match most of the observations , which , at the time of writing , are most complete for orbital periods of 100  days or less . figure  [ snap1.5 ] shows the evolution for an intermediate slope ( @xmath28 ) . here we find that the most massive final planets are found between 0.2 and 0.5  au . although the original oligarchs still increase in mass with semi - major axis , there is sufficient mass on smaller scales so that the faster dynamical evolution results in the accumulation of mass into a few large , comparable - mass objects by 10  myr . figure  [ snap2.5 ] shows the evolution of a disk with a steep density profile ( @xmath29 ) . in this case , the oligarch masses decrease with semi - major axis , and the final result retains many of the characteristics of the original mass profile , with the largest planets having the shortest orbital periods . the three histories shown here demonstrate the characteristic behaviour of the models  steep density profiles result in a distribution of planet masses that decreases with increasing orbital period , while flatter profiles move the maximum planet mass to larger scales . the stochastic nature of n - body systems means that there is some variation in the total number and division of mass from one simulation to the next , but the general distribution of mass with semi - major axis is reasonably constant .    for a fixed density profile , the final planet masses will also , obviously , scale with the total mass budget . figure  [ param2_4 ] shows the mass and orbital period of the largest final planet for several simulations with @xmath28 ( the choice which comes closest to matching the observations ) . the open circles indicate a total disk mass of 100@xmath17 , the open triangles indicate @xmath26 and the crosses indicate a total mass of @xmath30 . the filled circles are the observed periods and masses of the largest bodies in systems detected by howard et al . ( 2010 ) ( jupiter mass planets will lie off the top of this plot ) . we see that the known super - earth candidates require total disk masses in the range @xmath31 .    as noted in  [ inventory ] , even the optimistic estimates of the static nebula solid inventory are not sufficient to meet this requirement of @xmath0@xmath1 needed to match the observed planets . such massive planetesimal disks are more plausible if we allow for the possibility of migration of solids from several au to small scales . another generic consequence of in situ accumulation of rocky planets is that the resulting systems are almost always multiple , and the degree of multiplicity and clustering can be characterised by a variety of statistics ( see chambers 2001 for a summary , and appendix  [ stats ] ) that were initially developed to analyse the formation of the terrestrial planets of our own solar system . this proves a useful test of the in situ accumulation hypothesis , as we would then expect the observed systems to be primarily multiple , and to be consistent with the statistical distribution of the simulated systems . a complication for such a comparison is that the observations are likely to be still severely incomplete ( howard et al . 2010 ) and there may be many additional planets lurking in the observed systems , but we can hope for this to improve with time .    however , there are some systems for which the noise levels are low enough to demonstrate the existence of multiple systems already . the howard et al sample contains several multiple systems , although some contain planets beyond 1  au , and so are not directly comparable to our calculations . however , there are two systems , 61 vir ( originally announced by vogt et al . 2009 ) and hd69830 ( originally announced by lovis et al . 2006 ) which contain multiple planets within the period and mass range under discussion .. ]    figures  [ param2 ] and [ param3 ] show the comparison of these two systems with the corresponding simulation results . in order to ensure a proper comparison , we do not count planets in the simulations whose radial velocity amplitude would be @xmath32 . figure  [ param2 ] plots two statistical measures , @xmath33 and @xmath34 , of how the mass in the system is distributed with semi - major axis . the observed systems hd69830 and 61 vir are shown as well , where the range indicates an uncertainty in the global system mass because of the unknown inclination angle of the system to the line of sight , assuming coplanarity to within 5 degrees . again , we see that the @xmath35 model produces systems with similar values of both mass concentration and distribution in eccentricity and semi - major axis , whereas the shallower and steeper slopes can not do match all observed statistics simultaneously . to some extent this is not surprising , given the qualitative features of figures  [ snap0][snap2.5 ] , although it is a more quantitative statement . figure  [ param3 ] compares the measures of angular momentum deficit , amd , ( a measure of eccentricity and inclination spread ) and mass - weighted mean semi - major axis . the spread on the observed systems account for randomly varying the observed eccentricities according to the observed error bars , and allowing the unknown mutual inclinations to fluctuate randomly by 5 degrees . again , we see that the @xmath35 models bracket the observations . it is striking that both the 61 virginis and hd69830 parameters overlap with the simulated values , suggesting that it is possible to realise these observed systems within this model . the spread in the model statistics derives primarily from two sources . the variation in the mass concentration statistics ( figure  [ param2 ] ) reflects the degree to which the mass is spread evenly amongst the surviving planets , or concentrated in one or two more massive bodies . the spread in the amd ( see appendix  [ stats ] ) seen in figure  [ param3 ] is caused primarily by one of the surviving planets having significant ( @xmath36 ) inclination with respect to the others . if such conditions apply to the observed systems , then their parameter range will be expanded from that shown here , in which we limited the amount of assumed mutual inclinations . the outlier systems in the simulation results usually result from an incomplete accounting due to one or more of the final bodies being too small to meet our radial velocity criterion . the result of these accumulation models suggests that both the 61 vir and hd69830 systems are consistent with in situ accumulation of rocky planets from disks of sufficient mass . it also suggests that many of the other systems in the howard et al sample should harbour planets not far below the current detection limits , as the ratio of the two most massive planets in the simulated systems ranges from 0.40.98 , with a median value of 0.78 . the ability of the in situ accumulation model to match the properties of the observed low mass planet population is encouraging , but the masses of the rocky bodies that result suggest that a reconsideration of the issue of gaseous accretion is warranted . in the traditional model of gaseous giant planet formation , the first stage is the collisional assembly of a rocky core , which eventually gravitationally captures gas from the surrounding nebula ( e.g. pollack et al . 1996 ) when it reaches as mass @xmath3 . the final planet masses in   [ accum ] exceed this threshold , albeit in a much different region of the disk than the original cores of jupiter and saturn . therefore , the assumption that these planets are wholly rocky is likely to be in error . the question of gas giant formation on small scales was widely reconsidered upon the discovery of 51  pegasi and other ` hot jupiters ' . however , a consensus quickly emerged that these planets most likely formed on larger scales and migrated to their currently observed positions . part of this consensus was the finding , from detailed calculations of the core accretion process , that in situ formation was unlikely ( bodenheimer , hubickyj & lissauer 2000 ; ikoma , emori & nakazawa 2001 , rafikov 2006 ) . however , these findings rest largely on the assumed implausibility of assembling a massive enough rocky core on small scales , given traditional estimates of the rocky material reservoir . if the solid mass available is substantially enhanced by inward radial transport of material , then such massive cores are no longer an outlandish idea . furthermore , on these scales , the assembly timescales of the planets in   [ accum ] is short , often @xmath37 years , and the system configurations are largely in place by @xmath38  years . thus , it is quite conceivable that cores of several earth masses are in place before the gas is removed from the system . as a demonstration of this , figure  [ acc ] shows the temporal history of the accumulation in a representative simulation with total disk mass in solids of 50 @xmath17 interior to 1  au . the final configuration has 4 planets in the ` super - earth ' mass range interior to 1 au . the assembly of all four planets is largely completed within 1  myr and so it is quite likely that the planets achieve a significant mass while there is still gas present in the system .    on scales @xmath39  au , detailed calculations suggest that the energy transport in a hydrostatic atmosphere above a rocky core is regulated by convection ( ikoma et al . 2001 ; rafikov 2006 ) . however , these authors also note that traditional estimates of the critical core mass on these scales become uncertain because the hill sphere radii of the resulting cores begin to approach the size of the disk scale height  i.e. a core could potentially open a gap in the disk before the threshold mass is reached . for a disk temperature profile given by @xmath40 this is @xmath41 . the criterion for the planet to open a gap in the gas disk is to set the hill diameter of a planet of mass @xmath42 equal to the gas disk scale height @xmath43 , which yields @xmath44 this threshold mass is less than many of the simulated planets in   [ accum ] and most of the observed planets in the howard et al thus , if the assembly of these planets takes place in the presence of gas , accretion is likely to occur . however , the final outcome will not be a traditional gas giant planet . the normal core accretion sequence of events proceeds through a sequence of hydrostatically supported atmospheres until runaway gas accretion occurs ( see rafikov 2006 for a review ) . however , for planets on these small scales , the hydrostatic atmosphere becomes tidally limited before instability sets in , and the opening of the gap limits the amount of gas accretion . therefore , tidally limited accretion from the disk will produce a different , gas - starved , branch of the core accretion family than is found on larger scales .    for a gas disk of surface density @xmath45 , the gas mass within an annulus of half width equal to the hill sphere of a core of mass @xmath42 is @xmath46 where @xmath47 is the gas disk mass inside 1 au , and @xmath48 is the orbital semimajor axis of the core . the value @xmath49 corresponds to the gas mass in the standard hayashi minimum mass solar nebula . the relation between this mass and the solid mass will depend on whether or not solids migrate radially . if the gas is eventually accreted by the planet , then the final mass of the planet will result from the combination of the two estimates . in the case where both the gas and solids follow a profile with @xmath50 , this yields a constant core mass fraction for the final planet , @xmath51 suggesting core mass fractions @xmath52 , for ranges of disk mass from 13 the mmsn value . if the gas disk had a different profile , the core mass fraction will vary with orbital period . for a gas disk profile @xmath53 , the mass - period relation would be @xmath54 it is striking that this estimate seems to naturally produce planets in the right range of masses and orbital periods . for example , at the fiducial period of 50 days , we derive a gap opening core mass of 7 @xmath17 which results in a total mass of 2149@xmath17 , for gas disks within a factor of three of the mmsn . this naturally covers the range of the observed masses , with most of the currently observed planets falling in this intermediate class with comparable amounts of rock and gas ( see figure [ pm ] ) . there is potentially also an intermediate regime , where the planet core is large enough to maintain a gravitationally bound atmosphere but not big enough to open a gap in the gas disk . as the disk gas dissipates , the planet may retain enough of this gas to have a significant observational effect . this is the case treated by rafikov ( 2006 ) , whose estimate of the core instability threshold is really a criterion for the gravitionally bound gaseous atmosphere to be some fraction of the core mass . if we adopt an atmosphere fraction of 10% , rafikov s estimate ( in the case of the inner disk  his equation  60 ) , amounts to @xmath55 this turns out to be very similar to equation  ( [ gapmaker ] ) , so that there are probably few planets in the intermediate regime . if a planet accretes a non - negligible mass fraction of gas , the accretion is likely limited by the opening of a gap . thus , we envision two different outcomes for in situ rocky accretion , depending on the amount of mass available . traditional terrestrial class planets , such as those in the solar system , emerge from disks below some mass threshold so that the time to assemble bodies of several earth masses is longer than the gas dispersal timescale ( @xmath56  years ) . for larger mass disks , cores form that are massive enough to bind gaseous atmospheres of significant mass before the dispersal of the disk . this process is also self - limiting , as they will accrete enough gas to actually open a gap in the disk , which limits their masses to the observed values . it does not depend sensitively on the evolution of the background gas disk  a difference from the traditional ` failed core ' model , often invoked to explain these objects ( e.g. mordasini et al . 2009 ; ida & lin 2010 ) .      to determine the threshold of disk mass that divides the gas - poor and gas - rich systems , we have simulated the assembly of disks of various masses for an @xmath35 profile , and find that the threshold disk mass is @xmath57 of solids inside 1  au . for disk masses above this value , one or more planets can assemble , within 1 myr , to masses that satisfy equation  ( [ gapmaker ] ) , and therefore should accrete gas . for lower mass disks , planets either never reach the threshold value or do so when the gas has dispersed , leaving behind rocky cores . figure  [ pm ] shows the overall picture that emerges , along with the observed planets in the howard et al . sample . the solid and open points indicate planets and planet candidates . the dotted line indicates the threshold mass for gas accretion ( eqn . [ gapmaker ] ) . planets above this dotted line will frequently harbour non - negligible gaseous envelopes . the exact size of the envelope depends on the mass in the gaseous disk . the long - dashed line indicates the final mass if the gas accreted is given by equation  ( [ gaplimit ] ) , and assuming a gas disk that contains @xmath58 inside 1  au ( equivalent to four times the mmsn value ) . the region between these two lines is approximately the domain of the neptune - like planets born in situ . the short - dashed line running diagonally from upper left to lower right divides the parameter space into regions where accretion of solid bodies is efficient ( to the left of the line ) and not efficient ( to the right ) . the in situ solid body accumulation process is largely conservative of mass on scales smaller than this , whereas planets that lie above the line are more likely to eject smaller bodies than to accrete them . these planets will not grow significantly , but rather will migrate inwards until they reach the accretion / ejection boundary given by the short - dashed line in figure  [ pm ] ( murray et al . further accretion will drive planets inwards approximately along the boundary ( i.e. they will grow and migrate inwards ) . the fact that there are few objects located near the intersection of these lines suggests that concurrent migration and accretion in this fashion is not a significant contributor to the demographics of the observed populations . the accretion of gas will change the assembly history , because one or more bodies will increase in mass and therefore the gravitational perturbations on other bodies in the system will be stronger . to account for the effects of this , we repeated the in situ assembly simulations using the @xmath35 profile , but we stop the simulation at 0.1 , 0.3 and 1  myr to compare the masses of all bodies to the threshold value of equation  ( [ gapmaker ] ) . if any body has crossed the threshold in the interim , we instantaneously increase the mass according to a model for gas accretion , and then continue the simulation . this model is not fully self - consistent as it does not treat the gravitational influence of the gas disk before it is accreted , nor does it account for any difference in the specific angular momentum of the accreted gas , but it does address the issue of how the rapid mass growth of some bodies affects the gravitational interactions that regulate the assembly of the n - body system into a final system of planets . after the last gas accretion episode the systems are run to ages of 10  myr . we examined two models for gas accretion . in the first , we started with a rocky disk of mass @xmath26 , @xmath35 and assumed gas was accreted according to equation  ( [ gaplimit ] ) , i.e. , from a gas disk that also obeys @xmath59 , and with a mass three times the mmsn value . this gas disk mass is chosen to produce planets that trace the upper edge of the observed lower mass planet population . in the second model , we start with a lower mass rocky disk ( @xmath60 , @xmath35 ) and assume gas is accreted according to equation  ( [ gasscale ] ) , i.e. , from a gas profile @xmath53 , with a mmsn normalisation . this is chosen to capture the median parameters of the observed low mass population . note that the gas and solid disks have different power law slopes , which is possible if the solid inventory was assembled by radial inward migration , and not in situ deposition . figure  [ param4 ] shows the statistical comparison of these simulations with the observed parameters of 61  vir and hd69830 . both sets of simulations show an inward shift in terms of the mass - weighted semi - major axis . this is a result of the tendency of planets with shorter orbital periods to assemble more rapidly and therefore they are more likely accrete gas . the end result is that including gas accretion in the model tends to steepen the mass profile . we have therefore repeated the experiments using a slightly shallower rocky profile of @xmath61 , using two different normalisations ( total masses of 50 @xmath17 and 80 @xmath17 ) . in both cases , the gas accretion is given by equation  ( [ gasscale ] ) . the results are shown in figure  [ param7 ] . with the use of a steeper profile , both 61 virginis and hd69830 remain within the parameter space of the simulations , although the simulations are now biased to larger semi - major axis ( because there is a larger fraction of mass at large radii and so it is easier to assemble the critical core mass ) . the overall conclusion to draw is that simulations can produce a wide range of systems , depending on the original slope and the amount of gas accreted , bracketing those observed . the dynamical effect of the gas accretion is to accelerate the orbital instabilities and dynamical interactions that characterise the final stages of terrestrial planet assembly . figure  [ ta22 ] shows an example of one system ( drawn from the @xmath35 case ) in which three different planets cross the threshold given by equation  ( [ gapmaker ] ) , two at 0.1  myr , and a third at 0.3  myr . close examination of the various curves shows that a rapid increase in mass leads to orbital instability in neighbouring , less massive bodies . the ultimate result is frequent collisions , and occasionally accretion onto the star or ejection from the system . this example also shows a common occurrence in which two of the planets which accreted gas eventually merge . the end result of this simulation is that we are left with three planets , two of which possess substantial gas fractions , although somewhat less than indicated by equation  ( [ gasscale ] ) because of the accretion of substantial amounts of rocky material during the final clearing stages . this plot also shows why the mass - weighted semi - major axis @xmath62 moves inwards for these simulations , relative to the purely rocky ones  the innermost planets grow the fastest and accrete gas earlier . this shifts the mass distribution to shorter radii . nevertheless , accretion occurs for planets at larger radii in some simulations as well , and we also find collisions between planets on larger scales too . this offers an interesting avenue to get the upper envelope of the masses in figure  [ pm ] even if the gas disk is not quite as massive as assumed in equation  ( [ gaplimit ] ) . our model thus far has addressed the nature of the low mass population at short periods . the howard et al . ( 2010 ) sample also contains a handful of genuine jovian mass planets with orbital periods @xmath6350  days . these must arise in some other fashion . the original model for such planets was migration through a gas disk , but recent observations suggest the situation may be different . for planets with orbital periods of a few days , tidal capture from eccentric orbits might be an alternative model ( rasio & ford 1996 ; matsumura , peale & rasio 2010 ) . however , this seems less likely to produce a planet with orbital periods @xmath64  days , because tidal forces weaken rapidly with distance . the planetary system around the star 55  cancri may offer insights because it appears to contain members of both the low mass and high mass population . we have investigated several hypotheses to see which are consistent with the existence of the 55  cancri system . * can the disk of rocky material we invoke play the role of a dissipative medium in capturing an eccentric jupiter ? the model in which a high eccentricity planet is tidally captured was originally motivated by the onset of dynamical instability , which can potentially occur late in a planetary system s lifetime ( rasio & ford 1996 ) after the removal of the gas disk . we have therefore examined whether a rocky planet system as described in previous sections can capture and circularise a jovian mass planet injected into the system on a highly eccentric orbit . the short answer is that this model is wildly unsuccessful , because the mass of the injected planet is comparable to , or larger than , the entire mass budget of the disks we have considered . thus , the end result of this model is that the injected planet remains on an orbit with a slightly decreased semi - major axis , and partially circularised eccentricity , having ejected essentially all of the mass of the initial disk . to reduce the semi - major axis and eccentricity to the level required by observation would require more mass , and would leave very little surviving mass to make the other planets in the 55  cancri system . + * it is possible to generate rocky planets of jovian mass , given a large enough mass reservoir . in a model where the solid material is laid down from a marginally gravitationally stable @xmath65 surface density profile ( appendix  [ gi ] ) , the resulting innermost planet is quite massive , @xmath66 , with orbital periods @xmath64  days ] . this opens the possibility for the formation of rocky planets of saturn mass in short period orbits . at present there is little evidence for the existence of such objects in the sample discovered in transit surveys , although that sample is heavily biased to shorter orbital periods , and is possibly dominated by the tidal capture of genuine gas giants from eccentric orbits . these simulations do not , however , faithfully reproduce the lower mass planets found in the 55  cancri system . nevertheless , if some fraction of protoplanetary disks are significantly sculpted by gravitational instability , this class of planets may eventually emerge in a sample like that of the kepler mission . * a hybrid version of the model in   [ accum ] is to simply assume the giant planet was in place in the beginning through some other , unspecified , process ) in the case of the q=1.5 gas disk , which is much more massive . however , the rest of the 55 cancri system is consistent with a less massive disk . ] , and then to repeat our calculation of in situ accumulation with this additional component . the end result of such a simulation actually produces systems qualitatively similar to that of 55  cancri , as shown in figure  [ 55match ] . we are able to produce both an inner planet and a system of planets further out . in the event of no gas accretion , the resulting planets are too small compared to the observed planets , but they do cross the threshold of equation  ( [ gapmaker ] ) and should grow to appropriate mass . figure  [ 55match ] shows this for the analogue of 55  cancri  c. the conclusion is that the presence of a jovian mass planet in the disk does not impede the in situ assembly of planets along the lines we have investigated , and that the resulting dynamical environment can produce systems with similar organisation to 55  cancri . jovian mass planets can also potentially influence the assembly on sub - au scales even if they orbit farther out ( such perturbations are often advanced as potential causes for the low mass of the asteroid belt , and possibly mars  e.g. chambers 2001 ) . there are several systems in the howard et al . sample which contain a gas giant at long orbital periods and a lower mass planet at shorter orbital periods . to investigate the degree to which an external jupiter could influence the assembly histories considered here , we have repeated the @xmath35 , 50 @xmath17 runs in the presence of a planet of mass 351  @xmath17 , orbital period 493  days , and @xmath67 . this is chosen based on the planet hd114783b , as an example . we find that planet assembly is largely unaffected interior to @xmath68  au ( orbital periods @xmath69  days ) . the bulk of our analysis applies to smaller scales and shorter periods , so that the general character of the solutions is largely unaffected by the presence or absence of jupiters on scales @xmath70  au . this is consistent with the approximate separation between jupiter and the asteroid belt in the solar system and the rule - of - thumb that a gas giant does not significantly reduce planet assembly interior to it s 3:1 mean motion resonance . it has been suggested by numerous authors that short period planets may experience significant atmospheric mass loss driven by irradiation from the central star ( e.g. lammer et al . 2003 ; yelle 2004 ; lecavelier des etangs et al . 2004 ; tian et al . 2005 ; murray - clay , chiang & murray 2009 ) . as a result , we need to consider how much of a limit this will place on the accretion or retention of gas by planets at short orbital periods . empirical estimates of the evaporation rate for exoplanets are based on the so - called ` energy diagram ' , in which the gravitational binding energy of the planet ( @xmath71 ) is compared to the energy input ( @xmath72 ) from x - ray and euv radiation ( lecavelier des etangs 2007 , ehrenreich & desert 2011 ) , after accounting for various uncertain efficiency factors . using the estimates of ehrenreich & desert for the life - time averaged euv flux for the host star , we can rewrite the formula for the characteristic evaporation time , @xmath73 , as a mass - period relation @xmath74 where @xmath75 is the efficiency factor for converting incident radiation into mass loss ( set equal to unity here , so that this is a very generous estimate ) . the meaning of this equation is that planets with masses above this threshold are too large to experience significant dimunition through the evaporation process , under even very optimistic assumptions . we see that this relation has the wrong slope to explain the edge in figure  [ pm ] , and will likely affect only the planets with periods @xmath76  days . the slope may be changed if we postulate a particular relationship between mass and radius , but not by much , as the mass - radius relationship for planets with gaseous envelopes is relatively flat in this mass range . therefore , evaporation is unlikely to be responsible for the shape of the mass - period relation in figure  [ pm ] . nevertheless , it is likely to be important for the lowest mass , shortest period planets . in cases where the mass given by equation  ( [ evaplim ] ) is greater than that given by equation  ( [ gaplimit ] ) , the planet may very well accrete some gas , but it will likely lose it again over the course of time . thus , we may expect to see a transition from significant gaseous envelopes to thin , remnant envelopes below some period . estimating the period is somewhat uncertain because of the rather unconstrained nature of @xmath75 , but it likely lies in the range of 3.7  days ( @xmath77 ) to 8.5  days ( @xmath78 ) . motivated by recent observational challenges to the prevailing paradigm of exoplanet migration , we have re - examined the role that in situ accumulation might play in the assembly of planets on small scales , with masses ranging up to several tens of earth masses . we have demonstrated that the observed properties of planets in the  hot neptune``/super - earth '' mass range , orbiting solar mass stars with periods of days to several months , can be reproduced with a simple model of in - situ accumulation of rocky planets from traditional oligarchic conditions , without any additional radial migration . however , in order for this to be the case , the original amount of rocky material required is @xmath79 within 1  au . this is somewhat in excess of that expected from traditional estimates of the mass profiles of protoplanetary disks , but could be achieved if there was significant radial migration of solid material from the outer parts of disks to the scales of interest , before the migrating bodies achieved the mass of planetary embryos . thus , our calculations argue that migration of solid material did indeed occur , but that it served to rearrange the mass budget of the nebula before the final stages of planet assembly . such a position is not without physical foundation , as the timescale for small bodies to migrate over au scales through the action of aerodynamic drag forces can be as small as @xmath80  years ( weidenschilling 1977b ; nakagawa et al . 1986 ; youdin & shu 2002 ; bai & stone 2010 ) . a further consequence of such a mass concentration is that the collisional assembly of planets is fast enough that bodies of several earth masses can form within 1  myr of the original material being emplaced . if this growth occurs in the vicinity of a gas disk , accretion and gravitational retention of gas from the nebula is possible . the accretion of this gas is limited , however , by the fact that it is easier to open a gap in the gas disk on smaller scales , and the accretion may be quenched before the accretion flow enters the traditional runaway regime in the core instability scenario ( e.g. rafikov 2006 ) . this process provides a natural mechanism for the production of planets with comparable amounts of gas and rock , intermediate between the traditional classifications of rocky terrestrial planets , and gas - dominated giants . figure  [ mx ] shows the heavy element fraction of planets of various mass formed in our simulations , compared to published estimates for observed hot neptunes around solar mass stars . it should be noted that many of the published systems do not strictly fall within the bounds of our model ( e.g. hd149026b ) in terms of the parameter space of figure  [ pm ] . furthermore , our crude estimate of a tidally - limited gas accretion ignores a lot of potentially interesting physics that determines the final mass of a planet embedded in a gas disk . in particular , planets that open gaps in disks can continue to accrete gas if the disk viscosity is high enough to allow flow into the gap ( e.g. artymowicz & lubow 1996 ; bryden et al . 1999 ; kley 1999 ) . based on the formulae of lin & papaloizou ( 1979 ) or taniguchi & ikoma ( 2007 ) , we estimate that disks with @xmath81 ( shakura & sunyaev 1973 ) may result in further gas accretion across the gaps we have estimated here . nevertheless , it is interesting that the simple considerations outlined above provide a natural framework for producing planets of this observed class . it is particularly encouraging that the objects on this plot most likely to fall within our model are the kepler-9 planets , which match the results of our simulations quite well . the quantitative agreement between observations and model is in striking contrast to the failures of traditional migration - based models ( alibert et al . 2006 , terquem & papaloizou 2007 ; ida & lin 2010 ) to match the same observations ( mcneil & nelson 2010 ) . the neptune mass analogs in these calculations formed further out and represent ` failed cores ' of giant planets , whose inward migration due to disk torques quenches the gas accretion that would produce a traditional gas giant . multiplicity is often important in the survival of such objects , as strong mutual gravitational interaction can serve to slow the migration and avoid spiralling into the central star . however , the effects of such strong interaction imprints itself on the final system , even after the loss of gas and further dynamical evolution ( ida & lin 2010 ; mcneil & nelson 2010 ) . recent observations of hot jupiters have cast doubt on the migration of gas giants as well , suggesting that perhaps the very idea of planetary migration needs to be extensively reviewed ( although we stress that the inventory of solids invoked in our model most likely still requires significant radial migration of solid material ) . we should be clear , however , that the model presented here is not complete , as we have not addressed the details of how these systems are related to those with hot jupiters or jupiters on more traditional scales . our aim is to better understand the individual physical processes at work before attempting a complete synthesis of all known planetary systems . furthermore , it is now evident that the lower planet mass systems are the most generic kind , and that stars which host hot jupiters are a small minority . it may thus be useful to consider our models independently of those developed for hot jupiters , which could represent only one very special case amongst the pathways to planetary system formation . our original discussion was purposely limited to the sample of howard et al . ( 2010 ) because of the homogeneity of selection . mayor et al . ( 2011 ) recently released a catalog of detections from a homogeneously selected sample of their own . in addition to 61  vir and hd  69830 , their sample contains seven more stars with at least three sub - jupiter mass companions . figure  [ param6 ] shows a comparison of these systems with the same models as discussed before . the model range encompasses nearly all of these , with hd  10180 ( lovis et al . 2011 ) , hd  31527 and hd  136352 located closest to the bulk of the simulation points . several systems , such as hd  40307 and hd  39194 , could be brought closer to the main locus of points if additional planets are discovered in these systems with longer orbital period ( which seems likely in the case of hd  40307 , since mayor et al . 2009a also report the detection of a linear trend ) . conversely , hd  134606 has an anomalously large @xmath82 , which could be made consistent if another planet is found in the system with orbital period @xmath83  days .    as a check on the howard et al . sample , figure  [ pmall ] shows the harps ( mayor et al . 2011 ) and kepler - based exoplanet samples in the orbital period - mass plane , along with the proposed physical classification based on our model . at the lower right , at masses below those given by equation  ( [ gapmaker ] ) , planets are not large enough to accrete significant gas and thus are likely to be truly rocky planets . this is a regime that is poorly explored as yet , although current radial velocity limits are slowly inching lower . above this limit , and at orbital periods @xmath84  days , we expect to find cores which have accreted substantial , but not overwhelming , gaseous envelopes . within our model , the bulk of the planets in this period and mass range are expected to be more like neptunes than scaled - up rocky planets ( a.k.a . super - earths ) . this is consistent with the newly emerging constraints from kepler , which suggests planetary radii for this class of object that require substantial gaseous envelopes ( holman et al . 2010 ; lissauer et al . further encouragment comes from the analysis of howard et al . ( 2011 ) , who estimate that planets with masses @xmath85 and periods @xmath7  days have substantial gas fractions on average , while those at lower masses are less gas - rich . evaluating our equation  ( [ gapmaker ] ) at a mean period for the howard sample @xmath86  days , yields a mass of 5@xmath17 . as the observations and analysis improve , we expect that this division will be found to rise to higher masses at longer orbital periods . although the criterion is based on the ability of a core to gravitationally retain gas ( e.g. ikoma et al . 2001 ; rafikov 2006 ) , if the ordering of mass with semi - major axis is preserved to the present day , it would represent strong evidence that there was little orbital evolution post assembly . for orbital periods @xmath4  days , the division between the two classes is less uncertain because of the possibility that a gaseous envelope could be accreted but then removed over time . the fact that the observations tend to group near the line in figure  [ pmall ] is a tantalising suggestion that evaporation could be very efficient . if true , it would suggest that most of the planets with @xmath87  days , and @xmath88 are neptune - like planets that have experienced significant erosion over the course of their lifetime . the essential point of these calculations is that many of the quantitative features of the observations can be modelled within a relatively simple in situ accumulation formalism . as such , several aspects of the more general astrophysical context need to be addressed , but the agreement with the observations argues that the basic model is correct . the mechanism(s ) by which heavy elements are concentrated in the inner parts of the nebula remain to be properly quantified . although the mobility of small bodies in protoplanetary nebulae is well established , the precise nature of the material rearrangement will depend on the properties of the background nebula . similarly , the effects of a residual gas disk on the damping of planetary eccentricities and inclinations is something to be considered . however , moeckel & armitage ( 2012 ) have shown that the statistical outcome of planetary scattering remains largely unaffected by the presence of gas . they do identify a subpopulation of low eccentricity orbits in their population , but presently the observations are still in excellent agreement with the outcomes of pure n - body calculations ( e.g. figure  [ pmall ] ) . the model also predicts a period - dependent threshold mass that separates true rocky planets from those with comparable amounts of rock and gas . this will be quantitatively tested as more low mass planets are found , but already finds some statistical support in the planet candidates uncovered by kepler ( howard et al . 2011 ; wolfgang & laughlin 2011 ; youdin 2011 ) . the authors thank phil armitage , ravit helled , kristen menou , hilke schlichting and chris thompson for discussions and comments . is supported by the nasa atp program . is supported in part by the canada research chair program and by nserc of canada . adams , f. c. & bloch , a. m. , 2009 , apj , 701 , 1381 agnor , c. b. , canup , r. m. & levison , h. f. , 1999 , icarus 142 , 219 alibert , y. et al . , 2006 , a&a , 455 , l25 artymowicz , p. & lubow , s. h. , 1996 , apj , 467 , l77 bai , x .- & stone , j. m. , 2010 , apj , 722 , 1437 batalha , m. , et al . , 2011 , apj , 729 , 27 bodenheimer , p. , hubickyj , o. & lissauer , j. j. , 2000 , icarus , 143 , 2 borucki , w. b. et al . , 2011 , apj , 736 , 19 bryden , g. , chen , x. , lin , d. n. c. , nelson , r. p. & papaloizou , j. c. b. , 1999 , apj , 514 , 344 burrows , a. et al . , 2007 , apj , 661 , 502 carter , j. a. , winn , j. n , gilliland , r. & holman , m. j. , 2009 , apj , 696 , 241 chambers , j. e. , 1999 , mnras , 304 , 793 chambers , j. e. , 2001 , icarus , 152 , 205 chambers , j. e. & wetherill , g. w. , 1998 , icarus , 136 , 304 chambers , j. e. , wetherill , g. w. & boss , a. p. , 1996 , icarus , 119 , 261 dawson , r. i. & fabrycky , d. c. , 2010 , apj , 722 , 937 ehrenreich , d. , & dsert , j .- m . , 2011 , arxiv:1103.0011 fischer , d. a. & valenti , j. , 2005 , apj , 622 , 1102 gammie , c. , 2001 , apj , 553 , 174 goldreich , p. & tremaine , s. , 1980 , apj , 241 , 425 gonzalez , g. , 1998 , a&a , 334 , 221 hansen , b. m. s. , 2009 , apj , 703 , 1131 hartman , j. d. , et al . , 2011 , apj , 728 , 138 havel , m. , guillot , t. , valencia , d. & crida , a. , 2011 , arxiv:1103.6020 hayashi , c. , 1981 , prog . suppl . , 70 , 35 holman , m. et al . , 2010 , science , 330 , 51 howard , a. w. et al . , 2010 , science , 330 , 653 howard , a. w. et al . , 2011 , apj , 726 , 73 ida , s. & lin , d. n. c. , 2008 , apj , 685 , 84 ida , s. & lin , d. n. c. , 2010 , apj , 719 , 810 ikoma , m. , emori , h. & nakazawa , k , 2001 , apj , 553 , 999 kley , w. , 1999 , mnras , 303 , 696 kley , w. & nelson , r. p. , 2012 , arxiv:1203.1184 kokubo , e. , kominami , j. & ida , s. , 2006 , apj , 642 , 1131 kokubo , e. & ida , s. , 1998 , icarus , 131 , 171 kokubo , e. & ida , s. , 2002 , apj , 581 , 666 lammer , h. , selsis , f. , ribas , i. , guinan , e. , bauer , s. j. , & weiss , w. , 2003 , apj , 598 , l121 laughlin , g. , & rozyczka , m. , 1996 , apj , 456 , 279 leger , a. , et al . , 2009 , a&a , 506 , 287 lecavelier des etangs , a. , vidal - madjar , a. , mcconnell , j. c. & hebrard , g. , 2004 , apj , 621 , 1049 lecavelier des etangs , a. , 2007 , a&a , 461 , 1185 lin , d. n. c. , bodenheimer , p. & richardson , d. c. , 1996 , nature , 380 , 606 lin , d. n. c. & papaloizou , j. , 1979 , mnras , 186 , 799 lissauer , j. j. , et al . , 2011b , nature , 470 , 53 lissauer , j. j. et al . , 2011a , apj , 197 , 8 lovis , c. et al . , 2006 , nature , 441 , 305 lovis , c. et al . , 2011 , a&a , 528 , a112 matsumura , s. , peale , s. j. & rasio , f. a. , 2010 , apj , 725 , 1995 marcy , g. , butler , p. , fischer , d. , vogt , s. , lissuaer , j. & rivera , e. , 2001 , apj , 556 , 296 mayor , m. & queloz , d. , 1995 , nature , 378 , 355 mayor , m. et al . , 2009a , a&a , 493 , 639 mayor , m. et al . , 2009b , a&a , 507 , 487 mayor , m. et al . , 2011 , arxiv:1109.2497 mcneil , d. s. & nelson , r. p. , 2010 , mnras , 401 , 1691 moeckel , n. & armitage , p. j. , 2012 , mnras , 419 , 366 mordasini , c. , alibert , y. & benz , w. , 2009 , a&a , 501 , 1139 murray , n. , b. hansen , m. holman & s. tremaine , 1998 , science , 279 , 69 murray - clay , r. a. , chiang , e. i. & murray , n. , 2009 , apj , 693 , 23 nakagawa , y. , sekiya , m. & hayashi , c. , 1986 , icarus , 67 , 375 oishi , j. s. , mac  low , m .- & menou , k. , 2007 , apj , 670 , 805 paardekooper , s. j. & mellema , g. , 2006 , a&a , 459 , l17 pollack , j. b. , hubickyj , o. , bodenheimer , p. , lissauer , j. j. , podolak , m. & greenzweig , y. , 1996 , icarus , 124 , 62 rafikov , r. , 2006 , apj , 648 , 666 rasio , f. a. & ford , e. b. , 1996 , science , 274 , 954 raymond , s. n. , barnes , r. & mandell , a. , 2008 , mnras , 384 , 663 raymond , s. n. , quinn , t. & lunine , j. i. , 2004 , icarus , 168 , 1 raymond , s. n. , quinn , t. & lunine , j. i. , 2005 , apj , 632 , 670 rice , w. k. m. & armitage , p. j. , 2009 , mnras , 396 , 2228 rogers , l. a. & seager , s. , 2010 , apj , 712 , 974 santos , n. c. , israelian , g. & mayor , m. , 2001 , a&a , 373 , 1019 shakura , n. i. & sunyaev , r. a. , 1973 , a&a , 24 , 337 taniguchi , t. & ikoma , m. , 2007 , apj , 667 , 55 terquem , c. & papaloizou , j. c. b. , 2007 , apj , 654 , 1110 tian , f. , toon , o. b. , pavlov , a. a. & de sterck , h. , 2005 , apj , 621 , 1049 udry , s. & santos , n. c. , 2007 , ara&a , 45 , 397 vogt , s. et al . , 2009 , apj , 708 , 1366 ward , w. r. , 1997 , icarus , 126 , 261 weidenschilling , s. j. , 1977a , ap&ss , 51 , 153 weidenschilling , s. j. , 1977b , mnras , 180 , 57 winn , j.n . , et al . , 2011a , aj , 141 , 63 winn , j. n. et al . , 2011b , apj , 737 , l18 wolfgang , a. & laughlin , g. , 2011 , arxiv:1108.5842 yelle , r. v. , 2004 , icarus , 170 , 167 youdin , a. & shu , f. , 2002 , apj , 580 , 494 youdin , a. , 2011 , apj , 742 , 38 for the kinds of disk masses considered in this paper , the oligarchic formalism of kokubo & ida ( 1998 ) can produce bodies with masses @xmath89 or more . with masses this large , one might wonder whether the long - range gravitational perturbations between neighbouring annuli become important sooner , and perhaps the assumption of initially independent annuli no longer holds . to examine this , we have repeated our calculations but have restricted the maximum size of an oligarch to @xmath25 . the bodies are still spaced by annuli with width ten times the hill radius , and so the overall surface density distribution is preserved . in table [ democheck ] we show the statistics for 5 realisations of this model , for a disk with power law @xmath28 and total mass @xmath26 . these can be compared with the results of the runs run@xmath9012 in table  [ simparam ] . we see that there is little difference in the final statistical measures . the only noticeable difference is that the runs for this section lose a few percent of material in scattered and ejected bodies , which leads to a slightly higher probability that one of the surviving bodies will have a radial velocity signature @xmath32 . as a result we quote two sets of statistics for some of the runs in table  [ democheck ] . the first line represents the direct comparison to other simulations . in those cases where we show a second set of values for a simulation , it means that one of the surviving planets would not have been detected with a 1 m / s radial velocity threshold . the second line shows the statistics calculated if we include those planets that would likely have been missed . the equations collected here offer a brief description of each of the statistical measures used in the text . further information can be found in chambers ( 2001 ) and references therein . the angular momentum deficit is defined as @xmath91 } { \sum_j m_j \sqrt{a_j}}\ ] ] where @xmath92 , @xmath93 , @xmath94 and @xmath95 are the mass , semi - major axis , eccentricity and inclination of the j - th particle . it is a measure of how much the orbits in the system deviate from the case of circular , coplanar orbits , which represent the state of maximal angular momentum per unit energy . the mass - weighted semi - major axis is @xmath96 and is a measure of where the mass is eventually accumulated . the mass concentration statistic is given as the maximum value of the function @xmath97 ^ 2 } \right)\ ] ] as the semi - major axis @xmath48 is varied . the orbital spacing statistic is @xmath98 where @xmath99 is the mean mass of the planets , and @xmath100(@xmath101 ) is the semi - major axis of the shortest ( longest ) period of the n planets in the system . this quantity is often considered a measure of the long - term dynamical stability of the system , as it bears a qualitative similarity to the average spacing of the planets measured in terms of their hill radii , @xmath102 . dynamical simulations can determine the threshold value of these parameters for long term stability ( e.g. chambers , wetherill & boss 1996 ) , but require some choices with regards to the masses of the various components . thus , we have repeated this calculation for the specific case of interest here . we consider three - planet systems , with the middle planet at 0.35  au . the planet masses are given by equation  ( [ gasscale ] ) , and two other planets in the system were placed at equal separation @xmath103 on either side . the results suggest that stability is achieved for @xmath104  au , corresponding to @xmath105 . comparison to figure  [ param6 ] shows that nearly all the observed systems satisfy this criterion . the two that do not are gliese  876 ( which corresponds to a pair of jupiter mass planets orbiting an m - dwarf in resonance ) and kepler-11 ( which satisfies this criterion if the mass of the outermost planet , kepler-11 g , lies somewhat below the currently quoted upper limit ) . in the early stages of star formation , the mass fraction of material in the disk may be substantial and the disk may exhibit a variety of large scale gravitational instabilities . this can also lead to significant radial transport of gas and dust ( e.g. laughlin & rozyczka 1996 ; gammie 2001 ) and may result ina somewhat different radial profile than those considered above . thus , we will also consider a model in which the solid distribution is deposited in situ , but from a gas disk which follows a disk profile dictated by setting the toomre parameter to @xmath65 , which we take as the marginally stable end state of a gravitationally unstable disk . the temperature profile is also somewhat different because the disk is heated by the dissipation inher ent in the mass transport . based on the results of rice & armitage ( 2009 ) , we assume a constant temperature @xmath106  k interior to 1 au , and derive @xmath107 where we have used a conversion factor of @xmath108 between gas and solid material . this is a small fraction of the gas metallicity , but captures the likely inefficiency of sedimentation at these temperatures .    the assembly of embryos , drawn from this distribution , into planets proceeds in the same way as in the other simulations , except for the fact that there is considerably more mass available . the final planets produced have masses that range up to @xmath1 . this particular scenario may never be realised in practice , but is worth calculating if only to demonstrate that the gap in the observed masses is the result of the available reservoir of material , not of the physics of assembly itself . it is only when the masses begin to approach the ejection / accretion division shown in figure  [ pm ] that the process of collisional accumulation becomes self - limiting . lccccc|cccccc run@xmath109 & 230 & 1.75 & @xmath110 & @xmath110 & 57 & 2 & 99.3 & 36.5 & 0.073 & 8.5 & 0.168 + run@xmath109 & 230 & 1.75 & @xmath110 & @xmath110 & 57 & 2 & 127.5 & 31.1 & 0.012 & 10.3 & 0.291 + run@xmath109 & 230 & 1.75 & @xmath110 & @xmath110 & 57 & 2 & 97.9 & 29.4 & 0.027 & 13.5 & 0.161 + run@xmath111 & 25 & 0 & @xmath110 & @xmath110 & 79 & 1 & 16.3 & @xmath110 & 0.034 & @xmath112 & 0.545 + run@xmath113 & 25 & 0 & @xmath110 & @xmath110 & 79 & 1 & 12.5 & @xmath110 & 0.1011 & @xmath112 & 0.485 + run@xmath114 & 25 & 0 & @xmath110 & @xmath110 & 79 & 1 & 9.9 & @xmath115 & 0.0005 & @xmath112 & 0.906 + run@xmath116 & 50 & 0 & @xmath110 & @xmath110 & 57 & 2 & 46.1 & 38.6 & 0.0136 & 22.8 & 0.740 + run@xmath117 & 50 & 0 & @xmath110 & @xmath110 & 57 & 3 & 27.1 & 21.2 & 0.0084 & 17.3 & 0.738 + run@xmath118 & 50 & 0 & @xmath110 & @xmath110 & 57 & 3 & 27.1 & 17.5 & 0.0031 & 45.0 & 0.705 + run@xmath119 & 50 & 0 & @xmath110 & @xmath110 & 57 & 3 & 42.7 & 26.3 & 0.0081 & 15.8 & 0.717 + run@xmath120 & 50 & 0 & @xmath110 & @xmath110 & 57 & 3 & 32.1 & 15.6 & 0.0013 & 59.9 & 0.734 + run@xmath121 & 50 & 0 & @xmath110 & @xmath110 & 57 & 3 & 28.0 & 11.4 & 0.0006 & 87.1 & 0.750 + run@xmath122 & 50 & 0 & @xmath110 & @xmath110 & 57 & 2 & 28.0 & 21.5 & 0.0014 & 43.3 & 0.767 + run@xmath123 & 50 & 0 & @xmath110 & @xmath110 & 57 & 1 & 40.5 & @xmath110 & 0.0194 & @xmath112 & 0.559 + run@xmath124 & 50 & 0 & @xmath110 & @xmath110 & 57 & 3 & 28.0 & 21.9 & 0.0148 & 16.4 & 0.755 + run@xmath125 & 100 & 0 & @xmath110 & @xmath110 & 41 & 2 & 61.3 & 23.4 & 0.0125 & 46.4 & 0.664 + run@xmath126 & 100 & 0 & @xmath110 & @xmath110 & 41 & 3 & 51.3 & 19.3 & 0.0020 & 37.2 & 0.628 + run@xmath127 & 100 & 0 & @xmath110 & @xmath110 & 41 & 4 & 50.3 & 18.5 & 0.0042 & 23.9 & 0.614 + run@xmath128 & 25 & 1.5 & @xmath110 & @xmath110 & 34 & 5 & 6.3 & 24.1 & 0.0114 & 6.2 & 0.401 + run@xmath129 & 25 & 1.5 & @xmath110 & @xmath110 & 34 & 3 & 10.2 & 36.6 & 0.0189 & 11.7 & 0.234 + run@xmath130 & 25 & 1.5 & @xmath110 & @xmath110 & 34 & 4 & 9.3 & 29.4 & 0.0136 & 6.9 & 0.380 + run@xmath90 & 50 & 1.5 & @xmath110 & @xmath110 & 25 & 3 & 17.3 & 32.1 & 0.0102 & 6.2 & 0.383 + run@xmath131 & 50 & 1.5 & @xmath110 & @xmath110 & 25 & 2 & 16.5 & 48.5 & 0.0417 & 11.3 & 0.156 + run@xmath132 & 50 & 1.5 & @xmath110 & @xmath110 & 25 & 2 & 23.3 & 48.7 & 0.0484 & 12.8 & 0.199 + run@xmath133 & 50 & 1.5 & @xmath110 & @xmath110 & 25 & 3 & 22.6 & 33.7 & 0.0033 & 5.1 & 0.513 + run@xmath134 & 50 & 1.5 & @xmath110 & @xmath110 & 25 & 4 & 19.9 & 24.6 & 0.0100 & 5.9 & 0.474 + run@xmath135 & 50 & 1.5 & @xmath110 & @xmath110 & 25 & 3 & 17.8 & 33.3 & 0.0059 & 8.4 & 0.304 + run@xmath136 & 50 & 1.5 & @xmath110 & @xmath110 & 25 & 4 & 11.7 & 30.4 & 0.0569 & 10.0 & 0.208 + run@xmath137 & 50 & 1.5 & @xmath110 & @xmath110 & 25 & 3 & 26.4 & 31.9 & 0.0395 & 7.4 & 0.341 + run@xmath138 & 50 & 1.5 & @xmath110 & @xmath110 & 25 & 4 & 15.9 & 34.0 & 0.0087 & 5.2 & 0.428 + run@xmath139 & 50 & 1.5 & @xmath110 & @xmath110 & 25 & 4 & 12.8 & 24.6 & 0.0094 & 8.6 & 0.286 + run@xmath140 & 50 & 1.5 & @xmath110 & @xmath110 & 25 & 3 & 23.3 & 33.9 & 0.0121 & 8.0 & 0.316 + run@xmath141 & 50 & 1.5 & @xmath110 & @xmath110 & 25 & 4 & 14.5 & 24.2 & 0.0076 & 6.4 & 0.384 + run@xmath142 & 100 & 1.5 & @xmath110 & @xmath110 & 18 & 5 & 28.3 & 16.8 & 0.0034 & 6.5 & 0.426 + run@xmath143 & 100 & 1.5 & @xmath110 & @xmath110 & 18 & 2 & 53.3 & 54.5 & 0.0067 & 13.6 & 0.588 + run@xmath144 & 100 & 1.5 & @xmath110 & @xmath110 & 18 & 4 & 37.3 & 20.5 & 0.0159 & 5.9 & 0.442 + run@xmath145 & 50 & 2.5 & @xmath110 & @xmath110 & 20 & 3 & 24.0 & 26.5 & 0.0093 & 13.8 & 0.132 + run@xmath146 & 50 & 2.5 & @xmath110 & @xmath110 & 20 & 3 & 19.4 & 30.9 & 0.0063 & 9.9 & 0.155 + run@xmath147 & 50 & 2.5 & @xmath110 & @xmath110 & 20 & 3 & 22.6 & 28.0 & 0.0125 & 16.1 & 0.123 + run@xmath148 & 50 & 2.5 & @xmath110 & @xmath110 & 20 & 4 & 19.8 & 19.4 & 0.0057 & 12.2 & 0.116 + run@xmath149 & 50 & 2.5 & @xmath110 & @xmath110 & 20 & 4 & 19.8 & 18.7 & 0.0040 & 14.8 & 0.102 + run@xmath150 & 50 & 2.5 & @xmath110 & @xmath110 & 20 & 5 & 16.6 & 17.3 & 0.0082 & 10.7 & 0.126 + run@xmath151 & 50 & 2.5 & @xmath110 & @xmath110 & 20 & 4 & 14.5 & 20.7 & 0.0023 & 13.8 & 0.110 + run@xmath152 & 50 & 2.5 & @xmath110 & @xmath110 & 20 & 4 & 14.5 & 21.2 & 0.0068 & 10.2 & 0.131 + run@xmath153 & 50 & 2.5 & @xmath110 & @xmath110 & 20 & 4 & 13.4 & 21.3 & 0.0039 & 11.8 & 0.133 + run@xmath154 & 100 & 2.5 & @xmath110 & @xmath110 & 14 & 4 & 40.5 & 21.0 & 0.0132 & 6.9 & 0.198 + run@xmath155 & 100 & 2.5 & @xmath110 & @xmath110 & 14 & 3 & 60.3 & 25.2 & 0.0099 & 9.3 & 0.165 + run@xmath156 & 100 & 2.5 & @xmath110 & @xmath110 & 14 & 4 & 49.7 & 18.6 & 0.0133 & 9.6 & 0.164 +   + run@xmath157 & 50 & 1.5 & 2300 & 1.5 & 25 & 4 & 62.0 & 19.6 & 0.023 & 11.1 & 0.311 + run@xmath158 & 50 & 1.5 & 2300 & 1.5 & 25 & 4 & 60.7 & 20.6 & 0.011 & 10.9 & 0.267 + run@xmath159 & 50 & 1.5 & 2300 & 1.5 & 25 & 4 & 101.4 & 17.6 & 0.0047 & 9.6 & 0.354 + run@xmath160 & 50 & 1.5 & 2300 & 1.5 & 25 & 3 & 40.8 & 27.9 & 0.0230 & 6.8 & 0.315 + run@xmath161 & 50 & 1.5 & 2300 & 1.5 & 25 & 3 & 67.3 & 23.5 & 0.0091 & 12.8 & 0.230 + run@xmath162 & 50 & 1.5 & 2300 & 1.5 & 25 & 3 & 52.0 & 27.6 & 0.0146 & 6.8 & 0.301 + run@xmath163 & 50 & 1.5 & 2300 & 1.5 & 25 & 3 & 42.6 & 27.9 & 0.0145 & 8.1 & 0.259 + run@xmath164 & 50 & 1.5 & 2300 & 1.5 & 25 & 2 & 100.0 & 19.9 & 0.0059 & 34.3 & 0.410 + run@xmath165 & 40 & 1.5 & 767 & 1.0 & 30 & 4 & 22.1 & 23.9 & 0.0025 & 9.3 & 0.217 + run@xmath166 & 40 & 1.5 & 767 & 1.0 & 30 & 3 & 15.7 & 32.3 & 0.011 & 12.6 & 0.150 + run@xmath167 & 40 & 1.5 & 767 & 1.0 & 30 & 2 & 30.6 & 44.0 & 0.0025 & 11.3 & 0.196 + run@xmath168 & 40 & 1.5 & 767 & 1.0 & 30 & 4 & 28.9 & 24.0 & 0.0121 & 6.2 & 0.341 + run@xmath169 & 40 & 1.5 & 767 & 1.0 & 30 & 5 & 26.4 & 18.5 & 0.0891 & 7.4 & 0.313 + run@xmath170 & 40 & 1.5 & 767 & 1.0 & 30 & 5 & 13.9 & 20.2 & 0.0057 & 7.3 & 0.280 + run@xmath171 & 40 & 1.5 & 767 & 1.0 & 30 & 4 & 19.2 & 22.8 & 0.0097 & 6.7 & 0.317 + run@xmath172 & 40 & 1.5 & 767 & 1.0 & 30 & 3 & 18.7 & 35.5 & 0.0225 & 6.9 & 0.242 + run@xmath173 & 40 & 1.5 & 767 & 1.0 & 30 & 5 & 23.3 & 18.3 & 0.0063 & 5.6 & 0.350 + run@xmath174 & 40 & 1.5 & 767 & 1.0 & 30 & 3 & 32.9 & 34.1 & 0.0277 & 8.0 & 0.241 + run@xmath175 & 40 & 1.5 & 767 & 1.0 & 30 & 11 & 7.2 & 15.4 & 0.0007 & 6.1 & 0.442 + run@xmath176 & 40 & 1.5 & 767 & 1.0 & 30 & 1 & 43.2 & @xmath110 & 0.163 & @xmath112 & 0.106 + run@xmath177 & 50 & 1.0 & 767 & 1.5 & 30 & 4 & 28.2 & 20.9 & 0.0083 & 7.9 & 0.511 + run@xmath178 & 50 & 1.0 & 767 & 1.5 & 30 & 3 & 37.8 & 28.8 & 0.0011 & 10.6 & 0.344 + run@xmath179 & 50 & 1.0 & 767 & 1.5 & 30 & 3 & 51.1 & 27.9 & 0.0243 & 8.3 & 0.549 + run@xmath180 & 50 & 1.0 & 767 & 1.5 & 30 & 4 & 39.8 & 21.2 & 0.0042 & 8.6 & 0.511 + run@xmath181 & 50 & 1.0 & 767 & 1.5 & 30 & 3 & 55.4 & 25.9 & 0.0151 & 13.8 & 0.491 + run@xmath182 & 80 & 1.0 & 767 & 1.5 & 22 & 2 & 45.2 & 28.4 & 0.101 & 21.9 & 0.368 + run@xmath183 & 80 & 1.0 & 767 & 1.5 & 22 & 4 & 57.3 & 18.1 & 0.0078 & 9.1 & 0.580 + run@xmath184 & 80 & 1.0 & 767 & 1.5 & 22 & 3 & 49.5 & 25.7 & 0.0290 & 10.2 & 0.400 + run@xmath185 & 80 & 1.0 & 767 & 1.5 & 22 & 3 & 50.9 & 26.6 & 0.0352 & 10.6 & 0.339 + run@xmath186 & 80 & 1.0 & 767 & 1.5 & 22 & 2 & 57.3 & 31.2 & 0.0041 & 14.3 & 0.433 + run@xmath187 & 80 & 1.0 & 767 & 1.5 & 22 & 3 & 48.9 & 25.8 & 0.0086 & 7.4 & 0.500 + run@xmath188 & 50 & 1.5 & @xmath110 & @xmath110 & 25 & 3 & 11.8 & 33.3 & 0.0200 & 9.4 & 0.215 + run@xmath189 & 50 & 1.5 & @xmath110 & @xmath110 & 25 & 3 & 10.6 & 30.9 & 0.0305 & 10.4 & 0.194 + run@xmath190 & 50 & 1.5 & @xmath110 & @xmath110 & 25 & 3 & 14.5 & 30.8 & 0.0173 & 8.2 & 0.264 +    lccccc|cccccc demo@xmath191 & 50 & 1.5 & @xmath110 & @xmath110 & 52 & 3 & 23.7 & 31.1 & 0.0254 & 6.5 & 0.379 + demo@xmath192 & 50 & 1.5 & @xmath110 & @xmath110 & 52 & 3 & 15.0 & 28.4 & 0.0214 & 13.1 & 0.291 + & 50 & 1.5 & @xmath110 & @xmath110 & 52 & 4 & 15.0 & 25.3 & 0.0199 & 9.7 & 0.356 + demo@xmath193 & 50 & 1.5 & @xmath110 & @xmath110 & 52 & 2 & 17.0 & 45.3 & 0.0467 & 13.2 & 0.169 + & 50 & 1.5 & @xmath110 & @xmath110 & 52 & 3 & 17.0 & 33.4 & 0.0471 & 7.6 & 0.272 + demo@xmath194 & 50 & 1.5 & @xmath110 & @xmath110 & 52 & 2 & 22.0 & 45.0 & 0.0495 & 13.9 & 0.174 + & 50 & 1.5 & @xmath110 & @xmath110 & 52 & 3 & 22.0 & 36.0 & 0.0498 & 5.0 & 0.422 + demo@xmath195 & 50 & 1.5 & @xmath110 & @xmath110 & 52 & 3 & 21.0 & 34.2 & 0.0218 & 8.4 & 0.291
Abstract:

Answer: we demonstrate that the observed distribution of ` hot neptune'/`super - earth ' systems is well reproduced by a model in which planet assembly occurs in situ , with no significant migration post - assembly . this is achieved only if the amount of mass in rocky material is @xmath0@xmath1 interior to 1  au . such a reservoir of material implies that significant radial migration of solid material takes place , and that it occur before the stage of final planet assembly . the model not only reproduces the general distribution of mass versus period , but also the detailed statistics of multiple planet systems in the sample . we furthermore demonstrate that cores of this size are also likely to meet the criterion to gravitationally capture gas from the nebula , although accretion is rapidly limited by the opening of gaps in the gas disk . if the mass growth is limited by this tidal truncation , then the scenario sketched here naturally produces neptune - mass objects with substantial components of both rock and gas , as is observed . the quantitative expectations of this scenario are that most planets in the ` hot neptune / super - earth ' class inhabit multiple - planet systems , with characteristic orbital spacings . the model also provides a natural division into gas - rich ( hot neptune ) and gas - poor ( super - earth ) classes at fixed period . the dividing mass ranges from @xmath2 at 10  day orbital periods to @xmath3 at 100  day orbital periods . for orbital periods @xmath4  days , the division is less clear because a gas atmosphere may be significantly eroded by stellar radiation .