Source author record

J. -B. Le Bouquin

J. -B. Le Bouquin 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

48works
7topics
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

48 published item(s)

preprint2026arXiv

The ExoGRAVITY survey: A K-band spectral library of giant exoplanet and brown dwarf companions

Direct observations of exoplanet and brown dwarf companions with near-infrared interferometry, first enabled by the dual-field mode of VLTI/GRAVITY, provide unique measurements of the objects' orbital motions and atmospheric compositions. Here, we compile a homogeneous library of all exoplanet and brown dwarf K-band spectra observed by GRAVITY thus far. We re-reduced all the available GRAVITY dual-field high-contrast data and, where companions are detected, extract their ~2.0-2.4 $μ$m K-band contrast spectra. We then derived stellar model atmospheres for all employed flux references, which we used to convert the companion contrast into companion flux spectra. Solely from the resulting GRAVITY spectra, we extracted spectral types, spectral indices, and bulk physical properties for all companions. Finally, and with the help of age constraints from the literature, we also derived isochronal masses for most companions using evolutionary models. The resulting library contains R ~ 500 GRAVITY spectra of 39 substellar companions from late M to late T spectral types, including the entire L-T transition. Throughout this transition, a shift from CO-dominated late M- and L-type dwarfs to CH4-dominated T-type dwarfs can be observed in the K-band. The GRAVITY spectra alone constrain the objects' bolometric luminosity to typically within $\pm$0.15 dex. The derived isochronal masses agree with dynamical masses from the literature where available, except for HD 4113 c for which we confirm its previously reported potential underluminosity. Medium-resolution spectroscopy of substellar companions with GRAVITY provides insight into the carbon chemistry and the cloudiness of these objects' atmospheres. It also constrains these objects' bolometric luminosities which can yield measurements of their formation entropy if combined with dynamical masses, for instance from Gaia and GRAVITY astrometry.

preprint2022arXiv

A dusty veil shading Betelgeuse during its Great Dimming

Red supergiants are the most common final evolutionary stage of stars that have initial masses between 8 and 35 times that of the Sun. During this stage, which lasts roughly 100,000 years1, red supergiants experience substantial mass loss. However, the mechanism for this mass loss is unknown. Mass loss may affect the evolutionary path, collapse and future supernova light curve of a red supergiant, and its ultimate fate as either a neutron star or a black hole. From November 2019 to March 2020, Betelgeuse - the second-closest red supergiant to Earth (roughly 220 parsecs, or 724 light years, away) - experienced a historic dimming of its visible brightness. Usually having an apparent magnitude between 0.1 and 1.0, its visual brightness decreased to 1.614 +/- 0.008 magnitudes around 7-13 February 2020 - an event referred to as Betelgeuse's Great Dimming. Here we report high-angular-resolution observations showing that the southern hemisphere of Betelgeuse was ten times darker than usual in the visible spectrum during its Great Dimming. Observations and modelling support a scenario in which a dust clump formed recently in the vicinity of the star, owing to a local temperature decrease in a cool patch that appeared on the photosphere. The directly imaged brightness variations of Betelgeuse evolved on a timescale of weeks. Our findings suggest that a component of mass loss from red supergiants is inhomogeneous, linked to a very contrasted and rapidly changing photosphere

preprint2022arXiv

First Light for GRAVITY Wide: Large Separation Fringe Tracking for the Very Large Telescope Interferometer

GRAVITY+ is the upgrade of GRAVITY and the Very Large Telescope Interferometer (VLTI) with wide-separation fringe tracking, new adaptive optics, and laser guide stars on all four 8~m Unit Telescopes (UTs), for ever fainter, all-sky, high contrast, milliarcsecond interferometry. Here we present the design and first results of the first phase of GRAVITY+, called GRAVITY Wide. GRAVITY Wide combines the dual-beam capabilities of the VLTI and the GRAVITY instrument to increase the maximum separation between the science target and the reference star from 2 arcseconds with the 8 m UTs up to several 10 arcseconds, limited only by the Earth's turbulent atmosphere. This increases the sky-coverage of GRAVITY by two orders of magnitude, opening up milliarcsecond resolution observations of faint objects, and in particular the extragalactic sky. The first observations in 2019 - 2022 include first infrared interferometry of two redshift $z\sim2$ quasars, interferometric imaging on the binary system HD 105913A, and repeated observations of multiple star systems in the Orion Trapezium Cluster. We find the coherence loss between the science object and fringe-tracking reference star well described by the turbulence of the Earth's atmosphere. We confirm that the larger apertures of the UTs result in higher visibilities for a given separation due to larger overlap of the projected pupils on sky and give predictions for visibility loss as a function of separation to be used for future planning.

preprint2022arXiv

GRAVITY faint: reducing noise sources in GRAVITY$^+$ with a fast metrology attenuation system

With the upgrade from GRAVITY to GRAVITY$^+$ the instrument will evolve into an all-sky interferometer that can observe faint targets, such as high redshift AGN. Observing the faintest targets requires reducing the noise sources in GRAVITY as much as possible. The dominant noise source, especially in the blue part of the spectrum, is the backscattering of the metrology laser light onto the detector. To reduce this noise we introduce two new metrology modes. With a combination of small hardware changes and software adaptations, we can dim the metrology laser during the observation without losing the phase referencing. For single beam targets, we can even turn off the metrology laser for the maximum SNR on the detector. These changes lead to an SNR improvement of over a factor of two averaged over the whole spectrum and up to a factor of eight in the part of the spectrum currently dominated by laser noise.

preprint2021arXiv

Constraining the Nature of the PDS 70 Protoplanets with VLTI/GRAVITY

We present K-band interferometric observations of the PDS 70 protoplanets along with their host star using VLTI/GRAVITY. We obtained K-band spectra and 100 $μ$as precision astrometry of both PDS 70 b and c in two epochs, as well as spatially resolving the hot inner disk around the star. Rejecting unstable orbits, we found a nonzero eccentricity for PDS 70 b of $0.17 \pm 0.06$, a near-circular orbit for PDS 70 c, and an orbital configuration that is consistent with the planets migrating into a 2:1 mean motion resonance. Enforcing dynamical stability, we obtained a 95% upper limit on the mass of PDS 70 b of 10 $M_\textrm{Jup}$, while the mass of PDS 70 c was unconstrained. The GRAVITY K-band spectra rules out pure blackbody models for the photospheres of both planets. Instead, the models with the most support from the data are planetary atmospheres that are dusty, but the nature of the dust is unclear. Any circumplanetary dust around these planets is not well constrained by the planets' 1-5 $μ$m spectral energy distributions (SEDs) and requires longer wavelength data to probe with SED analysis. However with VLTI/GRAVITY, we made the first observations of a circumplanetary environment with sub-au spatial resolution, placing an upper limit of 0.3~au on the size of a bright disk around PDS 70 b.

preprint2021arXiv

Deep Images of the Galactic Center with GRAVITY

Stellar orbits at the Galactic Center provide a very clean probe of the gravitational potential of the supermassive black hole. They can be studied with unique precision, beyond the confusion limit of a single telescope, with the near-infrared interferometer GRAVITY. Imaging is essential to search the field for faint, unknown stars on short orbits which potentially could constrain the black hole spin. Furthermore, it provides the starting point for astrometric fitting to derive highly accurate stellar positions. Here, we present $\mathrm{G^R}$, a new imaging tool specifically designed for Galactic Center observations with GRAVITY. The algorithm is based on a Bayesian interpretation of the imaging problem, formulated in the framework of information field theory and building upon existing works in radio-interferometric imaging. Its application to GRAVITY observations from 2021 yields the deepest images to date of the Galactic Center on scales of a few milliarcseconds. The images reveal the complicated source structure within the central $100\,\mathrm{mas}$ around Sgr A*, where we detected the stars S29 and S55 and confirm S62 on its trajectory, slowly approaching Sgr A*. Furthermore, we were able to detect S38, S42, S60, and S63 in a series of exposures for which we offset the fiber from Sgr A*. We provide an update on the orbits of all aforementioned stars. In addition to these known sources, the images also reveal a faint star moving to the west at a high angular velocity. We cannot find any coincidence with any known source and, thus, we refer to the new star as S300. From the flux ratio with S29, we estimate its K-band magnitude as $m_\mathrm{K}\left(\mathrm{S300}\right)\simeq 19.0 - 19.3$. Images obtained with CLEAN confirm the detection.

preprint2021arXiv

Detection of faint stars near SgrA* with GRAVITY

The spin of the supermassive black hole that resides at the Galactic Centre can in principle be measured by accurate measurements of the orbits of stars that are much closer to SgrA* than S2, the orbit of which recently provided the measurement of the gravitational redshift and the Schwarzschild precession. The GRAVITY near-infrared interferometric instrument combining the four 8m telescopes of the VLT provides a spatial resolution of 2-4 mas, breaking the confusion barrier for adaptive-optics-assisted imaging with a single 8-10m telescope. We used GRAVITY to observe SgrA* over a period of six months in 2019 and employed interferometric reconstruction methods developed in radio astronomy to search for faint objects near SgrA*. This revealed a slowly moving star of magnitude 18.9 in K band within 30mas of SgrA*. The position and proper motion of the star are consistent with the previously known star S62, which is at a substantially larger physical distance, but in projection passes close to SgrA*. Observations in August and September 2019 easily detected S29, with K-magnitude of 16.6, at approximately 130 mas from SgrA*. The planned upgrades of GRAVITY, and further improvements in the calibration, hold the promise of finding stars fainter than magnitude 19 at K.

preprint2021arXiv

Improved GRAVITY astrometric accuracy from modeling of optical aberrations

The GRAVITY instrument on the ESO VLTI pioneers the field of high-precision near-infrared interferometry by providing astrometry at the $10 - 100\,μ$as level. Measurements at such high precision crucially depend on the control of systematic effects. Here, we investigate how aberrations introduced by small optical imperfections along the path from the telescope to the detector affect the astrometry. We develop an analytical model that describes the impact of such aberrations on the measurement of complex visibilities. Our formalism accounts for pupil-plane and focal-plane aberrations, as well as for the interplay between static and turbulent aberrations, and successfully reproduces calibration measurements of a binary star. The Galactic Center observations with GRAVITY in 2017 and 2018, when both Sgr A* and the star S2 were targeted in a single fiber pointing, are affected by these aberrations at a level of less than 0.5 mas. Removal of these effects brings the measurement in harmony with the dual beam observations of 2019 and 2020, which are not affected by these aberrations. This also resolves the small systematic discrepancies between the derived distance $R_0$ to the Galactic Center reported previously.

preprint2021arXiv

The ExoGRAVITY project: using single mode interferometry to characterize exoplanets

Combining adaptive optics and interferometric observations results in a considerable contrast gain compared to single-telescope, extreme AO systems. Taking advantage of this, the ExoGRAVITY project is a survey of known young giant exoplanets located in the range of 0.1'' to 2'' from their stars. The observations provide astrometric data of unprecedented accuracy, being crucial for refining the orbital parameters of planets and illuminating their dynamical histories. Furthermore, GRAVITY will measure non-Keplerian perturbations due to planet-planet interactions in multi-planet systems and measure dynamical masses. Over time, repetitive observations of the exoplanets at medium resolution ($R=500$) will provide a catalogue of K-band spectra of unprecedented quality, for a number of exoplanets. The K-band has the unique properties that it contains many molecular signatures (CO, H$_2$O, CH$_4$, CO$_2$). This allows constraining precisely surface gravity, metallicity, and temperature, if used in conjunction with self-consistent models like Exo-REM. Further, we will use the parameter-retrieval algorithm petitRADTRANS to constrain the C/O ratio of the planets. Ultimately, we plan to produce the first C/O survey of exoplanets, kick-starting the difficult process of linking planetary formation with measured atomic abundances.

preprint2020arXiv

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

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.

preprint2020arXiv

Detection of the Schwarzschild precession in the orbit of the star S2 near the Galactic centre massive black hole

The star S2 orbiting the compact radio source Sgr A* is a precision probe of the gravitational field around the closest massive black hole (candidate). Over the last 2.7 decades we have monitored the star's radial velocity and motion on the sky, mainly with the SINFONI and NACO adaptive optics (AO) instruments on the ESO VLT, and since 2017, with the four-telescope interferometric beam combiner instrument GRAVITY. In this paper we report the first detection of the General Relativity (GR) Schwarzschild Precession (SP) in S2's orbit. Owing to its highly elliptical orbit (e = 0.88), S2's SP is mainly a kink between the pre-and post-pericentre directions of motion ~ +- 1 year around pericentre passage, relative to the corresponding Kepler orbit. The superb 2017-2019 astrometry of GRAVITY defines the pericentre passage and outgoing direction. The incoming direction is anchored by 118 NACO-AO measurements of S2's position in the infrared reference frame, with an additional 75 direct measurements of the S2-Sgr A* separation during bright states ('flares') of Sgr A*. Our 14-parameter model fits for the distance, central mass, the position and motion of the reference frame of the AO astrometry relative to the mass, the six parameters of the orbit, as well as a dimensionless parameter f_SP for the SP (f_SP = 0 for Newton and 1 for GR). From data up to the end of 2019 we robustly detect the SP of S2, del phi = 12' per orbital period. From posterior fitting and MCMC Bayesian analysis with different weighting schemes and bootstrapping we find f_SP = 1.10 +- 0.19. The S2 data are fully consistent with GR. Any extended mass inside S2's orbit cannot exceed ~ 0.1% of the central mass. Any compact third mass inside the central arcsecond must be less than about 1000 M_sun.

preprint2020arXiv

Dynamically important magnetic fields near the event horizon of Sgr A*

We study the time-variable linear polarisation of Sgr A* during a bright NIR flare observed with the GRAVITY instrument on July 28, 2018. Motivated by the time evolution of both the observed astrometric and polarimetric signatures, we interpret the data in terms of the polarised emission of a compact region ('hotspot') orbiting a black hole in a fixed, background magnetic field geometry. We calculated a grid of general relativistic ray-tracing models, created mock observations by simulating the instrumental response, and compared predicted polarimetric quantities directly to the measurements. We take into account an improved instrument calibration that now includes the instrument's response as a function of time, and we explore a variety of idealised magnetic field configurations. We find that the linear polarisation angle rotates during the flare, which is consistent with previous results. The hotspot model can explain the observed evolution of the linear polarisation. In order to match the astrometric period of this flare, the near horizon magnetic field is required to have a significant poloidal component, which is associated with strong and dynamically important fields. The observed linear polarisation fraction of $\simeq 30\%$ is smaller than the one predicted by our model ($\simeq 50\%$). The emission is likely beam depolarised, indicating that the flaring emission region resolves the magnetic field structure close to the black hole.

preprint2020arXiv

Optical interferometry and Gaia measurement uncertainties reveal the physics of asymptotic giant branch stars

Context. Asymptotic giant branch stars are cool luminous evolved stars that are well observable across the Galaxy and populating Gaia data. They have complex stellar surface dynamics Aims. On the AGB star CL Lac, it has been shown that the convection-related variability accounts for a substantial part of the Gaia DR2 parallax error. We observed this star with the MIRC-X beam combiner installed at the CHARA interferometer to detect the presence of stellar surface inhomogeneities. Methods. We performed the reconstruction of aperture synthesis images from the interferometric observations at different wavelengths. Then, we used 3D radiative hydrodynamics simulations of stellar convection with CO5BOLD and the post-processing radiative transfer code Optim3D to compute intensity maps in the spectral channels of MIRC-X observations. Then, we determined the stellar radius and compared the 3D synthetic maps to the reconstructed ones focusing on matching the intensity contrast, the morphology of stellar surface structures, and the photocentre position at two different spectral channels, 1.52 and 1.70 micron, simultaneously. Results. We measured the apparent diameter of CL Lac at two wavelengths and recovered the radius using a Gaia parallax. In addition to this, the reconstructed images are characterised by the presence of a brighter area that largely affects the position of the photocentre. The comparison with 3D simulation shows good agreement with the observations both in terms of contrast and surface structure morphology, meaning that our model is adequate for explaining the observed inhomogenities. Conclusions. This work confirms the presence of convection-related surface structures on an AGB star of Gaia DR2. Our result will help us to take a step forward in exploiting Gaia measurement uncertainties to extract the fundamental properties of AGB stars using appropriate RHD simulations.

preprint2020arXiv

Probing the magnetospheric accretion region of the young pre-transitional disk system DoAr 44 using VLTI/GRAVITY

Young stellar objects are thought to accrete material from their circumstellar disks through their strong stellar magnetospheres. We aim to directly probe the magnetospheric accretion region on a scale of a few 0.01 au in a young stellar system using long-baseline optical interferometry. We observed the pre-transitional disk system DoAr 44 with VLTI/GRAVITY on two consecutive nights in the K-band. We computed interferometric visibilities and phases in the continuum and in the BrG line in order to constrain the extent and geometry of the emitting regions. We resolve the continuum emission of the inner dusty disk and measure a half-flux radius of 0.14 au. We derive the inclination and position angle of the inner disk, which provides direct evidence that the inner and outer disks are misaligned in this pre-transitional system. This may account for the shadows previously detected in the outer disk. We show that BrG emission arises from an even more compact region than the inner disk, with an upper limit of 0.047 au (5 Rstar). Differential phase measurements between the BrG line and the continuum allow us to measure the astrometric displacement of the BrG line-emitting region relative to the continuum on a scale of a few tens of microarcsec, corresponding to a fraction of the stellar radius. Our results can be accounted for by a simple geometric model where the BrG line emission arises from a compact region interior to the inner disk edge, on a scale of a few stellar radii, fully consistent with the concept of magnetospheric accretion process in low-mass young stellar systems.

preprint2020arXiv

The Flux Distribution of Sgr A*

The Galactic Center black hole Sagittarius A* is a variable NIR source that exhibits bright flux excursions called flares. The low-flux density turnover of the flux distribution is below the sensitivity of current single-aperture telescopes. We use the unprecedented resolution of the GRAVITY instrument at the VLTI. Our light curves are unconfused, overcoming the confusion limit of previous photometric studies. We analyze the light curves using standard statistical methods and obtain the flux distribution. We find that the flux distribution of SgrA* turns over at a median flux density of (1.1\pm0.3)mJy. We measure the percentiles of the flux distribution and use them to constrain the NIR K-band SED. Furthermore, we find that the flux distribution is intrinsically right-skewed to higher flux density in log space. Flux densities below 0.1mJy are hardly ever observed. In consequence, a single powerlaw or lognormal distribution does not suffice to describe the observed flux distribution in its entirety. However, if one takes into account a power law component at high flux densities, a lognormal distribution can describe the lower end of the observed flux distribution. We confirm the RMS-flux relation for Sgr~A* and find it to be linear for all flux densities in our observation. We conclude that Sgr~A* has two states: the bulk of the emission is generated in a lognormal process with a well-defined median flux density and this quiescent emission is supplemented by sporadic flares that create the observed power law extension of the flux distribution.

preprint2020arXiv

The GRAVITY Young Stellar Object survey IV. The CO overtone emission in 51 Oph at sub-au scales

51 Oph is a Herbig Ae/Be star that exhibits strong near-infrared CO ro-vibrational emission at 2.3 micron, most likely originating in the innermost regions of a circumstellar disc. We aim to obtain the physical and geometrical properties of the system by spatially resolving the circumstellar environment of the inner gaseous disc. We used the second-generation VLTI/GRAVITY to spatially resolve the continuum and the CO overtone emission. We obtained data over 12 baselines with the auxiliary telescopes and derive visibilities, and the differential and closure phases as a function of wavelength. We used a simple LTE ring model of the CO emission to reproduce the spectrum and CO line displacements. Our interferometric data show that the star is marginally resolved at our spatial resolution, with a radius of 10.58+-2.65 Rsun.The K-band continuum emission from the disc is inclined by 63+-1 deg, with a position angle of 116+-1 deg, and 4+-0.8 mas (0.5+-0.1 au) across. The visibilities increase within the CO line emission, indicating that the CO is emitted within the dust-sublimation radius.By modelling the CO bandhead spectrum, we derive that the CO is emitted from a hot (T=1900-2800 K) and dense (NCO=(0.9-9)x10^21 cm^-2) gas. The analysis of the CO line displacement with respect to the continuum allows us to infer that the CO is emitted from a region 0.10+-0.02 au across, well within the dust-sublimation radius. The inclination and position angle of the CO line emitting region is consistent with that of the dusty disc. Our spatially resolved interferometric observations confirm the CO ro-vibrational emission within the dust-free region of the inner disc. Conventional disc models exclude the presence of CO in the dust-depleted regions of Herbig AeBe stars. Ad hoc models of the innermost disc regions, that can compute the properties of the dust-free inner disc, are therefore required.

preprint2020arXiv

The GRAVITY young stellar object survey. II. First spatially resolved observations of the CO bandhead emission in a high-mass YSO

The inner regions of the discs of high-mass young stellar objects (HMYSOs) are still poorly known due to the small angular scales and the high visual extinction involved. We deploy near-infrared (NIR) spectro-interferometry to probe the inner gaseous disc in HMYSOs and investigate the origin and physical characteristics of the CO bandhead emission (2.3-2.4 $μ$m). We present the first GRAVITY/VLTI observations at high spectral (R=4000) and spatial (mas) resolution of the CO overtone transitions in NGC 2024 IRS2. The continuum emission is resolved in all baselines and is slightly asymmetric, displaying small closure phases ($\leq$8$^{\circ}$). Our best ellipsoid model provides a disc inclination of 34$^{\circ}$$\pm$1$^{\circ}$, a disc major axis position angle of 166$^{\circ}$$\pm$1$^{\circ}$, and a disc diameter of 3.99$\pm$0.09 mas (or 1.69$\pm$0.04 au, at a distance of 423 pc). The small closure phase signals in the continuum are modelled with a skewed rim, originating from a pure inclination effect. For the first time, our observations spatially and spectrally resolve the first four CO bandheads. Changes in visibility, as well as differential and closure phases across the bandheads are detected. Both the size and geometry of the CO-emitting region are determined by fitting a bidimensional Gaussian to the continuum-compensated CO bandhead visibilities. The CO-emitting region has a diameter of 2.74$\pm^{0.08}_{0.07}$ mas (1.16$\pm$0.03 au), and is located in the inner gaseous disc, well within the dusty rim, with inclination and $PA$ matching the dusty disc geometry, which indicates that both dusty and gaseous discs are coplanar. Physical and dynamical gas conditions are inferred by modelling the CO spectrum. Finally, we derive a direct measurement of the stellar mass of $M_*\sim$14.7$^{+2}_{-3.6}$ M$_{\odot}$ by combining our interferometric and spectral modelling results.

preprint2019arXiv

Peering into the formation history of beta Pictoris b with VLTI/GRAVITY long baseline interferometry

Our objective is to estimate the C/O ratio in the atmosphere of beta Pictoris b and obtain an estimate of the dynamical mass of the planet, as well as to refine its orbital parameters using high-precision astrometry. We used the GRAVITY instrument with the four 8.2 m telescopes of the Very Large Telescope Interferometer to obtain K-band spectro-interferometric data on $β$ Pic b. We extracted a medium resolution (R=500) K-band spectrum of the planet and a high-precision astrometric position. We estimated the planetary C/O ratio using two different approaches (forward modeling and free retrieval) from two different codes (ExoREM and petitRADTRANS, respectively). Finally, we used a simplified model of two formation scenarios (gravitational collapse and core-accretion) to determine which can best explain the measured C/O ratio. Our new astrometry disfavors a circular orbit for $β$ Pic b ($e=0.15^{+0.05}_{-0.04}$). Combined with previous results and with Hipparcos/GAIA measurements, this astrometry points to a planet mass of $M = 12.7\pm{}2.2\,M_\mathrm{Jup}$. This value is compatible with the mass derived with the free-retrieval code petitRADTRANS using spectral data only. The forward modeling and free-retrieval approches yield very similar results regarding the atmosphere of beta Pic b. In particular, the C/O ratios derived with the two codes are identical ($0.43\pm{}0.05$ vs $0.43^{+0.04}_{-0.03}$). We argue that if the stellar C/O in $β$ Pic is Solar, then this combination of a very high mass and a low C/O ratio for the planet suggests a formation through core-accretion, with strong planetesimal enrichment.

preprint2016arXiv

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

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.

preprint2016arXiv

Imaging the dust sublimation front of a circumbinary disk

We present the first near-IR milli-arcsecond-scale image of a post-AGB binary that is surrounded by hot circumbinary dust. A very rich interferometric data set in six spectral channels was acquired of IRAS08544-4431 with the new RAPID camera on the PIONIER beam combiner at the Very Large Telescope Interferometer (VLTI). A broadband image in the \textit{H} band was reconstructed by combining the data of all spectral channels using the SPARCO method. We spatially separate all the building blocks of the IRAS08544-4431 system in our milliarcsecond-resolution image. Our dissection reveals a dust sublimation front that is strikingly similar to that expected in early-stage protoplanetary disks, as well as an unexpected flux signal of $\sim$4\% from the secondary star. The energy output from this companion indicates the presence of a compact circum-companion accretion disk, which is likely the origin of the fast outflow detected in H$α$. Our image provides the most detailed view into the heart of a dusty circumstellar disk to date. Our results demonstrate that binary evolution processes and circumstellar disk evolution can be studied in detail in space and over time.

preprint2016arXiv

The red dwarf pair GJ65AB: inflated, spinning twins of Proxima

The nearby red dwarf binary GJ65AB (UB+BL Ceti, M5.5Ve+M6Ve) is a cornerstone system to probe the physics of very low mass stars. We present new interferometric measurements of the angular diameters of the two components of GJ65AB with the PIONIER instrument in the H band: thetaUD (A) = 0.558 +/- 0.008 +/- 0.020 mas and thetaUD (B) = 0.539 +/- 0.009 +/- 0.020 mas. They translate into limb-darkened angular diameters of LD (A) = 0.573 +/- 0.021 mas and LD (B) = 0.554 +/- 0.022 mas. From the known parallax, the radii are R(A) = 0.165 +/- 0.006 Rsun and R(B) = 0.159 +/- 0.006 Rsun (sigma(R)/R = 4%). We also observed GJ65AB with the VLT/NACO adaptive optics and refine the orbital parameters and infrared magnitudes of the system. We derive masses for the two components of m(A) = 0.1225 +/- 0.0043 Msun and m(B) = 0.1195 +/- 0.0043 Msun (sigma(m)/m = 4%). To derive the radial and rotational velocities of the two stars, as well as their relative metallicity with respect to Proxima, we also present new individual UVES high resolution spectra. Placing GJ65 A&B in the mass-radius diagram shows that their radii are larger respectively by 14 +/- 4% and 12 +/- 4% compared to expectations from the models by Baraffe et al. (2015). Following Chabrier et al. (2007), we propose that this discrepancy is caused by the inhibition of convective energy transport by a strong internal magnetic field generated by dynamo effect in these two fast rotating stars. A comparison with the almost identical twin but slowly-rotating Proxima strengthens this hypothesis, as the radius of Proxima does not appear inflated compared to models.

preprint2015arXiv

Inner disk clearing around the Herbig Ae star HD\,139614: Evidence for a planet-induced gap ?

Spatially resolving the inner dust cavity of the transitional disks is a key to understanding the connection between planetary formation and disk dispersal. The disk around the Herbig star HD 139614 is of particular interest since it presents a pretransitional nature with an au-sized gap, in the dust, that was spatially resolved by mid-IR interferometry. Using new NIR interferometric observations, we aim to characterize the 0.1-10~au region of the HD~139614 disk further and identify viable mechanisms for the inner disk clearing. We report the first multiwavelength radiative transfer modeling of the interferometric data acquired on HD~139614 with PIONIER, AMBER, and MIDI, complemented by Herschel/PACS photometries. We confirm a gap structure in the um-sized dust, extending from about 2.5 au to 6 au, and constrained the properties of the inner dust component: e.g., a radially increasing surface density profile, and a depletion of 10^3 relative to the outer disk. Since self-shadowing and photoevaporation appears unlikely to be responsible for the au-sized gap of HD~139614, we thus tested if dynamical clearing could be a viable mechanism using hydrodynamical simulations to predict the gaseous disk structure. Indeed, a narrow au-sized gap is expected when a single giant planet interacts with the disk. Assuming that small dust grains are well coupled to the gas, we found that a ~ 3~Mjup planet located at 4.5 au from the star could, in less than 1 Myr, reproduce most of the aspects of the dust surface density profile, while no significant depletion in gas occurred in the inner disk, in contrast to the dust. However, the dust-depleted inner disk could be explained by the expected dust filtration by the gap and the efficient dust growth/fragmentation in the inner disk regions. Our results support the hypothesis of a giant planet opening a gap and shaping the inner region of the HD~139614 disk.

preprint2015arXiv

Masses of the components of SB2 binaries observed with Gaia. II. Masses derived from PIONIER interferometric observations for Gaia validation

In anticipation of the Gaia astrometric mission, a sample of spectroscopic binaries is being observed since 2010 with the Sophie spectrograph at the Haute--Provence Observatory. Our aim is to derive the orbital elements of double-lined spectroscopic binaries (SB2s) with an accuracy sufficient to finally obtain the masses of the components with relative errors as small as 1 % when combined with Gaia astrometric measurements. In order to validate the masses derived from Gaia, interferometric observations are obtained for three SB2s in our sample with F-K components: HIP 14157, HIP 20601 and HIP 117186. The masses of the six stellar components are derived. Due to its edge-on orientation, HIP 14157 is probably an eclipsing binary. We note that almost all the derived masses are a few percent larger than the expectations from the standard spectral-type-mass calibration and mass-luminosity relation. Our calculation also leads to accurate parallaxes for the three binaries, and the Hipparcos parallaxes are confirmed.

preprint2014arXiv

Binarity of the LBV HR Car

VLTI/AMBER and VLTI/PIONIER observations of the LBV HR Car show an interferometric signature that could not possibly be explained by an extended wind, more or less symmetrically distributed around a single object. Instead, observations both in the Br$γ$ line and the H-band continuum are best explained by two point sources (or alternatively one point source and one slightly extended source) at about 2 mas separation and a contrast ratio of about 1:5. These observations establish that HR Car is a binary, but further interpretation will only be possible with future observations to constrain the orbit. Under the assumption that the current separation is close to the maximum one, the orbital period can be estimated to be of the order of 5 years, similar as in the $η$ Car system. This would make HR Car the second such LBV binary.

preprint2014arXiv

GG Tau: the fifth element

We aim at unveiling the observational imprint of physical mechanisms that govern planetary formation in young, multiple systems. In particular, we investigate the impact of tidal truncation on the inner circumstellar disks. We observed the emblematic system GG Tau at high-angular resolution: a hierarchical quadruple system composed of low-mass T Tauri binary stars surrounded by a well-studied, massive circumbinary disk in Keplerian rotation. We used the near-IR 4-telescope combiner PIONIER on the VLTI and sparse-aperture-masking techniques on VLT/NaCo to probe this proto-planetary system at sub-AU scales. We report the discovery of a significant closure-phase signal in H and Ks bands that can be reproduced with an additional low-mass companion orbiting GG Tau Ab, at a (projected) separation rho = 31.7 +/- 0.2mas (4.4 au) and PA = 219.6 +/- 0.3deg. This finding offers a simple explanation for several key questions in this system, including the missing-stellar-mass problem and the asymmetry of continuum emission from the inner dust disks observed at millimeter wavelengths. Composed of now five co-eval stars with 0.02 <= Mstar <= 0.7 Msun, the quintuple system GG Tau has become an ideal test case to constrain stellar evolution models at young ages (few 10^6yr).

preprint2014arXiv

HD 152246 - a new high-mass triple system and its basic properties

Analyses of multi-epoch, high-resolution (R ~ 50.000) optical spectra of the O-type star HD 152246 (O9 IV according to the most recent classification), complemented by a limited number of earlier published radial velocities, led to the finding that the object is a hierarchical triple system, where a close inner pair (Ba-Bb) with a slightly eccentric orbit (e = 0.11) and a period of 6.0049 days revolves in a 470-day highly eccentric orbit (e = 0.865) with another massive and brighter component A. The mass ratio of the inner system must be low since we were unable to find any traces of the secondary spectrum. The mass ratio A/(Ba+Bb) is 0.89. The outer system has recently been resolved using long-baseline interferometry on three occasions. The interferometry confirms the spectroscopic results and specifies elements of the system. Our orbital solutions, including the combined radial-velocity and interferometric solution indicate an orbital inclination of the outer orbit of 112° and stellar masses of 20.4 and 22.8 solar masses. We also disentangled the spectra of components A and Ba and compare them to synthetic spectra from two independent programmes, TLUSTY and FASTWIND. In either case, the fit was not satisfactory and we postpone a better determination of the system properties for a future study, after obtaining observations during the periastron passage of the outer orbit (the nearest chance being March 2015). For the moment, we can only conclude that component A is an O9 IV star with v*sin(i) = 210 +\- 10 km/s and effective temperature of 33000 +\- 500 K, while component Ba is an O9 V object with v*sin(i) = 65 +/- 3 km/s and T_eff = 33600 +\- 600 K.

preprint2014arXiv

On the structure of the transition disk around TW Hya

For over a decade, the structure of the inner cavity in the transition disk of TW Hydrae has been a subject of debate. Modeling the disk with data obtained at different wavelengths has led to a variety of proposed disk structures. Rather than being inconsistent, the individual models might point to the different faces of physical processes going on in disks, such as dust growth and planet formation. Our aim is to investigate the structure of the transition disk again and to find to what extent we can reconcile apparent model differences. A large set of high-angular-resolution data was collected from near-infrared to centimeter wavelengths. We investigated the existing disk models and established a new self-consistent radiative-transfer model. A genetic fitting algorithm was used to automatize the parameter fitting. Simple disk models with a vertical inner rim and a radially homogeneous dust composition from small to large grains cannot reproduce the combined data set. Two modifications are applied to this simple disk model: (1) the inner rim is smoothed by exponentially decreasing the surface density in the inner ~3 AU, and (2) the largest grains (>100 um) are concentrated towards the inner disk region. Both properties can be linked to fundamental processes that determine the evolution of protoplanetary disks: the shaping by a possible companion and the different regimes of dust-grain growth, respectively. The full interferometric data set from near-infrared to centimeter wavelengths requires a revision of existing models for the TW Hya disk. We present a new model that incorporates the characteristic structures of previous models but deviates in two key aspects: it does not have a sharp edge at 4 AU, and the surface density of large grains differs from that of smaller grains. This is the first successful radiative-transfer-based model for a full set of interferometric data.

preprint2014arXiv

Southern Massive Stars at High Angular Resolution: Observational Campaign and Companion Detection

Multiplicity is one of the most fundamental observable properties of massive O-type stars and offers a promising way to discriminate between massive star formation theories. Nevertheless, companions at separations between 1 and 100 mas remain mostly unknown due to intrinsic observational limitations. [...] The Southern MAssive Stars at High angular resolution survey (SMASH+) was designed to fill this gap by providing the first systematic interferometric survey of Galactic massive stars. We observed 117 O-type stars with VLTI/PIONIER and 162 O-type stars with NACO/SAM, respectively probing the separation ranges 1-45 and 30-250mas and brightness contrasts of Delta H < 4 and Delta H < 5. Taking advantage of NACO's field-of-view, we further uniformly searched for visual companions in an 8''-radius down to Delta H = 8. This paper describes the observations and data analysis, reports the discovery of almost 200 new companions in the separation range from 1mas to 8'' and presents the catalog of detections, including the first resolved measurements of over a dozen known long-period spectroscopic binaries. Excluding known runaway stars for which no companions are detected, 96 objects in our main sample (DEC < 0 deg; H<7.5) were observed both with PIONIER and NACO/SAM. The fraction of these stars with at least one resolved companion within 200mas is 0.53. Accounting for known but unresolved spectroscopic or eclipsing companions, the multiplicity fraction at separation < 8'' increases to f_m = 0.91 +/- 0.03. The fraction of luminosity class V stars that have a bound companion reaches 100% at 30mas while their average number of physically connected companions within 8'' is f_c = 2.2 +/- 0.3. This demonstrates that massive stars form nearly exclusively in multiple systems. Additionally, the nine non-thermal (NT) radio emitters observed by SMASH+ are all resolved [...]

preprint2014arXiv

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

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.

preprint2014arXiv

The BinaMIcS project: understanding the origin of magnetic fields in massive stars through close binary systems

It is now well established that a fraction of the massive (M>8 Msun) star population hosts strong, organised magnetic fields, most likely of fossil origin. The details of the generation and evolution of these fields are still poorly understood. The BinaMIcS project takes an important step towards the understanding of the interplay between binarity and magnetism during the stellar formation and evolution, and in particular the genesis of fossil fields, by studying the magnetic properties of close binary systems. The components of such systems are most likely formed together, at the same time and in the same environment, and can therefore help us to disentangle the role of initial conditions on the magnetic properties of the massive stars from other competing effects such as age or rotation. We present here the main scientific objectives of the BinaMIcS project, as well as preliminary results from the first year of observations from the associated ESPaDOnS and Narval spectropolarimetric surveys.

preprint2014arXiv

The Three body system delta Circini

Delta Cir is known as an O7.5 III eclipsing and spectroscopic binary with an eccentric orbit. Penny et al. discovered the presence of a third component in the IUE spectra. The eclipsing binary and the third body revolve around a common centre of gravity with a period of 1644 days in an eccentric orbit with the semimajor axis of 10 AU. We demonstrate the presence of the apsidal-line rotation with a period of about 141 years, which is considerably longer than its theoretically predicted value, based on the published radii of the binary components derived from the Hipparchos Hp light curve. However, our new solution of the same light curve resulted in smaller radii and a better agreement between the observed and predicted period of the apsidal- line advance. There are indications that the third body is a binary. The object was resolved by VLTI with the PIONIER combiner; in June 2012 the separation was 3.78 mas, magnitude difference in the H region 1.75 mag. This result means that (assuming the distance 770 pc) the inclination of the long orbit is 87.7 degrees.

preprint2013arXiv

Multiplicity of Galactic Cepheids from long-baseline interferometry. II. The Companion of AX Circini revealed with VLTI/PIONIER

Aims: We aim at detecting and characterizing the main-sequence companion of the Cepheid AX Cir ($P_\mathrm{orb} \sim $ 18 yrs). The long-term objective is to estimate the mass of both components and the distance to the system. Methods: We used the PIONIER combiner at the VLT Interferometer to obtain the first interferometric measurements of the short-period Cepheid AX Cir and its orbiting component. Results: The companion is resolved by PIONIER at a projected separation $ρ= 29.2 \pm 0.2$ mas and projection angle $PA = 167.6 \pm 0.3^{\circ}$. We measured $H$-band flux ratios between the companion and the Cepheid of $0.90 \pm 0.10$ % and $0.75 \pm 0.17$ %, respectively at a pulsation phase for the Cepheid $ϕ= 0.24$ and 0.48. The lower contrast at $ϕ= 0.48$ is due to increased brightness of the Cepheid compared to the $ϕ= 0.24$. This gives an average apparent magnitude $m\mathrm{_H (comp)} = 9.06 \pm 0.24$ mag. The limb-darkened angular diameter of the Cepheid at the two pulsation phases was measured to be $θ_\mathrm{LD} = 0.839 \pm 0.023$ mas and $θ_\mathrm{LD} = 0.742 \pm 0.020$ mas, respectively at $ϕ= 0.24$ and 0.48. A lower limit on the total mass of the system was also derived based on our measured separation, we found $M_\mathrm{T} \geq 9.7 \pm 0.6 M_\odot$.

preprint2013arXiv

The close environment of high-mass X-ray binaries at high angular resolution I. VLTI/AMBER and VLTI/PIONIER near-infrared interferometric observations of Vela X-1

Recent improvements on the sensitivity and spectral resolution of X-ray observations have led to a better understanding of the properties of matter in the vicinity of High Mass X-ray Binaries hosting a supergiant star and a compact object. However the geometry and physical properties of their environment at larger scales are currently only predicted by simulations. We aim at exploring the environment of Vela X-1 at a few stellar radii of the supergiant using spatially resolved observations in the near-infrared and at studying its dynamical evolution along the 9-day orbital period of the system. We observed Vela X-1 in 2010 and 2012 using long baseline interferometry at VLTI, respectively with the AMBER instrument in the K band and the PIONIER instrument in the H band. The PIONIER observations span through one orbital period to monitor possible evolutions in the geometry of the system. We resolved a structure of $8\pm3~R_\star$ from the AMBER data and $2.0\,_{-1.2}^{+0.7}~R_\star$ from the PIONIER data. From the closure phase we found that the environment of Vela X-1 is symmetrical. We observed comparable measurements between the continuum and the spectral lines in the K band, meaning that both emissions originate from the same forming region. From the monitoring of the system over one period in 2012, we found the signal to be constant with the orbital phase within the error bars. We propose three scenarios for the discrepancy between the two measurements: either there is a strong temperature gradient in the supergiant wind leading to a hot component much more compact than the cool part of the wind observed in the K band, or we observed a diffuse shell in 2010 possibly triggered by an off-state in the accretion rate of the pulsar that was dissolved in the interstellar medium in 2012, or the structure observed in the H band was the stellar photosphere instead of the supergiant wind.

preprint2013arXiv

Three-dimensional orbits of the triple-O stellar system HD 150136

Context. HD 150136 is a triple hierarchical system and a non-thermal radio emitter. It is formed by an O3-3.5 V + O5.5-6 V close binary and a more distant O6.5-7 V tertiary. So far, only the inner orbital properties have been reliably constrained. Aims. To quantitatively understand the non-thermal emission process, accurate knowledge of the physical and orbital properties of the object is crucial. Here, we aim to investigate the orbital properties of the wide system and to constrain the inclinations of the inner and outer binaries, and with these the absolute masses of the system components. Methods. We used the PIONIER combiner at the Very Large Telescope Interferometer to obtain the very first interferometric measurements of HD 150136. We combined the interferometric observations with new and existing high resolution spectroscopic data to derive the orbital solution of the outer companion in the three-dimensional space. Results. The wide system is clearly resolved by PIONIER, with a projected separation on the plane of the sky of about 9 milli-arcsec. The best-fit orbital period, eccentricity, and inclination are 8.2 yr, 0.73 and 108 degr. We constrain the masses of the three stars of the system to 63 +/- 10, 40 +/- 6, and 33 +/- 12 Msun for the O3-3.5 V, O5.5-6 V and O6.5-7 V components. Conclusions. The dynamical masses agree within errors with the evolutionary masses of the components. Future interferometric and spectroscopic monitoring of HD 150136 should allow one to reduce the uncertainties to a few per cent only and to accurately constrain the distance to the system. This makes HD 150136 an ideal system to quantitatively test evolutionary models of high-mass stars as well as the physics of non-thermal processes occurring in O-type systems.

preprint2012arXiv

Binary interaction dominates the evolution of massive stars

The presence of a nearby companion alters the evolution of massive stars in binary systems, leading to phenomena such as stellar mergers, X-ray binaries and gamma-ray bursts. Unambiguous constraints on the fraction of massive stars affected by binary interaction were lacking. We simultaneously measured all relevant binary characteristics in a sample of Galactic massive O stars and quantified the frequency and nature of binary interactions. Over seventy per cent of all massive stars will exchange mass with a companion, leading to a binary merger in one third of the cases. These numbers greatly exceed previous estimates and imply that binary interaction dominates the evolution of massive stars, with implications for populations of massive stars and their supernovae.

preprint2012arXiv

Hot circumstellar material resolved around beta Pic with VLTI/PIONIER

We aim at resolving the circumstellar environment around beta Pic in the near-infrared in order to study the inner planetary system (< 200 mas, i.e., ~4 AU). Precise interferometric fringe visibility measurements were obtained over seven spectral channels dispersed across the H band with the four-telescope VLTI/PIONIER interferometer. Thorough analysis of interferometric data was performed to measure the stellar angular diameter and to search for circumstellar material. We detected near-infrared circumstellar emission around beta Pic that accounts for 1.37% +/- 0.16% of the near-infrared stellar flux and that is located within the field-of-view of PIONIER (i.e., ~200 mas in radius). The flux ratio between this excess and the photosphere emission is shown to be stable over a period of 1 year and to vary only weakly across the H band, suggesting that the source is either very hot (> 1500 K) or dominated by the scattering of the stellar flux. In addition, we derived the limb-darkened angular diameter of beta Pic with an unprecedented accuracy (theta_LD= 0.736 +/- 0.019 mas). The presence of a small H-band excess originating in the vicinity of beta Pic is revealed for the first time thanks to the high-precision visibilities enabled by VLTI/PIONIER. This excess emission is likely due to the scattering of stellar light by circumstellar dust and/or the thermal emission from a yet unknown population of hot dust, although hot gas emitting in the continuum cannot be firmly excluded.

preprint2012arXiv

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

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.

preprint2012arXiv

Multiplicity of massive O stars and evolutionary implications

Nearby companions alter the evolution of massive stars in binary systems. Using a sample of Galactic massive stars in nearby young clusters, we simultaneously measure all intrinsic binary characteristics relevant to quantify the frequency and nature of binary interactions. We find a large intrinsic binary fraction, a strong preference for short orbital periods and a flat distribution for the mass-ratios. Our results do not support the presence of a significant peak of equal-mass `twin' binaries. As a result of the measured distributions, we find that over seventy per cent of all massive stars exchange mass with a companion. Such a rate greatly exceeds previous estimates and implies that the majority of massive stars have their evolution strongly affected by interaction with a nearby companion.

preprint2012arXiv

New activity in the large circumstellar disk of the Be-shell star 48 Lib

Spectroscopic, polarimetric, and high spectral resolution interferometric data covering the period 1995-2011 are analyzed to document the transition into a new phase of circumstellar disk activity in the classical Be-shell star 48 Lib. The long-term disk evolution is described using the V/R ratio of the violet and red emission components of H alpha and Br gamma, radial velocities and profiles of He I and optical metal shell lines, multi-band BVRI polarimetry, broad-band and high-resolution interferometric visibilities and phases. Spectroscopic signatures of disk asymmetries in 48 Lib vanished in the late nineties but recovered some time between 2004 and 2007. Variations in the radial velocity and line profile of conventional shell lines correlate with the V/R behavior. They are shared by narrow absorption cores superimposed on otherwise seemingly photospheric He I lines, which may form in high-density gas at the inner disk close to the photosphere. Large radial velocity variations continued also during the V/R-quiet years. The broad H-band interferometry gives a disk diameter of (1.72+-0.2) mas (equivalent to 15 stellar radii), position angle of the disk (50+-9) deg and a relatively low disk flattening of 1.66+-0.3. Within the errors the same disk position angle is derived from polarimetric observations and from photocenter shifts across Br gamma. The high-resolution interferometric visibility and phase profiles show a double or even multiple-component structure. A preliminary estimate based on the size of the Br gamma emitting region indicates a large diameter for the disk (tens of stellar radii).

preprint2012arXiv

PIONIER: a status report

The visitor instrument PIONIER provides VLTI with improved imaging capabilities and sensitivity. The instrument started routinely delivering scientific data in November 2010, that is less than 12 months after being approved by the ESO Science and Technical Committee. We recall the challenges that had to be tackled to design, built and commission PIONIER. We summarize the typical performances and some astrophysical results obtained so far. We conclude this paper by summarizing lessons learned.

preprint2012arXiv

The long period eccentric orbit of the particle accelerator HD167971 revealed by long baseline interferometry

Using optical long baseline interferometry, we resolved for the first time the two wide components of HD167971, a candidate hierarchical triple system known to efficiently accelerate particles. Our multi-epoch VLTI observations provide direct evidence for a gravitational link between the O8 supergiant and the close eclipsing O + O binary. The separation varies from 8 to 15 mas over the three-year baseline of our observations, suggesting that the components evolve on a wide and very eccentric orbit (most probably e>0.5). These results provide evidence that the wide orbit revealed by our study is not coplanar with the orbit of the inner eclipsing binary. From our measurements of the near-infrared luminosity ratio, we constrain the spectral classification of the components in the close binary to be O6-O7, and confirm that these stars are likely main-sequence objects. Our results are discussed in the context of the bright non-thermal radio emission already reported for this system, and we provide arguments in favour of a maximum radio emission coincident with periastron passage. HD167971 turns out to be an efficient O-type particle accelerator that constitutes a valuable target for future high angular resolution radio imaging using VLBI facilities.

preprint2011arXiv

An incisive look at the symbiotic star SS Leporis -- Milli-arcsecond imaging with PIONIER/VLTI

Context. Determining the mass transfer in a close binary system is of prime importance for understanding its evolution. SS Leporis, a symbiotic star showing the Algol paradox and presenting clear evidence of ongoing mass transfer, in which the donor has been thought to fill its Roche lobe, is a target particularly suited to this kind of study. Aims. Since previous spectroscopic and interferometric observations have not been able to fully constrain the system morphology and characteristics, we go one step further to determine its orbital parameters, for which we need new interferometric observations directly probing the inner parts of the system with a much higher number of spatial frequencies. Methods. We use data obtained at eight different epochs with the VLTI instruments AMBER and PIONIER in the H- and K-bands. We performed aperture synthesis imaging to obtain the first model-independent view of this system. We then modelled it as a binary (whose giant is spatially resolved) that is surrounded by a circumbinary disc. Results. Combining these interferometric measurements with previous radial velocities, we fully constrain the orbit of the system. We then determine the mass of each star and significantly revise the mass ratio. The M giant also appears to be almost twice smaller than previously thought. Additionally, the low spectral resolution of the data allows the flux of both stars and of the dusty disc to be determined along the H and K bands, and thereby extracting their temperatures. Conclusions. We find that the M giant actually does not stricto sensus fill its Roche lobe. The mass transfer is more likely to occur through the accretion of an important part of the giant wind. We finally rise the possibility for an enhanced mass loss from the giant, and we show that an accretion disc should have formed around the A star.

preprint2011arXiv

Hot exozodiacal dust resolved around Vega with IOTA/IONIC

Although debris discs have been detected around a significant number of main-sequence stars, only a few of them are known to harbour hot dust in their inner part where terrestrial planets may have formed. Thanks to infrared interferometric observations, it is possible to obtain a direct measurement of these regions, which are of prime importance for preparing future exo-Earth characterisation missions. In this context, we have resolved the exozodiacal dust disc around Vega with the help of infrared stellar interferometry and estimated the integrated H-band flux originating from the first few AUs of the debris disc. Using precise H-band interferometric measurements obtained with the 3-telescope IOTA/IONIC interferometer (Mount Hopkins, Arizona), thorough modelling of both interferometric data (squared visibility and closure phase) and spectral energy distribution was performed to constrain the nature of the near-infrared excess emission. The most straightforward scenario consists in a compact dust disc producing a thermal emission that is largely dominated by small grains located between 0.1 and 0.3 AU from Vega and accounting for 1.23 +/- 0.45% of the near-infrared stellar flux for our best-fit model. This flux ratio is shown to vary slightly with the geometry of the model used to fit our interferometric data (variations within +/-0.19%). Initially revealed by K-band CHARA/FLUOR observations, the presence of hot exozodiacal dust in the vicinity of Vega is confirmed by our H-band IOTA/IONIC measurements at the 3-sigma level. Whereas the origin of the dust is still uncertain, its presence and the possible connection with the outer disc suggest that the Vega system is currently undergoing major dynamical perturbations.

preprint2011arXiv

Sparse Aperture Masking of Massive Stars

We present the earliest results of our NACO/VLT sparse aperture masking (SAM) campaign to search for binarity in a sample of 60 O-type stars. We detect Delta Ks < 5 mag companions for 20-25% of our targets with separations in the range 30-100 mas (typically, 40 - 200 A.U.). Most of these companions were unknown, shedding thus new light on the multiplicity properties of massive stars in a separation and brightness regime that has been difficult to explore so far. Adding detections from other techniques (spectroscopy, interferometry, speckle, lucky imaging, AO), the fraction of O stars with at least one companion is 85% (51/60 targets). This is the largest multiplicity fraction ever found.

preprint2011arXiv

The non-thermal radio emitter HD 93250 resolved by long baseline interferometry

As the brightest O-type X-ray source in the Carina nebula, HD 93250 (O4 III(fc)) is X-ray overluminous for its spectral type and has an unusually hard X-ray spectrum. Two different scenarios have been invoked to explain its X-ray properties: wind-wind interaction and magnetic wind confinement. Yet, HD 93250 shows absolutely constant radial velocities over time scales of years suggesting either a single star, a binary system seen pole-one or a very long period and/or highly eccentric system. Using the ESO Very Large Telescope Interferometer, we resolved HD 93250 as a close pair with similar components. We measured a near-infrared flux ratio of 0.8+/-0.1 and a separation of 1.5+/-0.2 x 10E-03 arcsec. At the distance of Carina, this corresponds to a projected physical distance of 3.5 A.U. While a quantitative investigation would require a full characterization of the orbit, the binary nature of HD 93250 allows us to qualitatively explain both its X-ray flux and hardness and its non-thermal radio emission in the framework of a colliding wind scenario. We also discuss various observational biases. We show that, due to line-blending of two similar spectral components, HD 93250 could have a period as short as one to several years despite the lack of measurable radial velocity variations.

preprint2011arXiv

The size of ABDorA from VLTI/AMBER interferometry

The pre-main sequence (PMS) star ABDorA is the main component of the quadruple system ABDoradus. The precise determination of the mass and photometry of the close companion to ABDorA, ABDorC, has provided an important benchmark for calibration of theoretical evolutionary models of low-mass stars. The limiting factor to the precision of this calibration is the age of the system, as both the mass and luminosity of ABDorA and C are well monitored by other ongoing programs. In this paper we present VLTI/AMBER observations of ABDorA which provide a direct measurement of the size of this star, 0.96+/-0.06 Rsun. The latter estimate, combined with other fundamental parameters also measured for this star, allows a precise test of PMS evolutionary models using both H-R diagrams and mass-radius relationships. We have found that our radius measurement is larger than that predicted by the models, which we interpret as an evidence of the oversizing produced by the strong magnetic activity of ABDorA. Considering, at least partially, this magnetic effect, theoretical isochrones have been used to derive constraints to the age of ABDorA, favouring an age about 40-50 Myr for this system. Older ages are not completely excluded by our data.

preprint2009arXiv

Mass-radius relation of low and very low-mass stars revisited with the VLTI

We measured the radii of 7 low and very low-mass stars using long baseline interferometry with the VLTI interferometer and its VINCI and AMBER near-infrared recombiners. We use these new data, together with literature measurements, to examine the luminosity- radius and mass-radius relations for K and M dwarfs. The precision of the new interferometric radii now competes with what can be obtained for double-lined eclipsing binaries. Interferometry provides access to much less active stars, as well as to stars with much better measured distances and luminosities, and therefore complements the information obtained from eclipsing systems. The radii of magnetically quiet late-K to M dwarfs match the predictions of stellar evolution models very well, providing direct confirmation that magnetic activity explains the discrepancy that was recently found for magnetically active eclipsing systems. The radii of the early K dwarfs are well reproduced for a mixing length parameter that approaches the solar value, as qualitatively expected.