arxiv: v2 [astro-ph.sr] 28 Jul 2010

Transcription

1 Astronomy & Astrophysics manuscript no. srenard 1491 accepted final c ESO 218 November 8, 218 arxiv:17.293v2 [astro-ph.sr] 28 Jul 21 Milli-arcsecond images of the Herbig Ae star HD S. Renard 1, F. Malbet 1, M. Benisty 2, E. Thiébaut 3, and J-P. Berger 4,1 1 Laboratoire d Astrophysique de Grenoble, CNRS-UJF UMR5571, BP 53, F-3841 Grenoble, France Stephanie.Renard, Fabien.Malbet@obs.ujf-grenoble.fr 2 INAF-Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 5125 Firenze, Italy benisty@arcetri.astro.it 3 Centre de Recherche Astrophysique de Lyon, CNRS-UCBL-ENSL UMR5574, F St-Genis-Laval, France thiebaut@obs.univ-lyon1.fr 4 European Southern Observatory, Alonso de Cordova, 317, Vitacura, Chile jpberger@eso.org Received ; accepted ABSTRACT Context. The very close environments of young stars are the hosts of fundamental physical processes, such as planet formation, stardisk 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. Aims. We gathered numerous interferometric measurements for the young star HD with various interferometers (, IOTA, KeckI and CHARA), allowing for the first time an image independent of any a priori model to be reconstructed. Methods. Using the Multi-aperture image Reconstruction Algorithm (MiRA), we reconstruct images of HD 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 Results. We obtain model-independent H and K-band images of the surroundings of HD The images present several significant features that we can relate to an inclined asymmetric flared disk around HD 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. Conclusions. 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. Key words. Instrumentation: interferometers - Techniques: image processing - Stars: pre-main sequence - Stars: individual: HD Introduction The process of star formation triggered by the collapse and fragmentation of a molecular cloud leads to the birth of a young star surrounded by a circumstellar disk and outflows. The disks are believed to be the place where the planets form (Boss 1997; Mayer et al. 22), and are composed of a mixture of gas and dust with a wide range of grain composition (Henning & Meeus 29). The standard picture of the close environment of premain sequence stars is so far limited to a quasi-stationary accreting disk partially reprocessing the irradiation from the central protostar with potential planetary gaps opened by newly formed planets. Some of the accreted material is ejected by means of bipolar outflows but their precise origin has not yet been identified. Most of the models are presently assumed to be symmetric around the star rotation axis and stationary. However hydrodynamical turbulence, gravitational waves, magneto-rotational instabilities, vortices and thermal instabilities occurring on the AU-scale are known to play a major role in star and planet formation (e.g., Balbus & Hawley 1991). The close environments of young stars are therefore not expected to be as simple as they are currently modeled, but the observational measurements do not provide tight enough constraints to unambiguously identify strong departures from axisymmetrical models. The photometric and spectroscopic observations obtained with very modest spatial resolution are usually integrated over a subarcsecond field of view which corresponds to several tens of AUs at the distance of the closest star formation regions (14pc). Observations at the scale of 1 AU and below correspond to 7 milli-arcseconds (mas) angular scale and requires therefore optical long-baseline interferometry. However, even if long-baseline optical interferometry is capable of reaching very high angular resolution, the observations are usually limited to a small number of measurements (see Millan-Gabet et al. 27, for a review). With the advent of interferometers with more than two telescopes and with significantly high spectral resolution (e.g., /AMBER, CHARA/MIRC), the number of interferometric observations has significantly increased allowing the first images to be reconstructed with aperture synthesis techniques. Images of stellar surfaces (e.g., Monnier et al. 27; Haubois et al. 29; Zhao et al. 29), binaries (e.g., Zhao et al. 28; Kraus et al. 29), or circumstellar shells around evolved stars (e.g., Le Bouquin et al. 29) have been obtained mostly for objects brighter than the brightest young stars. In the young stellar object field, we are at a comparable stage to that reached 4 years ago when the first 3 antennas of the VLA became operational producing the first radio-interferometry images (Hogg et al. 1969). We report the first attempt to reconstruct an image of a circumstellar disk around a young star using mostly AMBER/ interferometric data. In this study, we present reconstructed images of the young star HD (MWC 275), an isolated Herbig Ae star (HAe) of spectral type A1, with a 3 L luminosity, and a 2.3 M mass located at pc (van den Ancker et al. 1998; Natta et al. 24; Montesinos et al. 29). In scattered light (Grady et al. 2) and at millimeter wavelengths (Mannings & Sargent

2 2 S. Renard et al.: Milli-arcsecond images of the Herbig Ae star HD ), a disk has been detected on large scales, traced out to 54 AU. The CO millimeter line observations have revealed a large-scale inclined disk in Keplerian rotation probably evolving towards a debris disk phase (Isella et al. 27). HD also exhibits an asymmetric outflow perpendicular to the disk, with a chain of six Herbig-Haro knots (HH49) tracing the history of mass loss (Devine et al. 2; Wassell et al. 26). The emission of the innermost regions, observed in far-uv emission lines, have been attributed to optically thin gas accreting onto the stellar surface, a magnetically confined wind, or shocks at the base of the jet (Deleuil et al. 25; Swartz et al. 25). HD has been observed with several interferometers (see Benisty et al. 21, B1 hereafter). Using the largest set of interferometric data of a young star, we present here the first reconstructed images of a complex young stellar object. The article is organized as follow. Sect. 2 describes the image reconstruction method and the methodology employed to extract the best reconstructed image. Sect. 3 presents the reconstructed images obtained in the H and K bands and their analysis using simulated data generated from inner disk models. In Sect. 4, we describe the choices made during the image reconstruction process, and discuss the physical meaning of the features seen in the image, as well as the consequences for the models commonly used. Finally, Sect. 5 summarizes our results and provides some perspectives for the future. 2. Image reconstruction 2.1. Image reconstruction by MiRA The principle of interferometry is to interfere coherently the light coming from a single astronomical source from two or more independent telescopes (Lawson 2; Malbet & Perrin 27). An interferometer measures a complex number referred to as the visibility. According to the Van Cittert-Zernicke theorem, this complex visibility, V e i φ, is the Fourier transform of the object brightness distribution at the spatial frequency of the observations, given by the projected baseline in units of wavelength (B/λ). The visibility amplitude, V, is related to the spatial extent of the emission, while the phase, φ, provides the location of the photocenter. However, the requested infrastructure to carry out optical interferometry is complex and has led to a limitation in the number of telescopes (Baldwin & Haniff 22). The main consequence is to provide a sparse sampling of the spatial frequencies, the socalled (u, v) plane. Moreover, the absolute value of the phase φ is lost due to atmospheric turbulence that randomly modifies it. However, by adding the phases of the fringes measured for each baseline over a 3-telescope configuration, one can measure an additional quantity, the closure phase, which is insensitive to the atmospheric disturbance (Monnier 23; Monnier et al. 26). The closure phase includes part of the Fourier phase information, and is related to the global asymmetry of the emission: a point-symmetric object has a zero closure phase. The main observables are therefore the squared visibility amplitudes, V 2, and the closure phases (CP). The objective of the image reconstruction is to numerically retrieve an approximation of the true brightness distribution of the source given the interferometric measurements. To account for the data, the Fourier transform of the sought image should fit the measured complex visibilities. However, due to the sparse (u, v) coverage, the image reconstruction problem is ill-posed as there are more unknowns, e.g., the pixels of the image, than measurements. Additional prior constraints are therefore required to supplement the available data and retrieve a unique and stable solution. A very general solution is to define the optimal image to be the solution to the optimization problem (Thiébaut 25; Thiébaut & Giovannelli 29): min x { fdata (x) + µ f prior (x) } with x and x n = 1 (1) where x R N are the pixel values in the discretized image and N the number of pixels. In problem (1) the strict constraints x and n x n = 1 account for the non-negativity and the normalization of the brightness distribution. The objective function in the optimization problem (1) is a joint criterion of two components: Minimization of the likelihood term f data (x), which enforces agreement of the sought image with the data. In practice, this term is derived from the noise statistics. For instance, it is the χ 2 of the data for Gaussian statistics. Minimization of the regularization term f prior (x), which favors images that are the simplest or the smoothest according to a priori assumptions. The parameter µ > is used to tune the relative importance of the two terms so as to select the most regular image among all those compatible with the data. Many different algorithms have been developed to solve the image reconstruction problem from optical interferometry data (e.g., Cotton et al. 28; Thiébaut & Giovannelli 29). In this paper, we use the Multi-Aperture Image Reconstruction Algorithm (MiRA; by Thiébaut 28). MiRA is capable of dealing with any available interferometric data (complex visibilities, squared visibilities, closure phases, etc.), and has been successfully used to process real data (e.g., Lacour et al. 28; Lacour et al. 29; Le Bouquin et al. 29; Haubois et al. 29). When phases are missed because of the atmospheric turbulence, the MiRA algorithm can directly fit the available interferometric observables without explicitly rebuilding the missing phases. Finally, one can choose the most effective of various regularization methods integrated into MiRA for the type of object observed and check the effect of the regularization on the resulting image. MiRA directly attempts to solve the problem in Eq. (1) using an iterative non-linear optimization algorithm (Thiébaut 22). Because of the missing phases, the function to be minimized is not however convex thus has several local minima. The final image is therefore determined by the data, the choice of the regularization (and its level), and the initial image Methodology used in this work Systematic tests have been performed on the MiRA algorithm by Renard et al. (in prep. ; RTM1 hereafter) in which images of ten astrophysical objects were reconstructed, for different (u, v) coverages and signal-to-noise ratios. Twelve regularizations were tested and images were reconstructed for a set of weight factors µ. The optimal solution is the one that minimizes the meansquared distance between the model and the reconstructed image. The RTM1 tests led to the following conclusions: The total variation regularization (Strong & Chan 23), which minimizes the total norm of the image gradient, is the best regularization method in most of the cases. The weight factor µ depends, within one order of magnitude, on neither the amount of data, the signal-to-noise ratio, nor the object type, but only on the type of the regularization used. Each regularization has its best value for µ. n

3 S. Renard et al.: Milli-arcsecond images of the Herbig Ae star HD The chosen regularization for this work is therefore the total variation and the weight factor µ is set to be 1 2. The effect of changing µ is discussed in Appendix B. As explained above, because part of the phase information is lost, the image reconstruction from squared visibilities and closure phases yields to a criterion with multiple local minima. The solution depends on the initial image used to start the iterative process. Different starting images, such as a Dirac or a random flux distribution among pixels, were tested and, thanks to a good regularization term and a sufficient number of data points, the solutions obtained were always found to be the same. Here, we choose a symmetric Gaussian as a starting image. Once the parameters of the image reconstruction process are set, the images can be reconstructed. The solutions are not straightforward to analyze because of artifacts caused by the image reconstruction process or the quality of the data set, e.g., voids in the (u, v) plane, error bars. There are no objective criteria to distinguish between the actual structures from the object and the artifacts caused only by the data structure. We therefore performed a comparative analysis between our results and the results obtained from simulated data from the B1 model for HD To do so, we simulated fake data sets, using the B1 model and the same (u, v) plane and errors as in the real data set. Image reconstruction was also performed for the simulated data in the same conditions and compared to the image model. This comparative method is important to understand what could be trusted in the actual reconstructed images and what could not Interferometry data set The principal characteristics of the data set is summarized in this section. A more detailed description is given in the Appendix A. HD was observed with several interferometers. The large data set comprises H and K-band data from (Benisty et al. 21), IOTA (Monnier et al. 26), Keck-I (Monnier et al. 25), and CHARA (Tannirkulam et al. 28). This data set represents the largest set of interferometric data available so far for a young stellar object. Since a trade-off has to be found between enough data to pave the (u, v) plane and the wavelength dependency of the observed object, we decided to reconstruct two different images, one in H and one in K band, using all the spectral channels available in each band i.e. assuming the object to be grey in each band (see Sect. 4.2 for further discussion). 3. Results and first image analysis 3.1. Reconstructed images in the H and K bands Using the methodology described in Sect. 2.2, we reconstruct the images of HD in the H and K bands. To produce a rendu similar to that usually used in radio-interferometry, we convolve all the resulting images with a Gaussian beam at the interferometer resolution defined by the (u, v) plane. The resulting images are plotted in Fig. 1 with a color scale and in Fig. 2 with linear contours. As the minimum cut, we use the level corresponding to the expected dynamic range (see Sect. 4.1 for the discussion on how to compute this value). At first look, the spot representing the star is unambiguous and corresponds to the maximum of the images in the two bands. Around this central spot, one can see many secondary blobs. In the next section, we will discuss the level of confidence in these blobs by comparing them to the results obtained on simulated data from the B1 model. The main secondary blobs are concentrated in the center of the image around the brightest spot at a distance of smaller than 4-5 mas. To show the location of these secondary blobs, an ellipse is drawn using a dashed blue line in Figs. 1 and 2 with a semi-major axis of 4.5 mas, an inclination of 55, and a position angle of 155. When considering their repartition along the ellipse, we find that they are less numerous at the bottom of the ellipse. The percentage of flux in these blobs is around 3% in the K band and 24% in the H band. Inside the ellipse, the central spot is not point-like. In the H band, the central spot is extended, while in the K band the energy is spread into two close separated spots. The emission does not decrease slowly from the central spot towards the exterior but instead shows a rapid decay before increasing again when crossing the ellipse to finally decrease at large distances. We also note that the ellipse and the central spots are not exactly centered, although this result may not be relevant Comparison with simulated images from models To analyze the artifacts in the images, we use the model of HD presented in B1 (see the bottom right squares of Fig. 3 for the H and K-band models) to simulate images with the same conditions as the actual ones. The model is composed of a star (producing 3 % and 14 % of the flux in the H and K bands, respectively, estimated from spectral energy distribution fitting), a dust rim located at.45 AU ( 3.6 mas) representing 16 % (36 %) in the H (K) band and a bright inner disk, from.1 to.45 AU, contributes to the remaining emission. The reconstructed images of the model from simulated data with the same (u, v) coverage and the same error bars as the real data are shown in Fig. 3 for the H and K bands. The analysis of the reconstructed images of the model in the K band illustrates that the following structures are well retrieved by the image reconstruction process: The dust rim, which is clearly visible at the right location and appears as a somewhat blobby ellipse. We checked that by changing the (u, v) plane filling, the blob location changes but remains aligned along the ellipse. The energy inside the disk, which is more spread than if there was only the star inside the dust rim in the model (see below). This emission between the star and the rim represents the bright inner disk. The skewness of the dust rim is visible in the reconstructed images: the bottom part of the blobby ellipse has less flux than the top part. The star, the dust rim, and the inner disk provide 15%, 3%, and 55% of the flux, respectively. We emphasize that these values are in close agreement with the model. The only structure that is not retrieved in the image is the hole between the star and the bright inner disk, inside.1 AU, which is too small to be resolved by the interferometer. In the H band, the dust rim disappears in the reconstructed image as a ring of blobs, but still seems to define the outer boundary of the object. The main reasons are that the rim represents only 16% of the flux to be compared with 36% in the K band, and, that the angular resolution is not high enough in the H band data compared to the K band data (with the CHARA very long baselines). For the dust rim to be seen in the H band, we would need data on longer baselines and higher dynamics in the reconstructed image (1 at least, 2 to be unambiguously seen). As in the K band, the bright inner disk is also present in the H band image as a large spot in the middle, which would not have existed with the star only (see below), and represents 86% of the total flux.

4 4 S. Renard et al.: Milli-arcsecond images of the Herbig Ae star HD Fig. 1. Reconstructed images of HD in the H (left) and K bands (right), after a convolution with a Gaussian beam at the interferometer resolution. The colors are scaled to the squared root of the intensity with a cut corresponding to the maximum expected dynamic range (see text for details). The blue ellipse traces the location of the main secondary blobs, and the green dotdashed ellipse corresponds to the location of the rim in the B1 model, with its width given by the green dashed ellipses. North is up and east is left. The sub-panel in the right corner of each plot indicates the Gaussian beam at the interferometer resolution, applicable to Figs. 3, 4, and 6. H band K band 1 1 relative δ (mas) relative δ (mas) relative α (mas) relative α (mas) Fig. 2. Contours of the reconstructed images of HD in the H (left) and K bands (right), after a convolution with a Gaussian beam at the interferometer resolution. The contours vary linearly between the minimum cut corresponding to the maximum expected dynamic range and the image maximum, with a step around.1. To demonstrate that the bright inner disk is clearly seen in the reconstructed images, i.e., that the central spot includes more energy than that from the star alone, reconstruction of a simpler model is performed. This model is the same as the B1 model but without the bright inner disk, i.e. a star surrounded by a Gaussian ring. The star fluxes in the H and K bands remain the same and the Gaussian ring accounts for 7% of the flux in the H band and 86% in the K band. Figure 4 clearly indicates that the star alone does not spread across more than over 4 pixels in the K band and 7 in the H band, which is less than in Fig. 3. This analysis performed on existing models allows us to state which features in the reconstructed images from Fig. 1 can be trusted. We argue that the main secondary blobs present around the main central spot are real. Their spatial distribution along an ellipse and the intensity present between these peaks and the central spot are also real. However, the clumpy structure of the ring is probably not representative of the reality, but only of the actual (u, v) plane. More observations at different spatial frequencies will probably change the actual position of these peaks along the ellipse, which may be smoothed. However, we conclude that the inclination and orientation of the observed distribution of peaks along an ellipse are real. This orientation and inclination are indeed very close to the ones fitted by B1 and Tannirkulam et al. (28), and are consistent with previous estimates at different wavelengths (Isella et al. 27). The second feature that we think is representative of the reconstructed image is that the central spot is extended and not reduced to an unresolved point as a point-like star would be. The shape of this central spot is certainly dependent on the filling of the (u, v) plane, although the position of the centroid is certainly representative of reality.

5 S. Renard et al.: Milli-arcsecond images of the Herbig Ae star HD Fig. 3. Reconstructed images of the B1 model of HD in the H (left) and K bands (right). The dashed green ellipse corresponds to the location of the rim in this model. The models used are presented in the upper left corner. Same conventions as in Fig. 1. Fig. 4. Reconstructed images of a geometrical model of HD with only a star plus a Gaussian ring in the H (left) and K bands (right). The dashed green ellipse corresponds to the location of the rim in this model. The models are presented in the upper left corner. Same conventions as in Fig Discussion In this section, we discuss the reconstructed images Dynamic range To compute the theoretical dynamical range of our image, we use an estimator based on the one proposed by Baldwin & Haniff (22) r n, (2) 2 (δv/v) + (δcp)2 where n is the total number of measurements, δv/v is the relative error in the visibilities, and δcp is the error in the closure phases (in radian). This number indicates the maximum contrast that can be reached in the image given the data, i.e. the ratio of the maximum of the image to the minimum value that can be trusted. Applying Eq. (2) to the data, a dynamic range of 78 is found in the K band. In the H band, for which there is fewer data and the error bars are slightly larger, a dynamic range of 4 is com- puted. Minimum cut levels of 1/78 and 1/4, respectively, are applied to all the K and H figures Use of the spectral information The decision to present one reconstructed image in the H and one in the K band results from different tests made on the B1 model. A trade-off has to be made between ensuring that we have enough data to pave the (u, v) plane and the wavelength dependency expected from circumstellar disks. Because of the intrinsic chromaticity of the object, we prefer to reconstruct two separated images in the H and K bands, otherwise two separated visibilities, sampling different emitting regions, may correspond to the same spatial frequency. We show in Fig. 5 the combination of both reconstructed images that illustrate complementary features. For each band, we used all the data points in the different spectral channels, assuming implicitly the object to be grey in each separated band. Compared to the model, this method pro-

6 6 S. Renard et al.: Milli-arcsecond images of the Herbig Ae star HD Fig. 5. Combination of the reconstructed images of HD in a two-color image (H band in green, K band in red). The blue ellipse traces the location of the main secondary blobs of the K- band emission. The sub-panels indicate the Gaussian beam at the interferometer resolution used in the convolution. duces images with fewer artifacts than when reconstructing a single broadband image with the average data for all spectral channels. The (u, v) plane is indeed far more filled taking into account all the wavelengths, because there are more spatial frequencies Physical consequences for the models We emphasize that we have obtained a new type of data in the form of images reconstructed with limited assumptions. The analysis of the results described in Sect. 3.2 is useful for distinguishing artifacts caused mainly by the shape of the (u, v) plane from what we infer to be true features. In this section, we highlight the new information provided by these images, but are also aware that these new pieces of evidence have to be handled with great care. The ellipse described by the successive blobs in the reconstructed image at a distance of 4-5 mas from the center of the image and underlined by a blue dot-dashed line in Figs. 1 and 2 certainly traces an external ring. The characteristics of this ring are not exactly the same as the rim found in the B1 model plotted with a green dashed line in Figs. 1 and 2. The radius of the ring is.55 AU instead of.45 ±.5 AU in the B1 model, the inclination is 55 instead of 48 ± 2, and the position angle is 155 instead of 136 ± 2. The intensity is close to that of the model with 3% of the total flux instead of 36%. The location of the blobs differs a little bit between the H and K bands: those in the H-band image are closer to the center of the image than those in the K-band image. This behavior could be explained by a temperature gradient in the disk, which has a tendency to move the peak of intensity closer to the center for the shorter wavelengths. We probably also need data on longer baselines (equivalent to the K-band CHARA data) to resolve unambiguously the external ring in the H-band image. This bright ring is not clear evidence of a physical rim, which was proposed by Dullemond et al. (21). We propose that this feature instead traces an enhancement of the intensity caused by a change in the opacity in the disk probably due to sublimation of dust. We assume that the number of blobs along this ellipse is representative of the azimuthal distribution of the intensity. Following B1 who in their analysis proposed a model of a rim with a skewed distribution of intensity along the ring (see Figs. 4), we assume that the distribution of light along the blue ellipse indicates that the actual ring of light is less luminous in the south-west direction than in the north-east direction at least in the K-band image. In the H-band image, there are blobs around the central spot that do not appear in the reconstructed image of the model. Does this mean that there is more flux in the external ring than expected? Is the skewness factor more important because the blobs appeared at one side only? Does the presence of blobs in the south part of the ellipse rather than in the north in the K-band image have some significance? To remain on solid ground, we assume that the intensity varies along the distribution of blobs, but do not have definitive data to determine the magnitude of this effect. The departure from the axisymmetry may be caused by an inclined surface of the disk or even a strong dust opacity change. The energy detected within the ellipse is certainly real: as explained in Sect. 3.2, the central region, which contains about 7% of the flux, certainly does not originate in a single unresolved star but from an extended source. These results independently confirm the conclusions of B1, without using any model. Indeed, B1 were unable to fit the visibilities at higher spatial frequencies without introducing some continuous emission in the space between the rim and the star. They called this region the inner disk. We do not know whether the shape of this central source in the K-band image in the form of two spots is real or not, since in the simulated image this central source also seems to be decomposed into 3 single sources, which were not present in the model. It might be only the effect of the (u, v) plane coverage, but we cannot exclude too that it might be due to a hot spot in the disk, although we should in principle then see it in the two images. The reconstructed images cannot help us to determine the origin of the inner disk emission, and ascertain whether it comes from hot gas or very refractory grains. These unknowns may be solved by combining high resolution spectroscopic observations in the near-infrared with advanced models that selfconsistently compute the emission of both dust and gas. Finally, it remains unclear whether the non-zero value of the closure phases found by B1 is caused by the contribution from the inner disk not being exactly centered with the external ring. This might be the case when the disk surface is flaring and the system is seen with non-zero inclination. If the curvature of the surface of the disk probed by our images is large enough, then the ellipses tracing equal distances to the star will be shifted in the polar direction Consequences on the image reconstruction Since image reconstruction in optical interferometry remains in its infancy, only a too sparse (u, v) coverage is available to reconstruct an unambiguous image of a complex object and analyze it without the help of the model-fitting technique. Repeated comparisons between the model and the reconstructed image have to be performed to avoid over-interpreting the structures in the images. However, in all these cases, the image reconstruction technique remains the only technique able to perform modelindependent analysis of the data and is a powerful tool to give

7 S. Renard et al.: Milli-arcsecond images of the Herbig Ae star HD Fig. 6. Reconstructed image of the B1 model in the K band, using a synthetic (u, v) plane obtained with 3 quadruplets (A-K- G-I1, D-H-G1-I1, E-G-H-I1) at and the CHARA baselines (S1-W1, W1-W2, S2-W2, E1-W1, E2-S2, S2-W1). more credits to the models, derive tighter constraints, or even reveal unexpected structures. In the future, two aspects of the image reconstruction will be important: 1. The homogeneity of the (u, v) plane: thanks to several tests on the B1 model, we found that the global blobby aspect of the reconstructed image is almost certainly caused by the non-homogeneity of the (u, v) plane. Indeed the reconstructed images from a homogeneous (u, v) plane are smoother. The holes in the (u, v) plane correspond to blobby or point-like structures in the reconstructed images, and the baselines for the observations have to be carefully chosen to map the (u, v) plane as homogeneously as possible. 2. The quality of the measurements: the number of data points and their associated errors can clearly determine the dynamic range of the image. A larger number of data points and smaller error bars are needed to improve the dynamics of the image. These aspects are illustrated in Fig. 6 where data on the B1 model are simulated with a (u, v) plane for 3 quadruplets of telescopes at, as offered in the coming period P86 for the AMBER instrument and later for the PIONIER visitor instrument. The reconstructed image from these data contains fewer artifacts than those from the actual data set. The advent of the next generation imaging instruments at the, such as GRAVITY, MATISSE, and VSI, with up to 8 telescopes, will allow us to study in detail the inner disk structure, to detect planetary gaps and map outflows bringing unambiguous constraints on crucial mechanisms for star and planet formation. 5. Conclusion Since the renewal of optical interferometry in the mid-7 s, imaging using aperture synthesis between different telescopes has been attempted for all possible astrophysical topics following the path opened by radio-interferometry with the Very Large Array. Unfortunately, optical interferometry requires precisions in the optical path delay 1 3 to 1 5 smaller than in the radio domain and therefore the development of large optical arrays has been slower. Although images of stellar surfaces or binaries have been obtained, no synthetic images have so far been possible at the milli-arcsecond scale of complex young stellar objects. This paper is the first one attempting to reconstruct images of the close environment of a young star independent of any a priori model by gathering many interferometric measurements on the young star HD with various interferometers (, IOTA, KeckI, and CHARA). Using the Multi-Aperture Image Reconstruction Algorithm (MiRA), we have reconstructed images of HD in the H and K bands. The quality of the images compares well to those obtained 4 years ago with the first three antennas of the Very Large Array. To assess the reality of the features that are present, we have compared these images with images obtained from simulated data using the best known physical model of the environment of HD and reconstructed with the same conditions of (u, v) plane and noise. Those model-independent H and K-band images of the surrounding of HD reveal several significant features that we can relate to the presence of an inclined asymmetric disk around HD We also detect a pre-eminent intensity at the location of dust sublimation above a structure extending from the central source with some differential effect along the azimuth that is not necessarily related to a puffed-up inner rim. Together with a slight offset of the central source compared to this bright ring, we proposed this to be the signature of the flared surface of a disk with a discontinuity caused by a opacity change. These images confirm that the morphology of the close environment of young stars is more complex than the simple models used in the literature. We have also shown that having a more uniform (u, v) plane coverage and higher measurement accuracies should help us to obtain higher quality images in the future. Acknowledgements. The authors are grateful to the different institutions that trusted them in giving them guaranteed observing time with the to perform the first image of a complex young stellar object: CNRS through its national programs (ASHRA and PNPS), INAF and the AMBER consortium. This research has made use of NASA s Astrophysics Data System service, of the Jean-Marie Mariotti Center (JMMC) resources of the SIMBAD database, operated at CDS, Strasbourg, France, and of the Yorick, a free data processing language written by D. Munro ( M.B. acknowledges fundings from INAF (grant ASI-INAF I/16/7/). References Balbus, S. A. & Hawley, J. F. 1991, ApJ, 376, 214 Baldwin, J. E. & Haniff, C. A. 22, Royal Society of London Philosophical Transactions Series A, 36, 969 Benisty, M., Natta, A., Isella, A., et al. 21, A&A, 511, A74+ Boss, A. P. 1997, Science, 276, 1836 Cotton, W., Monnier, J., Baron, F., et al. 28, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 713, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series Deleuil, M., Bouret, J., Catala, C., et al. 25, A&A, 429, 247 Devine, D., Grady, C. A., Kimble, R. A., et al. 2, ApJ, 542, L115 Dullemond, C. P., Dominik, C., & Natta, A. 21, ApJ, 56, 957 Grady, C. A., Devine, D., Woodgate, B., et al. 2, ApJ, 544, 895 Haubois, X., Perrin, G., Lacour, S., et al. 29, A&A, 58, 923 Henning, T. & Meeus, G. 29, in Physical Processes in Circumstellar Disks around Young Stars, Garcia, PJV (Ed.), Theoretical Astrophysics Series, Chicago University Press, Vol. in press (arxiv:911.11) Hogg, D. E., MacDonald, G. H., Conway, R. G., & Wade, C. M. 1969, AJ, 74, 126 Isella, A., Testi, L., Natta, A., et al. 27, A&A, 469, 213 Kraus, S., Weigelt, G., Balega, Y. Y., et al. 29, A&A, 497, 195 Lacour, S., Meimon, S., Thiébaut, E., et al. 28, A&A, 485, 561 Lacour, S., Thibaut, E., Perrin, G., et al. 29, ApJ, 77, 632 Lawson, P. R., ed. 2, Principles of Long Baseline Stellar Interferometry Le Bouquin, J., Lacour, S., Renard, S., et al. 29, A&A, 496, L1 Malbet, F. & Perrin, G. 27, New Astronomy Review, 51, 563 Mannings, V. & Sargent, A. I. 1997, ApJ, 49, 792

8 8 S. Renard et al.: Milli-arcsecond images of the Herbig Ae star HD Mayer, L., Quinn, T., Wadsley, J., & Stadel, J. 22, Science, 298, 1756 Millan-Gabet, R., Malbet, F., Akeson, R., et al. 27, Protostars and Planets V, 539 Monnier, J. D. 23, Reports on Progress in Physics, 66, 789 Monnier, J. D., Berger, J., Millan-Gabet, R., et al. 26, ApJ, 647, 444 Monnier, J. D., Millan-Gabet, R., Billmeier, R., et al. 25, ApJ, 624, 832 Monnier, J. D., Zhao, M., Pedretti, E., et al. 27, Science, 317, 342 Montesinos, B., Eiroa, C., Mora, A., & Merín, B. 29, A&A, 495, 91 Natta, A., Testi, L., Neri, R., Shepherd, D. S., & Wilner, D. J. 24, A&A, 416, 179 Strong, D. & Chan, T. 23, Inverse Problems, 19, S165 Swartz, D. A., Drake, J. J., Elsner, R. F., et al. 25, ApJ, 628, 811 Tannirkulam, A., Monnier, J. D., Millan-Gabet, R., et al. 28, ApJ, 677, L51 Thiébaut, E. 22, in Presented at the Society of Photo-Optical Instrumentation Engineers (SPIE) Conference, Vol. 4847, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, ed. J.-L. Starck & F. D. Murtagh, Thiébaut, E. 25, in NATO ASIB Proc. 198: Optics in astrophysics, ed. R. Foy & F. C. Foy, Thiébaut, E. 28, in Presented at the Society of Photo-Optical Instrumentation Engineers (SPIE) Conference, Vol. 713, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series Thiébaut, É. & Giovannelli, J.-F. 29, ArXiv e-prints van den Ancker, M. E., de Winter, D., & Tjin A Djie, H. R. E. 1998, A&A, 33, 145 Wassell, E. J., Grady, C. A., Woodgate, B., Kimble, R. A., & Bruhweiler, F. C. 26, ApJ, 65, 985 Zhao, M., Gies, D., Monnier, J. D., et al. 28, ApJ, 684, L95 Zhao, M., Monnier, J. D., Pedretti, E., et al. 29, ApJ, 71, 29

9 Table A.1. Log of the data used for the image reconstruction. S. Renard et al.: Milli-arcsecond images of the Herbig Ae star HD Interferom. Obs. B P (m) # meas. Band Disp. Refs. KI K BB M5 IOTA H BB M6 CHARA K BB T8 / K 35 B1 AMBER 544 H Notes: B P is the projected baselines; Disp. is the spectral dispersion; Refs. are M5 (Monnier et al. 25), M6 (Monnier et al. 26), T8 (Tannirkulam et al. 28), B1 (Benisty et al. 21) Appendix A: The data set in detail The data used in this article comes from several interferometers. The name of these facilities, the observation date, the projected baselines, the number of measurements, the spectral dispersion, and the references have been summarized in Table A.1. The (u, v) plane in the H and K bands is plotted in Fig. A.1. The data (squared visibilities and closure phases) is displayed in Fig. A.2. For a better view, a zoom on the data is shown in Fig. A.3, with various colors for different ranges of baseline position angles. Appendix B: The effect of the regularization weight The systematic tests performed in RTM1 shows that a weight factor µ can be associated with a regularization term to within an order of magnitude. In this appendix, we illustrate the effect of the µ factor on the reconstructed images, therefore determining the influence of µ and demonstrating that µ = 1 2 is the best value. The upper part of Fig. B.1 presents the reconstructed images from the B1 model in the K band for 3 different values of the weight factor µ (1, 1, 1). When comparing them to the model (see the upper left corner of Fig. 3), the following information can be extracted: In the left part of Fig. B.1, the image is too regularized, meaning that too much weight is put on the regularization term and not enough on the data. The reconstructed image is too smooth and several structure are not visible, because of the lack of distinction between the external ring and the internal disk. In the right part of Fig. B.1, the image is insufficiently regularized, meaning that too much weight has been placed on the data. The reconstructed image is far more blobby than the original and the flux in the internal disk starts to disappear. The ideal weight factor is shown in the middle part of Fig. B.1. In this figure, all the different characteristics of the model are presented and it seems to provide the best reconstruction between the three. The same analysis can be performed on the reconstructed images from the real data and are illustrated in the bottom part of Fig. B.1.

10 1 S. Renard et al.: Milli-arcsecond images of the Herbig Ae star HD H band K band 5 IOTA 1 KeckI CHARA V/λ [m/µ m] V/λ [m/µ m] U/λ [m/µ m] U/λ [m/µ m] Fig. A.1. (u, v) plane coverage of the data used for the image reconstruction in spatial frequencies for the H (left) and K (right) bands. The different interferometers are plotted in different colors and symbols. 1. H band IOTA 1. K band Keck CHARA Squared Visibility V 2.5 Squared Visibility V Spatial frequency B/λ [m/µ m] H band Spatial frequency B/λ [m/µ m] K band IOTA 5 5 Closure Phase (deg) Closure Phase (deg) Maximum spatial frequency [m/µ m] Maximum spatial frequency [m/µ m] Fig. A.2. Squared visibilities (up) and closure phases (bottom) in the H (left) and K (right) bands. The different interferometers are plotted in different colors and symbols.

11 S. Renard et al.: Milli-arcsecond images of the Herbig Ae star HD H band < PA B < 3 1. K band < PA B < 3 3 < PA B < 6 3 < PA B < 6 6 < PA B < 9 6 < PA B < 9 Squared Visibility V < PA B < < PA B < < PA B < 18 Squared Visibility V < PA B < < PA B < < PA B < Spatial frequency B/λ [m/µ m] Spatial frequency B/λ [m/µ m] Fig. A.3. Squared visibilities in the H (left) and K (right) bands for the /AMBER data. The different colors and symbols correspond to different baseline position angle ranges. Fig. B.1. Reconstructed images of HD in the K band from simulated data of the B1 model (up) and from the real data (bottom), for 3 different values of the weight factor µ. conventions as in Fig. 1.