Source author record

J. Milli

J. Milli appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.EP astro-ph.SR astro-ph.IM astro-ph.GA

Catalog footprint

What is connected

45works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

preprint2026arXiv

JWST/MIRI coronagraphic search for planets in systems with gapped exoKuiper belts and proper motion anomalies

Over the past decade, ALMA has uncovered a range of substructures within exoKuiper belts, pointing to a population of undetected planets. With JWST's sensitivity, we now have the opportunity to identify these planets thought to be responsible for the observed substructures in debris discs. We present Cycle 1 JWST/MIRI 11.4 μm coronagraphic observations of three exoKuiper belts that exhibit gaps in their radial structures: HD 92945, HD 107146, and HD 206893, to determine whether planets are responsible for carving these structures, as seen in our Solar System with the gas giants. We reduce the JWST/MIRI data using spaceKLIP, and introduce new routines to mitigate the Brighter-Fatter effect and persistence. We do not detect any planet candidates, and all detected objects in the field-of-view are consistent with background stars or galaxies. However, by combining JWST mass limits, archival observational constraints, and astrometric accelerations, we rule out a significant portion of planet parameter space, placing tight constraints on the planets possibly responsible for these gaps. To interpret these results, we explore multiple gap-carving scenarios in discs, either massless or with non-zero mass, including clearing by in-situ planet(s), as well as shaping by inner planets through mean-motion or secular apsidal resonances. Finally, we conclude that the planets causing the proper motion anomaly in these systems must reside within the inner 20 au.

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) I: Motivation, sample, data reduction, and results overview

The outer regions of planetary systems host dusty debris discs analogous to the Kuiper belt (exoKuiper belts), which provide crucial constraints on planet formation and evolution processes. ALMA dust observations have revealed a great diversity, and that some belts contain CO gas, whose origin and implications are uncertain. Most of this progress, however, has been limited by low-resolution observations. We conducted the first ALMA large programme dedicated to debris discs: the ALMA survey to Resolve exoKuiper belt Substructures (ARKS). We selected the 24 most promising belts to constrain their detailed radial and vertical structure, and to characterise the gas content. We constrained the radial and vertical distribution of dust, as well as the presence of asymmetries. For a subset of six belts with CO gas, we constrained the gas distribution and kinematics. To interpret these observations, we used a wide range of dynamical models. The first ARKS results are presented as a series of ten papers. We discovered that up to 33% of our sample exhibits multiple dusty rings. For highly inclined belts, we found that non-Gaussian vertical distributions are common and are indicative of multiple dynamical populations. We also found that 10 of the 24 belts present asymmetries. We find that the CO gas is radially broader than the dust, but this could be an effect of optical depth. At least one system shows non-Keplerian kinematics due to strong pressure gradients, which may have triggered a vortex that trapped dust in an arc. Finally, we find evidence that the micron-sized grains may be affected by gas drag in gas rich systems. ARKS has revealed a great diversity of structures in exoKuiper belts that may arise when they are formed in protoplanetary discs or subsequently via interactions with planets and/or gas. We encourage the community to explore the reduced data and data products.

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) IV: CO gas imaging and overview

CO gas is detected in a significant number of debris discs, but its origin and evolution remains unclear. Key constraints are its mass and spectro-spatial distribution, which are coupled through optical depth and have only been analysed at low to moderate resolution so far. The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) is the first ALMA large program to target debris discs at high spectro-spatial resolution. We used $^{12}$CO and $^{13}$CO J=3-2 line data of 18 ARKS debris belts, 5 of which were already known to host gas, to analyse the spectro-spatial distribution of CO, constrain the gas masses, and to search for gas in the remaining systems. We developed a line-imaging pipeline and produced line cubes for each disc, with a spatial resolution down to $\sim$70 mas and spectral resolution of 26 m s$^{-1}$. Using spectro-spatial shifting and stacking, we produced high signal-to-noise maps, and radial and spectral profiles that reveal the distribution and kinematics of gas in 5 gas-bearing discs. For these discs, we constrained the inner radius of the $^{12}$CO, and found the radial brightness profile of CO peaked interior to the dust ring, but that CO was more radially extended than the dust. We present the first radially resolved $^{12}$CO/$^{13}$CO isotopologue flux ratios in gas-bearing debris discs, which are constant with radius for the majority of systems, indicating $^{12}$CO and $^{13}$CO are both optically thick or thin throughout the discs. We report CO line fluxes/upper limits for all systems and optical depth dependant masses for the 5 gas-bearing systems. Finally, we analysed the $^{12}$CO J=3-2 line luminosities for the ARKS debris discs and discs from the literature. We confirm that gas is mostly detected in young systems. However, the high scatter seen in young/high fractional luminosity systems indicates no trend within the systems with detected gas.

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) V: Comparison between scattered light and thermal emission

Debris discs are analogues to our own Kuiper belt around main-sequence stars and are therefore referred to as exoKuiper belts. They have been resolved at high angular resolution at wavelengths spanning the optical to the submillimetre-millimetre regime. Short wavelengths probe the light scattered by such discs, which is dominated by micron-sized dust particles, while millimetre wavelengths probe the thermal emission of millimetre-sized particles. Determining differences in the dust distribution between millimetre- and micron-sized dust is fundamental to revealing the dynamical processes affecting the dust in debris discs. We aim to compare the scattered light from the discs of the ALMA survey to Resolve exoKuiper belt Substructures (ARKS) with the thermal emission probed by ALMA. We focus on the radial distribution of the dust. We used high-contrast scattered light observations obtained with VLT/SPHERE, GPI, and the HST to uniformly study the dust distribution in those systems and compare it to the dust distribution extracted from the ALMA observations carried out in the course of the ARKS project. We also set constraints on the presence of planets by using these high-contrast images combined with exoplanet evolutionary models. 15 of the 24 discs comprising the ARKS sample are detected in scattered light, with TYC9340-437-1 being imaged for the first time at near-infrared wavelengths. For 6 of those 15 discs, the dust surface density seen in scattered light peaks farther out compared to that observed with ALMA. These 6 discs except one are known to also host cold CO gas. Conversely, the systems without significant offsets are not known to host gas, except one. This observational study suggests that the presence of gas in debris discs may affect the small and large grains differently, pushing the small dust to greater distances where the gas is less abundant.

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) VII: Optically thick gas with broad CO gaussian local line profiles in the HD 121617 disc

CO gas has been detected in $\sim$20 debris discs. We present ALMA observations of the CO-rich HD 121617 debris disc from the ARKS survey. Using high-resolution Band 7 observations of $^{12}CO \ J=3-2$, we analyse local CO line profiles to investigate optical depth, CO mass, and temperature. Spectra are aligned and stacked in concentric annuli to produce local line profiles. The resulting profiles are Gaussian-shaped and broadened by Keplerian shear. The line profiles are modelled using both a simplified toy model and a RADMC-3D model including projection effects and Keplerian shear. Fitting the RADMC-3D model to the $^{13}$CO data, we find that an optically thick model with a temperature of 38 K and a CO mass of $2 \times 10^{-3}$ M$_{\oplus}$ reproduces the observations. The model reproduces the enhanced emission at orbital azimuths of $\sim \pm45^{\circ}$ and $\pm135^{\circ}$, forming an X-shaped structure in the velocity-integrated intensity map, as well as the broader $^{12}$CO linewidth relative to $^{13}$CO. Scaling the model by the ISM abundance ratio ($\sim$77) also reproduces the $^{12}$CO data, though high optical depths and model assumptions limit mass constraints. We find that azimuthally averaged local line profiles appear Gaussian regardless of optical depth, cautioning against their use for distinguishing optically thin and thick emission. We constrain the mean molecular weight to $12.6_{-1.1}^{+1.3}$, dependent on model assumptions. Our $^{13}$CO results suggest that C$^{18}$O may also be optically thick in CO-rich debris discs and that the mean molecular weight is significantly higher than if H$_2$ were the dominant gas species, suggesting a non-primordial composition.

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) VIII: A dust arc and non-Keplerian gas kinematics in HD 121617

ExoKuiper belts around young A-type stars often host CO gas, whose origin is still unclear. The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) includes 6 of these gas-bearing belts, to characterise their dust and gas distributions and investigate the gas origin. As part of ARKS, we observed the gas-rich system HD121617 and discovered an arc of enhanced dust density. In this paper, we analyse in detail the dust and gas distributions and the gas kinematics of this system. We extracted radial and azimuthal profiles of the dust (in the millimetre and near-infrared) and gas emission ($^{12}$CO and $^{13}$CO) from reconstructed images. To constrain the morphology of the arc, we fitted an asymmetric model to the dust emission. To characterise the gas kinematics, we fitted a Keplerian model to the velocity map and extracted the azimuthal velocity profile by deprojecting the data. We find that the dust arc is narrow (1-5 au wide at a radius of 75 au), azimuthally extended, and asymmetric; the emission is more azimuthally compact in the direction of the system's rotation, and represents 13% of the total dust mass (0.2$M_\oplus$). The arc is much less pronounced or absent for small grains and gas. Finally, we find strong non-Keplerian azimuthal velocities at the inner and outer wings of the ring, as was expected due to strong pressure gradients. The dust arc resembles the asymmetries found in protoplanetary discs, often interpreted as the result of dust trapping in vortices. If the gas disc mass is high enough ($\gtrsim20M_\oplus$, requiring a primordial gas origin), both the radial confinement of the ring and the azimuthal arc may result from dust grains responding to gas drag. Alternatively, it could result from planet-disc interactions via mean motion resonances. Further studies should test these hypotheses and may provide a dynamical gas mass estimate in this CO-rich exoKuiper belt.

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS): VI. Asymmetries and offsets

Asymmetries in debris discs provide unique clues to understand the evolution and architecture of planetary systems.** The aim of the ALMA survey to Resolve exoKuiper belt Substructures (ARKS) is to expand our understanding of radial and vertical dust structures, as well as gas distributions and kinematics, in debris discs.** Here, in ARKS~VI, we present a systematic analysis of the asymmetries and stellocentric offsets present in the ALMA continuum data for the ARKS survey. Our aims are to identify asymmetries in debris disc dust distributions, quantify debris disc asymmetry properties, and discuss the potential origins of debris disc asymmetries.** We utilised empirical methods to identify emission asymmetries** and the presence of offset emission between disc centres and the locations of the host stars, via an analysis of their calibration procedures and disc properties. We associated observational asymmetry types** and plausible physical classes** associated with each source. We show that there are ten systems, almost half of the ARKS sample, that host either a continuum emission asymmetry or offset emission. Three systems host offsets (HD15115, HD32297, and HD109573 (HR4796)), four host azimuthal asymmetries (HD9672 (49Ceti), HD92945, HD107146, and HD121617), two host an asymmetry in their major axis (HD10647 (q$^1$ Eri), and HD39060 ($β$ Pic)), and one hosts an asymmetry in their minor axis (HD61005). We attribute the offset asymmetries to non-zero eccentricities, and three of the azimuthal asymmetries to arcs. The presence of an asymmetry or offset in the ARKS sample appears to be correlated with the fractional luminosity of cold dust.** Conclusions: This study demonstrates that debris disc asymmetries in the ARKS sample are common, and plausibly so in the wider population of debris discs at (sub)-millimetre wavelengths.** ** = ABRIDGED FOR ARXIV: FULL ABSTRACT IN PAPER

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS). X. Interpreting the peculiar dust rings around HD 131835

Dusty discs detected around main-sequence stars are thought to be signs of planetesimal belts in which the dust distribution is shaped by collisional and dynamical processes, including interactions with gas if present. The debris disc around the young A-type star HD 131835 is composed of two dust rings at ~65 au and ~100 au, a third unconstrained innermost component, and a gaseous component centred at ~65 au. New ALMA observations show that the inner of the two dust rings is brighter than the outer one, in contrast with previous observations in scattered light. We explore two scenarios that could explain these observations: the two dust rings might represent distinct planetesimal belts with different collisional properties, or only the inner ring might contain planetesimals while the outer ring consists entirely of dust that has migrated outwards due to gas drag. To explore the first scenario, we employed a state-of-the-art collisional evolution code. To test the second scenario, we used a simple dynamical model of dust grain evolution in an optically thin gaseous disc. Collisional models of two planetesimal belts cannot fully reproduce the observations by only varying their dynamical excitation, and matching the data through a different material strength requires an extreme difference in dust composition. The gas-driven scenario can reproduce the location of the outer ring and the brightness ratio of the two rings from scattered light observations, but the resulting outer ring is too faint overall in both scattered light and sub-millimetre emission. The dust rings in HD 131835 could be produced from two planetesimal belts, although how these belts would attain the required extremely different properties needs to be explained. The dust-gas interaction is a plausible alternative explanation and deserves further study using a more comprehensive model.

preprint2024arXiv

Dynamics of the Beta Pictoris planetary system and its falling evaporating bodies

For decades, the spectral variations of Beta Pictoris have been modelled as the result of the evaporation of exocomets close to the star, termed falling evaporating bodies (FEBs). Resonant perturbations by a giant planet have been proposed to explain the dynamical origin of these stargrazers. The disk is now known to harbour two giant planets, Beta Pic b and c, orbiting the star at 9.9 au and 2.7 au. While the former almost matches the planet formerly suspected, the discovery of the latter complicates the picture. We first question the stability of the two-planet system. Then we investigate the dynamics of a disk of planetesimals orbiting the star with both planets to check the validity of the FEB generation mechanism. Symplectic N-body simulations are used to determine which regions of the planetesimal disk are dynamically stable. Then we focus on regions where disk particles are able to reach high eccentricities thanks to resonant mechanisms. The first result is that the system is dynamically stable. Both planets may temporarily fall in 7:1 mean motion resonance (MMR). Then, simulations reveal that the whole region extending between ~1.5 au and ~25 au is unstable to planetary perturbations. However, a disk below 1.5 au survives, which appears to constitute an active source of FEBs via high-order MMRs with Beta Pic c. Beta Pic b acts as a distant perturber that helps sustain the whole process. These simulations rule out the preceding FEB generation mechanism model, which placed their origin at around 4-5 au. Conversely, FEBs are likely to originate from a region much further in and related to MMRs with Beta Pic c. That mechanism also appears to last longer, as new planetesimals are able to continuously enter the MMRs and evolve towards the FEB state. Subsequently, the physical nature of the FEBs may differ from that previously thought, and presumably may not be icy.

preprint2023arXiv

Probing the innermost region of the AU~Microscopii debris disk

AU Mic is a young and nearby M-dwarf star harbouring a circumstellar debris disk and one recently discovered planet on an 8d orbit. Large-scale structures within the disk were also discovered and are moving outward at high velocity. We aim at studying this system with the highest spatial resolution in order to probe the innermost regions and to search for additional low-mass companion or set detection limits. The star was observed with two different techniques probing complementary spatial scales. We obtained new SAM observations with SPHERE, which we combined with data from NACO, PIONIER and GRAVITY. We did not detect additional companions within 0.02-7au from the star. We determined magnitude upper limits for companions of H~9.8mag within 0.02-0.5au, Ks~11.2mag within 0.4-2.4au and L'~10.7mag within 0.7-7au. Using theoretical isochrones, we converted into mass upper limits of ~17Mjup, ~12Mjup and ~9jup, respectively. The PIONIER observations allowed us to determine the angular diameter of AU Mic, 0.825+/-0.050mas, which converts to R = 0.862+/-0.052Rsun. We did not detect the newly discovered planets, but we derived upper limit masses for the innermost region of AU Mic. We do not have any detection with a significance beyond 3sigma, the most significant signal with PIONIER being 2.9sigma and with SPHERE being 1.6σ. We applied the pyMESS2 code to estimate the detection probability of companions by combining radial velocities, SPHERE imaging and our interferometric detection maps. We show that 99% of the companions down to ~0.5Mjup can be detected within 0.02au or 1Mjup down to 0.2au. The low-mass planets orbiting at <0.11au will not be directly detectable with the current AO and interferometric instruments due to its close orbit and very high contrast (~10e-10 in K). It will be also below the angular resolution and contrast limit of the next ELT IR imaging instruments.

preprint2022arXiv

Calibration of quasi-static aberrations in exoplanet direct-imaging instruments with a Zernike phase-mask sensor. IV. Temporal stability of non-common path aberrations in VLT/SPHERE

Coronagraphic imaging of exoplanets using ground-based instruments on large telescopes is intrinsically limited by speckles induced by uncorrected aberrations. These aberrations originate from the imperfect correction of the atmosphere by an extreme adaptive optics system; from static optical defects; or from small opto-mechanical variations due to changes in temperature, pressure, or gravity vector. More than the speckles themselves, the performance of high-contrast imagers is ultimately limited by their temporal stability, since most post-processing techniques rely on difference of images acquired at different points in time. Identifying the origin of the aberrations and the timescales involved is therefore crucial to understanding the fundamental limits of dedicated high-contrast instruments. We previously demonstrated the use of a Zernike wavefront sensor called ZELDA for sensing non-common path aberrations (NCPA) in VLT/SPHERE. We now use ZELDA to investigate the stability of the instrumental aberrations using 5 long sequences of measurements obtained at high cadence on the internal source. Our study reveals two regimes of decorrelation of the NCPA. The first, with a characteristic timescale of a few seconds and an amplitude of a few nanometers, is induced by a fast internal turbulence within the enclosure. The second is a slow quasi-linear decorrelation on the order of a few $10^{-3}$ nm rms/s that acts on timescales from minutes to hours. We use coronagraphic image reconstruction to demonstrate that these two NCPA contributions have a measurable impact on differences of images, and that the fast internal turbulence is a dominating term over to the slow linear decorrelation. We also use dedicated sequences where the derotator and atmospheric dispersion compensators emulate a real observation to demonstrate the importance of performing observations symmetric around the meridian.

preprint2022arXiv

In-depth direct imaging and spectroscopic characterization of the young Solar System analog HD 95086

Context. HD 95086 is a young nearby Solar System analog hosting a giant exoplanet orbiting at 57 au from the star between an inner and outer debris belt. The existence of additional planets has been suggested as the mechanism that maintains the broad cavity between the two belts. Aims. We present a dedicated monitoring of HD 95086 with the VLT/SPHERE instrument to refine the orbital and atmospheric properties of HD 95086 b, and to search for additional planets in this system. Methods. SPHERE observations, spread over ten epochs from 2015 to 2019 and including five new datasets, were used. Combined with archival observations, from VLT/NaCo (2012-2013) and Gemini/GPI (2013-2016), the extended set of astrometric measurements allowed us to refine the orbital properties of HD 95086 b. We also investigated the spectral properties and the presence of a circumplanetary disk around HD 95086 b by using the special fitting tool exploring the diversity of several atmospheric models. In addition, we improved our detection limits in order to search for a putative planet c via the K-Stacker algorithm. Results. We extracted for the first time the JH low-resolution spectrum of HD 95086 b by stacking the six best epochs, and confirm its very red spectral energy distribution. Combined with additional datasets from GPI and NaCo, our analysis indicates that this very red color can be explained by the presence of a circumplanetary disk around planet b, with a range of high-temperature solutions (1400-1600 K) and significant extinction (Av > 10 mag), or by a super-solar metallicity atmosphere with lower temperatures (800-1300 K), and small to medium amount of extinction (Av < 10 mag). We do not find any robust candidates for planet c, but give updated constraints on its potential mass and location.

preprint2021arXiv

A SPHERE survey of self-shadowed planet-forming disks

To date, nearly two hundred planet-forming disks have been imaged with high resolution. Our propensity to study bright and extended objects is however biasing our view of the disk demography. In this work, we contribute to alleviate this bias by analyzing fifteen disks targeted with VLT/SPHERE that look faint in scattered light. Sources were selected based on a low far-IR excess from the spectral energy distribution. The comparison with the ALMA images available for a few sources shows that the scattered light surveyed by these datasets is only detected from a small portion of the disk extent. The mild anti-correlation between the disk brightness and the near-IR excess demonstrates that these disks are self-shadowed: the inner disk rim intercepts much starlight and leaves the outer disk in penumbra. Based on the uniform distribution of the disk brightness in scattered light across all spectral types, self-shadowing would act similarly for inner rims at a different distance from the star. We discuss how the illumination pattern of the outer disk may evolve with time. Some objects in the sample are proposed to be at an intermediate stage toward bright disks from the literature with either no shadow or with sign of azimuthally confined shadows.

preprint2021arXiv

A survey of the linear polarization of directly imaged exoplanets and brown dwarf companions with SPHERE-IRDIS. First polarimetric detections revealing disks around DH Tau B and GSC 6214-210 B

Young giant planets and brown dwarf companions emit near-infrared radiation that can be linearly polarized up to several percent. This polarization can reveal the presence of a circumsubstellar accretion disk, rotation-induced oblateness of the atmosphere, or an inhomogeneous distribution of atmospheric dust clouds. We measured the near-infrared linear polarization of 20 known directly imaged exoplanets and brown dwarf companions with the high-contrast imager SPHERE-IRDIS at the VLT. We reduced the data using the IRDAP pipeline to correct for the instrumental polarization and crosstalk with an absolute polarimetric accuracy <0.1% in the degree of polarization. We report the first detection of polarization originating from substellar companions, with a polarization of several tenths of a percent for DH Tau B and GSC 6214-210 B in H-band. By comparing the measured polarization with that of nearby stars, we find that the polarization is unlikely to be caused by interstellar dust. Because the companions have previously measured hydrogen emission lines and red colors, the polarization most likely originates from circumsubstellar disks. Through radiative transfer modeling, we constrain the position angles of the disks and find that the disks must have high inclinations. The presence of these disks as well as the misalignment of the disk of DH Tau B with the disk around its primary star suggest in situ formation of the companions. For the 18 other companions, we do not detect significant polarization and place subpercent upper limits on their degree of polarization. These non-detections may indicate the absence of circumsubstellar disks, a slow rotation rate of young companions, the upper atmospheres containing primarily submicron-sized dust grains, and/or limited cloud inhomogeneity. Finally, we present images of the circumstellar disks of DH Tau, GQ Lup, PDS 70, Beta Pic, and HD 106906.

preprint2021arXiv

An extended scattered light disk around AT Pyx -- Possible planet formation in a cometary globule

To understand how the multitude of planetary systems that have been discovered come to be, we need to study systems at different evolutionary stages, with different central stars but also in different environments. The most challenging environment for planet formation may be the harsh UV radiation field of nearby massive stars which quickly erodes disks by external photo-evaporation. We have observed the AT Pyx system, located in the head of a cometary globule in the Gum Nebula, to search for signs of ongoing planet formation. We used the extreme adaptive optics imager VLT/SPHERE to observe AT Pyx in polarized light as well as total intensity in the J, H and K-band. Additionally we employed VLT/NACO to observe the system in the L-band. We resolve the disk around AT Pyx in scattered light across multiple wavelengths. We find an extended (>126 au) disk, with an intermediate inclination between 35 deg and 42 deg. The disk shows complex sub-structure and we identify 2 and possibly 3 spiral-like features. Depending on the precise geometry of the disk (which we can not unambiguously infer from our data) the disk may be eccentric with an eccentricity of ~0.16 or partially self-shadowed. The spiral features and possible eccentricity are both consistent with signatures of an embedded gas giant planet equal in mass to Jupiter. Our own observations can rule out brown dwarf companions embedded in the resolved disk, but are not sensitive enough to detect gas giants. AT Pyx is the first disk in a cometray globule in the Gum Nebula which is spatially resolved. By comparison with disks in the Orion Nebula Cluster we note that the extension of the disk may be exceptional for this environment if the external UV radiation field is comparable to other cometary globules in the region. The signposts of ongoing planet formation are intriguing and need to be followed up with higher sensitivity.

preprint2021arXiv

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

The planetary system discovered around the young A-type HR8799 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. We present and exploit new near-infrared images and integral-field spectra of the four gas giants surrounding HR8799 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 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. 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\sim30$), near-infrared (0.94-1.64 $μ$m) spectra of the two innermost planets HR8799d and e. The data were reduced with dedicated algorithms, such as the Karhunen-Loève 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.

preprint2020arXiv

CS Cha B: A disc-obscured M-type star mimicking a polarised planetary companion

Context. Direct imaging provides a steady flow of newly discovered giant planets and brown dwarf companions. These multi-object systems can provide information about the formation of low-mass companions in wide orbits and/or help us to speculate about possible migration scenarios. Accurate classification of companions is crucial for testing formation pathways. Aims. In this work we further characterise the recently discovered candidate for a planetary-mass companion CS Cha b and determine if it is still accreting. Methods. MUSE is a four-laser-adaptive-optics-assisted medium-resolution integral-field spectrograph in the optical part of the spectrum. We observed the CS Cha system to obtain the first spectrum of CS Cha b. The companion is characterised by modelling both the spectrum from 6300 $\unicode{x212B}$ to 9300 $\unicode{x212B}$ and the photometry using archival data from the visible to the near-infrared (NIR). Results. We find evidence of accretion and outflow signatures in H$\mathrmα$ and OI emission. The atmospheric models with the highest likelihood indicate an effective temperature of $3450\pm50$ K with a $\log{g}$ of $3.6\pm0.5$ dex. Based on evolutionary models, we find that the majority of the object is obscured. We determine the mass of the faint companion with several methods to be between 0.07 $M_{\odot}$ and 0.71 $M_{\odot}$ with an accretion rate of $\dot{M} = 4 \times 10^{-11 \pm 0.4}$ Myr$^{-1}$. Conclusions. Our results show that CS Cha B is most likely a mid-M-type star that is obscured by a highly inclined disc, which has led to its previous classification using broadband NIR photometry as a planetary-mass companion. This shows that it is important and necessary to observe over a broad spectral range to constrain the nature of faint companions

preprint2020arXiv

HD 117214 debris disk: scattered-light images and constraints on the presence of planets

We performed observations of the Sco-Cen F star HD 117214 aiming at a search for planetary companions and the characterization of the debris disk structure. HD 117214 was observed with the SPHERE subsystems IRDIS, IFS and ZIMPOL at optical and near-IR wavelengths using angular and polarimetric differential imaging techniques. This provided the first images of scattered light from the debris disk with a spatial resolution reaching 25 mas and an inner working angle $< 0.1''$. With the observations with IRDIS and IFS we derive detection limits for substellar companions. The geometrical parameters of the detected disk are constrained by fitting 3D models for the scattering of an optically thin dust disk. Investigating the possible origin of the disk gap, we introduced putative planets therein and modeled the planet-disk and planet-planet dynamical interactions. The obtained planetary architectures are compared with the detection limit curves. The debris disk has an axisymmetric ring structure with a radius of $0.42(\pm 0.01)''$ or $\sim45$ au and an inclination of $71(\pm 2.5)^\circ$ and exhibits a $0.4''$ ($\sim40$ au) wide inner cavity. From the polarimetric data, we derive a polarized flux contrast for the disk of $(F_{\rm pol})_{\rm disk}/F_{\rm \ast}> (3.1 \pm 1.2)\cdot 10^{-4}$ in the RI band. The fractional scattered polarized flux of the disk is eight times smaller than the fractional infrared flux excess. This ratio is similar to the one obtained for the debris disk HIP 79977 indicating that dust radiation properties are not very different between these two disks. Inside the disk cavity we achieve the high sensitivity limits on planetary companions with a mass down to $\sim 4 M_{\rm J}$ at projected radial separations between $0.2''$ and $0.4''$. We can exclude the stellar companions at a radial separation larger than 75 mas from the star.

preprint2020arXiv

Possible evidence of ongoing planet formation in AB Aurigae. A showcase of the SPHERE/ALMA synergy

Context. Planet formation is expected to take place in the first million years of a planetary system through various processes, which remain to be tested through observations. Aims. With the recent discovery, using ALMA, of two gaseous spiral arms inside the 120 au cavity and connected to dusty spirals, the famous protoplanetary disk around AB Aurigae presents a strong incentive for investigating the mechanisms that lead to giant planet formation. A candidate protoplanet located inside a spiral arm has already been claimed in an earlier study based on the same ALMA data. Methods. We used SPHERE at the Very Large Telescope (VLT) to perform near-infrared (IR) high-contrast imaging of AB Aur in polarized and unpolarized light in order to study the morphology of the disk and search for signs of planet formation. Results. SPHERE has delivered the deepest images ever obtained for AB Aur in scattered light. Among the many structures that are yet to be understood, we identified not only the inner spiral arms, but we also resolved a feature in the form of a twist in the eastern spiral at a separation of about 30 au. The twist of the spiral is perfectly reproduced with a planet-driven density wave model when projection effects are accounted for. We measured an azimuthal displacement with respect to the counterpart of this feature in the ALMA data, which is consistent with Keplerian motion on a 4-yr baseline. Another point sxce is detected near the edge of the inner ring, which is likely the result of scattering as opposed to the direct emission from a planet photosphere. We tentatively derived mass constraints for these two features. Conclusions. The twist and its apparent orbital motion could well be the first direct evidence of a connection between a protoplanet candidate and its manifestation as a spiral imprinted in the gas and dust distributions.

preprint2020arXiv

Searching for the near infrared counterpart of Proxima c using multi-epoch high contrast SPHERE data at VLT

Proxima Centauri is known to host an earth-like planet in its habitable zone; very recently a second candidate planet was proposed based on radial velocities. At quadrature, the expected projected separation of this new candidate is larger than 1 arcsec, making it a potentially interesting target for direct imaging. While difficult, identification of the optical counterpart of this planet would allow detailed characterization of the closest planetary system. We searched for a counterpart in SPHERE images acquired during four years through the SHINE survey. In order to account for the large orbital motion of the planet, we used a method that assumes the circular orbit obtained from radial velocities and exploits the sequence of observations acquired close to quadrature in the orbit. We checked this with a more general approach that considers keplerian motion, K-stacker. We did not obtain a clear detection. The best candidate has S/N=6.1 in the combined image. A statistical test suggests that the probability that this detection is due to random fluctuation of noise is < 1% but this result depends on the assumption that distribution of noise is uniform over the image. The position of this candidate and the orientation of its orbital plane fit well with observations in the ALMA 12m array image. However, the astrometric signal expected from the orbit of the candidate we detected is 3-sigma away from the astrometric motion of Proxima as measured from early Gaia data. This, together with the unexpectedly high flux associated with our direct imaging detection, means we cannot confirm that our candidate is indeed Proxima c. On the other hand, if confirmed, this would be the first observation in imaging of a planet discovered from radial velocities and the second one (after Fomalhaut b) of reflecting circumplanetary material. Further confirmation observations should be done as soon as possible.

preprint2020arXiv

SPHERE+: Imaging young Jupiters down to the snowline

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

preprint2020arXiv

Successful application of PSF-R techniques to the case of the globular cluster NGC6121 (M4)

Precise photometric and astrometric measurements on astronomical images require an accurate knowledge of the Point Spread Function (PSF). When the PSF cannot be modelled directly from the image, PSF-reconstruction techniques become the only viable solution. So far, however, their performance on real observations has rarely been quantified. Aims. In this Letter, we test the performance of a novel hybrid technique, called PRIME, on Adaptive Optics-assisted SPHERE/ZIMPOL observations of the Galactic globular cluster NGC6121. Methods. PRIME couples PSF-reconstruction techniques, based on control-loop data and direct image fitting performed on the only bright point-like source available in the field of view of the ZIMPOL exposures, with the aim of building the PSF model. Results. By exploiting this model, the magnitudes and positions of the stars in the field can be measured with an unprecedented precision, which surpasses that obtained by more standard methods by at least a factor of four for on-axis stars and by up to a factor of two on fainter, off-axis stars. Conclusions. Our results demonstrate the power of PRIME in recovering precise magnitudes and positions when the information directly coming from astronomical images is limited to only a few point-like sources and, thus, paving the way for a proper analysis of future Extremely Large Telescope observations of sparse stellar fields or individual extragalactic objects.

preprint2020arXiv

The challenge of measuring the phase function of debris disks. Application to HR\,4796

Abridged: Debris disks are valuable systems to study dust properties. Because they are optically thin at all wavelengths, we have direct access to the properties of dust grains. One very promising technique to study them is to measure their phase function. Disks that are highly inclined are promising targets as a wider range of scattering angles can be probed. The phase function is usually either inferred by comparing the observations to synthetic disk models assuming a parametrized phase function, or estimating it from the surface brightness of the disk. We argue here that the latter approach can be biased due to projection effects leading to an increase in column density along the major axis of a non flat disk. We present a novel approach to account for those column density effects. The method remains model dependent, as one still requires a disk model to estimate the density variations as a function of the scattering angle. This method allows us however to estimate the shape of the phase function without having to invoke any parametrized form. We apply our method to SPHERE/ZIMPOL observations of HR\,4796 and highlight the differences with previous measurements. Our modelling results suggest that the disk is not vertically flat at optical wavelengths. We discuss some of the caveats of the approach, mostly that our method remains blind to real local increase of the dust density, and that it cannot yet be readily applied to angular differential imaging observations. Similarly to previous studies on HR\,4796, we still cannot reconcile the full picture using a given scattering theory to explain the shape of the phase function, a long lasting problem for debris disks. Nonetheless, we argue that similar effects as the ones highlighted in this study can also bias the determination of the phase function in total intensity.

preprint2019arXiv

RefPlanets: Search for reflected light from extra-solar planets with SPHERE/ZIMPOL

RefPlanets is a guaranteed time observation (GTO) programme that uses the Zurich IMaging POLarimeter (ZIMPOL) of SPHERE/VLT for a blind search for exoplanets in wavelengths from 600-900 nm. The goals of this study are the characterization of the unprecedented high polarimetic contrast and polarimetric precision capabilities of ZIMPOL for bright targets, the search for polarized reflected light around some of the closest bright stars to the Sun and potentially the direct detection of an evolved cold exoplanet for the first time. For our observations of Alpha Cen A and B, Sirius A, Altair, Eps Eri and Tau Ceti we used the polarimetric differential imaging (PDI) mode of ZIMPOL which removes the speckle noise down to the photon noise limit for angular separations >0.6". We describe some of the instrumental effects that dominate the noise for smaller separations and explain how to remove these additional noise effects in post-processing. We then combine PDI with angular differential imaging (ADI) as a final layer of post-processing to further improve the contrast limits of our data at these separations. For good observing conditions we achieve polarimetric contrast limits of 15.0-16.3 mag at the effective inner working angle of about 0.13", 16.3-18.3 mag at 0.5" and 18.8-20.4 mag at 1.5". The contrast limits closer in (<0.6") depend significantly on the observing conditions, while in the photon noise dominated regime (>0.6"), the limits mainly depend on the brightness of the star and the total integration time. We compare our results with contrast limits from other surveys and review the exoplanet detection limits obtained with different detection methods. For all our targets we achieve unprecedented contrast limits. Despite the high polarimetric contrasts we are not able to find any additional companions or extended polarized light sources in the data that has been taken so far.

preprint2019arXiv

The polarimetric imaging mode of VLT/SPHERE/IRDIS I: Description, data reduction and observing strategy

Context. Polarimetric imaging is one of the most effective techniques for high-contrast imaging and characterization of protoplanetary disks, and has the potential to be instrumental in characterizing exoplanets. VLT/SPHERE contains the InfraRed Dual-band Imager and Spectrograph (IRDIS) with a dual-beam polarimetric imaging (DPI) mode, which offers the capability to obtain linear polarization images at high contrast and resolution. Aims. We aim to provide an overview of IRDIS/DPI and study its optical design to improve observing strategies and data reduction. Methods. For H-band observations of TW Hya, we compare two data reduction methods that correct for instrumental polarization effects in different ways: a minimization of the noise image, and a polarimetric-model-based correction method that we present in Paper II of this study. Results. We use observations of TW Hya to illustrate the data reduction. In the images of the protoplanetary disk around this star we detect variability in the polarized intensity and angle of linear polarization with pointing-dependent instrument configuration. We explain these variations as instrumental polarization effects and correct for these effects using our model-based correction method. Conclusions. IRDIS/DPI has proven to be a very successful and productive high-contrast polarimetric imaging system. However, the instrument performance depends on the specific instrument configuration. We suggest adjustments to future observing strategies to optimize polarimetric efficiency in field tracking mode by avoiding unfavourable derotator angles. We recommend reducing on-sky data with the pipeline called IRDAP that includes the model-based correction method (described in Paper II) to optimally account for the remaining telescope and instrumental polarization effects and to retrieve the true polarization state of the incident light.

preprint2019arXiv

The polarimetric imaging mode of VLT/SPHERE/IRDIS II: Characterization and correction of instrumental polarization effects

Context. Circumstellar disks and self-luminous giant exoplanets or companion brown dwarfs can be characterized through direct-imaging polarimetry at near-infrared wavelengths. SPHERE/IRDIS at the Very Large Telescope has the capabilities to perform such measurements, but uncalibrated instrumental polarization effects limit the attainable polarimetric accuracy. Aims. We aim to characterize and correct the instrumental polarization effects of the complete optical system, i.e. the telescope and SPHERE/IRDIS. Methods. We create a detailed Mueller matrix model in the broadband filters Y-, J-, H- and Ks, and calibrate it using measurements with SPHERE's internal light source and observations of two unpolarized stars. We develop a data-reduction method that uses the model to correct for the instrumental polarization effects, and apply it to observations of the circumstellar disk of T Cha. Results. The instrumental polarization is almost exclusively produced by the telescope and SPHERE's first mirror and varies with telescope altitude angle. The crosstalk primarily originates from the image derotator (K-mirror). At some orientations, the derotator causes severe loss of signal (>90% loss in H- and Ks-band) and strongly offsets the angle of linear polarization. With our correction method we reach in all filters a total polarimetric accuracy of <0.1% in the degree of linear polarization and an accuracy of a few degrees in angle of linear polarization. Conclusions. The correction method enables us to accurately measure the polarized intensity and angle of linear polarization of circumstellar disks, and is a vital tool for detecting unresolved (inner) disks and measuring the polarization of substellar companions. We have incorporated the correction method in a highly-automatic end-to-end data-reduction pipeline called IRDAP which is publicly available at https://irdap.readthedocs.io.

preprint2018arXiv

Resolving faint structures in the debris disk around TWA7

Debris disks are the intrinsic by-products of the star and planet formation processes. Most likely due to instrumental limitations and their natural faintness, little is known about debris disks around low-mass stars, especially when it comes to spatially resolved observations. We present new VLT/SPHERE IRDIS Dual-Polarization Imaging (DPI) observations in which we detect the dust ring around the M2 spectral type star TWA\,7. Combined with additional Angular Differential Imaging observations we aim at a fine characterization of the debris disk and setting constraints on the presence of low-mass planets. We model the SPHERE DPI observations and constrain the location of the small dust grains, as well as the spectral energy distribution of the debris disk, using the results inferred from the observations, and perform simple N-body simulations. We find that the dust density distribution peaks at 25 au, with a very shallow outer power-law slope, and that the disk has an inclination of 13 degrees with a position angle of 90 degrees East of North. We also report low signal-to-noise detections of an outer belt at a distance of ~52 au from the star, of a spiral arm in the Southern side of the star, and of a possible dusty clump at 3.9 au. These findings seem to persist over timescales of at least a year. Using the intensity images, we do not detect any planets in the close vicinity of the star, but the sensitivity reaches Jovian planet mass upper limits. We find that the SED is best reproduced with an inner disk at 7 au and another belt at 25 au. We report the detections of several unexpected features in the disk around TWA\,7. A yet undetected 100 M$_\oplus$ planet with a semi-major axis at 20-30 au could possibly explain the outer belt as well as the spiral arm. We conclude that stellar winds are unlikely to be responsible for the spiral arm.

preprint2016arXiv

Compositional characterisation of the Themis family

Context. It has recently been proposed that the surface composition of icy main-belt asteroids (B-,C-,Cb-,Cg-,P-,and D-types) may be consistent with that of chondritic porous interplanetary dust particles (CPIDPs). Aims. In the light of this new association, we re-examine the surface composition of a sample of asteroids belonging to the Themis family in order to place new constraints on the formation and evolution of its parent body. Methods. We acquired NIR spectral data for 15 members of the Themis family and complemented this dataset with existing spectra in the visible and mid-infrared ranges to perform a thorough analysis of the composition of the family. Assuming end-member minerals and particle sizes (<2\mum) similar to those found in CPIDPs, we used a radiative transfer code adapted for light scattering by small particles to model the spectral properties of these asteroids. Results. Our best-matching models indicate that most objects in our sample possess a surface composition that is consistent with the composition of CP IDPs.We find ultra-fine grained Fe-bearing olivine glasses to be among the dominant constituents. We further detect the presence of minor fractions of Mg-rich crystalline silicates. The few unsuccessfully matched asteroids may indicate the presence of interlopers in the family or objects sampling a distinct compositional layer of the parent body. Conclusions. The composition inferred for the Themis family members suggests that the parent body accreted from a mixture of ice and anhydrous silicates (mainly amorphous) and subsequently underwent limited heating. By comparison with existing thermal models that assume a 400km diameter progenitor, the accretion process of the Themis parent body must have occurred relatively late (>4Myr after CAIs) so that only moderate internal heating occurred in its interior, preventing aqueous alteration of the outer shell.

preprint2016arXiv

Direct detection of scattered light gaps in the transitional disk around HD 97048 with VLT/SPHERE

We studied the well known circumstellar disk around the Herbig Ae/Be star HD 97048 with high angular resolution to reveal undetected structures in the disk, which may be indicative of disk evolutionary processes such as planet formation. We used the IRDIS near-IR subsystem of the extreme adaptive optics imager SPHERE at the ESO/VLT to study the scattered light from the circumstellar disk via high resolution polarimetry and angular differential imaging. We imaged the disk in unprecedented detail and revealed four ring-like brightness enhancements and corresponding gaps in the scattered light from the disk surface with radii between 39 au and 341 au. We derived the inclination and position angle as well as the height of the scattering surface of the disk from our observational data. We found that the surface height profile can be described by a single power law up to a separation ~270 au. Using the surface height profile we measured the scattering phase function of the disk and found that it is well consistent with theoretical models of compact dust aggregates. We discuss the origin of the detected features and find that low mass (< 1 M_Jup) nascent planets are a possible explanation.

preprint2016arXiv

Discovery of concentric broken rings at sub-arcsec separations in the HD 141569A gas-rich, debris disk with VLT/SPHERE

Transition disks correspond to a short stage between the young protoplanetary phase and older debris phase. Along this evolutionary sequence, the gas component disappears leaving room for a dust-dominated environment where already-formed planets signpost their gravitational perturbations. We endeavor to study the very inner region of the well-known and complex debris, but still gas-rich disk, around HD 141569A using the exquisite high-contrast capability of SPHERE at the VLT. Recent near-infrared (IR) images suggest a relatively depleted cavity within ~200 au, while former mid-IR data indicate the presence of dust at separations shorter than ~100 au. We obtained multi-wavelength images in the near-IR in J, H2, H3 and Ks bands with the IRDIS camera and a 0.95-1.35 micrometers spectral data cube with the IFS. Data were acquired in pupil-tracking mode, thus allowing for angular differential imaging. We discovered several new structures inside 1", of which the most prominent is a bright ring with sharp edges (semi-major axis: 0.4") featuring a strong north-south brightness asymmetry. Other faint structures are also detected from 0.4" to 1" in the form of concentric ringlets and at least one spiral arm. Finally, the VISIR data at 8.6 micrometers suggests the presence of an additional dust population closer in. Besides, we do not detect companions more massive than 1-3 mass of Jupiter. The performance of SPHERE allows us to resolve the extended dust component, which was previously detected at thermal and visible wavelengths, into very complex patterns with strong asymmetries ; the nature of these asymmetries remains to be understood. Scenarios involving shepherding by planets or dust-gas interactions will have to be tested against these observations.

preprint2015arXiv

A narrow, edge-on disk resolved around HD 106906 with SPHERE

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.

preprint2015arXiv

An Extreme-AO Search for Giant Planets around a White Dwarf --VLT/SPHERE performance on a faint target GD 50

CONTEXT. Little is known about the planetary systems around single white dwarfs although there is strong evidence that they do exist. AIMS. We performed a pilot study with the extreme-AO system on the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) on the Very Large Telescopes (VLT) to look for giant planets around a young white dwarf, GD 50. METHODS. We were awarded science verification time on the new ESO instrument SPHERE. Observations were made with the InfraRed Dual-band Imager and Spectrograph in classical imaging mode in H band. RESULTS. Despite the faintness of the target (14.2 mag in R band), the AO loop was closed and a strehl of 37\% was reached in H band. No objects were detected around GD 50. We achieved a 5-sigma contrast of 6.2, 8.0 and 8.25 mags at 0{\farcs}2, 0{\farcs}4 and 0{\farcs}6 and beyond, respectively. We exclude any substellar objects more massive than 4.0 M$_\textrm{J}$ at 6.2 AU, 2.9 M$_\textrm{J}$ at 12.4 AU and 2.8 M$_\textrm{J}$ at 18.6 AU and beyond. This rivals the previous upper limit set by Spitzer. We further show that SPHERE is the most promising instrument available to search for close-in substellar objects around nearby white dwarfs.

preprint2015arXiv

Asymmetric features in the protoplanetary disk MWC758

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. In this paper, we aim to derive new constraints on the structure of the transition disk MWC 758, to detect non-axisymmetric features and understand their origin. We obtained infrared polarized intensity observations of the protoplanetary disk MWC 758 with SPHERE/VLT at 1.04 microns to resolve scattered light at a smaller inner working angle (0.093") and a higher angular resolution (0.027") than previously achieved. We observe polarized scattered light within 0.53" (148 au) down to the inner working angle (26 au) and detect distinct non-axisymmetric 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. Our observations reveal the complex morphology of the disk MWC758. To understand the origin of the detected features, the combination of high-resolution observations in the submillimeter with ALMA and detailed modeling is needed.

preprint2015arXiv

Direct exoplanet detection and characterization using the ANDROMEDA method: Performance on VLT/NaCo data

Context. The direct detection of exoplanets with high-contrast imaging requires advanced data processing methods to disentangle potential planetary signals from bright quasi-static speckles. Among them, angular differential imaging (ADI) permits potential planetary signals with a known rotation rate to be separated from instrumental speckles that are either statics or slowly variable. The method presented in this paper, called ANDROMEDA for ANgular Differential OptiMal Exoplanet Detection Algorithm is based on a maximum likelihood approach to ADI and is used to estimate the position and the flux of any point source present in the field of view. Aims. In order to optimize and experimentally validate this previously proposed method, we applied ANDROMEDA to real VLT/NaCo data. In addition to its pure detection capability, we investigated the possibility of defining simple and efficient criteria for automatic point source extraction able to support the processing of large surveys. Methods. To assess the performance of the method, we applied ANDROMEDA on VLT/NaCo data of TYC-8979-1683-1 which is surrounded by numerous bright stars and on which we added synthetic planets of known position and flux in the field. In order to accommodate the real data properties, it was necessary to develop additional pre-processing and post-processing steps to the initially proposed algorithm. We then investigated its skill in the challenging case of a well-known target, $β$ Pictoris, whose companion is close to the detection limit and we compared our results to those obtained by another method based on principal component analysis (PCA). Results. Application on VLT/NaCo data demonstrates the ability of ANDROMEDA to automatically detect and characterize point sources present in the image field. We end up with a robust method bringing consistent results with a sensitivity similar to the recently published algorithms, with only two parameters to be fine tuned. Moreover, the companion flux estimates are not biased by the algorithm parameters and do not require a posteriori corrections. Conclusions. ANDROMEDA is an attractive alternative to current standard image processing methods that can be readily applied to on-sky data.

preprint2015arXiv

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

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]

preprint2015arXiv

First light of the VLT planet finder SPHERE. II. The physical properties and the architecture of the young systems PZ Tel and HD 1160 revisited

[Abridged] Context. The young systems PZ Tel and HD 1160, hosting known low-mass companions, were observed during the commissioning of the new planet finder SPHERE with several imaging and spectroscopic modes. Aims. We aim to refine the physical properties and architecture of both systems. Methods. We use SPHERE commissioning data and REM observations, as well as literature and unpublished data from VLT/SINFONI, VLT/NaCo, Gemini/NICI, and Keck/NIRC2. Results. We derive new photometry and confirm the nearly daily photometric variability of PZ Tel A. Using literature data spanning 38 yr, we show that the star also exhibits a long-term variability trend. The 0.63-3.8 mic SED of PZ Tel B allows us to revise its properties: spectral type M7+/-1, Teff=2700+/-100 K, log(g)<4.5 dex, log(L/L_Sun)=-2.51+/-0.10 dex, and mass 38-72 MJ. The 1-3.8 mic SED of HD 1160 B suggests a massive brown dwarf or a low-mass star with spectral type M5.5-7.0, Teff=3000+/-100 K, [M/H]=-0.5-0.0 dex, log(L/L_Sun)=-2.81+/-0.10 dex, and mass 39-168 MJ. We confirm the deceleration and high eccentricity (e>0.66) of PZ Tel B. For e<0.9, the inclination, longitude of the ascending node, and time of periastron passage are well constrained. The system is seen close to an edge-on geometry. We reject other brown dwarf candidates outside 0.25" for both systems, and massive giant planets (>4 MJ) outside 0.5" for the PZ Tel system. We also show that K1-K2 color can be used with YJH low-resolution spectra to identify young L-type companions, provided high photometric accuracy (<0.05 mag) is achieved. Conclusions. SPHERE opens new horizons in the study of young brown dwarfs and giant exoplanets thanks to high-contrast imaging capabilities at optical and near-infrared wavelengths, as well as high signal-to-noise spectroscopy in the near-infrared from low (R~30-50) to medium resolutions (R~350).

preprint2015arXiv

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

The system of four planets around HR8799 offers a unique opportunity 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). In this paper (paper IV), we compare the available spectra and photometry of the planets to known objects and atmospheric models (BT-SETTL14, Cloud-AE60, Exo-REM) to characterize the atmospheric properties of the planets. 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 two 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 Teff = 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. We conclude that an enhanced content in dust and decreased CIA of H2 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.

preprint2015arXiv

Variability and dust filtration in the transition disk J160421.7-213028 observed in optical scattered light

Context. Some of transition disks show asymmetric structures in thermal sub-millimetre emission and optical scattered light. These structures can be the result of planet(s) or companions embedded in the disk. Aims. We aim to detect and analyse the scattered light of the transition disk J160421.7-213028, identify disk structures, and compare the results with previous observations of this disk at other wavelengths. Methods. We obtained and analysed new polarised intensity observations of the transition disk J160421.7-213028 with VLT/SPHERE using the visible light instrument ZIMPOL at $R'$ band (0.626$μ$m). We probe the disk gap down to a radius of confidence of 0.1'' (${\sim}15$ AU at 145 pc). We interpret the results in the context of dust evolution when planets interact with the parental disk. Results. We observe a gap from 0.1 to 0.3'' (${\sim}15$ to 40AU) and a bright annulus as previously detected by HiCIAO $H$ band observations at $1.65μ$m. The radial width of the annulus is around $40$AU, and its centre is at ${\sim}61$AU from the central star. The peak of the reflected light at 0.626$μ$m is located 20 AU inward of the cavity detected in the submillimetre. In addition, we detect a dip at a position angle of ${\sim}46.2 \pm 5.4^\circ$. A dip was also detected with HiCIAO, but located at ${\sim}85^\circ$. If the dip observed with HiCIAO is the same, this suggests an average dip rotation of ${\sim}12^\circ/$year, which is inconsistent with the local Keplerian angular velocity of $\sim$0.8$^\circ$/yr at $\sim$61AU. Conclusions. The spatial discrepancy in the radial emission in J160421.7-213028 at different wavelengths is consistent with dust filtration at the outer edge of a gap carved by a massive planet. The dip rotation can be interpreted by fast variability of the inner disk and/or the presence of a warp or circumplanetary material of a planet at ${\sim}9.6$ AU.

preprint2014arXiv

Very deep images of the innermost regions of the beta Pictoris debris disc at Lp

Very few debris discs have been imaged in scattered light at wavelengths beyond 3 microns because the thermal emission from both the sky and the telescope is generally too strong with respect to the faint emission of a debris disc. We present here the first analysis of a high angular resolution image of the disc of beta Pictoris at 3.8 microns. Our primary objective is to probe the innermost parts of the beta Pictoris debris disc and describe its morphology. We performed extensive forward modelling to correct for the biases induced by angular differential imaging on extended objects and derive the physical parameters of the disc. This work relies on a new analysis of seven archival datasets of beta Pictoris observed with VLT/NaCo in the Lp band, associated with disc forward modelling to correct for the biases induced by that technique. The disc is detected above a 5 sigma level between 0.4" and 3.8". The two extensions have a similar brightness within error bars. We confirm an asymmetry previously observed at larger distances from the star and at shorter wavelengths: the isophotes are more widely spaced on the north-west side (above the disc apparent midplane) than on the south-east side. This is interpreted as a small inclination of the disc combined with anisotropic scattering. Our best-fit model has an inclination of 86 degrees with an anisotropic Henyey- Greenstein coefficient of 0.36. This interpretation is supported by a new asymmetry detected in the disc: the disc is significantly bowed towards the north-west within 3" (above the apparent midplane). We also detect a possible new asymmetry within 1", but at this stage we cannot discern between a real feature and an underlying speckle.

preprint2013arXiv

Can eccentric debris disks be long-lived? A first numerical investigation and application to $ζ^2$ Reticuli

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 $ζ^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 $ζ^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 $ζ^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 $ζ^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 $ζ^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$).

preprint2013arXiv

L'-band AGPM vector vortex coronagraph's first light on VLT/NACO: Discovery of a late-type companion at two beamwidths from an F0V star

Context. High contrast imaging has thoroughly combed through the limited search space accessible with first-generation ground-based adaptive optics instruments and the Hubble Space Telescope. Only a few objects were discovered, and many non-detections reported and statistically interpreted. The field is now in need of a technological breakthrough. Aim. Our aim is to open a new search space with first-generation systems such as NACO at the Very Large Telescope, by providing ground-breaking inner working angle (IWA) capabilities in the L' band. The L' band is a sweet spot for high contrast coronagraphy since the planet-to-star brightness ratio is favorable, while the Strehl ratio is naturally higher. Methods. An annular groove phase mask (AGPM) vector vortex coronagraph optimized for the L' band, made from diamond subwavelength gratings was manufactured and qualified in the lab. The AGPM enables high contrast imaging at very small IWA, potentially being the key to unexplored discovery space. Results. Here we present the installation and successful on-sky tests of an L'-band AGPM coronagraph on NACO. Using angular differential imaging, which is well suited to the rotational symmetry of the AGPM, we demonstrated a ΔL' > 7.5 mag contrast from an IWA ~ 0".09 onwards, during average seeing conditions, and for total integration times of a few hundred seconds.

preprint2012arXiv

An insight in the surroundings of HR4796

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

preprint2012arXiv

beta Pic b position relative to the Debris Disk

Context. We detected in 2009 a giant, close-by planet orbiting β Pic, a young star surrounded with a disk, extensively studied for more than 20 years. We showed that if located on an inclined orbit, the planet could explain several peculiarities of β Pictoris system. However, the available data did not permit to measure the inclination of β Pic b with respect to the disk, and in particular to establish in which component of the disk - 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 on 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 degrades the precision of the disk and planet relative position measurements. Methods. We gathered new NaCo data at Ks band, with a set-up optimized to derive simultaneously the orientation(s) of the disk(s) and that of the planet. Results. We show that the projected position of β Pic b is above the midplane of the main disk. With the current data and knowledge on the system, this implies that β Pic b cannot be located in the main disk. The data rather suggest the planet being located in the inclined component.

preprint2012arXiv

High Contrast Imaging of the Close Environment of HD 142527 -

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

preprint2012arXiv

Impact of angular differential imaging on circumstellar disk images

Direct imaging of circumstellar disks requires high-contrast and high-resolution techniques. The angular differential imaging (ADI) technique is one of them, initially developed for point-like sources but now increasingly applied to extended objects. This new field of application raises many questions because the disk images reduced with ADI depend strongly on the amplitude of field rotation and the ADI data reduction strategy. Both of them directly affect the disk observable properties. Our aim is to characterize the applicability and biases of some ADI data reduction strategies for different disk morphologies. A particular emphasis is placed on parameters mostly used for disks: their surface brightness, their width for a ring, and local features such as gaps or asymmetries. We first present a general method for predicting and quantifying those biases. In a second step we illustrate them for some widely used ADI algorithms applied to typical debris disk morphologies: inclined rings with various inner/outer slopes and width. Last, our aim is also to propose improvements of classical ADI to limit the biases on extended objects. Simulated disks seen under various observing conditions were used to reduce ADI data and quantify the resulting biases. These conclusions complements previous results from NaCo L' real-disk images of HR4796A. ADI induces flux losses on disks. This makes this technique appropriate only for low- to medium-inclination disks. A theoretical criterion is derived to predict the amount of flux loss for a given disk morphology, and quantitative estimates of the biases are given in some specific configurations. These biases alter the disk observable properties, such as the slopes of their surface brightness or their radial/azimuthal extent. Additionally, this work demonstrates that ADI can very easily create artificial features without involving astrophysical processes.

J. Milli

What is connected

Connect this record

See the researcher in context

Building this map preview

45 published item(s)

JWST/MIRI coronagraphic search for planets in systems with gapped exoKuiper belts and proper motion anomalies

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) I: Motivation, sample, data reduction, and results overview

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) IV: CO gas imaging and overview

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) V: Comparison between scattered light and thermal emission

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) VII: Optically thick gas with broad CO gaussian local line profiles in the HD 121617 disc

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) VIII: A dust arc and non-Keplerian gas kinematics in HD 121617

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS): VI. Asymmetries and offsets

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS). X. Interpreting the peculiar dust rings around HD 131835

Dynamics of the Beta Pictoris planetary system and its falling evaporating bodies

Probing the innermost region of the AU~Microscopii debris disk

Calibration of quasi-static aberrations in exoplanet direct-imaging instruments with a Zernike phase-mask sensor. IV. Temporal stability of non-common path aberrations in VLT/SPHERE

In-depth direct imaging and spectroscopic characterization of the young Solar System analog HD 95086

A SPHERE survey of self-shadowed planet-forming disks

A survey of the linear polarization of directly imaged exoplanets and brown dwarf companions with SPHERE-IRDIS. First polarimetric detections revealing disks around DH Tau B and GSC 6214-210 B

An extended scattered light disk around AT Pyx -- Possible planet formation in a cometary globule

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

CS Cha B: A disc-obscured M-type star mimicking a polarised planetary companion

HD 117214 debris disk: scattered-light images and constraints on the presence of planets

Possible evidence of ongoing planet formation in AB Aurigae. A showcase of the SPHERE/ALMA synergy

Searching for the near infrared counterpart of Proxima c using multi-epoch high contrast SPHERE data at VLT

SPHERE+: Imaging young Jupiters down to the snowline

Successful application of PSF-R techniques to the case of the globular cluster NGC6121 (M4)

The challenge of measuring the phase function of debris disks. Application to HR\,4796

RefPlanets: Search for reflected light from extra-solar planets with SPHERE/ZIMPOL

The polarimetric imaging mode of VLT/SPHERE/IRDIS I: Description, data reduction and observing strategy

The polarimetric imaging mode of VLT/SPHERE/IRDIS II: Characterization and correction of instrumental polarization effects

Resolving faint structures in the debris disk around TWA7

Compositional characterisation of the Themis family

Direct detection of scattered light gaps in the transitional disk around HD 97048 with VLT/SPHERE

Discovery of concentric broken rings at sub-arcsec separations in the HD 141569A gas-rich, debris disk with VLT/SPHERE

A narrow, edge-on disk resolved around HD 106906 with SPHERE

An Extreme-AO Search for Giant Planets around a White Dwarf --VLT/SPHERE performance on a faint target GD 50

Asymmetric features in the protoplanetary disk MWC758

Direct exoplanet detection and characterization using the ANDROMEDA method: Performance on VLT/NaCo data

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

First light of the VLT planet finder SPHERE. II. The physical properties and the architecture of the young systems PZ Tel and HD 1160 revisited

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

Variability and dust filtration in the transition disk J160421.7-213028 observed in optical scattered light

Very deep images of the innermost regions of the beta Pictoris debris disc at Lp

Can eccentric debris disks be long-lived? A first numerical investigation and application to $ζ^2$ Reticuli

L'-band AGPM vector vortex coronagraph's first light on VLT/NACO: Discovery of a late-type companion at two beamwidths from an F0V star

An insight in the surroundings of HR4796

beta Pic b position relative to the Debris Disk

High Contrast Imaging of the Close Environment of HD 142527 -

Impact of angular differential imaging on circumstellar disk images