A large fraction of stars is found to be part of binary or higher-order multiple systems. The ubiquity of planets found around single stars raises the question if and how planets in binary systems may form. Protoplanetary disks are the birthplaces of planets, and their characterization is crucial in order to understand the planet formation process. Our aim is to characterize the morphology of the GG Tau A disk, one of the largest and most massive circumbinary disks, and trace evidence for binary-disk interactions. We obtained observations in polarized scattered light of GG Tau A using the SPHERE/IRDIS instrument in the H-band filter. We analyze the observed disk morphology and substructures. We run 2D hydrodynamical models simulating the evolution of the circumbinary ring over the lifetime of the disk. The disk, as well as the cavity and the inner region are highly structured with several shadowed regions, spiral structures, and streamer-like filaments, some of them detected for the first time. The streamer-like filaments appear to connect the outer ring with the northern arc. Their azimuthal spacing suggests that they may be generated by periodic perturbations by the binary, tearing off material from the inner edge of the outer disk once during each orbit. By comparing observations to hydrodynamical simulations we find that the main features, in particular the gap size, as well as the spiral and streamer filaments, can be qualitatively explained by the gravitational interactions of a binary with semi-major axis of $\sim$35 au on an orbit coplanar with the circumbinary ring.

SPHERE (Beuzit et al,. 2019) has now been in operation at the VLT for more than 5 years, demonstrating a high level of performance. SPHERE has produced outstanding results using a variety of operating modes, primarily in the field of direct imaging of exoplanetary systems, focusing on exoplanets as point sources and circumstellar disks as extended objects. The achievements obtained thus far with SPHERE (~200 refereed publications) in different areas (exoplanets, disks, solar system, stellar physics...) have motivated a large consortium to propose an even more ambitious set of science cases, and its corresponding technical implementation in the form of an upgrade. The SPHERE+ project capitalizes on the expertise and lessons learned from SPHERE to push high contrast imaging performance to its limits on the VLT 8m-telescope. The scientific program of SPHERE+ described in this document will open a new and compelling scientific window for the upcoming decade in strong synergy with ground-based facilities (VLT/I, ELT, ALMA, and SKA) and space missions (Gaia, JWST, PLATO and WFIRST). While SPHERE has sampled the outer parts of planetary systems beyond a few tens of AU, SPHERE+ will dig into the inner regions around stars to reveal and characterize by mean of spectroscopy the giant planet population down to the snow line. Building on SPHERE's scientific heritage and resounding success, SPHERE+ will be a dedicated survey instrument which will strengthen the leadership of ESO and the European community in the very competitive field of direct imaging of exoplanetary systems. With enhanced capabilities, it will enable an even broader diversity of science cases including the study of the solar system, the birth and death of stars and the exploration of the inner regions of active galactic nuclei.

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.

The new NIKA2 camera at the IRAM 30m radiotelescope was used to observe three known debris disks in order to constrain the SED of their dust emission in the millimeter wavelength domain. We have found that the spectral index between the two NIKA2 bands (1mm and 2mm) is consistent with the Rayleigh-Jeans regime (lambda^{-2}), unlike the steeper spectra (lambda^{-3}) measured in the submillimeter-wavelength domain for two of the three disks $-$ around the stars Vega and HD107146. We provide a succesful proof of concept to model this spectral inversion in using two populations of dust grains, those smaller and those larger than a grain radius a0 of 0.5mm. This is obtained in breaking the slope of the size distribution and the functional form of the absorption coefficient of the standard model at a0. The third disk - around the star HR8799 - does not exhibit this spectral inversion but is also the youngest.

Context. Direct imaging of debris discs gives important information about their nature, their global morphology, and allows us to identify specific structures possibly in connection with the presence of gravitational perturbers. It is the most straightforward technique to observe planetary systems as a whole. Aims. We present the first resolved images of the debris disc around the young F-type star HD 160305, detected in scattered light using the VLT/SPHERE instrument in the near infrared. Methods. We used a post-processing method based on angular differential imaging and synthetic images of debris discs produced with a disc modelling code (GRaTer) to constrain the main characteristics of the disc around HD 160305. All of the point sources in the field of the IRDIS camera were analysed with an astrometric tool to determine whether they are bound objects or background stars. Results. We detect a very inclined (~ 82{\deg}) ring-like debris disc located at a stellocentric distance of about 86au (deprojected width ~27 au). The disc displays a brightness asymmetry between the two sides of the major axis, as can be expected from scattering properties of dust grains. We derive an anisotropic scattering factor g>0.5. A second right-left asymmetry is also observed with respect to the minor axis. We measure a surface brightness ratio of 0.73 $\pm$ 0.18 between the bright and the faint sides. Because of the low signal-to-noise ratio (S/N) of the images we cannot easily discriminate between several possible explanations for this left-right asymmetry, such as perturbations by an unseen planet, the aftermath of the breakup of a massive planetesimal, or the pericenter glow effect due to an eccentric ring. Two epochs of observations allow us to reject the companionship hypothesis for the 15 point sources present in the field.

The scattering properties of the dust originating from debris discs are still poorly known. The analysis of scattered light is however a powerful remote-sensing tool to understand the physical properties of dust particles orbiting other stars. Scattered light is indeed widely used to characterise the properties of cometary dust in the solar system. We aim to measure the morphology and scattering properties of the dust from the debris ring around HR4796A in polarised optical light. We obtained high-contrast polarimetric images of HR4796A in the wavelength range 600-900nm with the SPHERE / ZIMPOL instrument on the Very Large Telescope. We measured for the first time the polarised phase function of the dust in a debris system over a wide range of scattering angles in the optical. We confirm that it is incompatible with dust particles being compact spheres under the assumption of the Mie theory, and propose alternative scenarios compatible with the observations, such as particles with irregular surface roughness or aggregate particles.

The presence of submicron grains has been inferred in several debris discs, despite the fact that these particles should be blown out by stellar radiation pressure on very short timescales. So far, no fully satisfying explanation has been found for this apparent paradox. We investigate the possibility that the observed abundances of submicron grains could be "naturally" produced in bright debris discs, where the high collisional activity produces them at a rate high enough to partially compensate for their rapid removal. We also investigate to what extent this potential presence of small grains can affect our understanding of some debris disc characteristics. We use a code following the collisional evolution of a debris disc down to submicron grains far below the limiting blow-out size $s_{blow}$. We explore different configurations: A and G stars, cold and warm discs, "bright" and "very bright" systems. We find that, in bright discs (fractional luminosity $>10^{-3}$) around A stars, there is always a high-enough amount of submicron grains to leave detectable signatures, both in scattered-light, where the disc's color becomes blue, and in the mid-IR ($10<\lambda<20\mu$m), where it boosts the disc's luminosity by at least a factor of 2 and induces a pronounced silicate solid-state band around $10\mu$m. We also show that, with this additional contribution of submicron grains, the SED can mimic that of two debris belts separated by a factor of 2 in radial distance. For G stars, the effect of $s<s_{blow}$ grains remains limited in the spectra, in spite of the fact that they dominate the system's geometrical cross section. We also find that, for all considered cases, the halo of small (bound and unbound) grains that extends far beyond the main disc contributes to $\sim50$% of the flux up to $\lambda\sim50\mu$m wavelengths.

We detect the accreting low-mass companion HD142527B at a separation of 73 mas (11.4 au) from the star. No other companions with mass greater than 10 MJ are visible in the field of view of IFS (\sim 100 au centered on the star) or in the IRDIS field of view (\sim 400 au centered on the star). Measurements from IFS, SAM IFS, and IRDIS suggest an M6 spectral type for HD142527B, with an uncertainty of one spectral subtype, compatible with an object of M=0.11 \pm 0.06 MSun and R=0.15 \pm 0.07 RSun. The determination of the mass remains a challenge using contemporary evolutionary models, as they do not account for the energy input due to accretion from infalling material. We consider that the spectral type of the secondary may also be earlier than the type we derived from IFS spectra. From dynamical considerations, we further constrain the mass to 0.26^{+0.16}_{-0.14} MSun , which is consistent with both our spectroscopic analysis and the values reported in the literature. Following previous methods, the lower and upper dynamical mass values correspond to a spectral type between M2.5 and M5.5 for the companion. By fitting the astrometric points, we find the following orbital parameters: a period of P=35-137 yr; an inclination of i=121-130 deg.; , a value of Omega=124-135 deg for the longitude of node, and an 68% confidence interval of \sim 18 - 57 au for the separation at periapsis. Eccentricity and time at periapsis passage exhibit two groups of values: \sim0.2-0.45 and \sim0.45-0.7 for e, and \sim 2015-2020 and \sim2020-2022 for T_0. While these orbital parameters might at first suggest that HD142527B is not the companion responsible for the outer disk truncation, a previous hydrodynamical analysis of this system showed that they are compatible with a companion that is able to produce the large cavity and other observed features.

For a comprehensive understanding of planetary formation and evolution, we need to investigate the environment in which planets form: circumstellar disks. Here we present high-contrast imaging observations of V4046 Sagittarii, a 20-Myr-old close binary known to host a circumbinary disk. We have discovered the presence of rotating shadows in the disk, caused by mutual occultations of the central binary. Shadow-like features are often observed in disks\cite{garufi,marino15}, but those found thus far have not been due to eclipsing phenomena. We have used the phase difference due to light travel time to measure the flaring of the disk and the geometrical distance of the system. We calculate a distance that is in very good agreement with the value obtained from the Gaia mission's Data Release 2 (DR2), and flaring angles of $\alpha = 6.2 \pm 0.6 $ deg and $\alpha = 8.5 \pm 1.0 $ deg for the inner and outer disk rings, respectively. Our technique opens up a path to explore other binary systems, providing an independent estimate of distance and the flaring angle, a crucial parameter for disk modelling.

The recent discovery of hot dust grains in the vicinity of main-sequence stars has become a hot issue among the scientific community of debris disks. Hot grains must have been enormously accumulated near their sublimation zones, but it is a mystery how such a high concentration of hot grains is sustained. The most difficult conundrum is that the size of hot dust grains is estimated to lie in the submicrometer range, while submicrometer-sized grains are instantly swept away from near-stellar environments by stellar radiation pressure. One and only mechanism proposed for prolonging the residence time of hot grains in the near-stellar environments is trapping of charged nanoparticles by stellar magnetic fields. We revisit the model of magnetic grain trapping around main-sequence stars of various spectral classes by taking into account sublimation and electric charging of the grains. The model of magnetic grain trapping predicts that hot dust grains are present in the vicinity of main-sequence stars with high rotation velocities and intermediate magnetic-field strengths. On the contrary, we find that the detection of hot dust grains has no correlation with the rotation velocities of central stars nor the magnetic field strengths of the stars. Our numerical evaluation of electric grain charging indicates that the surface potential of submicrometer-sized grains in the vicinity of main-sequence stars is typically 4-5 V, which is one order of magnitude smaller than the value assumed by the model of magnetic grain trapping. On the basis of our numerical simulation on sublimation of dust grains in the vicinity of a star, it turns out that their lives end due to sublimation in a timescale much shorter than the period of one revolution at the gyroradius. It is, therefore, infeasible to dynamically extend the dwell time of hot grains inside the sublimation zone by magnetic trapping, ...

The TRAPPIST-1 system is unique in that it has a chain of seven terrestrial Earth-like planets located close to or in its habitable zone. In this paper, we study the effect of potential cometary impacts on the TRAPPIST-1 planets and how they would affect the primordial atmospheres of these planets. We consider both atmospheric mass loss and volatile delivery with a view to assessing whether any sort of life has a chance to develop. We ran N-body simulations to investigate the orbital evolution of potential impacting comets, to determine which planets are more likely to be impacted and the distributions of impact velocities. We consider three scenarios that could potentially throw comets into the inner region (i.e within 0.1au where the seven planets are located) from an (as yet undetected) outer belt similar to the Kuiper belt or an Oort cloud: Planet scattering, the Kozai-Lidov mechanism and Galactic tides. For the different scenarios, we quantify, for each planet, how much atmospheric mass is lost and what mass of volatiles can be delivered over the age of the system depending on the mass scattered out of the outer belt. We find that the resulting high velocity impacts can easily destroy the primordial atmospheres of all seven planets, even if the mass scattered from the outer belt is as low as that of the Kuiper belt. However, we find that the atmospheres of the outermost planets f, g and h can also easily be replenished with cometary volatiles (e.g. $\sim$ an Earth ocean mass of water could be delivered). These scenarios would thus imply that the atmospheres of these outermost planets could be more massive than those of the innermost planets, and have volatiles-enriched composition.

A warm/hot dust component (at temperature $>$ 300K) has been detected around $\sim$ 20% of stars. This component is called "exozodiacal dust" as it presents similarities with the zodiacal dust detected in our Solar System, even though its physical properties and spatial distribution can be significantly different. Understanding the origin and evolution of this dust is of crucial importance, not only because its presence could hamper future detections of Earth-like planets in their habitable zones, but also because it can provide invaluable information about the inner regions of planetary systems. In this review, we present a detailed overview of the observational techniques used in the detection and characterisation of exozodiacal dust clouds ("exozodis") and the results they have yielded so far, in particular regarding the incidence rate of exozodis as a function of crucial parameters such as stellar type and age, or the presence of an outer cold debris disc. We also present the important constraints that have been obtained, on dust size distribution and spatial location, by using state-of-the-art radiation transfer models on some of these systems. Finally, we investigate the crucial issue of how to explain the presence of exozodiacal dust around so many stars (regardless of their ages) despite the fact that such dust so close to its host star should disappear rapidly due to the coupled effect of collisions and stellar radiation pressure. Several potential mechanisms have been proposed to solve this paradox and are reviewed in detail in this paper. The review finishes by presenting the future of this growing field.

A collisional avalanche is set off by the breakup of a large planetesimal, releasing small unbound grains that enter a debris disc located further away from the star, triggering there a collisional chain reaction that can potentially create detectable transient structures. We explore this mechanism, using for the first time a code coupling dynamical and collisional evolutions, and investigate if avalanches could explain the short-term luminosity variations observed in some extremely bright discs. We consider two set-ups: a cold disc case, with a dust release at 10au and an outer disc extending from 50 to 120au, and a warm disc case with the release at 1au and a 5-12au outer disc. We find that avalanches could leave detectable structures on resolved images, for both cold and warm disc cases, in discs with optical depth $\tau$ of a few $10^{-3}$, provided that large dust masses ($\gtrsim$10$^{20}$-5$\times$10$^{22}$g) are initially released. The integrated photometric excess due to an avalanche is limited, less than 10% for these released dust masses, peaking in the mid-IR and becoming insignificant beyond $\sim$40-50$\mu$m. Contrary to earlier studies, we do not obtain stronger avalanches when increasing $\tau$ to higher values. Likewise, we do not observe a significant luminosity deficit, as compared to the pre-avalanche level, after the passage of the avalanche. These two results concur to make avalanches an unlikely explanation for the sharp luminosity drops observed in some extremely bright debris discs. The ideal configuration for observing an avalanche would be a two-belt structure, with an inner belt of fractional luminosity >10$^{-4}$ where breakups of massive planetesimals occur, and a more massive outer belt, with $\tau$ of a few $10^{-3}$, into which the avalanche chain reaction develops and propagates.

We explore the multiplicity of exoplanet host stars with high-resolution images obtained with VLT/SPHERE. Two different samples of systems were observed: one containing low-eccentricity outer planets, and the other containing high-eccentricity outer planets. We find that 10 out of 34 stars in the high-eccentricity systems are members of a binary, while the proportion is 3 out of 27 for circular systems. Eccentric-exoplanet hosts are, therefore, significantly more likely to have a stellar companion than circular-exoplanet hosts. The median magnitude contrast over the 68 data sets is 11.26 and 9.25, in H and K, respectively, at 0.30 arcsec. The derived detection limits reveal that binaries with separations of less than 50au are rarer for exoplanet hosts than for field stars. Our results also imply that the majority of high-eccentricity planets are not embedded in multiple stellar systems (24 out of 34), since our detection limits exclude the presence of a stellar companion. We detect the low-mass stellar companions of HD 7449 and HD 211847, both members of our high-eccentricity sample. HD 7449B was already detected by Rodigas et al (2016) and our independent observation is in agreement with this earlier work. HD 211847's substellar companion, previously detected by the radial velocity method, is actually a low-mass star seen face-on. The role of stellar multiplicity in shaping planetary systems is confirmed by this work, although it does not appear as the only source of dynamical excitation.

Spurred by the discovery of numerous exoplanets in multiple systems, binaries have become in recent years one of the main topics in planet formation research. Numerous studies have investigated to what extent the presence of a stellar companion can affect the planet formation process. Such studies have implications that can reach beyond the sole context of binaries, as they allow to test certain aspects of the planet formation scenario by submitting them to extreme environments. We review here the current understanding on this complex problem. We show in particular how each of the different stages of the planet-formation process is affected differently by binary perturbations. We focus especially on the intermediate stage of kilometre-sized planetesimal accretion, which has proven to be the most sensitive to binarity and for which the presence of some exoplanets observed in tight binaries is difficult to explain by in-situ formation following the "standard" planet-formation scenario. Some tentative solutions to this apparent paradox are presented. The last part of our review presents a thorough description of the problem of planet habitability, for which the binary environment creates a complex situation because of the presence of two irradation sources of varying distance.

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 $\mu$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).

The precise mechanisms that provide the non-radiative energy for heating the chromosphere and the corona of the Sun and those of other stars constitute an active field of research. By studying stellar chromospheres one aims at identifying the relevant physical processes. Defining the permittable extent of the parameter space can also serve as a template for the Sun-as-a-star. Earlier observations with Herschel and APEX have revealed the temperature minimum of alpha Cen, but these were unable to spatially resolve the binary into individual components. With the data reported here, we aim at remedying this shortcoming. Furthermore, these earlier data were limited to the wavelength region between 100 and 870mu. In the present context, we intend to extend the spectral mapping to longer wavelengths, where the contrast between stellar photospheric and chromospheric emission becomes increasingly evident. ALMA is particularly suited to point sources, such as unresolved stars. ALMA provides the means to achieve our objectives with both its high sensitivity of the collecting area for the detection of weak signals and the high spatial resolving power of its adaptable interferometer for imaging close multiple stars. This is the first detection of main-sequence stars at a wavelength of 3mm. Furthermore, the individual components of the binary alpha CenAB are clearly detected and spatially well resolved at all ALMA wavelengths. The high S/N of these data permit accurate determination of their relative flux ratios. The previously obtained flux ratio of 0.44, which was based on measurements in the optical and at 70mu, is consistent with the present ALMA results, albeit with a large error bar. Given the distinct difference in their cyclic activity, the similarity of their submm SEDs appears surprising.

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$\mu$m (marginally), 11.4 and 15.5$\mu$m would show the inner disc structures while 23$\mu$m images would display the outer disc asymmetries.

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.

[abridged] Aims. Our Herschel Open Time Key Programme DUNES aims at detecting and characterizing debris disks around nearby, sun-like stars. In addition to the statistical analysis of the data, the detailed study of single objects through spatially resolving the disk and detailed modeling of the data is a main goal of the project. Methods. We obtained the first observations spatially resolving the debris disk around the sun-like star HIP 17439 (HD23484) using the instruments PACS and SPIRE on board the Herschel Space Observatory. Simultaneous multi-wavelength modeling of these data together with ancillary data from the literature is presented. Results. A standard single component disk model fails to reproduce the major axis radial profiles at 70 um, 100 um, and 160 um simultaneously. Moreover, the best-fit parameters derived from such a model suggest a very broad disk extending from few au up to few hundreds of au from the star with a nearly constant surface density which seems physically unlikely. However, the constraints from both the data and our limited theoretical investigation are not strong enough to completely rule out this model. An alternative, more plausible, and better fitting model of the system consists of two rings of dust at approx. 30 au and 90 au, respectively, while the constraints on the parameters of this model are weak due to its complexity and intrinsic degeneracies. Conclusions. The disk is probably composed of at least two components with different spatial locations (but not necessarily detached), while a single, broad disk is possible, but less likely. The two spatially well-separated rings of dust in our best-fit model suggest the presence of at least one high mass planet or several low-mass planets clearing the region between the two rings from planetesimals and dust.

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.

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 "at rest" 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.

The first attempt at developing a fully self-consistent code coupling dynamics and collisions to study debris discs (Kral, Thebault, Charnoz, 2013) is presented. So far, these two crucial mechanisms were studied separately, with N-body and statistical collisional codes respectively, because of stringent computational constraints. In particular, incorporating collisional effects (especially destructive collisions) into an N-body scheme was deemed an impossible task because of the exponential increase of particles it would imply. We present here a new model named LIDT-DD which is able to follow over long timescales the coupled evolution of dynamics (including radiation pressure) and collisions in a self-consistent way.

Recent observations from NASA'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.

We report on high contrast mid-infrared observations of Fomalhaut obtained with the Keck Interferometer Nuller (KIN) showing a small resolved excess over the level expected from the stellar photosphere. The measured null excess has a mean value of 0.35% +/- 0.10% between 8 and 11 microns and increases from 8 to 13 microns. Given the small field of view of the instrument, the source of this marginal excess must be contained within 2AU of Fomalhaut. This result is reminiscent of previous VLTI K-band observations, which implied the presence of a ~ 0.88% excess, and argued that thermal emission from hot dusty grains located within 6 AU from Fomalhaut was the most plausible explanation. Using a parametric 2D radiative transfer code and a Bayesian analysis, we examine different dust disk structures to reproduce both the near and mid-infrared data simultaneously. While not a definitive explanation of the hot excess of Fomalhaut, our model suggests that the most likely inner few AU disk geometry consists of a two-component structure, with two different and spatially distinct grain populations. The 2 to 11 microns data are consistent with an inner hot ring of very small (~ 10 to 300 nm) carbon-rich grains concentrating around 0.1AU. The second dust population consists of larger grains (size of a few microns to a few tens of microns) located further out in a colder region where regular astronomical silicates could survive, with an inner edge around 1AU. From a dynamical point of view, the presence of the inner concentration of sub-micron sized grains is surprising, as such grains should be expelled from the inner planetary system by radiation pressure within only a few years. This could either point to some inordinate replenishment rates (e.g. many grazing comets coming from an outer reservoir) or to the existence of some braking mechanism preventing the grains from moving out.

Spectral modeling of the large infrared excess in the Spitzer IRS spectra of HD 172555 suggests that there is more than 10^19 kg of sub-micron dust in the system. Using physical arguments and constraints from observations, we rule out the possibility of the infrared excess being created by a magma ocean planet or a circumplanetary disk or torus. We show that the infrared excess is consistent with a circumstellar debris disk or torus, located at approximately 6 AU, that was created by a planetary scale hypervelocity impact. We find that radiation pressure should remove submicron dust from the debris disk in less than one year. However, the system's mid-infrared photometric flux, dominated by submicron grains, has been stable within 4 percent over the last 27 years, from IRAS (1983) to WISE (2010). Our new spectral modeling work and calculations of the radiation pressure on fine dust in HD 172555 provide a self-consistent explanation for this apparent contradiction. We also explore the unconfirmed claim that 10^47 molecules of SiO vapor are needed to explain an emission feature at 8 um in the Spitzer IRS spectrum of HD 172555. We find that unless there are 10^48 atoms or 0.05 Earth masses of atomic Si and O vapor in the system, SiO vapor should be destroyed by photo-dissociation in less than 0.2 years. We argue that a second plausible explanation for the 8 um feature can be emission from solid SiO, which naturally occurs in submicron silicate "smokes" created by quickly condensing vaporized silicate.

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.

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's mass and orbit, as well as the disc'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'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.

Context. It has long been suggested that circumstellar disks surrounding young stars may be the signposts of planets, and still more since the recent discoveries of embedded substellar companions. The planet-disk interaction may create, according to models, large structures, gaps, rings or spirals, in the disk. In that sense, the Herbig star HD 142527 is particularly compelling as, its massive disk displays intriguing asymmetries that suggest the existence of a dynamical peturber of unknown nature. Aims. Our goal was to obtain deep thermal images of the close circumstellar environment of HD 142527 to re-image the reported close-in structures (cavity, spiral arms) of the disk and to search for stellar and substellar companions that could be connected to their presence. Results. The circumstellar environment of HD 142527 is revealed at an unprecedented spatial resolution down to the sub arcsecond level for the first time at 3.8 microns. Our images reveal important radial and azimuthal asymmetries which invalidate an elliptical shape for the disk as previously proposed. It rather suggests a bright inhomogeneous spiral arm plus various fainter spiral arms. We also confirm an inner cavity down to 30 AU and two important dips at position angles of 0 and 135 deg. The detection performance in angular differential imaging enables the exploration of the planetary mass regime for projected physical separations as close as 40 AU. The use of our detection map together with Monte Carlo simulations sets stringent constraints on the presence of planetary mass, brown dwarf or stellar companions as a function of the semi-major axis. They severely constrain the presence of massive giant planets with semi-major axis beyond 50AU, i.e. probably within the large disk's cavity that radially extends up to 145 AU or even further outside.

HR4796 is a young, early A-type star harbouring a well structured debris disk, shaped as a ring with sharp inner edges. It forms with the M-type star HR4796B a binary system, with a proj. sep. ~560 AU. Our aim is to explore the surroundings of HR4796A and B, both in terms of extended or point-like structures. Adaptive optics images at L'-band were obtained with NaCo in Angular Differential Mode and with Sparse Aperture Masking (SAM). We analyse the data as well as the artefacts that can be produced by ADI reduction on an extended structure with a shape similar to that of HR4796A dust ring. We determine constraints on the presence of companions using SAM and ADI on HR4796A, and ADI on HR4796B. We also performed dynamical simulations of a disk of planetesimals and dust produced by collisions, perturbed by a planet located close to the disk outer edge. The disk ring around HR4796A is well resolved. We highlight the potential effects of ADI reduction of the observed disk shape and surface brightness distribution, and side-to-side asymmetries. No planet is detected around the star, with masses as low as 3.5 M_Jup at 0.5" (58 AU) and less than 3 M_Jup in the 0.8-1" range along the semi-major axis. We exclude massive brown dwarfs at separations as close as 60 mas (4.5 AU) from the star thanks to SAM data. The detection limits obtained allow us to exclude a possible close companion to HR4796A as the origin of the offset of the ring center with respect to the star; they also allow to put interesting constraints on the (mass, separation) of any planet possibly responsible for the inner disk steep edge. Using detailed dynamical simulations, we show that a giant planet orbiting outside the ring could sharpen the disk outer edge and reproduce the STIS images published by Schneider et al. (2009).

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.

With the uniquely high contrast within 0.1" (\Delta mag(L') = 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' of 0.016+/-0.007, 0.012+/-0.008, 0.0086+/-0.0011 respectively (3-\sigma 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.

Discs in binaries have a complex behavior because of the perturbations of the companion star. Planet formation in binary-star systems both depend on the companion star parameters and on the properties of the circumstellar disc. An eccentric disc may increase the impact velocity of planetesimals and therefore jeopardize the accumulation process. We model the evolution of discs in close binaries including the effects of self-gravity and adopting different prescriptions to model the disc's radiative properties. We focus on the dynamical properties and evolutionary tracks of the discs. We use the hydrodynamical code FARGO and we include in the energy equation heating and cooling effects. Radiative discs have a lower disc eccentricity compared to locally isothermal discs with same temperature profile. As a consequence, we do not observe the formation of an internal elliptical low density region as in locally isothermal disc models. However, the disc eccentricity depends on the disc mass through the opacities. Akin to locally isothermal disc models, self-gravity forces the disc's longitude of pericenter to librate about a fixed orientation with respect to the binary apsidal line ($\pi$). The disc's radiative properties play an important role in the evolution of discs in binaries. A radiative disc has an overall shape and internal structure that are significantly different compared to a locally isothermal disc with same temperature profile. This is an important finding both for describing the evolutionary track of the disc during its progressive mass loss, and for planet formation since the internal structure of the disc is relevant for planetesimals growth in binary systems. The non-symmetrical distribution of mass in these discs causes large eccentricities for planetesimals that may affect their growth.

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's structure is time varying, with spiral arms in the dynamically "forbidden" 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's periastron direction. (better resolution version of the paper at http://lesia.obspm.fr/perso/philippe-thebault/theb2011.pdf)

We present Herschel PACS 100 and 160 micron observations of the solar-type stars alpha Men, HD 88230 and HD 210277, which form part of the FGK stars sample of the Herschel Open Time Key Programme (OTKP) DUNES (DUst around NEarby S tars). Our observations show small infrared excesses at 160 micron for all three stars. HD 210277 also shows a small excess at 100 micron, while the 100 micron fluxes of alpha Men and HD 88230 agree with the stellar photospheric predictions. We attribute these infrared excesses to a new class of cold, faint debris discs. alpha Men and HD 88230 are spatially resolved in the PACS 160 micron images, while HD 210277 is point-like at that wavelength. The projected linear sizes of the extended emission lie in the range from ~ 115 to ~ 250 AU. The estimated black body temperatures from the 100 and 160 micron fluxes are $\lesssim$ 22 K, while the fractional luminosity of the cold dust is Ldust/Lstar ~ 10E-6, close to the luminosity of the Solar-System's Kuiper belt. These debris discs are the coldest and faintest discs discovered so far around mature stars and cannot easily be explained by invoking "classical" debris disc models.

We have analyzed Spitzer and NASA/IRTF 2 - 35 \mum spectra of the warm, ~350 K circumstellar dust around the nearby MS star {\eta} Corvi (F2V, 1.4 \pm 0.3 Gyr). The spectra show clear evidence for warm, water- and carbon-rich dust at ~3 AU from the central star, in the system's Terrestrial Habitability Zone. Spectral features due to ultra-primitive cometary material were found, in addition to features due to impact produced silica and high temperature carbonaceous phases. At least 9 x 10^18 kg of 0.1 - 100 \mum warm dust is present in a collisional equilibrium distribution with dn/da ~ a^-3.5, the equivalent of a 130 km radius KBO of 1.0 g/cm^3 density and similar to recent estimates of the mass delivered to the Earth at 0.6 - 0.8 Gyr during the Late Heavy Bombardment. We conclude that the parent body was a Kuiper-Belt body or bodies which captured a large amount of early primitive material in the first Myrs of the system's lifetime and preserved it in deep freeze at ~150 AU. At ~1.4 Gyr they were prompted by dynamical stirring of their parent Kuiper Belt into spiraling into the inner system, eventually colliding at 5-10 km/sec with a rocky planetary body of mass \leq MEarth at ~3 AU, delivering large amounts of water (>0.1% of MEarth's Oceans) and carbon-rich material. The Spitzer spectrum also closely matches spectra reported for the Ureilite meteorites of the Sudan Almahata Sitta fall in 2008, suggesting that one of the Ureilite parent bodies was a KBO.

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'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.

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' orbital evolution, in particular on the evolution of the planetesimals' 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.

HD196885Ab is the most "extreme" 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'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 "snowball" 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.

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 "snowball phase"), 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 "room at the bottom" for theories of planet formation.

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'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' 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.

We study the evolution of circumstellar massive disks around the primary star of a binary system focusing on the computation of disk eccentricity. In particular, we concentrate on its dependence on the binary eccentricity. Self-gravity is included in our numerical simulations. Our standard model assumes a semimajor axis for the binary of 30 AU, the most probable value according to the present binary statistics.

We explore planet formation in binary systems around the central star where the protoplanetary disk plane is highly inclined with respect to the companion star orbit. This might be the most frequent scenario for binary separations larger than 40 AU, according to Hale (1994). We focus on planetesimal accretion and compute average impact velocities in the habitable region and up to 6 AU from the primary.

The vertical thickness of debris discs is often used as a measure of these systems' 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 "natural" 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 "natural" 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

Recent studies have shown that alpha Centauri B might be, from an observational point of view, an ideal candidate for the detection of an Earth-like planet in or near its habitable zone (0.5-0.9AU). We study here if such habitable planets can form, by numerically investigating the planet-formation stage which is probably the most sensitive to binarity effects: the mutual accretion of km-sized planetesimals. Using a state-of-the-art algorithm for computing the impact velocities within a test planetesimal population, we find that planetesimal growth is only possible, although marginally, in the innermost part of the HZ around 0.5AU. Beyond this point, the combination of secular perturbations by the binary companion and gas drag drive the mutual velocities beyond the erosion limit. Impact velocities might later decrease during the gas removal phase, but this probably happens too late for preventing most km-sized objects to be removed by inward drift, thus preventing accretion from starting anew. A more promising hypothesis is that the binary formed in a crowded cluster, where it might have been wider in its initial stages, when planetary formation was ongoing. We explore this scenario and find that a starting separation roughly 15 AU wider, or an eccentricity 2.5 times lower than the present ones are required to have an accretion-friendly environment in the whole HZ.

We numerically explore planet formation around alpha Cen A by focusing on the crucial planetesimals-to-embryos phase. Our code computes the relative velocity distribution, and thus the accretion vs. fragmentation trend, of planetesimal populations having any given size distribution. This is a critical aspect of planet formation in binaries since the pericenter alignment of planetesimal orbits due to the gravitational perturbations of the companion star and to gas friction strongly depends on size. We find that, for the nominal case of a MMSN gas disc, the region beyond 0.5AU from the primary is hostile to planetesimal accretion. In this area, impact velocities between different-size bodies are increased, by the differential orbital phasing, to values too high to allow mutual accretion. For any realistic size distribution for the planetesimal population, this accretion-inhibiting effect is the dominant collision outcome and the accretion process is halted. Results are robust with respect to the profile and density of the gas disc: except for an unrealistic almost gas-free case, the inner accretion safe area never extends beyond 0.75AU. We conclude that planet formation is very difficult in the terrestrial region around alpha Cen A, unless it started from fast-formed very large (>30km) planetesimals. Notwithstanding these unlikely initial conditions, the only possible explanation for the presence of planets around 1 AU from the star would be the hypothetical outward migration of planets formed closer to the star or a different orbital configuration in the binary's early history. Our conclusions differ from those of several studies focusing on the later embryos-to-planets stage, confirming that the planetesimals-to-embryos phase is more affected by binary perturbations.

Only very few solar-type stars exhibiting an infrared excess and harbouring planets are known to date. Indeed, merely a single case of a star-planet-disk system has previously been detected at submillimeter (submm) wavelengths. Consequently, one of our aims is to understand the reasons for these poor statistics, i.e., whether these results reflected the composition and/or the physics of the planetary disks or were simply due to observational bias and selection effects. Finding more examples would be very significant. The selected target, q1 Eri, is a solar-type star, which was known to possess a planet, q1 Eri b, and to exhibit excess emission at IRAS wavelengths, but had remained undetected in the millimeter regime. Therefore, submm flux densities would be needed to better constrain the physical characteristics of the planetary disk. Consequently, we performed submm imaging observations of q1 Eri. The detected dust toward q1 Eri at 870 micron exhibits the remarkable fact that the entire SED, from the IR to mm-wavelengths, is fit by a single temperature blackbody function (60 K). This would imply that the emitting regions are confined to a narrow region (ring) at radial distances much larger than the orbital distance of q1 Eri b, and that the emitting particles are considerably larger than some hundred micron. However, the 870 micron source is extended, with a full-width-half-maximum of roughly 600 AU. Therefore, a physically more compelling model also invokes a belt of cold dust (17 K), located at 300 AU from the star and about 60 AU wide. The minimum mass of 0.04 Mearth (3 Mmoon) of 1 mm-size icy ring-particles is considerable, given the stellar age of about 1 Gyr. These big grains form an inner edge at about 25 AU, which may suggest the presence of an unseen outer planet (q1 Eri c).

Observations of extrasolar planets reveal that planets can be found in close binary systems, where the semi-major axis of the binary orbit is less than 20 AU. The existence of these planets challenges planet formation theory, because the strong gravitational perturbations due to the companion increase encounter velocities between planetesimals and make it difficult for them to grow through accreting collisions. We study planetesimal encounter velocities in binary systems, where the planetesimals are embedded in a circumprimary gas disc that is allowed to evolve under influence of the gravitational perturbations of the companion star. We find that the encounter velocities between planetesimals of different size strongly depend on the gas disc eccentricity. In all cases studied, inclusion of the full gas dynamics increases the encounter velocity compared to the case of a static, circular gas disc. Full numerical parameter exploration is still impossible, but we derive analytical formulae to estimate encounter velocities between bodies of different sizes given the gas disc eccentricity. The gas dynamical evolution of a protoplanetary disc in a binary system tends to make planetesimal accretion even more difficult than in a static, axisymmetric gas disc.

Ring-like features have been observed in several debris discs. Outside the main ring, while some systems exhibit smooth surface brightness profiles (SB) that fall off roughly as r**-3.5, others display large luminosity drops at the ring's outer edge and steeper radial SB profiles. We seek to understand this diversity of outer edge profiles under the ``natural'' collisional evolution of the system, without invoking external agents such as planets or gas. We use a statistical code to follow the evolution of a collisional population, ranging from dust grains (submitted to radiation pressure) to planetesimals and initially confined within a belt (the 'birth ring'). The system typically evolves toward a "standard" steady state, with no sharp edge and SB \propto r**-3.5 outside the birth ring. Deviations from this standard profile, in the form of a sharp outer edge and a steeper fall-off, occur only when two parameters take their extreme values: 1) When the birth ring is so massive that it becomes radially optically thick for the smallest grains. However, the required disc mass is here probably too high to be realistic. 2) When the dynamical excitation of the dust-producing planetesimals is so low (<e> <0.01) that the smallest grains, which otherwise dominate the total optical depth, are preferentially depleted. This low-excitation case, although possibly not generic, cannot be ruled out by observations. Our "standard" profile provides a satisfactory explanation for a large group of debris discs with outer edges and SB falling as r**-3.5. Systems with sharper outer edges, barring other confining agents, could still be explained by ``natural'' collisional evolution if their dynamical excitation is very low. We show that such a dynamically-cold case provides a satisfactory fit for HR4796A

Numerical simulations of planetesimal accretion in circumprimary and circumbinary orbits are described. The secular perturbations by the com- panion star and gas drag are included in our models. We derive limits on the parameters of the binary system for which accretion and then planetary forma- tion are possible. In the circumbinary case we also outline the radial distance from the baricenter of the stars beyond which accumulation always occurs. Hy- drodynamical simulations are also presented to validate our N-body approach based on the axisymmetric approximation for the gas of the disk.

We put theoretical constraints on the presence and survival of icy grains in debris discs. Particular attention is paid to UV sputtering of water ice, which has so far not been studied in detail in this context. We present a photosputtering model based on available experimental and theoretical studies. We quantitatively estimate the erosion rate of icy and ice-silicate grains, under the influence of both sublimation and photosputtering, as a function of grain size, composition and distance from the star. The effect of erosion on the grain's location is investigated through numerical simulations coupling the grain size to its dynamical evolution. Our model predicts that photodesorption efficiently destroy ice in optically thin discs, even far beyond the sublimation snow line. For the reference case of beta Pictoris, we find that only > 5mm grains can keep their icy component for the age of the system in the 50-150AU region. When taking into account the collisional reprocessing of grains, we show that the water ice survival on grains improves (grains down to ~ 20 um might be partially icy). However, estimates of the amount of gas photosputtering would produce on such a hypothetical population of big icy grains lead to values for the OI column density that strongly exceed observational constraints for beta Pic, thus ruling out the presence of a significant amount of icy grains in this system. Erosion rates and icy grains survival timescales are also given for a set of 11 other debris disc systems. We show that, with the possible exception of M stars, photosputtering cannot be neglected in calculations of icy grain lifetimes.

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.

The inner Solar System contains a cloud of small (1-100 micron) dust grains created when small bodies-asteroids, comets, and Kuiper belt objects-collide and outgas. This dust cloud, the zodiacal cloud probably has extrasolar analogs, exozodiacal clouds. Exozodiacal clouds are related to debris disks, clouds of rocks and dust orbiting main sequence stars thought to represent the debris left over from planet formation. Some debris disks appear to contain distinct inner clouds that could be considered massive exozodiacal clouds (e.g. Koerner et al. 1998, Absil et al. 2006). This white paper addresses the need for future theoretical work on the dynamics of exozodiacal clouds. This theoretical work should help us discover new planets and understand exozodiacal clouds as astrophysical noise. So far, observations of nearby stars have not provided good constraints on the structures of exozodiacal clouds. But future observations probably will demand a better theoretical understanding of these systems.

We numerically investigate the possibility of planetesimal accretion in circumbinary disks, under the coupled influence of both stars' secular perturbations and friction due to the gaseous component of the protoplanetary disk. We focus on one crucial parameter: the distribution of encounter velocities between planetesimals in the 0.5 to 100km size range. An extended range of binary systems with differing orbital parameters is explored. The resulting encounter velocities are compared to the threshold velocities below which the net outcome of a collision is accumulation into a larger body instead of mass erosion. For each binary configuration, we derive the critical radial distance from the binary barycenter beyond which planetesimal accretion is possible. This critical radial distance is smallest for equal-mass binaries on almost circular orbits. It shifts to larger values for increasing eccentricities and decreasing mass ratio. The importance of the planetesimals' orbital alignments of planetesimals due to gas drag effects is discussed.

We present a new multi-annulus code for the study of collisionally evolving extended debris discs. We first aim to confirm results obtained for a single-annulus system, namely that the size distribution in "real" debris discs always departs from the theoretical collisional equilibrium $dN\proptoR^{-3.5}dR$ power law, especially in the crucial size range of observable particles (<1cm), where it displays a characteristic wavy pattern. We also aim at studying how debris discs density distributions, scattered light luminosity profiles, and SEDs are affected by the coupled effect of collisions and radial mixing due to radiation pressure affected small grains. The size distribution evolution is modeled from micron-sized grains to 50km-sized bodies. The model takes into account the crucial influence of radiation pressure-affected small grains. We consider the collisional evolution of a fiducial a=120AU radius disc with an initial surface density in $\Sigma(a)\propto a^{\alpha}$. We show that the system's radial extension plays a crucial role: in most regions the collisional and size evolution of the dust is imposed by small particles on eccentric or unbound orbits produced further inside the disc. The spatial distribution of small grains strongly departs from the initial profile, while the bigger objects, containing most of the system's mass, still follow the initial distribution. This has consequences on the scattered--light radial profiles which get significantly flatter, and we propose an empirical law to trace back the distribution of large unseen parent bodies from the observed profiles. We finally provide empirical formula for the collisional size distribution and collision timescale that can be used for future debris disc modeling.

We study to which extent planetesimal accretion is affected by the perturbing presence of a compagnon star. We concentrate on one crucial parameter: the distribution of encounter velocities within the planetesimal swarm. We numerically explore the evolution of this parameter taking into account the secular perturbations of the binary and friction due to the very likely presence of gas in the disk. $<\Delta v>$ maps are derived, for planetesimals of different sizes, for a total of 120 binary configurations (eccentricity eb and separation ab). We identify for each case 3 different accreting behaviours. 1) In regions where no significant dV increase is observed, ``classical'' (i.e., single-star) accretion is possible. 2) In regions where $dV>v_{ero}$, the threshold velocity above which all impacts are eroding, no accretion is possible and planet growth is stopped. 3) In between these 2 limiting behaviours, a large fraction of binary configurations leads to significant dV increase, but still below the erosion threshold. In this intermediate case, planetesimal growth can occur, but proceeds slower than in the single-star case, possibly following the so-called type II runaway groth mode.

We investigate classical planetesimal accretion in a binary star system of separation ab<50AU by numerical simulations, with particular focus on the region at a distance of 1 AU from the primary. The planetesimals orbit the primary, are perturbed by the companion and are in addition subjected to a gas drag force. We concentrate on the problem of relative velocities dv among planetesimals of different sizes. For various stellar mass ratios and binary orbital parameters we determine regions where dv exceed planetesimal escape velocities v_esc (thus preventing runaway accretion) or even the threshold velocity v_ero for which erosion dominates accretion. Gaseous friction has two crucial effects on the velocity distribution: it damps secular perturbations by forcing periastron alignment of orbits, but at the same time the size--dependence of this orbital alignment induces a significant dv increase between bodies of different sizes. This differential phasing effect proves very efficient and almost always increases dv to values preventing runaway accretion, except in a narrow domain of almost circular companion orbits. The erosion threshold dv>v_ero is reached in a wide (ab,eb) space for small (<10km) planetesimals, but in a much more limited region for bigger ~50km objects. In the intermediate v_esc<dv < v_ero domain, a possible growth mode would be the type II runaway growth identified by Kortenkmap et al.(2001)

We investigate in this paper the effect of gas drag on the dynamics of the dust particles in the edge-on Beta-Pictoris disc in order to derive an upper limit on the mass of gas in this system. Our study is motivated by the large uncertainties on the amount of gas in the Beta-Pictoris disc currently found in the literature. The dust particles are assumed to originate from a colliding annulus of planetesimals peaked around 100AU from the central star as proposed by Augereau et al.(2001). We consider the various gas densities that have been inferred from independent observing techniques and we discuss their impact on the dust dynamics and on the disc profile in scattered light along the midplane. We show that the observed scattered light profile of the disc cannot be properly reproduced if hydrogen gas number density at 117AU exceeds 10**4 cm**-3. This corresponds to an upper limit on the total gas mass of about 0.4 Mearth assuming the gas density profile inferred by Brandeker et al.(2004) and thus to a gas to dust mass ratio smaller than 1. Our approach therefore provides an independent diagnostic for a gas depletion in the Beta-Pictoris system relative to the dust disc. Such an approach could also be used to constrain the gas content of recently identified systems like the edge-on disc around AUmic.

We numerically investigate under which conditions the planet detected at 2.1 AU of Gamma-Cephei could form through the core-accretion scenario despite the perturbing presence of the highly eccentric companion star. We first show that the initial stage of runaway accretion of kilometer-sized planetesimals is possible within 2.5 AU from the central star only if large amounts of gas are present. In this case, gaseous friction induces periastron alignment of the orbits which reduces the otherwise high mutual impact velocities due to the companion's secular perturbations. The following stage of mutual accretion of large embryos is also modeled. According to our simulations, the giant impacts among the embryos always lead to a core of 10 Mearth within 10 Myr, the average lifetime of gaseous discs. However, the core always ends up within 1.5 AU from the central star. Either the core grows more quickly in the inner region of the disc, or it migrates inside by scattering the residual embryos

Dust particles observed in extrasolar planetary discs originate from undetectable km-sized bodies but this valuable information remains uninteresting if the theoretical link between grains and planetesimals is not properly known. We outline in this paper a numerical approach we developed in order to address this issue for the case of dust producing collisional cascades. The model is based on a particle-in-a-box method. We follow the size distribution of particles over eight orders of magnitude in radius taking into account fragmentation and cratering according to different prescriptions. A very particular attention is paid to the smallest particles, close to the radiation pressure induced cut-off size $R_{pr}$, which are placed on highly eccentric orbits by the stellar radiation pressure....(abstract continued in the uploaded paper)

We have started since 1997 the Meudon Multicolor Survey of Outer Solar System Objects with the aim of collecting a large and homogeneous set of color data for Trans-Neptunian and Centaurs objects [...] We have a combined sample of 52 B-R color measurements for 8 Centaurs, 22 Classicals, 13 Plutinos, 8 Scattered objects and 1 object with unidentified dynamical class. This dataset is the largest single and homogeneous published dataset to date [...]. A strong (color) correlation with mean excitation velocity points toward a space weathering/impact origin for the color diversity. However, thorough modeling of the collisional/dynamical environment in the Edgeworth-Kuiper belt needs to be done in order to confirm this scenario. We found also that the Classical TNOs consist in the superposition of two distinct populations: the dynamically Cold Classical TNOs (red colors, low i, small sizes) and the dynamically Hot Classical TNOs (diverse colors, moderate and high i, larger sizes). [...] Our specific observation strategy [...] permitted us to highlight a few objects suspected to have true compositional and/or texture variation on their surfaces. These are 1998 HK151, 1999 DF9, 1999 OY3, 2000 GP183, 2000 OK67, and 2001 KA77 and should be prime targets for further observations [...]. Our survey has also highlighted 1998 SN165 whose colors and dynamical properties puts it in a new dynamical class distinct from the Classicals, its previously assigned dynamical class.

Observations of the November 1995 Sun crossing of the Saturn's ring-plane made with the 3.6m CFH telescope, using the UHAO adaptive optics system, are presented here. We report the detection of four arcs located in the vicinity of the F ring. They can be seen one day later in HST images. The combination of both data sets gives accurate determinations of their orbits. Semi-major axes range from 140020 km to 140080 km, with a mean of 140060 +- 60 km. This is about 150 km smaller than previous estimates of the F ring radius from Voyager 1 and 2 data, but close to the orbit of another arc observed at the same epoch in HST images.