VIGAN: Missing View Imputation with Generative Adversarial Networks

In an era when big data are becoming the norm, there is less concern with the quantity but more with the quality and completeness of the data. In many disciplines, data are collected from heterogeneous sources, resulting in multi-view or multi-modal datasets. The missing data problem has been challenging to address in multi-view data analysis. Especially, when certain samples miss an entire view of data, it creates the missing view problem. Classic multiple imputations or matrix completion methods are hardly effective here when no information can be based on in the specific view to impute data for such samples. The commonly-used simple method of removing samples with a missing view can dramatically reduce sample size, thus diminishing the statistical power of a subsequent analysis. In this paper, we propose a novel approach for view imputation via generative adversarial networks (GANs), which we name by VIGAN. This approach first treats each view as a separate domain and identifies domain-to-domain mappings via a GAN using randomly-sampled data from each view, and then employs a multi-modal denoising autoencoder (DAE) to reconstruct the missing view from the GAN outputs based on paired data across the views. Then, by optimizing the GAN and DAE jointly, our model enables the knowledge integration for domain mappings and view correspondences to effectively recover the missing view. Empirical results on benchmark datasets validate the VIGAN approach by comparing against the state of the art. The evaluation of VIGAN in a genetic study of substance use disorders further proves the effectiveness and usability of this approach in life science.Comment: 10 pages, 8 figures, conferenc

Characterizing HR3549B using SPHERE

Aims. In this work, we characterize the low mass companion of the A0 field star HR3549. Methods. We observed HR3549AB in imaging mode with the the NIR branch (IFS and IRDIS) of SPHERE@VLT, with IFS in YJ mode and IRDIS in the H band. We also acquired a medium resolution spectrum with the IRDIS long slit spectroscopy mode. The data were reduced using the dedicated SPHERE GTO pipeline, purposely designed for this instrument. We employed algorithms such as PCA and TLOCI to reduce the speckle noise. Results. The companion was clearly visible both with IRDIS and IFS.We obtained photometry in four different bands as well as the astrometric position for the companion. Based on our astrometry, we confirm that it is a bound object and put constraints on its orbit. Although several uncertainties are still present, we estimate an age of ~100-150 Myr for this system, yielding a most probable mass for the companion of 40-50MJup and T_eff ~300-2400 K. Comparing with template spectra points to a spectral type between M9 and L0 for the companion, commensurate with its position on the color-magnitude diagram.Comment: Accepted by A&A, 13 pages, 10 Figures (Figures 9 and 10 degraded to reduce the dimension

High-contrast spectroscopy of SCR J1845-6357 B

Spectral characterization of sub-stellar companions is essential to understand their composition and formation processes. However, the large contrast ratio of the brightness of each object to that of its parent star limits our ability to extract a clean spectrum, free from any significant contribution from the star. During the development of the long slit spectroscopy (LSS) mode of IRDIS, the dual-band imager and spectrograph of SPHERE, we proposed a data analysis method to estimate and remove the contributions of the stellar spectrum. This method has never been tested on real data because of the lack of instrumentation capable of combining adaptive optics (AO), coronagraphy, and LSS. Nonetheless, a similar attenuation of the star can be obtained using a particular observing configuration. Test data were acquired using the AO-assisted spectrograph VLT/NACO. We obtained new J- and H-band spectra of SCR J1845-6357 B, a T6 companion to a nearby (3.85\pm0.02 pc) M8 star. This system is a well-suited benchmark as it is relatively wide (~1.0") with a modest contrast ratio (~4 mag), and a previously published JHK spectrum is available for reference. We demonstrate that (1) our method is efficient at estimating and removing the stellar contribution, (2) it allows to properly recover the spectral shape of the companion, and (3) it is essential to obtain an unbiased estimation of physical parameters. We also show that the slit configuration associated with this method allows us to use long exposure times with high throughput producing high signal-to-noise ratio data. However, the signal of the companion gets over-subtracted, particularly in our J-band data, compelling us to use a fake companion spectrum to estimate and compensate for the loss of flux. Finally, we report a new astrometric measurement of the position of the companion (sep = 0.817", PA = 227.92 deg).Comment: 11 pages, 8 figures, 4 tables. Accepted for publication in A&

Apodization in high-contrast long-slit spectroscopy. Closer, deeper, fainter, cooler

The spectroscopy of faint planetary-mass companions to nearby stars is one of the main challenges that new-generation high-contrast spectro-imagers are going to face. In a previous work we presented a long slit coronagraph (LSC), for which the presence of a slit in the coronagraphic focal plane induces a complex distribution of energy in the Lyot pupil-plane that cannot be easily masked with a binary Lyot stop. To alleviate this concern, we propose to use a pupil apodization to suppress diffraction, creating an apodized long slit coronagraph (ALSC). After describing how the apodization is optimized, we demonstrate its advantages with respect to the CLC in the context of SPHERE/IRDIS long slit spectroscopy (LSS) mode at low-resolution with a 0.12" slit and 0.18" coronagraphic mask. We perform different sets of simulations with and without aberrations, and with and without a slit to demonstrate that the apodization is a more appropriate concept for LSS, at the expense of a significantly reduced throughput (37%) compared to the LSC. Then we perform detailed end-to-end simulations of the LSC and the ALSC that include realistic levels of aberrations to obtain datasets representing 1h of integration time on stars of spectral types A0 to M0 located at 10 pc. We insert spectra of planetary companions at different effective temperatures (Teff) and surface gravities (log g) into the data at angular separations of 0.3" to 1.5" and with contrast ratios from 6 to 18 mag. Using the SD method to subtract the speckles, we show that the ALSC brings a gain in sensitivity of up to 3 mag at 0.3" with respect to the LSC, which leads to a much better spectral extraction below 0.5". In terms of Teff, we demonstrate that at small angular separations the limit with the ALSC is always lower by at least 100K, inducing an increase of sensitivity of a factor up to 1.8 in objects' masses at young ages. [Abridged]Comment: 15 pages, 17 figures. Accepted for publication in A&

First light of the VLT planet finder SPHERE. I. Detection and characterization of the sub-stellar companion GJ 758 B

GJ758 B is a brown dwarf companion to a nearby (15.76 pc) solar-type, metal-rich (M/H = +0.2 dex) main-sequence star (G9V) that was discovered with Subaru/HiCIAO in 2009. From previous studies, it has drawn attention as being the coldest (~600K) companion ever directly imaged around a neighboring star. We present new high-contrast data obtained during the commissioning of the SPHERE instrument at the VLT. The data was obtained in Y-, J-, H-, and Ks-bands with the dual-band imaging (DBI) mode of IRDIS, providing a broad coverage of the full near-infrared (near-IR) range at higher contrast and better spectral sampling than previously reported. In this new set of high-quality data, we report the re-detection of the companion, as well as the first detection of a new candidate closer-in to the star. We use the new 8 photometric points for an extended comparison of GJ758 B with empirical objects and 4 families of atmospheric models. From comparison to empirical object, we estimate a T8 spectral type, but none of the comparison object can accurately represent the observed near-IR fluxes of GJ758 B. From comparison to atmospheric models, we attribute a Teff = 600K $\pm$ 100K, but we find that no atmospheric model can adequately fit all the fluxes of GJ758 B. The photometry of the new candidate companion is broadly consistent with L-type objects, but a second epoch with improved photometry is necessary to clarify its status. The new astrometry of GJ758 B shows a significant proper motion since the last epoch. We use this result to improve the determination of the orbital characteristics using two fitting approaches, Least-Square Monte Carlo and Markov Chain Monte Carlo. Finally, we analyze the sensitivity of our data to additional closer-in companions and reject the possibility of other massive brown dwarf companions down to 4-5 AU. [abridged]Comment: 20 pages, 15 figures. Accepted for publication in A&

High-contrast imaging of Sirius~A with VLT/SPHERE: Looking for giant planets down to one astronomical unit

Sirius has always attracted a lot of scientific interest, especially after the discovery of a companion white dwarf at the end of the 19th century. Very early on, the existence of a potential third body was put forward to explain some of the observed properties of the system. We present new coronagraphic observations obtained with VLT/SPHERE that explore, for the very first time, the innermost regions of the system down to 0.2" (0.5 AU) from Sirius A. Our observations cover the near-infrared from 0.95 to 2.3 $\mu$m and they offer the best on-sky contrast ever reached at these angular separations. After detailing the steps of our SPHERE/IRDIFS data analysis, we present a robust method to derive detection limits for multi-spectral data from high-contrast imagers and spectrographs. In terms of raw performance, we report contrasts of 14.3 mag at 0.2", ~16.3 mag in the 0.4-1.0" range and down to 19 mag at 3.7". In physical units, our observations are sensitive to giant planets down to 11 $M_{Jup}$ at 0.5 AU, 6-7 $M_{Jup}$ in the 1-2 AU range and ~4 $M_{Jup}$ at 10 AU. Despite the exceptional sensitivity of our observations, we do not report the detection of additional companions around Sirius A. Using a Monte Carlo orbital analysis, we show that we can reject, with about 50% probability, the existence of an 8 $M_{Jup}$ planet orbiting at 1 AU. In addition to the results presented in the paper, we provide our SPHERE/IFS data reduction pipeline at http://people.lam.fr/vigan.arthur/ under the MIT license.Comment: 16 pages, 10 figures, accepted for publication in MNRA

Calibration of quasi-static aberrations in exoplanet direct-imaging instruments with a Zernike phase-mask sensor. II. Concept validation with ZELDA on VLT/SPHERE

Warm or massive gas giant planets, brown dwarfs, and debris disks around nearby stars are now routinely observed by dedicated high-contrast imaging instruments on large, ground-based observatories. These facilities include extreme adaptive optics (ExAO) and state-of-the-art coronagraphy to achieve unprecedented sensitivities for exoplanet detection and spectral characterization. However, differential aberrations between the ExAO sensing path and the science path represent a critical limitation for the detection of giant planets with a contrast lower than a few $10^{-6}$ at very small separations (<0.3\as) from their host star. In our previous work, we proposed a wavefront sensor based on Zernike phase contrast methods to circumvent this issue and measure these quasi-static aberrations at a nanometric level. We present the design, manufacturing and testing of ZELDA, a prototype that was installed on VLT/SPHERE during its reintegration in Chile. Using the internal light source of the instrument, we performed measurements in the presence of Zernike or Fourier modes introduced with the deformable mirror. Our experimental and simulation results are consistent, confirming the ability of our sensor to measure small aberrations (<50 nm rms) with nanometric accuracy. We then corrected the long-lived non-common path aberrations in SPHERE based on ZELDA measurements. We estimated a contrast gain of 10 in the coronagraphic image at 0.2\as, reaching the raw contrast limit set by the coronagraph in the instrument. The simplicity of the design and its phase reconstruction algorithm makes ZELDA an excellent candidate for the on-line measurements of quasi-static aberrations during the observations. The implementation of a ZELDA-based sensing path on the current and future facilities (ELTs, future space missions) could ease the observation of the cold gaseous or massive rocky planets around nearby stars.Comment: 13 pages, 12 figures, A&A accepted on June 3rd, 2016. v2 after language editin

SPHERE: the exoplanet imager for the Very Large Telescope

Observations of circumstellar environments to look for the direct signal of exoplanets and the scattered light from disks has significant instrumental implications. In the past 15 years, major developments in adaptive optics, coronagraphy, optical manufacturing, wavefront sensing and data processing, together with a consistent global system analysis have enabled a new generation of high-contrast imagers and spectrographs on large ground-based telescopes with much better performance. One of the most productive is the Spectro-Polarimetic High contrast imager for Exoplanets REsearch (SPHERE) designed and built for the ESO Very Large Telescope (VLT) in Chile. SPHERE includes an extreme adaptive optics system, a highly stable common path interface, several types of coronagraphs and three science instruments. Two of them, the Integral Field Spectrograph (IFS) and the Infra-Red Dual-band Imager and Spectrograph (IRDIS), are designed to efficiently cover the near-infrared (NIR) range in a single observation for efficient young planet search. The third one, ZIMPOL, is designed for visible (VIR) polarimetric observation to look for the reflected light of exoplanets and the light scattered by debris disks. This suite of three science instruments enables to study circumstellar environments at unprecedented angular resolution both in the visible and the near-infrared. In this work, we present the complete instrument and its on-sky performance after 4 years of operations at the VLT.Comment: Final version accepted for publication in A&

The SPHERE data center: a reference for high contrast imaging processing

The objective of the SPHERE Data Center is to optimize the scientific return of SPHERE at the VLT, by providing optimized reduction procedures, services to users and publicly available reduced data. This paper describes our motivation, the implementation of the service (partners, infrastructure and developments), services, description of the on-line data, and future developments. The SPHERE Data Center is operational and has already provided reduced data with a good reactivity to many observers. The first public reduced data have been made available in 2017. The SPHERE Data Center is gathering a strong expertise on SPHERE data and is in a very good position to propose new reduced data in the future, as well as improved reduction procedures.Comment: SF2A proceeding

Spectroscopy across the brown dwarf/planetary mass boundary - I. Near-infrared JHK spectra

With a uniform VLT SINFONI data set of nine targets, we have developed an empirical grid of J,H,K spectra of the atmospheres of objects estimated to have very low substellar masses of \sim5-20 MJup and young ages of \sim1-50 Myr. Most of the targets are companions, objects which are especially valuable for comparison with atmosphere and evolutionary models, as they present rare cases in which the age is accurately known from the primary. Based on the sample youth, all objects are expected to have low surface gravity, and this study investigates the critical early phases of the evolution of substellar objects. The spectra are compared with grids of five different theoretical atmosphere models. This analysis represents the first systematic model comparison with infrared spectra of young brown dwarfs. The fits to the full JHK spectra of each object result in a range of best fit effective temperatures of +/-150-300K whether or not the full model grid or a subset restricted to lower log(g) values is used. This effective temperature range is significantly larger than the uncertainty typically assigned when using a single model grid. Fits to a single wavelength band can vary by up to 1000K using the different models. Since the overall shape of these spectra is governed more by the temperature than surface gravity, unconstrained model fits did not find matches with low surface gravity or a trend in log(g) with age. This suggests that empirical comparison with spectra of unambiguously young objects targets (such as these SINFONI data) may be the most reliable method to search for indications of low surface gravity and youth. For two targets, the SINFONI data are a second epoch and the data show no variations in morphology over time. The analysis of two other targets, AB Pic B and CT Cha B, suggests that these objects may have lower temperatures, and consequently lower masses, than previously estimated.Comment: 15 pages, 13 figure