Issue |
A&A
Volume 651, July 2021
|
|
---|---|---|
Article Number | A92 | |
Number of page(s) | 13 | |
Section | Interstellar and circumstellar matter | |
DOI | https://doi.org/10.1051/0004-6361/202140626 | |
Published online | 21 July 2021 |
Mid-infrared circumstellar emission of the long-period Cepheid ℓ Carinae resolved with VLTI/MATISSE★
1
Université Côte d’Azur, Observatoire de la Côte d’Azur, CNRS, Laboratoire Lagrange,
France,
e-mail: vincent.hocde@oca.eu
2
European Southern Observatory,
Karl-Schwarzschild-Str. 2,
85748
Garching,
Germany
3
Max-Planck-Institut für Radioastronomie,
Auf dem Hügel 69,
53121
Bonn,
Germany
4
Leiden Observatory, Leiden University,
Niels Bohrweg 2,
2333
CA
Leiden,
The Netherlands
5
LESIA, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, Université de Paris,
5 place Jules Janssen,
92195
Meudon,
France
6
Nicolaus Copernicus Astronomical Centre, Polish Academy of Sciences,
Bartycka 18,
00-716
Warszawa,
Poland
7
Unidad Mixta Internacional Franco-Chilena de Astronomía (CNRS UMI 3386), Departamento de Astronomía, Universidad de Chile,
Camino El Observatorio 1515, Las Condes,
Santiago,
Chile
8
Departamento de Astronomía, Universidad de Concepción,
Casilla 160-C,
Concepción,
Chile
9
Univ. Grenoble Alpes, CNRS, IPAG,
38000,
Grenoble,
France
10
NASA Goddard Space Flight Center, Astrophysics Division,
Greenbelt,
MD
20771,
USA
11
Anton Pannekoek Institute for Astronomy, University of Amsterdam,
Science Park 904,
1090
GE Amsterdam,
The Netherlands
12
Department of Astrophysics, University of Vienna,
Türkenschanzstrasse 17,
Vienna,
Austria
13
AIM, CEA, CNRS, Université Paris-Saclay, Université Paris Diderot,
Sorbonne Paris Cité,
91191
Gif-sur-Yvette,
France
14
Institute for Mathematics, Astrophysics and Particle Physics, Radboud University,
PO Box 9010,
MC 62, 6500
GL Nijmegen,
the Netherlands
15
SRON Netherlands Institute for Space Research,
Sorbonnelaan 2,
3584
CA Utrecht,
The Netherlands
16
Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel,
Leibnizstraße 15,
24118,
Kiel,
Germany
17
Konkoly Observatory, Research Centre for Astronomy and Earth Sciences, Eötvös Loránd Research Network (ELKH),
Konkoly-Thege Miklós út 15-17,
1121
Budapest,
Hungary
18
CSFK Lendület Near-Field Cosmology Research Group,
Budapest,
Hungary
19
Max Planck Institute for Astronomy,
Königstuhl 17,
69117
Heidelberg,
Germany
20
Institute for Astrophysics, University of Vienna,
Türkenschanzstrasse 17,
1180
Vienna,
Austria
21
I. Physikalisches Institut, Universität zu Köln,
Zülpicher Str. 77,
50937,
Köln,
Germany
22
Zselic Park of Stars,
064/2 hrsz.,
7477
Zselickisfalud,
Hungary
23
Sydney Institute for Astronomy (SIfA), School of Physics, The University of Sydney,
NSW
2006,
Australia
24
NOVA Optical IR Instrumentation Group at ASTRON,
Oude Hoogeveensedijk 4 7991 PD Dwingeloo,
The Netherlands
25
ASTRON (Netherlands),
Oude Hoogeveensedijk 4, 7991 PD Dwingeloo,
The Netherlands
Received:
22
February
2021
Accepted:
24
March
2021
Context. The nature of circumstellar envelopes (CSEs) around Cepheids is a matter of ongoing debate. The physical origin of their infrared (IR) excess could be shown to either be made up of a shell of ionized gas, a dust envelope, or a combination of both.
Aims. This study is aimed at constraining the geometry and the IR excess of the environment of the bright long-period Cepheid ℓ Car (P = 35.5 days) at mid-IR wavelengths in order to understand its physical nature.
Methods. We first used photometric observations in various bands (from the visible domain to the infrared) and Spitzer Space Telescope spectroscopy to constrain the IR excess of ℓ Car. Then we analyzed the VLTI/MATISSE measurements at a specific phase of observation in order to determine the flux contribution as well as the size and shape of the environment of the star in the L band. Finally, we tested the hypothesis of a shell of ionized gas in order to model the IR excess.
Results. We report the first detection in the L band of a centro-symmetric extended emission around ℓ Car, of about 1.7 R⋆ in full width at half maximum, producing an excess of about 7.0% in this band.This latter value is used to calibrate the IR excess found when comparing the photometric observations in various bands and quasi-static atmosphere models. In the N band, there is no clear evidence for dust emission from VLTI/MATISSE correlated flux and Spitzer data. On the other side, the modeled shell of ionized gas implies a more compact CSE (1.13 ± 0.02 R⋆) that is also fainter (IR excess of 1% in the L band).
Conclusions. We provide new evidence supporting a compact CSE for ℓ Car and we demonstrate the capabilities of VLTI/MATISSE for determining common properties of CSEs. While the compact CSE of ℓ Car is likely to be of a gaseous nature, the tested model of a shell of ionized gas is not able to simultaneously reproduce the IR excess and the interferometric observations. Further Galactic Cepheid observations with VLTI/MATISSE are necessary for determining the properties of CSEs, which may also depend on both the pulsation period and the evolutionary state of the stars.
Key words: instrumentation: interferometers / circumstellar matter / infrared: stars / stars: variables: Cepheids / stars: atmospheres
© V. Hocdé et al. 2021
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Circumstellar envelopes (CSEs) around Cepheids have been spatially resolved by long-baseline interferometry in the K band with the Very Large Telescope Interferometer (VLTI) and the Center for High Angular Resolution Astronomy (CHARA; Kervella et al. 2006; Mérand et al. 2006). These structures have been detected around four Cepheids, including ℓ Car, in the N band with VLTI/VISIR and VLTI/MIDI (Kervella et al. 2009; Gallenne et al. 2013). In the K band, the diameter of the envelope appears to be at least two stellar radii and the flux contribution is up to 5% of the continuum (Mérand et al. 2007). Both the size and flux contribution are expected to be larger in the N band, as demonstrated by Kervella et al. (2009) for ℓ Car, Gallenne et al. (2013) for X Sgr and T Mon, and Gallenne et al. (2021) for other Cepheids. Moreover, a CSE has also been discovered in the visible domain with the CHARA/VEGA instrument around δ Cep (Nardetto et al. 2016). The systematic presence of a CSE is still a matter of debate. While some studies have found no significant observational evidence for circumstellar dust envelopes in a large number of Cepheids (Schmidt 2015; Groenewegen 2020), Gallenne et al. (2021) found a significant infrared (IR) excess attributed to a CSE for 10 out of 45 Cepheids.
Recently, Hocdé et al. (2020b; hereafter, Paper I) used an analytical model of free-free and bound-free emission from a thin shell of ionized gas to explain the near and mid-IR excess of Cepheids. They found a typical radius for this shell of ionized gas of about Rshell = 1.15 R⋆.
This shell of ionized gas could be due to periodic shocks occurring in both the atmosphere and the chromosphere, which heat up and ionize the gas. Using VLT/UVES data Hocdé et al. (2020a), found a radius for the chromosphere of at least Rchromo = 1.5 R⋆ in long-period Cepheids. In addition, they found a motionless Hα feature in UVES high-resolution spectra obtained for several long-period Cepheids, including ℓ Car. This absorption feature was attributed to a static CSE surrounding the chromosphere above at least 1.5 R⋆, and reported by various authors around ℓ Car (Rodgers & Bell 1968; Bohm-Vitense & Love 1994; Nardetto et al. 2008).
Determining the IR excess and the size of the CSE is a key for understanding the physical processes at play. In this paper, we study the long-period Cepheid ℓ Car. Our aim is to determine its IR excess from photometric measurements in various bands and Spitzer spectroscopy, while inferring the size of the CSE and its flux contribution thanks to the unique capabilities provided by the Multi AperTure mid-Infrared SpectroScopic Experiment (VLTI/MATISSE; Lopez et al. 2014; Allouche et al. 2016; Robbe-Dubois et al. 2018) in the L (2.8–4.0 μm), M (4.5–5 μm), and N bands (8–13 μm).
We first reconstruct the IR excess in Sect. 2. Then we present the data reduction and calibration process of VLTI/MATISSE data in Sect. 3. We use a simple model of the CSE to reproduce the visibility measurement of VLTI/MATISSE in Sect. 4. In Sect. 5, we discuss the physical origin of this envelope and we present a model of a shell of ionized gas. We summarize our conclusions in Sect. 6.
2 Deriving the infrared excess of ℓ Car
2.1 Near-IR excess modeling with SPIPS
Due to their intrinsic variability, the photospheres of the Cepheids are difficult to model along the pulsation cycle. However, it is an essential prerequisite for deriving both the IR excess in a given photometric band and the expected angular stellar diameter at a specific phase of interferometric observations (see Sect. 3). In order to model the photosphere, we used the SpectroPhoto-Interferometric modeling of Pulsating Stars (SPIPS), which is a model-based parallax-of-pulsation code that includes photometric, interferometric, effective temperature, and radial velocity measurements in a robust model-fitting process (Mérand et al. 2015).
The SPIPS uses a grid of ATLAS9 atmospheric models1 (Castelli & Kurucz 2003) with solar metallicity and a standard turbulent velocity of 2 km s−1. SPIPS was already extensively described and used in several studies (Mérand et al. 2015; Breitfelder et al. 2016; Kervella et al. 2017; Gallenne et al. 2017, 2021; Javanmardi et al. 2021), and Paper I. SPIPS also takes into account the possible presence of an IR excess by fitting an ad hoc analytic power law IRex (see green curve in Fig. 1a). We note that Gallenne et al. (2021) used SPIPS to show that ℓ Car may possibly exhibit an IR excess, but its low detection level (2.9 ± 3.2% in the K band) is not significant compared to a single-star model (i.e., without an IR excess model). Breitfelder et al. (2016) also provided a SPIPS fitting of ℓ Car that did not present any significant IR excess. However, this work was carried out by considering only the V, J, H, and K photometric bands, that is, without the mid-infrared bands that are critical for detecting such an IR excess. Finally, since a CSE around ℓ Car has been detected by interferometry in the infrared (Kervella et al. 2006, 2009), we enabled the SPIPS to fit an IR excess in the following analysis. The SPIPS provides a fit to the photometry along the pulsation cycle, which is in agreement with the observational data (description in Appendix A.1). Based on Fig. A.1, we emphasize that the SPIPS fitted model is well constrained by the numerous observations and, thus, the physical parameters of the photosphere are accurately derived. However, as noted in the SPIPS original paper, the uncertainties on the derived parameters are purely statistical and do not take systematics into account. While the fitting is satisfactory from the visible domain to the near infrared, the SPIPS fits an IR excess for wavelengths above about 1.2 μm. Indeed the observed brightness (mobs) significantly exceeds the predictions (mkurucz). The averaged IR excess over a pulsation cycle is presented in Fig. 1a. In the next section, we combine this average IR excess with Spitzer mid-infrared spectroscopy.
2.2 Mid-IR excess from Spitzer observations
Space-based observations such as those of obtained with the Spitzer space telescope are useful in avoiding perturbations caused by the Earth’s atmosphere, which is essential for the study of the spectral features of dust. We use spectroscopic observations made with the InfraRed Spectrograph IRS (Houck et al. 2004) on board the Spitzer telescope (Werner et al. 2004). High spectral resolution (R = 600) from 10 μm to 38 μm was used with Short-High (SH) and Long-High (LH) modules placed in the focal plane instrument. The full spectra were retrieved from the CASSIS atlas (Lebouteiller et al. 2011) and the best flux calibrated spectrum was obtained from the optimal extraction using differential method, which eliminates low-level rogue pixels (Lebouteiller et al. 2015). Table 1 provides an overview of the Spitzer observation. The reference epoch at maximum light (MJD0 = 50583.742), pulsation period (P =35.557 days) and the rate of period change (Ṗ = 33.4 ± 1.0 s yr−1) used to derive the pulsation phase corresponding to the Spitzer observation are those computed by the SPIPS fitting of ℓ Car. We note that the most up-to-date and independent value of rate of period change, derived by Neilson et al. (2016) (Ṗ = 20.23 ± 1.38 s yr−1) is slightly lower than the value estimated by SPIPS. However, adopting this latter value would not change the result of this paper.
In order to derive the IR excess of Spitzer observation and correct for the absorption of the interstellar matter (ISM), we followed the method developed in Paper I. First, we derived the IR excess of ℓ Car at the specific Spitzer epoch using:
(1)
where mSpitzer is the magnitude of the Spitzer observation and mkurucz[ϕSpitzer] is the magnitude of the ATLAS9 atmospheric model interpolated at the phase of Spitzer observations (ϕSpitzer), using the parameters derived by SPIPS and given in Table 1. We discarded wavelengths of Spitzer observations longer than 30 μm due to extremely large uncertainties. Secondly, we corrected the spectrum for ISM extinction by subtracting a synthetic ISM composed of silicates from Draine & Lee (1984) (see orange curve in Fig. 1b). This calculation assumes a relation between the extinction E(B − V) derived by SPIPS (i.e., E(B − V) = 0.148 mag) and the silicate absorption from diffuse interstellar medium (see Eq. (4) in Paper I). Finally, we combine this result with SPIPS-averaged IR excess in order to reconstruct the IR excess from the visible to the mid-infrared domain. The determined IR excess is presented in Fig. 1b (red dots and blue line).
Similarly to the five Cepheids presented in Paper I, we observe a continuum IR excess increasing up to − 0.1 mag at 10 μm. In particular, within the specific bands of MATISSE (see vertical strips in Fig. 1b), the IR excess is between − 0.05 and − 0.1 mag at 3.5 μm and 4.5 μm (L and M bands), respectively, which represents an excess of 5–10% above the stellar continuum. In addition, we find no silicate emission feature in the Spitzer spectrum around 10 μm (N band) within the uncertainties, which points toward an absence of significant amount of circumstellar dust with amorphous silicate components.
Spitzer data set of ℓ Car (program number 40968).
![]() |
Fig. 1 (a) Averaged IR excess of ℓ Car as derived with the SPIPS algorithm (red dots) and presented with the fitted ad hoc analytic law |
3 VLTI/MATISSE interferometric observations
MATISSE is the four-telescope beam combiner in the L, M, and N bands of the Very Large Telescope Interferometer (VLTI). The VLTI array consists of four 1.8-m auxiliary telescopes (ATs) and four 8-m unit telescopes (UTs), and provides baseline lengths from 11 m to 150 m. As a spectro-interferometer, MATISSE provides dispersed fringes and photometries. The standard observing mode of MATISSE (so-called hybrid) uses the two following photometric measurement modes: SIPHOT for the L and M bands, in which the photometry is measured simultaneously with the dispersed interference fringes, and HIGH SENS for the N band, in which the photometric flux is measured separately after the interferometric observations. The observations were carried out during the nights of 27 and 28 February 2020 with the so-called large configuration of the ATs quadruplet (A0-G1-J2-J3, with ground baseline lengths from 58 to 132 m) in low spectral resolution (R = λ∕Δλ ≈ 30). The log of the MATISSE observations is given in Table 2. The raw data were processed using the version (1.5.5) of the MATISSE data reduction software2. The steps of the data reduction process are described in Millour et al. (2016).
The MATISSE absolute visibilities are estimated by dividing the measured correlated flux by the photometric flux. Thermal background effects affect the measurement in the M and N bands significantly more than in the L band. Thus, we use only the L band visibilities in our modeling and analysis, while the L and M photometries are used for deriving the observed total flux (see Sect. 3.2). In the N band, the total flux of ℓ Car (≈ 17 Jy) is at the lower limit of the ATs sensitivity with MATISSE for accurate visibility measurements (~ 20 Jy, as stated on the ESO webpage of the MATISSE instrument). Thus, we rather use the correlated flux measurements (not normalized by the photometric flux) as an estimate of the N band photometry of ℓ Car. In the following of the paper, we discarded the 4.1–4.5 μm spectral region, where the atmosphere is not transmissive, and also the noisy edges of the atmospheric spectral bands. As a result, we analyze the following spectral region in L (3.1–3.75 μm) for absolute visibility and flux, M (4.75–4.9 μm) for flux, and in N (8.2–12 μm) for correlated flux.
Log of the observations for nights of 27 February (observations 1 and 2) and 28 February 2020 (observation 3).
3.1 Calibration of the squared visibility in the L band
Calibrators were selected using the SearchCal tool of the JMMC3, with a high confidence level on the derived angular diameter (χ2 ≤ 5, see Appendix A.2 in Chelli et al. 2016). Moreover, calibrator fluxes in the L and N bands have to be of the same order than for the science target (i.e., ≈10 to 100 Jy in L and N) to ensure a reliable calibration (Cruzalèbes et al. 2019). These flux requirements induce the use of partially resolved calibrators at the 132-m longest baseline. The L band uniform-disk (UD) angular diameters for the standard stars as well as the corresponding L and N band fluxes and the spectral type are given in Table 3.
We analyzed these calibration stars with a special care. The first calibrator of this night, q Car, was observed in the frame of a different program and was not used in the interferometric visibility calibration process. Indeed, it is a supergiant star of spectral type KII, which makes its angular diameter of 5.2 mas too large (compared to ℓ Car), and rather uncertain in the L band. However, q Car appears to be a suitable total flux calibrator for ℓ Car in both the L and M bands (see Sect. 3.2). Indeed, although q Car is an irregular long-period variable, its photometric variation in the visible is Δmag ≈ 0.14 (GCVS, Samus’ et al. 2017); thus, its flux variation in the Rayleigh-Jeans domain, that is, in the L and M bands, only represents about 1%. Moreover, q Car was observed at an air mass very close to the one of ℓ Car (see Table 2). In the analysis, we discarded the calibrator B Cen since it has inconsistent visibilities, as explained in Appendix B. Thus, we used ϵ Ant to calibrate observation 1. For observation 2, we also used ϵ Ant. For observation 3, we bracketed the science target with the standard CAL-SCI-CAL strategy calibration using β Vol. The calibrated squared visibilities (V2) in the L band are presented in Fig. 2. The visibility curve associated with the limb-darkened angular diameter of the star (without CSE), which was derived from the SPIPS analysis at the specific phase of VLTI/MATISSE (ϕ = 0.07), is shown for comparison.
MATISSE also provides closure phase measurements, which contain information about the spatial centro-symmetry of the brightness distribution of the source. For all the closure phase measurements, we find an average closure phase of about 0° in the L, M, and the N band from 8.2 to 9.25 μm (see Fig. 3). In the N band, we indicated the typical peak-to-valley dispersion of the MATISSE closure phase measurements with an error bar (≈ 10 degrees) and we discarded the spectral region beyond 10 μm due to extremely large uncertainties (see Fig. 3). Down to a sub-degree level (resp. 10 degree level), our closure phase measurements are consistent with the absence of significant brightness spatial asymmetries in the environment around ℓ Car in the L and M bands (resp. the N band). That justifies our use of centro-symmetric models in the following of the paper.
Calibrator properties.
![]() |
Fig. 2 Calibrated squared visibilities of ℓ Car as a function of the spatial frequencies in the L band for observations 1, 2, and 3. The theoretical visibility curve corresponding to a uniform disk of θUD = 2.887 ± 0.003 mas in the L band, as derived from the SPIPS analysis (i.e., without CSE) is indicated for comparison. |
3.2 Flux calibration
In the L and M bands, we calibrate the total flux Ftot,sci of ℓ Car using
(2)
where Itot,sci and Itot,cal are the observed total raw flux of the science target and the calibrator respectively. Here, Ftot,cal is the known flux of the standard star. We calibrate ℓ Car using standard stars with the closest air mass which are q Car and β Vol for the observations of the first and second night respectively; Ftot,cal is given by their atmospheric templates from Cohen et al. (1999). Since the air mass of both target and calibrator are comparable, we do not correct for the air mass. Moreover, while a chromatic correction exists for the N band (Schütz & Sterzik 2005), it is not calibrated (to our knowledge) for the L and M bands. Then we average the three observations, adding uncertainties and systematics between each measurement in quadrature.
Since ℓ Car is too faint in N band for accurate photometry (and thus absolute visibility) measurements, we use the correlated flux for the N band, which is the flux contribution from the spatially unresolved structures of the source. We calibrated the correlated flux of the science target Fcorr,sci following:
(3)
where Icorr,sci and Icorr,cal are the observed raw correlated flux of the science target and the calibrator respectively, and Vcal the calibrator visibility. For such an absolute calibration of the correlated flux, we need a robust interferometric calibrator with an atmospheric template given by Cohen et al. (1999). Only β Vol, observed during the second night, meets these two requirements. Thus, we only calibrate the correlated flux measurements from observation 3 using the known flux of β Vol. The calibrated correlated fluxes from the different baselines do not present any resolved features, allowing us to average the correlated flux from the different baselines. We note that we have discarded the spectral region between 9.3 and 10 μm due to the presence of a telluric ozone absorption feature around 9.6 μm that strongly impacts the data quality. The calibrated total fluxes in the L and M bands and the averaged calibrated N band correlated flux are presented in Fig. 4.
At this point, we can draw several intermediate conclusions. First, the agreement between the measured flux with MATISSE, Spitzer, and the SPIPS atmosphere model of ℓ Car is qualitatively satisfactory within the uncertainties. However, the MATISSE flux measurements have rather large uncertainties (≥10%), making it difficult to determine the IR excess of ℓ Car with a precision at the few percent level. Secondly, we can see that the N band correlated flux follows a pure Rayleigh-Jeans slope. This indicates the absence of silicate emission, which is consistent with the Spitzer data (Sect. 2.2). This is also in agreement with previous MIDI spectrum given by Kervella et al. (2009). Third, a resolved environment around ℓ Car is seen in the visibility measurements in the L band. Indeed, the calibrated visibilities are significantly lower than the visibilities corresponding to a model of the star without CSE, as derived from the SPIPS algorithm at the specific phase of VLTI/MATISSE observations (see Fig. 2). Since the interferometric closure phase is zero the resolved structure around the star is centro-symmetric. In the next section, we apply a centro-symmetric model of envelope on the observed V2 measurements in the L band in order to investigate the diameter and flux contribution of the envelope.
4 Gaussian envelope model
In this section, we fit a geometrical model on the measured L band visibilities, with ℓ Car and its CSE modeled with a uniform disk (UD) and a surimposed Gaussian distribution, as in previous studies on ℓ Car CSE (Kervella et al. 2006, 2009). The UD diameter corresponds to the one at the specific phase of the observation derived by SPIPS that is θUD= 2.887 ± 0.003 mas in the L band. To model the CSE, we superimposed a Gaussian intensity distribution centered on the stellar disk. This CSE Gaussian model has two parameters which are the CSE diameter θCSE taken as the Full Width at Half Maximum (FWHM), and the CSE flux contribution normalized to the total flux. We then performed a reduced fitting to adjust the total squared visibility of the model
over the MATISSE observation in the L band. We computed the total squared visibility of the star plus CSE model using
(4)
where f = Bp/λ is the spatial frequency (Bp the length of the projected baseline), F⋆ and FCSE are the normalized stellar and CSE flux contribution to the total flux respectively, normalized to unity F⋆ + FCSE = 1 × VUD(f) is the visibility derived from the UD diameter of the star given by
(5)
where J1 is the Bessel function of the first order. VCSE(f) is the visibility of a Gaussian intensity distribution:
(6)
We performed the fit for observations 1, 2, and 3 simultaneously, since ℓ Car is not expected to vary significantly between ϕ = 0.05 and 0.07. The derived visibility model is shown in Fig. 5 and we present the best-fit parameters in Table 4. The visibility model turns out to be consistent, within the error bars, with the three observations. However, the reduced is low (i.e., 0.2), thus, the data are overfitted by the model and the uncertainties on the best-fit parameters are underestimated. In that case, toobtain more reliable uncertainties on the best-fit parameters, we independently fit the three observations and we took the uncertainty as the standard deviation as of the resulted parameters. We resolved a CSE with a radius of 1.76 ± 0.28 R⋆ accounting for about 7.0 ± 1.4% of the flux contribution in the L band, which gives an IR excess of − 0.07 mag. Both the size and the flux contribution are in agreement with previous studies which have found a compact environment of size 1.9 ± 1.4 R⋆ in the K band with VINCI and 3.0 ± 1.1 R⋆ in the N band with MIDI (Kervella et al. 2006), with a few percents in flux contribution in both cases. As noted by Kervella et al. (2006), the large uncertainties in the K band diameter are due to a lack of interferometric data for baselines between 15 to 75 m. We also note that cycle-to-cycle amplitude variations discovered in ℓ Car (Anderson 2014; Anderson et al. 2016) could slightly affect the diameter in the K band. On the other hand, Kervella et al. (2009) found extended emission between 100 and 1000 AU (≈ 100–1000 R⋆), using MIDI and VISIR observations in the N band. Thus, it is possible that the CSE have both a compact and an extended component which are observable at different wavelengths. The CSE flux contribution of 7% is also in agreement with the MATISSE total flux within the uncertainties (see Fig. 4), as well as the determined IR excess of about − 0.05 to − 0.1 mag we derived from Fig. 1b in Sect. 2.2.
![]() |
Fig. 3 Closure phase for observations 1, 2, and 3 for LM and the N bands for each ATs triplet. The double arrows in the N band panels represent the typical peak-to-valley variation. |
![]() |
Fig. 4 Averaged calibrated total flux for observations 1, 2, and 3 in LM bands plus the correlated flux in N band for observation 3. Panels (a) and (b) refer to LM and N bands respectively. SPIPS photosphere model is interpolated at the phase corresponding to the MATISSE the observation with parameters: Teff(ϕobs) = 5014 ± 15 K; log g(ϕobs) = 0.96 ± 0.10; and θUD(ϕobs) = 2.887 ± 0.003 mas. The Spitzer spectrum used in Sect. 2.2 is plotted for comparison. |
![]() |
Fig. 5 Fitted Gaussian CSE around ℓ Car for the combined observations 1, 2, and 3 in the L band. The error on the visibility model is obtained using the covariance matrix of the fitting result. |
Fitted parameters of a Gaussian CSE with flux contribution FCSE(%) and FWHM θCSE(mas) which is also expressed in stellar radius RCSE(R⋆).
5 Physical origin of the IR excess
5.1 Dust envelope
The absence of emission features in the N band in both MATISSE and Spitzer spectra rules out the presence of dust from typical oxygen rich star mineralogy with amorphous silicates or aluminum oxide, which present a characteristic spectral shape in the N band. A pure iron envelope, which presents a continuum emission, is also unlikely to be created. We initially considered and then rejected these possibilities in Paper I (see Figs. 8 and 9 in that paper). We note that the presence of large grains of silicate, with sizes of about 1 μm and more, would lead to a broader emission which could completely disappear (Henning 2010). However, the reason why large grains would preferentially be created in the envelope of Cepheids remains to be explained. In addition, considering the best-fit of the Gaussian CSE model, at a distance of 1.76 R⋆ from the stellar center (0.76 R⋆ from the photosphere), the temperature would exceed 2000 K; therefore, dust cannot survive since it would sublimate (Gail & Sedlmayr 1999). These physical difficulties to reproduce the IR excess with dust envelopes were also recently pointed out by Groenewegen (2020), who fitted the SEDs of 477 Cepheids (including ℓ Car) with a dust radiative transfer code.
5.2 Thin shell of ionized gas
In Paper I, we suggested the presence of a thin shell of ionized gas, with a radius of about 1.15 R⋆, to explain the reconstructed IR excess of five Cepheids. Here, we use the same model for ℓ Car, with a shell of ionized gas, to test its consistency with IR excess and interferometric observations. In this part of the study, we performed a reduced fitting of the shell parameters on the IR excess reconstructed in Sect. 2.2. This simple model has three physical parameters: the shell radius, the temperature, and the mass of ionized gas. An additional parameter is applied as an offset to the whole reconstructed IR excess from Fig. 1b, allowing the model to present either a deficit or an excess in the visible domain if necessary. Indeed, the key point is that when considering a shell of ionized gas, this shell is expected to absorb the light coming from the star in the visible domain, which is currently not considered in the SPIPS algorithm when reconstructing the IR excess (the V magnitude is forced to be zero, see Fig. 1a). However, thanks to the VLTI/MATISSE measurements of flux contribution of the CSE in the L band (7%), we can fix this offset parameter consistently to a 0.03 magnitude in order to force the IR excess in the L band from SPIPS (that is ≈10%) to also reach 7%, which is within the uncertainties. We present the results of the modeling (
) constrained by the spectrum of ℓ Car in Fig. 6. We find a thin shell radius of Rshell = 1.13 ± 0.02 R⋆ which is lower than the Gaussian CSE of 1.76 ± 0.28 R⋆ constrained by MATISSE observations. On the other hand, this model also reproduces the Spitzer spectrum better than the SPIPS predicted photometries. In particular, its value in the L band is about 1% (see Fig. 6) which is not consistent with the value derived using L band visibilities (~ 7%). This model has also derived weak near-infrared excess for several stars in Paper I (see its Fig. 11). This behaviour is physically explained by the free-bound absorption of hydrogen before the free-free emission dominates at longer wavelength. We explored the spatial parameters to test if a larger ionized envelope has the potential to match MATISSE observations with the IR excess. However, our model is not able to accurately reproduce all the observations. Also, even if our model is complex when considering the physics involved (free-free and bound-free opacities description), its geometrical description is rather simple, with a constant density and temperature distribution. Thus, it might not be adapted to the modeling of a larger envelope.
![]() |
Fig. 6 ℓ Car IR excess fitting result of ionized gas shell following the method described in Paper I. Yellow region is the error on the magnitude obtained using the covariance matrix of the fitting result. We find a thin shell of ionized gas with radius Rshell = 1.13 ± 0.02 R⋆; temperature Tshell = 3791 ± 85 K; mass of ionized gas M = 9.10 × 10−8 ± 8.5 × 10−9 M⊙, using an ad hoc offset of +0.03 mag on the data in order to match the IR excess found by VLTI/MATISSE and the one of SPIPS (see the text); and |
5.3 Limit of the model
The derived parameters of the shell of ionized gas depend on the initial IR excess derived by SPIPS. Any systematics in this calculation could in turn affect the IR excess and the parameters of the shell of ionized gas. As noted by Hocdé et al. (2020b), the derived distance of SPIPS would be systematically affected if the shell of ionized gas is not taken into account because of its absorption in the visible. This systematic should be on the level of only a few percents (10Δm∕2.5), with all parameters remaining unchanged. The direct implementation of the gas model into the SPIPS fitting has to be tested in forthcoming studies to correct this uncertainty. Moreover, Gallenne et al. (2021) studied the ad hoc IR excess derived by SPIPS for a larger sample of stars and found that angular diameter and the colour excess are the most impacted parameters when no ad hoc IR excess model is included. However, in the case of ℓ Car, these parameters appear to be well constrained by angular diameter measurements, the reddening and the distance are also in agreement with values found in the literature (see Appendix A).
5.4 Perspectives
An interesting physical alternative is to consider free-free emission produced by negative hydrogen ion H−, with free electrons provided by metals. Indeed for an envelope of 2 R⋆ the average temperature is about 3000 K and the hydrogen is neutral. Thus, in that case, most of the free electrons would be provided by silicon, iron,and magnesium which have a mean first ionization potential of 7.89 eV. As neutral hydrogen is able to form negative hydrogen ion, H− could generate an infrared excess with free-free emission. This phenomenon is known to produce a significant IR excess compared to dust emission in the extended chromosphere of cool supergiant stars such as Betelgeuse (Gilman 1974; Humphreys 1974; Lambert & Snell 1975; Altenhoff et al. 1979; Skinner & Whitmore 1987). Moreover, the photo-detachment potential occurs at 1.6 μm for the H− mechanism, thus, free-free emission should dominate the IR excess above this wavelength, as it is suggested by the SPIPS cycle-averaged IR excess. This single photo-detachment has also the advantage to avoid the important bound-free absorption in the near-infrared produced in the preceding model. We thus suggest that the compact structure we resolved around ℓ Car has the potential to produce IR excess through H− free-free mechanism via an analogy with extended chromosphere of cool supergiant stars. Further investigations are necessary to confirm this hypothesis.
6 Conclusions
Determining the nature and the occurrence of the CSEs of Cepheids is of considerable interest in quantifying their impact on the period-luminosity relation, and also for understanding the mass-loss mechanisms at play. In this paper, we constrain both the IR excess and the geometry of the CSE of ℓ Car using both Spitzer low-resolution spectroscopy and MATISSE interferometric observations in the mid-infrared. Assuming an IR excess ad hoc model, we used SPIPS to derive the photospheric parameters of the star in order to compare with Spitzer and MATISSE at their respective phase of observation. This analysis allows to derive the IR excess from Spitzer and also to derive the CSE properties in the L band spatially resolved with MATISSE. These observations lead to the following conclusions with regard to the physical nature of ℓ Car’s CSE:
- 1.
We resolved a centro-symmetric and compact structure in the L band with VLTI/MATISSE that has a radius of about 1.76 R⋆. The flux contribution is about 7%.
- 2.
We find no clear evidence for dust emission features around 10 μm (N band) in MATISSE and Spitzer spectra, which suggests an absence of circumstellar dust.
- 3.
Our dedicated model of a shell of ionized gas better reproduces the mid-infrared Spitzer spectrum and implies a size for the CSE that is lower than the one derived from VLTI/MATISSE observations (1.13 ± 0.02 R⋆ versus 1.76 ± 0.28 R⋆), as well as a lower flux in the L band (1% versus 7%).
- 4.
We suggest that improving our model of shell of ionized gas by including the free-free emission from negative hydrogen ion is likely to help in reproducing the observations, that is, with regard to the size of the CSE and the IR excess, particularly in the L band.
While the compact CSE of ℓ Car is likely gaseous, the exact physical origin of the IR excess remains uncertain. Further observations of Cepheids that depend on both their pulsation period and their location in the HR diagram are necessary for understanding the CSE’s IR excess. This could be a key to unbias period-luminosity relation from Cepheids IR excess.
Acknowledgements
The authors acknowledge the support of the French Agence Nationale de la Recherche (ANR), under grant ANR-15-CE31-0012-01 (project UnlockCepheids). The research leading to these results has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program under grant agreement no. 695099 (project CepBin). This research was supported by the LP2018-7/2019 grant of the Hungarian Academy of Sciences. This research has been supported by the Hungarian NKFIH grant K132406. This research made use of the SIMBAD and VIZIER (http://cdsweb.u-strasbg.fr/) databases at CDS, Strasbourg (France) and the electronic bibliography maintained by the NASA/ADS system. This research also made use of Astropy, a community-developed corePython package for Astronomy (Astropy Collaboration 2018). Based on observations made with ESO telescopes at Paranal observatory under program IDs: 0104.D-0554(A). This research has benefited from the help of SUV, the VLTIuser support service of the Jean-Marie Mariotti Center (http://www.jmmc.fr/suv.htm). This research has also made use of the Jean-Marie Mariotti Center Aspro service (http://www.jmmc.fr/aspro).
Appendix A SPIPS data set and fitted pulsational model of the star sample
Figure A.1 is organized as follows: pulsational velocity, effective temperature, and angular diameter curves according to the pulsation phase are shown on the left panels, while the right panels display photometric data in various bands. Above the figure, the projection factor set to p = 1.270 is indicated, along with the fitted distance d using parallax-of-pulsation method, the fitted color excess E(B − V), and the ad hoc IR excess law.
In this work, we arbitrarily set the p-factor to 1.27 following the general agreement around this value for several stars in the literature (Nardetto et al. 2004, 2009; Mérand et al. 2005). However, there is no agreement for the optimum value to use and this is also the case for ℓ Car. For this star, Nardetto et al. (2007, 2009) found a p-factor of 1.28 ± 0.02 and 1.22 ± 0.04 using different observational techniques that are both in agreement with 1.27 within about 1 sigma. Thus, a possible systematic cannot be excluded although the distance derived by SPIPS (d = 519.2 ± 4.5 pc) is in agreement, for example, with the distance obtained using Gaia parallaxes (505 ± 28 pc; Gaia Collaboration 2018) or the one derived using a recent period-luminosity calibrated for Cepheids in the Milky Way (544 ± 32 pc) (Breuval et al. 2020).
Gallenne et al. (2021) have shown that the absence of IR excess model in the SPIPS fit can affect the derivation of the angular diameter and the colour excess. In our fit, based on an IR excess model, we emphasize the good agreement of the angular diameter model with the data (panel b) and also the consistency of the derived colour excess with the one given by Stilism 3D extinction map (Capitanio et al. 2017), that is, 0.148 ± 0.006 versus 0.116 ± 0.26. In the photometric panels, the gray dashed line corresponds to the magnitude of the SPIPS model without CSE. It actually corresponds to the magnitude of a Kurucz atmosphere model, mkurucz, obtained with the ATLAS9 simulation code from Castelli & Kurucz (2003), with solar metallicity and a standard turbulent velocity of 2 km s−1. These models are indeed suitable for deriving precise synthetic photometry (Castelli 1999), in comparison with MARCS models giving rough flux estimates5. The gray line corresponds to the best SPIPS model, which is composed of the latter model without CSE plus an IR excess model. In the angular diameter panels, the gray curve corresponds to limb-darkened (LD) angular diameters. For stars with solar metallicity, when effective temperature is low enough, CO molecules can form in the photosphere and absorb light in the CO band-head at 4.6 μm (Scowcroft et al. 2016). This effect is observed in the Spitzer I2 IRAC data set. In this case, these data were ignored during the fitting of SPIPS. Panels presenting a horizontal blue bar contain only one data with undetermined phase of observation, thus these photometric mid-infrared bands were not used in the IR excess reconstruction in Sect. 2.1 and 2.2
![]() |
Fig. A.1 SPIPS results of ℓ Car. Velocity: Anderson (2014). Effective temperature: no data. Angular diameter: Kervella et al. (2006). Photometry: Madore (1975), Pel (1976), Neugebauer et al. (1984), Laney & Stobie (1992), ESA (1997), Berdnikov & Turner (2002), Smith et al. (2004), Ammons et al. (2006), Wright et al. (2010), Ishihara et al. (2010), Monson et al. (2012), Gaia Collaboration (2018). |
Appendix B Verification of the diameters of the calibrators
In order to verify the consistency of the standard stars of the first night (27 February), we consider all calibrators available during the night (B Cen, e Cen, ϵ Ant) and weinter-calibrate them. Thus, each calibrator is calibrated using the two others. We performed six different calibrations in total, following Table B.1 and Fig. B.1. The results are also compared with the known values from the literature.
From these results, the diameters of ϵ Ant and e Cen are in excellent agreement with those given by the JSDC when they are calibrated from one to another, and the visibilities are well fitted by an uniform disk (UD) model (see Figs. B.1a and b). On the contrary, we find important inconsistencies when deriving the diameter of B Cen or using it as a calibrator (see Figs. B.1c to f). In particular this standard star is suspected to be 20% larger than the diameter given by the JSDC (≈ 3 mas instead of 2.5 mas). Moreover, a simple UD model seems to be unsuitable for fitting the observed visibilities (see black curve in Figs. B.1e and f). On the other hand, the diameter of B Cen is also well established by interferometry in the K band (Richichi et al. 2009) and was used many times as an interferometric calibrator (Kervella et al. 2006; Kervella 2007). We emphasize that the transfer function was very stable during this night, thus this result cannot be attributed to the atmospheric conditions (see coherence time in Table 2). Since the origin of this discrepancy remains unknown, this analysis has prevented us from using B Cen to calibrate the Cepheid ℓ Car.
Diameters in milliarcseconds obtained by the inter-calibration between B Cen, ϵ Ant, and e Cen.
![]() |
Fig. B.1 Intercalibration of the different calibrators and comparison with JSDC diameters. (a) and (b) Diameters of ϵ Ant and e Cen are in agreement with the JSDC when they are calibrated from one another. (c) and (d) On the contrary, the diameters of ϵ Ant and e Cen are inconsistent with the JSDC catalogue when they are calibrated with B Cen. (e) and (f) B Cen’s diameter is inconsistent with the JSDC catalogue. |
References
- Allouche, F., Robbe-Dubois, S., Lagarde, S., et al. 2016, Proc. SPIE, 9907, 99070C [Google Scholar]
- Altenhoff, W. J., Oster, L., & Wendker, H. J. 1979, A&A, 73, L21 [NASA ADS] [Google Scholar]
- Ammons, S. M., Robinson, S. E., Strader, J., et al. 2006, ApJ, 638, 1004 [NASA ADS] [CrossRef] [Google Scholar]
- Anderson, R. I. 2014, A&A, 566, L10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Anderson, R. I., Mérand, A., Kervella, P., et al. 2016, MNRAS, 455, 4231 [NASA ADS] [CrossRef] [Google Scholar]
- Astropy Collaboration (Price-Whelan, A. M., et al.) 2018, AJ, 156, 123 [Google Scholar]
- Berdnikov, L. N., & Turner, D. G. 2002, VizieR On-line Data Catalog: J/ApJS/137/209 [Google Scholar]
- Bohm-Vitense, E., & Love, S. G. 1994, ApJ, 420, 401 [NASA ADS] [CrossRef] [Google Scholar]
- Bourgés, L., Lafrasse, S., Mella, G., et al. 2014, in ASP Conf. Ser., 485, Astronomical Data Anaylsis Softward and Systems XXIII, eds. N. Manset, & P. Forshay, 223 [Google Scholar]
- Breitfelder, J., Mérand, A., Kervella, P., et al. 2016, A&A, 587, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Breuval, L., Kervella, P., Anderson, R. I., et al. 2020, A&A, 643, A115 [EDP Sciences] [Google Scholar]
- Capitanio, L., Lallement, R., Vergely, J. L., Elyajouri, M., & Monreal-Ibero, A. 2017, A&A, 606, A65 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Castelli, F. 1999, A&A, 346, 564 [NASA ADS] [Google Scholar]
- Castelli, F., & Kurucz, R. L. 2003, in IAU Symp., 210, Modelling of Stellar Atmospheres, eds. N. Piskunov, W. W. Weiss, & D. F. Gray, A20 [Google Scholar]
- Chelli, A., Duvert, G., Bourgès, L., et al. 2016, A&A, 589, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cohen, M., Walker, R. G., Carter, B., et al. 1999, AJ, 117, 1864 [Google Scholar]
- Cruzalèbes, P., Petrov, R. G., Robbe-Dubois, S., et al. 2019, MNRAS, 490, 3158 [Google Scholar]
- Draine, B. T.,& Lee, H. M. 1984, ApJ, 285, 89 [Google Scholar]
- ESA, 1997, ESA SP, 1200, The HIPPARCOS and TYCHO Catalogues. Astrometric and Photometric Star Catalogues Derived from the ESA HIPPARCOS Space Astrometry Mission [Google Scholar]
- Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gail, H.-P., & Sedlmayr, E. 1999, A&A, 347, 594 [Google Scholar]
- Gallenne, A., Mérand, A., Kervella, P., et al. 2013, A&A, 558, A140 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gallenne, A., Kervella, P., Mérand, A., et al. 2017, A&A, 608, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gallenne, A., Mérand, A., Kervella, P., et al. 2021, A&A, in press https://doi.org/10.1051/0004-6361/202140350 [Google Scholar]
- Gilman, R. C. 1974, ApJ, 188, 87 [CrossRef] [Google Scholar]
- Groenewegen, M. A. T. 2020, A&A, 635, A33 [EDP Sciences] [Google Scholar]
- Henning, T. 2010, ARA&A, 48, 21 [NASA ADS] [CrossRef] [Google Scholar]
- Hocdé, V., Nardetto, N., Borgniet, S., et al. 2020a, A&A, 641, A74 [EDP Sciences] [Google Scholar]
- Hocdé, V., Nardetto, N., Lagadec, E., et al. 2020b, A&A, 633, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Houck, J. R., Roellig, T. L., van Cleve, J., et al. 2004, ApJS, 154, 18 [NASA ADS] [CrossRef] [Google Scholar]
- Humphreys, R. M. 1974, ApJ, 188, 75 [Google Scholar]
- Ishihara, D., Onaka, T., Kataza, H., et al. 2010, A&A, 514, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Javanmardi, B., Mérand, A., Kervella, P., et al. 2021, ApJ, 911, 12 [Google Scholar]
- Kervella, P. 2007, A&A, 464, 1045 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kervella, P., Mérand, A., Perrin, G., & Coudé du Foresto, V. 2006, A&A, 448, 623 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kervella, P., Mérand, A., & Gallenne, A. 2009, A&A, 498, 425 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kervella, P., Trahin, B., Bond, H. E., et al. 2017, A&A, 600, A127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Lambert, D. L., & Snell, R. L. 1975, MNRAS, 172, 277 [NASA ADS] [CrossRef] [Google Scholar]
- Laney, C. D., & Stobie, R. S. 1992, A&AS, 93, 93 [NASA ADS] [Google Scholar]
- Lebouteiller, V., Barry, D. J., Spoon, H. W. W., et al. 2011, ApJS, 196, 8 [NASA ADS] [CrossRef] [Google Scholar]
- Lebouteiller, V., Barry, D. J., Goes, C., et al. 2015, ApJS, 218, 21 [NASA ADS] [CrossRef] [Google Scholar]
- Lopez, B., Lagarde, S., Jaffe, W., et al. 2014, The Messenger, 157, 5 [NASA ADS] [Google Scholar]
- Madore, B. F. 1975, ApJS, 29, 219 [NASA ADS] [CrossRef] [Google Scholar]
- Mérand, A., Kervella, P., Coudé du Foresto, V., et al. 2005, A&A, 438, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Mérand, A., Kervella, P., Coudé du Foresto, V., et al. 2006, A&A, 453, 155 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Mérand, A., Aufdenberg, J. P., Kervella, P., et al. 2007, ApJ, 664, 1093 [NASA ADS] [CrossRef] [Google Scholar]
- Mérand, A., Kervella, P., Breitfelder, J., et al. 2015, A&A, 584, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Millour, F., Berio, P., Heininger, M., et al. 2016, Proc. SPIE, 9907, 990723 [Google Scholar]
- Monson, A. J., Freedman, W. L., Madore, B. F., et al. 2012, ApJ, 759, 146 [NASA ADS] [CrossRef] [Google Scholar]
- Nardetto, N., Fokin, A., Mourard, D., et al. 2004, A&A, 428, 131 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Nardetto, N., Mourard, D., Mathias, P., Fokin, A., & Gillet, D. 2007, A&A, 471, 661 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Nardetto, N., Groh, J. H., Kraus, S., Millour, F., & Gillet, D. 2008, A&A, 489, 1263 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Nardetto, N., Gieren, W., Kervella, P., et al. 2009, A&A, 502, 951 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Nardetto, N., Mérand, A., Mourard, D., et al. 2016, A&A, 593, A45 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Neilson, H. R., Engle, S. G., Guinan, E. F., Bisol, A. C., & Butterworth, N. 2016, ApJ, 824, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Neugebauer, G., Habing, H. J., van Duinen, R., et al. 1984, ApJ, 278, L1 [Google Scholar]
- Pel, J. W. 1976, A&AS, 24, 413 [NASA ADS] [Google Scholar]
- Richichi, A., Percheron, I., & Davis, J. 2009, MNRAS, 399, 399 [NASA ADS] [CrossRef] [Google Scholar]
- Robbe-Dubois, S., Lagarde, S., Antonelli, P., et al. 2018, Proc. SPIE, 10701, 107010H [Google Scholar]
- Rodgers, A. W., & Bell, R. A. 1968, MNRAS, 138, 23 [NASA ADS] [CrossRef] [Google Scholar]
- Samus’, N. N., Kazarovets, E. V., Durlevich, O. V., Kireeva, N. N., & Pastukhova, E. N. 2017, Astron. Rep., 61, 80 [NASA ADS] [CrossRef] [Google Scholar]
- Schmidt, E. G. 2015, ApJ, 813, 29 [NASA ADS] [CrossRef] [Google Scholar]
- Schütz, O., & Sterzik, M. 2005, in High Resolution Infrared Spectroscopy in Astronomy, 104 [CrossRef] [Google Scholar]
- Scowcroft, V., Seibert, M., Freedman, W. L., et al. 2016, MNRAS, 459, 1170 [Google Scholar]
- Skinner, C. J., & Whitmore, B. 1987, MNRAS, 224, 335 [NASA ADS] [Google Scholar]
- Smith, B. J., Price, S. D., & Baker, R. I. 2004, ApJS, 154, 673 [NASA ADS] [CrossRef] [Google Scholar]
- Werner, M. W., Roellig, T. L., Low, F. J., et al. 2004, ApJS, 154, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868 [Google Scholar]
The MATISSE reduction pipeline is publicly available at http://www.eso.org/sci/software/pipelines/matisse/
SearchCal is publicly available at https://www.jmmc.fr/english/tools/proposal-preparation/search-cal/
All Tables
Log of the observations for nights of 27 February (observations 1 and 2) and 28 February 2020 (observation 3).
Fitted parameters of a Gaussian CSE with flux contribution FCSE(%) and FWHM θCSE(mas) which is also expressed in stellar radius RCSE(R⋆).
Diameters in milliarcseconds obtained by the inter-calibration between B Cen, ϵ Ant, and e Cen.
All Figures
![]() |
Fig. 1 (a) Averaged IR excess of ℓ Car as derived with the SPIPS algorithm (red dots) and presented with the fitted ad hoc analytic law |
In the text |
![]() |
Fig. 2 Calibrated squared visibilities of ℓ Car as a function of the spatial frequencies in the L band for observations 1, 2, and 3. The theoretical visibility curve corresponding to a uniform disk of θUD = 2.887 ± 0.003 mas in the L band, as derived from the SPIPS analysis (i.e., without CSE) is indicated for comparison. |
In the text |
![]() |
Fig. 3 Closure phase for observations 1, 2, and 3 for LM and the N bands for each ATs triplet. The double arrows in the N band panels represent the typical peak-to-valley variation. |
In the text |
![]() |
Fig. 4 Averaged calibrated total flux for observations 1, 2, and 3 in LM bands plus the correlated flux in N band for observation 3. Panels (a) and (b) refer to LM and N bands respectively. SPIPS photosphere model is interpolated at the phase corresponding to the MATISSE the observation with parameters: Teff(ϕobs) = 5014 ± 15 K; log g(ϕobs) = 0.96 ± 0.10; and θUD(ϕobs) = 2.887 ± 0.003 mas. The Spitzer spectrum used in Sect. 2.2 is plotted for comparison. |
In the text |
![]() |
Fig. 5 Fitted Gaussian CSE around ℓ Car for the combined observations 1, 2, and 3 in the L band. The error on the visibility model is obtained using the covariance matrix of the fitting result. |
In the text |
![]() |
Fig. 6 ℓ Car IR excess fitting result of ionized gas shell following the method described in Paper I. Yellow region is the error on the magnitude obtained using the covariance matrix of the fitting result. We find a thin shell of ionized gas with radius Rshell = 1.13 ± 0.02 R⋆; temperature Tshell = 3791 ± 85 K; mass of ionized gas M = 9.10 × 10−8 ± 8.5 × 10−9 M⊙, using an ad hoc offset of +0.03 mag on the data in order to match the IR excess found by VLTI/MATISSE and the one of SPIPS (see the text); and |
In the text |
![]() |
Fig. A.1 SPIPS results of ℓ Car. Velocity: Anderson (2014). Effective temperature: no data. Angular diameter: Kervella et al. (2006). Photometry: Madore (1975), Pel (1976), Neugebauer et al. (1984), Laney & Stobie (1992), ESA (1997), Berdnikov & Turner (2002), Smith et al. (2004), Ammons et al. (2006), Wright et al. (2010), Ishihara et al. (2010), Monson et al. (2012), Gaia Collaboration (2018). |
In the text |
![]() |
Fig. B.1 Intercalibration of the different calibrators and comparison with JSDC diameters. (a) and (b) Diameters of ϵ Ant and e Cen are in agreement with the JSDC when they are calibrated from one another. (c) and (d) On the contrary, the diameters of ϵ Ant and e Cen are inconsistent with the JSDC catalogue when they are calibrated with B Cen. (e) and (f) B Cen’s diameter is inconsistent with the JSDC catalogue. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.