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.

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.

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.

This paper reports on H-band interferometric observations of Betelgeuse made at the three-telescope interferometer IOTA. We image Betelgeuse and its asymmetries to understand the spatial variation of the photosphere, including its diameter, limb darkening, effective temperature, surrounding brightness, and bright (or dark) star spots. We used different theoretical simulations of the photosphere and dusty environment to model the visibility data. We made images with parametric modeling and two image reconstruction algorithms: MIRA and WISARD. We measure an average limb-darkened diameter of 44.28 +/- 0.15 mas with linear and quadratic models and a Rosseland diameter of 45.03 +/- 0.12 mas with a MARCS model. These measurements lead us to derive an updated effective temperature of 3600 +/- 66 K. We detect a fully-resolved environment to which the silicate dust shell is likely to contribute. By using two imaging reconstruction algorithms, we unveiled two bright spots on the surface of Betelgeuse. One spot has a diameter of about 11 mas and accounts for about 8.5% of the total flux. The second one is unresolved (diameter < 9 mas) with 4.5% of the total flux. Resolved images of Betelgeuse in the H band are asymmetric at the level of a few percent. The MOLsphere is not detected in this wavelength range. The amount of measured limb-darkening is in good agreement with model predictions. The two spots imaged at the surface of the star are potential signatures of convective cells.

Within the context of upcoming full-sky lensing surveys, the edge-preserving non- linear algorithm Aski is presented. Using the framework of Maximum A Posteriori inversion, it aims at recovering the full-sky convergence map from surveys with masks. It proceeds in two steps: CCD images of crowded galactic fields are deblurred using automated edge-preserving deconvolution; once the reduced shear is estimated, the convergence map is also inverted via an edge- preserving method. For the deblurring, it is found that when the observed field is crowded, this gain can be quite significant for realistic ground-based surveys when both positivity and edge-preserving penalties are imposed during the iterative deconvolution. For the convergence inversion, the quality of the reconstruction is investigated on noisy maps derived from the horizon N-body simulation, with and without Galactic cuts, and quantified using one-point statistics, power spectra, cluster counts, peak patches and the skeleton. It is found that the reconstruction is able to interpolate and extrapolate within the Galactic cuts/non-uniform noise; its sharpness-preserving penalization avoids strong biasing near the clusters of the map; it reconstructs well the shape of the PDF as traced by its skewness and kurtosis; the geometry and topology of the reconstructed map is close to the initial map as traced by the peak patch distribution and the skeleton's differential length; the two-points statistics of the recovered map is consistent with the corresponding smoothed version of the initial map; the distribution of point sources is also consistent with the corresponding smoothing, with a significant improvement when edge preserving prior is applied. The contamination of B-modes when realistic Galactic cuts are present is also investigated. Leakage mainly occurs on large scales.

Diffuse envelopes around Mira variables are among the most important sources influencing the chemical evolution of galaxies. However they represent an observational challenge because of their complex spectral features and their rapid temporal variability. We constrained the exact brightness distribution of the Mira star TLep with a model-independent analysis. We obtained single-epoch interferometric observations with a dataset continuous in the spectral domain (1.5-2.4mum) and in the spatial domain (baselines ranging from 11 to 96m). We performed a model independent image reconstruction for each spectral bin using the MIRA software. We completed the analysis by modeling the data with a simple star+layer model inspired from the images. Reconstructed images confirm the general picture of a central star partially obscured by the surrounding molecular shell of changing opacity. At 1.7mum, the shell becomes optically thin, with corresponding emission appearing as a ring circling the star. This is the first direct evidence of the spherical morphology of the molecular shell. Model fitting confirmed a spherical layer of constant size and changing opacity over the wavelengths. Rough modeling points to a continuum opacity within the shell, in addition to the CO and H2O features. Accordingly, it appeared impossible to model the data by a photosphere alone in any of the spectral bins.

Classically, optical and near-infrared interferometry have relied on closure phase techniques to produce images. Such techniques allow us to achieve modest dynamic ranges. In order to test the feasibility of next generation optical interferometers in the context of the VLTI-spectro-imager (VSI), we have embarked on a study of image reconstruction and analysis. Our main aim was to test the influence of the number of telescopes, observing nights and distribution of the visibility points on the quality of the reconstructed images. Our results show that observations using six Auxiliary Telescopes (ATs) during one complete night yield the best results in general and is critical in most science cases; the number of telescopes is the determining factor in the image reconstruction outcome. In terms of imaging capabilities, an optical, six telescope VLTI-type configuration and ~200 meter baseline will achieve 4 mas spatial resolution, which is comparable to ALMA and almost 50 times better than JWST will achieve at 2.2 microns. Our results show that such an instrument will be capable of imaging, with unprecedented detail, a plethora of sources, ranging from complex stellar surfaces to microlensing events.

One of the aims of next generation optical interferometric instrumentation is to be able to make use of information contained in the visibility phase to construct high dynamic range images. Radio and optical interferometry are at the two extremes of phase corruption by the atmosphere. While in radio it is possible to obtain calibrated phases for the science objects, in the optical this is currently not possible. Instead, optical interferometry has relied on closure phase techniques to produce images. Such techniques allow only to achieve modest dynamic ranges. However, with high contrast objects, for faint targets or when structure detail is needed, phase referencing techniques as used in radio interferometry, should theoretically achieve higher dynamic ranges for the same number of telescopes. Our approach is not to provide evidence either for or against the hypothesis that phase referenced imaging gives better dynamic range than closure phase imaging. Instead we wish to explore the potential of this technique for future optical interferometry and also because image reconstruction in the optical using phase referencing techniques has only been performed with limited success. We have generated simulated, noisy, complex visibility data, analogous to the signal produced in radio interferometers, using the VLTI as a template. We proceeded with image reconstruction using the radio image reconstruction algorithms contained in AIPS IMAGR (CLEAN algorithm). Our results show that image reconstruction is successful in most of our science cases, yielding images with a 4 milliarcsecond resolution in K band. (abridged)

The study of protoplanetary disks, where the planets are believed to form, will certainly allow the formation of our Solar System to be understood. To conduct observations of these objects at the milli-arcsecond scale, infrared interferometry provides the right performances for T Tauri, FU Ori or Herbig Ae/Be stars. However, the only information obtained so far are scarce visibility measurements which are directly tested with models. With the outcome of recent interferometers, one can foresee obtaining images reconstructed independently of the models. In fact, several interferometers including IOTA and AMBER on the VLTI already provide the possibility to recombine three telescopes at once and thus to obtain the data necessary to reconstruct images. In this paper, we describe the use of MIRA, an image reconstruction algorithm developed for optical inter- ferometry data (squared visibilities and closure phases) by E. Thiebaut. We foresee also to use the spectral information given by AMBER data to constrain even better the reconstructed images. We describe the use of MIRA to reconstruct images of young stellar objects out of actual data, in particular the multiple system GW Orionis (IOTA, 2004), and discuss the encountered difficulties.

The VLTI Spectro Imager (VSI) was proposed as a second-generation instrument of the Very Large Telescope Interferometer providing the ESO community with spectrally-resolved, near-infrared images at angular resolutions down to 1.1 milliarcsecond and spectral resolutions up to R=12000. Targets as faint as K=13 will be imaged without requiring a brighter nearby reference object. The unique combination of high-dynamic-range imaging at high angular resolution and high spectral resolution enables a scientific program which serves a broad user community and at the same time provides the opportunity for breakthroughs in many areas of astrophysic including: probing the initial conditions for planet formation in the AU-scale environments of young stars; imaging convective cells and other phenomena on the surfaces of stars; mapping the chemical and physical environments of evolved stars, stellar remnants, and stellar winds; and disentangling the central regions of active galactic nuclei and supermassive black holes. VSI will provide these new capabilities using technologies which have been extensively tested in the past and VSI requires little in terms of new infrastructure on the VLTI. At the same time, VSI will be able to make maximum use of new infrastructure as it becomes available; for example, by combining 4, 6 and eventually 8 telescopes, enabling rapid imaging through the measurement of up to 28 visibilities in every wavelength channel within a few minutes. The current studies are focused on a 4-telescope version with an upgrade to a 6-telescope one. The instrument contains its own fringe tracker and tip-tilt control in order to reduce the constraints on the VLTI infrastructure and maximize the scientific return.

The discovery of extra-solar planets is one of the greatest achievements of modern astronomy. The detection of planets with a wide range of masses demonstrates that extra-solar planets of low mass exist. In this paper we describe a mission, called Darwin, whose primary goal is the search for, and characterization of, terrestrial extrasolar planets and the search for life. Accomplishing the mission objectives will require collaborative science across disciplines including astrophysics, planetary sciences, chemistry and microbiology. Darwin is designed to detect and perform spectroscopic analysis of rocky planets similar to the Earth at mid-infrared wavelengths (6 - 20 micron), where an advantageous contrast ratio between star and planet occurs. The baseline mission lasts 5 years and consists of approximately 200 individual target stars. Among these, 25 to 50 planetary systems can be studied spectroscopically, searching for gases such as CO2, H2O, CH4 and O3. Many of the key technologies required for the construction of Darwin have already been demonstrated and the remainder are estimated to be mature in the near future. Darwin is a mission that will ignite intense interest in both the research community and the wider public.

This paper is an H band interferometric examination of Arcturus, a star frequently used as a spatial and spectral calibrator. Using the IOTA 3 telescope interferometer, we performed spectro-interferometric observation (R~35) of Arcturus. Atmospheric models and prescriptions were fitted to the data to derive the brightness distribution of the photosphere. Image reconstruction was also obtained using two software algorithms: Wisard and Mira. An achromatic power law proved to be a good model of the brightness distribution, with a limb darkening compatible with the one derived from atmospheric model simulations using our Marcs model. A Rosseland diameter of 21.05+/-0.21 was derived, corresponding to an effective temperature of T_eff=4295+/-26 K. No companion was detected from the closure phases, with an upper limit on the brightness ratio of 8e-4 at 1AU. Dynamic range at such distance from the photosphere was established at 1.5e-4 (1sigma rms). An upper limit of 1.7e-3 was also derived for the level of brightness asymmetries present on the photosphere.

Because of atmospheric turbulence, obtaining high angular resolution images with a high dynamic range is difficult even in the near infrared domain of wavelengths. We propose a novel technique to overcome this issue. The fundamental idea is to apply techniques developed for long baseline interferometry to the case of a single-aperture telescope. The pupil of the telescope is broken down into coherent sub-apertures each feeding a single-mode fiber. A remapping of the exit pupil allows interfering all sub-apertures non-redundantly. A diffraction-limited image with very high dynamic range is reconstructed from the fringe pattern analysis with aperture synthesis techniques, free of speckle noise. The performances of the technique are demonstrated with simulations in the visible range with an 8 meter telescope. Raw dynamic ranges of 1:$10^6$ can be obtained in only a few tens of seconds of integration time for bright objects.

We introduce STECKMAP (STEllar Content and Kinematics via Maximum A Posteriori), a method to recover the kinematical properties of a galaxy simultaneously with its stellar content from integrated light spectra. It is an extension of STECMAP (astro-ph/0505209) to the general case where the velocity distribution of the underlying stars is also unknown. %and can be used as is for the analysis of large sets of data. The reconstructions of the stellar age distribution, the age-metallicity relation, and the Line-Of-Sight Velocity Distribution (LOSVD) are all non-parametric, i.e. no specific shape is assumed. The only a propri we use are positivity and the requirement that the solution is smooth enough. The smoothness parameter can be set by GCV according to the level of noise in the data in order to avoid overinterpretation. We use single stellar populations (SSP) from PEGASE-HR (R=10000, lambda lambda = 4000-6800 Angstrom, Le Borgne et al. 2004) to test the method through realistic simulations. Non-Gaussianities in LOSVDs are reliably recovered with SNR as low as 20 per 0.2 Angstrom pixel. It turns out that the recovery of the stellar content is not degraded by the simultaneous recovery of the kinematic distribution, so that the resolution in age and error estimates given in Ocvirk et al. 2005 remain appropriate when used with STECKMAP. We also explore the case of age-dependent kinematics (i.e. when each stellar component has its own LOSVD). We separate the bulge and disk components of an idealized simplified spiral galaxy in integrated light from high quality pseudo data (SNR=100 per pixel, R=10000), and constrain the kinematics (mean projected velocity, projected velocity dispersion) and age of both components.

This paper describes STECMAP (STEllar Content via Maximum A Posteriori), a flexible, non-parametric inversion method for the interpretation of the integrated light spectra of galaxies, based on synthetic spectra of single stellar populations (SSPs). We focus on the recovery of a galaxy's star formation history and stellar age-metallicity relation. We use the high resolution SSPs produced by PEGASE-HR to quantify the informational content of the wavelength range 4000 - 6800 Angstroms. A detailed investigation of the properties of the corresponding simplified linear problem is performed using singular value decomposition. It turns out to be a powerful tool for explaining and predicting the behaviour of the inversion. We provide means of quantifying the fundamental limitations of the problem considering the intrinsic properties of the SSPs in the spectral range of interest, as well as the noise in these models and in the data. We performed a systematic simulation campaign and found that, when the time elapsed between two bursts of star formation is larger than 0.8 dex, the properties of each episode can be constrained with a precision of 0.04 dex in age and 0.02 dex in metallicity from high quality data (R=10 000, signal-to-noise ratio SNR=100 per pixel), not taking model errors into account. The described methods and error estimates will be useful in the design and in the analysis of extragalactic spectroscopic surveys.

A general inversion technique for the recovery of the underlying distribution function for observed galactic disks is presented and illustrated. Under the assumption that these disks are axi-symmetric and thin, the proposed method yields the unique distribution compatible with all the observables available. The derivation may be carried out from the measurement of the azimuthal velocity distribution arising from positioning the slit of a spectrograph along the major axis of the galaxy. More generally, it may account for the simultaneous measurements of velocity distributions corresponding to slits presenting arbitrary orientations with respect to the major axis. The approach is non-parametric, i.e. it does not rely on a particular algebraic model for the distribution function. Special care is taken to account for the fraction of counter-rotating stars which strongly affects the stability of the disk. An optimisation algorithm is devised -- generalising the work of Skilling & Bryan (1984) -- to carry this truly two-dimensional ill-conditioned inversion efficiently. The performances of the overall inversion technique with respect to the noise level and truncation in the data set is investigated with simulated data. Reliable results are obtained up to a mean signal to noise ratio of~5 and when measurements are available up to $4 R_{e}$. A discussion of the residual biases involved in non parametric inversions is presented. Prospects of application to observed galaxies and other inversion problems are discussed.