J. Milli

Asymmetric features in the protoplanetary disk MWC 758

Context. The study of dynamical processes in protoplanetary disks is essential to understand planet formation. In this context, transition disks are prime targets because they are at an advanced stage of disk clearing and may harbor direct signatures of disk evolution. Aims. We aim to derive new constraints on the structure of the transition disk MWC 758, to detect non-axisymmetric features and understand their origin. Methods. We obtained infrared polarized intensity observations of the protoplanetary disk MWC 758 with SPHERE/VLT at 1.04 m to resolve scattered light at a smaller inner working angle (0.093 00 ) and a higher angular resolution (0.027 00 ) than previously achieved. Results. We observe polarized scattered light within 0.53 00 (148 au) down to the inner working angle (26 au) and detect distinct nonaxisymmetric features but no fully depleted cavity. The two small-scale spiral features that were previously detected with HiCIAO are resolved more clearly, and new features are identified, including two that are located at previously inaccessible radii close to the star. We present a model based on the spiral density wave theory with two planetary companions in circular orbits. The best model requires a high disk aspect ratio (H=r 0.20 at the planet locations) to account for the large pitch angles which implies a very warm disk. Conclusions. Our observations reveal the complex morphology of the disk MWC 758. To understand the origin of the detected features, the combination of high-resolution observations in the submillimeter with ALMA and detailed modeling is needed.

Fundamental Limitations of High Contrast Imaging Set by Small Sample Statistics

In this paper, we review the impact of small sample statistics on detection thresholds and corresponding confidence levels (CLs) in high-contrast imaging at small angles. When looking close to the star, the number of resolution elements decreases rapidly toward small angles. This reduction of the number of degrees of freedom dramatically affects CLs and false alarm probabilities. Naively using the same ideal hypothesis and methods as for larger separations, which are well understood and commonly assume Gaussian noise, can yield up to one order of magnitude error in contrast estimations at fixed CL. The statistical penalty exponentially increases toward very small inner working angles. Even at 5-10 resolution elements from the star, false alarm probabilities can be significantly higher than expected. Here we present a rigorous statistical analysis that ensures robustness of the CL, but also imposes a substantial limitation on corresponding achievable detection limits (thus contrast) at small angles. This unavoidable fundamental statistical effect has a significant impact on current coronagraphic and future high-contrast imagers. Finally, the paper concludes with practical recommendations to account for small number statistics when computing the sensitivity to companions at small angles and when exploiting the results of direct imaging planet surveys.

The position of β Pictoris b position relative to the debris disk

Context. We detected in 2009 a giant, close-by planet orbiting β Pic, a young star surrounded by a disk that has been extensively studied for more than 20 years. We showed that if the planet were located on an inclined orbit, this could account for several peculiarities of the β Pictoris system. However, the available data did not permit us to measure the inclination of β Pic b with respect to the disk, and in particular to establish in which component of the disk – either the main, extended disk or the inner inclined component/disk – the planet was located. Comparison between the observed planet position and the disk orientation measured using previous imaging data was not an option because of potential biases in the measurements. Aims. Our aim is to measure precisely the planet location with respect to the dust disk using a single high-resolution image, and correcting for systematics or errors that degrade the precision of the disk and planet relative-position measurements. Methods. We gathered new NaCo data in the Ks band, with a set-up optimized to derive simultaneously the orientation(s) of the disk(s) and the planet projected position. Results. We show that the projected position of β Pic b is above the midplane of the main disk. With the current data and knowledge of the system, this implies that β Pic b cannot be located in the main disk. The data instead suggest that the planet is located in the inclined component.

First light of the VLT planet finder SPHERE III. New spectrophotometry and astrometry of the HR 8799 exoplanetary system

Context. The planetary system discovered around the young A-type HR 8799 provides a unique laboratory to: a) test planet formation theories; b) probe the diversity of system architectures at these separations, and c) perform comparative (exo)planetology. Aims. We present and exploit new near-infrared images and integral-field spectra of the four gas giants surrounding HR 8799 obtained with SPHERE, the new planet finder instrument at the Very Large Telescope, during the commissioning and science verification phase of the instrument (July–December 2014). With these new data, we contribute to completing the spectral energy distribution (SED) of these bodies in the 1.0–2.5 μm range. We also provide new astrometric data, in particular for planet e, to further constrain the orbits. Methods. We used the infrared dual-band imager and spectrograph (IRDIS) subsystem to obtain pupil-stabilized, dual-band H2H3 (1.593 μm, 1.667 μm), K1K2 (2.110 μm, 2.251 μm), and broadband J (1.245 μm) images of the four planets. IRDIS was operated in parallel with the integral field spectrograph (IFS) of SPHERE to collect low-resolution (R ~ 30), near-infrared (0.94–1.64 μm) spectra of the two innermost planets HR 8799 d and e. The data were reduced with dedicated algorithms, such as the Karhunen-Loeve image projection (KLIP), to reveal the planets. We used the so-called negative planets injection technique to extract their photometry, spectra, and measure their positions. We illustrate the astrometric performance of SPHERE through sample orbital fits compatible with SPHERE and literature data. Results. We demonstrated the ability of SPHERE to detect and characterize planets in this kind of systems, providing spectra and photometry of its components. The spectra improve upon the signal-to-noise ratio of previously obtained data and increase the spectral coverage down to the Y band. In addition, we provide the first detection of planet e in the J band. Astrometric positions for planets HR 8799 bcde are reported for the epochs of July, August, and December 2014. We measured the photometric values in J, H2H3, K1K2 bands for the four planets with a mean accuracy of 0.13 mag. We found upper limit constraints on the mass of a possible planet f of 3–7 MJup . Our new measurements are more consistent with the two inner planets d and e being in a 2d:1e or 3d:2e resonance. The spectra of HR 8799 d and e are well matched by those of L6-8 field dwarfs. However, the SEDs of these objects are redder than field L dwarfs longward of 1.6 μm.

Review of small-angle coronagraphic techniques in the wake of ground-based second-generation adaptive optics systems

Small-angle coronagraphy is technically and scientifically appealing because it enables the use of smaller telescopes, allows covering wider wavelength ranges, and potentially increases the yield and completeness of circumstellar environment – exoplanets and disks – detection and characterization campaigns. However, opening up this new parameter space is challenging. Here we will review the four posts of high contrast imaging and their intricate interactions at very small angles (within the first 4 resolution elements from the star). The four posts are: choice of coronagraph, optimized wavefront control, observing strategy, and post-processing methods. After detailing each of the four foundations, we will present the lessons learned from the 10+ years of operations of zeroth and first-generation adaptive optics systems. We will then tentatively show how informative the current integration of second-generation adaptive optics system is, and which lessons can already be drawn from this fresh experience. Then, we will review the current state of the art, by presenting world record contrasts obtained in the framework of technological demonstrations for space-based exoplanet imaging and characterization mission concepts. Finally, we will conclude by emphasizing the importance of the cross-breeding between techniques developed for both ground-based and space-based projects, which is relevant for future high contrast imaging instruments and facilities in space or on the ground.

First light of the VLT planet finder SPHERE I. Detection and characterization of the substellar companion GJ 758 B

GJ 758 B is a brown dwarf companion to a nearby (15.76%) 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 (similar to 600 K) companion ever directly imaged around a neighboring star. We present new high-contrast data obtained during the commissioning of the SPHERE instrument at the Very Large Telescope (VLT). The data was obtained in Y-, J-, H-, and K-s-bands with the dual-band imaging (DBI) mode of IRDIS, thus 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 eight photometric points for an extended comparison of GJ 758 B with empirical objects and four families of atmospheric models. From comparison to empirical object, we estimate a T8 spectral type, but none of the comparison objects can accurately represent the observed near-IR fluxes of GJ 758 B. From comparison to atmospheric models, we attribute a T-eff = 600 +/- 100 K, but we find that no atmospheric model can adequately fit all the fluxes of GJ 758 B. The lack of exploration of metal enrichment in model grids appears as a major limitation that prevents an accurate estimation of the companion physical parameters. 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 GJ 758 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-Squares Monte Carlo and Markov chain Monte Carlo. We confirm the high-eccentricity of the orbit (peak at 0.5), and find a most likely semi-major axis of 46.05 AU. We also use our imaging data, as well as archival radial velocity data, to reject the possibility that this is a false positive effect created by an unseen, closer-in, companion. 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.

A near-infrared interferometric survey of debris-disk stars. IV. An unbiased sample of 92 southern stars observed in H band with VLTI/PIONIER ⋆

Context. Detecting and characterizing circumstellar dust is a way to study the architecture and evolution of planetary systems. Cold dust in debris disks only traces the outer regions. Warm and hot exozodiacal dust needs to be studied in order to trace regions close to the habitable zone. Aims. We aim to determine the prevalence and to constrain the properties of hot exozodiacal dust around nearby main-sequence stars. Methods. We searched a magnitude-limited (H ≤ 5) sample of 92 stars for bright exozodiacal dust using our VLTI visitor instrument PIONIER in the H band. We derived statistics of the detection rate with respect to parameters, such as the stellar spectral type and age or the presence of a debris disk in the outer regions of the systems. We derived more robust statistics by combining our sample with the results from our CHARA/FLUOR survey in the K band. In addition, our spectrally dispersed data allowed us to put constraints on the emission mechanism and the dust properties in the detected systems. Results. We find an overall detection rate of bright exozodiacal dust in the H band of 11% (9 out of 85 targets) and three tentative detections. The detection rate decreases from early type to late type stars and increases with the age of the host star. We do not confirm the tentative correlation between the presence of cold and hot dust found in our earlier analysis of the FLUOR sample alone. Our spectrally dispersed data suggest that either the dust is extremely hot or the emission is dominated by the scattered light in most cases. The implications of our results for the target selection of future terrestrial planet-finding missions using direct imaging are discussed.

First light of the VLT planet finder SPHERE IV : Physical and chemical properties of the planets around HR8799

Context. The system of fourplanets discovered around the intermediate-mass star HR8799 offers a unique opportunity to test planet formation theories at large orbital radii and to probe the physics and chemistry at play in the atmospheres of self-luminous young (~30 Myr) planets. We recently obtained new photometry of the four planets and low-resolution (R ~ 30) spectra of HR8799 d and e with the SPHERE instrument (Paper III). Aims. In this paper (Paper IV), we aim to use these spectra and available photometry to determine how they compare to known objects, what the planet physical properties are, and how their atmospheres work. Methods. We compare the available spectra, photometry, and spectral energy distribution (SED) of the planets to field dwarfs and young companions. In addition, we use the extinction from corundum, silicate (enstatite and forsterite), or iron grains likely to form in the atmosphere of the planets to try to better understand empirically the peculiarity of their spectrophotometric properties. To conclude, we use three sets of atmospheric models (BT-SETTL14, Cloud-AE60, Exo-REM) to determine which ingredients are critically needed in the models to represent the SED of the objects, and to constrain their atmospheric parameters (T_(eff), log g, M/H). Results. We find that HR8799d and e properties are well reproduced by those of L6-L8 dusty dwarfs discovered in the field, among which some are candidate members of young nearby associations. No known object reproduces well the properties of planets b and c. Nevertheless, we find that the spectra and WISE photometry of peculiar and/or young early-T dwarfs reddened by submicron grains made of corundum, iron, enstatite, or forsterite successfully reproduce the SED of these planets. Our analysis confirms that only the Exo-REM models with thick clouds fit (within 2σ) the whole set of spectrophotometric datapoints available for HR8799 d and e for T_(eff) = 1200 K, log g in the range 3.0−4.5, and M/H = +0.5. The models still fail to reproduce the SED of HR8799c and b. The determination of the metallicity, log g, and cloud thickness are degenerate. Conclusions. Our empirical analysis and atmospheric modelling show that an enhanced content in dust and decreased CIA of H_2 is certainly responsible for the deviation of the properties of the planet with respect to field dwarfs. The analysis suggests in addition that HR8799c and b have later spectral types than the two other planets, and therefore could both have lower masses.

Fast-moving features in the debris disk around AU Microscopii

In the 1980s, excess infrared emission was discovered around main-sequence stars; subsequent direct-imaging observations revealed orbiting disks of cold dust to be the source. These ‘debris disks’ were thought to be by-products of planet formation because they often exhibited morphological and brightness asymmetries that may result from gravitational perturbation by planets. This was proved to be true for the β Pictoris system, in which the known planet generates an observable warp in the disk. The nearby, young, unusually active late-type star AU Microscopii hosts a well-studied edge-on debris disk; earlier observations in the visible and near-infrared found asymmetric localized structures in the form of intensity variations along the midplane of the disk beyond a distance of 20 astronomical units. Here we report high-contrast imaging that reveals a series of five large-scale features in the southeast side of the disk, at projected separations of 10–60 astronomical units, persisting over intervals of 1–4 years. All these features appear to move away from the star at projected speeds of 4–10 kilometres per second, suggesting highly eccentric or unbound trajectories if they are associated with physical entities. The origin, localization, morphology and rapid evolution of these features are difficult to reconcile with current theories.

A near-infrared interferometric survey of debris-disc stars. V. PIONIER search for variability

Context. Extended circumstellar emission has been detected within a few 100 milli-arcsec around ≳10% of nearby main sequence stars using near-infrared interferometry. Follow-up observations using other techniques, should they yield similar results or non-detections, can provide strong constraints on the origin of the emission. They can also reveal the variability of the phenomenon. Aims. We aim to demonstrate the persistence of the phenomenon over the timescale of a few years and to search for variability of our previously detected excesses. Methods. Using Very Large Telescope Interferometer (VLTI)/Precision Integrated Optics Near Infrared ExpeRiment (PIONIER) in H band we have carried out multi-epoch observations of the stars for which a near-infrared excess was previously detected using the same observation technique and instrument. The detection rates and distribution of the excesses from our original survey and the follow-up observations are compared statistically. A search for variability of the excesses in our time series is carried out based on the level of the broadband excesses. Results. In 12 of 16 follow-up observations, an excess is re-detected with a significance of > 2σ, and in 7 of 16 follow-up observations significant excess (> 3σ) is re-detected. We statistically demonstrate with very high confidence that the phenomenon persists for the majority of the systems. We also present the first detection of potential variability in two sources. Conclusions. We conclude that the phenomenon responsible for the excesses persists over the timescale of a few years for the majority of the systems. However, we also find that variability intrinsic to a target can cause it to have no significant excess at the time of a specific observation.