Source author record

David Mary

David Mary appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2022arXiv

3-beam self-calibrated Kernel nulling photonic interferometer

The use of interferometric nulling for the direct characterization of extrasolar planets is an exciting prospect, but one that faces many practical challenges when deployed on telescopes. The largest limitation is the extreme sensitivity of nullers to any residual optical path differences between the incoming telescope beams even after adaptive optics or fringe-tracker correction. The recently proposed kernel-nulling architecture attempts to alleviate this by producing the destructive interference required for nulling, in a scheme whereby self-calibrated observables can be created efficiently, in effect canceling out residual atmospheric piston terms. Here we experimentally demonstrate for the first time a successful creation of self-calibrated kernel-null observables for nulling interferometry in the laboratory. We achieved this through the use of a purpose-built photonic integrated device, containing a multimode interference coupler that creates one bright, and two nulled outputs when injected with three co-phased telescope beams. The device produces the nulled outputs in a way that, by the subtraction of the measured output flux, create a single self-calibrated kernel-null. We experimentally demonstrate the extraction of kernel-nulls for up to 200 nm induced piston error using a laboratory test-bench at a wavelength of 1.55 μm. Further, we empirically demonstrate the kernel-null behaviour when injected with a binary companion analogue equivalent to a 2.32 mas separation at a contrast of 10^{-2}, under 100 nm RMS upstream piston residuals.

preprint2022arXiv

Scaling relations of convective granulation noise across the HR diagram from 3D stellar atmosphere models

High-precision photometric data from space missions have improved our understanding of stellar granulation. These observations have shown with precision the stochastic brightness fluctuations of stars across the HR diagram, allowing us to better understand how stellar surface convection reacts to a change in stellar parameters. These fluctuations need to be understood and quantified in order to improve the detection and characterization of exoplanets. In this work, we provide new scaling relations of two characteristic properties of the brightness fluctuations time series, the standard deviation ($σ$) and the auto-correlation time ($τ\rm_{eff}$). This was done by using long time series of 3D stellar atmosphere models at different metallicities and across the HR diagram, generated with a 3D radiative hydrodynamical code: the STAGGER code. We compared our synthetic granulation properties with the values of a large sample of Kepler stars, and analyzed selected stars with accurate stellar parameters from the Kepler LEGACY sample. Our 3D models showed that $σ\proptoν\rm_{max}^{-0.567\pm0.012}$ and $τ\rm_{eff}\proptoν\rm_{max}^{-0.997\pm0.018}$ for stars at solar metallicity. We showed that both $σ$ and $τ\rm_{eff}$ decrease with metallicity, although the metallicity dependence is more significant on $σ$. Unlike previous studies, we found very good agreement between $σ$ from Kepler targets and the 3D models at $\log{g}\leq3.5$, and a good correlation between the stars and models with $\log{g}\geq3.5$. For $τ\rm_{eff}$, we found that the 3D models reproduced well the Kepler LEGACY star values. Overall, this study shows that 3D stellar atmosphere models reproduce the granulation properties of stars across the HR diagram.

preprint2022arXiv

Statistical tests with multi-wavelength Kernel-phase analysis for the detection and characterization of planetary companions

Kernel phase is a method to interpret stellar point source images by considering their formation as the analytical result of an interferometric process. Using Fourier formalism, this method allows for observing planetary companions around nearby stars at separations down to half a telescope resolution element, typically 20\,mas for a 8\,m class telescope in H band. The Kernel-phase analysis has so far been mainly focused on working with a single monochromatic light image, recently providing theoretical contrast detection limits down to $10^{-4}$ at 200\,mas with JWST/NIRISS in the mid-infrared by using hypothesis testing theory. In this communication, we propose to extend this approach to data cubes provided by integral field spectrographs (IFS) on ground-based telescopes with adaptive optics to enhance the detection of planetary companions and explore the spectral characterization of their atmosphere by making use of the Kernel-phase multi-spectral information. Using ground-based IFS data cube with a spectral resolution R=20, we explore different statistical tests based on kernel phases at three wavelengths to estimate the detection limits for planetary companions. Our tests are first conducted with synthetic data before extending their use to real images from ground-based exoplanet imagers such as Subaru/SCExAO and VLT/SPHERE in the near future. Future applications to multi-wavelength data from space telescopes are also discussed for the observation of planetary companions with JWST.

preprint2020arXiv

Angular differential kernel phases

To reach its optimal performance, Fizeau interferometry requires that we work to resolve instrumental biases through calibration. One common technique used in high contrast imaging is angular differential imaging, which calibrates the point spread function and flux leakage using a rotation in the focal plane. Our aim is to experimentally demonstrate and validate the efficacy of an angular differential kernel-phase approach, a new method for self-calibrating interferometric observables that operates similarly to angular differential imaging, while retaining their statistical properties. We used linear algebra to construct new observables that evolve outside of the subspace spanned by static biases. On-sky observations of a binary star with the SCExAO instrument at the Subaru telescope were used to demonstrate the practicality of this technique. We used a classical approach on the same data to compare the effectiveness of this method. The proposed method shows smaller and more Gaussian residuals compared to classical calibration methods, while retaining compatibility with the statistical tools available. We also provide a measurement of the stability of the SCExAO instrument that is relevant to the application of the technique. Angular differential kernel phases provide a reliable method for calibrating biased observables. Although the sensitivity at small separations is reduced for small field rotations, the calibration is effectively improved and the number of subjective choices is reduced.

preprint2020arXiv

Kernel-phase analysis: aperture modeling prescriptions that minimize calibration errors

Kernel-phase is a data analysis method based on a generalization of the notion of closure-phase invented in the context of interferometry, but that applies to well corrected diffraction dominated images produced by an arbitrary aperture. The linear model upon which it relies theoretically leads to the formation of observable quantities robust against residual aberrations. In practice, detection limits reported thus far seem to be dominated by systematic errors induced by calibration biases not sufficiently filtered out by the kernel projection operator. This paper focuses on the impact the initial modeling of the aperture has on these errors and introduces a strategy to mitigate them, using a more accurate aperture transmission model. The paper first uses idealized monochromatic simulations of a non trivial aperture to illustrate the impact modeling choices have on calibration errors. It then applies the outlined prescription to two distinct data-sets of images whose analysis has previously been published. The use of a transmission model to describe the aperture results in a significant improvement over the previous type of analysis. The thus reprocessed data-sets generally lead to more accurate results, less affected by systematic errors. As kernel-phase observing programs are becoming more ambitious, accuracy in the aperture description is becoming paramount to avoid situations where contrast detection limits are dominated by systematic errors. Prescriptions outlined in this paper will benefit any attempt at exploiting kernel-phase for high-contrast detection.

preprint2020arXiv

ORIGIN: Blind detection of faint emission line galaxies in MUSE datacubes

One of the major science cases of the MUSE integral field spectrograph is the detection of Lyman-alpha emitters at high redshifts. The on-going and planned deep fields observations will allow for one large sample of these sources. An efficient tool to perform blind detection of faint emitters in MUSE datacubes is a prerequisite of such an endeavor. Several line detection algorithms exist but their performance during the deepest MUSE exposures is hard to quantify, in particular with respect to their actual false detection rate, or purity. {The aim of this work is to design and validate} an algorithm that efficiently detects faint spatial-spectral emission signatures, while allowing for a stable false detection rate over the data cube and providing in the same time an automated and reliable estimation of the purity. Results on simulated data cubes providing ground truth show that the method reaches its aims in terms of purity and completeness. When applied to the deep 30-hour exposure MUSE datacube in the Hubble Ultra Deep Field, the algorithms allows for the confirmed detection of 133 intermediate redshifts galaxies and 248 Lyman Alpha Emitters, including 86 sources with no HST counterpart. The algorithm fulfills its aims in terms of detection power and reliability. It is consequently implemented as a Python package whose code and documentation are available on GitHub and readthedocs.

preprint2016arXiv

Distributed multi-frequency image reconstruction for radio-interferometry

The advent of enhanced technologies in radio interferometry and the perspective of the SKA telescope bring new challenges in image reconstruction. One of these challenges is the spatio-spectral reconstruction of large (Terabytes) data cubes with high fidelity. This contribution proposes an alternative implementation of one such 3D prototype algorithm, MUFFIN (MUlti-Frequency image reconstruction For radio INterferometry), which combines spatial and spectral analysis priors. Using a recently proposed primal dual algorithm, this new version of MUFFIN allows a parallel implementation where computationally intensive steps are split by spectral channels. This parallelization allows to implement computationally demanding translation invariant wavelet transforms (IUWT), as opposed to the union of bases used previously. This alternative implementation is important as it opens the possibility of comparing these efficient dictionaries, and others, in spatio-spectral reconstruction. Numerical results show that the IUWT-based version can be successfully implemented at large scale with performances comparable to union of bases.

preprint2016arXiv

Multi-frequency image reconstruction for radio interferometry. A regularized inverse problem approach

We describe a "spatio-spectral" deconvolution algorithm for wide-band imaging in radio interferometry. In contrast with the existing multi-frequency reconstruction algorithms, the proposed method does not rely on a model of the sky-brightness spectral distribution. This non-parametric approach can be of particular interest for the new generation of low frequency radiotelescopes. The proposed solution formalizes the reconstruction problem as a convex optimization problem with spatial and spectral regularizations. The efficiency of this approach has been already proven for narrow-band image reconstruction and the present contribution can be considered as its extension to the multi-frequency case. Because the number of frequency bands multiplies the size of the inverse problem, particular attention is devoted to the derivation of an iterative large scale optimization algorithm. It is shown that the main computational bottleneck of the approach, which lies in the resolution of a linear system, can be efficiently overcome by a fully parallel implementation w.r.t. the frequencies, where each processor reconstructs a narrow-band image. All the other optimization steps are extremely fast. A parallel implementation of the algorithm in Julia is publicly available at https://github.com/andferrari. Preliminary simulations illustrate the performances of the method and its ability to reconstruct complex spatio-spectral structures.

preprint2015arXiv

Distributed image reconstruction for very large arrays in radio astronomy

Current and future radio interferometric arrays such as LOFAR and SKA are characterized by a paradox. Their large number of receptors (up to millions) allow theoretically unprecedented high imaging resolution. In the same time, the ultra massive amounts of samples makes the data transfer and computational loads (correlation and calibration) order of magnitudes too high to allow any currently existing image reconstruction algorithm to achieve, or even approach, the theoretical resolution. We investigate here decentralized and distributed image reconstruction strategies which select, transfer and process only a fraction of the total data. The loss in MSE incurred by the proposed approach is evaluated theoretically and numerically on simple test cases.

preprint2015arXiv

Large Scale 3D Image Reconstruction in Optical Interferometry

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.

preprint2014arXiv

Blind and fully constrained unmixing of hyperspectral images

This paper addresses the problem of blind and fully constrained unmixing of hyperspectral images. Unmixing is performed without the use of any dictionary, and assumes that the number of constituent materials in the scene and their spectral signatures are unknown. The estimated abundances satisfy the desired sum-to-one and nonnegativity constraints. Two models with increasing complexity are developed to achieve this challenging task, depending on how noise interacts with hyperspectral data. The first one leads to a convex optimization problem, and is solved with the Alternating Direction Method of Multipliers. The second one accounts for signal-dependent noise, and is addressed with a Reweighted Least Squares algorithm. Experiments on synthetic and real data demonstrate the effectiveness of our approach.

preprint2014arXiv

MORESANE: MOdel REconstruction by Synthesis-ANalysis Estimators. A sparse deconvolution algorithm for radio interferometric imaging

(arXiv abridged abstract) The current years are seeing huge developments of radio telescopes and a tremendous increase of their capabilities. Such systems make mandatory the design of more sophisticated techniques not only for transporting, storing and processing this new generation of radio interferometric data, but also for restoring the astrophysical information contained in such data. In this paper we present a new radio deconvolution algorithm named MORESANE and its application to fully realistic simulated data of MeerKAT, one of the SKA precursors. This method has been designed for the difficult case of restoring diffuse astronomical sources which are faint in brightness, complex in morphology and possibly buried in the dirty beam's side lobes of bright radio sources in the field. MORESANE is a greedy algorithm which combines complementary types of sparse recovery methods in order to reconstruct the most appropriate sky model from observed radio visibilities. A synthesis approach is used for the reconstruction of images, in which the synthesis atoms representing the unknown sources are learned using analysis priors. We apply this new deconvolution method to fully realistic simulations of radio observations of a galaxy cluster and of an HII region in M31. We show that MORESANE is able to efficiently reconstruct images composed from a wide variety of sources from radio interferometric data. Comparisons with other available algorithms, which include multi-scale CLEAN and the recently proposed methods by Li et al. (2011) and Carrillo et al. (2012), indicate that MORESANE provides competitive results in terms of both total flux/surface brightness conservation and fidelity of the reconstructed model. MORESANE seems particularly well suited for the recovery of diffuse and extended sources, as well as bright and compact radio sources known to be hosted in galaxy clusters.

preprint2014arXiv

Non-thermal emission from galaxy clusters: feasibility study with SKA1

Galaxy clusters are known to host a variety of extended radio sources: tailed radio galaxies whose shape is modelled by the interaction with the intra-cluster medium (ICM); radio bubbles filling cavities in the ICM distribution and rising buoyantly through the thermal gas; diffuse giant radio sources ("halos" and "relics") revealing the presence of relativistic electrons and magnetic fields in the intra-cluster volume. It is currently the subject of an active debate how the non-thermal components that we observe at radio wavelengths affect the physical properties of the ICM and depend on the dynamical state of galaxy clusters. In this work we start our SKA1 feasibility study of the "radio cluster zoo" through simulations of a typical radio-loud cluster, hosting several bright tailed radio galaxies and a diffuse radio halo. Realistic simulations of SKA1 observations are obtained through the MeqTrees software. A new deconvolution algorithm, based on sparse representations and optimised for the detection of faint diffuse astronomical sources, is tested and compared to the classical CLEAN method.

preprint2014arXiv

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

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.

preprint2014arXiv

Statistical characterization of polychromatic absolute and differential squared visibilities obtained from AMBER/VLTI instrument

In optical interferometry, the visibility squared modulus are generally assumed to follow a Gaussian distribution and to be independent of each other. A quantitative analysis of the relevance of such assumptions is important to help improving the exploitation of existing and upcoming multi-wavelength interferometric instruments. Analyze the statistical behaviour of both the absolute and the colour-differential squared visibilities: distribution laws, correlations and cross-correlations between different baselines. We use observations of stellar calibrators obtained with AMBER instrument on VLTI in different instrumental and observing configurations, from which we extract the frame-by-frame transfer function. Statistical hypotheses tests and diagnostics are then systematically applied. For both absolute and differential squared visibilities and under all instrumental and observing conditions, we find a better fit for the Student distribution than for the Gaussian, log-normal and Cauchy distributions. We find and analyze clear correlation effects caused by atmospheric perturbations. The differential squared visibilities allow to keep a larger fraction of data with respect to selected absolute squared visibilities and thus benefit from reduced temporal dispersion, while their distribution is more clearly characterized. The frame selection based on the criterion of a fixed SNR value might result in either a biased sample of frames or in a too severe selection.

preprint2012arXiv

Radio Continuum Surveys with Square Kilometre Array Pathfinders

In the lead-up to the Square Kilometre Array (SKA) project, several next-generation radio telescopes and upgrades are already being built around the world. These include APERTIF (The Netherlands), ASKAP (Australia), eMERLIN (UK), VLA (USA), e-EVN (based in Europe), LOFAR (The Netherlands), Meerkat (South Africa), and the Murchison Widefield Array (MWA). Each of these new instruments has different strengths, and coordination of surveys between them can help maximise the science from each of them. A radio continuum survey is being planned on each of them with the primary science objective of understanding the formation and evolution of galaxies over cosmic time, and the cosmological parameters and large-scale structures which drive it. In pursuit of this objective, the different teams are developing a variety of new techniques, and refining existing ones. Here we describe these projects, their science goals, and the technical challenges which are being addressed to maximise the science return.

preprint2012arXiv

The 2012 Interferometric Imaging Beauty Contest

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.