Dusty debris disks around pre- and main-sequence stars are potential signposts for the existence of planetesimals and exoplanets. Giant planet formation is therefore expected to play a key role in the evolution of the disk. This is indirectly confirmed by extant sub-millimeter near-infrared images of young protoplanetary and cool dusty debris disks around main sequence stars usually showing substantial spatial structures. A majority of recent discoveries of imaged giant planets have been obtained around young, early-type stars hosting a circumstellar disk. In this context, we have carried out a direct imaging program designed to maximize our chances of giant planet discovery and targeting twenty-two young, early-type stars. About half of them show indication of multi-belt architectures. Using the IRDIS dual-band imager and the IFS integral field spectrograph of SPHERE to acquire high-constrast coronagraphic differential near-infrared images, we have conducted a systematic search in the close environment of these young, dusty and early-type stars. We confirmed that companions detected around HIP 34276, HIP 101800 and HIP 117452 are stationary background sources and binary companions. The companion candidates around HIP 8832, HIP 16095 and HIP 95619 are determined as background contamination. For stars for which we infer the presence of debris belts, a theoretical minimum mass for planets required to clear the debris gaps can be calculated . The dynamical mass limit is at least $0.1 M_J$ and can exceed $1 M_J$. Direct imaging data is typically sensitive to planets down to $\sim 3.6 M_J$ at 1 $''$, and $1.7 M_J$ in the best case. These two limits tightly constrain the possible planetary systems present around each target. These systems will be probably detectable with the next generation of planet imagers.

Recent high-contrast imaging surveys, looking for planets in young, nearby systems showed evidence of a small number of giant planets at relatively large separation beyond typically 20 au where those surveys are the most sensitive. Access to smaller physical separations between 5 and 20 au is the next step for future planet imagers on 10 m telescopes and ELTs in order to bridge the gap with indirect techniques (radial velocity, transit, astrometry with Gaia). In that context, we recently proposed a new algorithm, Keplerian-Stacker, combining multiple observations acquired at different epochs and taking into account the orbital motion of a potential planet present in the images to boost the ultimate detection limit. We showed that this algorithm is able to find planets in time series of simulated images of SPHERE even when a planet remains undetected at one epoch. Here, we validate the K-Stacker algorithm performances on real SPHERE datasets, to demonstrate its resilience to instrumental speckles and the gain offered in terms of true detection. This will motivate future dedicated multi-epoch observation campaigns in high-contrast imaging to search for planets in emitted and reflected light. Results. We show that K-Stacker achieves high success rate when the SNR of the planet in the stacked image reaches 7. The improvement of the SNR ratio goes as the square root of the total exposure time. During the blind test and the redetection of HD 95086 b, and betaPic b, we highlight the ability of K-Stacker to find orbital solutions consistent with the ones derived by the state of the art MCMC orbital fitting techniques, confirming that in addition to the detection gain, K-Stacker offers the opportunity to characterize the most probable orbital solutions of the exoplanets recovered at low signal to noise.

Our objective is to estimate the C/O ratio in the atmosphere of beta Pictoris b and obtain an estimate of the dynamical mass of the planet, as well as to refine its orbital parameters using high-precision astrometry. We used the GRAVITY instrument with the four 8.2 m telescopes of the Very Large Telescope Interferometer to obtain K-band spectro-interferometric data on $\beta$ Pic b. We extracted a medium resolution (R=500) K-band spectrum of the planet and a high-precision astrometric position. We estimated the planetary C/O ratio using two different approaches (forward modeling and free retrieval) from two different codes (ExoREM and petitRADTRANS, respectively). Finally, we used a simplified model of two formation scenarios (gravitational collapse and core-accretion) to determine which can best explain the measured C/O ratio. Our new astrometry disfavors a circular orbit for $\beta$ Pic b ($e=0.15^{+0.05}_{-0.04}$). Combined with previous results and with Hipparcos/GAIA measurements, this astrometry points to a planet mass of $M = 12.7\pm{}2.2\,M_\mathrm{Jup}$. This value is compatible with the mass derived with the free-retrieval code petitRADTRANS using spectral data only. The forward modeling and free-retrieval approches yield very similar results regarding the atmosphere of beta Pic b. In particular, the C/O ratios derived with the two codes are identical ($0.43\pm{}0.05$ vs $0.43^{+0.04}_{-0.03}$). We argue that if the stellar C/O in $\beta$ Pic is Solar, then this combination of a very high mass and a low C/O ratio for the planet suggests a formation through core-accretion, with strong planetesimal enrichment.

Mixed-variable symplectic integrators are widely used in orbital dynamics. However, they have been developed for Solar system-type architectures, and can not handle evolving hierarchy, in particular in systems with two or more stellar components. Such configuration may have occurred in the history of HD 106906, a tight pair of F-type stars surrounded by a debris disk and a planetary-mass companion on a wide orbit. We present the new algorithm ODEA, based on the symplectic algorithm SWIFT HJS, that can model any system (binary,...) with unstable architecture. We study the peculiar system HD 106906 as a testcase for the code. We define and compute a criterion based on acceleration ratios to indicate when the initial hierarchy is not relevant anymore. A new hierarchy is then computed. The code is applied to study the two fly-bys that occurred on system HD 106906, recently evidenced by De Rosa & Kalas (2019), to determine if they could account for the wide orbit of the planet. Thousands of simulations have been performed to account for the uncertainty on the perturbers coordinates and velocities. The algorithm is able to handle any change of hierarchy, temporary or not. We used it to fully model HD 106906 encounters. The simulations confirm that the fly-bys could have stabilized the planet orbit, and show that it can account for the planet probable misalignment with respect to the disk plane as well as the disk morphology. However, that requires a small distance at closest approach (< 0.05 pc), and this configuration is not guaranteed. ODEA is a very good choice for the study of non-Solar type architecture. It can now adapt to an evolving hierarchy, and is thus suitable to study capture of planets and dust. Further observations of the perturbers, in particular their radial velocity, are required to conclude on the effects of the fly-by on system HD 106906.

High contrast imaging enables the determination of orbital parameters for substellar companions (planets, brown dwarfs) from the observed relative astrometry and the estimation of model and age-dependent masses from their observed magnitudes or spectra. Combining astrometric positions with radial velocity gives direct constraints on the orbit and on the dynamical masses of companions. A brown dwarf was discovered with the VLT/SPHERE instrument in 2017, which orbits at $\sim$ 11 au around HD 206893. Its mass was estimated between 12 and 50 $M_{Jup}$ from evolutionary models and its photometry. However, given the significant uncertainty on the age of the system and the peculiar spectrophotometric properties of the companion, this mass is not well constrained. We aim at constraining the orbit and dynamical mass of HD 206893 B. We combined radial velocity data obtained with HARPS spectra and astrometric data obtained with the high contrast imaging VLT/SPHERE and VLT/NaCo instruments, with a time baseline less than three years. We then combined those data with astrometry data obtained by Hipparcos and Gaia with a time baseline of 24 years. We used a MCMC approach to estimate the orbital parameters and dynamical mass of the brown dwarf from those data. We infer a period between 21 and 33{\deg} and an inclination in the range 20-41{\deg} from pole-on from HD 206893 B relative astrometry. The RV data show a significant RV drift over 1.6 yrs. We show that HD 206893 B cannot be the source of this observed RV drift as it would lead to a dynamical mass inconsistent with its photometry and spectra and with Hipparcos and Gaia data. An additional inner (semimajor axis in the range 1.4-2.6 au) and massive ($\sim$ 15 $M_{Jup}$) companion is needed to explain the RV drift, which is compatible with the available astrometric data of the star, as well as with the VLT/SPHERE and VLT/NaCo nondetection.

Planetary systems hold the imprint of the formation and of the evolution of planets especially at young ages, and in particular at the stage when the gas has dissipated leaving mostly secondary dust grains. The dynamical perturbation of planets in the dust distribution can be revealed with high-contrast imaging in a variety of structures. SPHERE, the high-contrast imaging device installed at the VLT, was designed to search for young giant planets in long period, but is also able to resolve fine details of planetary systems at the scale of astronomical units in the scattered-light regime. As a young and nearby star, NZ Lup was observed in the course of the SPHERE survey. A debris disk had been formerly identified with HST/NICMOS. We observed this system in the near-infrared with the camera in narrow and broad band filters and with the integral field spectrograph. High contrasts are achieved by the mean of pupil tracking combined with angular differential imaging algorithms. The high angular resolution provided by SPHERE allows us to reveal a new feature in the disk which is interpreted as a superimposition of two belts of planetesimals located at stellocentric distances of $\sim$85 and $\sim$115\,au, and with a mutual inclination of about 5$\degb$. Despite the very high inclination of the disk with respect to the line of sight, we conclude that the presence of a gap, that is, a void in the dust distribution between the belts, is likely. We discuss the implication of the existence of two belts and their relative inclination with respect to the presence of planets.

Context. The 51 Eridani system harbors a complex architecture with its primary star forming a hierarchical system with the binary GJ 3305AB at a projected separation of 2000 au, a giant planet orbiting the primary star at 13 au, and a low-mass debris disk around the primary star with possibly a cold component and a warm component inferred from the spectral energy distribution. Aims. We aim to better constrain the orbital parameters of the known giant planet. Methods. We monitored the system over three years from 2015 to 2018 with the VLT/SPHERE exoplanet imaging instrument. Results. We measure an orbital motion for the planet of ~130 mas with a slightly decreasing separation (~10 mas) and find a hint of curvature. This potential curvature is further supported at 3$\sigma$ significance when including literature GPI astrometry corrected for calibration systematics. Fits of the SPHERE and GPI data using three complementary approaches provide broadly similar results. The data suggest an orbital period of 32$^{+17}_{-9}$ yr (i.e. 12$^{+4}_{-2}$ au in semi-major axis), an inclination of 133$^{+14}_{-7}$ deg, an eccentricity of 0.45$^{+0.10}_{-0.15}$, and an argument of periastron passage of 87$^{+34}_{-30}$ deg [mod 180 deg]. The time at periastron passage and the longitude of node exhibit bimodal distributions because we do not detect yet if the planet is accelerating or decelerating along its orbit. Given the inclinations of the planet's orbit and of the stellar rotation axis (134-144 deg), we infer alignment or misalignment within 18 deg for the star-planet spin-orbit. Further astrometric monitoring in the next 3-4 years is required to confirm at a higher significance the curvature in the planet's motion, determine if the planet is accelerating or decelerating on its orbit, and further constrain its orbital parameters and the star-planet spin-orbit.

To date, infrared interferometry at best achieved contrast ratios of a few times $10^{-4}$ on bright targets. GRAVITY, with its dual-field mode, is now capable of high contrast observations, enabling the direct observation of exoplanets. We demonstrate the technique on HR8799, a young planetary system composed of four known giant exoplanets. We used the GRAVITY fringe tracker to lock the fringes on the central star, and integrated off-axis on the HR8799e planet situated at 390 mas from the star. Data reduction included post-processing to remove the flux leaking from the central star and to extract the coherent flux of the planet. The inferred K band spectrum of the planet has a spectral resolution of 500. We also derive the astrometric position of the planet relative to the star with a precision on the order of 100$\,\mu$as. The GRAVITY astrometric measurement disfavors perfectly coplanar stable orbital solutions. A small adjustment of a few degrees to the orbital inclination of HR 8799 e can resolve the tension, implying that the orbits are close to, but not strictly coplanar. The spectrum, with a signal-to-noise ratio of $\approx 5$ per spectral channel, is compatible with a late-type L brown dwarf. Using Exo-REM synthetic spectra, we derive a temperature of $1150\pm50$\,K and a surface gravity of $10^{4.3\pm0.3}\,$cm/s$^{2}$. This corresponds to a radius of $1.17^{+0.13}_{-0.11}\,R_{\rm Jup}$ and a mass of $10^{+7}_{-4}\,M_{\rm Jup}$, which is an independent confirmation of mass estimates from evolutionary models. Our results demonstrate the power of interferometry for the direct detection and spectroscopic study of exoplanets at close angular separations from their stars.

With an orbital distance comparable to that of Saturn in the solar system, \bpic b is the closest (semi-major axis $\simeq$\,9\,au) exoplanet that has been imaged to orbit a star. Thus it offers unique opportunities for detailed studies of its orbital, physical, and atmospheric properties, and of disk-planet interactions. With the exception of the discovery observations in 2003 with NaCo at the Very Large Telescope (VLT), all following astrometric measurements relative to \bpic have been obtained in the southwestern part of the orbit, which severely limits the determination of the planet's orbital parameters. We aimed at further constraining \bpic b orbital properties using more data, and, in particular, data taken in the northeastern part of the orbit. We used SPHERE at the VLT to precisely monitor the orbital motion of beta \bpic b since first light of the instrument in 2014. We were able to monitor the planet until November 2016, when its angular separation became too small (125 mas, i.e., 1.6\,au) and prevented further detection. We redetected \bpic b on the northeast side of the disk at a separation of 139\,mas and a PA of 30$^{\circ}$ in September 2018. The planetary orbit is now well constrained. With a semi-major axis (sma) of $a = 9.0 \pm 0.5$ au (1 $\sigma $), it definitely excludes previously reported possible long orbital periods, and excludes \bpic b as the origin of photometric variations that took place in 1981. We also refine the eccentricity and inclination of the planet. From an instrumental point of view, these data demonstrate that it is possible to detect, if they exist, young massive Jupiters that orbit at less than 2 au from a star that is 20 pc away.

Context: Numerous theoretical studies of the stellar dynamics of triple systems have been carried out, but fewer purely empirical studies that have addressed planetary orbits within these systems. Most of these empirical studies have been for coplanar orbits and with a limited number of orbital parameters. Aims: Our objective is to provide a more generalized empirical mapping of the regions of planetary stability in triples by considering both prograde and retrograde motion of planets and the outer star; investigating highly inclined orbits of the outer star; extending the parameters used to all relevant orbital elements of the triple's stars and expanding these elements and mass ratios to wider ranges that will accommodate recent and possibly future observational discoveries. Methods: Using N-body simulations, we integrated numerically the various four-body configurations over the parameter space, using a symplectic integrator designed specifically for the integration of hierarchical multiple stellar systems. The triples were then reduced to binaries and the integrations repeated to highlight the differences between these two types of system. Results: This established the regions of secular stability and resulted in 24 semi-empirical models describing the stability bounds for planets in each type of triple orbital configuration. The results were then compared with the observational extremes discovered to date to identify regions that may contain undiscovered planets.

The presence of a debris disc around the Gyr-old solar-type star $\zeta^2\,\mathrm{Reticuli}$ was suggested by the $\mathit{Spitzer}$ infrared excess detection. Follow-up observations with $\mathit{Herschel}$/PACS revealed a double-lobed feature, that displayed asymmetries both in brightness and position. Therefore, the disc was thought to be edge-on and significantly eccentric. Here we present ALMA/ACA observations in Band 6 and 7 which unambiguously reveal that these lobes show no common proper motion with $\zeta^2\,\mathrm{Reticuli}$. In these observations, no flux has been detected around $\zeta^2\,\mathrm{Reticuli}$ that exceeds the $3\sigma$ levels. We conclude that surface brightness upper limits of a debris disc around $\zeta^2\,\mathrm{Reticuli}$ are $5.7\,\mathrm{\mu Jy/arcsec^2}$ at 1.3 mm, and $26\,\mathrm{\mu Jy/arcsec^2}$ at 870 microns. Our results overall demonstrate the capability of the ALMA/ACA to follow-up $\mathit{Herschel}$ observations of debris discs and clarify the effects of background confusion.

Context. A low-mass brown dwarf has been recently imaged around HR 2562 (HD 50571), a star hosting a debris disk resolved in the far infrared. Interestingly, the companion location is compatible with an orbit coplanar with the disk and interior to the debris belt. This feature makes the system a valuable laboratory to analyze the formation of substellar companions in a circumstellar disk and potential disk-companion dynamical interactions. Aims. We aim to further characterize the orbital motion of HR 2562 B and its interactions with the host star debris disk. Methods. We performed a monitoring of the system over ~10 months in 2016 and 2017 with the VLT/SPHERE exoplanet imager. Results. We confirm that the companion is comoving with the star and detect for the first time an orbital motion at high significance, with a current orbital motion projected in the plane of the sky of 25 mas (~0.85 au) per year. No orbital curvature is seen in the measurements. An orbital fit of the SPHERE and literature astrometry of the companion without priors on the orbital plane clearly indicates that its orbit is (quasi-)coplanar with the disk. To further constrain the other orbital parameters, we used empirical laws for a companion chaotic zone validated by N-body simulations to test the orbital solutions that are compatible with the estimated disk cavity size. Non-zero eccentricities (>0.15) are allowed for orbital periods shorter than 100 yr, while only moderate eccentricities up to ~0.3 for orbital periods longer than 200 yr are compatible with the disk observations. A comparison of synthetic Herschel images to the real data does not allow us to constrain the upper eccentricity of the companion.

The G-type star GJ504A is known to host a 3 to 35 MJup companion whose temperature, mass, and projected separation all contribute to make it a test case for the planet formation theories and for atmospheric models of giant planets and light brown dwarfs. We collected data from the CHARA interferometer, SOPHIE spectrograph, and VLT/SPHERE high contrast imager to revisit the properties of the system. We measure a radius of 1.35+/- 0.04Rsun for GJ504A which yields isochronal ages of 21+/-2Myr or 4.0+/-1.8Gyr for the system and line-of-sight stellar rotation axis inclination of $162.4_{-4.3}^{+3.8}$ degrees or $18.6_{-3.8}^{+4.3}$ degrees. We re-detect the companion in the Y2, Y3, J3, H2, and K1 dual band SPHERE images. The complete 1-4 $\mu$m SED shape of GJ504b is best reproduced by T8-T9.5 objects with intermediate ages ($\leq1.5$Gyr), and/or unusual dusty atmospheres and/or super-solar metallicities. All six atmospheric models used yield $\mathrm{T_{eff}=550 \pm 50}$K for GJ504b and point toward a low surface gravity (3.5-4.0 dex). The accuracy on the metallicity value is limited by model-to-model systematics. It is not degenerate with the C/O ratio. We derive $\mathrm{log\:L/L_{\odot}=-6.15\pm0.15}$ dex for the companion compatible with masses of $\mathrm{M=1.3^{+0.6}_{-0.3}M_{Jup}}$ and $\mathrm{M=23^{+10}_{-9} M_{Jup}}$ for the young and old age ranges, respectively. The semi-major axis (sma) is above 27.8 au and the eccentricity lower than 0.55. The posterior on GJ~504b's orbital inclination suggests a misalignment with GJ~504A rotation axis. We combine the radial velocity and multi-epoch imaging data to exclude additional objects (90\% prob.) more massive than 2.5 and 30 $\mathrm{M_{Jup}}$ with sma in the range 0.01-80 au for the young and old system ages, respectively. The companion is in the envelope of the population of planets synthetized with our core-accretion model.

Aims: We aim to characterize the orbital and atmospheric properties of PDS 70 b, which was first identified on May 2015 in the course of the SHINE survey with SPHERE, the extreme adaptive-optics instrument at the VLT. Methods: We obtained new deep SPHERE/IRDIS imaging and SPHERE/IFS spectroscopic observations of PDS 70 b. The astrometric baseline now covers 6 years which allows us to perform an orbital analysis. For the first time, we present spectrophotometry of the young planet which covers almost the entire near-infrared range (0.96 to 3.8 micrometer). We use different atmospheric models covering a large parameter space in temperature, log(g), chemical composition, and cloud properties to characterize the properties of the atmosphere of PDS 70 b. Results: PDS 70 b is most likely orbiting the star on a circular and disk coplanar orbit at ~22 au inside the gap of the disk. We find a range of models that can describe the spectrophotometric data reasonably well in the temperature range between 1000-1600 K and log(g) no larger than 3.5 dex. The planet radius covers a relatively large range between 1.4 and 3.7 R_jupiter with the larger radii being higher than expected from planet evolution models for the age of the planet of 5.4 Myr. Conclusions: This study provides a comprehensive dataset on the orbital motion of PDS 70 b, indicating a circular orbit and a motion coplanar with the disk. The first detailed spectral energy distribution of PDS 70 b indicates a temperature typical for young giant planets. The detailed atmospheric analysis indicates that a circumplanetary disk may contribute to the total planet flux.

Evolutionary models are widely used to infer the mass of stars, brown dwarfs, and giant planets. Their predictions are thought to be less reliable at young ages ($<$ 200 Myr) and in the low-mass regime ($\mathrm{<1~M_{\odot}}$). GJ 2060 AB and TWA 22 AB are two rare astrometric M-dwarf binaries respectively members of the AB Doradus and Beta Pictoris moving groups. As their dynamical mass can be measured within a few years, they can be used to calibrate the evolutionary tracks and set new constraints on the age of young moving groups. We find a total mass of $\mathrm{0.18\pm 0.02~M_\odot}$ for TWA 22. That mass is in good agreement with model predictions at the age of the Beta Pic moving group. We obtain a total mass of $\mathrm{1.09 \pm 0.10~M_{\odot}}$ for GJ 2060. We estimate a spectral type of M$1\pm0.5$, $\mathrm{L/L_{\odot}=-1.20\pm0.05}$ dex, and $\mathrm{T_{eff}=3700\pm100}$ K for GJ 2060 A. The B component is a M$3\pm0.5$ dwarf with $\mathrm{L/L_{\odot}=-1.63\pm0.05}$ dex and $\mathrm{T_{eff}=3400\pm100}$ K. The dynamical mass of GJ 2060 AB is inconsistent with the most recent models predictions (BCAH15, PARSEC) for an ABDor age in the range 50-150 Myr. It is 10 to 20\% (1-2 sigma, depending on the assumed age) above the models predictions, corresponding to an underestimation of $0.10$ to $0.20~\mathrm{M_\odot}$. Coevality suggests a young age for the system ($\sim$ 50 Myr) according to most evolutionary models. TWA 22 validates the predictions of recent evolutionary tracks at $\sim$20 Myr. On the other hand, we evidence a 1-2 sigma mismatch between the predicted and observed mass of GJ 2060 AB. This slight departure may indicate that one of the star hosts a tight companion. Alternatively, this would confirm the models tendency to underestimate the mass of young low-mass stars.

HD95086 (A8V, 17Myr) hosts a rare planetary system for which a multi-belt debris disk and a giant planet of 4-5MJup have been directly imaged. Our study aims to characterize the physical and orbital properties of HD95086b, search for additional planets at short and wide orbits and image the cold outer debris belt in scattered light. We used HARPS at the ESO 3.6m telescope to monitor the radial velocity of HD95086 over 2 years and investigate the existence of giant planets at less than 3au orbital distance. With the IRDIS dual-band imager and the IFS integral field spectrograph of SPHERE at VLT, we imaged the faint circumstellar environment beyond 10au at six epochs between 2015 and 2017. We do not detect additional giant planets around HD95086. We identified the nature (bound companion or background contaminant) of all point-like sources detected in the IRDIS field of view. None of them correspond to the ones recently discovered near the edge of the cold outer belt by ALMA. HD95086b is resolved for the first time in J-band with IFS. Its near-infrared spectral energy distribution is well fitted by a few dusty and/or young L7-L9 dwarf spectral templates. The extremely red 1-4um spectral distribution is typical of low-gravity objects at the L/T spectral type transition. The planet's orbital motion is resolved between January 2015 and May 2017. Together with past NaCo measurements properly re-calibrated, our orbital fitting solutions favor a retrograde low to moderate-eccentricity orbit e=0.2 (0.0 to 0.5), with a semi-major axis 52au corresponding to orbital periods of 288$ yrs and an inclination that peaks at i = 141deg, which is compatible with a planet-disk coplanar configuration. Finally, we report the detection in polarimetric differential imaging of the cold outer debris belt between 100 and 300au, consistent in radial extent with recent ALMA 1.3mm resolved observations.

The angle between the spin of a star and its planets' orbital planes traces the history of the planetary system. Exoplanets orbiting close to cool stars are expected to be on circular, aligned orbits because of strong tidal interactions with the stellar convective envelope. Spin-orbit alignment can be measured when the planet transits its star, but such ground-based spectroscopic measurements are challenging for cool, slowly-rotating stars. Here we report the characterization of a planet three-dimensional trajectory around an M dwarf star, derived by mapping the spectrum of the stellar photosphere along the chord transited by the planet. We find that the eccentric orbit of the Neptune-mass exoplanet GJ 436b is nearly perpendicular to the stellar equator. Both eccentricity and misalignment, surprising around a cool star, can result from dynamical interactions (via Kozai migration) with a yet-undetected outer companion. This inward migration of GJ 436b could have triggered the atmospheric escape that now sustains its giant exosphere. Eccentric, misaligned exoplanets orbiting close to cool stars might thus hint at the presence of unseen perturbers and illustrate the diversity of orbital architectures seen in exoplanetary systems.

A giant planet has been recently resolved at a projected distance of 730 au from the tight pair of young ($\sim$ 13 Myr) intermediate-mass stars HD 106906AB in the Lower Centaurus Crux (LCC) group. The stars are surrounded by a debris disk which displays a ring-like morphology and strong asymmetries at multiple scales. We aim at studying the likelihood of a scenario where the planet formed closer to the stars in the disk, underwent inward disk-induced migration, and got scattered away by the binary star before being stabilized by a close encounter (fly-by). We performed semi-analytical calculations and numerical simulations (Swift_HJS package) to model the interactions between the planet and the two stars. We accounted for the migration as a simple force. We studied the LCC kinematics to set constraints on the local density of stars, and therefore on the fly-by likelihood. We performed N-body simulations to determine the effects of the planet trajectories (ejection and secular effects) onto the disk morphology. The combination of the migration and mean-motion resonances with the binary star (often 1:6) can eject the planet. Nonetheless, we estimate that the fly-by hypothesis decreases the scenario probability to less than $10^{-7}$ for a derived local density of stars of 0.11 stars/pc$^{3}$. We show that the concomitant effect of the planet and stars trajectories induce spiral-features in the disk which may correspond to the observed asymmetries. Moreover, the present disk shape suggests that the planet is on an eccentric orbit. The scenario we explored is a natural hypothesis if the planet formed within a disk. Conversely, its low probability of occurrence and the fact that HD 106906 b shares some characteristics with other systems in Sco-Cen (e.g. HIP 78530, in terms of mass ratio and separation) may indicate an alternative formation pathway for those objects.

The young and nearby star \beta\ Pictoris (\beta\ Pic) is surrounded by a debris disk composed of dust and gas known to host a myriad evaporating exocomets, planetesimals and at least one planet. At an edge-on inclination, as seen from Earth, this system is ideal for debris disk studies providing an excellent opportunity to use absorption spectroscopy to study the planet forming environment. Using the Cosmic Origins Spectrograph (COS) instrument on the Hubble Space Telescope (HST) we observe the most abundant element in the disk, hydrogen, through the HI Lyman \alpha\ (Ly-\alpha\) line. We present a new technique to decrease the contamination of the Ly-\alpha\ line by geocoronal airglow in COS spectra. This Airglow Virtual Motion (AVM) technique allows us to shift the Ly-\alpha\ line of the astrophysical target away from the contaminating airglow emission revealing more of the astrophysical line profile. The column density of hydrogen in the \beta\ Pic stable gas disk at the stellar radial velocity is measured to be $\log(N_{\mathrm{H}}/1 \mathrm{cm}^2) \ll 18.5$. The Ly-\alpha\ emission line profile is found to be asymmetric and we propose that this is caused by HI falling in towards the star with a bulk radial velocity of $41\pm6$ km/s relative to \beta\ Pic and a column density of $\log(N_{\mathrm{H}}/1 \mathrm{cm}^2) = 18.6\pm0.1$. The high column density of hydrogen relative to the hydrogen content of CI chondrite meteorites indicates that the bulk of the hydrogen gas does not come from the dust in the disk. This column density reveals a hydrogen abundance much lower than solar, which excludes the possibility that the detected hydrogen could be a remnant of the protoplanetary disk or gas expelled by the star. We hypothesise that the hydrogen gas observed falling towards the star arises from the dissociation of water originating from evaporating exocomets.

In 2015, we initiated a survey of Scorpius-Centaurus A-F stars that are predicted to host warm-inner and cold-outer belts of debris similar to the case of the system HR~8799. The survey aims to resolve the disks and detect planets responsible for the disk morphology. In this paper, we study the F-type star HIP~67497 and present a first-order modelisation of the disk in order to derive its main properties. We used the near-infrared integral field spectrograph (IFS) and dual-band imager IRDIS of VLT/SPHERE to obtain angular-differential imaging observations of the circumstellar environnement of HIP~67497. We removed the stellar halo with PCA and TLOCI algorithms. We modeled the disk emission with the GRaTeR code. We resolve a ring-like structure that extends up to $\sim$450 mas ($\sim$50 au) from the star in the IRDIS and IFS data. It is best reproduced by models of a non-eccentric ring with an inclination of $80\pm1^{\circ}$, a position angle of $-93\pm1^{\circ}$, and a semi-major axis of $59\pm3$ au. We also detect an additional, but fainter, arc-like structure with a larger extension (0.65 arcsec) South of the ring that we model as a second belt of debris at $\sim$130 au. We detect 10 candidate companions at separations $\geq$1''. We estimate the mass of putative perturbers responsible for the disk morphology and compare it to our detection limits. Additional data are needed to find those perturbers, and to relate our images to large-scale structures seen with HST/STIS.

Statistical analysis of the orbits of distant Kuiper Belt Objects (KBOs) have led to suggest that an additional planet should reside in the Solar System. According to recent models, the secular action of this body should cause orbital alignment of the KBOs. It was recently claimed that the KBOs concerned by this dynamics are presumably trapped in mean motion resonances with the suspected planet. I reinvestigate here the secular model underlying this idea. The original analysis was done expanding and truncating the secular Hamiltonian. I show that this is inappropriate here, as the series expansion is not convergent. I present a study based on numerical computation of the Hamiltonian with no expansion. I show in phase-space diagrams the existence of apsidally anti-aligned, high eccentricity libration islands that were not present in the original modelling, but that match numerical simulations. These island were claimed to correspond to bodies trapped in mean-motion resonances with the hypothetical planet, and match the characteristics of the distant KBOs observed. My main result is that regular secular dynamics can account for the anti-aligned particles itself as well as mean-motion resonances. I also perform a semi-analytical study of resonant motion and show that some resonance are actually capable of producing the same libration islands. I discuss then the relative importance of both mechanisms.

We aim to unveil the observational imprint of physical mechanisms that govern planetary formation in the young, multiple system GG Tau A. We present ALMA observations of $^{12}$CO and $^{13}$CO 3-2 and 0.9 mm continuum emission with 0.35" resolution. The $^{12}$CO 3-2 emission, found within the cavity of the circumternary dust ring (at radius $< 180$ AU) where no $^{13}$CO emission is detected, confirms the presence of CO gas near the circumstellar disk of GG Tau Aa. The outer disk and the recently detected hot spot lying at the outer edge of the dust ring are mapped both in $^{12}$CO and $^{13}$CO. The gas emission in the outer disk can be radially decomposed as a series of slightly overlapping Gaussian rings, suggesting the presence of unresolved gaps or dips. The dip closest to the disk center lies at a radius very close to the hot spot location at $\sim250-260$~AU. The CO excitation conditions indicate that the outer disk remains in the shadow of the ring. The hot spot probably results from local heating processes. The two latter points reinforce the hypothesis that the hot spot is created by an embedded proto-planet shepherding the outer disk.

Imaging companions to main-sequence stars often allows to detect a projected orbital motion. MCMC has become very popular in for fitting their orbits. Some of these companions appear to move on very eccentric, possibly unbound orbits. This is the case for the exoplanet Fomalhaut b and the brown dwarf companion PZ Tel B. For such orbits, standard MCMC codes assuming only bound orbits may be inappropriate. We develop a new MCMC implementation able to handle bound and unbound orbits as well in a continuous manner, and we apply it to the cases of Fomalhaut b and PZ Tel B. This code is based on universal Keplerian variables and Stumpff functions formalism. We present two versions of this code, the second one using a different set of angular variables designed to avoid degeneracies arising when the projected orbital motion is quasi-radial, as it is the case for PZ Tel B. We also present additional observations of PZ Tel B. The code is applied to Fomalhaut b and PZ Tel B. Concerning Fomalhaut b, we confirm previous results, but we show that open orbital solutions are also possible. The eccentricity distribution nevertheless peaks around ~0.9 in the bound regime. We present a first successful orbital fit of PZ Tel B, showing in particular that the eccentricity distribution presents a sharp peak very close to e=1, meaning a quasi-parabolic orbit. It was recently suggested that unseen inner companions may lead orbital fitting algorithms to artificially give high eccentricities. We show that this caveat is unlikely to apply to Fomalhaut b. Concerning PZ Tel B, an inner ~12 MJup companion would mimic a e=1 orbit despite a real eccentricity around 0.7, but a dynamical analysis reveals that such a system would not be stable. We conclude that our orbital fit is robust.

GJ758 B is a brown dwarf companion to a nearby (15.76 pc) solar-type, metal-rich (M/H = +0.2 dex) main-sequence star (G9V) that was discovered with Subaru/HiCIAO in 2009. From previous studies, it has drawn attention as being the coldest (~600K) companion ever directly imaged around a neighboring star. We present new high-contrast data obtained during the commissioning of the SPHERE instrument at the VLT. The data was obtained in Y-, J-, H-, and Ks-bands with the dual-band imaging (DBI) mode of IRDIS, providing a broad coverage of the full near-infrared (near-IR) range at higher contrast and better spectral sampling than previously reported. In this new set of high-quality data, we report the re-detection of the companion, as well as the first detection of a new candidate closer-in to the star. We use the new 8 photometric points for an extended comparison of GJ758 B with empirical objects and 4 families of atmospheric models. From comparison to empirical object, we estimate a T8 spectral type, but none of the comparison object can accurately represent the observed near-IR fluxes of GJ758 B. From comparison to atmospheric models, we attribute a Teff = 600K $\pm$ 100K, but we find that no atmospheric model can adequately fit all the fluxes of GJ758 B. The photometry of the new candidate companion is broadly consistent with L-type objects, but a second epoch with improved photometry is necessary to clarify its status. The new astrometry of GJ758 B shows a significant proper motion since the last epoch. We use this result to improve the determination of the orbital characteristics using two fitting approaches, Least-Square Monte Carlo and Markov Chain Monte Carlo. Finally, we analyze the sensitivity of our data to additional closer-in companions and reject the possibility of other massive brown dwarf companions down to 4-5 AU. [abridged]

HD~106906AB is so far the only young binary system around which a planet has been imaged and a debris disk evidenced thanks to a strong IR excess. As such, it represents a unique opportunity to study the dynamics of young planetary systems. We aim at further investigating the close (tens of au scales) environment of the HD~106906AB system. We used the extreme AO fed, high contrast imager SPHERE recently installed on the VLT to observe HD~106906. Both the IRDIS imager and the Integral Field Spectrometer were used. We discovered a very inclined, ring-like disk at a distance of 65~au from the star. The disk shows a strong brightness asymmetry with respect to its semi-major axis. It shows a smooth outer edge, compatible with ejection of small grains by the stellar radiation pressure. We show furthermore that the planet's projected position is significantly above the disk's PA. Given the determined disk inclination, it is not excluded though that the planet could still orbit within the disk plane if at a large separation (2000--3000 au). We identified several additional point sources in the SPHERE/IRDIS field-of-view, that appear to be background objects. We compare this system with other debris disks sharing similarities, and we briefly discuss the present results in the framework of dynamical evolution.

The eccentric shape of the debris disk observed around Fomalhaut was first attributed to Fom b, a companion detected near the belt inner-edge, but new constraints on its orbit revealed that it is belt-crossing, highly eccentric $(e \sim 0.6-0.9)$, and can hardly account for the shape of the belt. The best scenario to explain this paradox is that there is another massive body in this system, Fom c, which drives the debris disk shape. The resulting planetary system is highly unstable, which hints at a dynamical scenario involving a recent scattering of Fom b on its current orbit, potentially with the putative Fom c. Our goal is to give insights on the probability for Fom b to have been set on its highly eccentric orbit by a close-encounter with the putative Fom c. We aim to study in particular the part played by mean-motion resonances with Fom c, which could have brought Fom b sufficiently close to Fom c for it to be scattered on its current orbit, but also delay this scattering event. Using N-body simulations, we found that the generation of orbits similar to that of Fom b, either in term of dimensions or orientation, is a robust process involving a scattering event and a further secular evolution of inner material with an eccentric massive body such as the putative Fom c. We found in particular that mean-motion resonances can delay scattering events, and thus the production of Fom b-like orbits, on timescales comparable to the age of the system, thus explaining the witnessing of an unstable configuration. We conclude that Fom b probably originated from an inner resonance with Fom c, which is at least Neptune-Saturn size, and was set on its current orbit by a scattering event with Fom c. Since Fom b could not have formed from material in resonance, our scenario also hints at former migration processes in this planetary system.

The intermediate-mass star Beta Pictoris is known to be surrounded by a structured edge-on debris disk within which a gas giant planet was discovered orbiting at 8-10 AU. The physical properties of Beta Pic b were previously inferred from broad and narrow-band 0.9-4.8 microns photometry. We used commissioning data of the Gemini Planet Imager (GPI) to obtain new astrometry and a low-resolution (R=35-39) J-band (1.12-1.35 microns) spectrum of the planet. We find that the planet has passed the quadrature. We constrain its semi-major axis to $\leq$ 10 AU (90 % prob.) with a peak at 8.9+0.4-0.6 AU. The joint fit of the planet astrometry and the most recent radial velocity measurements of the star yields a planet's dynamical mass $\leq$ 20 MJup (greater than 96 % prob.). The extracted spectrum of Beta Pic b is similar to those of young L1-1.5+1 dwarfs. We use the spectral type estimate to revise the planet luminosity to log(L/Lsun)=-3.90+-0.07. The 0.9-4.8 microns photometry and spectrum are reproduced for Teff=1650+-150 K and a log g lower than 4.7 dex by 12 grids of PHOENIX-based and LESIA atmospheric models. If we adopt the most recent system age estimate (21+-4 Myr), the bolometric luminosity and the constraints on the dynamical mass of Beta Pic b are only reproduced by warm- and hot-start tracks with initial entropies Si greater than 10.5 kB/baryon. Such initial conditions may result from an inefficient accretion shock and/or a planetesimal density at formation higher than in the classical core accretion model. Considering a younger age for the system or a conservative formation time for Beta Pic b does not change these conclusions.

We aim at unveiling the observational imprint of physical mechanisms that govern planetary formation in young, multiple systems. In particular, we investigate the impact of tidal truncation on the inner circumstellar disks. We observed the emblematic system GG Tau at high-angular resolution: a hierarchical quadruple system composed of low-mass T Tauri binary stars surrounded by a well-studied, massive circumbinary disk in Keplerian rotation. We used the near-IR 4-telescope combiner PIONIER on the VLTI and sparse-aperture-masking techniques on VLT/NaCo to probe this proto-planetary system at sub-AU scales. We report the discovery of a significant closure-phase signal in H and Ks bands that can be reproduced with an additional low-mass companion orbiting GG Tau Ab, at a (projected) separation rho = 31.7 +/- 0.2mas (4.4 au) and PA = 219.6 +/- 0.3deg. This finding offers a simple explanation for several key questions in this system, including the missing-stellar-mass problem and the asymmetry of continuum emission from the inner dust disks observed at millimeter wavelengths. Composed of now five co-eval stars with 0.02 <= Mstar <= 0.7 Msun, the quintuple system GG Tau has become an ideal test case to constrain stellar evolution models at young ages (few 10^6yr).

The source HD172555 is a young A7V star surrounded by a debris disk with a gaseous component. Here, we present the detection of variable absorption features detected simultaneously in the Ca II K and H doublet lines (at 3,933 and 3,968 Angstrom). We identified the presence of these absorption signatures at four different epochs in the 129 HARPS high-resolution spectra gathered between 2004 and 2011. These transient absorption features are most likely due to Falling Evaporating Bodies (FEBs, or exocomets) that produce absorbing gas observed transiting in front of the central star. We also detect a stable Ca II absorption component at the star's radial velocity. With no corresponding detection in the Na I line, the resulting very low upper limit for the NaI/CaII ratio suggests that this absorption is due to circumstellar gas.

Imaging of debris disks has found evidence for both eccentric and offset disks. One hypothesis is that these provide evidence for massive perturbers that sculpt the observed structures. One such disk was recently observed in the far-IR by the Herschel Space Observatory around $\zeta^2$ Ret. In contrast with previously reported systems, the disk is significantly eccentric, and the system is Gyr-old. We aim to investigate the long-term evolution of eccentric structures in debris disks caused by a perturber on an eccentric orbit. Both analytical predictions and numerical N-body simulations are used to investigate the observable structures that could be produced by eccentric perturbers. The long-term evolution of the disk geometry is examined, with particular application to the $\zeta^2$ Ret system. In addition, synthetic images of the disk are produced for comparison with Herschel observations. We show that an eccentric companion can produce both the observed offsets and eccentric disks. Such effects are not immediate and we characterise the timescale required for the disk to develop to an eccentric state. For the case of $\zeta^2$ Ret, we place limits on the mass and orbit of the companion required to produce the observations. Synthetic images show that the pattern observed around $\zeta^2$ Ret can be produced by an eccentric disk seen close to edge-on, and allow us to bring additional constraints on the disk parameters of our model (disk flux, extent). We determine that eccentric planets or stellar companions can induce long-lived eccentric structures in debris disks. Observations of such eccentric structures provide potential evidence of the presence of such a companion in a planetary system. We consider the example of $\zeta^2$ Ret, whose observed eccentric disk can be explained by a distant companion at tens of AU, on an eccentric orbit ($e_p\gtrsim 0.3$).

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.

[Abridged] Debris disks are extrasolar analogs to the solar system planetesimal belts. The star Fomalhaut harbors a cold debris belt at 140 AU as well as evidence of a warm dust component, which is suspected of being a bright analog to the solar system's zodiacal dust. Interferometric observations obtained with the VLTI and the KIN have identified near- and mid-infrared excesses attributed to hot and warm exozodiacal dust in the inner few AU of the star. We performed parametric modeling of the exozodiacal disk using the GRaTeR radiative transfer code to reproduce the interferometric data, complemented by mid- to far-infrared measurements. A detailed treatment of sublimation temperatures was introduced to explore the hot population at the sublimation rim. We then used an analytical approach to successively testing several source mechanisms. A good fit to the data is found by two distinct dust populations: (1) very small, hence unbound, hot dust grains confined in a narrow region at the sublimation rim of carbonaceous material; (2) bound grains at 2 AU that are protected from sublimation and have a higher mass despite their fainter flux level. We propose that the hot dust is produced by the release of small carbon grains following the disruption of aggregates that originate from the warm component. A mechanism, such as gas braking, is required to further confine the small grains for a long enough time. In situ dust production could hardly be ensured for the age of the star, so the observed amount of dust must be triggered by intense dynamical activity. Fomalhaut may be representative of exozodis that are currently being surveyed worldwide. We propose a framework for reconciling the hot exozodi phenomenon with theoretical constraints: the hot component of Fomalhaut is likely the tip of the iceberg since it could originate from a warm counterpart residing near the ice line.

We aim at measuring the stellar parameters of the two Chemically Peculiar components of the B9.5Vp HgMn + A2 Vm double-lined spectroscopic binary HD141556, whose period is 15.25 days. We combined historical radial velocity measurements with new spatially resolved astrometric observations from PIONIER/VLTI to reconstruct the three-dimensional orbit of the binary, and thus obtained the individual masses. We fit the available photometric points together with the flux ratios provided by interferometry to constrain the individual sizes, which we compared to predictions from evolutionary models.The individual masses of the components are $\Ma = 2.84 \pm 0.12\ \Msun$ and $\Mb = 1.94 \pm 0.09\ \Msun$. The dynamical distance is compatible with the Hipparcos parallax. We find linear stellar radii of $\Ra=2.85 \pm 0.15\ \Rsun$ and $\Rb=1.75 \pm 0.18\ \Rsun$. This result validates a posteriori the flux ratio used in previous detailed abundance studies. We determine a sub-solar initial metallicity $Z=0.012\pm0.003$ and an age of $(2.8\pm0.3)\times10^8\ $years. Our results imply that the primary rotates more slowly than its synchronous velocity, while the secondary is probably synchronous. We show that strong tidal coupling during the pre-main sequence evolution followed by a full decoupling at zero-age main sequence provides a plausible explanation for these very low rotation rates.

A gas giant planet has previously been directly seen orbiting at 8-10 AU within the debris disk of the ~12 Myr old star {\beta} Pictoris. The {\beta} Pictoris system offers the rare opportunity to study the physical and atmospheric properties of an exoplanet placed on a wide orbit and to establish its formation scenario. We obtained J (1.265 {\mu}m), H (1.66 {\mu}m), and M' (4.78 {\mu}m) band angular differential imaging of the system between 2011 and 2012. We detect the planetary companion in our four-epoch observations. We estimate J = 14.0 +- 0.3, H = 13.5 +- 0.2, and M' = 11.0 +- 0.3 mag. Our new astrometry consolidates previous semi-major axis (sma=8-10 AU) and excentricity (e <= 0.15) estimates of the planet. These constraints, and those derived from radial velocities of the star provides independent upper limits on the mass of {\beta} Pictoris b of 12 and 15.5 MJup for semi-major axis of 9 and 10 AU. The location of {\beta} Pictoris b in color-magnitude diagrams suggests it has spectroscopic properties similar to L0-L4 dwarfs. This enables to derive Log10(L/Lsun) = -3.87 +- 0.08 for the companion. The analysis with 7 PHOENIX-based atmospheric models reveals the planet has a dusty atmosphere with Teff = 1700 +- 100 K and log g = 4.0+- 0.5. "Hot-start" evolutionary models give a new mass of 10+3-2 MJup from Teff and 9+3-2 MJup from luminosity. Predictions of "cold-start" models are inconsistent with independent constraints on the planet mass. "Warm-start" models constrain the mass to M >= 6MJup and the initial entropies to values (Sinit >= 9.3Kb/baryon), intermediate between those considered for cold/hot-start models, but likely closer to those of hot-start models.

The close-in planet orbiting GJ 436 presents a puzzling orbital eccentricity considering its very short orbital period. Given the age of the system, this planet should have been tidally circularized a long time ago. Many attempts to explain this were proposed in recent years, either involving abnormally weak tides, or the perturbing action of a distant companion. We address here the latter issue based on Kozai migration. We propose that GJ 436b was formerly located further away from the star and that it underwent a migration induced by a massive, inclined perturber via Kozai mechanism. In this context, the perturbations by the companion trigger high amplitude variations to GJ 436b that cause tides to act at periastron. Then the orbit tidally shrinks to reach its present day location. We numerically integrate the 3-body system including tides and General Relativity correction. We first show that starting from the present-day location of GJ 436b inevitably leads to damping the Kozai oscillations and to rapidly circularizing the planet. Conversely, starting from 5-10 times further away allows the onset of Kozai cycles. The tides act in peak eccentricity phases and reduce the semi-major axis of the planet. The net result is an evolution characterized by two phases: a first one with Kozai cycles and a slowly shrinking semi-major axis, and a second one once the planet gets out of the Kozai resonance characterized by a more rapid decrease. The timescale of this process appears in most cases much longer than the standard circularization time of the planet by a factor larger than 50. This model can provide a solution to the eccentricity paradox of GJ 436b. Depending on the various orbital configurations, it can take several Gyrs to GJ 436b to achieve a full orbital decrease and circularization. According to this scenario, we could be witnessing today the second phase of the scenario where the semi-major axis is already reduced while the eccentricity is still significant. We then explore the parameter space and derive in which conditions this model can be realistic given the age of the system. This yields constraints on the characteristics of the putative companion.

In June 2010, we confirmed the existence of a giant planet in the disk of the young star Beta Pictoris, located between 8 AU and 15 AU from the star. This young planet offers the rare opportunity to monitor a large fraction of the orbit using the imaging technique over a reasonably short timescale. Using the NAOS-CONICA adaptive-optics instrument (NACO) at the Very Large Telescope (VLT), we obtained repeated follow-up images of the Bpic system in the Ks and L' filters at four new epochs in 2010 and 2011. Complementing these data with previous measurements, we conduct a homogeneous analysis, which covers more than eight yrs, to accurately monitor the Bpic b position relative to the star. On the basis of the evolution of the planet's relative position with time, we derive the best-fit orbital solutions for our measurements. More reliable results are found with a Markov-chain Monte Carlo approach. The solutions favor a low-eccentricity orbit e < 0.17, with semi-major axis in the range 8--9 AU corresponding to orbital periods of 17--21 yrs. Our solutions favor a highly inclined solution with a peak around i=88.5+-1.7 deg, and a longitude of ascending node tightly constrained at Omega = -147.5+-1.5 deg. These results indicate that the orbital plane of the planet is likely to be above the midplane of the main disk, and compatible with the warp component of the disk being tilted between 3.5 deg and 4.0 deg. This suggests that the planet plays a key role in the origin of the inner warped-disk morphology of the Bpic disk. Finally, these orbital parameters are consistent with the hypothesis that the planet is responsible for the transit-like event observed in November 1981, and also linked to the cometary activity observed in the Bpic system.

Context. The {\beta} Pictoris system with its debris disk and a massive giant planet orbiting at \simeq 9 AU represents an ideal laboratory to study giant planet formation and evolution as well as planet-disk interactions. {\beta} Pic b can also help testing brightness-mass relations at young ages. Other planets, yet undetected, may of course be present in the system. Aims. We aim at putting direct constrains on the mass of {\beta} Pic b and at searching for additional jovian planets on orbits closer than typically 2 AU. Methods. We use high precision Harps data collected over 8 years since 2003 to measure and analyse {\beta} Pic radial velocities. Results. We show that the true mass of {\beta} Pic b is less than 10, 12, 15.5, 20 and 25 MJup if orbiting respectively at 8, 9, 10, 11 and 12 AU. This is the first direct constraint on the mass of an imaged planet. The upper mass found is well in the range predicted by brightness-mass relations provided by current "hot start" models. We also exclude the presence of giant planets more massive than 2.5 MJup with periods less than 100 days (hot Jupiters), more massive than 9 MJup for periods in the range 100-500 days. In the 500-1000 day range, the detection limit is in the brown dwarf domain. Beyond the intrinsic interest for {\beta} Pic, these results show the possibilities of precise RV measurements of early type, rapidly rotating stars.

Context. High contrast imaging is a powerful technique to search for gas giant planets and brown dwarfs orbiting at separation larger than several AU. Around solar-type stars, giant planets are expected to form by core accretion or by gravitational instability, but since core accretion is increasingly difficult as the primary star becomes lighter, gravitational instability would be the a probable formation scenario for yet-to-be-found distant giant planets around a low-mass star. A systematic survey for such planets around M dwarfs would therefore provide a direct test of the efficiency of gravitational instability. Aims. We search for gas giant planets orbiting around late-type stars and brown dwarfs of the solar neighborhood. Methods. We obtained deep high resolution images of 16 targets with the adaptive optic system of VLT-NACO in the Lp band, using direct imaging and angular differential imaging. This is currently the largest and deepest survey for Jupiter-mass planets around Mdwarfs. We developed and used an integrated reduction and analysis pipeline to reduce the images and derive our 2D detection limits for each target. The typical contrast achieved is about 9 magnitudes at 0.5" and 11 magnitudes beyond 1". For each target we also determine the probability of detecting a planet of a given mass at a given separation in our images. Results. We derived accurate detection probabilities for planetary companions, taking into account orbital projection effects, with in average more than 50% probability to detect a 3MJup companion at 10AU and a 1.5MJup companion at 20AU, bringing strong constraints on the existence of Jupiter-mass planets around this sample of young M-dwarfs.

The dynamical features of the irregular satellites of the giant planets argue against an in-situ formation and are strongly suggestive of a capture origin. Since the last detailed investigations of their dynamics, the total number of satellites have doubled, increasing from 50 to 109, and almost tripled in the case of Saturn system. We have performed a new dynamical exploration of Saturn system to test whether the larger sample of bodies could improve our understanding of which dynamical features are primordial and which are the outcome of the secular evolution of the system. We have performed detailed N--Body simulations using the best orbital data available and analysed the frequencies of motion to search for resonances and other possible perturbing effects. We took advantage of the Hierarchical Jacobian Symplectic algorithm to include in the dynamical model of the system also the gravitational effects of the two outermost massive satellites, Titan and Iapetus. Our results suggest that Saturn's irregular satellites have been significantly altered and shaped by the gravitational perturbations of Jupiter, Titan, Iapetus and the Sun and by the collisional sweeping effect of Phoebe. In particular, the effects on the dynamical evolution of the system of the two massive satellites appear to be non-negligible. Jupiter perturbs the satellites through its direct gravitational pull and, indirectly, via the effects of the Great Inequality, i.e. its almost resonance with Saturn. Finally, by using the Hierarchical Clustering Method we found hints to the existence of collisional families and compared them with the available observational data.

Planets can form and survive in close binaries, although dynamical interactions with the secondary component can actually significantly impact the giant planet formation and evolution. Rare close binaries hosting giant planets offer therefore an ideal laboratory to explore the properties and the stability of such extreme planetary systems. In the course of our CFHT and VLT coronographic imaging survey dedicated to the search for faint companions of exoplanet host stars, a close (about 20 AU) secondary stellar companion to the exoplanet host HD196885 A was discovered. For more than 4 years, we have used the NaCo near-infrared adaptive optics instrument to monitor the astrometric position of HD196885 B relative to A. The system was observed at five different epochs from August 2005 to August 2009 and accurate relative positions were determined. Our observations fully reject the stationary background hypothesis for HD196885 B. The two components are found to be comoving. The orbital motion of HD196885 B is well resolved and the orbital curvature is even detected. From our imaging data combined with published radial velocity measurements, we refine the complete orbital parameters of the stellar component. We derive for the first time its orbital inclination and its accurate mass. We find also solutions for the inner giant planet HD196885 Ab compatible with previous independent radial velocity studies. Finally, we investigate the stability of the inner giant planet HD196885 Ab due to the binary companion proximity. Our dynamical simulations show that the system is currently and surprisingly more stable in a high mutual inclination configuration that falls in the Kozai resonance regime. If confirmed, this system would constitute one of the most compact non-coplanar systems known so far. It would raise several questions about its formation and stability

The Fourier-Kelvin Stellar Interferometer (FKSI) is a mission concept for a nulling interferometer for the near-to-mid-infrared spectral region. FKSI is conceived as a mid-sized strategic or Probe class mission. FKSI has been endorsed by the Exoplanet Community Forum 2008 as such a mission and has been costed to be within the expected budget. The current design of FKSI is a two-element nulling interferometer. The two telescopes, separated by 12.5 m, are precisely pointed (by small steering mirrors) on the target star. The two path lengths are accurately controlled to be the same to within a few nanometers. A phase shifter/beam combiner (Mach-Zehnder interferometer) produces an output beam consisting of the nulled sum of the target planet's light and the host star's light. When properly oriented, the starlight is nulled by a factor of 10^-4, and the planet light is undiminished. Accurate modeling of the signal is used to subtract the residual starlight, permitting the detection of planets much fainter than the host star. The current version of FKSI with 0.5-m apertures and waveband 3-8 microns has the following main capabilities: (1) detect exozodiacal emission levels to that of our own solar system (1 Solar System Zodi) around nearby F, G, and K, stars; (2) characterize spectroscopically the atmospheres of a large number of known non-transiting planets; (3) survey and characterize nearby stars for planets down to 2 Earth radii from just inside the habitable zone and inward. An enhanced version of FKSI with 1-m apertures separated by 20 m and cooled to 40 K, with science waveband 5-15 microns, allows for the detection and characterization of 2 Earth-radius super-Earths and smaller planets in the habitable zone around stars within about 30 pc.

With the development of high contrast imaging techniques and infrared detectors, vast efforts have been devoted during the past decade to detect and characterize lighter, cooler and closer companions to nearby stars, and ultimately image new planetary systems. Complementary to other observing techniques (radial velocity, transit, micro-lensing, pulsar-timing), this approach has opened a new astrophysical window to study the physical properties and the formation mechanisms of brown dwarfs and planets. I here will briefly present the observing challenge, the different observing techniques, strategies and samples of current exoplanet imaging searches that have been selected in the context of the LAOG-Planet Imaging Surveys. I will finally describe the most recent results that led to the discovery of giant planets probably formed like the ones of our solar system, offering exciting and attractive perspectives for the future generation of deep imaging instruments.

(abridged) In the frame of the search for extrasolar planets and brown dwarfs around early-type main-sequence stars, we present the results obtained on the early F-type star Theta Cygni. Elodie and Sophie at OHP were used to obtain the spectra. Our dedicated radial-velocity measurement method was used to monitor the star's radial velocities over five years. We also use complementary, high angular resolution and high-contrast images taken with PUEO at CFHT. We show that Theta Cygni radial velocities are quasi-periodically variable, with a ~150-day period. These variations are not due to the ~0.35-Msun stellar companion that we detected in imaging at more than 46 AU from the star. The absence of correlation between the bisector velocity span variations and the radial velocity variations for this 7 km/s vsini star, as well as other criteria indicate that the observed radial velocity variations are not due to stellar spots. The observed amplitude of the bisector velocity span variations also seems to rule out stellar pulsations. However, we observe a peak in the bisector velocity span periodogram at the same period as the one found in the radial velocity periodogram, which indicates a probable link between these radial velocity variations and the low amplitude lineshape variations which are of stellar origin. Long-period variations are not expected from this type of star to our knowledge. If a stellar origin (hence of new type) was to be confirmed for these long-period radial velocity variations, this would have several consequences on the search for planets around main-sequence stars, both in terms of observational strategy and data analysis. An alternative explanation for these variable radial velocities is the presence of at least one planet of a few Jupiter masses orbiting at less than 1 AU. (abridged)

The GJ 581 planetary system was already known to harbour three planets, including two super-Earths planets which straddle its habitable zone. We report here the detection of an additional planet -- GJ 581e -- with a minimum mass of 1.9 M_earth. With a period of 3.15 days, it is the innermost planet of the system and has a ~5% transit probability. We also correct our previous confusion of the orbital period of GJ 581d (the outermost planet) with a one-year alias, thanks to an extended time span and many more measurements. The revised period is 66.8 days, and locates the semi-major axis inside the habitable zone of the low mass star. The dynamical stability of the 4-planet system imposes an upper bound on the orbital plane inclination. The planets cannot be more massive than approximately 1.6 times their minimum mass.

Tight binaries discovered in young, nearby associations, with known distances, are ideal targets to provide dynamical mass measurements to test the physics of evolutionary models at young ages and very low masses. We report for the first time the binarity of TWA22, possible new dynamical calibrator for evolutionary models at young ages. Based on an accurate trigonometric distance (17.53 +- 0.21 pc) determination, we infer a total dynamical mass of 220 +- 21 MJup for the system. From the resolved near-infrared integral-field spectroscopy, we find an effective temperature Teff=2900+200-200 K for TWA22 A and Teff=2900+200-100 K for TWA22 B and surface gravities between 4.0 and 5.5 dex. From our photometry and a M6 +- 1 spectral type for both components, we find luminosities of log(L/Lsun)=-2.11 +- 0.13 dex and log(L/Lsun)=-2.30 +- 0.16 dex for TWA22 A and B respectively. By comparing these parameters with evolutionary models, we question the age and the multiplicity of this system. We also discuss a possible underestimation of the mass predicted by evolutionary models for young stars close to the substellar boundary.

HD 141569, a triple star system, has been intensively observed and studied for its massive debris disk. It was rather regarded as a gravitationally bound triple system but recent measurements of the HD 141569A radial velocity seem to invalidate this hypothesis. The flyby scenario has therefore to be investigated to test its compatibility with the observations. We present a study of the flyby scenario for the HD141569 system, by considering 3 variants: a sole flyby, a flyby associated with one planet and a flyby with two planets. We use analytical calculations and perform N-body numerical simulations of the flyby encounter. The binary orbit is found to be almost fixed by the observational constraint on a edge-on plane with respect to the observers. If the binary has had an influence on the disk structure, it should have a passing time at the periapsis between 5000 and 8000 years ago and a distance at periapsis between 600 and 900 AU. The best scenario for reproducing the disk morphology is a flyby with only 1 planet. For a 2 Mj (resp. 8 Mj) planet, its eccentricity must be around 0.2 (resp. below 0.1). In the two cases, its apoapsis is about 130 AU. Although the global disk shape is reasonably well reproduced, some features cannot be explain by the present model and the likehood of the flyby event remains an issue. Dynamically speaking, HD 141569 is still a puzzling system.

Aims: In the frame of the search for extrasolar planets and brown dwarfs around early-type stars, we present the results obtained for the F-type main-sequence star HD 60532 (F6V) with HARPS. Methods: Using 147 spectra obtained with HARPS at La Silla on a time baseline of two years, we study the radial velocities of this star. Results: HD 60532 radial velocities are periodically variable, and the variations have a Keplerian origin. This star is surrounded by a planetary system of two planets with minimum masses of 1 and 2.5 Mjup and orbital separations of 0.76 and 1.58 AU respectively. We also detect high-frequency, low-amplitude (10 m/s peak-to-peak) pulsations. Dynamical studies of the system point toward a possible 3:1 mean-motion resonance which should be confirmed within the next decade.

We present a thorough study of the impact of a migrating planet on a planetesimal disk, by exploring a broad range of masses and eccentricities for the planet. We discuss the sensitivity of the structures generated in debris disks to the basic planet parameters. We perform many N-body numerical simulations, using the symplectic integrator SWIFT, taking into account the gravitational influence of the star and the planet on massless test particles. A constant migration rate is assumed for the planet. The effect of planetary migration on the trapping of particles in mean motion resonances is found to be very sensitive to the initial eccentricity of the planet and of the planetesimals. A planetary eccentricity as low as 0.05 is enough to smear out all the resonant structures, except for the most massive planets. The planetesimals also initially have to be on orbits with a mean eccentricity of less than than 0.1 in order to keep the resonant clumps visible. This numerical work extends previous analytical studies and provides a collection of disk images that may help in interpreting the observations of structures in debris disks. Overall, it shows that stringent conditions must be fulfilled to obtain observable resonant structures in debris disks. Theoretical models of the origin of planetary migration will therefore have to explain how planetary systems remain in a suitable configuration to reproduce the observed structures.

We address the issue of the dynamical evolution of the Gliese 581 planetary system. It is crucial when considering the planets' habitability because the secular evolution of the orbits may regulate their climate, even in the case where the system is stable. We have numerically integrated the planetary system over 10^8 yrs, starting from the present fitted solution. In all cases, the system appears dynamically stable, even in close to pole-on configurations. Only a limited range of inclinations can be excluded. The climate on the planets is expected to be secularly stable, thus not precluding the development of life. Gl 581 remains the best candidate for a planetary system with planets that potentially bear primitive forms of life.

The presence of off-plane Ca II ions in the Beta Pictoris disk, and the non-detection of off-plane Na I atoms, can be explained as a consequence of the evaporation process of Falling Evaporating Bodies (FEBs). In the star-grazing regime, the FEBs are subject to inclination oscillations up to 30 - 40 degrees that causes most metallic species released by sublimation to move off plane The ions are be stopped at about 100 AU from the star. We show that collisions with a neutral medium can stop the ions. The required H I column density is reduced to 10^17 cm^-2, one order of magnitude below present detection limits. We also investigate the possibility that the ions are slowed down magnetically. While the sole action of a magnetic field of the order of 1 microGauss is not effective, the combined effect of magnetic and collisional deceleration processes lead to an additional lowering of the required H I column density.

The goal of this study is to obtain accurate estimates for the individual masses of the components of the tight binary system T Tau S in order to settle the ongoing debate on the nature of T Tau Sa, a so-called infrared companion. We take advantage of the fact that T Tau S belongs to a triple system composed of two hierarchical orbits to simultaneously analyze the motion of T Tau Sb in the rest frames of T Tau Sa and T Tau N. With this method, it is possible to pinpoint the location of the center of mass of T Tau S and, thereby, to determine individual masses for T Tau Sa and T Tau Sb with no prior assumption about the mass/flux ratio of the system. This improvement over previous studies of the system results in much better constraints on orbital parameters. We find individual masses of 2.73+/-0.31 Msun for T Tau Sa and of 0.61+/-0.17 Msun for T Tau Sb (in agreement with its early-M spectral type), including the uncertainty on the distance to the system. These are among the most precise estimates of the mass of any Pre-Main Sequence star, a remarkable result since this is the first system in which individual masses of T Tauri stars can be determined from astrometry only. This model-independent analysis confirms that T Tau Sa is an intermediate-mass star, presumably a very young Herbig Ae star, that may possess an almost edge-on disk.

We present the results of multi epochs imaging observations of the companion to the planetary host Gliese 86. Associated to radial velocity measurements, this study aimed at characterizing dynamically the orbital properties and the mass of this companion (here after Gliese 86 B), but also at investigating the possible history of this particular system. We used the adaptive optics instrument NACO at the ESO Very Large Telescope to obtain deep coronographic imaging in order to determine new photometric and astrometric measurements of Gliese 86 B. Part of the orbit is resolved. The photometry of Gliese B indicates colors compatible with a ~70 Jupiter mass brown dwarf or a white dwarf. Both types of objects allow to fit the available, still limited astrometric data. Besides, if we attribute the long term radial velocity residual drift observed for Gliese A to B, then the mass of the latter object is ~0.5 Msun. We analyse both astrometric and radial velocity data to propose first orbital parameters for Gliese B. Assuming Gliese B is a ~0.5 Msun white dwarf, we explore the constraints induced by this hypothesis and refine the parameters of the system.

We present the first comprehensive analysis of the AU Mic debris disk properties since the system was discovered by Kalas et al. (2004), and we explore whether the dynamical model, successful to reproduce the Beta Pic brightness profile could apply to AU Mic. We calculate the surface density profile of the AU Mic disk by performing the inversion of the near-IR and visible scattered light brightness profiles measured by Liu (2004a) and Krist et al. (2005), respectively. We discuss the grain properties by analysing the blue color of the disk in the visible (Krist et al. 2005) and by fitting the disk spectral energy distribution. We show that irrespective of the mean scattering asymmetry factor of the grains, most of the emission arises from an asymmetric, collisionally-dominated region that peaks close to the surface brightness break around 35 AU. The elementary scatterers at visible wavelengths are found to be sub-micronic, but the inferred size distribution underestimates the amount of large grains, resulting in too low sub-millimeter emissions compared to the observations. From our inversion procedure, we find that the V- to H-band scattering cross sections ratio increases outside 40 AU, in line with the observed color gradient of the disk. We show that a standard, solar-like stellar wind generates a pressure force onto the dust particles that behaves much like a radiation pressure force. With an assumed Mdot ~ 300 Mdot_sun, the wind pressure overcomes the radiation pressure and this effect is enhanced by the stellar flares. This explains the similarity between the Beta Pic and AU Mic brightness profiles. In both cases, the color gradient beyond 120 AU for Beta Pic and 35 AU for AU Mic, is believed to be a direct consequence of the dust dynamics.

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 report the detection of a possible subsolar mass companion to the triple young system TY CrA using the NACO instrument at the VLT UT4 during its commissioning. Assuming for TY CrA a distance similar to that of the close binary system HD 176386, the photometric spectral type of this fourth stellar component candidate is consistent with an ~M4 star. We discuss the dynamical stability of this possible quadruple system as well as the possible location of dusty particles inside or outside the system.

We present the results of exploratory HST observations to detect the emission lines of ions in the Beta Pic disk with the spectrograph slit placed off the star. Possible emission lines from FeII have been detected at $\sim 10^{-14}$erg cm-2 s-1 at 0.5 arcsec from Beta Pic, which would suggest a total FeII density of > 2.$10^{-2}$ cm-3 at 10 AU. This detection needs confirmation. If real, it requires a large production rate of gas and dust equivalent to the total disruption of ten bodies 30 kilometers in radius per year.