The number of publications of aperture-synthesis images based on optical long-baseline interferometry measurements has recently increased due to easier access to visible and infrared interferometers. The interferometry technique has now reached a technical maturity level that opens new avenues for numerous astrophysical topics requiring milli-arcsecond model-independent imaging. In writing this paper our motivation was twofold: 1) review and publicize emblematic excerpts of the impressive corpus accumulated in the field of optical interferometry image reconstruction; 2) discuss future prospects for this technique by selecting four representative astrophysical science cases in order to review the potential benefits of using optical long baseline interferometers. For this second goal we have simulated interferometric data from those selected astrophysical environments and used state-of-the-art codes to provide the reconstructed images that are reachable with current or soon-to-be facilities. The image reconstruction process was "blind" in the sense that reconstructors had no knowledge of the input brightness distributions. We discuss the impact of optical interferometry in those four astrophysical fields. We show that image reconstruction software successfully provides accurate morphological information on a variety of astrophysical topics and review the current strengths and weaknesses of such reconstructions. We investigate how to improve image reconstruction and the quality of the image possibly by upgrading the current facilities. We finally argue that optical interferometers and their corresponding instrumentation, existing or to come, with 6 to 10 telescopes, should be well suited to provide images of complex sceneries.

With the advent of infrared long-baseline interferometers with more than two telescopes, both the size and the completeness of interferometric data sets have significantly increased, allowing images based on models with no a priori assumptions to be reconstructed. Our main objective is to analyze the multiple parameters of the image reconstruction process with particular attention to the regularization term and the study of their behavior in different situations. The secondary goal is to derive practical rules for the users. Using the Multi-aperture image Reconstruction Algorithm (MiRA), we performed multiple systematic tests, analyzing 11 regularization terms commonly used. The tests are made on different astrophysical objects, different (u,v) plane coverages and several signal-to-noise ratios to determine the minimal configuration needed to reconstruct an image. We establish a methodology and we introduce the mean-square errors (MSE) to discuss the results. From the ~24000 simulations performed for the benchmarking of image reconstruction with MiRA, we are able to classify the different regularizations in the context of the observations. We find typical values of the regularization weight. A minimal (u,v) coverage is required to reconstruct an acceptable image, whereas no limits are found for the studied values of the signal-to-noise ratio. We also show that super-resolution can be achieved with increasing performance with the (u,v) coverage filling. Using image reconstruction with a sufficient (u,v) coverage is shown to be reliable. The choice of the main parameters of the reconstruction is tightly constrained. We recommend that efforts to develop interferometric infrastructures should first concentrate on the number of telescopes to combine, and secondly on improving the accuracy and sensitivity of the arrays.

Circumstellar disks surrounding young stars are known to be the birthplaces of planets, and the innermost astronomical unit is of particular interest. We present new long-baseline spectro-interferometric observations of the HerbigAe star, HR5999, obtained in the H and K bands with the AMBER instrument at the VLTI, and aim to produce near-infrared images at the sub-AU spatial scale. We spatially resolve the circumstellar material and reconstruct images using the MiRA algorithm. In addition, we interpret the interferometric observations using models that assume that the near-infrared excess is dominated by the emission of a circumstellar disk. We compare the images reconstructed from the VLTI measurements to images obtained using simulated model data. The K-band image reveals three main elements: a ring-like feature located at ~0.65 AU, a low surface brightness region inside, and a central spot. At the maximum angular resolution of our observations (1.3 mas), the ring is resolved while the central spot is only marginally resolved, preventing us from revealing the exact morphology of the circumstellar environment. We suggest that the ring traces silicate condensation, i.e., an opacity change, in a circumstellar disk around HR 5999. We build a model that includes a ring at the silicate sublimation radius and an inner disk of low surface brightness responsible for a large amount of the near-infrared continuum emission. The model successfully fits the SED, visibilities, and closure phases, and provides evidence of a low surface brightness region inside the silicate sublimation radius. This study provides additional evidence that in HerbigAe stars, there is material in a low surface brightness region, probably a low optical depth region, located inside the silicate sublimation radius and of unknown nature.

Young and close multiple systems are unique laboratories to probe the initial dynamical interactions between forming stellar systems and their dust and gas environment. Their study is a key building block to understanding the high frequency of main-sequence multiple systems. However, the number of detected spectroscopic young multiple systems that allow dynamical studies is limited. GW Orionis is one such system. It is one of the brightest young T Tauri stars and is surrounded by a massive disk. Our goal is to probe the GW Orionis multiplicity at angular scales at which we can spatially resolve the orbit. We used the IOTA/IONIC3 interferometer to probe the environment of GW Orionis with an astronomical unit resolution in 2003, 2004, and 2005. By measuring squared visibilities and closure phases with a good UV coverage we carry out the first image reconstruction of GW Ori from infrared long-baseline interferometry. We obtain the first infrared image of a T Tauri multiple system with astronomical unit resolution. We show that GW Orionis is a triple system, resolve for the first time the previously known inner pair (separation $\rho\sim$1.4 AU) and reveal a new more distant component (GW Ori C) with a projected separation of $\sim$8 AU with direct evidence of motion. Furthermore, the nearly equal (2:1) H-band flux ratio of the inner components suggests that either GW Ori B is undergoing a preferential accretion event that increases its disk luminosity or that the estimate of the masses has to be revisited in favour of a more equal mass-ratio system that is seen at lower inclination. Accretion disk models of GW Ori will need to be completely reconsidered because of this outer companion C and the unexpected brightness of companion B.

The very close environments of young stars are the hosts of fundamental physical processes, such as planet formation, star-disk interactions, mass accretion, and ejection. The complex morphological structure of these environments has been confirmed by the now quite rich data sets obtained for a few objects by near-infrared long-baseline interferometry. We gathered numerous interferometric measurements for the young star HD163296 with various interferometers (VLTI, IOTA, KeckI and CHARA), allowing for the first time an image independent of any a priori model to be reconstructed. Using the Multi-aperture image Reconstruction Algorithm (MiRA), we reconstruct images of HD 163296 in the H and K bands. We compare these images with reconstructed images obtained from simulated data using a physical model of the environment of HD 163296. We obtain model-independent $H$ and $K$-band images of the surroundings of HD 163296. The images present several significant features that we can relate to an inclined asymmetric flared disk around HD 163296 with the strongest intensity at about 4-5 mas. Because of the incomplete spatial frequency coverage, we cannot state whether each of them individually is peculiar in any way. For the first time, milli-arcsecond images of the environment of a young star are produced. These images confirm that the morphology of the close environment of young stars is more complex than the simple models used in the literature so far.

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.

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.

In this paper, we study the feasibility of obtaining near-infrared spectra of bright extrasolar planets with the 2nd generation VLTI Spectro-Imager instrument (VSI), which has the required angular resolution to resolve nearby hot Extrasolar Giant Planets (EGPs) from their host stars. Taking into account fundamental noises, we simulate closure phase measurements of several extrasolar systems using four 8-m telescopes at the VLT and a low spectral resolution (R = 100). Synthetic planetary spectra from T. Barman are used as an input. Standard chi2-fitting methods are then used to reconstruct planetary spectra from the simulated data. These simulations show that low-resolution spectra in the H and K bands can be retrieved with a good fidelity for half a dozen targets in a reasonable observing time (about 10 hours, spread over a few nights). Such observations would strongly constrain the planetary temperature and albedo, the energy redistribution mechanisms, as well as the chemical composition of their atmospheres. Systematic errors, not included in our simulations, could be a serious limitation to these performance estimations. The use of integrated optics is however expected to provide the required instrumental stability (around 10^-4 on the closure phase) to enable the first thorough characterisation of extrasolar planetary emission spectra in the near-infrared.