Showing 1–50 of 64 results for author: Thiebaut, E

Shrinkage MMSE estimators of covariances beyond the zero-mean and stationary variance assumptions

Authors: Olivier Flasseur, Eric Thiébaut, Loïc Denis, Maud Langlois

Abstract: We tackle covariance estimation in low-sample scenarios, employing a structured covariance matrix with shrinkage methods. These involve convexly combining a low-bias/high-variance empirical estimate with a biased regularization estimator, striking a bias-variance trade-off. Literature provides optimal settings of the regularization amount through risk minimization between the true covariance and i… ▽ More We tackle covariance estimation in low-sample scenarios, employing a structured covariance matrix with shrinkage methods. These involve convexly combining a low-bias/high-variance empirical estimate with a biased regularization estimator, striking a bias-variance trade-off. Literature provides optimal settings of the regularization amount through risk minimization between the true covariance and its shrunk counterpart. Such estimators were derived for zero-mean statistics with i.i.d. diagonal regularization matrices accounting for the average sample variance solely. We extend these results to regularization matrices accounting for the sample variances both for centered and non-centered samples. In the latter case, the empirical estimate of the true mean is incorporated into our shrinkage estimators. Introducing confidence weights into the statistics also enhance estimator robustness against outliers. We compare our estimators to other shrinkage methods both on numerical simulations and on real data to solve a detection problem in astronomy. △ Less

Submitted 26 June, 2024; v1 submitted 11 March, 2024; originally announced March 2024.

Comments: Accepted to EUSIPCO 2024

Closing the loop as an inverse problem: the real-time control of THEMIS adaptive optics

Authors: Éric Thiébaut, Michel Tallon, Isabelle Tallon-Bosc, Bernard Gelly, Richard Douet, Maud Langlois, Gil Moretto

Abstract: We have taken advantage of the implementation of an adaptive optics system on the Themis solar telescope to implement innovative strategies based on an inverse problem formulation for the control loop. Such an approach encompassing the whole system implies the estimation of the pixel variances of the Shack-Hartmann wavefront sensor, a novel real-time method to extract the wavefront slopes as well… ▽ More We have taken advantage of the implementation of an adaptive optics system on the Themis solar telescope to implement innovative strategies based on an inverse problem formulation for the control loop. Such an approach encompassing the whole system implies the estimation of the pixel variances of the Shack-Hartmann wavefront sensor, a novel real-time method to extract the wavefront slopes as well as their associated noise covariance, and the computation of pseudo-open loop data. The optimal commands are computed by iteratively solving a regularized inverse problem with spatio-temporal constraints including Kolmogorov statistics. The latency of the dedicated real-time control software with conventional CPU is shorter than 300 $μ$s from the acquisition of the raw 400 x 400 pixel wavefront sensor image to the sending of the commands. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Journal ref: Adaptive Optics Systems VIII, Jul 2022, Montr{é}al, Canada. pp.1218507

Beyond FRiM, ASAP: a family of sparse approximation for covariance matrices and preconditioners

Authors: Éric Thiébaut, Michel Tallon, Samuel Thé, Loïc Denis

Abstract: The FRiM fractal operator belongs to a family of operators, called ASAP, defined by an ordered selection of nearest neighbors. This generalization provides means to improve upon the good properties of FRiM. We propose a fast algorithm to build an ASAP operator mimicking the fractal structure of FRiM for pupils of any size and geometry and to learn the sparse coefficients from empirical data. We em… ▽ More The FRiM fractal operator belongs to a family of operators, called ASAP, defined by an ordered selection of nearest neighbors. This generalization provides means to improve upon the good properties of FRiM. We propose a fast algorithm to build an ASAP operator mimicking the fractal structure of FRiM for pupils of any size and geometry and to learn the sparse coefficients from empirical data. We empirically show the good approximation by ASAP of correlated statistics and the benefits of ASAP for solving phase restoration problems. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Journal ref: Adaptive Optics Systems VIII, Jul 2022, Montr{é}al, Canada. pp.121852Q

PACOME: Optimal multi-epoch combination of direct imaging observations for joint exoplanet detection and orbit estimation

Authors: Jules Dallant, Maud Langlois, Olivier Flasseur, Éric Thiébaut

Abstract: Exoplanet detections and characterizations via direct imaging require high contrast and high angular resolution. These requirements typically require (i) cutting-edge instrumental facilities, (ii) optimized differential imaging to introduce a diversity in the signals of the sought-for objects, and (iii) dedicated processing algorithms to further eliminate the residual stellar leakages. Substanti… ▽ More Exoplanet detections and characterizations via direct imaging require high contrast and high angular resolution. These requirements typically require (i) cutting-edge instrumental facilities, (ii) optimized differential imaging to introduce a diversity in the signals of the sought-for objects, and (iii) dedicated processing algorithms to further eliminate the residual stellar leakages. Substantial efforts have been undertaken on the design of more efficient post-processing algorithms but their performance remains upper-bounded at shorter angular separations due to the the lack of diversity induced by the processing of each epoch of observations individually. We propose a new algorithm that is able to combine several observations of the same star by accounting for the Keplerian orbital motion across epochs of the sought-for sources in order to constructively co-add their weak signals. The proposed algorithm, PACOME, integrates an exploration of the plausible orbits within a statistical detection and estimation formalism. It is extended to a multi-epoch combination of the maximum likelihood framework of PACO, which is a mono-epoch post-processing algorithm. We derive a reliable multi-epoch detection criterion, interpretable both in terms of probability of detection and of false alarm. We tested the proposed algorithm on several datasets obtained from the VLT/SPHERE instrument with IRDIS and IFS. By resorting to injections of synthetic exoplanets, we show that PACOME is able to detect sources remaining undetectable in mono-epoch frameworks. The gain in detection sensitivity scales as high as the square root of the number of epochs. We also applied PACOME on a set of observations from the HR 8799 star hosting four known exoplanets, which are detected with very high signal-to-noise ratios. In addition, its implementation is efficient, fast, and fully automatized. △ Less

Submitted 15 September, 2023; originally announced September 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 679, A38 (2023)

Analytical modelling of adaptive optics systems: Role of the influence function

Authors: Anthony Berdeu, Michel Tallon, Éric Thiébaut, Maud Langlois

Abstract: Context. Adaptive optics (AO) is now a tool commonly deployed in astronomy. The real time correction of the atmospheric turbulence that AO enables allows telescopes to perform close to the diffraction limit at the core of their point spread function (PSF). Among other factors, AO-corrected PSFs depend on the ability of the wavefront corrector (WFC), generally a deformable mirror, to fit the incide… ▽ More Context. Adaptive optics (AO) is now a tool commonly deployed in astronomy. The real time correction of the atmospheric turbulence that AO enables allows telescopes to perform close to the diffraction limit at the core of their point spread function (PSF). Among other factors, AO-corrected PSFs depend on the ability of the wavefront corrector (WFC), generally a deformable mirror, to fit the incident wavefront corrugations. Aims. In this work, we focus on this error introduced by the WFC, the so-called fitting error. To date, analytical models only depend on the WFC cut-off frequency, and Monte Carlo simulations are the only solution for studying the impact of the WFC influence function shape on the AO-corrected PSF. We aim to develop an analytical model accounting for the influence function shape. Methods. We first obtain a general analytical model of the fitting error structure function. With additional hypotheses, we then derive an analytical model of the AO-corrected power spectral density. These two analytical solutions are compared with Monte Carlo simulations on different ideal profiles (piston, pyramid, Gaussian) as well as with real hardware (DM192 from ALPAO). Results. Our analytical predictions show a very good agreement with the Monte Carlo simulations. We show that in the image plane, the depth of the correction as well as the transition profile between the AO-corrected area and the remaining turbulent halo depend on the influence functions of the WFC. We also show that the generally assumed hypothesis of stationarity of the AO correction is actually not met. Conclusions. As the fitting error is the intrinsic optimal limit of an AO system, our analytical model allows for the assessment of the theoretical limits of extreme AO systems limited by the WFC in high-contrast imaging through a context where other errors become comparable. △ Less

Submitted 19 June, 2023; originally announced June 2023.

Journal ref: Astronomy and Astrophysics - A&A, 2023, 674, pp.A112

Characterization of stellar companion from high-contrast long-slit spectroscopy data: The EXtraction Of SPEctrum of COmpanion (EXOSPECO) algorithm

Authors: Samuel Thé, Éric Thiébaut, Loïc Denis, Thibault Wanner, Rémi Thiébaut, Maud Langlois, Ferréol Soulez

Abstract: High-contrast long-slit spectrographs can be used to characterize exoplanets. High-contrast long-slit spectroscopic data are however corrupted by stellar leakages which largely dominate other signals and make the process of extracting the companion spectrum very challenging. This paper presents a complete method to calibrate the spectrograph and extract the signal of interest. The proposed metho… ▽ More High-contrast long-slit spectrographs can be used to characterize exoplanets. High-contrast long-slit spectroscopic data are however corrupted by stellar leakages which largely dominate other signals and make the process of extracting the companion spectrum very challenging. This paper presents a complete method to calibrate the spectrograph and extract the signal of interest. The proposed method is based on a flexible direct model of the high-contrast long-slit spectroscopic data. This model explicitly accounts for the instrumental response and for the contributions of both the star and the companion. The contributions of these two components and the calibration parameters are jointly estimated by solving a regularized inverse problem. This problem having no closed-form solution, we propose an alternating minimization strategy to effectively find the solution. We have tested our method on empirical long-slit spectroscopic data and by injecting synthetic companion signals in these data. The proposed initialization and the alternating strategy effectively avoid the self-subtraction bias, even for companions observed very close to the coronagraphic mask. Careful modeling and calibration of the angular and spectral dispersion laws of the instrument clearly reduce the contamination by the stellar leakages. In practice, the outputs of the method are mostly driven by a single hyper-parameter which tunes the level of regularization of the companion SED. △ Less

Submitted 6 June, 2023; originally announced June 2023.

Comments: Paper under review by Astronomy & Astrophysics

Journal ref: A&A 678, A77 (2023)

Inverse problem approach in Extreme Adaptive Optics: analytical model of the fitting error and lowering of the aliasing

Authors: Anthony Berdeu, Michel Tallon, Éric Thiébaut, Mary Angelie Alagao, Sitthichat Sukpholtham, Maud Langlois, Adithep Kawinkij, Puttiwat Kongkaew

Abstract: We present the results obtained with an end-to-end simulator of an Extreme Adaptive Optics (XAO) system control loop. It is used to predict its on-sky performances and to optimise the AO loop algorithms. It was first used to validate a novel analytical model of the fitting error, a limit due to the Deformable Mirror (DM) shape. Standard analytical models assume a sharp correction under the DM cuto… ▽ More We present the results obtained with an end-to-end simulator of an Extreme Adaptive Optics (XAO) system control loop. It is used to predict its on-sky performances and to optimise the AO loop algorithms. It was first used to validate a novel analytical model of the fitting error, a limit due to the Deformable Mirror (DM) shape. Standard analytical models assume a sharp correction under the DM cutoff frequency, disregarding the transition between the AO corrected and turbulence dominated domains. Our model account for the influence function shape in this smooth transition. Then, it is well-known that Shack-Hartmann wavefront sensors (SH-WFS) have a limited spatial bandwidth, the high frequencies of the wavefront being seen as low frequencies. We show that this aliasing error can be partially compensated (both in terms of Strehl ratio and contrast) by adding priors on the turbulence statistics in the framework of an inverse problem approach. This represents an alternative to the standard additional optical filter used in XAO systems. In parallel to this numerical work, a bench was aligned to experimentally test the AO system and these new algorithms comprising a DM192 ALPAO deformable mirror and a 15x15 SH-WFS. We present the predicted performances of the AO loop based on end-to-end simulations. △ Less

Submitted 9 January, 2023; originally announced January 2023.

Journal ref: Adaptive Optics Systems VIII, Jul 2022, Montr{é}al, Canada. pp.121850K

Adaptive Optics system of the Evanescent Wave Coronagraph (EvWaCo): optimised phase plate and DM characterisation

Authors: Anthony Berdeu, Sitthichat Sukpholtham, Puttiwat Kongkaew, Adithep Kawinkij, Matthew Ridsdill-Smith, Michel Tallon, Éric Thiébaut, Maud Langlois, Mary Angelie Alagao

Abstract: The Evanescent Wave Coronagraph (EvWaCo) is an achromatic coronagraph mask with adjustable size over the spectral domain [600nm, 900nm] that will be installed at the Thai National Observatory. We present in this work the development of a bench to characterise its Extreme Adaptive Optics system (XAO) comprising a DM192 ALPAO deformable mirror (DM) and a 15x15 Shack-Hartmann wavefront sensor (SH-WFS… ▽ More The Evanescent Wave Coronagraph (EvWaCo) is an achromatic coronagraph mask with adjustable size over the spectral domain [600nm, 900nm] that will be installed at the Thai National Observatory. We present in this work the development of a bench to characterise its Extreme Adaptive Optics system (XAO) comprising a DM192 ALPAO deformable mirror (DM) and a 15x15 Shack-Hartmann wavefront sensor (SH-WFS). In this bench, the turbulence is simulated using a rotating phase plate in a pupil plane. In general, such components are designed using a randomly generated phase screen. Such single realisation does not necessarily provide the wanted structure function. We present a solution to design the printed pattern to ensure that the beam sees a strict and controlled Kolmogorov statistics with the correct 2D structure function. This is essential to control the experimental conditions in order to compare the bench results with the numerical simulations and predictions. This bench is further used to deeply characterise the full 27 mm pupil of the ALPAO DM using a 54x54 ALPAO SH-WFS. We measure the average shape of its influence functions as well as the influence function of each single actuator to study their dispersion. We study the linearity of the actuator amplitude with the command as well as the linearity of the influence function profile. We also study the actuator offsets as well as the membrane shape at 0-command. This knowledge is critical to get a forward model of the DM for the XAO control loop. △ Less

Submitted 28 December, 2022; originally announced December 2022.

IMAGE-OI: an OIFITS extension and its application in OImaging to compare image reconstruction algorithms

Authors: Ferréol Soulez, Laurent Bourgès, Antoine Kaszczyc, Guillaume Mella, Martin Pratoussy, Gilles Duvert, Jacques Kluska, Eric Thiébaut, John Young

Abstract: In interferometry, the quality of the reconstructed image depends on the algorithm used and its parameters, and users often need to compare the results of several algorithms to disentangle artifacts from actual features of the astrophysical object. Such comparisons can rapidly become cumbersome, as these software packages are very different. OImaging is a graphical interface intended to be a commo… ▽ More In interferometry, the quality of the reconstructed image depends on the algorithm used and its parameters, and users often need to compare the results of several algorithms to disentangle artifacts from actual features of the astrophysical object. Such comparisons can rapidly become cumbersome, as these software packages are very different. OImaging is a graphical interface intended to be a common frontend to image reconstruction software packages. With OImaging, the user can now perform multiple reconstructions within a single interface. From a given dataset, OImaging allows benchmarking of different image reconstruction algorithms and assessment of the reliability of the image reconstruction process. To that end, OImaging uses the IMAGE-OI OIFITS extension proposed to standardize communication with image reconstruction algorithms. △ Less

Submitted 16 September, 2022; originally announced September 2022.

Journal ref: Proc. SPIE 12183, Optical and Infrared Interferometry and Imaging VIII, 121831H (26 August 2022);

Upgrading the high contrast imaging facility SPHERE: science drivers and instrument choices

Authors: A. Boccaletti, G. Chauvin, F. Wildi, J. Milli, E. Stadler, E. Diolaiti, R. Gratton, F. Vidal, M. Loupias, M. Langlois, F. Cantalloube, M. N'Diaye, D. Gratadour, F. Ferreira, M. Tallon, J. Mazoyer, D. Segransan, D. Mouillet, J. -L. Beuzit, M. Bonnefoy, R. Galicher, A. Vigan, I. Snellen, M. Feldt, S. Desidera , et al. (49 additional authors not shown)

Abstract: SPHERE+ is a proposed upgrade of the SPHERE instrument at the VLT, which is intended to boost the current performances of detection and characterization for exoplanets and disks. SPHERE+ will also serve as a demonstrator for the future planet finder (PCS) of the European ELT. The main science drivers for SPHERE+ are 1/ to access the bulk of the young giant planet population down to the snow line (… ▽ More SPHERE+ is a proposed upgrade of the SPHERE instrument at the VLT, which is intended to boost the current performances of detection and characterization for exoplanets and disks. SPHERE+ will also serve as a demonstrator for the future planet finder (PCS) of the European ELT. The main science drivers for SPHERE+ are 1/ to access the bulk of the young giant planet population down to the snow line ($3-10$ au), to bridge the gap with complementary techniques (radial velocity, astrometry); 2/ to observe fainter and redder targets in the youngest ($1-10$\,Myr) associations compared to those observed with SPHERE to directly study the formation of giant planets in their birth environment; 3/ to improve the level of characterization of exoplanetary atmospheres by increasing the spectral resolution in order to break degeneracies in giant planet atmosphere models. Achieving these objectives requires to increase the bandwidth of the xAO system (from $\sim$1 to 3\,kHz) as well as the sensitivity in the infrared (2 to 3\,mag). These features will be brought by a second stage AO system optimized in the infrared with a pyramid wavefront sensor. As a new science instrument, a medium resolution integral field spectrograph will provide a spectral resolution from 1000 to 5000 in the J and H bands. This paper gives an overview of the science drivers, requirements and key instrumental trade-off that were done for SPHERE+ to reach the final selected baseline concept. △ Less

Submitted 5 September, 2022; originally announced September 2022.

Comments: To appear in the Proceedings of the SPIE Astronomical Telescopes + Instrumentation (2022), 13 pages, 6 figure

REXPACO ASDI: Joint unmixing and deconvolution of the circumstellar environment by angular and spectral differential imaging

Authors: Olivier Flasseur, Loïc Denis, Éric Thiébaut, Maud Langlois

Abstract: Angular and spectral differential imaging is an observational technique of choice to investigate the immediate vicinity of stars. The relative angular motion and spectral scaling between on-axis and off-axis sources are exploited by post-processing techniques to separate two components: the residual star light and the light coming from surrounding objects such as circumstellar disks or point-like… ▽ More Angular and spectral differential imaging is an observational technique of choice to investigate the immediate vicinity of stars. The relative angular motion and spectral scaling between on-axis and off-axis sources are exploited by post-processing techniques to separate two components: the residual star light and the light coming from surrounding objects such as circumstellar disks or point-like objects. This paper introduces a new algorithm to jointly unmix these two components and deconvolve disk images. The proposed algorithm is based on a statistical modeling of the residual star light, accounting for its spatial and spectral correlations. While critical, these correlations are not modeled and compensated for by existing reconstruction algorithms. We leverage dedicated shrinkage techniques to estimate the large number of parameters of our model of correlations in a data-driven fashion. We show that the resulting separable model of the spatial and spectral covariances captures very accurately the star light, thus allowing its efficient suppression. We apply our method on several datasets from the VLT/SPHERE instrument and compare the performances against standard algorithms of the field (cADI, PCA). We show that accounting for the multiple correlations of the data significantly enhances the reconstruction quality, leading to better preservation of both the disk morphology and photometry. Thanks to its unique joint spectral modeling of the data, the proposed algorithm can reconstruct disks having a circular symmetry (e.g., rings, spirals) at intensities one million times fainter than the star, without resorting to additional reference datasets free from the off-axis objects. △ Less

Submitted 14 May, 2024; v1 submitted 26 September, 2021; originally announced September 2021.

Comments: Submitted to MNRAS

RHAPSODIE: Reconstruction of High-contrAst Polarized SOurces and Deconvolution for cIrcumstellar Environments

Authors: Laurence Denneulin, Maud Langlois, Éric Thiébaut, Nelly Pustelnik

Abstract: Polarimetric imaging is one of the most effective techniques for high-contrast imaging and characterization of circumstellar environments. These environments can be characterized through direct-imaging polarimetry at near-infrared wavelengths. The SPHERE/IRDIS instrument installed on the Very Large Telescope in its dual-beam polarimetric imaging (DPI) mode, offers the capability to acquire polarim… ▽ More Polarimetric imaging is one of the most effective techniques for high-contrast imaging and characterization of circumstellar environments. These environments can be characterized through direct-imaging polarimetry at near-infrared wavelengths. The SPHERE/IRDIS instrument installed on the Very Large Telescope in its dual-beam polarimetric imaging (DPI) mode, offers the capability to acquire polarimetric images at high contrast and high angular resolution. However dedicated image processing is needed to get rid of the contamination by the stellar light, of instrumental polarization effects, and of the blurring by the instrumental point spread function. We aim to reconstruct and deconvolve the near-infrared polarization signal from circumstellar environments. We use observations of these environments obtained with the high-contrast imaging infrared polarimeter SPHERE-IRDIS at the VLT. We developed a new method to extract the polarimetric signal using an inverse approach method that benefits from the added knowledge of the detected signal formation process. The method includes weighted data fidelity term, smooth penalization, and takes into account instrumental polarization. The method enables to accurately measure the polarized intensity and angle of linear polarization of circumstellar disks by taking into account the noise statistics and the convolution of the observed objects by the instrumental point spread function. It has the capability to use incomplete polarimetry cycles which enhance the sensitivity of the observations. The method improves the overall performances in particular for low SNR/small polarized flux compared to standard methods. △ Less

Submitted 12 May, 2021; originally announced May 2021.

Journal ref: A&A 653, A138 (2021)

REXPACO: an algorithm for high contrast reconstruction of the circumstellar environment by angular differential imaging

Authors: Olivier Flasseur, Samuel Thé, Loïc Denis, Éric Thiébaut, Maud Langlois

Abstract: Aims. The purpose of this paper is to describe a new post-processing algorithm dedicated to the reconstruction of the spatial distribution of light received from off-axis sources, in particular from circumstellar disks. Methods. Built on the recent PACO algorithm dedicated to the detection of point-like sources, the proposed method is based on the local learning of patch covariances capturing th… ▽ More Aims. The purpose of this paper is to describe a new post-processing algorithm dedicated to the reconstruction of the spatial distribution of light received from off-axis sources, in particular from circumstellar disks. Methods. Built on the recent PACO algorithm dedicated to the detection of point-like sources, the proposed method is based on the local learning of patch covariances capturing the spatial fluctuations of the stellar leakages. From this statistical modeling, we develop a regularized image reconstruction algorithm (REXPACO) following an inverse problem approach based on a forward image formation model of the off-axis sources in the ADI sequences. Results. Injections of fake circumstellar disks in ADI sequences from the VLT/SPHERE-IRDIS instrument show that both the morphology and the photometry of the disks are better preserved by REXPACO compared to standard postprocessing methods like cADI. In particular, the modeling of the spatial covariances proves usefull in reducing typical ADI artifacts and in better disentangling the signal of these sources from the residual stellar contamination. The application to stars hosting circumstellar disks with various morphologies confirms the ability of REXPACO to produce images of the light distribution with reduced artifacts. Finally, we show how REXPACO can be combined with PACO to disentangle the signal of circumstellar disks from the signal of candidate point-like sources. Conclusions. REXPACO is a novel post-processing algorithm producing numerically deblurred images of the circumstellar environment. It exploits the spatial covariances of the stellar leakages and of the noise to efficiently eliminate this nuisance term. △ Less

Submitted 13 May, 2021; v1 submitted 19 April, 2021; originally announced April 2021.

Journal ref: A&A 651, A62 (2021)

Deep Contrast and Companion Detection Using the EvWaCo Testbed Equipped with an Achromatic Focal Plane Mask and an Adjustable Inner Working Angle

Authors: M. A. Alagao, A. Kawinkij, C. Buisset, A. Prasit, T. Lépine, Y. Rabbia, A. Berdeu, É. Thiébaut, M. Langlois, M. Tallon, S. Awiphan, E. Semenko, P. Rittipruk, D. Mkrtichian, A. Leckngam, G. Thuammasorn, P. Kaewsamoeta, A. Inpan, T. Kuha, A. Laoyang, W. Somboonchai, S. Kanthum, S. Poshyachinda, B. Soonthornthum

Abstract: The evanescent wave coronagraph uses the principle of frustrated total internal reflection (FTIR) to suppress the light coming from the star and study its close environment. Its focal plane mask is composed of a lens and a prism placed in contact with each other to produce the coronagraphic effect. In this paper, we present the experimental results obtained using an upgraded focal plane mask of th… ▽ More The evanescent wave coronagraph uses the principle of frustrated total internal reflection (FTIR) to suppress the light coming from the star and study its close environment. Its focal plane mask is composed of a lens and a prism placed in contact with each other to produce the coronagraphic effect. In this paper, we present the experimental results obtained using an upgraded focal plane mask of the Evanescent Wave Coronagraph (EvWaCo). These experimental results are also compared to the theoretical performance of the coronagraph obtained through simulations. Experimentally, we reach a raw contrast equal to a few $10^{-4}$ at a distance equal to 3 $λ/D$ over the full I-band ($λ_c =$ $800$ $nm$, $Δλ/λ \approx 20\%$) and equal to 4 $λ/D$ over the full R-band ($λ_c =$ $650$ $nm$, $Δλ/λ \approx 23\%$) in unpolarized light. However, our simulations show a raw contrast close to $10^{-4}$ over the full I-band and R-band at the same distance, thus, confirming the theoretical achromatic advantage of the coronagraph. We also verify the stability of the mask through a series of contrast measurements over a period of 8 months. Furthermore, we measure the sensitivity of the coronagraph to the lateral and longitudinal misalignment of the focal plane mask, and to the lateral misalignment of the Lyot stop. △ Less

Submitted 19 March, 2021; v1 submitted 18 March, 2021; originally announced March 2021.

Comments: 27 pages, 19 figures, accepted by AJ

The SPHERE infrared survey for exoplanets (SHINE) -- II. Observations, Data reduction and analysis Detection performances and early-results

Authors: M. Langlois, R. Gratton, A. -M. Lagrange, P. Delorme, A. Boccaletti, M. Bonnefoy, A. -L. Maire, D. Mesa, G. Chauvin, S. Desidera, A. Vigan, A. Cheetham, J. Hagelberg, M. Feldt, M. Meyer, P. Rubini, H. Le Coroller, F. Cantalloube, B. Biller, M. Bonavita, T. Bhowmik, W. Brandner, S. Daemgen, V. D'Orazi, O. Flasseur , et al. (96 additional authors not shown)

Abstract: Over the past decades, direct imaging has confirmed the existence of substellar companions (exoplanets or brown dwarfs) on wide orbits (>10 au) from their host stars. To understand their formation and evolution mechanisms, we have initiated in 2015 the SPHERE infrared survey for exoplanets (SHINE), a systematic direct imaging survey of young, nearby stars to explore their demographics.} {We aim to… ▽ More Over the past decades, direct imaging has confirmed the existence of substellar companions (exoplanets or brown dwarfs) on wide orbits (>10 au) from their host stars. To understand their formation and evolution mechanisms, we have initiated in 2015 the SPHERE infrared survey for exoplanets (SHINE), a systematic direct imaging survey of young, nearby stars to explore their demographics.} {We aim to detect and characterize the population of giant planets and brown dwarfs beyond the snow line around young, nearby stars. Combined with the survey completeness, our observations offer the opportunity to constrain the statistical properties (occurrence, mass and orbital distributions, dependency on the stellar mass) of these young giant planets.} {In this study, we present the observing and data analysis strategy, the ranking process of the detected candidates, and the survey performances for a subsample of 150 stars, which are representative of the full SHINE sample. The observations were conducted in an homogeneous way from February 2015 to February 2017 with the dedicated ground-based VLT/SPHERE instrument equipped with the IFS integral field spectrograph and the IRDIS dual-band imager covering a spectral range between 0.9 and 2.3 $μ$m. We used coronographic, angular and spectral differential imaging techniques to reach the best detection performances for this study down to the planetary mass regime.} △ Less

Submitted 5 March, 2021; originally announced March 2021.

Journal ref: A&A 651, A71 (2021)

The wavefront sensing making-of for THEMIS solar telescope

Authors: Michel Tallon, Éric Thiébaut, Maud Langlois, Bernard Gelly, Richard Douet, Clémentine Béchet, Loïc Denis, Gil Moretto

Abstract: An adaptive optics system with a single deformable mirror is being implemented on the THEMIS 90cm solar telescope. This system is designed to operate in the visible and is required to be as robust as possible in order to deliver the best possible correction in any atmospheric conditions, even if wavefronts are sensed on some low-contrast solar granulation. In extreme conditions, the images given b… ▽ More An adaptive optics system with a single deformable mirror is being implemented on the THEMIS 90cm solar telescope. This system is designed to operate in the visible and is required to be as robust as possible in order to deliver the best possible correction in any atmospheric conditions, even if wavefronts are sensed on some low-contrast solar granulation. In extreme conditions, the images given by the subapertures of the Shack-Hartmann wavefront sensor get randomly blurred in space, in the set of subapertures, and the distribution of blurred images is rapidly changing in time, some of them possibly fading away. The algorithms we have developed for such harsh conditions rely on inverse problem approach. As an example, with the gradients of the wavefronts, the wavefront sensor also estimates their errors, including their covariance. This information allows the control loop to promptly optimize itself to the fast varying conditions, both in space (wavefront reconstruction) and in time. A major constraint is to fit the calculations in a low-cost multi-core CPU. An overview of the algorithms in charge of implementing this strategy is presented, focusing on wavefront sensing. △ Less

Submitted 8 January, 2021; originally announced January 2021.

Journal ref: Adaptive Optics for Extremely Large Telescopes 6, Jun 2019, Qu{é}bec CIty, Canada

Evidence for localized onset of episodic mass loss in Mira

Authors: G. Perrin, S. T. Ridgway, S. Lacour, X. Haubois, E. Thiebaut, J. P. Berger, M. G. Lacasse, R. Millan-Gabet, J. D. Monnier, E. Pedretti, S. Ragland, W. Traub

Abstract: We report Multi-telescope interferometric measurements taken with the Interferometric Optical Telescope Array (IOTA) to provide imagery of the LPV Mira in the H-band. This wavelength region is well suited to studying mass loss given the low continuum opacity, which allows for emission to be observed over a very long path in the stellar atmosphere and envelope. The observed visibilities are consist… ▽ More We report Multi-telescope interferometric measurements taken with the Interferometric Optical Telescope Array (IOTA) to provide imagery of the LPV Mira in the H-band. This wavelength region is well suited to studying mass loss given the low continuum opacity, which allows for emission to be observed over a very long path in the stellar atmosphere and envelope. The observed visibilities are consistent with a simple core-halo model to represent the central object and the extended molecular layers but, in addition, they demonstrate a substantial asymmetry. An analysis with image reconstruction software shows that the asymmetry is consistent with a localized absorbing patch. The observed opacity is tentatively associated with small dust grains, which will grow substantially during a multi-year ejection process. Spatial information along with a deduced dust content of the cloud, known mass loss rates, and ejection velocities provide evidence for the pulsational pumping of the extended molecular layers. The cloud may be understood as a spatially local zone of enhanced dust formation, very near to the pulsating halo. The observed mass loss could be provided by several such active regions around the star. This result provides an additional clue for better understanding the clumpiness of dust production in the atmosphere of AGB stars. It is compatible with scenarios where the combination of pulsation and convection play a key role in the process of mass loss. △ Less

Submitted 22 August, 2020; originally announced August 2020.

Comments: Accepted by A&A, 12 pages, 9 figures

Journal ref: A&A 642, A82 (2020)

A family portrait of disk inner rims around Herbig Ae/Be stars: Hunting for warps, rings, self shadowing, and misalignments in the inner astronomical units

Authors: J. Kluska, J. -P. Berger, F. Malbet, B. Lazareff, M. Benisty, J. -B. Le Bouquin, O. Absil, F. Baron, A. Delboulbé, G. Duvert, A. Isella, L. Jocou, A. Juhasz, S. Kraus, R. Lachaume, F. Ménard, R. Millan-Gabet, J. D. Monnier, T. Moulin, K. Perraut, C. Pinte, S. Rochat, F. Soulez, M. Tallon, W. -F. Thi , et al. (3 additional authors not shown)

Abstract: 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 inter… ▽ More 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. △ Less

Submitted 3 April, 2020; originally announced April 2020.

Comments: Accepted in A&A, 24pages, 22 figures

Journal ref: A&A 636, A116 (2020)

SPHERE+: Imaging young Jupiters down to the snowline

Authors: A. Boccaletti, G. Chauvin, D. Mouillet, O. Absil, F. Allard, S. Antoniucci, J. -C. Augereau, P. Barge, A. Baruffolo, J. -L. Baudino, P. Baudoz, M. Beaulieu, M. Benisty, J. -L. Beuzit, A. Bianco, B. Biller, B. Bonavita, M. Bonnefoy, S. Bos, J. -C. Bouret, W. Brandner, N. Buchschache, B. Carry, F. Cantalloube, E. Cascone , et al. (108 additional authors not shown)

Abstract: 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 S… ▽ More 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. △ Less

Submitted 13 March, 2020; v1 submitted 12 March, 2020; originally announced March 2020.

Comments: White paper submitted to ESO on Feb. 20th, 2020

Why Chromatic Imaging Matters

Authors: Joel Sanchez-Bermudez, Florentin Millour, Fabien Baron, Roy van Boekel, Laurent Bourgès, Gilles Duvert, Paulo J. V. Garcia, Nuno Gomes, Karl-Heinz Hofmann, Thomas Henning, Jacob W. Isbell, Bruno Lopez, Alexis Matter, J-Uwe Pott, Dieter Schertl, Eric Thiébaut, Gerd Weigelt, John Young

Abstract: During the last two decades, the first generation of beam combiners at the Very Large Telescope Interferometer has proved the importance of optical interferometry for high-angular resolution astrophysical studies in the near- and mid-infrared. With the advent of 4-beam combiners at the VLTI, the u-v coverage per pointing increases significantly, providing an opportunity to use reconstructed images… ▽ More During the last two decades, the first generation of beam combiners at the Very Large Telescope Interferometer has proved the importance of optical interferometry for high-angular resolution astrophysical studies in the near- and mid-infrared. With the advent of 4-beam combiners at the VLTI, the u-v coverage per pointing increases significantly, providing an opportunity to use reconstructed images as powerful scientific tools. Therefore, interferometric imaging is already a key feature of the new generation of VLTI instruments, as well as for other interferometric facilities like CHARA and JWST. It is thus imperative to account for the current image reconstruction capabilities and their expected evolutions in the coming years. Here, we present a general overview of the current situation of optical interferometric image reconstruction with a focus on new wavelength-dependent information, highlighting its main advantages and limitations. As an Appendix we include several cookbooks describing the usage and installation of several state-of-the art image reconstruction packages. To illustrate the current capabilities of the software available to the community, we recovered chromatic images, from simulated MATISSE data, using the MCMC software SQUEEZE. With these images, we aim at showing the importance of selecting good regularization functions and their impact on the reconstruction. △ Less

Submitted 14 December, 2018; originally announced December 2018.

Comments: Accepted for publication in Experimental Astronomy as part of the topical collection: Future of Optical-infrared Interferometry in Europe

Proximity Operators for Phase Retrieval

Authors: Ferréol Soulez, Éric Thiébaut, Antony Schutz, André Ferrari, Frédéric Courbin, Michael Unser

Abstract: 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 a… ▽ More 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. △ Less

Submitted 27 October, 2017; originally announced October 2017.

Journal ref: Applied Optics Vol. 55, Issue 26, pp. 7412-7421 (2016)

Principles of image reconstruction in optical interferometry: tutorial

Authors: Éric Thiébaut, John Young

Abstract: This paper provides a general introduction to the problem of image reconstruction from interferometric data. A simple model of the interferometric observables is given and the issues arising from sparse Fourier data are discussed. The effects of various regularizations are described. In the proposed general framework, most existing algorithms can be understood. For an astronomer, such an understan… ▽ More This paper provides a general introduction to the problem of image reconstruction from interferometric data. A simple model of the interferometric observables is given and the issues arising from sparse Fourier data are discussed. The effects of various regularizations are described. In the proposed general framework, most existing algorithms can be understood. For an astronomer, such an understanding is crucial not only for selecting and using an algorithm but also to ensure correct interpretation of the resulting image. △ Less

Submitted 28 August, 2017; originally announced August 2017.

Comments: 20 pages, 5 figures

Journal ref: Journal of the Optical Society of America A, vol. 34, pp. 904-923 (2017)

PeX 1. Multi-spectral expansion of residual speckles for planet detection

Authors: Nicholas Devaney, Éric Thiébaut

Abstract: The detection of exoplanets in coronographic images is severely limited by residual starlight speckles. Dedicated post-processing can drastically reduce this "stellar leakage" and thereby increase the faintness of detectable exoplanets. Based on a multi-spectral series expansion of the diffraction pattern, we derive a multi-mode model of the residuals which can be exploited to estimate and thus re… ▽ More The detection of exoplanets in coronographic images is severely limited by residual starlight speckles. Dedicated post-processing can drastically reduce this "stellar leakage" and thereby increase the faintness of detectable exoplanets. Based on a multi-spectral series expansion of the diffraction pattern, we derive a multi-mode model of the residuals which can be exploited to estimate and thus remove the residual speckles in multi-spectral coronographic images. Compared to other multi-spectral processing methods, our model is physically grounded and is suitable for use in an (optimal) inverse approach. We demonstrate the ability of our model to correctly estimate the speckles in simulated data and demonstrate that very high contrasts can be achieved. We further apply our method to removing speckles from a real data cube obtained with the SPHERE IFS instrument. △ Less

Submitted 29 August, 2017; v1 submitted 28 August, 2017; originally announced August 2017.

Comments: accepted for publication in MNRAS on 25th of August 2017, 17 pages, 15 figures

Science cases for a visible interferometer

Authors: Philippe Stee, France Allard, Myriam Benisty, Lionel Bigot, Nicolas Blind, Henri Boffin, Marcelo Borges Fernandes, Alex Carciofi, Andrea Chiavassa, Orlagh Creevey, Pierre Cruzalebes, Willem-Jan de Wit, Armando Domiciano de Souza, Martin Elvis, Nicolas Fabas, Daniel Faes, Alexandre Gallenne, Carlos Guerrero Pena, Michel Hillen, Sebastian Hoenig, Michael Ireland, Pierre Kervella, Makoto Kishimoto, Nadia Kostogryz, Stefan Kraus , et al. (32 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 21 March, 2017; v1 submitted 7 March, 2017; originally announced March 2017.

Comments: White Paper prospective. This book is dedicated to the memory of our colleague Olivier Chesneau who passed away at the age of 41. v2 includes some corrections to text

Structure of Herbig AeBe disks at the milliarcsecond scale A statistical survey in the H band using PIONIER-VLTI

Authors: B. Lazareff, J. -P. Berger, J. Kluska, J. -B. Le Bouquin, M. Benisty, F. Malbet, C. Koen, C. Pinte, W. -F. Thi, O. Absil, F. Baron, A. Delboulbé, G. Duvert, A. Isella, L. Jocou, A. Juhasz, S. Kraus, R. Lachaume, F. Ménard, R. Millan-Gabet, J. D. Monnier, T. Moulin, K. Perraut, S. Rochat, F. Soulez , et al. (4 additional authors not shown)

Abstract: 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 interf… ▽ More 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. △ Less

Submitted 25 November, 2016; originally announced November 2016.

Comments: 41 pages. 23 figures. Four appendices

Journal ref: A&A 599, A85 (2017)

Assessing the quality of restored images in optical long-baseline interferometry

Authors: Nuno Gomes, Paulo J. V. Garcia, Éric Thiébaut

Abstract: Assessing the quality of aperture synthesis maps is relevant for benchmarking image reconstruction algorithms, for the scientific exploitation of data from optical long-baseline interferometers, and for the design/upgrade of new/existing interferometric imaging facilities. Although metrics have been proposed in these contexts, no systematic study has been conducted on the selection of a robust met… ▽ More Assessing the quality of aperture synthesis maps is relevant for benchmarking image reconstruction algorithms, for the scientific exploitation of data from optical long-baseline interferometers, and for the design/upgrade of new/existing interferometric imaging facilities. Although metrics have been proposed in these contexts, no systematic study has been conducted on the selection of a robust metric for quality assessment. This article addresses the question: what is the best metric to assess the quality of a reconstructed image? It starts by considering several metrics, and selecting a few based on general properties. Then, a variety of image reconstruction cases is considered. The observational scenarios are phase closure and phase referencing at the Very Large Telescope Interferometer (VLTI), for a combination of two, three, four and six telescopes. End-to-end image reconstruction is accomplished with the MiRA software, and several merit functions are put to test. It is found that convolution by an effective point spread function is required for proper image quality assessment. The effective angular resolution of the images is superior to naive expectation based on the maximum frequency sampled by the array. This is due to the prior information used in the aperture synthesis algorithm and to the nature of the objects considered. The l1 norm is the most robust of all considered metrics, because being linear it is less sensitive to image smoothing by high regularisation levels. For the cases considered, this metric allows the implementation of automatic quality assessment of reconstructed images, with a performance similar to human selection. △ Less

Submitted 20 February, 2017; v1 submitted 18 November, 2016; originally announced November 2016.

Comments: Key words: instrumentation: high angular resolution - instrumentation: interferosmeters - methods: data analysis - techniques: high angular resolution - techniques: image processing - techniques: interferometric. 18 pages, 8 figures, 4 tables

A disk asymmetry in motion around the B[e] star MWC158

Authors: J. Kluska, M. Benisty, F. Soulez, J. -P. Berger, J. -B. Le Bouquin, F. Malbet, B. Lazareff, E. Thiébaut

Abstract: 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 th… ▽ More 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. △ Less

Submitted 17 May, 2016; originally announced May 2016.

Comments: Accepted for publication in A&A

Journal ref: A&A 591, A82 (2016)

Development of an ELT XAO testbed using a Mach-Zehnder wavefront sensor: calibration of the deformable mirror

Authors: Christian Delacroix, Maud Langlois, Magali Loupias, Eric Thiébaut, Louisa Adjali, Jonathan Leger, Michel Tallon

Abstract: (abridged) Extreme adaptive optics (XAO) encounters severe difficulties to cope with the high speed (>1kHz), high accuracy and high order requirements for future extremely large telescopes. An innovative high order adaptive optics system using a self-referenced Mach-Zehnder wavefront sensor (MZWFS) allows counteracting these limitations. This sensor estimates very accurately the wavefront phase at… ▽ More (abridged) Extreme adaptive optics (XAO) encounters severe difficulties to cope with the high speed (>1kHz), high accuracy and high order requirements for future extremely large telescopes. An innovative high order adaptive optics system using a self-referenced Mach-Zehnder wavefront sensor (MZWFS) allows counteracting these limitations. This sensor estimates very accurately the wavefront phase at small spatial scale by measuring intensity differences between two outputs, with a $λ/4$ path length difference between its two legs, but is limited in dynamic range due to phase ambiguity. During the past few years, such an XAO system has been studied by our team in the framework of 8-meter class telescopes. In this work, we report on our latest results with the XAO testbed recently installed in our lab, and dedicated to high contrast imaging with 30m-class telescopes (such as the E-ELT or the TMT). After reminding the principle of a MZWFS and describing the optical layout of our experiment, we will show the results of the assessment of the woofer-tweeter phase correctors, i.e., a Boston Micromachine continuous membrane deformable mirror (DM) and a Boulder Nonlinear Systems liquid crystal spatial light modulator (SLM). In particular, we will detail the calibration of the DM using Zygo interferometer metrology. Our method consists in the precise measurement of the membrane deformation while applying a constant deformation to 9 out of 140 actuators at the same time. By varying the poke voltage across the DM operating range, we propose a simple but efficient way of modeling the DM influence function using a Gaussian model. Finally, we show the DM flattening on the MZWFS allowing to compensate for low order aberrations. △ Less

Submitted 28 September, 2015; originally announced September 2015.

Comments: To appear in SPIE proceedings vol. 9617

Resolving asymmetries along the pulsation cycle of the Mira star X Hya

Authors: X. Haubois, M. Wittkowski, G. Perrin, P. Kervella, A. Mérand, E. Thiébaut, S. T. Ridgway, M. Ireland, M. Scholz

Abstract: The mass-loss process in Mira stars probably occurs in an asymmetric way where dust can form in inhomogeneous circumstellar molecular clumps. Following asymmetries along the pulsation cycle can give us clues about these mass-loss processes. We imaged the Mira star X Hya and its environnement at different epochs to follow the evolution of the morphology in the continuum and in the molecular bands.… ▽ More The mass-loss process in Mira stars probably occurs in an asymmetric way where dust can form in inhomogeneous circumstellar molecular clumps. Following asymmetries along the pulsation cycle can give us clues about these mass-loss processes. We imaged the Mira star X Hya and its environnement at different epochs to follow the evolution of the morphology in the continuum and in the molecular bands. We observed X Hya with AMBER in J-H-K at low resolution at two epochs. We modelled squared visibilities with geometrical and physical models. We also present imaging reconstruction results obtained with MiRA and based on the physical a priori images. We report on the angular scale change of X Hya between the two epochs. 1D CODEX profiles allowed us to understand and model the spectral variation of squared visibilities and constrain the stellar parameters. Reconstructed model-dependent images enabled us to reproduce closure phase signals and the azimuthal dependence of squared visibilities. They show evidence for material inhomogeneities located in the immediate environment of the star. △ Less

Submitted 13 August, 2015; originally announced August 2015.

Comments: Accepted for publication in A&A, 17 pages, 16 figures

Journal ref: A&A 582, A71 (2015)

Large Scale 3D Image Reconstruction in Optical Interferometry

Authors: Antony Schutz, André Ferrari, David Mary, Eric Thiébaut, Ferréol Soulez

Abstract: 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, specifi… ▽ More 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. △ Less

Submitted 5 March, 2015; originally announced March 2015.

Comments: EUSIPCO, Aug 2015, NICE, France

Interferometric imaging diagnostics of X Hya's circumstellar environment

Authors: X. Haubois, M. Wittkowski, G. Perrin, P. Kervella, S. T. Ridgway, E. Thiébaut

Abstract: Optical interferometry is a powerful tool to investigate the close environment of AGB stars. With a spatial resolution of a few milli-arcseconds, it is even possible to image directly the surface of angularly large objects. This is of special interest forMira stars and red supergiants for which the dust-wind is initiated from or very close to the photosphere by an interplay between pulsation and c… ▽ More Optical interferometry is a powerful tool to investigate the close environment of AGB stars. With a spatial resolution of a few milli-arcseconds, it is even possible to image directly the surface of angularly large objects. This is of special interest forMira stars and red supergiants for which the dust-wind is initiated from or very close to the photosphere by an interplay between pulsation and convection. Based on two-epoch interferometric observations of the Mira star X Hya, we present how the variation of the angular size with wavelength challenges pulsation models and how reconstructed images can reveal the evolution of the object shape and of its asymmetric structures. △ Less

Submitted 4 March, 2015; originally announced March 2015.

Comments: Proceedings of the conference "Why Galaxies care about AGB III"

PAINTER: a spatio-spectral image reconstruction algorithm for optical interferometry

Authors: Antony Schutz, André Ferrari, David Mary, Férréol Soulez, Éric Thiébaut, Martin Vannier

Abstract: 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 polychr… ▽ More 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. △ Less

Submitted 27 September, 2014; v1 submitted 29 June, 2014; originally announced July 2014.

Comments: 12 pages, 10 figures, http://www.opticsinfobase.org/submit/review/copyright_permissions.cfm

SPARCO : a semi-parametric approach for image reconstruction of chromatic objects

Authors: J. Kluska, F. Malbet, J. -P. Berger, F. Baron, B. Lazareff, J. -B. Le Bouquin, J. D. Monnier, F. Soulez, E. Thiébaut

Abstract: 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… ▽ More 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. △ Less

Submitted 13 March, 2014; originally announced March 2014.

Comments: 11 pages, 11 figures, accepted in A&A

Exploiting spatial sparsity for multi-wavelength imaging in optical interferometry

Authors: Éric Thiébaut, Ferréol Soulez, Loïc Denis

Abstract: 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. Th… ▽ More 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. △ Less

Submitted 18 December, 2012; v1 submitted 11 September, 2012; originally announced September 2012.

Comments: This version has been accepted for publication in J. Opt. Soc. Am. A

The 2012 Interferometric Imaging Beauty Contest

Authors: Fabien Baron, William D. Cotton, Peter R. Lawson, Steve T. Ridgway, Alicia Aarnio, John D. Monnier, Karl-Heinz Hofmann, Dieter Schertl, Gerd Weigelt, Eric Thiébaut, Férréol Soulez, David Mary, Florentin Millour, Martin Vannier, John Young, Nicholas M. Elias II, Henrique R. Schmitt, Sridharan Rengaswamy

Abstract: 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 wit… ▽ More 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ébaut 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. △ Less

Submitted 1 August, 2012; v1 submitted 30 July, 2012; originally announced July 2012.

Imaging the heart of astrophysical objects with optical long-baseline interferometry

Authors: J. -P. Berger, F. Malbet, F. Baron, A. Chiavassa, G. Duvert, M. Elitzur, B. Freytag, F. Gueth, S. Hönig, J. Hron, H. Jang-Condell, J. -B. Le Bouquin, J. -L Monin, J. D. Monnier, G. Perrin, B. Plez, T. Ratzka, S. Renard, S. Stefl, E. Thiébaut, K. Tristram, T. Verhoelst, S. Wolf, J. Young

Abstract: 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 pa… ▽ More 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. △ Less

Submitted 19 April, 2012; originally announced April 2012.

Comments: Acccepted to Astronomy and Astrophysics Review

3-D deconvolution of hyper-spectral astronomical data

Authors: S. Bongard, F. Soulez, E. Thiebaut, E. Pécontal

Abstract: 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 subtractio… ▽ More 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. △ Less

Submitted 20 July, 2011; originally announced July 2011.

Comments: 14 pages, 12 figures, 3 tables. Accepted for publication in MNRAS

Image reconstruction in optical interferometry: Benchmarking the regularization

Authors: Stéphanie Renard, Eric Thiébaut, Fabien Malbet

Abstract: 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 regulariz… ▽ More 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. △ Less

Submitted 22 June, 2011; originally announced June 2011.

Comments: 15 pages, 16 figures; accepted in A&A

A low optical depth region in the inner disk of the HerbigAe star HR5999

Authors: M. Benisty, S. Renard, A. Natta, J. P. Berger, F. Massi, F. Malbet, P. J. V. Garcia, A. Isella, A. Mérand, J. L. Monin, L. Testi, E. Thiébaut, M. Vannier, G. Weigelt

Abstract: 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 re… ▽ More 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. △ Less

Submitted 21 June, 2011; originally announced June 2011.

Comments: 11 pages, 10 figures

Journal ref: 2011, A&A, 531A, 84B

First astronomical unit scale image of the GW Ori triple. Direct detection of a new stellar companion

Authors: J. -P. Berger, J. D. Monnier, R. Millan-Gabet, S. Renard, E. Pedretti, W. Traub, C. Bechet, M. Benisty, N. Carleton, P. Haguenauer, P. Kern, P. Labeye, F. Longa, M. Lacasse, F. Malbet, K. Perraut, S. Ragland, P. Schloerb, P. A. Schuller, E. Thiébaut

Abstract: 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… ▽ More 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 $ρ\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. △ Less

Submitted 20 March, 2011; originally announced March 2011.

Comments: 5 pages, 9 figures, accepted Astronomy and Astrophysics Letters. 2011

Milli-arcsecond images of the Herbig Ae star HD 163296

Authors: S. Renard, F. Malbet, M. Benisty, E. Thiébaut, J. -P. Berger

Abstract: 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 measurem… ▽ More 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. △ Less

Submitted 28 July, 2010; v1 submitted 17 July, 2010; originally announced July 2010.

Comments: 11 pages, 10 figures, accepted A&A paper

Imaging the spotty surface of Betelgeuse in the H band

Authors: X. Haubois, G. Perrin, S. Lacour, T. Verhoelst, S. Meimon, L. Mugnier, E. Thiebaut, J. P. Berger, S. T. Ridgway, J. D. Monnier, R. Millan-Gabet, W. Traub

Abstract: 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… ▽ More 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. △ Less

Submitted 9 November, 2009; v1 submitted 21 October, 2009; originally announced October 2009.

Comments: 10 pages, 10 figures, accepted for publication in A&A, references added

The Pulsation of Chi Cygni Imaged by Optical Interferometry; a Novel Technique to Derive Distance and Mass of Mira Stars

Authors: S. Lacour, E. Thiébaut, G. Perrin, S. Meimon, X. Haubois, E. Pedretti, S. Ridgway, J. D. Monnier, J. P. Berger, P. A. Schuller, H. Woodruff, A. Poncelet, H. Le Coroller, R. Millan-Gabet, M. Lacasse, W. Traub

Abstract: We present infrared interferometric imaging of the S-type Mira star Chi Cygni. The object was observed at four different epochs in 2005-2006 with the IOTA optical interferometer (H band). Images show up to 40% variation in the stellar diameter, as well as significant changes in the limb darkening and stellar inhomogeneities. Model fitting gave precise time-dependent values of the stellar diamete… ▽ More We present infrared interferometric imaging of the S-type Mira star Chi Cygni. The object was observed at four different epochs in 2005-2006 with the IOTA optical interferometer (H band). Images show up to 40% variation in the stellar diameter, as well as significant changes in the limb darkening and stellar inhomogeneities. Model fitting gave precise time-dependent values of the stellar diameter, and reveals presence and displacement of a warm molecular layer. The star radius, corrected for limb darkening, has a mean value of 12.1 mas and shows a 5.1mas amplitude pulsation. Minimum diameter was observed at phase 0.94+/-0.01. Maximum temperature was observed several days later at phase 1.02+/-0.02. We also show that combining the angular acceleration of the molecular layer with CO (Delta v = 3) radial velocity measurements yields a 5.9+/-1.5 mas parallax. The constant acceleration of the CO molecules -- during 80% of the pulsation cycle -- lead us to argument for a free-falling layer. The acceleration is compatible with a gravitational field produced by a 2.1(+1.5/-0.7) solar mass star. This last value is in agreement with fundamental mode pulsator models. We foresee increased development of techniques consisting in combining radial velocity with interferometric angular measurements, ultimately allowing total mapping of the speed, density, and position of the diverse species in pulsation driven atmospheres. △ Less

Submitted 20 October, 2009; originally announced October 2009.

Comments: 36 pages, accepted in ApJ

Journal ref: Astrophys.J.707:632-643,2009

Image Reconstruction in Optical Interferometry

Authors: Éric Thiébaut, Jean-François Giovannelli

Abstract: This tutorial paper describes the problem of image reconstruction from interferometric data with a particular focus on the specific problems encountered at optical (visible/IR) wavelengths. The challenging issues in image reconstruction from interferometric data are introduced in the general framework of inverse problem approach. This framework is then used to describe existing image reconstruct… ▽ More This tutorial paper describes the problem of image reconstruction from interferometric data with a particular focus on the specific problems encountered at optical (visible/IR) wavelengths. The challenging issues in image reconstruction from interferometric data are introduced in the general framework of inverse problem approach. This framework is then used to describe existing image reconstruction algorithms in radio interferometry and the new methods specifically developed for optical interferometry. △ Less

Submitted 12 September, 2009; v1 submitted 11 September, 2009; originally announced September 2009.

Comments: accepted for publication in IEEE Signal Processing Magazine

Journal ref: IEEE Sig.Proc.Mag.27:97-109, 2010

Probing magnetic fields with multi-frequency polarized synchrotron emission

Authors: Jerome Thiebaut, Simon Prunet, Christophe Pichon, Eric Thiébaut

Abstract: We investigate the problem of probing the local spatial structure of the magnetic field of the interstellar medium using multi-frequency polarized maps of the synchrotron emission at radio wavelengths. We focus in this paper on the three-dimensional reconstruction of the largest scales of the magnetic field, relying on the internal depolarization (due to differential Faraday rotation) of the emi… ▽ More We investigate the problem of probing the local spatial structure of the magnetic field of the interstellar medium using multi-frequency polarized maps of the synchrotron emission at radio wavelengths. We focus in this paper on the three-dimensional reconstruction of the largest scales of the magnetic field, relying on the internal depolarization (due to differential Faraday rotation) of the emitting medium as a function of electromagnetic frequency. We argue that multi-band spectroscopy in the radio wavelengths, developed in the context of high-redshift extragalactic HI lines, can be a very useful probe of the 3D magnetic field structure of our Galaxy when combined with a Maximum A Posteriori reconstruction technique. When starting from a fair approximation of the magnetic field, we are able to recover the true one by using a linearized version of the corresponding inverse problem. The spectral analysis of this problem allows us to specify the best sampling strategy in electromagnetic frequency and predicts a spatially anisotropic distribution of posterior errors. The reconstruction method is illustrated for reference fields extracted from realistic magneto-hydrodynamical simulations. △ Less

Submitted 30 November, 2009; v1 submitted 11 September, 2009; originally announced September 2009.

A binary engine fuelling HD87643' s complex circumstellar environment, using AMBER/VLTI

Authors: Florentin Millour, Olivier Chesneau, Marcelo Borges Fernandes, Anthony Meilland, Gilbert Mars, C. Benoist, E. Thiébaut, Philippe Stee, K. -H. Hofmann, Fabien Baron, John R. Young, Philippe Bendjoya, A. C. Carciofi, Armando Domiciano De Souza, Thomas Driebe, Slobodan Jankov, Pierre Kervella, R. G. Petrov, Sylvie Robbe-Dubois, Farrokh Vakili, L. B. F. M. Waters, Gerd Weigelt

Abstract: Context. The star HD 87643, exhibiting the "B[e] phenomenon", has one of the most extreme infrared excesses for this object class. It harbours a large amount of both hot and cold dust, and is surrounded by an extended reflection nebula. Aims. One of our major goals was to investigate the presence of a companion in HD87643. In addition, the presence of close dusty material was tested through a co… ▽ More Context. The star HD 87643, exhibiting the "B[e] phenomenon", has one of the most extreme infrared excesses for this object class. It harbours a large amount of both hot and cold dust, and is surrounded by an extended reflection nebula. Aims. One of our major goals was to investigate the presence of a companion in HD87643. In addition, the presence of close dusty material was tested through a combination of multi-wavelength high spatial 5Aresolution observations. Methods. We observed HD 87643 with high spatial resolution techniques, using the near-IR AMBER/VLTI interferometer with baselines ranging from 60 m to 130 m and the mid-IR MIDI/VLTI interferometer with baselines ranging from 25 m to 65 m. These observations are complemented by NACO/VLT adaptive-optics-corrected images in the K and L-bands, ESO-2.2m optical Wide-Field Imager large-scale images in the B, V and R-bands, Results. We report the direct detection of a companion to HD 87643 by means of image synthesis using the AMBER/VLTI instrument. The presence of the companion is confirmed by the MIDI and NACO data, although with a lower confidence. The companion is separated by ~ 34 mas with a roughly north-south orientation. The period must be large (several tens of years) and hence the orbital parameters are not determined yet. Binarity with high eccentricity might be the key to interpreting the extreme characteristics of this system, namely a dusty circumstellar envelope around the primary, a compact dust nebulosity around the binary system and a complex extended nebula witnessing past violent ejections. △ Less

Submitted 22 September, 2009; v1 submitted 3 August, 2009; originally announced August 2009.

Journal ref: Astronomy and Astrophysics (2009) xxx

Pre-maximum spectro-imaging of the Mira star T Lep with AMBER/VLTI

Authors: Jean-Baptiste Le Bouquin, Sylvestre Lacour, Stephanie Renard, Eric Thiebaut, Antoine Merand, Tijl Verhoelst

Abstract: 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 interferome… ▽ More 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. △ Less

Submitted 23 February, 2009; originally announced February 2009.

Comments: Accepted in A&A

ASKI: full-sky lensing map making algorithms

Authors: C. Pichon, E. Thiebaut, S. Prunet, K. Benabed, S. Colombi, T. Sousbie, R. Teyssier

Abstract: 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 i… ▽ More 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. △ Less

Submitted 25 August, 2009; v1 submitted 14 January, 2009; originally announced January 2009.

Comments: 24 pages, 21 figures accepted for publication to MNRAS;

Phase Closure Image Reconstruction for Future VLTI Instrumentation

Authors: Mercedes E. Filho, Stephanie Renard, Paulo Garcia, Gilles Duvert, Gaspard Duchene, Eric Thiebaut, John Young, Olivier Absil, Jean-Phillipe Berger, Thomas Beckert, Sebastian Hoenig, Dieter Schertl, Gerd Weigelt, Leonardo Testi, Eric Tatuli, Virginie Borkowski, Michael de Becker, Jean Surdej, Bernard Aringer, Joseph Hron, Thomas Lebzelter, Andrea Chiavassa, Romano Corradi, Tim Harries

Abstract: 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 influe… ▽ More 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. △ Less

Submitted 2 October, 2008; originally announced October 2008.

Comments: 11 pages, 26 figures

Phase Referencing in Optical Interferometry

Authors: Mercedes E. Filho, Paulo Garcia, Gilles Duvert, Gaspard Duchene, Eric Thiebaut, John Young, Olivier Absil, Jean-Phillipe Berger, Thomas Beckert, Sebastian Hoenig, Dieter Schertl, Gerd Weigelt, Leonardo Testi, Eric Tatuli, Virginie Borkowski, Michael de Becker, Jean Surdej, Bernard Aringer, Joseph Hron, Thomas Lebzelter, Andrea Chiavassa, Romano Corradi, Tim Harries

Abstract: 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 curren… ▽ More 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) △ Less

Submitted 2 October, 2008; originally announced October 2008.

Comments: 11 pages, 36 figures