Researcher profile

Philippe Thebault

Philippe Thebault contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
20works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2020arXiv

Debris discs in binaries: morphology and photometric signatures

We aim to see whether debris belts evolving in between two stars may be impacted by the presence of the companion and whether this leaves any detectable signature that could be observed with current or future instruments. We consider a circumprimary parent body (PB) planetesimal belt that is placed just inside the stability limit between the 2 stars and use the DyCoSS code to follow the coupled dynamical and collisional evolution of the dust produced by this PB belt. We explore several free parameters such as the belt's mass or the binary's mass ratio and orbit. We use the GraTeR package to produce 2-D luminosity maps and system-integrated SEDs. We confirm a preliminary result obtained by earlier DyCoSS studies, which is that the coupled effect of collisional activity, binary perturbations and stellar radiation pressure maintains a halo of small grains in the dynamically unstable region between the 2 stars. In addition, several spatial structures are identified, notably a single spiral arm stretching all the way from the PB belt to the companion star. We also identify a fainter and more compact disc around the secondary star, which is non-native and feeds off small grains from the unstable halo. Both the halo, spiral arm and secondary disc should be detectable on resolved images by instruments with capacities similar to SPHERE. The system as a whole is depleted in small grains when compared to a companion-free case. This depletion leaves an imprint on the system's integrated SED, which appears colder than for the same parent body belt around a single star. This new finding could explain why the SED-derived location $r_{disc}$ of some unresolved discs-in-binaries places their primary belt in the dynamically "forbidden" region between the 2 stars: this apparent paradox could indeed be due to overestimating $r_{disc}$ when using empirical prescriptions valid for a single star case

preprint2020arXiv

Planets in Binaries: Formation and Dynamical Evolution

Binary systems are very common among field stars. While this relatively small number of planets in binaries is probably partly due to strong observational biases, there is, however, statistical evidence that planets are indeed less frequent in binaries with separations smaller than 100 au, strongly suggesting that the presence of a close in companion star has an adverse effect on planet formation. It is indeed possible for the gravitational pull of the second star to affect all the different stages of planet formation, from proto-planetary disk formation to dust accumulation into planetesimals, to the accretion of these planetesimals into large planetary embryos and, eventually, the final growth of these embryos into planets. For the crucial planetesimal accretion phase, the complex coupling between dynamical perturbations from the binary and friction due to gas in the protoplanetary disk suggests that planetesimal accretion might be hampered due to increased, accretion hostile impact velocities. Likewise, the interplay between the binary secular perturbations and mean motion resonances lead to unstable regions, where not only planet formation is inhibited, but where a massive body would be ejected from the system on a hyperbolic orbit. The amplitude of these two main effects is different for S and P type planets, so that a comparison between the two populations might outline the influence of the companion star on the planet formation process. Unfortunately, at present the two populations (circumstellar or circumbinary) are not known equally well and different biases and uncertainties prevent a quantitative comparison. We also highlight the long term dynamical evolution of both S and P type systems and focus on how these different evolutions influence the final architecture of planetary systems in binaries.

preprint2016arXiv

Dust production in debris discs: constraints on the smallest grains

The surface energy constraint puts a limit on the smallest fragment $s_{surf}$ that can be produced after a collision. Based on analytical considerations, this mechanism has been recently identified as been potentially able to prevent the production of small dust grains in debris discs and cut off their size distribution at sizes larger than the blow-out size. We numerically investigate the importance of this effect to find under which conditions it can leave a signature in the small-size end of a disc's particle size distribution (PSD). An important part of this work is to map out, in a disc at steady-state, what is the most likely collisional origin for micron-sized grains, in terms of the sizes of their collisional progenitors. We implement, for the first time, the surface energy constraint into a collisional evolution code. We consider a debris disc extending from 50 to 100AU and 2 different stellar types. We also consider two levels of stirring in the disc: dynamically "hot" (e=0.075) and "cold" (e=0.01). For all cases, we derive $s_{surf}$ maps as a function of target and projectile sizes, $s_t$ and $s_p$, and compare them to equivalent maps for the dust-production rate. We then compute disc-integrated PSDs and estimate the imprint of the surface energy constraint. We find that the ($s_p$,$s_t$) regions of high $s_{surf}$ values do not coincide with those of high dust production rate. As a consequence, the surface energy constraint has generally a weak effect on the system's PSD. The maximum $s_{surf}$-induced depletion of $μ$m-sized grains is $\sim 30$% and is obtained for a sun-like star and a dynamically hot case. For the e=0.01 cases, the surface energy effect is negligible compared to the massive small grain depletion induced by another mechanism: the natural imbalance between dust production and destruction rates in low-stirring discs identified by Thebault&Wu(2008).

preprint2014arXiv

Forming Circumbinary Planets: N-body Simulations of Kepler-34

Observations of circumbinary planets orbiting very close to the central stars have shown that planet formation may occur in a very hostile environment, where the gravitational pull from the binary should be very strong on the primordial protoplanetary disk. Elevated impact velocities and orbit crossings from eccentricity oscillations are the primary contributors towards high energy, potentially destructive collisions that inhibit the growth of aspiring planets. In this work, we conduct high resolution, inter-particle gravity enabled N-body simulations to investigate the feasibility of planetesimal growth in the Kepler-34 system. We improve upon previous work by including planetesimal disk self-gravity and an extensive collision model to accurately handle inter-planetesimal interactions. We find that super-catastrophic erosion events are the dominant mechanism up to and including the orbital radius of Kepler-34(AB)b, making in-situ growth unlikely. It is more plausible that Kepler-34(AB)b migrated from a region beyond 1.5 AU. Based on the conclusions that we have made for Kepler-34 it seems likely that all of the currently known circumbinary planets have also migrated significantly from their formation location with the possible exception of Kepler-47(AB)c.

preprint2014arXiv

Signatures of massive collisions in debris discs

Violent stochastic collisional events have been invoked as a possible explanation for some debris discs displaying pronounced asymmetries or having a great luminosity excess. So far, no thorough modelling of the consequences of such events has been carried out, mainly because of the extreme numerical challenge of coupling the dynamical and collisional evolution of dust. We perform the first fully self-consistent modelling of the aftermath of massive breakups in debris discs. We follow the collisional and dynamical evolution of dust released after the breakup of a Ceres-sized body at 6 AU from its central star. We investigate the duration, magnitude and spatial structure of the signature left by such a violent event, as well as its observational detectability. We use the recently developed LIDT-DD code (Kral et al., 2013), which handles the coupled collisional and dynamical evolution of debris discs. The main focus is placed on the complex interplay between destructive collisions, Keplerian dynamics and radiation pressure forces. We use the GRaTer package to estimate the system's luminosity at different wavelengths. The breakup of a Ceres-sized body at 6 AU creates an asymmetric dust disc that is homogenized, by the coupled action of collisions and dynamics, on a timescale of a few $10^5$ years. The luminosity excess in the breakup's aftermath should be detectable by mid-IR photometry, from a 30 pc distance, over a period of $\sim 10^6$ years that exceeds the duration of the asymmetric phase of the disc (a few $10^5$ years). As for the asymmetric structures, we derive synthetic images for the SPHERE/VLT and MIRI/JWST instruments, showing that they should be clearly visible and resolved from a 10 pc distance. Images at 1.6$μ$m (marginally), 11.4 and 15.5$μ$m would show the inner disc structures while 23$μ$m images would display the outer disc asymmetries.

preprint2013arXiv

An independent determination of Fomalhaut b's orbit and the dynamical effects on the outer dust belt

The nearby star Fomalhaut harbours a cold, moderately eccentric dust belt with a sharp inner edge near 133 au. A low-mass, common proper motion companion (Fom b), was discovered near the inner edge and was identified as a planet candidate that could account for the belt morphology. However, the most recent orbit determination based on four epochs of astrometry over eight years reveals a highly eccentric orbit that appears to cross the belt in the sky plane projection. We perform here a full orbital determination based on the available astrometric data to independently validate the orbit estimates previously presented. Adopting our values for the orbital elements and their associated uncertainties, we then study the dynamical interaction between the planet and the dust ring, to check whether the proposed disk sculpting scenario by Fom b is plausible. We used a dedicated MCMC code to derive the statistical distributions of the orbital elements of Fom b. Then we used symplectic N-body integration to investigate the dynamics of the dust belt, as perturbed by a single planet. Different attempts were made assuming different masses for Fom b. We also performed a semi-analytical study to explain our results. Our results are in good agreement with others regarding the orbit of Fom b. We find that the orbit is highly eccentric, is close to apsidally aligned with the belt, and has a moderate mutual inclination relative to the belt plane of. If coplanar, this orbit crosses the disk. Our dynamical study then reveals that the observed planet could sculpt a transient belt configuration with a similar eccentricity to what is observed, but it would not be simultaneously apsidally aligned with the planet. This transient configuration only occurs a short time after the planet is placed on such an orbit (assuming an initially circular disk), a time that is inversely proportional to the planet's mass, and that is in any case much less than the 440 Myr age of the star. We constrain how long the observed dust belt could have survived with Fom b on its current orbit, as a function of its possible mass. This analysis leads us to conclude that Fom b is likely to have low mass, that it is unlikely to be responsible for the sculpting of the belt, and that it supports the hypothesis of a more massive, less eccentric planet companion Fom c.

preprint2013arXiv

Grain Size segregation in debris discs

In most debris discs, dust grain dynamics is strongly affected by stellar radiation pressure. As this mechanism is size-dependent, we expect dust grains to be spatially segregated according to their sizes. However, because of the complex interplay between radiation pressure, collisions and dynamical perturbations, this spatial segregation of the particle size distribution (PSD) has proven difficult to investigate with numerical models. We propose to explore this issue using a new-generation code that can handle some of the coupling between dynamical and collisional effects. We investigate how PSDs behave in both unperturbed discs &#34;at rest&#34; and in discs pertubed by planetary objects. We use the DyCoSS code of Thebault(2012) to investigate the coupled effect of collisions, radiation pressure and dynamical perturbations in systems having reached a steady state. We consider 2 setups: a narrow ring perturbed by an exterior planet, and an extended disc into which a planet is embedded. For both setups we consider an additional unperturbed case with no planet. We also investigate how possible spatial size segregation affect disc images at different wavelengths. We find that PSDs are always strongly spatially segregated. The only case for which they follow a standard dn/dr = C.r**(-3.5) law is for an unperturbed narrow ring, but only within the parent body ring itself. For all other configurations, the PSD can strongly depart from such power laws and have strong spatial gradients. As an example, the geometrical cross section of the disc is rarely dominated by the smallest grains on bound orbits, as it is expected to be in standard PSDs in s**q with q<-3. Although the exact profiles and spatial variations of PSDs are a complex function of the considered set-up, we are however able to derive some robust results that should be useful for image-or-SED-fitting models of observed discs.

preprint2013arXiv

Influence of the circumbinary disk gravity on planetesimal accumulation in the Kepler 16 system

Recent observations from NASA&#39;s Kepler mission detected the first planets in circumbinary orbits. The question we try to answer is where these planets formed in the circumbinary disk and how far inside they migrated to reach their present location. We investigate the first and more delicate phase of planet formation when planetesimals accumulate to form planetary embryos. We use the hydrodynamical code FARGO to study the evolution of the disk and of a test population of planetesimals embedded in it. With this hybrid hydrodynamical--N--body code we can properly account for the gas drag force on the planetesimals and for the gravitational force of the disk on them. The numerical simulations show that the gravity of the eccentric disk on the planetesimal swarm excites their eccentricities to values much larger than those induced by the binary perturbations only within 10 AU from the stars. Moreover, the disk gravity prevents a full alignment of the planetesimal pericenters. Both these effects lead to large impact velocities, beyond the critical value for erosion. Planetesimals accumulation in circumbinary disks appears to be prevented close to the stellar pair by the gravitational perturbations of the circumbinary disk. The observed planets possibly formed in the outer regions of the disk and then migrated inside by tidal interaction with the disk.

preprint2012arXiv

A Likely Close-In Low-Mass Stellar Companion to the Transitional Disk Star HD 142527

With the uniquely high contrast within 0.1&#34; (Δmag(L&#39;) = 5-6.5 magnitudes) available using Sparse Aperture Masking (SAM) with NACO at VLT, we detected asymmetry in the flux from the Herbig Fe star HD 142527 with a barycenter emission situated at a projected separation of 88+/-5 mas (12.8+/-1.5 AU at 145 pc) and flux ratios in H, K, and L&#39; of 0.016+/-0.007, 0.012+/-0.008, 0.0086+/-0.0011 respectively (3-σerrors) relative to the primary star and disk. After extensive closure-phase modeling, we interpret this detection as a close-in, low-mass stellar companion with an estimated mass of ~0.1-0.4 M_Sun. HD 142527 has a complex disk structure, with an inner gap imaged in both the near and mid-IR as well as a spiral feature in the outer disk in the near-IR. This newly detected low-mass stellar companion may provide a critical explanation of the observed disk structure.

preprint2012arXiv

How not to build Tatooine: the difficulty of in situ formation of circumbinary planets Kepler 16b, Kepler 34b and Kepler 35b

We study planetesimal evolution in circumbinary disks, focusing on the three systems Kepler 16, 34 and 35 where planets have been discovered recently. We show that for circumbinary planetesimals, in addition to secular forcing, eccentricities evolve on a dynamical timescale, which leads to orbital crossings even in the presence of gas drag. This makes the current locations of the circumbinary Kepler planets hostile to planetesimal accretion. We then present results from simulations including planetesimal formation and dust accretion, and show that even in the most favourable case of 100% efficient dust accretion, in situ growth starting from planetesimals smaller than ~10 km is difficult for Kepler 16b, Kepler 34b and Kepler 35b. These planets were likely assembled further out in the disk, and migrated inward to their current location.

preprint2012arXiv

Planet Signatures in Collisionally Active Debris Discs: scattered light images

Planet perturbations are often invoked as a potential explanation for many spatial structures that have been imaged in debris discs. So far this issue has been mostly investigated with collisionless N-body numerical models. We numerically investigate how the coupled effect of collisions and radiation pressure can affect the formation and survival of radial and azimutal structures in a disc perturbed by a planet. We consider two set-ups: a planet embedded within an extended disc and a planet exterior to an inner debris ring. We use the DyCoSS code of Thebault(2012) and derive synthetic images of the system in scattered light. The planet&#39;s mass and orbit, as well as the disc&#39;s collisional activity are explored as free parameters. We find that collisions always significantly damp planet-induced structures. For the case of an embedded planet, the planet&#39;s signature, mostly a density gap around its radial position, should remain detectable in head-on images if M_planet > M_Saturn. If the system is seen edge-on, however, inferring the presence of the planet is much more difficult, although some planet-induced signatures might be observable under favourable conditions. For the inner-ring/external-planet case, planetary perturbations cannot prevent collision-produced small fragments from populating the regions beyond the ring: The radial luminosity profile exterior to the ring is close to the one it should have in the absence of the planet. However, a Jovian planet on a circular orbit leaves precessing azimutal structures that can be used to indirectly infer its presence. For a planet on an eccentric orbit, the ring is elliptic and the pericentre glow effect is visible despite of collisions and radiation pressure, but detecting such features in real discs is not an unambiguous indicator of the presence of an outer planet.

preprint2012arXiv

Scattering of small bodies by planets: a potential origin for exozodiacal dust ?

High levels of exozodiacal dust are observed around a growing number of main sequence stars. The origin of such dust is not clear, given that it has a short lifetime against both collisions and radiative forces. Even a collisional cascade with km-sized parent bodies, as suggested to explain outer debris discs, cannot survive sufficiently long. In this work we investigate whether the observed exozodiacal dust could originate from an outer planetesimal belt. We investigate the scattering processes in stable planetary systems in order to determine whether sufficient material could be scattered inwards in order to retain the exozodiacal dust at its currently observed levels. We use N-body simulations to investigate the efficiency of this scattering and its dependence on the architecture of the planetary system. The results of these simulations can be used to assess the ability of hypothetical chains of planets to produce exozodi in observed systems. We find that for older (>100Myr) stars with exozodiacal dust, a massive, large radii (>20AU) outer belt and a chain of tightly packed, low-mass planets would be required in order to retain the dust at its currently observed levels. This brings into question how many, if any, real systems possess such a contrived architecture and are therefore capable of scattering at sufficiently high rates to retain exozodi dust on long timescales.

preprint2011arXiv

A new code to study structures in collisionally active, perturbed debris discs. Application to binaries

Debris discs are traditionally studied using two distinct types of numerical models: statistical particle-in-a-box codes to study their collisional and size distribution evolution, and dynamical N-body models to study their spatial structure. The absence of collisions from N-body codes is in particular a major shortcoming, as collisional processes are expected to significantly alter the results obtained from pure N-body runs. We present a new numerical model, to study the spatial structure of perturbed debris discs at dynamical and collisional steady-state. We focus on the competing effects between gravitational perturbations by a massive body (planet or star), collisional production of small grains, and radiation pressure placing these grains in possibly dynamically unstable regions. We consider a disc of parent bodies at dynamical steady-state, from which small radiation-pressure-affected grains are released in a series of runs, each corresponding to a different orbital position of the perturber, where particles are assigned a collisional destruction probability. These collisional runs produce successive position maps that are then recombined, following a complex procedure, to produce surface density profiles for each orbital position of the perturbing body. We apply our code to the case of a circumprimary disc in a binary. We find pronounced structures inside and outside the dynamical stability regions. For low $e_B$, the disc&#39;s structure is time varying, with spiral arms in the dynamically &#34;forbidden&#34; region precessing with the companion star. For high $e_B$, the disc is strongly asymmetric but time invariant, with a pronounced density drop in the binary&#39;s periastron direction. (better resolution version of the paper at http://lesia.obspm.fr/perso/philippe-thebault/theb2011.pdf)

preprint2011arXiv

Against all odds? Forming the planet of the HD196885 binary

HD196885Ab is the most &#34;extreme&#34; planet-in-a-binary discovered to date, whose orbit places it at the limit for orbital stability. The presence of a planet in such a highly perturbed region poses a clear challenge to planet-formation scenarios. We investigate this issue by focusing on the planet-formation stage that is arguably the most sensitive to binary perturbations: the mutual accretion of kilometre-sized planetesimals. To this effect we numerically estimate the impact velocities $dv$ amongst a population of circumprimary planetesimals. We find that most of the circumprimary disc is strongly hostile to planetesimal accretion, especially the region around 2.6AU (the planet&#39;s location) where binary perturbations induce planetesimal-shattering $dv$ of more than 1km/s. Possible solutions to the paradox of having a planet in such accretion-hostile regions are 1) that initial planetesimals were very big, at least 250km, 2) that the binary had an initial orbit at least twice the present one, and was later compacted due to early stellar encounters, 3) that planetesimals did not grow by mutual impacts but by sweeping of dust (the &#34;snowball&#34; growth mode identified by Xie et al., 2010b), or 4) that HD196885Ab was formed not by core-accretion but by the concurent disc instability mechanism. All of these 4 scenarios remain however highly conjectural.

preprint2011arXiv

On how the optical depth tunes the effects of ISM neutral atom flow on debris disks

The flux of ISM neutral atoms surrounding stars and their environment affects the motion of dust particles in debris disks, causing a significant dynamical evolution. Large values of eccentricity and inclination can be excited and strong correlations settle in among the orbital angles. This dynamical behaviour, in particular for bound dust grains, can potentially cause significant asymmetries in dusty disks around solar type stars which might be detected by observations. However, the amount of orbital changes due to this non--gravitational perturbation is strongly limited by the collisional lifetime of dust particles. We show that for large values of the disk&#39;s optical depth the influence of ISM flow on the disk shape is almost negligible because the grains are collisionally destroyed before they can accumulate enough orbital changes due to the ISM perturbations. On the other hand, for values smaller than $10^{-3}$, peculiar asymmetric patterns appear in the density profile of the disk when we consider 1-10 mum grains, just above the blow-out threshold. The extent and relevance of these asymmetries grow for lower values of the optical depth. An additional sink mechanism, which may prevent the formation of large clumps and warping in the disks is related to the fast inward migration due to the drag component of the forces. When a significant eccentricity is pumped up by the ISM perturbations, the drag forces (Poynting-Robertson and in particular ISM drag) drive the disk particles on fast migrating tracks leading them into the star on a short timescale. It is then expected that disks with small optical depth expand inside the parent body ring all the way towards the star while disks with large optical depth would not significantly extend inside.

preprint2011arXiv

Planet Formation In Highly Inclined Binary Systems I. Planetesimals Jump Inwards And Pile Up

Most detected planet-bearing binaries are in wide orbits, for which a high inclination, $i_B$, between the binary orbital plane and the plane of the planetary disk around the primary is likely to be common. In this paper, we investigate the intermediate stages - from planetesimals to planetary embryos/cores - of planet formation in such highly inclined cases. Our focus is on the effects of gas drag on the planetesimals&#39; orbital evolution, in particular on the evolution of the planetesimals&#39; semimajor axis distribution and their mutual relative velocities. We first demonstrate that a non-evolving axisymmetric disk model is a good approximation for studying the effects of gas drag on a planetesimal in the highly inclined case ($30^\circ <i_B<150^\circ$). We then find that gas drag plays a crucial role, and the results can be generally divided into two categories, i.e., the Kozai-on regime and the Kozai-off regime, depending on the specific value of $i_B$. For both regimes, a robust outcome over a wide range of parameters is that, planetesimals migrate/jump inwards and pile up, leading to a severely truncated and dense planetesimal disk around the primary. In this compact and dense disk, collision rates are high but relative velocities are low, providing conditions which are favorable for planetesimal growth, and potentially allow for the subsequent formation of planets.

preprint2010arXiv

Debris discs in binaries: a numerical study

Debris disc analysis and modelling provide crucial information about the structure and the processes at play in extrasolar planetary systems. In binary systems, this issue is more complex because the disc should in addition respond to the companion star&#39;s perturbations. We explore the dynamical evolution of a collisionally active debris disc for different initial parent body populations, diverse binary configurations and optical depths. We focus on the radial extent and size distribution of the disc at a stationary state. We numerically follow the evolution of $10^{5}$ massless small grains, initially produced from a circumprimary disc of parent bodies following a size distribution in $dN \propto s^{-3.5}$ds . Grains are submitted to both stars&#39; gravity as well as radiation pressure. In addition, particles are assigned an empirically derived collisional lifetime. For all the binary configurations the disc extends far beyond the critical semimajor axis $a_crit$ for orbital stability. This is due to the steady production of small grains, placed on eccentric orbits reaching beyond $a_crit$ by radiation pressure. The amount of matter beyond acrit depends on the balance between collisional production and dynamical removal rates: it increases for more massive discs as well as for eccentric binaries. Another important effect is that, in the dynamically stable region, the disc is depleted from its smallest grains. Both results could lead to observable signatures. We have shown that a companion star can never fully truncate a collisionally active disc. For eccentric companions, grains in the unstable regions can significantly contribute to the thermal emission in the mid-IR. Discs with sharp outer edges, especially bright ones such as HR4796A, are probably shaped by other mechanisms.

preprint2010arXiv

From Dust To Planetesimal: The Snowball Phase ?

The standard model of planet formation considers an initial phase in which planetesimals form from a dust disk, followed by a phase of mutual planetesimal-planetesimal collisions, leading eventually to the formation of planetary embryos. However, there is a potential transition phase (which we call the &#34;snowball phase&#34;), between the formation of the first planetesimals and the onset of mutual collisions amongst them, which has often been either ignored or underestimated in previous studies. In this snowball phase, isolated planetesimals move on Keplerian orbits and grow solely via the direct accretion of sub-cm sized dust entrained with the gas in the protoplanetary disk. Using a simplified model in which planetesimals are progressively produced from the dust, we consider the expected sizes to which the planetesimals can grow before mutual collisions commence and derive the dependence of this size on a number of critical parameters, including the degree of disk turbulence, the planetesimal size at birth and the rate of planetesimal creation. For systems in which turbulence is weak and the planetesimals are created at a low rate and with relatively small birth size, we show that the snowball growth phase can be very important, allowing planetesimals to grow by a factor of 10^6 in mass before mutual collisions take over. In such cases, the snowball growth phase can be the dominant mode to transfer mass from the dust to planetesimals. Moreover, such growth can take place within the typical lifetime of a protoplanetary gas disk. A noteworthy result is that ... ...(see the paper). For the specific case of close binaries such as Alpha Centauri ... ... (see the paper). From a more general perspective, these preliminary results suggest that an efficient snowball growth phase provides a large amount of &#34;room at the bottom&#34; for theories of planet formation.

preprint2009arXiv

Vertical structure of debris discs

The vertical thickness of debris discs is often used as a measure of these systems&#39; dynamical excitation and as clues to the presence of hidden massive perturbers such as planetary embryos. However, this argument could be flawed because the observed dust should be naturally placed on inclined orbits by the combined effect of radiation pressure and mutual collisions. We critically reinvestigate this issue and numerically estimate what the &#34;natural&#34; vertical thickness of a collisionally evolving disc is, in the absence of any additional perturbing body. We use a deterministic collisional code, following the dynamical evolution of a population of indestructible test grains suffering mutual inelastic impacts. Grain differential sizes as well as the effect of radiation pressure are taken into account. We find that, under the coupled effect of radiation pressure and collisions, grains naturally acquire inclinations of a few degrees. The disc is stratified with respect to grain sizes, with the smallest grains having the largest vertical dispersion and the bigger ones clustered closer to the midplane. Debris discs should have a minimum &#34;natural&#34; observed aspect ratio $h_{min}\sim 0.04\pm0.02$ at visible to mid-IR wavelengths where the flux is dominated by the smallest bound grains. These values are comparable to the estimated thicknesses of many vertically resolved debris discs, as is illustrated with the specific example of AU Mic. For all systems with $h \sim h_{min}$, the presence (or absence) of embedded perturbing bodies cannot be inferred from the vertical dispersion of the disc

preprint2007arXiv

Collisional Velocities and Rates in Resonant Planetesimal Belts

We consider a belt of small bodies around a star, captured in one of the external or 1:1 mean-motion resonances with a massive perturber. The objects in the belt collide with each other. Combining methods of celestial mechanics and statistical physics, we calculate mean collisional velocities and collisional rates, averaged over the belt. The results are compared to collisional velocities and rates in a similar, but non-resonant belt, as predicted by the particle-in-a-box method. It is found that the effect of the resonant lock on the velocities is rather small, while on the rates more substantial. The collisional rates between objects in an external resonance are by about a factor of two higher than those in a similar belt of objects not locked in a resonance. For Trojans under the same conditions, the collisional rates may be enhanced by up to an order of magnitude. Our results imply, in particular, shorter collisional lifetimes of resonant Kuiper belt objects in the solar system and higher efficiency of dust production by resonant planetesimals in debris disks around other stars.