The innermost astronomical unit in protoplanetary disks is a key region for stellar and planet formation, as exoplanet searches have shown a large occurrence of close-in planets that are located within the first au around their host star. We aim to reveal the morphology of the disk inner rim using near-infrared interferometric observations with milli-arcsecond resolution provided by infrared interferometry. We provide reconstructed images of 15 objects selected from the Herbig AeBe survey carried out with PIONIER at the VLTI, using SPARCO. We find that 40% of the systems are centrosymmetric at the angular resolution of the observations. For the rest of the objects, we find evidence for asymmetric emission due to moderate-to-strong inclination of a disk-like structure for 30% of the objects and noncentrosymmetric morphology due to a nonaxisymmetric and possibly variable environment (30%). Among the systems with a disk-like structure, 20% show a resolved dust-free cavity. The image reconstruction process is a powerful tool to reveal complex disk inner rim morphologies. At the angular resolution reached by near-infrared interferometric observations, most of the images are compatible with a centrally peaked emission (no cavity). For the most resolved targets, image reconstruction reveals morphologies that cannot be reproduced by generic parametric models. Moreover, the nonaxisymmetric disks show that the spatial resolution probed by optical interferometers makes the observations of the near-infrared emission sensitive to temporal evolution with a time-scale down to a few weeks. The evidence of nonaxisymmetric emission that cannot be explained by simple inclination and radiative transfer effects requires alternative explanations, such as a warping of the inner disks. Interferometric observations can, therefore, be used to follow the evolution of the asymmetry of those disks at a sub-au scale.

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.

We present a new formulation of a family of proximity operators that generalize the projector step for phase retrieval. These proximity operators for noisy intensity measurements can replace the classical "noise free" projection in any projection-based algorithm. They are derived from a maximum likelihood formulation and admit closed form solutions for both the Gaussian and the Poisson cases. In addition, we extend these proximity operators to undersampled intensity measurements. To assess their performance, these operators are exploited in a classical Gerchberg Saxton algorithm. We present numerical experiments showing that the reconstructed complex amplitudes with these proximity operators perform always better than using the classical intensity projector while their computational overhead is moderate.

High spatial resolution is the key for the understanding various astrophysical phenomena. But even with the future E-ELT, single dish instruments are limited to a spatial resolution of about 4 mas in the visible. For the closest objects within our Galaxy most of the stellar photosphere remains smaller than 1 mas. With the success of long baseline interferometry these limitations were soom overcome. Today low and high resolution interferometric instruments on the VLTI and CHARA offer an immense range of astrophysical studies. Combining more telescopes and moving to visible wavelengths broadens the science cases even more. With the idea of developing strong science cases for a future visible interferometer, we organized a science group around the following topics: pre-main sequence and main sequence stars, fundamental parameters, asteroseismology and classical pulsating stars, evolved stars, massive stars, active galactic nuclei (AGNs) and imaging techniques. A meeting was organized on the 15th and 16th of January, 2015 in Nice with the support of the Action Specific in Haute Resolution Angulaire (ASHRA), the Programme National en Physique Stellaire (PNPS), the Lagrange Laboratory and the Observatoire de la Cote d'Azur, in order to present these cases and to discuss them further for future visible interferometers. This White Paper presents the outcome of the exchanges. This book is dedicated to the memory of our colleague Olivier Chesneau who passed away at the age of 41.

Context. It is now generally accepted that the near-infrared excess of Herbig AeBe stars originates in the dust of a circumstellar disk. Aims. The aims of this article are to infer the radial and vertical structure of these disks at scales of order one au, and the properties of the dust grains. Methods. The program objects (51 in total) were observed with the H-band (1.6micron) PIONIER/VLTI interferometer. The largest baselines allowed us to resolve (at least partially) structures of a few tenths of an au at typical distances of a few hundred parsecs. Dedicated UBVRIJHK photometric measurements were also obtained. Spectral and 2D geometrical parameters are extracted via fits of a few simple models: ellipsoids and broadened rings with azimuthal modulation. Model bias is mitigated by parallel fits of physical disk models. Sample statistics were evaluated against similar statistics for the physical disk models to infer properties of the sample objects as a group. Results. We find that dust at the inner rim of the disk has a sublimation temperature Tsub~1800K. A ring morphology is confirmed for approximately half the resolved objects; these rings are wide delta_r>=0.5. A wide ring favors a rim that, on the star-facing side, looks more like a knife edge than a doughnut. The data are also compatible with a the combination of a narrow ring and an inner disk of unspecified nature inside the dust sublimation radius. The disk inner part has a thickness z/r~0.2, flaring to z/r~0.5 in the outer part. We confirm the known luminosity-radius relation; a simple physical model is consistent with both the mean luminosity-radius relation and the ring relative width; however, a significant spread around the mean relation is present. In some of the objects we find a halo component, fully resolved at the shortest interferometer spacing, that is related to the HAeBe class.

We propose a wavefront-based method to estimate the PSF over the whole field of view. This method estimate the aberrations of all the mirrors of the telescope using only field stars. In this proof of concept paper, we described the method and present some qualitative results.

MWC158 is a star with the B[e] phenomenon that shows strong spectrophotometric variability (in lines and in UV and visible continuum) attributed to phases of shell ejection. The evolutionary stage of this star was never clearly determined. Previous interferometric, spectropolarimetric and spectro-interferometric studies suggest a disk morphology for its environment. We investigate the origin of the variability within the inner astronomical unit of the central star using near-infrared interferometric observations with PIONIER at the VLTI over a two-year period. We performed an image reconstruction of the circumstellar environment using the SPARCO method. We discovered that the morphology of the circumstellar environment could vary on timescales of weeks or days. We carried out a parametric fit of the data with a model consisting of a star, a disk and a bright spot that represents a brighter emission in the disk. We detect strong morphological changes in the first astronomical unit around the star, that happen on a timescale of few months. We cannot account for such variability well with a binary model. Our parametric model fits the data well and allows us to extract the location of the asymmetry for different epochs. For the first time, we detect a morphological variability in the environment of MWC158. This variability is reproduced by a model of a disk and a bright spot. The locations of the bright spot suggest that it is located in the disk, but its precise motion is not determined. The origin of the asymmetry in the disk is complex and may be related to asymmetric shell ejections.

Astronomical optical interferometers (OI) sample the Fourier transform of the intensity distribution of a source at the observation wavelength. Because of rapid atmospheric perturbations, the phases of the complex Fourier samples (visibilities) cannot be directly exploited , and instead linear relationships between the phases are used (phase closures and differential phases). Consequently, specific image reconstruction methods have been devised in the last few decades. Modern polychromatic OI instruments are now paving the way to multiwavelength imaging. This paper presents the derivation of a spatio-spectral ("3D") image reconstruction algorithm called PAINTER (Polychromatic opticAl INTErferometric Reconstruction software). The algorithm is able to solve large scale problems. It relies on an iterative process, which alternates estimation of polychromatic images and of complex visibilities. The complex visibilities are not only estimated from squared moduli and closure phases, but also from differential phases, which help to better constrain the polychromatic reconstruction. Simulations on synthetic data illustrate the efficiency of the algorithm.

Astronomical optical interferometers sample the Fourier transform of the intensity distribution of a source at the observation wavelength. Because of rapid perturbations caused by atmospheric turbulence, the phases of the complex Fourier samples (visibilities) cannot be directly exploited. Consequently, specific image reconstruction methods have been devised in the last few decades. Modern polychromatic optical interferometric instruments are now paving the way to multiwavelength imaging. This paper is devoted to the derivation of a spatio-spectral (3D) image reconstruction algorithm, coined PAINTER (Polychromatic opticAl INTErferometric Reconstruction software). The algorithm relies on an iterative process, which alternates estimation of polychromatic images and of complex visibilities. The complex visibilities are not only estimated from squared moduli and closure phases, but also differential phases, which helps to better constrain the polychromatic reconstruction. Simulations on synthetic data illustrate the efficiency of the algorithm and in particular the relevance of injecting a differential phases model in the reconstruction.

The emergence of optical interferometers with three and more telescopes allows image reconstruction of astronomical objects at the milliarcsecond scale. However, some objects contain components with very different spectral energy distributions (SED; i.e. different temperatures), which produces strong chromatic effects on the interferograms that have to be managed with care by image reconstruction algorithms. For example, the gray approximation for the image reconstruction process results in a degraded image if the total (u, v)-coverage given by the spectral supersynthesis is used. The relative flux contribution of the central object and an extended structure changes with wavelength for different temperatures. For young stellar objects, the known characteristics of the central object (i.e., stellar SED), or even the fit of the spectral index and the relative flux ratio, can be used to model the central star while reconstructing the image of the extended structure separately. Methods. We present a new method, called SPARCO (semi-parametric algorithm for the image reconstruction of chromatic objects), which describes the spectral characteristics of both the central object and the extended structure to consider them properly when reconstructing the image of the surrounding environment. We adapted two image-reconstruction codes (Macim, Squeeze, and MiRA) to implement this new prescription. SPARCO is applied using Macim, Squeeze and MiRA on a young stellar object model and also on literature data on HR5999 in the near-infrared with the VLTI. This method paves the way to improved aperture-synthesis imaging of several young stellar objects with existing datasets. More generally, the approach can be used on astrophysical sources with similar features such as active galactic nuclei, planetary nebulae, and asymptotic giant branch stars.

Optical interferometers provide multiple wavelength measurements. In order to fully exploit the spectral and spatial resolution of these instruments, new algorithms for image reconstruction have to be developed. Early attempts to deal with multi-chromatic interferometric data have consisted in recovering a gray image of the object or independent monochromatic images in some spectral bandwidths. The main challenge is now to recover the full 3-D (spatio-spectral) brightness distribution of the astronomical target given all the available data. We describe a new approach to implement multi-wavelength image reconstruction in the case where the observed scene is a collection of point-like sources. We show the gain in image quality (both spatially and spectrally) achieved by globally taking into account all the data instead of dealing with independent spectral slices. This is achieved thanks to a regularization which favors spatial sparsity and spectral grouping of the sources. Since the objective function is not differentiable, we had to develop a specialized optimization algorithm which also accounts for non-negativity of the brightness distribution.

We present the results of the fifth Interferometric Imaging Beauty Contest. The contest consists in blind imaging of test data sets derived from model sources and distributed in the OIFITS format. Two scenarios of imaging with CHARA/MIRC-6T were offered for reconstruction: imaging a T Tauri disc and imaging a spotted red supergiant. There were eight different teams competing this time: Monnier with the software package MACIM; Hofmann, Schertl and Weigelt with IRS; Thi\'ebaut and Soulez with MiRA ; Young with BSMEM; Mary and Vannier with MIROIRS; Millour and Vannier with independent BSMEM and MiRA entries; Rengaswamy with an original method; and Elias with the radio-astronomy package CASA. The contest model images, the data delivered to the contestants and the rules are described as well as the results of the image reconstruction obtained by each method. These results are discussed as well as the strengths and limitations of each algorithm.

In this paper we present a general forward fitting method for multichannel image restoration based on regularized chi2. We introduce separable regularizations that account for the dynamic of the model and take advantage of the continuities present in the data, leaving only two hyper-parameters to tune. We illustrate a practical implementation of this method in the context of host galaxy subtraction for the Nearby SuperNova factory. We show that the image restoration obtained fulfills the stringent requirements on bias and photometricity needed by this program. The reconstruction yields sub-percent integrated residuals in all the synthetic filters considered both on real and simulated data. Even though our implementation is tied to the SNfactory data, the method translates to any hyper-spectral data. As such, it is of direct relevance to several new generation instruments like MUSE. Also, this technique could be applied to multi-band astronomical imaging for which image reconstruction is important, for example to increase image resolution for weak lensing surveys.