Search

SAXO+ is a planned enhancement of the existing SAXO, the VLT/ SPHERE adaptive optics system, deployed on ESO's Very Large Telescope. This upgrade is designed to significantly enhance the instrument's capacity to detect and analyze young Jupiter-like planets. The pivotal addition in SAXO+ is a second-stage adaptive optics system featuring a dedicated near-infrared pyramid wavefront sensor and a second deformable mirror. This secondary stage is strategically integrated to address any residual wavefront errors persisting after the initial correction performed by the current primary AO loop, SAXO. However, several recent studies clearly showed that in good conditions, even in the current system SAXO, non-common path aberrations (NCPAs) are the limiting factor of the final normalized intensity in focal plane, which is the final metric for ground-based high-contrast instruments. This is likely to be even more so the case with the new AO system, with which the AO residuals will be minimized. Several techniques have already been extensively tested on SPHERE in internal source and/or on-sky and will be presented in this paper. However, the use of a new type of sensor for the second stage, a pyramid wavefront sensor, will likely complicate the correction of these aberrations. Using an end-to-end AO simulation tool, we conducted simulations to gauge the effect of measured SPHERE NCPAs in the coronagraphic image on the second loop system and their correction using focal plane wavefront sensing systems. We finally analyzed how the chosen position of SAXO+ in the beam will impact the evolution of the NCPAs in the new instrument., Comment: 16 pages, 10 figures, submitted to the proceedings of SPIE Astronomical Telescopes + Instrumentation 2024, 13096-357

SPHERE, operating at the VLT since 2014, is currently one of the high-contrast instruments with a higher performance. Its adaptive optics system, known as SAXO, will be upgraded to SAXO+, which features the addition of a second stage of adaptive optics. This stage will use a near-infrared pyramid wavefront sensor to record images of fainter exoplanets around redder stars. In this work, we compare the performance of SAXO and SAXO+. We look for the optimal values of the key system parameters of SAXO+ for various science cases and turbulence conditions. We performed numerical simulations using COMPASS, an end-to-end adaptive optics simulation tool. We simulated perfect coronagraph images of an on-axis point source, and we minimized the residual starlight intensity between 3 and 5 ${\lambda}/D$ as a performance criterion. The explored parameter space includes science cases, turbulence conditions, and key system parameters. In every science case and turbulence condition, SAXO+ reduces the residual starlight intensity inside the correction zone of the second stage by a factor of ten compared to SAXO. The optimal first stage gain is lower for SAXO+ than for SAXO alone. We quantified the gain in performance of SAXO+ when changing the second stage frequency from 2 kHz to 3 kHz, and we conclude that 2 kHz may be sufficient for most realistic conditions. We give the optimal first stage gain as well as the first and second stage frequencies for every seeing, coherence time, and science case. Finally, we find that a 2 ${\lambda_{\mathrm{WFS}}}/D$ pyramid modulation radius is a good trade-off between performance and robustness against varying turbulence conditions. This study shows that the future SAXO+ system will outperform the current SAXO system in all studied cases., Comment: 14 pages, 11 figures

High-contrast imaging (HCI) is a technique designed to observe faint signals near bright sources, such as exoplanets and circumstellar disks. The primary challenge in revealing the faint circumstellar signal near a star is the presence of quasi-static speckles, which can produce patterns on the science images that are as bright, or even brighter, than the signal of interest. Strategies such as angular differential imaging (ADI) or reference-star differential imaging (RDI) aim to provide a means of removing the quasi-static speckles in post-processing. In this paper, we present and discuss the adaptation of state-of-the-art algorithms, initially designed for ADI, to jointly leverage angular and reference-star differential imaging (ARDI) for direct high-contrast imaging of circumstellar disks. Using a collection of high-contrast imaging data sets, we assess the performance of ARDI in comparison to ADI and RDI based on iterative principal component analysis (IPCA). These diverse data sets are acquired under various observing conditions and include the injection of synthetic disk models at various contrast levels. Our results demonstrate that ARDI with IPCA improves the quality of recovered disk images and the sensitivity to planets embedded in disks, compared to ADI or RDI individually. This enhancement is particularly pronounced when dealing with extended sources exhibiting highly ambiguous structures that cannot be accurately retrieved using ADI alone, and when the quality of the reference frames is suboptimal, leading to an underperformance of RDI. We finally apply our method to a sample of real observations of protoplanetary disks taken in star-hopping mode, and propose to revisit the protoplanetary claims associated with these disks., Comment: 21 pages, 18 figures, Accepted to A&A

Debris disks give us the unique opportunity to probe the properties of small $\mu$m-sized particles, allowing us to peer into the constituents of their parent bodies, young analogs of comets and asteroids of our solar system. In the past, studies of the total intensity phase function have proven powerful to constrain the main characteristics of the dust particles in debris disks. Nonetheless, there can remain some degeneracies in the modeling that can be alleviated when considering polarized intensity observations. We obtained new near-IR scattered light observations of four young debris disks which allow us to constrain the degree of linear polarization as a function of the scattering angle. All four debris disks are detected in polarized intensity, and three are also recovered in total intensity. We measured peak degree of polarization of $\lesssim 40$\% for all three disks. We find that the particles must consist of highly refractive and absorbing material. For HD129590, by measuring the polarization fraction beyond the birth ring, we constrain the width of the size distribution to be smaller and smaller, compatible with the effect of radiation pressure. We put these findings to the test and present a self-consistent approach to produce synthetic images, assuming different profiles for the radiation pressure strength, and accounting for the presence of unbound grains. We find the contribution of these grains to be especially critical to reproduce the increasing degree of polarization with stellocentric distances. Some of our results might seem difficult to reconcile with our understanding of cosmic dust but since similar results have been obtained for other disks, we discuss the current limitation of available light scattering models as well as possible avenues to alleviate these unfortunate limitations., Comment: Accepted for publication in A&A, abstract shortened

Extreme adaptive optics instruments have revealed exquisite details on debris discs, allowing to extract the optical properties of the dust particles such as the phase function, the degree of polarisation and the spectral reflectance. These are three powerful diagnostic tools to understand the physical properties of the dust : the size, shape and composition of the dust particles. This can inform us on the population of parent bodies, also called planetesimals, which generate those particles through collisions. It is however very rare to be able to combine all those three observables for the same system, as this requires different high-contrast imaging techniques to suppress the starlight and reveal the faint scattered light emission from the dust. Due to its brightness, the ring detected around the A-type star HR 4796 is a notable exception, with both unpolarised and polarised images covering near-infrared wavelengths. Here, we show how measurements of dust particles in the laboratory can reproduce the observed near-infrared photo-polarimetric properties of the HR 4796 disc. Experimental characterisation of dust allows to bypass the current limitations of dust models to reproduce simultaneously the phase function, the degree of polarisation and the spectral reflectance., Comment: 4 pages, 4 figures, Proceedings of the Annual meeting of the French Society of Astronomy and Astrophysics

Polarisation is a powerful remote-sensing tool to study the nature of particles scattering the starlight. It is widely used to characterise interplanetary dust particles in the Solar System and increasingly employed to investigate extrasolar dust in debris discs' systems. We aim to measure the scattering properties of the dust from the debris ring around HD 181327 at near-infrared wavelengths. We obtained high-contrast polarimetric images of HD 181327 in the H band with the SPHERE / IRDIS instrument on the Very Large Telescope (ESO). We complemented them with archival data from HST / NICMOS in the F110W filter reprocessed in the context of the Archival Legacy Investigations of Circumstellar Environments (ALICE) project. We developed a combined forward-modelling framework to simultaneously retrieve the scattering phase function in polarisation and intensity. We detected the debris disc around HD 181327 in polarised light and total intensity. We measured the scattering phase function and the degree of linear polarisation of the dust at 1.6 micron in the birth ring. The maximum polarisation is 23.6% +/- 2.6% and occurs between a scattering angle of 70 deg and 82 deg. We show that compact spherical particles made of a highly refractive and relatively absorbing material in a differential power-law size distribution of exponent $-3.5$ can simultaneously reproduce the polarimetric and total intensity scattering properties of the dust. This type of material cannot be obtained with a mixture of silicates, amorphous carbon, water ice, and porosity, and requires a more refracting component such as iron-bearing minerals. We reveal a striking analogy between the near-infrared polarisation of comets and that of HD 181327. The methodology developed here combining VLT/SPHERE and HST/NICMOS may be applicable in the future to combine the polarimetric capabilities of SPHERE with the sensitivity of JWST., Comment: 21 pages, 22 figures, accepted in Astronomy & Astrophysics

To keep current global warming below 1.5{\deg}C compared with the pre-industrial era, measures must be taken as quickly as possible in all spheres of society. Astronomy must also make its contribution. In this proceeding, and during the workshop to which it refers, different levers of actions are discussed through various examples: individual efforts, laboratory-level actions, impact evaluation and mitigation in major projects, institutional level, and involvement through collectives., Comment: 5 pages, 0 figures, SF2A conference 2023

The Nancy Grace Roman Space Telescope Coronagraph Instrument will enable the polarimetric imaging of debris disks and inner dust belts in the optical and near-infrared wavelengths, in addition to the high-contrast polarimetric imaging and spectroscopy of exoplanets. The Coronagraph uses two Wollaston prisms to produce four orthogonally polarized images and is expected to measure the polarization fraction with measurement errors < 3% per spatial resolution element. To simulate the polarization observations through the Hybrid Lyot Coronagraph (HLC) and Shaped Pupil Coronagraph (SPC), we model disk scattering, the coronagraphic point-response function, detector noise, speckles, jitter, and instrumental polarization and calculate the Stokes parameters. To illustrate the potential for discovery and a better understanding of known systems with both the HLC and SPC modes, we model the debris disks around Epsilon Eridani and HR 4796A, respectively. For Epsilon Eridani, using astrosilicates with 0.37+/-0.01 as the peak input polarization fraction in one resolution element, we recover the peak disk polarization fraction of 0.33+/-0.01. Similarly, for HR 4796A, for a peak input polarization fraction of 0.92+/-0.01, we obtain the peak output polarization fraction as 0.80+/-0.03. The Coronagraph design meets the required precision, and forward modeling is needed to accurately estimate the polarization fraction., Comment: 21 pages, 13 figures, Accepted in Publications of the Astronomical Society of the Pacific

We show the first SPHERE/IRDIS and IFS data of the CO-rich debris disc around HD 131488. We use N-body simulations to model both the scattered light images and the SED of the disc in a self-consistent way. We apply the Henyey-Greenstein approximation, Mie theory, and the Discrete Dipole Approximation to model the emission of individual dust grains. Our study shows that only when gas drag is taken into account can we find a model that is consistent with scattered light as well as thermal emission data of the disc. The models suggest a gas surface density of $2\times10^{-5}\,M_\oplus/$au$^2$ which is in agreement with estimates from ALMA observations. Thus, our modelling procedure allows us to roughly constrain the expected amount of gas in a debris disc without actual gas measurements. We also show that the shallow size distribution of the dust leads to a significant contribution of large particles to the overall amount of scattered light. The scattering phase function indicates a dust porosity of $\sim0.2\ldots 0.6$ which is in agreement with a pebble pile scenario for planetesimal growth., Comment: 27 pages, 25 figures, 2 tables, accepted for publication in MNRAS

SAXO+ is a proposed upgrade to SAXO, the AO system of the SPHERE instrument on the ESO Very Large Telescope. It will improve the capabilities of the instrument for the detection and characterization of young giant planets. It includes a second stage adaptive optics system composed of a dedicated near-infrared wavefront sensor and a deformable mirror. This second stage will remove the residual wavefront errors left by the current primary AO loop (SAXO). This paper focuses on the numerical simulations of the second stage (SAXO+) and concludes on the impact of the main AO parameters used to build the design strategy. Using an end-to-end AO simulation tool (COMPASS), we investigate the impact of several parameters on the performance of the AO system. We measure the performance in minimizing the star residuals in the coronagraphic image. The parameters that we study are : the second stage frequency, the photon flux on each WFS, the first stage gain and the DM number of actuators of the second stage. We show that the performance is improved by a factor 10 with respect to the current AO system (SAXO). The optimal second stage frequency is between 1 and 2 kHz under good observing conditions. In a red star case, the best SAXO+ performance is achieved with a low first stage gain of 0.05, which reduces the first stage rejection., Comment: 10 pages, 8 figures. Submitted to AO4ELT7 conference proceedings

Diverse protoplanetary disk morphology can result from planet-disk interaction, suggesting planetary presence. To date, most scattered light imaging campaigns have probed polarized light, which is only a fraction of the total light and not very sensitive to planets. To observe and characterize protoplanetary disk systems in the near-infrared in both polarized and total intensity light, we carried out an unprecedented study of scattering properties of disks, as well as of any planetary companions. Using SPHERE with star-hopping at the Very Large Telescope, we observed 29 disk hosts and their reference stars in $K_s$-band polarized light. We extracted disks in total intensity by adopting the data imputation concept with sequential non-negative matrix factorization (DI-sNMF). We obtained high-quality disk images in total intensity for 15 systems and in polarized light for 23. For well-recovered disks in polarized light and total intensity, we parameterized the polarization fraction phase functions using scaled beta distribution: the peak of polarization fraction tentatively correlates with the peak scattering angle, which could be reproduced using certain compact dust, yet more detailed modeling studies are needed. We investigated the empirical DI-sNMF detectability of disks using logistic regression: total intensity detectability of disks primarily depends on host star brightness. For disks with SPHERE data in $Y$-/$J$-/$H$-band, we summarized their polarized color at ~90 deg scattering angle: most of disks are blue in polarized $J-K_s$ color, and they are relatively redder as stellar luminosity increases, indicating larger scatterers. High-quality disk imagery in both total intensity and polarized light thus allows for disk characterization in polarization fraction, and reduces the confusion between disk and planetary signals., Comment: 25 pages, 16 figures, 3 tables, A&A accepted. Data files in FITS format will be publicly available

The discovery of planets orbiting at less than 1 au from their host star and less massive than Saturn in various exoplanetary systems revolutionized our theories of planetary formation. The fundamental question is whether these close-in low-mass planets could have formed in the inner disk interior to 1 au, or whether they formed further out in the planet-forming disk and migrated inward. Exploring the role of additional giant planets in these systems may help us to pinpoint their global formation and evolution. We searched for additional substellar companions by using direct imaging in systems known to host close-in small planets. The use of direct imaging complemented by radial velocity and astrometric detection limits enabled us to explore the giant planet and brown dwarf demographics around these hosts to investigate the potential connection between both populations. We carried out a direct imaging survey with VLT/SPHERE to look for outer giant planets and brown dwarf companions in 27 systems hosting close-in low-mass planets discovered by radial velocity. Our sample is composed of very nearby (<20pc) planetary systems, orbiting G-, K-, and M-type mature (0.5-10Gyr) stellar hosts. We performed homogeneous direct imaging data reduction and analysis to search for and characterize point sources, and derived robust statistical detection limits. Of 337 point-source detections, we do not find any new bound companions. We recovered the emblematic very cool T-type brown dwarf GJ229B. Our typical sensitivities in direct imaging range from 5 to 30 MJup beyond 2 au. The non-detection of massive companions is consistent with predictions based on models of planet formation by core accretion. Our pilot study opens the way to a multi-technique approach for the exploration of very nearby exoplanetary systems with future ground-based and space observatories., Comment: 49 pages including 31 pages of appendices and references, 31 figures, A&A, accepted

The ability to accurately predict the contrast achieved from high contrast imagers is important for efficient scheduling and quality control measures in modern observatories. We aim to consistently predict and measure the raw contrast achieved by SPHERE/IRDIS on a frame by frame basis to improve the efficiency and scientific yield with SPHERE at the Very Large Telescope (VLT).Contrast curves were calculated for over 5 years of archival data using the most common SPHERE/IRDIS coronagraphic mode in the H2/H3 dual band filter, consisting of approximately 80,000 individual frames. These were merged and interpolated with atmospheric data to create a large data-base of contrast curves with associated features. An empirical power law model for contrast, motivated by physical considerations, was then trained and finally tested on an out-of-sample test data set. At an angular separation of 300 mas, the contrast model achieved a mean (out-of-sample) test error of 0.13 magnitudes with the residual 5-95% percentiles between -0.23 and 0.64 magnitude respectively. The models test set root mean square error (RMSE) between 250-600 mas was between 0.31 - 0.40 magnitudes which is equivalent with other state-of-the-art contrast models presented in the literature. In general, the model performed best for targets between 5-9 G-band magnitude, with degraded performance for targets outside this range. This model is currently being incorporated into the Paranal SCUBA software for first level quality control and real time scheduling support.

Context: Circumstellar debris disks provide insight into the formation and early evolution of planetary systems. Resolved belts in particular help to locate planetesimals in exosystems, and can hint at the presence of disk-sculpting exoplanets. Aims: We study the circumstellar environment of HD 112810 (HIP 63439), a mid-F type star in the Sco-Cen association with a significant infrared excess indicating the presence of a circumstellar debris disk. Methods: We collected five high-contrast observations of HD 112810 with VLT/SPHERE. We identified a debris disk in scattered light, and found that the debris signature is robust over a number of epochs and a variety of reduction techniques. We modelled the disk, accounting for self-subtraction and assuming that it is optically thin. Results: We find a single-belt debris disk, with a radius of 118$\pm$9au and an inclination angle of ${75.7}^{+1.1}_{-1.3}$$\deg$. This is in good agreement with the constraints from SED modelling and from a partially-resolved ALMA image of the system. No planets are detected, though planets below the detection limit ($\sim$2.6M$_\textrm{J}$ at a projected separation of 118au) could be present and could have contributed to sculpting the ring of debris. Conclusions: HD 112810 adds to the growing inventory of debris disks imaged in scattered light. The disk is faint, but the radius and the inclination of the disk are promising for follow-up studies of the dust properties., Comment: A&A accepted. 13 pages, 6 figures + appendix

Observations of debris disks have significantly improved over the past decades, both in terms of sensitivity and spatial resolution. At near-infrared wavelengths, new observing strategies and post-processing algorithms allow us to drastically improve the final images, revealing faint structures in the disks. These structures inform us about the properties and spatial distribution of the small dust particles. We present new $H$-band observations of the disk around HD 129590, which display an intriguing arc-like structure in total intensity but not in polarimetry, and propose an explanation for the origin of this arc. Assuming geometric parameters for the birth ring of planetesimals, our model provides the positions of millions of particles of different sizes to compute scattered light images. We demonstrate that if the grain size distribution is truncated or strongly peaks at a size larger than the radiation pressure blow-out size we are able to produce an arc quite similar to the observed one. If the birth ring is radially narrow, given that particles of a given size have similar eccentricities, they will have their apocenters at the same distance from the star. Since this is where the particles will spend most of their time, this results in a "apocenter pile-up" that can look like a ring. Due to more efficient forward scattering this arc only appears in total intensity observations and remains undetected in polarimetric data. This scenario requires sharp variations either in the grain size distribution or for the scattering efficiencies $Q_\mathrm{sca}$. Alternative possibilities such as a wavy size distribution and a size-dependent phase function are interesting candidates to strengthen the apocenter pile-up. We also discuss why such arcs are not commonly detected in other systems, which can mainly be explained by the fact that most parent belts are usually broad., Comment: Accepted for publication in A&A, abstract shortened

Planetesimals -- asteroids and comets -- are the building blocks of planets in protoplanetary discs and the source of dust, ice and gas in debris discs. Along with planets they comprise the left-over material after star formation that constitutes a planetary system. Planets influence the dynamics of planetesimals, sculpting the orbits of debris belts to produce asymmetries or gaps. We can constrain the architecture of planetary systems, and infer the presence of unseen planetary companions, by high spatial resolution imaging of debris discs. HD~16743 is a relatively young F-type star that hosts a bright edge-on debris disc. Based on far-infrared \textit{Herschel} observations its disc was thought to be stirred by a planetary companion. Here we present the first spatially resolved observations at near-infrared and millimetre wavelengths with \textit{HST} and ALMA, revealing the disc to be highly inclined at $87\fdg3~^{+1\fdg9}_{-2\fdg5}$ with a radial extent of 157.7$^{+2.6}_{-1.5}$~au and a FWHM of 79.4$^{+8.1}_{-7.8}$~au ($\Delta R/R = 0.5$). The vertical scale height of the disc is $0.13~\pm~0.02$, significantly greater than typically assumed unstirred value of 0.05, and could be indicative of stirring of the dust-producing planetesimals within the disc by bodies at least a few times the mass of Pluto up to 18.3~$M_{\oplus}$ in the single object limit., Comment: 12 pages, 7 figures, 4 tables, accepted for publication in MNRAS

[abridged] We analyzed mid-infrared high-contrast coronagraphic images of the beta Pictoris system, taking advantage of the NEAR experiment using the VLT/VISIR instrument. The goal of our analysis is to investigate both the detection of the planet beta Pictoris b and of the disk features at mid-IR wavelengths. In addition, by combining several epochs of observation, we expect to constrain the position of the known clumps and improve our knowledge on the dynamics of the disk. To evaluate the planet b flux contribution, we extracted the photometry and compared it to the flux published in the literature. In addition, we used previous data from T-ReCS and VISIR, to study the evolution of the position of the southwest clump that was initially observed in the planetary disk back in 2003. While we did not detect the planet b, we were able to put constraints on the presence of circumplanetary material, ruling out the equivalent of a Saturn-like planetary ring around the planet. The disk presents several noticeable structures, including the known southwest clump. Using a 16-year baseline, sampled with five epochs of observations, we were able to examine the evolution of the clump: the clump orbits in a Keplerian motion with an sma of 56.1+-0.4 au. In addition to the known clump, the images clearly show the presence of a second clump on the northeast side of the disk and fainter and closer structures that are yet to be confirmed. We found correlations between the CO clumps detected with ALMA and the mid-IR images. If the circumplanetary material were located at the Roche radius, the maximum amount of dust determined from the flux upper limit around beta Pictoris b would correspond to the mass of an asteroid of 5 km in diameter. Finally, the Keplerian motion of the southwestern clump is possibly indicative of a yet-to-be-detected planet or signals the presence of a vortex., Comment: Accepted in Astronomy and Astrophysics

Over the past decade, hundreds of nights have been spent on the worlds largest telescopes to search for and directly detect new exoplanets using high-contrast imaging (HCI). Thereby, two scientific goals are of central interest: First, to study the characteristics of the underlying planet population and distinguish between different planet formation and evolution theories. Second, to find and characterize planets in our immediate Solar neighborhood. Both goals heavily rely on the metric used to quantify planet detections and non-detections. Current standards often rely on several explicit or implicit assumptions about the noise. For example, it is often assumed that the residual noise after data post-processing is Gaussian. While being an inseparable part of the metric, these assumptions are rarely verified. This is problematic as any violation of these assumptions can lead to systematic biases. This makes it hard, if not impossible, to compare results across datasets or instruments with different noise characteristics. We revisit the fundamental question of how to quantify detection limits in HCI. We focus our analysis on the error budget resulting from violated assumptions. To this end, we propose a new metric based on bootstrapping that generalizes current standards to non-Gaussian noise. We apply our method to archival HCI data from the NACO-VLT instrument and derive detection limits for different types of noise. Our analysis shows that current standards tend to give detection limit that are about one magnitude too optimistic in the speckle-dominated regime. That is, HCI surveys may have excluded planets that can still exist., Comment: After first iteration with the referee, resubmitted to AJ. Comments welcome!

Since 2005 ESO has been working with its community and industry to develop an extremely large optical/infrared telescope. ESO's Extremely Large Telescope, or ELT for short, is a revolutionary ground-based telescope that will have a 39-metre main mirror and will be the largest visible and infrared light telescope in the world. To address specific topics that are needed for the science operations and calibrations of the telescope, thirteen specific working groups were created to coordinate the effort between ESO, the instrument consortia, and the wider community. We describe here the goals of these working groups as well as their achievements so far., Comment: 8 pages, 4 figures

Multi-wavelength scattered light imaging of debris disks may inform dust properties including typical size and mineral composition. Existing studies have investigated a small set of individual systems across a variety of imaging instruments and filters, calling for uniform comparison studies to systematically investigate dust properties. We obtain the surface brightness of dust particles in debris disks by post-processing coronagraphic imaging observations, and compare the multi-wavelength reflectance of dust. For a sample of resolved debris disks, we perform a systematic analysis on the reflectance properties of their birth rings. We reduced the visible and near-infrared images of 23 debris disk systems hosted by A through M stars using two coronagraphs onboard the Hubble Space Telescope: the STIS instrument observations centering at 0.58 $\mu$m, and the NICMOS instrument at 1.12 $\mu$m or 1.60 $\mu$m. For proper recovery of debris disks, we used classical reference differential imaging for STIS, and adopted non-negative matrix factorization with forward modeling for NICMOS. By dividing disk signals by stellar signals to take into account of intrinsic stellar color effects, we systematically obtained and compared the reflectance of debris birth rings at ~90 deg scattering angle. Debris birth rings typically exhibit a blue color at ~90 deg scattering angle. As the stellar luminosity increases, the color tends to be more neutral. A likely L-shaped color-albedo distribution indicates a clustering of scatterer properties. The observed color trend correlates with the expected blow-out size of dust particles. The color-albedo clustering likely suggests different populations of dust in these systems. More detailed radiative transfer models with realistic dust morphology will contribute to explaining the observed color and color-albedo distribution of debris systems., Comment: 17 pages, 8 figures, 3 tables. A&A accepted

We present a new open-source data-reduction pipeline to reconstruct spectral data cubes from raw SPHERE integral-field spectrograph (IFS) data. The pipeline is written in Python and based on the pipeline that was developed for the CHARIS IFS. It introduces several improvements to SPHERE data analysis that ultimately produce significant improvements in postprocessing sensitivity. We first used new data to measure SPHERE lenslet point spread functions (PSFs) at the four laser calibration wavelengths. These lenslet PSFs enabled us to forward-model SPHERE data, to extract spectra using a least-squares fit, and to remove spectral crosstalk using the measured lenslet PSFs. Our approach also reduces the number of required interpolations, both spectral and spatial, and can preserve the original hexagonal lenslet geometry in the SPHERE IFS. In the case of least-squares extraction, no interpolation of the data is performed. We demonstrate this new pipeline on the directly imaged exoplanet 51 Eri b and on observations of the hot white dwarf companion to HD 2133. The extracted spectrum of HD 2133B matches theoretical models, demonstrating spectrophotometric calibration that is good to a few percent. Postprocessing on two 51 Eri b data sets demonstrates a median improvement in sensitivity of 80% and 30% for the 2015 and 2017 data, respectively, compared to the use of cubes reconstructed by the SPHERE Data Center. The largest improvements are seen for poorer observing conditions. The new SPHERE pipeline takes less than three minutes to produce a data cube on a modern laptop, making it practical to reprocess all SPHERE IFS data., Comment: 17 pages, 11 figures. Software available at: https://github.com/PrincetonUniversity/charis-dep

The imaging and characterization of a larger range of exoplanets, down to young Jupiters and exo-Earths will require accessing very high contrasts at small angular separations with an increased robustness to aberrations, three constraints that drive current instrumentation development. This goal relies on efficient coronagraphs set up on extremely large diameter telescopes such as the Thirty Meter Telescope (TMT), the Giant Magellan Telescope (GMT), or the Extremely Large Telescope (ELT). However, they tend to be subject to specific aberrations that drastically deteriorate the coronagraph performance: their primary mirror segmentation implies phasing errors or even missing segments, and the size of the telescope imposes large spiders, generating low-wind effect as already observed on the Very Large Telescope (VLT)/SPHERE instrument or at the Subaru telescope, or adaptive-optics-due petaling, studied in simulations in the ELT case. The ongoing development of coronagraphs has then to take into account their sensitivity to such errors. We propose an innovative method to generate coronagraphs robust to primary mirror phasing errors and low-wind and adaptive-optics-due petaling effect. This method is based on the apodization of the segment or petal instead of the entire pupil, this apodization being then repeated to mimic the pupil redundancy. We validate this so-called Redundant Apodized Pupil (RAP) method on a James Webb Space Telescope-like pupil composed of 18 hexagonal segments segments to align, and on the VLT architecture in the case of residual low-wind effect., Comment: 13 pages, 10 figures

Telescope pupil fragmentation from spiders generates specific aberrations observed at various telescopes and expected on the large telescopes under construction. This so-called island effect induces differential pistons, tips and tilts on the pupil petals, deforming the instrumental PSF, and is one of the main limitations to the detection of exoplanets with high-contrast imaging. These aberrations have different origins such as the low-wind effect or petaling errors in the adaptive-optics reconstruction. In this paper, we propose to alleviate the impact of the aberrations induced by island effects on high-contrast imaging by adapting the coronagraph design in order to increase its robustness to petal-level aberrations. Following a method first developed for errors due to primary mirror segmentation (segment phasing errors, missing segments...), we develop and test Redundant Apodized Pupils (RAP), i.e. apodizers designed at the petal-scale, then duplicated and rotated to mimic the pupil petal geometry. We apply this concept to the ELT architecture, made of six identical petals, to yield a 10^-6 contrast in a dark region from 8 to 40lambda/D. Both amplitude and phase apodizers proposed in this paper are robust to differential pistons between petals, with minimal degradation to their coronagraphic PSFs and contrast levels. In addition, they are also more robust to petal-level tip-tilt errors than apodizers designed for the whole pupil, with which the limit of contrast of 10^-6 in the coronagraph dark zone is achieved for constraints up to 2 rad RMS of these petal-level modes. The RAP concept proves its robustness to island effects (low-wind effect and post-adaptive optics petaling), with an application to the ELT architecture. It can also be considered for other 8- to 30-meter class ground-based units such as VLT/SPHERE, Subaru/SCExAO, GMT/GMagAO-X, or TMT/PSI., Comment: 11 pages, 15 figures, accepted in A&A

We present 870 $\mu$m ALMA polarization observations of thermal dust emission from the iconic, edge-on debris disk $\beta$ Pic. While the spatially resolved map does not exhibit detectable polarized dust emission, we detect polarization at the $\sim$3$\sigma$ level when averaging the emission across the entire disk. The corresponding polarization fraction is $P_\textrm{frac}$ = $0.51 \pm 0.19$%. The polarization position angle $\chi$ is aligned with the minor axis of the disk, as expected from models of dust grains aligned via radiative alignment torques (RAT) with respect to a toroidal magnetic field ($B$-RAT) or with respect to the anisotropy in the radiation field ($k$-RAT). When averaging the polarized emission across the outer versus inner thirds of the disk, we find that the polarization arises primarily from the SW third. We perform synthetic observations assuming grain alignment via both $k$-RAT and $B$-RAT. Both models produce polarization fractions close to our observed value when the emission is averaged across the entire disk. When we average the models in the inner versus outer thirds of the disk, we find that $k$-RAT is the likely mechanism producing the polarized emission in $\beta$ Pic. A comparison of timescales relevant to grain alignment also yields the same conclusion. For dust grains with realistic aspect ratios (i.e., $s > 1.1$), our models imply low grain-alignment efficiencies., Comment: 22 pages, 13 figures, 2 tables. Accepted for publication in the Astrophysical Journal. Materials accessible in the online version of the ApJ article include the FITS files used to make the ALMA images

The vertical structure of debris disks provides clues about their dynamical evolution and the collision rate of the unseen planetesimals. Thanks to the ever-increasing angular resolution of contemporary instruments and facilities, we are beginning to constrain the scale height of a handful of debris disks, either at near-infrared or millimeter wavelengths. Nonetheless, this is often done for individual targets only. We present here the geometric modeling of eight disks close to edge-on, all observed with the same instrument (SPHERE) and using the same mode (dual-beam polarimetric imaging). Motivated by the presence of CO gas in two out of the eight disks, we then investigate the impact that gas can have on the scale height by performing N-body simulations including gas drag and collisions. We show that gas can quickly alter the dynamics of particles (both in the radial and vertical directions), otherwise governed by gravity and radiation pressure. We find that, in the presence of gas, particles smaller than a few tens of microns can efficiently settle toward the midplane at the same time as they migrate outward beyond the birth ring. For second generation gas ($M_\mathrm{gas} \leq 0.1$ $M_\oplus$), the vertical settling should be best observed in scattered light images compared to observations at millimeter wavelengths. But if the gas has a primordial origin ($M_\mathrm{gas} \geq 1$ $M_\oplus$), the disk will appear very flat both at near-infrared and sub-mm wavelengths. Finally, far beyond the birth ring, our results suggest that the surface brightness profile can be as shallow as $\sim -2.25$., Comment: Accepted in MNRAS, 18 pages, 11 Figures

GQ Lup B is a young and accreting, substellar companion that appears to drive a spiral arm in the circumstellar disk of its host star. We report high-contrast imaging observations of GQ Lup B with VLT/NACO at 4-5 $\mu$m and medium-resolution integral field spectroscopy with VLT/MUSE. The optical spectrum is consistent with an M9 spectral type, shows characteristics of a low-gravity atmosphere, and exhibits strong H$\alpha$ emission. The $H-M'$ color is $\gtrsim$1 mag redder than field dwarfs with similar spectral types and a detailed analysis of the spectral energy distribution (SED) from optical to mid-infrared wavelengths reveals excess emission in the $L'$, NB4.05, and $M'$ bands. The excess flux is well described by a blackbody component with $T_\mathrm{disk} \approx 460$ K and $R_\mathrm{disk} \approx 65\,R_\mathrm{J}$ and is expected to trace continuum emission from small grains in a protolunar disk. We derive an extinction of $A_V \approx 2.3$ mag from the broadband SED with a suspected origin in the vicinity of the companion. We also combine 15 yr of astrometric measurements and constrain the mutual inclination with the circumstellar disk to $84 \pm 9$ deg, indicating a tumultuous dynamical evolution or a stellar-like formation pathway. From the measured H$\alpha$ flux and the estimated companion mass, $M_\mathrm{p} \approx 30\,M_\mathrm{J}$, we derive an accretion rate of $\dot{M} \approx 10^{-6.5}\,M_\mathrm{J}\,\mathrm{yr}^{-1}$. We speculate that the disk is in a transitional stage in which the assembly of satellites from a pebble reservoir has opened a central cavity while GQ Lup B is in the final stages of its formation., Comment: 28 pages, 11 figures, accepted for publication in AJ

We obtained polarimetric differential imaging of a gas-rich debris disk around HD 141569A with SPHERE in the H-band to compare the scattering properties of the innermost ring at 44 au with former observations in total intensity with the same instrument. In polarimetric imaging, we observed that the intensity of the ring peaks in the south-east, mostly in the forward direction, whereas in total intensity imaging, the ring is detected only at the south. This noticeable characteristic suggests a non-uniform dust density in the ring. We implemented a density function varying azimuthally along the ring and generated synthetic images both in polarimetry and in total intensity, which are then compared to the actual data. We find that the dust density peaks in the south-west at an azimuthal angle of $220^{\circ} \sim 238^{\circ}$ with a rather broad width of $61^{\circ} \sim 127^{\circ}$. Although there are still uncertainties that remain in the determination of the anisotropic scattering factor, the implementation of an azimuthal density variation to fit the data proved to be robust. Upon elaborating on the origin of this dust density distribution, we conclude that it could be the result of a massive collision when we account for the effect of the high gas mass that is present in the system on the dynamics of grains. Using the outcome of this modelization, we further measured the polarized scattering phase function for the observed scattering angle between 33$^{\circ}$ and 147$^{\circ}$ as well as the spectral reflectance of the southern part of the ring between 0.98 $\mu$m and 2.1 $\mu$m. We tentatively derived the grain properties by comparing these quantities with MCFOST models and assuming Mie scattering. Our preliminary interpretation indicates a mixture of porous sub-micron sized astro-silicate and carbonaceous grains., Comment: 21 pages, 18 figures, Accepted for publication in Astronomy & Astrophysics journal

We have obtained Hubble Space Telescope (HST) coronagraphic observations of the circumstellar disk around M star TWA 7 using the STIS instrument in visible light. Together with archival observations including HST/NICMOS using the F160W filter and Very Large Telescope/SPHERE at $H$-band in polarized light, we investigate the system in scattered light. By studying this nearly face-on system using geometric disk models and Henyey--Greenstein phase functions, we report new discovery of a tertiary ring and a clump. We identify a layered architecture: three rings, a spiral, and an ${\approx}150$ au$^2$ elliptical clump. The most extended ring peaks at $28$ au, and the other components are on its outskirts. Our point source detection limit calculations demonstrate the necessity of disk modeling in imaging fainter planets. Morphologically, we witness a clockwise spiral motion, and the motion pattern is consistent with both solid body and local Keplerian; we also observe underdensity regions for the secondary ring that might result from mean motion resonance or moving shadows: both call for re-observations to determine their nature. Comparing multi-instrument observations, we obtain blue STIS-NICMOS color, STIS-SPHERE radial distribution peak difference for the tertiary ring, and high SPHERE-NICMOS polarization fraction; these aspects indicate that TWA 7 could retain small dust particles. By viewing the debris disk around M star TWA 7 at a nearly face-on vantage point, our study allows for the understanding of such disks in scattered light in both system architecture and dust property., Comment: 20 pages, 10 figures, 3 tables, ApJ accepted. Data and animation available in the ancillary folder. Citation fixed, author list updated

Fomalhaut C (LP 876-10) is a low mass M4V star in the intriguing Fomalhaut triple system and, like Fomalhaut A, possesses a debris disc. It is one of very few nearby M-dwarfs known to host a debris disc and of these has by far the lowest stellar mass. We present new resolved observations of the debris disc around Fomalhaut C with the Atacama Large Millimetre Array which allow us to model its properties and investigate the system's unique history. The ring has a radius of 26 au and a narrow full width at half maximum of at most 4.2 au. We find a 3$\sigma$ upper limit on the eccentricity of 0.14, neither confirming nor ruling out previous dynamic interactions with Fomalhaut A that could have affected Fomalhaut C's disc. We detect no $^{12}$CO J=3-2 emission in the system and do not detect the disc in scattered light with HST/STIS or VLT/SPHERE. We find the original Herschel detection to be consistent with our ALMA model's radial size. We place the disc in the context of the wider debris disc population and find that its radius is as expected from previous disc radius-host luminosity trends. Higher signal-to-noise observations of the system would be required to further constrain the disc properties and provide further insight to the history of the Fomalhaut triple system as a whole., Comment: 15 pages, 10 figures, 2 tables, Accepted for publication in MNRAS

Optical turbulence affects significantly the quality of ground-based astronomical observations. An accurate and reliable forecast of optical turbulence can help to optimise the scheduling of the science observations and to improve both the quality of the data and the scientific productivity of the observatory. However, forecasts of the turbulence to a level of accuracy that is useful in the operations of large observatories are notoriously difficult to obtain. Several routes have been investigated, from detailed physical modelling of the atmosphere to empirical data-driven approaches. Here, we present an empirical approach exploiting spatial diversity and based on simultaneous measurements between two nearby sites, Cerro Paranal, host of the Very Large Telescope (VLT), and Cerro Armazones, future host of the Extremely Large Telescope (ELT) in Chile. We study the correlation between the high-altitude turbulence as measured between those two sites. This is part of the on-going efforts initiated by the European Southern Observatory (ESO), to obtain short-term forecasts of the turbulence to facilitate the operations of the VLT and prepare the ELT mode of operations., Comment: 13 pages, 9 figures, Proc. of SPIE Paper No. 11448-59

Climate change is affecting and will increasingly affect astronomical observations. In this paper, we investigated the role some key weather parameters play in the quality of astronomical observations, and analysed their long-term trends (longer than 30 years) in order to grasp the impact of climate change on current and future observations. In this preliminary study, we specifically analysed four parameters, the temperature, the surface layer turbulence, the wind speed at the jetstream layer and the humidity. The analyses is conducted with data from the Very Large Telescope (VLT), operated by the European Southern Observatory (ESO), located at Cerro Paranal in the Atacama desert, Chile, which is one of the driest places on Earth. To complete the data from the various sensors installed at Paranal, we used the fifth generation and 20th century European Centre Medium-Range Weather Forecasts (ECMWF) atmospheric reanalysis of the global climate, ERA5 (from 1980 to now) and ERA20C (from 1900 to 2010), which we interpolated at the Paranal observatory location. In addition, we also explored climate projections in this region, using the Coupled Model Intercomparison Project Phase 6 (CMIP6) multi-model ensemble, under the worst-case climate change Shared Socio-Economic Pathways (SSP5-8.5) scenario. Further investigation is needed to better understand the underlying mechanisms of change, as well as to assess the severity of the impact., Comment: Originated at EAS 2020 conference, sustainability session by https://astronomersforplanet.earth - published in Nature Astronomy, September 2020. Free access https://rdcu.be/b61Oh (courtesy of Nature Astronomy)

We have obtained Gemini Planet Imager (GPI) J-, H-, K1-, and K2-Spec observations of the iconic debris ring around the young, main-sequence star HR 4796A. We applied several point-spread function (PSF) subtraction techniques to the observations (Mask-and-Interpolate, RDI-NMF, RDI-KLIP, and ADI-KLIP) to measure the geometric parameters and the scattering phase function for the disk. To understand the systematic errors associated with PSF subtraction, we also forward-modeled the observations using a Markov Chain Monte Carlo framework and a simple model for the disk. We found that measurements of the disk geometric parameters were robust, with all of our analyses yielding consistent results; however, measurements of the scattering phase function were challenging to reconstruct from PSF-subtracted images, despite extensive testing. As a result, we estimated the scattering phase function using disk modeling. We searched for a dependence of the scattering phase function with respect to the GPI filters but found none. We compared the H-band scattering phase function with that measured by Hubble Space Telescope STIS at visual wavelengths and discovered a blue color at small scattering angles and a red color at large scattering angles, consistent with predictions and laboratory measurements of large grains. Finally, we successfully modeled the SPHERE H2 HR 4796A scattered phase function using a distribution of hollow spheres composed of silicates, carbon, and metallic iron., Comment: 34 pages, 30 figures, accepted in ApJ

RZ Psc is a young Sun-like star with a bright and warm infrared excess that is occasionally dimmed significantly by circumstellar dust structures. Optical depth arguments suggest that the dimming events do not probe a typical sight line through the circumstellar dust, and are instead caused by structures that appear above an optically thick mid-plane. This system may therefore be similar to systems where an outer disk is shadowed by material closer to the star. Here we report the discovery that RZ Psc hosts a 0.12$M_\odot$ companion at a projected separation of 23au. We conclude that the disk must orbit the primary star. While we do not detect orbital motion, comparison of the angle of linear polarization of the primary with the companion's on-sky position angle provides circumstantial evidence that the companion and disc may not share the same orbital plane. Whether the companion severely disrupts the disc, truncates it, or has little effect at all, will require further observations of both the companion and disc., Comment: MNRASL in press

Context. Since 1995 and the first discovery of an exoplanet orbiting a main-sequence star, 4000 exoplanets have been discovered using several techniques. However, only a few of these exoplanets were detected through direct imaging. Indeed, the imaging of circumstellar environments requires high-contrast imaging facilities and accurate control of wavefront aberrations. Ground-based planet imagers such as VLT/SPHERE or Gemini/GPI have already demonstrated great performance. However, their limit of detection is hampered by suboptimal correction of aberrations unseen by adaptive optics (AO). Aims. Instead of focusing on the phase minimization of the pupil plane as in standard AO, we aim to directly minimize the stellar residual light in the SPHERE science camera behind the coronagraph to improve the contrast as close as possible to the inner working angle. Methods. We propose a dark hole (DH) strategy optimized for SPHERE. We used a numerical simulation to predict the global improvement of such a strategy on the overall performance of the instrument for different AO capabilities and particularly in the context of a SPHERE upgrade. Then, we tested our algorithm on the internal source with the AO in closed loop. Results. We demonstrate that our DH strategy can correct for aberrations of phase and amplitude. Moreover, this approach has the ability to strongly reduce the diffraction pattern induced by the telescope pupil and the coronagraph, unlike methods operating at the pupil plane. Our strategy enables us to reach a contrast of 5e-7 at 150 mas from the optical axis in a few minutes using the SPHERE internal source. This experiment establishes the grounds for implementing the algorithm on sky in the near future., Comment: Accepted for publication in Astronomy & Astrophysics

Context. The wind driven halo is a feature observed within the images delivered by the latest generation of ground-based instruments equipped with an extreme adaptive optics system and a coronagraphic device, such as SPHERE at the VLT. This signature appears when the atmospheric turbulence conditions are varying faster than the adaptive optics loop can correct. The wind driven halo shows as a radial extension of the point spread function along a distinct direction (sometimes referred to as the butterfly pattern). When present, it significantly limits the contrast capabilities of the instrument and prevents the extraction of signals at close separation or extended signals such as circumstellar disks. This limitation is consequential because it contaminates the data a substantial fraction of the time: about 30% of the data produced by the VLT/SPHERE instrument are affected by the wind driven halo.Aims. This paper reviews the causes of the wind driven halo and presents a method to analyze its contribution directly from the scientific images. Its effect on the raw contrast and on the final contrast after post-processing is demonstrated.Methods. We used simulations and on-sky SPHERE data to verify that the parameters extracted with our method are capable of describing the wind driven halo present in the images. We studied the temporal, spatial and spectral variation of these parameters to point out its deleterious effect on the final contrast.Results. The data driven analysis we propose does provide information to accurately describe the wind driven halo contribution in the images. This analysis justifies why this is a fundamental limitation to the final contrast performance reached.Conclusions. With the established procedure, we will analyze a large sample of data delivered by SPHERE in order to propose, in the future, post-processing techniques tailored to remove the wind driven halo., Comment: Astronomy and Astrophysics - A&A, EDP Sciences, In press

Tiny meteoroids entering the Earth's atmosphere and inducing meteor showers have long been thought to originate partly from cometary dust. Together with other dust particles, they form a huge cloud around the Sun, the zodiacal cloud. From our previous studies of the zodiacal light, as well as other independent methods (dynamical studies, infrared observations, data related to Earth's environment), it is now established that a significant fraction of dust particles entering the Earth's atmosphere comes from Jupiter-family comets (JFCs). This paper relies on our understanding of key properties of the zodiacal cloud and of comet 67P/Churyumov-Gerasimenko, extensively studied by the Rosetta mission to a JFC. The interpretation, through numerical and experimental simulations of zodiacal light local polarimetric phase curves, has recently allowed us to establish that interplanetary dust is rich in absorbing organics and consists of fluffy particles. The ground-truth provided by Rosetta presently establishes that the cometary dust particles are rich in organic compounds and consist of quite fluffy and irregular aggregates. Our aims are as follows: (1) to make links, back in time, between peculiar micrometeorites, tiny meteoroids, interplanetary dust particles, cometary dust particles, and the early evolution of the Solar System, and (2) to show how detailed studies of such meteoroids and of cometary dust particles can improve the interpretation of observations of dust in protoplanetary and debris disks. Future modeling of dust in such disks should favor irregular porous particles instead of more conventional compact spherical particles., Comment: 15 pages, 5 figures, accepted for publication in Planetary and Space Science

Point spread function (PSF) reconstruction (PSF-R) is a well established technique to determine reliably and accurately the PSF from Adaptive Optics (AO) control loop data. We have successfully applied this technique to improve the precision on photometry and astrometry to observation of NGC6121 obtained with SPHERE/ZIMPOL as it will be presented in a forthcoming letter. Firstly, we present the methodology we followed to reconstruct the PSF combining pupil-plane and focal-plane measurements using using our PSF-R method PRIME (Beltramo-Martin et al. 2019), with upgrade of both the model and best-fitting steps compared to previous papers. Secondly, we highlight that PRIME allows to maintain the PSF fitting residual below 0.2% over 2 hours of observation and using only 30 s of AO telemetry, which may have important consequences for telemetry storage for PSF-R purpose on future 30-40 m class telescopes. Finally, we deploy PRIME in a more realistic regime using faint stars so as to identify the precision needed on the initial guess parameters to ensure the convergence towards the optimal solution.

METIS is one of the three first-light instruments planned for the ELT, mainly dedicated to high contrast imaging in the mid-infrared. On the SPHERE high-contrast instrument currently installed at the VLT, we observe that one of the main contrast limitations is the wind driven halo, due to the limited AO running speed with respect to the atmospheric turbulence temporal evolution. From this observation, we extrapolate this signature to the ELT/METIS instrument, which is equipped with a single conjugated adaptive optics system and with several coronagraphic devices. By making use of an analytic AO simulator, we compare the amount of wind driven halo observed with SPHERE and with METIS, under the same turbulence conditions., Comment: AO4ELT6 Conference Proceedings, 12 pages, 8 figures

Context. Adaptive optics (AO) systems greatly increase the resolution of large telescopes, but produce complex point spread function (PSF) shapes, varying in time and across the field of view. This PSF must be accurately known since it provides crucial information about optical systems for design, characterisation, diagnostics and image post processing. Aims. We develop here a model of the AO long exposure PSF, adapted to various seeing conditions and any AO system. This model is made to match accurately both the core of the PSF and its turbulent halo. Methods. The PSF model we develop is based on a parsimonious parameterization of the phase power spectral density with only five parameters to describe circularly symmetric PSFs and seven parameters for asymmetrical ones. Moreover, one of the parameters is directly the Fried parameter r0 of the turbulence s strength. This physical parameter is an asset in the PSF model since it can be correlated with external measurements of the r0, such as phase slopes from the AO real time computer (RTC) or site seeing monitoring. Results. We fit our model against endtoend simulated PSFs using OOMAO tool, and against on sky PSFs from the SPHERE ZIMPOL imager and the MUSE integral field spectrometer working in AO narrowfield mode. Our model matches the shape of the AO PSF both in the core and the halo, with a sub 1 percent relative error for simulated and experimental data. We also show that we retrieve the r0 parameter with subcentimeter precision on simulated data. For ZIMPOL data, we show a correlation of 97 percent between our r0 estimation and the RTC estimation. Finally, MUSE allows us to test the spectral dependency of the fitted r0 parameter. It follows the theoretical $\lambda^{6/5}$ evolution with a standard deviation of 0.3 cm. Evolution of other PSF parameters, such as residual phase variance or aliasing, is also discussed.

In this white paper, we recommend the European Space Agency plays a proactive role in developing a global collaborative effort to construct a large high-contrast imaging space telescope, e.g. as currently under study by NASA. Such a mission will be needed to characterize a sizable sample of temperate Earth-like planets in the habitable zones of nearby Sun-like stars and to search for extraterrestrial biological activity. We provide an overview of relevant European expertise, and advocate ESA to start a technology development program towards detecting life outside the Solar system., Comment: White paper for ESA Voyage 2050; 24 pages

We present new high fidelity optical coronagraphic imagery of the inner $\sim$50 au of AU Mic's edge-on debris disk using the BAR5 occulter of the Hubble Space Telescope Imaging Spectrograph (HST/STIS) obtained on 26-27 July 2018. This new imagery reveals that "feature A", residing at a projected stellocentric separation of 14.2 au on SE-side of the disk, exhibits an apparent "loop-like" morphology at the time of our observations. The loop has a projected width of 1.5 au and rises 2.3 au above the disk midplane. We also explored TESS photometric observations of AU Mic that are consistent with evidence of two starspot complexes in the system. The likely co-alignment of the stellar and disk rotational axes breaks degeneracies in detailed spot modeling, indicating that AU Mic's projected magnetic field axis is offset from its rotational axis. We speculate that small grains in AU Mic's disk could be sculpted by a time-dependent wind that is influenced by this offset magnetic field axis, analogous to co-rotating Solar interaction regions that sculpt and influence the inner and outer regions of our own Heliosphere. Alternatively, if the observed spot modulation is indicative of a significant mis-alignment of the stellar and disk rotational axes, we suggest the disk could still be sculpted by the differential equatorial versus polar wind that it sees with every stellar rotation., Comment: Accepted to ApJL

In a multi-wavelength study of thermal emission and scattered light images we analyse the dust properties and structure of the debris disc around the A1-type main sequence star 49~Cet. As a basis for this study, we present new scattered light images of the debris disc known to possess both a high amount of dust and gas. The outer region of the disc is revealed in former coronagraphic H-band and our new Y-band images from the Very Large Telescope SPHERE instrument. We use the knowledge of the disc's radial extent inferred from ALMA observations and the grain size distribution found by SED fitting to generate semi-dynamical dust models of the disc. We compare the models to scattered light and thermal emission data and find that a disc with a maximum of the surface density at 110~au and shallow edges can describe both thermal emission and scattered light observations. This suggests that grains close to the blow-out limit and large grains stem from the same planetesimal population and are mainly influenced by radiation pressure. The influence of inwards transport processes could not be analysed in this study., Comment: Accepted by MNRAS

Various structures are visible within Spectro-Polarimetric High-contrast Exoplanet REsearch instrument (SPHERE) images that are not always straightforward to interpret. In this article we present a review of these features and demonstrate their origin using simulations. We also identify which expected or unexpected features are limiting the contrast reached by the instrument and how they may be tackled. This vision paves the way to designing a future upgrade of the SPHERE instrument and the next generation of high-contrast instruments such as those planned for the Extremely Large Telescope (ELT).

Debris disks are extrasolar analogs to our own Kuiper Belt and they are detected around at least 17% of nearby Sun-like stars. The morphology and dynamics of a disk encode information about its history, as well as that of any exoplanets within the system. We used ALMA to obtain 1.3 mm observations of the debris disk around the nearby F5V star HD 170773. We image the face-on ring and determine its fundamental parameters by forward-modeling the interferometric visibilities through a Markov Chain Monte Carlo approach. Using a symmetric Gaussian surface density profile, we find a 71 $\pm$ 4 au wide belt with a radius of 193$^{+2}_{-3}$ au, a relatively large radius compared to most other millimeter-resolved belts around late A / early F type stars. This makes HD 170773 part of a group of four disks around A and F stars with radii larger than expected from the recently reported planetesimal belt radius - stellar luminosity relation. Two of these systems are known to host directly imaged giant planets, which may point to a connection between large belts and the presence of long-period giant planets. We also set upper limits on the presence of CO and CN gas in the system, which imply that the exocomets that constitute this belt have CO and HCN ice mass fractions of <77% and <3%, respectively, consistent with Solar System comets and other exocometary belts., Comment: 14 pages, 6 figures, accepted to ApJ

The inclination distribution of the Kuiper belt provides unique constraints on its origin and dynamical evolution, motivating vertically resolved observations of extrasolar planetesimal belts. We present ALMA observations of millimeter emission in the near edge-on planetesimal belt around $\beta$ Pictoris, finding that the vertical distribution is significantly better described by the sum of two Gaussians compared to a single Gaussian. This indicates that, as for the Kuiper belt, the inclination distribution of $\beta$ Pic's belt is better described by the sum of dynamically hot and cold populations rather than a single component. The hot and cold populations have RMS inclinations of 8.9$^{+0.7}_{-0.5}$ and 1.1$^{+0.5}_{-0.5}$ degrees. We also report that an axisymmetric belt model provides a good fit to new and archival ALMA visibilities, and confirm that the midplane is misaligned with respect to $\beta$ Pic b's orbital plane. However, we find no significant evidence for either the inner disk tilt observed in scattered light and CO emission or the South-West/North-East (SW/NE) asymmetry previously reported for millimeter emission. Finally, we consider the origin of the belt's inclination distribution. Secular perturbations from $\beta$ Pic b are unlikely to provide sufficient dynamical heating to explain the hot population throughout the belt's radial extent, and viscous stirring from large bodies within the belt alone cannot reproduce the two populations observed. This argues for an alternative or additional scenario, such as planetesimals being born with high inclinations, or the presence of a `$\beta$ Pic c' planet, potentially migrating outwards near the belt's inner edge., Comment: 19 pages, 12 figures; Accepted for publication in AJ

Nearly all young stars are initially surrounded by `protoplanetary' discs of gas and dust, and in the case of single stars at least 30\% of these discs go on to form planets. The process of protoplanetary disc formation can result in initial misalignments, where the disc orbital plane is different to the stellar equator in single star systems, or to the binary orbital plane in systems with two stars. A quirk of the dynamics means that initially misaligned `circumbinary' discs -- those that surround two stars -- are predicted to evolve to one of two possible stable configurations, one where the disc and binary orbital planes are coplanar, and one where they are perpendicular (a `polar' configuration). Prior work has found coplanar circumbinary discs, but no polar examples were known until now. Here we report the first discovery of a protoplanetary circumbinary disc in the polar configuration, supporting the predictions that such discs should exist. The disc shows some characteristics that are similar to discs around single stars, and that are attributed to dust growth. Thus, the first stages of planet formation appear able to proceed in polar circumbinary discs., Comment: Submitted version - do not read. Published in Nature Astronomy, free PDF at https://rdcu.be/bgNSO

The Coronagraph Instrument (CGI) for NASA's Wide Field Infrared Survey Telescope (WFIRST) will constitute a dramatic step forward for high-contrast imaging, integral field spectroscopy, and polarimetry of exoplanets and circumstellar disks, aiming to improve upon the sensitivity of current ground-based direct imaging facilities by 2-3 orders of magnitude. Furthermore, CGI will serve as a pathfinder for future exo-Earth imaging and characterization missions by demonstrating wavefront control, coronagraphy, and spectral retrieval in a new contrast regime, and by validating instrument and telescope models at unprecedented levels of precision. To achieve this jump in performance, it is critical to draw on the experience of ground-based high-contrast facilities. We discuss several areas of relevant commonalities, including: wavefront control, post-processing of integral field unit data, and calibration and observing strategies., Comment: Published in Proceedings of the SPIE; Space Telescopes and Instrumentation 2018: Optical, Infrared, and Millimeter Wave. 13 pages, 6 figures

The low wind effect is a phenomenon disturbing the phase of the wavefront in the pupil of a large telescope obstructed by spiders, in the absence of wind. It can be explained by the radiative cooling of the spiders, creating air temperature inhomogeneities across the pupil. Because it is unseen by traditional adaptive optics (AO) systems, thus uncorrected, it significantly degrades the quality of AO-corrected images. We provide a statistical analysis of the strength of this effect as seen by VLT/SPHERE after 4 years of operations. We analyse its dependence upon the wind and temperature conditions. We describe the mitigation strategy implemented in 2017: a specific coating with low thermal emissivity in the mid-infrared was applied on the spiders of Unit Telescope 3. We quantify the improvement in terms of image quality, contrast and wave front error using both focal plane images and measured phase maps., Comment: 20 pages, 21 figures, SPIE 2018 conference proceeding

We present multi-epoch VLT/SPHERE observations of the protoplanetary disk around HD 135344B (SAO 206462). The $J$-band scattered light imagery reveal, with high spatial resolution ($\sim$41 mas, 6.4 au), the disk surface beyond $\sim$20 au. Temporal variations are identified in the azimuthal brightness distributions of all epochs, presumably related to the asymmetrically shading dust distribution in the inner disk. These shadows manifest themselves as narrow lanes, cast by localized density enhancements, and broader features which possibly trace the larger scale dynamics of the inner disk. We acquired visible and near-infrared photometry which shows variations up to 10% in the $JHK$ bands, possibly correlated with the presence of the shadows. Analysis of archival VLTI/PIONIER $H$-band visibilities constrain the orientation of the inner disk to $i = 18.2\deg^{+3.4}_{-4.1}$ and ${\rm PA} = 57.3\deg \pm 5.7\deg$, consistent with an alignment with the outer disk or a minor disk warp of several degrees. The latter scenario could explain the broad, quasi-stationary shadowing in N-NW direction in case the inclination of the outer disk is slightly larger. The correlation between the shadowing and the near-infrared excess is quantified with a grid of radiative transfer models. The variability of the scattered light contrast requires extended variations in the inner disk atmosphere ($H/r \lesssim 0.2$). Possible mechanisms that may cause asymmetric variations in the optical depth ($\Delta\tau\lesssim1$) through the atmosphere of the inner disk include turbulent fluctuations, planetesimal collisions, or a dusty disk wind, possibly enhanced by a minor disk warp. A fine temporal sampling is required to follow day-to-day changes of the shadow patterns which may be a face-on variant of the UX Orionis phenomenon., Comment: 19 pages, 8 figures, accepted for publication in ApJ

The development of adaptive optics (AO) played a major role in modern astronomy over the last three decades. By compensating for the atmospheric turbulence, these systems enable to reach the diffraction limit on large telescopes. In this review, we will focus on high contrast applications of adaptive optics, namely, imaging the close vicinity of bright stellar objects and revealing regions otherwise hidden within the turbulent halo of the atmosphere to look for objects with a contrast ratio lower than 10^-4 with respect to the central star. Such high-contrast AO-corrected observations have led to fundamental results in our current understanding of planetary formation and evolution as well as stellar evolution. AO systems equipped three generations of instruments, from the first pioneering experiments in the nineties, to the first wave of instruments on 8m-class telescopes in the years 2000, and finally to the extreme AO systems that have recently started operations. Along with high-contrast techniques, AO enables to reveal the circumstellar environment: massive protoplanetary disks featuring spiral arms, gaps or other asymmetries hinting at on-going planet formation, young giant planets shining in thermal emission, or tenuous debris disks and micron-sized dust leftover from collisions in massive asteroid-belt analogs. After introducing the science case and technical requirements, we will review the architecture of standard and extreme AO systems, before presenting a few selected science highlights obtained with recent AO instruments., Comment: 24 pages, 14 figures