Source author record

P. Kervella

P. Kervella 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

79works
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

79 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.

preprint2025arXiv

ATOMIUM: Continuum emission and evidence of dust enhancement from binary motion

Low- and intermediate-mass stars on the asymptotic giant branch (AGB) account for a significant portion of the dust and chemical enrichment in their host galaxy. Here we present ALMA observations of the continuum emission at 1.24 mm around a sample of 17 stars from the ATOMIUM survey. From our analysis of the stellar contributions to the continuum flux, we find that the semi-regular variables all have smaller physical radii and fainter monochromatic luminosities than the Mira variables. Comparing these properties with pulsation periods, we find a positive trend between stellar radius and period only for the Mira variables with periods above 300 days and a positive trend between the period and the monochromatic luminosity only for the red supergiants and the most extreme AGB stars with periods above 500 days. We find that the continuum emission at 1.2 mm can be classified into four groups. "Featureless" continuum emission is confined to the (unresolved) regions close to the star for five stars in our sample, relatively uniform extended flux is seen for four stars, tentative elongated features are seen for three stars, and the remaining five stars have unique or unusual morphological features in their continuum maps. These features can be explained by binary companions to 10 out of the 14 AGB stars in our sample. Based on our results we conclude that there are two modes of dust formation: well established pulsation-enhanced dust formation and our newly proposed companion-enhanced dust formation. If the companion is located close to the AGB star, in the wind acceleration region, then additional dust formed in the wake of the companion can increase the mass lost through the dust driven wind. This explains the different dust morphologies seen around our stars and partly accounts for a large scatter in literature mass-loss rates, especially among semiregular stars with small pulsation periods.

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.

preprint2021arXiv

ATOMIUM: ALMA tracing the origins of molecules in dust forming oxygen rich M-type stars: Motivation, sample, calibration, and initial results

This overview paper presents ATOMIUM, a Large Programme in Cycle 6 with the Atacama Large Millimeter-submillimeter Array (ALMA). The goal of ATOMIUM is to understand the dynamics and the gas phase and dust formation chemistry in the winds of evolved asymptotic giant branch (AGB) and red supergiant (RSG) stars. A more general aim is to identify chemical processes applicable to other astrophysical environments. 17 oxygen-rich AGB and RSG stars spanning a range in (circum)stellar parameters and evolutionary phases were observed in a homogeneous observing strategy allowing for an unambiguous comparison. Data were obtained between 213.83 and 269.71 GHz at high (0.025-0.050 arcsec), medium (0.13-0.24 arcsec), and low (about 1 arcsec) angular resolution. The sensitivity per 1.3 km/s channel was 1.5-5 mJy/beam. 13 molecules were designated as primary molecules in the survey: CO, SiO, AlO, AlOH, TiO, TiO2, HCN, SO, SO2, SiS, CS, H2O, and NaCl. The scientific motivation, survey design, sample properties, data reduction, and an overview of the data products are described; and we highlight one scientific result - the wind kinematics of the ATOMIUM sources. The ATOMIUM sources often have a slow wind acceleration, and a fraction of the gas reaches a velocity which can be up to a factor of two times larger than previously reported terminal velocities assuming isotropic expansion, and the wind kinematic profiles establish that the radial velocity described by the momentum equation for a spherical wind structure cannot capture the complexity of the velocity field. In 15 sources, some molecular transitions other than 12CO v=0 J=2-1 reach a higher outflow velocity, with a spatial emission zone that is often greater than 30 stellar radii, but much less than the extent of CO. Binary interaction with a (sub)stellar companion might (partly) explain the non-monotonic behaviour of the projected velocity field.

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.

preprint2021arXiv

The inner circumstellar dust of the red supergiant Antares as seen with VLT/SPHERE/ZIMPOL

The processes by which red supergiants lose mass are not fully understood thus-far and their mass-loss rates lack theoretical constraints. The ambient surroundings of the nearby M0.5 Iab star Antares offers an ideal environment to obtain detailed empirical information on the outflow properties at its onset, and hence indirectly, on the mode(s) of mass loss. We present and analyse optical VLT/SPHERE/ZIMPOL polarimetric imaging with angular resolution down to 23 milli-arcsec, sufficient to spatially resolve both the stellar disk and its direct surroundings. We detect a conspicuous feature in polarised intensity that we identify as a clump containing dust, which we characterise through 3D radiative transfer modelling. The clump is positioned behind the plane of the sky, therefore has been released from the backside of the star, and its inner edge is only 0.3 stellar radii above the surface. The current dust mass in the clump is $1.3^{+0.2}_{-1.0} \times 10^{-8}$ M$_{\odot}$, though its proximity to the star implies that dust nucleation is probably still ongoing. The ejection of clumps of gas and dust makes a non-negligible contribution to the total mass lost from the star which could possibly be linked to localised surface activity such as convective motions or non-radial pulsations.

preprint2020arXiv

Calibrating the surface brightness - color relation for late-type red giants stars in the visible domain using VEGA/CHARA interferometric observations

The surface brightness - color relationship (SBCR) is a poweful tool for determining the angular diameter of stars from photometry. It was for instance used to derive the distance of eclipsing binaries in the Large Magellanic Cloud (LMC), which led to its distance determination with an accuracy of 1%. We calibrate the SBCR for red giant stars in the 2.1 < V-K < 2.5 color range using homogeneous VEGA/CHARA interferometric data secured in the visible domain, and compare it to the relation based on infrared interferometric observations, which were used to derive the distance to the LMC. Observations of eight G-K giants were obtained with the VEGA/CHARA instrument. The derived limb-darkened angular diameters were combined with a homogeneous set of infrared magnitudes in order to constrain the SBCR. The average precision we obtain on the limb-darkened angular diameters of the eight stars in our sample is 2.4%. For the four stars in common observed by both VEGA/CHARA and PIONIER/VLTI, we find a 1 sigma agreement for the angular diameters. The SBCR we obtain in the visible has a dispersion of 0.04 magnitude and is consistent with the one derived in the infrared (0.018 magnitude). The consistency of the infrared and visible angular diameters and SBCR reinforces the result of 1% precision and accuracy recently achieved on the distance of the LMC using the eclipsing-binary technique. It also indicates that it is possible to combine interferometric observations at different wavelengths when the SBCR is calibrated.

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

Imaging the expanding knotty structure in the close environment of the LBV star $η$ Carinae

$η$~Car is one of the most massive stars in the Galaxy. It underwent a massive eruption in the 19th century, which produced the impressive bipolar Homunculus nebula now surrounding it. The central star is an eccentric binary with a period of 5.54\,years. Although the companion has not been detected directly, it causes time-variable ionization and colliding-wind X-ray emission. By characterizing the complex structure and kinematics of the ejecta close to the star, we aim to constrain past and present mass loss of $η$~Car. $η$~Car is observed with the extreme adaptive optics instrument SPHERE at the Very Large Telescope, using its polarimetric mode in the optical with the ZIMPOL camera. A spatial resolution of 20\,mas was achieved, i.e. very close to the presumed 13 mas apastron separation of the companion star. We detect new structures within the inner arcsecond to the star (2\,300\,au at a 2.3\,kpc distance). We can relate these structures to the eruption near 1890 by tracking their proper motions derived from our new images and historical images over a 30\,years time span. Besides, we find a fan-shaped structure in the inner 200~au to the star in the H$α$ line, that could potentially be associated with the wind collision zone of the two stars.

preprint2020arXiv

Inner dusty envelope of the AGB stars W~Hydrae, SW\,Virginis, and R~Crateris using SPHERE/ZIMPOL

Stars with initial masses between $\sim0.8$ and 8~$M_\odot$ present copious mass loss during the asymptotic giant branch (AGB) at the end of their lives. The accepted mass-loss mechanism requires radiation pressure acting on dust grains that form in the extended AGB stellar atmospheres. The details of this process are not yet well understood, however. Using the extreme-adaptive-optics imager and polarimeter SPHERE/ZIMPOL, we observed light polarised by grains around W\,Hya, SW\,Vir, and R\,Crt, which have mass-loss rates between 10$^{-7}$ and 10$^{-6}~M_\odot~{\rm yr^{-1}}$. We find the distribution of dust to be asymmetric around the three targets. A biconical morphology is seen for R Crt, with a position angle that is very similar to those inferred from interferometric observations of maser emission and of mid-infrared continuum emission. The cause of the biconical outflow cannot be directly inferred from the ZIMPOL data. The dust grains polarise light more efficiently at 0.65~$μ$m for R\,Crt and SW\,Vir and at 0.82~$μ$m for W\,Hya. This indicates that at the time of the observations, the grains around SW\,Vir and R\,Crt had sizes $< 0.1~μ$m, while those around W\,Hya were larger, with sizes $\gtrsim 0.1~μ$m. The asymmetric distribution of dust around R\,Crt makes the interpretation more uncertain for this star, however. We find that polarised light is produced already from within the visible photosphere of W~Hya, which we reproduce using models with an inner dust shell that is optically thick to scattering. The radial profile of the polarised light observed around W\,Hya reveal a steep dust density profile. We find the wind-acceleration region of W\,Hya to extend to at least $\sim 7~R_\star$, in agreement with theoretical predictions of acceleration up to $\sim 12~R_\star$.

preprint2020arXiv

Pulsating chromosphere of classical Cepheids. Calcium infrared triplet and H$α$ profile variations

It has been shown recently that the infrared emission of Cepheids, which is constant over the pulsation cycle, might be due to a pulsating shell of ionized gas of about 15\% of the stellar radius, which could be attributed to the chromospheric activity of Cepheids. The aim of this paper is to investigate the dynamical structure of the chromosphere of Cepheids along the pulsation cycle and quantify its size. We present H$α$ and Calcium Near InfraRed triplet (Ca IR) profile variations using high-resolution spectroscopy with the UVES spectrograph of a sample of 24 Cepheids with a good period coverage from $\approx$ 3 to 60 days. After a qualitative analysis of the spectral lines profiles, we quantify the Van Hoof effect (velocity gradient between the H$α$ and Ca IR) as a function of the period of the Cepheids. Then, we use the Schwarzschild mechanism (a line doubling due to a shock wave) to quantify the size of the chromosphere. We find a significant Van Hoof effect for Cepheids with period larger than $P=10$ days, in particular H$α$ lines are delayed with a velocity gradient up to $Δv \approx$30 km/s compared to Ca IR. We find that the size of the chromosphere of long-period Cepheids is of at least $\approx$ 50\% of the stellar radius, which is consistent at first order with the size of the shell made of ionized gas previously found from the analysis of infrared excess. Last, for most of the long-period Cepheids in the sample, we report a motionless absorption feature in the H$α$ line that we attribute to a circumstellar envelope that surrounds the chromosphere. Analyzing the Ca~IR lines of Cepheids is of importance to potentially unbias the period-luminosity relation from their infrared excess, particularly in the context of forthcoming observations from the Radial Velocity Spectrometer (RVS) on board \textit{Gaia}, that could be sensitive to their chromosphere.

preprint2020arXiv

Retrieving scattering clouds and disequilibrium chemistry in the atmosphere of HR 8799e

Clouds are ubiquitous in exoplanet atmospheres and represent a challenge for the model interpretation of their spectra. Complex cloud models are too numerically costly for generating a large number of spectra, while more efficient models may be too strongly simplified. We aim to constrain the atmospheric properties of the directly imaged planet HR 8799e with a free retrieval approach. We use our radiative transfer code petitRADTRANS for generating spectra, which we couple to the PyMultiNest tool. We added the effect of multiple scattering which is important for treating clouds. Two cloud model parameterizations are tested: the first incorporates the mixing and settling of condensates, the second simply parameterizes the functional form of the opacity. In mock retrievals, using an inadequate cloud model may result in atmospheres that are more isothermal and less cloudy than the input. Applying our framework on observations of HR 8799e made with the GPI, SPHERE and GRAVITY, we find a cloudy atmosphere governed by disequilibrium chemistry, confirming previous analyses. We retrieve that ${\rm C/O}=0.60_{-0.08}^{+0.07}$. Other models have not yet produced a well constrained C/O value for this planet. The retrieved C/O values of both cloud models are consistent, while leading to different atmospheric structures: cloudy, or more isothermal and less cloudy. Fitting the observations with the self-consistent Exo-REM model leads to comparable results, while not constraining C/O. With data from the most sensitive instruments, retrieval analyses of directly imaged planets are possible. The inferred C/O ratio of HR 8799e is independent of the cloud model and thus appears to be a robust. This C/O is consistent with stellar, which could indicate that the HR 8799e formed outside the CO$_2$ or CO iceline. As it is the innermost planet of the system, this constraint could apply to all HR 8799 planets.

preprint2020arXiv

Searching for the near infrared counterpart of Proxima c using multi-epoch high contrast SPHERE data at VLT

Proxima Centauri is known to host an earth-like planet in its habitable zone; very recently a second candidate planet was proposed based on radial velocities. At quadrature, the expected projected separation of this new candidate is larger than 1 arcsec, making it a potentially interesting target for direct imaging. While difficult, identification of the optical counterpart of this planet would allow detailed characterization of the closest planetary system. We searched for a counterpart in SPHERE images acquired during four years through the SHINE survey. In order to account for the large orbital motion of the planet, we used a method that assumes the circular orbit obtained from radial velocities and exploits the sequence of observations acquired close to quadrature in the orbit. We checked this with a more general approach that considers keplerian motion, K-stacker. We did not obtain a clear detection. The best candidate has S/N=6.1 in the combined image. A statistical test suggests that the probability that this detection is due to random fluctuation of noise is < 1% but this result depends on the assumption that distribution of noise is uniform over the image. The position of this candidate and the orientation of its orbital plane fit well with observations in the ALMA 12m array image. However, the astrometric signal expected from the orbit of the candidate we detected is 3-sigma away from the astrometric motion of Proxima as measured from early Gaia data. This, together with the unexpectedly high flux associated with our direct imaging detection, means we cannot confirm that our candidate is indeed Proxima c. On the other hand, if confirmed, this would be the first observation in imaging of a planet discovered from radial velocities and the second one (after Fomalhaut b) of reflecting circumplanetary material. Further confirmation observations should be done as soon as possible.

preprint2020arXiv

SPHERE+: Imaging young Jupiters down to the snowline

SPHERE (Beuzit et al,. 2019) has now been in operation at the VLT for more than 5 years, demonstrating a high level of performance. SPHERE has produced outstanding results using a variety of operating modes, primarily in the field of direct imaging of exoplanetary systems, focusing on exoplanets as point sources and circumstellar disks as extended objects. The achievements obtained thus far with SPHERE (~200 refereed publications) in different areas (exoplanets, disks, solar system, stellar physics...) have motivated a large consortium to propose an even more ambitious set of science cases, and its corresponding technical implementation in the form of an upgrade. The SPHERE+ project capitalizes on the expertise and lessons learned from SPHERE to push high contrast imaging performance to its limits on the VLT 8m-telescope. The scientific program of SPHERE+ described in this document will open a new and compelling scientific window for the upcoming decade in strong synergy with ground-based facilities (VLT/I, ELT, ALMA, and SKA) and space missions (Gaia, JWST, PLATO and WFIRST). While SPHERE has sampled the outer parts of planetary systems beyond a few tens of AU, SPHERE+ will dig into the inner regions around stars to reveal and characterize by mean of spectroscopy the giant planet population down to the snow line. Building on SPHERE's scientific heritage and resounding success, SPHERE+ will be a dedicated survey instrument which will strengthen the leadership of ESO and the European community in the very competitive field of direct imaging of exoplanetary systems. With enhanced capabilities, it will enable an even broader diversity of science cases including the study of the solar system, the birth and death of stars and the exploration of the inner regions of active galactic nuclei.

preprint2020arXiv

TESS observations of Cepheid stars: first light results

We present the first analysis of Cepheid stars observed by the TESS space mission in Sectors 1 to 5. Our sample consists of 25 pulsators: ten fundamental mode, three overtone and two double-mode classical Cepheids, plus three Type II and seven anomalous Cepheids. The targets were chosen from fields with different stellar densities, both from the Galactic field and from the Magellanic System. Three targets have 2-minute cadence light curves available by the TESS Science Processing Operations Center: for the rest, we prepared custom light curves from the full-frame images with our own differential photometric FITSH pipeline. Our main goal was to explore the potential and the limitations of TESS concerning the various subtypes of Cepheids. We detected many low amplitude features: weak modulation, period jitter, and timing variations due to light-time effect. We also report signs of non-radial modes and the first discovery of such a mode in an anomalous Cepheid, the overtone star XZ Cet, which we then confirmed with ground-based multicolor photometric measurements. We prepared a custom photometric solution to minimize saturation effects in the bright fundamental-mode classical Cepheid, $β$ Dor with the lightkurve software, and we revealed strong evidence of cycle-to-cycle variations in the star. In several cases, however, fluctuations in the pulsation could not be distinguished from instrumental effects, such as contamination from nearby sources which also varies between sectors. Finally, we discuss how precise light curve shapes will be crucial not only for classification purposes but also to determine physical properties of these stars.

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

A realistic two-dimensional model of Altair

Fast rotation is responsible for important changes in the structure and evolution of stars. Optical long baseline interferometry now permits the study of its effects on the stellar surface, mainly gravity darkening and flattening. We aim to determine the fundamental parameters of the fast-rotating star Altair, in particular its evolutionary stage, mass, and differential rotation, using state-of-the-art stellar interior and atmosphere models together with interferometric, spectroscopic, and asteroseismic observations. We use ESTER 2D stellar models to produce the relevant surface parameters needed to create intensity maps from atmosphere models. Interferometric and spectroscopic observables are computed from these intensity maps and several stellar parameters are then adjusted using the MCMC algorithm Emcee. We determined Altair's equatorial radius to be 2.008 +/- 0.006 Rsun, the position angle 301.1 +/- 0.3 degrees, the inclination 50.7 +/- 1.2 degrees, and the equatorial angular velocity 0.74 +/- 0.01 times the Keplerian angular velocity. This angular velocity leads to a flattening of 0.220 +/- 0.003. We also deduce from the spectroscopically derived vsini ~ 243 km/s, a true equatorial velocity of ~314 km/s corresponding to a rotation period of 7h46m (~3 c/d). The data also impose a strong correlation between mass, metallicity, hydrogen abundance, and core evolution. Thanks to asteroseismic data, we constrain the mass of Altair to 1.86 +/- 0.03 Msun and further deduce its metallicity Z = 0.019 and its core hydrogen mass fraction Xc = 0.71, assuming an initial solar hydrogen mass fraction X = 0.739. These values suggest that Altair is ~100 Myrs old. Finally, the 2D ESTER model also gives the internal differential rotation of Altair, showing that its core rotates approximately 50% faster than the envelope, while the surface differential rotation does not exceed 6%.

preprint2019arXiv

A thin shell of ionized gas explaining the IR excess of classical Cepheids

Despite observational evidences, InfraRed (IR) excess of classical Cepheids are seldom studied and poorly understood, but probably induces systematics on the Period-Luminosity (PL) relation used in the calibration of the extragalactic distance scale. This study aims to understand the physical origin of the IR excess found in the spectral energy distribution (SED) of 5 Cepheids : RS Pup (P=41.46d), zeta Gem (P=10.15d), eta Aql (P=7.18d), V Cen (P=5.49d) and SU Cyg (P=3.85d). A time series of atmospheric models along the pulsation cycle are fitted to a compilation of data, including optical and near-IR photometry, Spitzer spectra (secured at a specific phase), interferometric angular diameters, effective temperature and radial velocity measurements. Herschel images in two bands are also analyzed qualitatively. In this fitting process, based on the SPIPS algorithm, a residual is found in the SED, whatever the pulsation phase, and for wavelengths larger than about $1.2μ$m, which corresponds to the so-determined infrared excess of Cepheids. This IR excess is then corrected from interstellar medium absorption in order to infer or not the presence of dust shells, and is finally used in order to fit a model of a shell of ionized gas. For all Cepheids, we find a continuum IR excess increasing up to about -0.1 magnitudes at 30$μ$m, which cannot be explained by a hot or cold dust model of CircumStellar Environment (CSE). We show, for the first time, that the IR excess of Cepheids can be explained by free-free emission from a thin shell of ionized gas, with a thickness of about 15% of the star radius, a mass of $10^{-9}-10^{-7}$ Msol and a temperature ranging from 3500 to 4500K. This result has to be tested with interferometers operating in visible, in the mid-IR or in the radio domain. The impact of such CSEs of ionized gas on the PL relation of Cepheids needs also more investigations.

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

ALMA observations of the nearby AGB star L2 Puppis - I. Mass of the central star and detection of a candidate planet

Six billion years from now, while evolving on the asymptotic giant branch (AGB), the Sun will metamorphose from a red giant into a beautiful planetary nebula. This spectacular evolution will impact the Solar System planets, but observational confirmations of the predictions of evolution models are still elusive as no planet orbiting an AGB star has yet been discovered. The nearby AGB red giant L2 Puppis (d = 64 pc) is surrounded by an almost edge-on circumstellar dust disk. We report new observations with ALMA at very high angular resolution (18 x 15 mas) in band 7 (f ~ 350 GHz) that allow us to resolve the velocity profile of the molecular disk. We establish that the gas velocity profile is Keplerian within the central cavity of the dust disk, allowing us to derive the mass of the central star L2 Pup A, mA = 0.659 +/- 0.011 +/- 0.041 Msun (+/- 6.6%). From evolutionary models, we determine that L2 Pup A had a near-solar main sequence mass, and is therefore a close analog of the future Sun in 5 to 6 Gyr. The continuum map reveals the presence of a secondary source (B) at a radius of 2 AU contributing fB/ fA = 1.3 +/- 0.1% of the flux of the AGB star. L2 Pup B is also detected in CO emission lines at a radial velocity of vB = 12.2 +/- 1.0 km/s. The close coincidence of the center of rotation of the gaseous disk with the position of the continuum emission from the AGB star allows us to constrain the mass of the companion to mB = 12 +/- 16 MJup. L2 Pup B is most likely a planet or low mass brown dwarf with an orbital period around 5 years. Its continuum brightness and molecular emission suggest that it may be surrounded by an extended molecular atmosphere or an accretion disk. L2 Pup therefore emerges as a promising vantage point on the distant future of our Solar System.

preprint2016arXiv

Close stellar conjunctions of alpha Centauri A and B until 2050 - An mK = 7.8 star may enter the Einstein ring of alpha Cen A in 2028

The rapid proper motion of the alpha Cen pair (3.7 arcsec/yr) and its location close to the galactic plane on a rich stellar background combine to make them excellent candidates for stellar conjunctions with distant stars. Adding new astrometry to archival data, we have refined the orbital parameters, barycentric proper motion and parallax of alpha Cen and compute its apparent trajectory on sky over the coming decades. We present a catalog of the expected close conjunctions until 2050. An exceptional event will take place in early May 2028, when alpha Cen A will come within 0.015 +/- 0.135 arcseconds of the mK = 7.8 star 2MASS 14392160-6049528 (hereafter S5). In terms of impact parameter and contrast, this is the most favorable stellar conjunction of alpha Cen within at least the next three decades. With an angular diameter of LD = 0.47 +/- 0.05 mas, it is likely that S5 is a red giant or supergiant located at several kiloparsecs. The approached stars will act as moving light probes in transmission through the environment of alpha Cen. The observation of these close conjunctions holds great promises to search for planets and other low mass objects in the alpha Cen system using photometry and astrometry. The relativistic deflection of the approached star images will be in the milliarcsecond range. The small impact parameter of the conjunction with S5 means that this star has a probability of 45% of entering the Einstein ring of alpha Cen A. The gravitational amplification of the flux of S5 could reach a factor five for the combination of the two lensed images. The proper motion, orbital parameters and parallax of alpha Cen will be measurable with an extreme accuracy from differential astrometry with the S stars. This will be valuable, for example to prepare the Breakthrough Starshot initiative to send interstellar nanocrafts to alpha Centauri.

preprint2016arXiv

Discovery of a complex linearly polarized spectrum of Betelgeuse dominated by depolarization of the continuum

Betelgeuse is an M supergiant that harbors spots and giant granules at its surface and presents linear polarization of its continuum. We have previously discovered linear polarization signatures associated with individual lines in the spectra of cool and evolved stars. Here, we investigate whether a similar linearly polarized spectrum exists for Betelgeuse. We used the spectropolarimeter Narval, combining multiple polarimetric sequences to obtain high signal-to-noise ratio spectra of individual lines, as well as the least-squares deconvolution (LSD) approach. We have discovered the existence of a linearly polarized spectrum for Betelgeuse, detecting a rather strong signal (at a few times 10$^{-4}$ of the continuum intensity level), both in individual lines and in the LSD profiles. Studying its properties and the signal observed for the resonant \ion{Na}{i}\,D lines, we conclude that we are mainly observing depolarization of the continuum by the absorption lines. The linear polarization of the Betelgeuse continuum is due to the anisotropy of the radiation field induced by brightness spots at the surface and Rayleigh scattering in the atmosphere. We have developed a geometrical model to interpret the observed polarization, from which we infer the presence of two brightness spots and their positions on the surface of Betelgeuse. We show that applying the model to each velocity bin along the Stokes Q and U profiles allows the derivation of a map of the bright spots. We use the Narval linear polarization observations of Betelgeuse obtained over a period of 1.4 years to study the evolution of the spots. Our study of the linearly polarized spectrum of Betelgeuse provides a novel method for studying the evolution of brightness spots at its surface and complements quasi-simultaneous observations obtained with PIONIER at the VLTI.

preprint2016arXiv

Multiplicity of Galactic Cepheids from long-baseline interferometry~III. Sub-percent limits on the relative brightness of a close companion of $δ$~Cephei

We report new CHARA/MIRC interferometric observations of the Cepheid archetype $δ$ Cep, which aimed at detecting the newly discovered spectroscopic companion. We reached a maximum dynamic range $ΔH $ = 6.4, 5.8, and 5.2 mag, respectively within the relative distance to the Cepheid $r < 25$ mas, $25 < r < 50$ mas and $50 < r < 100$ mas. Our observations did not show strong evidence of a companion. We have a marginal detection at $3σ$ with a flux ratio of 0.21\%, but nothing convincing as we found other possible probable locations. We ruled out the presence of companion with a spectral type earlier than F0V, A1V and B9V, respectively for the previously cited ranges $r$. From our estimated sensitivity limits and the Cepheid light curve, we derived lower-limit magnitudes in the $H$ band for this possible companion to be $H_\mathrm{comp} > 9.15, 8.31$ and 7.77 mag, respectively for $r < 25$ mas, $25 < r < 50$ mas and $50 < r < 100$ mas. We also found that to be consistent with the predicted orbital period, the companion has to be located at a projected separation $< 24$ mas with a spectral type later than a F0V star.

preprint2016arXiv

Pluto's atmosphere from the 29 June 2015 ground-based stellar occultation at the time of the New Horizons flyby

We present results from a multi-chord Pluto stellar occultation observed on 29 June 2015 from New Zealand and Australia. This occurred only two weeks before the NASA New Horizons flyby of the Pluto system and serves as a useful comparison between ground-based and space results. We find that Pluto's atmosphere is still expanding, with a significant pressure increase of 5$\pm$2\% since 2013 and a factor of almost three since 1988. This trend rules out, as of today, an atmospheric collapse associated with Pluto's recession from the Sun. A central flash, a rare occurrence, was observed from several sites in New Zealand. The flash shape and amplitude are compatible with a spherical and transparent atmospheric layer of roughly 3~km in thickness whose base lies at about 4~km above Pluto's surface, and where an average thermal gradient of about 5 K~km$^{-1}$ prevails. We discuss the possibility that small departures between the observed and modeled flash are caused by local topographic features (mountains) along Pluto's limb that block the stellar light. Finally, using two possible temperature profiles, and extrapolating our pressure profile from our deepest accessible level down to the surface, we obtain a possible range of 11.9-13.7~$μ$bar for the surface pressure.

preprint2016arXiv

Study of the inner dust envelope and stellar photosphere of the AGB star R Doradus using SPHERE/ZIMPOL

We use high-angular-resolution images obtained with SPHERE/ZIMPOL to study the photosphere, the warm molecular layer, and the inner wind of the close-by oxygen-rich AGB star R Doradus. We present observations in filters V, cntH$α$, and cnt820 and investigate the surface brightness distribution of the star and of the polarised light produced in the inner envelope. Thanks to second-epoch observations in cntH$α$, we are able to see variability on the stellar photosphere. We find that in the first epoch the surface brightness of R Dor is asymmetric in V and cntH$α$, the filters where molecular opacity is stronger, while in cnt820 the surface brightness is closer to being axisymmetric. The second-epoch observations in cntH$α$ show that the morphology of R Dor changes completely in a timespan of 48 days to a more axisymmetric and compact configuration. The polarised intensity is asymmetric in all epochs and varies by between a factor of 2.3 and 3.7 with azimuth for the different images. We fit the radial profile of the polarised intensity using a spherically symmetric model and a parametric description of the dust density profile, $ρ(r)=ρ_\circ r^{-n}$. On average, we find exponents of $- 4.5 \pm 0.5$ that correspond to a much steeper density profile than that of a wind expanding at constant velocity. The dust densities we derive imply an upper limit for the dust-to-gas ratio of $\sim 2\times10^{-4}$ at 5.0 $R_\star$. Given the uncertainties in observations and models, this value is consistent with the minimum values required by wind-driving models for the onset of a wind, of $\sim 3.3\times10^{-4}$. However, if the steep density profile we find extends to larger distances from the star, the dust-to-gas ratio will quickly become too small for the wind of R Dor to be driven by the grains that produce the scattered light.

preprint2016arXiv

The close circumstellar environment of Betelgeuse. IV. VLTI/PIONIER interferometric monitoring of the photosphere

Context. The mass-loss mechanism of cool massive evolved stars is poorly understood. The proximity of Betelgeuse makes it an appealing target to study its atmosphere, map the shape of its envelope, and follow the structure of its wind from the photosphere out to the interstellar medium. Aims. A link is suspected between the powerful convective motions in Betelgeuse and its mass loss. We aim to constrain the spatial structure and temporal evolution of the convective pattern on the photosphere and to search for evidence of this link. Methods. We report new interferometric observations in the infrared H band using the VLTI/PIONIER instrument. We monitored the photosphere of Betelgeuse between 2012 January and 2014 November to look for evolutions that may trigger the outflow. Results. Our interferometric observations at low spatial frequencies are compatible with the presence of a hot spot on the photosphere that has a characteristic width of one stellar radius. It appears to be superposed on the smaller scale convective pattern. In the higher spatial frequency domain, we observe a significant difference between the observations and the predictions of 3D hydrodynamical simulations. Conclusions. We bring new evidence for the presence of a convective pattern in the photosphere of red supergiants. The inferred hot spot is probably the top of a giant convection cell although an asymmetric extension of the star cannot be excluded by these interferometric observations alone. The properties of the observed surface features show a stronger contrast and inhomogeneity as predicted by 3D radiative hydrodynamical simulations. We propose that the large observed feature is modifying the signature of the convective pattern at the surface of the star in a way that simulations cannot reproduce.

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.

preprint2016arXiv

VEGA/CHARA interferometric observations of Cepheids. I. A resolved structure around the prototype classical Cepheid delta Cep in the visible spectral range

The B-W method is used to determine the distance of Cepheids and consists in combining the angular size variations of the star, as derived from infrared surface-brightness relations or interferometry, with its linear size variation, as deduced from visible spectroscopy using the projection factor. While many Cepheids have been intensively observed by infrared beam combiners, only a few have been observed in the visible. This paper is part of a project to observe Cepheids in the visible with interferometry as a counterpart to infrared observations already in hand. Observations of delta Cep itself were secured with the VEGA/CHARA instrument over the full pulsation cycle of the star. These visible interferometric data are consistent in first approximation with a quasi-hydrostatic model of pulsation surrounded by a static circumstellar environment (CSE) with a size of theta_cse=8.9 +/- 3.0 mas and a relative flux contribution of f_cse=0.07+/-0.01. A model of visible nebula (a background source filling the field of view of the interferometer) with the same relative flux contribution is also consistent with our data at small spatial frequencies. However, in both cases, we find discrepancies in the squared visibilities at high spatial frequencies (maximum 2sigma) with two different regimes over the pulsation cycle of the star, phi=0.0-0.8 and phi=0.8-1.0. We provide several hypotheses to explain these discrepancies, but more observations and theoretical investigations are necessary before a firm conclusion can be drawn. For the first time we have been able to detect in the visible domain a resolved structure around delta~Cep. We have also shown that a simple model cannot explain the observations, and more work will be necessary in the future, both on observations and modelling.

preprint2015arXiv

Interferometric imaging diagnostics of X Hya's circumstellar environment

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

preprint2015arXiv

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

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

preprint2015arXiv

Robust high-contrast companion detection from interferometric observations. The CANDID algorithm and an application to six binary Cepheids

Long-baseline interferometry is an important technique to spatially resolve binary or multiple systems in close orbits. By combining several telescopes together and spectrally dispersing the light, it is possible to detect faint components around bright stars. Aims. We provide a rigorous and detailed method to search for high-contrast companions around stars, determine the detection level, and estimate the dynamic range from interferometric observations. We developed the code CANDID (Companion Analysis and Non-Detection in Interferometric Data), a set of Python tools that allows us to search systematically for point-source, high-contrast companions and estimate the detection limit. The search pro- cedure is made on a N x N grid of fit, whose minimum needed resolution is estimated a posteriori. It includes a tool to estimate the detection level of the companion in the number of sigmas. The code CANDID also incorporates a robust method to set a 3σ detection limit on the flux ratio, which is based on an analytical injection of a fake companion at each point in the grid. We used CANDID to search for the companions around the binary Cepheids V1334 Cyg, AX Cir, RT Aur, AW Per, SU Cas, and T Vul. First, we showed that our previous discoveries of the components orbiting V1334 Cyg and AX Cir were detected at > 13 sigmas. The companion around AW Per is detected at more than 15 sigmas with a flux ratio of f = 1.22 +/- 0.30 %. We made a possible detection of the companion orbiting RT Aur with f = 0.22 +/- 0.11 %. It was detected at 3.8σ using the closure phases only, and so more observations are needed to confirm the detection. We also set the detection limit for possible undetected companions. We found that there is no companion with a spectral type earlier than B7V, A5V, F0V, B9V, A0V, and B9V orbiting V1334 Cyg, AX Cir, RT Aur, AW Per, SU Cas, and T Vul, respectively.

preprint2015arXiv

The Araucaria Project: High-precision orbital parallax and masses of the eclipsing binary TZ~Fornacis

Context: Independent distance estimates are particularly useful to check the precision of other distance indicators, while accurate and precise masses are necessary to constrain evolution models. Aim: The goal is to measure the masses and distance of the detached eclipsing-binary TZ~For with a precision level lower than 1\,\% using a fully geometrical and empirical method. Method: We obtained the first interferometric observations of TZ~For with the VLTI/PIONIER combiner, which we combined with new and precise radial velocity measurements to derive its three-dimensional orbit, masses, and distance. Results: The system is well resolved by PIONIER at each observing epoch, which allowed a combined fit with eleven astrometric positions. Our derived values are in a good agreement with previous work, but with an improved precision. We measured the mass of both components to be $M_1 = 2.057 \pm 0.001\,M_\odot$ and $M_2 = 1.958 \pm 0.001\,M_\odot$. The comparison with stellar evolution models gives an age of the system of $1.20 \pm 0.10$\,Gyr. We also derived the distance to the system with a precision level of 1.1\,\%: $d = 185.9 \pm 1.9$\,pc. Such precise and accurate geometrical distances to eclipsing binaries provide a unique opportunity to test the absolute calibration of the surface brightness-colour relation for late-type stars, and will also provide the best opportunity to check on the future Gaia measurements for possible systematic errors.

preprint2015arXiv

The close circumstellar environment of Betelgeuse - III. SPHERE/ZIMPOL visible polarimetry of the inner envelope and photosphere

The physical mechanism through which the outgoing material of massive red supergiants is accelerated above the escape velocity is unclear. Thanks to the transparency of its circumstellar envelope, the nearby red supergiant Betelgeuse gives the opportunity to probe the innermost layers of the envelope of a typical red supergiant down to the photosphere, i.e. where the acceleration of the wind is expected to occur. We took advantage of the SPHERE/ZIMPOL adaptive optics imaging polarimeter to resolve the visible photosphere and close envelope of Betelgeuse. We detect an asymmetric gaseous envelope inside a radius of 2 to 3 times the near-infrared photospheric radius of the star (R*), and a significant Halpha emission mostly contained within 3 R*. From the polarimetric signal, we also identify the signature of dust scattering in an asymmetric and incomplete dust shell located at a similar radius. The presence of dust so close to the star may have a significant impact on the wind acceleration through radiative pressure on the grains. The 3 R* radius emerges as a major interface between the hot gaseous and dusty envelopes. The detected asymmetries strengthen previous indications that the mass loss of Betelgeuse is likely tied to the vigorous convective motions in its atmosphere.

preprint2015arXiv

The dust disk and companion of the nearby AGB star L2 Puppis - SPHERE/ZIMPOL polarimetric imaging at visible wavelengths

The bright southern star L2 Pup is a particularly prominent asymptotic giant branch (AGB) star, as its distance of 64 pc makes it the nearest of its type. We report new adaptive optics observations of L2 Pup at visible wavelengths with the SPHERE/ZIMPOL instrument of the VLT that confirm the presence of the circumstellar dust disk at high inclination discovered recently by Kervella et al. (2014b). The signature of the three-dimensional structure of the disk is clearly observed in the map of the degree of linear polarization pL. We identify the inner rim of the disk through its polarimetric signature at a radius of 6 AU from the AGB star. The ZIMPOL intensity images in the V and R bands also reveal a close-in secondary source at a projected separation of 2 AU from the primary. The identification of the spectral type of this companion is uncertain due to the strong reddening from the disk, but its photometry suggests that it is a late K giant, of comparable mass to the AGB star. We present refined physical parameters for the dust disk derived using the RADMC-3D radiative transfer code. We also interpret the pL map using a simple polarization model to infer the three-dimensional structure of the envelope. Interactions between the inner binary system and the disk apparently form spiral structures that propagate along the orthogonal axis to the disk to form streamers. Two dust plumes propagating orthogonally to the disk are also detected. They originate in the inner stellar system, and are possibly related to the interaction of the wind of the two stars with the material in the disk. Based on the morphology of the envelope of L2 Pup, we propose that this star is at an early stage of the formation of a bipolar planetary nebula.

preprint2015arXiv

The EChO science case

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

preprint2015arXiv

The nearby AGB star L2 Puppis: the birth of a planetary nebula ?

Adaptive optics observations in the infrared (VLT/NACO, Kervella et al. 2014) and visible (VLT/SPHERE, Kervella et al. 2015) domains revealed that the nearby AGB star L2 Pup (d=64 pc) is surrounded by a dust disk seen almost edge-on. Thermal emission from a large dust "loop" is detected at 4 microns up to more than 10 AU from the star. We also detect a secondary source at a separation of 32 mas, whose nature is uncertain. L2 Pup is currently a relatively "young" AGB star, so we may witness the formation of a planetary nebula. The mechanism that breaks the spherical symmetry of mass loss is currently uncertain, but we propose that the dust disk and companion are key elements in the shaping of the bipolar structure. L2 Pup emerges as an important system to test this hypothesis.

preprint2014arXiv

A ring system detected around the Centaur (10199) Chariklo

Until now, rings have been detected in the Solar System exclusively around the four giant planets. Here we report the discovery of the first minor-body ring system around the Centaur object (10199) Chariklo, a body with equivalent radius 124$\pm$9 km. A multi-chord stellar occultation revealed the presence of two dense rings around Chariklo, with widths of about 7 km and 3 km, optical depths 0.4 and 0.06, and orbital radii 391 and 405 km, respectively. The present orientation of the ring is consistent with an edge-on geometry in 2008, thus providing a simple explanation for the dimming of Chariklo's system between 1997 and 2008, and for the gradual disappearance of ice and other absorption features in its spectrum over the same period. This implies that the rings are partially composed of water ice. These rings may be the remnants of a debris disk, which were possibly confined by embedded kilometre-sized satellites.

preprint2014arXiv

ALMA Observations of Anisotropic Dust Mass-loss in the Inner Circumstellar Environment of the Red Supergiant VY Canis Majoris

The processes leading to dust formation and the subsequent role it plays in driving mass loss in cool evolved stars is an area of intense study. Here we present high resolution ALMA Science Verification data of the continuum emission around the highly evolved oxygen-rich red supergiant VY CMa. These data enable us to study the dust in its inner circumstellar environment at a spatial resolution of 129 mas at 321 GHz and 59 mas at 658 GHz, thus allowing us to trace dust on spatial scales down to 11 R$_{\star}$ (71 AU). Two prominent dust components are detected and resolved. The brightest dust component, C, is located 334 mas (61 R$_{\star}$) South East of the star and has a dust mass of at least $2.5\times 10^{-4}$ M$_{\odot}$. It has a dust emissivity spectral index of $β=-0.1$ at its peak, implying that it is optically thick at these frequencies with a cool core of $T_{d}\lesssim 100$ K. Interestingly, not a single molecule in the ALMA data has emission close to the peak of this massive dust clump. The other main dust component, VY, is located at the position of the star and contains a total dust mass of $4.0 \times 10^{-5} $M$_{\odot}$. It also contains a weaker dust feature extending over $60$ R$_{\star}$ to the North with the total component having a typical dust emissivity spectral index of $β=0.7$. We find that at least $17\%$ of the dust mass around VY CMa is located in clumps ejected within a more quiescent roughly spherical stellar wind, with a quiescent dust mass loss rate of $5 \times 10^{-6}$ M$_{\odot} $yr$^{-1}$. The anisotropic morphology of the dust indicates a continuous, directed mass loss over a few decades, suggesting that this mass loss cannot be driven by large convection cells alone.

preprint2014arXiv

ALMA sub-mm maser and dust distribution of VY Canis Majoris

Cool, evolved stars have copious, enriched winds. The structure of these winds and the way they are accelerated is not well known. We need to improve our understanding by studying the dynamics from the pulsating stellar surface to about 10 stellar radii, where radiation pressure on dust is fully effective. Some red supergiants have highly asymmetric nebulae, implicating additional forces. We retrieved ALMA Science Verification data providing images of sub-mm line and continuum emission from VY CMa. This enables us to locate water masers with milli-arcsec precision and resolve the dusty continuum. The 658-, 321- and 325-GHz masers lie in irregular, thick shells at increasing distances from the centre of expansion. For the first time this is confirmed as the stellar position, coinciding with a compact peak offset to the NW of the brightest continuum emission. The maser shells (and dust formation zone) overlap but avoid each other on tens-au scales. Their distribution is broadly consistent with excitation models but the conditions and kinematics appear to be complicated by wind collisions, clumping and asymmetries.

preprint2014arXiv

An interferometric view on binarity and circumstellar envelopes of Cepheids

Optical interferometry is the only technique giving access to milli-arcsecond (mas) resolution at infrared wavelengths. For Cepheids, this is a powerful and unique tool to detect the orbiting companions and the circumstellar envelopes (CSE). CSEs are interesting because they might be used to trace the Cepheid evolution history, and more particularly they could impact the distance scale. Cepheids belonging to binary systems offer an unique opportunity to make progress in resolving the Cepheid mass discrepancy. The combination of spectroscopic and interferometric measurements will allow us to derive the orbital elements, distances, and dynamical masses. Here we focus on recent results using 2- to 6-telescopes beam combiners for the Cepheids X Sgr, T Mon and V1334 Cyg.

preprint2014arXiv

Binary Cepheids from optical interferometry

Classical Cepheid stars have been considered since more than a century as reliable tools to estimate distances in the universe thanks to their Period-Luminosity (P-L) relationship. Moreover, they are also powerful astrophysical laboratories, providing fundamental clues for studying the pulsation and evolution of intermediate-mass stars. When in binary systems, we can investigate the age and evolution of the Cepheid, estimate the mass and distance, and constrain theoretical models. However, most of the companions are located too close to the Cepheid (1-40 mas) to be spatially resolved with a 10-meter class telescope. The only way to spatially resolve such systems is to use long-baseline interferometry. Recently, we have started a unique and long-term interferometric program that aims at detecting and characterizing physical parameters of the Cepheid companions, with as main objectives the determination of accurate masses and geometric distances.

preprint2014arXiv

GCIRS 7, a pulsating M1 supergiant at the Galactic centre. Physical properties and age

The stellar population in the central parsec of the Galaxy is dominated by an old (several Gyr) population, but young, massive stars dominate the luminosity function. We have studied the most luminous of these stars, GCIRS 7, in order to constrain the age of the recent star formation event in the Galactic Centre and to characterise it as an interferometric reference for observations of the Galactic Centre with the instrument GRAVITY, which will equip the Very Large Telescope Interferometer in the near future. We present the first H-band interferometric observations of GCIRS 7, obtained using the PIONIER visitor instrument on the VLTI using the four 8.2-m unit telescopes. In addition, we present unpublished K-band VLTI/AMBER data, build JHKL light-curves based on data spanning 4 decades, and measured the star's effective temperature using SINFONI spectroscopy. GCIRS 7 is marginally resolved at H-band (in 2013: uniform-disk diameter=1.076+/-0.093mas, R=960+/-92Rsun at 8.33+/-0.35kpc). We detect a significant circumstellar contribution at K-band. The star and its environment are variable in brightness and in size. The photospheric H-band variations are well modelled with two periods: P0~470+/-10 days (amplitude ~0.64mag) and long secondary period LSP~2700-2850 days (~1.1mag). As measured from CO equivalent width, <Teff>=3600+/-195K. The size, periods, luminosity (<Mbol>=-8.44+/-0.22) and effective temperature are consistent with an M1 supergiant with an initial mass of 22.5+/-2.5Msun and an age of 6.5-10Myr (depending on rotation). This age is in remarkable agreement with most estimates for the recent star formation event in the central parsec. Caution should be taken when using this star as an interferometric reference as it is variable in size, is surrounded by a variable circumstellar environment and large convection cells may form on its photosphere.

preprint2014arXiv

Multiplicity of Galactic Cepheids from long-baseline interferometry I. CHARA/MIRC detection of the companion of V1334 Cygni

We aim at determining the masses of Cepheids in binary systems, as well as their geometric distances and the flux contribution of the companions. The combination of interferometry with spectroscopy will offer a unique and independent estimate of the Cepheid masses. Using long-baseline interferometry at visible and infrared wavelengths, it is possible to spatially resolve binary systems containing a Cepheid down to milliarcsecond separations. Based on the resulting visual orbit and radial velocities, we can then derive the fundamental parameters of these systems, particularly the masses of the components and the geometric distance. We therefore performed interferometric observations of the first-overtone mode Cepheid V1334 Cyg with the CHARA/MIRC combiner. We report the first detection of a Cepheid companion using long-baseline interferometry. We detect the signature of a companion orbiting V1334 Cyg at two epochs. We measure a flux ratio between the companion and the Cepheid f = 3.10+/-0.08%, giving an apparent magnitude mH = 8.47+/-0.15mag. The combination of interferometric and spectroscopic data have enabled the unique determination of the orbital elements: P = 1938.6+/-1.2 days, Tp = 2 443 616.1+/-7.3, a = 8.54+/-0.51mas, i = 124.7+/-1.8°, e = 0.190+/-0.013, ω = 228.7+/-1.6°, and Ω = 206.3+/-9.4°. We derive a minimal distance d ~ 691 pc, a minimum mass for both stars of 3.6 Msol, with a spectral type earlier than B5.5V for the companion star. Our measured flux ratio suggests that radial velocity detection of the companion using spectroscopy is within reach, and would provide an orbital parallax and model-free masses.

preprint2014arXiv

Reaching micro-arcsecond astrometry with long baseline optical interferometry; application to the GRAVITY instrument

A basic principle of long baseline interferometry is that an optical path difference (OPD) directly translates into an astrometric measurement. In the simplest case, the OPD is equal to the scalar product between the vector linking the two telescopes and the normalized vector pointing toward the star. However, a too simple interpretation of this scalar product leads to seemingly conflicting results, called here "the baseline paradox". For micro-arcsecond accuracy astrometry, we have to model in full the metrology measurement. It involves a complex system subject to many optical effects: from pure baseline errors to static, quasi-static and high order optical aberrations. The goal of this paper is to present the strategy used by the "General Relativity Analysis via VLT InTerferometrY" instrument (GRAVITY) to minimize the biases introduced by these defects. It is possible to give an analytical formula on how the baselines and tip-tilt errors affect the astrometric measurement. This formula depends on the limit-points of three type of baselines: the wide-angle baseline, the narrow-angle baseline, and the imaging baseline. We also, numerically, include non-common path higher-order aberrations, whose amplitude were measured during technical time at the Very Large Telescope Interferometer. We end by simulating the influence of high-order common-path aberrations due to atmospheric residuals calculated from a Monte-Carlo simulation tool for Adaptive optics systems. The result of this work is an error budget of the biases caused by the multiple optical imperfections, including optical dispersion. We show that the beam stabilization through both focal and pupil tracking is crucial to the GRAVITY system. Assuming the instrument pupil is stabilized at a 4 cm level on M1, and a field tracking below 0.2$λ/D$, we show that GRAVITY will be able to reach its objective of 10$μ$as accuracy.

preprint2014arXiv

Searching for visual companions of close Cepheids. VLT/NACO lucky imaging of Y~Oph, FF~Aql, X~Sgr, W~Sgr and $η$~Aql

Aims: High-resolution imaging in several photometric bands can provide color and astrometric information of the wide-orbit component of Cepheid stars. Such measurements are needed to understand the age and evolution of pulsating stars. In addition, binary Cepheids have the potential to provide direct and model-independent distances and masses. Methods: We used the NAOS-CONICA adaptive optics instrument (NACO) in the near-infrared to perform a deep search for wide components around the classical Cepheids, Y~Oph, FF~Aql, X~Sgr, W~Sgr, and $η$~Aql, within a field of view (FoV) of $1.7"\times 1.7"$ ($3.4"\times 3.4"$ for $η$~Aql). Results: We were able to reach contrast $ΔH = 5$-8\,mag and $ΔK_\mathrm{s} = 4$-7\,mag in the radius range $r > 0.2"$, which enabled us to constrain the presence of wide companions. For Y~Oph, FF~Aql, X~Sgr, W~Sgr, and $η$~Aql at $r > 0.2"$, we ruled out the presence of companions with a spectral type that is earlier than a B7V, A9V, A9V, A1V, and G5V star, respectively. For $0.1"< r < 0.2"$, no companions earlier than O9V, B3V, B4V, B2V, and B2V star, respectively, are detected. A component is detected close to $η$~Aql at projected separation $ρ= 654.7 \pm 0.9$\,mas and a position angle $PA = 92.8 \pm 0.1^\circ$. We estimated its dereddened apparent magnitude to be $m_H^0 = 9.34 \pm 0.04$ and derived a spectral type that ranges between an F1V and F6V star. Additional photometric and astrometric measurements are necessary to better constrain this star and check its physical association to the $η$~Aql system.

preprint2014arXiv

The long-period Galactic Cepheid RS Puppis - III. A geometric distance from HST polarimetric imaging of its light echoes

As one of the most luminous Cepheids in the Milky Way, the 41.5-day RS Puppis is an analog of the long-period Cepheids used to measure extragalactic distances. An accurate distance to this star would therefore help anchor the zero-point of the bright end of the period-luminosity relation. But, at a distance of about 2 kpc, RS Pup is too far away for measuring a direct trigonometric parallax with a precision of a few percent with existing instrumentation. RS Pup is unique in being surrounded by a reflection nebula, whose brightness varies as pulses of light from the Cepheid propagate outwards. We present new polarimetric imaging of the nebula obtained with HST/ACS. The derived map of the degree of linear polarization pL allows us to reconstruct the three-dimensional structure of the dust distribution. To retrieve the scattering angle from the pL value, we consider two different polarization models, one based on a Milky Way dust mixture and one assuming Rayleigh scattering. Considering the derived dust distribution in the nebula, we adjust a model of the phase lag of the photometric variations over selected nebular features to retrieve the distance of RS Pup. We obtain a distance of 1910 +/- 80 pc (4.2%), corresponding to a parallax of 0.524 +/- 0.022 mas. The agreement between the two polarization models we considered is good, but the final uncertainty is dominated by systematics in the adopted model parameters. The distance we obtain is consistent with existing measurements from the literature, but light echoes provide a distance estimate that is not subject to the same systematic uncertainties as other estimators (e.g. the Baade-Wesselink technique). RS Pup therefore provides an important fiducial for the calibration of systematic uncertainties of the long-period Cepheid distance scale.

preprint2013arXiv

A near-infrared interferometric survey of debris disc stars. III. First statistics based on 42 stars observed with CHARA/FLUOR

(Abridged) Dust is expected to be ubiquitous in extrasolar planetary systems owing to the dynamical activity of minor bodies. Inner dust populations are, however, still poorly known because of the high contrast and small angular separation with respect to their host star. We aim to determine the level of near-infrared exozodiacal dust emission around a sample of 42 nearby main sequence stars with spectral types ranging from A to K and to investigate its correlation with various stellar parameters and with the presence of cold dust belts. We use high-precision K-band visibilities obtained with the FLUOR interferometer on the shortest baseline of the CHARA array. The calibrated visibilities are compared with the expected visibility of the stellar photosphere to assess whether there is an additional, fully resolved circumstellar emission. Near-infrared circumstellar emission amounting to about 1% of the stellar flux is detected around 13 of our 42 target stars. Follow-up observations showed that one of them (eps Cep) is associated with a stellar companion, while another one was detected around what turned out to be a giant star (kap CrB). The remaining 11 excesses found around single main sequence stars are most probably associated with hot circumstellar dust, yielding an overall occurrence rate of 28+8-6% for our (biased) sample. We show that the occurrence rate of bright exozodiacal discs correlates with spectral type, K-band excesses being more frequent around A-type stars. It also correlates with the presence of detectable far-infrared excess emission in the case of solar-type stars. This study provides new insight into the phenomenon of bright exozodiacal discs, showing that hot dust populations are probably linked to outer dust reservoirs in the case of solar-type stars. For A-type stars, no clear conclusion can be made regarding the origin of the detected near-infrared excesses.

preprint2013arXiv

Extended envelopes around Galactic Cepheids IV. T Monocerotis and X Sagittarii from mid-infrared interferometry with VLTI/MIDI

Aims. We study the close environment of nearby Cepheids using high spatial resolution observations in the mid-infrared with the VLTI/MIDI instrument, a two-beam interferometric recombiner. Methods. We obtained spectra and visibilities for the classical Cepheids X Sgr and T Mon. We fitted the MIDI measurements, supplemented by B, V, J, H, K literature photometry, with the numerical transfer code DUSTY to determine the dust shell parameters. We used a typical dust composition for circumstellar environments. Results. We detect an extended dusty environment in the spectra and visibilities for both stars, although T Mon might suffer from thermal background contamination. We attribute this to the presence of a circumstellar envelope (CSE) surrounding the Cepheids. This is optically thin for X Sgr (tau(0.55microns) = 0.008), while it appears to be thicker for T Mon (tau(0.55micron) = 0.15). They are located at about 15-20 stellar radii. Following our previous work, we derived a likely period-excess relation in the VISIR PAH1 filter, f(8.6micron)[%]= 0.81(+/-0.04)P[day]. We argue that the impact of CSEs on the mid-IR period-luminosity (P-L) relation cannot be negligible because they can bias the Cepheid brightness by up to about 30 %. For the K-band P-L relation, the CSE contribution seems to be lower (< 5 %), but the sample needs to be enlarged to firmly conclude that the impact of the CSEs is negligible in this band.

preprint2013arXiv

Flares and variability from Sagittarius A*: five nights of simultaneous multi-wavelength observations

Aims. We report on simultaneous observations and modeling of mid-infrared (MIR), near-infrared (NIR), and submillimeter (submm) emission of the source Sgr A* associated with the supermassive black hole at the center of our Galaxy. Our goal was to monitor the activity of Sgr A* at different wavelengths in order to constrain the emitting processes and gain insight into the nature of the close environment of Sgr A*. Methods. We used the MIR instrument VISIR in the BURST imaging mode, the adaptive optics assisted NIR camera NACO, and the sub-mm antenna APEX to monitor Sgr A* over several nights in July 2007. Results. The observations reveal remarkable variability in the NIR and sub-mm during the five nights of observation. No source was detected in the MIR, but we derived the lowest upper limit for a flare at 8.59 microns (22.4 mJy with A_8.59mu = 1.6+/- 0.5). This observational constraint makes us discard the observed NIR emission as coming from a thermal component emitting at sub-mm frequencies. Moreover, comparison of the sub-mm and NIR variability shows that the highest NIR fluxes (flares) are coincident with the lowest sub-mm levels of our five-night campaign involving three flares. We explain this behavior by a loss of electrons to the system and/or by a decrease in the magnetic field, as might conceivably occur in scenarios involving fast outflows and/or magnetic reconnection.

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.

preprint2012arXiv

Fundamental properties of the Population II fiducial stars HD 122563 and Gmb 1830 from CHARA interferometric observations

We have determined the angular diameters of two metal-poor stars, HD 122563 and Gmb 1830, using CHARA and Palomar Testbed Interferometer observations. For the giant star HD 122563, we derive an angular diameter theta_3D = 0.940 +- 0.011 milliarcseconds (mas) using limb-darkening from 3D convection simulations and for the dwarf star Gmb 1830 (HD 103095) we obtain a 1D limb-darkened angular diameter theta_1D = 0.679 +- 0.007 mas. Coupling the angular diameters with photometry yields effective temperatures with precisions better than 55 K (Teff = 4598 +- 41 K and 4818 +- 54 K --- for the giant and the dwarf star, respectively). Including their distances results in very well-determined luminosities and radii (L = 230 +- 6 L_sun, R = 23.9 +- 1.9 R_sun and L = 0.213 +- 0.002 L_sun, R = 0.664 +- 0.015 R_sun, respectively). We used the CESAM2k stellar structure and evolution code in order to produce models that fit the observational data. We found values of the mixing-length parameter alpha (which describes 1D convection) that depend on the mass of the star. The masses were determined from the models with precisions of <3% and with the well-measured radii excellent constraints on the surface gravity are obtained (log g = 1.60 +- 0.04, 4.59 +- 0.02, respectively). The very small errors on both log g and Teff provide stringent constraints for spectroscopic analyses given the sensitivity of abundances to both of these values. The precise determination of Teff for the two stars brings into question the photometric scales for metal-poor stars.

preprint2012arXiv

Interferometry, spectroscopy and astrometry of the bright eclipsing system Delta Velorum

The bright southern star Delta Vel is a multiple system comprising at least three stars. Its brightest component, Delta Vel A, was identified in 2000 as one of the brightest eclipsing system in the sky. Its eclipses are easily observable with the unaided eye, a remarkable property shared only by Algol, Beta Aur, Alpha CrB and Psi Cen. We determined dynamical masses from a combination of spectroscopy, high-precision astrometry of the orbits of Aab-B and Aa-Ab using adaptive optics (VLT/NACO) and optical interferometry (VLTI/AMBER). The main eclipsing component is a pair of A-type stars in rapid rotation. We modeled the photometric and radial velocity measurements of the eclipsing pair Aa-Ab using a self consistent method based on physical parameters (mass, radius, luminosity, rotational velocity). From this modeling, we derive the fundamental parameters of the eclipsing stars with a typical accuracy of 1%. We find that they have similar masses, respectively 2.43 +/- 0.02 and 2.27 +/- 0.02 Msun. The physical parameters of the tertiary component (Delta Vel B) are also derived, although to a lower accuracy, as well as the parallax of the system, 39.8 +/- 0.4 mas. This value is in satisfactory agreement (-1.2 sigma) with the Hipparcos parallax of the system (pi_Hip =40.5 +/- 0.4 mas).

preprint2012arXiv

Mean angular diameters, distances and pulsation modes of the classical Cepheids FF Aql and T Vul - CHARA/FLUOR near-infrared interferometric observations

We report the first angular diameter measurements of two classical Cepheids, FF Aql and T Vul, that we have obtained with the FLUOR instrument installed at the CHARA interferometric array. We obtain average limb-darkened angular diameters of θ_LD = 0.878 +/- 0.013 mas and θ_LD = 0.629 +/- 0.013 mas, respectively for FF Aql and T Vul. Combining these angular diameters with the HST-FGS trigonometric parallaxes leads to linear radii R = 33.6 +/- 2.2 Rsol and R = 35.6 +/- 4.4 Rsol, respectively. The comparison with empirical and theoretical Period-Radius relations leads to the conclusion that these Cepheids are pulsating in their fundamental mode. The knowledge of the pulsation mode is of prime importance to calibrate the Period-Luminosity relation with a uniform sample of fundamental mode Cepheids.

preprint2012arXiv

The expanding dusty bipolar nebula around the nova V1280 Sco

V1280 Sco is one of the slowest dust-forming nova ever historically observed. We performed multi-epoch high-spatial resolution observations of the circumstellar dusty environment of V1280 Sco to investigate the level of asymmetry of the ejecta We observed V1280 Sco in 2009, 2010 and 2011 using unprecedented high angular resolution techniques. We used the NACO/VLT adaptive optics system in the J, H and K bands, together with contemporaneous VISIR/VLT mid-IR imaging that resolved the dust envelope of V1280 Sco, and SINFONI/VLT observations secured in 2011. We report the discovery of a dusty hourglass-shaped bipolar nebula. The apparent size of the nebula increased from 0.30" x 0.17" in July 2009 to 0.64" x 0.42" in July 2011. The aspect ratio suggests that the source is seen at high inclination. The central source shines efficiently in the K band and represents more than 56+/-5% of the total flux in 2009, and 87+/-6% in 2011. A mean expansion rate of 0.39+/-0.03 mas per day is inferred from the VISIR observations in the direction of the major axis, which represents a projected upper limit. Assuming that the dust shell expands in that direction as fast as the low-excitation slow ejecta detected in spectroscopy, this yields a lower limit distance to V1280 Sco of 1kpc; however, the systematic errors remain large due to the complex shape and velocity field of the dusty ejecta. The dust seems to reside essentially in the polar caps and no infrared flux is detected in the equatorial regions in the latest dataset. This may imply that the mass-loss was dominantly polar.

preprint2012arXiv

The long-period Galactic Cepheid RS Puppis - II. 3D structure and mass of the nebula from VLT/FORS polarimetry

The long-period Cepheid RS Pup is surrounded by a large dusty nebula reflecting the light from the central star. Due to the changing luminosity of the central source, light echoes propagate into the nebula. This remarkable phenomenon was the subject of Paper I.The origin and physical properties of the nebula are however uncertain: it may have been created through mass loss from the star itself, or it could be the remnant of a pre-existing interstellar cloud. Our goal is to determine the 3D structure of the nebula, and estimate its mass. Knowing the geometrical shape of the nebula will also allow us to retrieve the distance of RS Pup in an unambiguous manner using a model of its light echoes (in a forthcoming work). The scattering angle of the Cepheid light in the circumstellar nebula can be recovered from its degree of linear polarization. We thus observed the nebula surrounding RS Pup using the polarimetric imaging mode of the VLT/FORS instrument, and obtained a map of the degree and position angle of linear polarization. From our FORS observations, we derive a 3D map of the distribution of the dust, whose overall geometry is an irregular and thin layer. The nebula does not present a well-defined symmetry. Using a simple model, we derive a total dust mass of M(dust) = 2.9 +/- 0.9 Msun for the dust within 1.8 arcmin of the Cepheid. This translates into a total mass of M(gas+dust) = 290 +/- 120 Msun, assuming a dust-to-gas ratio of 1.0 +/- 0.3 %. The high mass of the dusty nebula excludes that it was created by mass-loss from the star. However, the thinness nebula is an indication that the Cepheid participated to its shaping, e.g. through its radiation pressure or stellar wind. RS Pup therefore appears as a regular long-period Cepheid located in an exceptionally dense interstellar environment.

preprint2012arXiv

Three-dimensional interferometric, spectrometric, and planetary views of Procyon

We used a new realistic 3D radiative-hydrodynamical model atmosphere of Procyon generated with the Stagger Code and synthetic spectra computed with the radiative transfer code Optim3D to re-analyze interferometric and spectroscopic data from the optical to the infrared of Procyon. We compute intensity maps in two optical filters centered at 500 and 800 nm (MARK III) and one infrared filter centered at 2200 nm (VINCI). We constructed stellar disk images accounting for the center-to-limb variations and used them to derive visibility amplitudes and closure phases. We provide 3D limb-darkening coefficients in the optical as well as in the infrared. We show that visibility curves and closure phases show clear deviations from circular symmetry from the 3rd lobe on. These deviations are detectable with current interferometers using closure phases. We derive new angular diameters at different wavelengths with two independent methods based on 3D simulations. We find a diameter_Vinci = 5.390 \pm 0.03 mas that this is confirmed by an independent asteroseismic estimation. The resulting Teff is 6591 K, which is consistent with the infrared flux method determinations. We find also a value of the surface gravity log g = 4.01 \pm 0.03 that is larger by 0.05 dex from literature values. Spectrophotometric comparisons with observations provide very good agreement with the spectral energy distribution and photometric colors, allowing us to conclude that the thermal gradient of the simulation matches fairly well Procyon. Finally, we show that the granulation pattern of a planet hosting Procyon-like star has a non-negligible impact on the detection of hot Jupiters in the infrared using interferometry closure phases. It is then crucial to have a comprehensive knowledge of the host star to directly detect and characterize hot Jupiters. In this respect, RHD simulations are very important to reach this aim.

preprint2011arXiv

EChO - Exoplanet Characterisation Observatory

A dedicated mission to investigate exoplanetary atmospheres represents a major milestone in our quest to understand our place in the universe by placing our Solar System in context and by addressing the suitability of planets for the presence of life. EChO -the Exoplanet Characterisation Observatory- is a mission concept specifically geared for this purpose. EChO will provide simultaneous, multi-wavelength spectroscopic observations on a stable platform that will allow very long exposures. EChO will build on observations by Hubble, Spitzer and groundbased telescopes, which discovered the first molecules and atoms in exoplanetary atmospheres. EChO will simultaneously observe a broad enough spectral region -from the visible to the mid-IR- to constrain from one single spectrum the temperature structure of the atmosphere and the abundances of the major molecular species. The spectral range and resolution are tailored to separate bands belonging to up to 30 molecules to retrieve the composition and temperature structure of planetary atmospheres. The target list for EChO includes planets ranging from Jupiter-sized with equilibrium temperatures Teq up to 2000 K, to those of a few Earth masses, with Teq ~300 K. We have baselined a dispersive spectrograph design covering continuously the 0.4-16 micron spectral range in 6 channels (1 in the VIS, 5 in the IR), which allows the spectral resolution to be adapted from several tens to several hundreds, depending on the target brightness. The instrument will be mounted behind a 1.5 m class telescope, passively cooled to 50 K, with the instrument structure and optics passively cooled to ~45 K. EChO will be placed in a grand halo orbit around L2. We have also undertaken a first-order cost and development plan analysis and find that EChO is easily compatible with the ESA M-class mission framework.

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

The Baade-Wesselink p-factor applicable to LMC Cepheids

Context. Recent observations of LMC Cepheids bring new constraints on the slope of the period-projection factor relation (hereafter Pp relation) that is currently used in the Baade-Wesselink (hereafter BW) method of distance determination. The discrepancy between observations and theoretical analysis is particularly significant for short period Cepheids Aims. We investigate three physical effects that might possibly explain this discrepancy: (1) the spectroscopic S/N that is systematically lower for LMC Cepheids (around 10) compared to Galactic ones (up to 300), (2) the impact of the metallicity on the dynamical structure of LMC Cepheids, and (3) the combination of infrared photometry/interferometry with optical spectroscopy. Methods. To study the S/N we use a very simple toy model of Cepheids. The impact of metallicity on the projection factor is based on the hydrodynamical model of delta Cep already described in previous studies. This model is also used to derive the position of the optical versus infrared photospheric layers. Results. We find no significant effect of S/N, metallicity, and optical-versus-infrared observations on the Pp relation. Conclusions. The Pp relation of Cepheids in the LMC does not differ from the Galactic relation. This allows its universal application to determine distances to extragalactic Cepheids via BW analysis.

preprint2011arXiv

The diameter of the CoRoT target HD 49933. Combining the 3D limb darkening, asteroseismology, and interferometry

Context. The interpretation of stellar pulsations in terms of internal structure depends on the knowledge of the fundamental stellar parameters. Long-base interferometers permit us to determine very accurate stellar radii, which are independent constraints for stellar models that help us to locate the star in the HR diagram. Aims: Using a direct interferometric determination of the angular diameter and advanced three-dimensional (3D) modeling, we derive the radius of the CoRoT target HD 49933 and reduce the global stellar parameter space compatible with seismic data. Methods: The VEGA/CHARA spectro-interferometer is used to measure the angular diameter of the star. A 3D radiative hydrodynamical simulation of the surface is performed to compute the limb darkening and derive a reliable diameter from visibility curves. The other fundamental stellar parameters (mass, age, and Teff) are found by fitting the large and small p-mode frequency separations using a stellar evolution model that includes microscopic diffusion. Results: We obtain a limb-darkened angular diameter of θLD = 0.445 \pm 0.012 mas. With the Hipparcos parallax, we obtain a radius of R = 1.42 \pm 0.04 Rsun. The corresponding stellar evolution model that fits both large and small frequency separations has a mass of 1.20 \pm 0.08 Msun and an age of 2.7 Gy. The atmospheric parameters are Teff = 6640 \pm 100 K, log g = 4.21 \pm 0.14, and [Fe/H] = -0.38.

preprint2011arXiv

The nearby eclipsing stellar system delta Velorum III. Self-consistent fundamental parameters and distance

The triple stellar system delta Vel (composed of two A-type and one F-type main sequence stars) is particularly interesting as it contains one of the nearest and brightest eclipsing binaries. It therefore presents a unique opportunity to determine independently the physical properties of the three components of the system, as well as its distance. We aim at determining the fundamental parameters (masses, radii, luminosities, rotational velocities) of the three components of delta Vel, as well as the parallax of the system, independently from the existing Hipparcos} measurement. We determined dynamical masses from high-precision astrometry of the orbits of Aab-B and Aa-Ab using adaptive optics (VLT/NACO) and optical interferometry (VLTI/AMBER). The main component is an eclipsing binary composed of two early A-type stars in rapid rotation. We modeled the photometric and radial velocity measurements of the eclipsing pair Aa-Ab using a self consistent method based on physical parameters (mass, radius, luminosity, rotational velocity). From our self-consistent modeling of the primary and secondary components of the delta Vel A eclipsing pair, we derive their fundamental parameters with a typical accuracy of 1%. We find that they have similar masses, respectively 2.43+/-0.02Msol and 2.27+/-0.02Msol. The physical parameters of the tertiary component (delta Vel B) are also estimated, although to a lower accuracy. We obtain a parallax of 39.8+/-0.4mas for the system, in satisfactory agreement (-1.2 sigma) with the Hipparcos value (40.5+/-0.4mas). The physical parameters we derive represent a consistent set of constraints for the evolutionary modeling of this system. The agreement of the parallax we measure with the Hipparcos value to a 1% accuracy is also an interesting confirmation of the true accuracy of these two independent measurements.

preprint2011arXiv

The nearby eclipsing stellar system delta Velorum. II. First reliable orbit for the eclipsing pair

Context. The nearby multiple system delta Velorum contains a widely detached eclipsing binary and a third component. Aims. We take advantage of this system offering the opportunity to determine the set of fundamental parameters (masses, luminosities, and radii) of three coeval stars with sufficient precision to test models of stellar evolution. Methods. Extensive high-resolution spectroscopy is analyzed by the broadening function technique to provide the first spectroscopic orbit of the eclipsing pair. Simultaneous analysis of the spectroscopic data and the SMEI satellite light curve is performed to provide astrophysical parameters for the components. We use a modified Roche model assuming an eccentric orbit and asynchronous rotation. Results. The observations show that components of the eclipsing pair rotate at about two-thirds of the break-up velocity, which excludes any chemical peculiarity and results in a non-uniform surface brightness. Although the inner orbit is eccentric, no apsidal motion is seen during the SMEI photometric observations. For the inner orbit, the orbital parameters are eccentricity e = 0.29, longitude of the periastron passage omega = 109 degrees, and inclination 89.0 degrees. Conclusions. Component's masses MAa = 2.53+/-0.11 Msun, MAb = 2.37+/-0.10 Msun and MB = 1.5 Msun combined with inferred radii of Aa and Ab components indicate that the eclipsing pair has already left the Main sequence and the estimated age of the system is about 400 Myr.

preprint2010arXiv

Surface convection and red giants radii measurements

The phenomenological models of convection use characteristic length scales they do not determine but that are chosen to fit solar or stellar observations. We investigate if changes of these length scales are required between the Sun and low mass stars on the red giant branch (RGB). The question is addressed jointly in the frameworks of the mixing length theory and of the full spectrum of turbulence model. For both models, the convective length scale is assumed to be a fixed fraction of the local pressure scale height. We use constraints coming from the observed effective temperatures and linear radii independently. We rely on a sample of 38 nearby giants and subgiants for which surface temperatures and luminosities are known accurately and the radii are determined through interferometry to better than 10%. For the few cases where the stellar masses were determined by asteroseismological measurements, we computed dedicated models. First we calibrate the solar models. Then, with the same physics, we compute RGB models for masses between 0.9 Mo and 2.5 Mo and metallicities ranging from $\rm [Fe/H]=-0.34$ to solar. The evolution is followed up to 1000 Lo. A special attention is given to the opacities and to the non grey atmosphere models used as boundary conditions for which the model of convection is the same as in the interior. We find that for both the mixing length theory and the full spectrum of turbulence model the characteristic solar length scale for convection has to be slightly reduced to fit the lower edge of the observed RGB. The corresponding models also better match the expected mass distribution on the RGB and are in better agreement to the seismic constraints. These results are robust whether effective temperatures determined spectroscopically or radii determined interferometrically are used.

preprint2009arXiv

Asteroseismology and interferometry of the red giant star epsilon Oph

The GIII red giant star epsilon Oph has been found to exhibit several modes of oscillation by the MOST mission. We interpret the observed frequencies of oscillation in terms of theoretical radial p-mode frequencies of stellar models. Evolutionary models of this star, in both shell H-burning and core He-burning phases of evolution, are constructed using as constraints a combination of measurements from classical ground-based observations (for luminosity, temperature, and chemical composition) and seismic observations from MOST. Radial frequencies of models in either evolutionary phase can reproduce the observed frequency spectrum of epsilon Oph almost equally well. The best-fit models indicate a mass in the range of 1.85 +/- 0.05 Msun with radius of 10.55 +/- 0.15 Rsun. We also obtain an independent estimate of the radius of epsilon Oph using high accuracy interferometric observations in the infrared K' band, using the CHARA/FLUOR instrument. The measured limb darkened disk angular diameter of epsilon Oph is 2.961 +/- 0.007 mas. Together with the Hipparcos parallax, this translates into a photospheric radius of 10.39 +/- 0.07 Rsun. The radius obtained from the asteroseismic analysis matches the interferometric value quite closely even though the radius was not constrained during the modelling.

preprint2009arXiv

High resolution spectroscopy for Cepheids distance determination. V. Impact of the cross-correlation method on the p-factor and the gamma-velocities

The cross correlation method (hereafter CC) is widely used to derive the radial velocity curve of Cepheids when the signal to noise of the spectra is low. However, if it is used with the wrong projection factor, it might introduce some biases in the Baade-Wesselink (hereafter BW) methods of determining the distance of Cepheids. In addition, it might affect the average value of the radial velocity curve (or gamma-velocity) important for Galactic structure studies. We aim to derive a period-projection factor relation (hereafter Pp) appropriate to be used together with the CC method. Moreover, we investigate whether the CC method can explain the misunderstood previous calculation of the K-term of Cepheids. We observed eight galactic Cepheids with the HARPS spectrograph. For each star, we derive an interpolated CC radial velocity curve using the HARPS pipeline. The amplitudes of these curves are used to determine the correction to be applied to the semi-theoretical projection factor derived in Nardetto et al. (2007). Their average value (or gamma-velocity) are also compared to the center-of-mass velocities derived in Nardetto et al. (2008). The correction in amplitudes allows us to derive a new Pp relation: p = [-0.08+-0.05] log P +[1.31+-0.06]. We also find a negligible wavelength dependence (over the optical range) of the Pp relation. We finally show that the gamma-velocity derived from the CC method is systematically blue-shifted by about 1.0 +- 0.2km/s compared to the center-of-mass velocity of the star. An additional blue-shift of 1.0km/s is thus needed to totally explain the previous calculation of the K-term of Cepheids (around 2km/s). The new Pp relation we derived is a solid tool for the distance scale calibration (abridged).

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.

preprint2009arXiv

The radius and effective temperature of the binary Ap star beta CrB from CHARA/FLUOR and VLT/NACO observations

The prospects for using asteroseismology of rapidly oscillating Ap (roAp) stars are hampered by the large uncertainty in fundamental stellar parameters. Results in the literature for the effective temperature (Teff) often span a range of 1000 K. Our goal is to reduce systematic errors and improve the Teff calibration of Ap stars based on new interferometric measurements. We obtained long-baseline interferometric observations of beta CrB using the CHARA/FLUOR instrument. To disentangle the flux contributions of the two components of this binary star, we obtained VLT/NACO adaptive optics images. We determined limb darkened angular diameters of 0.699+-0.017 mas for beta CrB A (from interferometry) and 0.415+-0.017 mas for beta CrB B (from surface brightness- color relations), corresponding to radii of 2.63+-0.09 Rsun (3.4 percent uncertainty) and 1.56+-0.07 Rsun (4.5 percent). The combined bolometric flux of the A and B components was determined from satellite UV data, spectrophotometry in the visible and broadband data in the infrared. The flux from the B component constitutes 16+-4 percent of the total flux and was determined by fitting an ATLAS9 model atmosphere to the broad-band NACO J and K magnitudes. Combining the flux of the A component with its measured angular diameter, we determine the effective temperature Teff(A) = 7980+-180 K (2.3 percent). Our new interferometric and imaging data enable a nearly model-independent determination of the effective temperature of beta CrB A. Including our recent study of alpha Cir, we now have direct Teff measurements of two of the brightest roAp stars, providing a strong benchmark for an improved calibration of the Teff scale for Ap stars. This will support the use of potentially strong constraints imposed by asteroseismic studies of roAp stars.

preprint2003arXiv

The diameter and evolutionary state of Procyon A

We report the angular diameter measurement obtained with the VINCI/VLTI instrument on the nearby star Procyon A (alpha CMi A, F5IV-V), at a relative precision of +/- 0.9%. We obtain a uniform disk angular diameter in the K band of UD = 5.376 +/- 0.047 mas and a limb darkened value of LD = 5.448 +/- 0.053 mas. Together with the Hipparcos parallax, this gives a linear diameter of 2.048 +/- 0.025 D_sun. We use this result in combination with spectroscopic, photometric and asteroseismic constraints to model this star with the CESAM code. One set of modeling parameters that reproduces the observations within their error bars are an age of 2314 Myr, an initial helium mass fraction Y_i = 0.301 and an initial mass ratio of heavy elements to hydrogen Z/X_i = 0.0314. We also computed the adiabatic oscillation spectrum of our model of Procyon A, giving a mean large frequency separation of Delta nu_0 = 54.7 microHz. This value is in agreement with the seismic observations by Martic et al. (1999, 2001). The interferometric diameter and the asteroseismic large frequency spacing together suggest a mass closer to 1.4 M_sun than to 1.5 M_sun. We conclude that Procyon is currently ending its life on the main sequence, as its luminosity class indicates.

preprint2003arXiv

The diameters of Alpha Centauri A and B - A comparison of the asteroseismic and VINCI/VLTI views

We compare the first direct angular diameter measurements obtained on our closest stellar neighbour, Alpha Centauri, to recent model diameters constrained by asteroseismic observations. Using the VINCI instrument installed at ESO's VLT Interferometer (VLTI), the angular diameters of the two main components of the system, Alpha Cen A and B, were measured with a relative precision of 0.2% and 0.6%, respectively. Particular care has been taken in the calibration of these measurements, considering that VINCI is estimating the fringe visibility using a broadband K filter. We obtain uniform disk angular diameters for Alpha Cen A and B of UD[A] = 8.314 +/- 0.016 mas and UD[B] = 5.856 +/- 0.027 mas, and limb darkened angular diameters of LD[A] = 8.511 +/- 0.020 mas and LD[B] = 6.001 +/- 0.034 mas. Combining these values with the parallax from Soderhjelm (1999), we derive linear diameters of D[A] = 1.224 +/- 0.003 Dsun and D[B] = 0.863 +/- 0.005 Dsun. These values are compatible with the masses published by Thevenin et al.(2002) for both stars.

preprint2003arXiv

The spinning-top Be star Achernar from VLTI-VINCI

We report here the first observations of a rapidly rotating Be star, $α$ Eridani, using Earth-rotation synthesis on the Very Large Telescope (VLT) Interferometer. Our measures correspond to a $2a/2b = 1.56\pm0.05$ apparent oblate star, $2a$ and $2b$ being the equivalent uniform disc angular diameters in the equatorial and polar direction. Considering the presence of a circumstellar envelope (CSE) we argue that our measurement corresponds to a truly distorted star since $α$ Eridani exhibited negligible H$α$ emission during the interferometric observations. In this framework we conclude that the commonly adopted Roche approximation (uniform rotation and centrally condensed mass) should not apply to $α$ Eridani. This result opens new perspectives to basic astrophysical problems, such as rotationally enhanced mass loss and internal angular momentum distribution. In addition to its intimate relation with magnetism and pulsation, rapid rotation thus provides a key to the Be phenomenon: one of the outstanding non-resolved problems in stellar physics.