Issue |
A&A
Volume 525, January 2011
|
|
---|---|---|
Article Number | A22 | |
Number of page(s) | 11 | |
Section | Interstellar and circumstellar matter | |
DOI | https://doi.org/10.1051/0004-6361/201015194 | |
Published online | 29 November 2010 |
Fast ray-tracing algorithm for circumstellar structures (FRACS)
II. Disc parameters of the B[e] supergiant CPD-57° 2874 from VLTI/MIDI data⋆,⋆⋆
1
Laboratoire H. Fizeau, UMR CNRS 6525, Université de Nice-Sophia Antipolis
(UNS), Observatoire de la Côte d’Azur (OCA),
Campus Valrose,
06108
Nice Cedex 2,
France
e-mail: Armando.Domiciano@obs-azur.fr
2
Observatório Nacional, Rua General José Cristino, 77, 20921-400
São Cristovão, Rio de Janeiro,
Brazil
3
Instituto de Astronomia, Geofísica e Ciências Atmosféricas,
Universidade de São Paulo (USP), Rua do Matão 1226, Cidade
Universitária, São
Paulo, SP - 05508-900,
Brazil
4
Max-Planck-Institut für Radioastronomie,
Auf dem Hügel 69,
53121
Bonn,
Germany
5
Deutsches Zentrum für Luft- und Raumfahrt, Raumfahrt-Agentur,
Königswinterer Strasse 522-524,
53227
Bonn,
Germany
Received: 10 June 2010
Accepted: 15 September 2010
Context. B[e] supergiants are luminous, massive post-main sequence stars exhibiting non-spherical winds, forbidden lines, and hot dust in a disc-like structure. The physical properties of their rich and complex circumstellar environment (CSE) are not well understood, partly because these CSE cannot be easily resolved at the large distances found for B[e] supergiants (typically ≳1 kpc).
Aims. From mid-IR spectro-interferometric observations obtained with VLTI/MIDI we seek to resolve and study the CSE of the Galactic B[e] supergiant CPD-57° 2874.
Methods. For a physical interpretation of the observables (visibilities and spectrum) we use our ray-tracing radiative transfer code (FRACS), which is optimised for thermal spectro-interferometric observations.
Results. Thanks to the short computing time required by FRACS (<10 s per monochromatic model), best-fit parameters and uncertainties for several physical quantities of CPD-57° 2874 were obtained, such as inner dust radius, relative flux contribution of the central source and of the dusty CSE, dust temperature profile, and disc inclination.
Conclusions. The analysis of VLTI/MIDI data with FRACS allowed one of the first direct determinations of physical parameters of the dusty CSE of a B[e] supergiant based on interferometric data and using a full model-fitting approach. In a larger context, the study of B[e] supergiants is important for a deeper understanding of the complex structure and evolution of hot, massive stars.
Key words: stars: individual: CPD-57° 2874 / methods: observational / methods: numerical / techniques: high angular resolution / techniques: interferometric / stars: emission-line, Be
Based on VLTI/MIDI observations collected at the European Southern Observatory (ESO), Paranal, Chile under ESO Programmes 074.D-0101 and 078.D-0213. Also based on observations at the ESO 2.2-m telescope, La Silla, Chile, under agreement with the Observatório Nacional-MCT (Brazil).
Figure 5 is only available in electronic form at http://www.aanda.org
© ESO, 2010
1. Introduction
The B[e] stars form a heterogeneous group composed of objects at different evolutionary stages, but with many similar observational signatures characterising the so-called “B[e] phenomenon”: simultaneous presence of hydrogen emission lines, low-excitation forbidden and permitted metallic lines in emission, and a significant infrared (IR) excess mainly caused by hot circumstellar dust. Another common property of B[e] stars is the presence of a non-spherical circumstellar environment (hereafter CSE, e.g. Zickgraf 2003).
Lamers et al. (1998) defined five B[e] sub-classes, one of which contains unclassified stars. Miroshnichenko (2007) proposed an additional B[e] sub-class (the FS CMa stars) to explain at least part of the unclassified B[e]-type stars as binaries at a phase of ongoing or recently ended rapid mass transfer and dust formation.
One of the B[e] sub-classes is composed of luminous (log (L ⋆ /L⊙) ≳ 4) post-main sequence objects: the B[e] supergiant stars (hereafter sgB[e]). Previous spectroscopic and polarimetric observations of sgB[e] (e.g. Zickgraf et al. 1985; Magalhães 1992) show that the wind of these massive and luminous stars is composed of two distinct components: (1) a wind of low density and high velocity and (2) a wind of high density and low velocity. Zickgraf et al. (1985) proposed a picture where the sgB[e] winds consist of a hot and fast radiation-driven polar wind and a slow, much cooler and denser (by a factor of 102 or 103) equatorial wind. This disc-like structure provides a natural explanation for the existence of dust around those objects, because the presence of dust requires regions of sufficiently high density and low kinetic temperatures. One possible explanation for this two-component CSE is that rapid rotation of the central star leads to the formation of an equatorial disc because of the combination of rotation-induced bi-stability and rotation-induced wind compression (Lamers & Pauldrach 1991; Bjorkman 1998; Pelupessy et al. 2000). Other mechanisms, such as binarity, are evoked to explain the disc-like CSE, but in any case, rapid rotation seems to play a key role in the origin of these discs (e.g. Meynet & Maeder 2006). Owing to their physical characteristics (fast rotation, disc-like CSE, high luminosity, evolutionary status), it has also been suggested (Vink et al. 2009) that sgB[e] might share evolutionary links with rapidly rotating O-stars and long-duration gamma-ray bursts (GRBs).
Because of the large distances of sgB[e] (≳1 kpc for the closest ones), the geometry and physical structure (e.g. density and temperature distribution) of their CSE could be only quite recently directly probed, thanks to modern high angular resolution (HAR) techniques. For example, Domiciano de Souza et al. (2008) used ESO’s VLT/VISIR instrument to directly measure the typical size of the dusty CSE of the sgB[e] MWC 300 from diffraction-limited mid-IR images. Among HAR techniques, optical/IR long baseline interferometry (OLBI) provides the highest resolving power, allowing the study of sgB[e] CSEs at angular resolutions ~1−10 milliarcsec (mas). For example, 3 combined adaptive optics (VLT/NACO) and OLBI (VLTI/AMBER and VLTI/MIDI) to detect a companion around the sgB[e] HD87643 and to perform an extensive study of the dusty CSE of the binary system.
In the examples above as well as in most works based on HAR data of sgB[e], two different strategies are commonly adopted to interpret the observations: (1) geometrical analytical modelling and (2) radiative transfer modelling (e.g. using the Monte Carlo approach). The geometrical analytical models have the advantage to be very fast to calculate, allowing a full model-fitting (for example a χ2 minimisation) and error estimate of the model parameters (mostly geometrical parameters of the CSE). However, these simple models do not give access to physical parameters of the targets such as temperature and density distributions, optical depths, etc. On the other hand, most radiative transfer models present a consistent description of the physical conditions of the CSE. However, because these models are quite complex, they demand a lot of computing time, which prevents one from exploring a large domain of the parameter space and also from obtaining a good estimate of the uncertainties on the fitted parameters. In this work we adopt a third approach for the data modelling, which tries to keep the advantages of the other approaches, without the drawbacks. To this aim we use our fast ray-tracing algorithm for circumstellar structures (FRACS), which is based on a parametrised CSE combined to a simplified radiative transfer (no scattering). A complete description of this algorithm is given in Niccolini et al. (2010, hereafter Paper I).
In the present paper we apply FRACS to study the CSE of the Galactic sgB[e] CPD-57° 2874 (also named Hen 3-394, WRAY 15-535) based on mid-IR spectro-interferometric observations performed with ESO’s VLTI/MIDI beam-combiner instrument. Previous near- and mid-IR interferometric observations of CPD-57° 2874 directly revealed an elongated CSE that is compatible with a disc-like structure formed by gas and dust (Domiciano de Souza et al. 2007, hereafter DS07). However, because only a limited number (four) of baselines was available and since the authors adopted simple analytical models, only geometrical parameters could be derived from this first analysis of CPD-57° 2874. As shown below, the use of FRACS allowed us to confirm the previous results and, most importantly, to derive physical parameters for this Galactic sgB[e].
In Sect. 2 we give the log of the VLTI/MIDI observations and describe the data reduction procedure. In Sect. 3 we provide a new distance estimate of CPD-57° 2874 obtained from spectroscopic observations with FEROS. A short reminder of the ray-tracing code FRACS is presented in Sect. 4, followed by the results obtained from a model-fitting analysis of the VLTI/MIDI observations (Sect. 5). A discussion of the results and the conclusions of this work are presented in Sects. 6 and 7, respectively.
2. VLTI/MIDI observations
The interferometric observations of CPD-57° 2874 were performed with MIDI, the mid-infrared 2-telescope beam-combiner instrument of ESO’s VLTI (Leinert et al. 2004). All four 8.2 m unit telescopes (UTs) were used. The N-band spectrum as well as spectrally dispersed fringes have been recorded between λ ≃ 7.5 μm and λ ≃ 13.5 μm with a spectral resolution of R ≃ 30 using a prism. In total, nB = 10 data sets have been obtained with projected baselines (Bproj) ranging from ≃40 m to ≃130 m, and baseline position angles (PA) between ≃8° and ≃105° (from North to East). A summary of the VLTI/MIDI observations of CPD-57° 2874 is given in Table 1, and the corresponding uv-plane coverage is shown in Fig. 1.
Summary of VLTI/MIDI observations of CPD-57° 2874: data set index, date, Coordinated Universal Time (UTC) of observation, baseline configuration, projected baseline length and position angle.
Fig. 1 uv-Plane coverage: projected baselines (length and position angle) for the VLTI/MIDI observations of CPD-57° 2874 (further details are given in Table 1). |
The MIDI data were reduced with the MIA+EWS data reduction package, which includes two different sub-packages: the MIA package developed at the Max-Planck-Institut für Astronomie, and the EWS package developed at the Leiden Observatory1.
While MIA is based on the power spectrum analysis, which measures the total power of observed fringes (Leinert et al. 2004), EWS coherently adds the fringes after correction for optical path differences (instrumental as well as atmospheric delays) in each scan (Jaffe et al. 2004). The data reduction results obtained with the MIA and EWS packages agree well within de uncertainties.
Observation log and angular diameters of calibrators (values from DS07) used to derive the calibrated N-band visibilities of CPD-57° 2874.
The instrumental transfer function at each spectral channel was obtained from the observations of calibrator stars with known uniform-disc diameters (⌀UD). The calibrators used in the data reduction and the adopted angular diameters and uncertainties are listed in Table 2. The calibrated visibilities were calculated from the ratio of the targets’ raw visibilities and the average transfer function derived from the calibrator measurements of the corresponding night (Fig. 3). The error of the calibrated MIDI visibilities is of the order of ≃5−10% and includes the raw visibility error as well as the error of the transfer function. The uncertainties on the calibrator angular diameter are negligible compared to the standard deviation of the transfer function. Usually, ≃3−5 calibrator measurements per night were available (Table 2). In the few cases where only one suitable calibrator observation was available we assumed a typical transfer function error of 5% to estimate the errors on the final calibrated visibilities.
VLTI/MIDI also provides spectral fluxes of CPD-57° 2874 in the N-band. On average, all fluxes are compatible within ≃10% with those previously presented by DS07 (MIDI and ISO-SWS spectra). We here considered the VLTI/MIDI spectrum used by DS07 with an uncertainty of 20% (Fig. 4). This larger uncertainty ensures a complete agreement between the observed MIDI and ISO fluxes at all wavelengths. We note that the mid-IR spectrum of CPD-57° 2874 does not show any clear evidence of an important silicate feature around 10 μm.
3. FEROS observations and distance estimate
In addition to our mid-IR interferometric observations, we obtained high-resolution optical spectra of CPD-57° 2874. The spectra were recorded with the high-resolution Fiber-fed Extended Range Optical Spectrograph (FEROS), attached to the 2.2-m telescope at ESO in La Silla (Chile). FEROS is a bench-mounted Echelle spectrograph with fibers, covering a sky area of ≃2″ × 2″ and a wavelength range from 3600 Å to 9200 Å. Its spectral resolution is R ≃ 55000 (around 6000 Å). We have adopted its complete automatic online reduction, which includes the heliocentric correction. The FEROS spectra were obtained on 2008 December 21. We recorded two exposures of 1000 s with S/N of ~60 in the 5500 Å region.
These observations were used to estimate the distance of CPD-57° 2874. A previous distance estimate of d = 2.5 kpc has been proposed by 2, assuming that this star belongs to the Carina OB association.
From our FEROS high-resolution spectra it is possible to estimate the distance of CPD-57° 2874, based on the statistical relation cited by Allen (1973). This relation uses the equivalent widths of interstellar Na i lines. In our data, each Na i line is composed of two absorption components. However, owing to the lack of data from different epochs, it is impossible to see any temporal changes, which would allow us to derive a possible circumstellar contamination. We have therefore assumed that both Na i components are of interstellar origin. The measured equivalent widths are given in Table 3. Our estimated distance for CPD-57° 2874 is d = 1.7 kpc with an uncertainty of 0.7 kpc. This large error is firstly due to a possible contamination from the circumstellar emission component and the saturation of the absorption one, and secondly to a systematic error caused by the statistical relation used. Within the error bars, our distance estimate is roughly compatible with the result of 2. We have considered both distances in our analysis: 1.7 and 2.5 kpc.
Equivalent widths (EW) of the Na i absorption lines (λλ 5890 Å, 5896 Å), obtained from our FEROS data.
4. Description of the model FRACS
Here we present a short description of our numerical model FRACS and the parametrisation adopted to describe a sgB[e]. A full description of FRACS is given in Paper I.
FRACS is based on the ray-tracing technique using quadtree meshes for the grid and the full symmetries of the problem to hand to significantly decrease the computing time necessary to obtain monochromatic images (calculated with 300 × 300 pixels) within seconds. Complex visibilities and fluxes can be directly derived from these monochromatic images. FRACS neglects scattering, an approximation well suited to interpret spectro-interferometric observations in the thermal IR. Indeed, compared to absorption of light from dust, scattering can be neglected in the IR and beyond (Paper I; Lamers & Cassinelli 1999).
To analyse the VLTI/MIDI data of CPD-57° 2874 we adopted the same parametrised description of a sgB[e] (central star and dusty CSE) as given in Paper I. Below we summarise the main equations of this description in axis-symmetric spherical coordinates, and define the free parameters used in the model-fitting.
We assume the specific intensity from the central regions of the star to be a power-law with spectral index α and level at a fiducial wavelength λ0 = 10 μm: (1)This emission from the central region includes a contribution from the stellar photosphere and from the continuum radiation (free-free and free-bound) of the circumstellar ionised gas. The spectral index α is sensitive to the nature of the central source. 4 and Felli & Panagia (1981) give theoretical values for the spectral index for spherical envelopes: α ≃ 4 for a blackbody and α ≃ 2.6 for a fully ionised gas (free-free emission) with an electron density proportional to r-2. Their estimates are valid within the Rayleigh-Jeans domain of the spectrum, which fits to our case when we consider the hot central parts of a sgB[e] in the mid-IR.
A radius Rs = 54 R⊙ was adopted for the central region. This value is used simply as a scaling factor and to convert to the observed 10 μm flux from the central region: (2)We note that at distances of a few kpc the central regions of a sgB[e] are not resolved by VLTI/MIDI and can thus be considered as point sources.
For the gas number density we adopt the bi-modal distribution used by Carciofi et al. (2010) to study another B[e] star. Similar density distribution descriptions were adopted for Be stars (e.g. Stee et al. 1995). The adopted distribution is motivated by the two-wind scenario proposed by Zickgraf et al. (1985) and assumes a fast polar wind, a slow equatorial outflow, and a latitude-dependent mass loss rate. The number density of dust grains is therefore given by (3)where (r,θ) are the radial coordinate and co-latitude, and nin is the dust grain number density at θ = 90° and at r = Rin, which is the inner dust radius where dust starts to survive. Dust is confined between Rin and the outer dust radius Rout. The value of Rout cannot be determined from the VLTI/MIDI data and has been fixed to a high value: 750 AU (the exact value does not affect our results).
The parameter A1 controls the ratio between the equatorial and polar mass loss rates. More precisely, (1 + A1) specifies the ratio between the mass loss rate per unit solid angle in the equator and the pole. In the bi-modal scenario, the poles are assumed to be much less dense than the equator, therefore A1 must have a large value (≳10). Our models indicate that for a wide range of values this parameter does not have a strong influence on the results because it can always be compensated by nin (see also discussion in Paper I). Therefore, we have arbitrarily fixed (1 + A1) to 50.
The parameter A2 indicates how much faster the polar wind is compared to the slow equatorial wind. (1 + A2) is the equatorial to polar terminal velocity, i.e., v∞(90°)/v∞(0°). This parameter is also quite uncertain and in principle can assume values ranging from ~1 to ~100. We kept A2 as a free parameter, although it is not well constrained from the observations as shown in the next sections.
Finally, parameter m controls how fast the mass loss (and consequently the density) drops from the equator to the pole. Defining the disc opening angle Δθd as the latitudinal range within which the mass loss rate is higher then half its equatorial value, we have (4)High m values correspond to thinner regions of high density around the equator. These dense, slowly flowing disc-like regions around the equatorial plane of sgB[e] stars provide favourable conditions for dust to form and survive. Different approaches exist to define regions of dust formation (e.g. Carciofi et al. 2010). Here we adopt the relatively simple assumption where dust is allowed to exist only within the disc opening angle, i.e., at co-latitudes between 90°−0.5Δθd and 90° + 0.5Δθd.
The dust grain opacity was calculated in the Mie theory (Mie 1908) for silicate dust and for a dust size distribution following the commonly adopted MRN (Mathis et al. 1977) power-law ∝ a-3.5, where a is the dust grain radius. The Mie absorption cross sections are computed from the optical indices of astronomical silicate (Draine & Lee 1984; see also Paper I). One possibility to reproduce the absence of a silicate feature in the N-band spectrum of CPD-57° 2874 is to have relatively large grain sizes. We thus used grain radii ranging from a = 0.5 to 50 μm, which can relatively well reproduce the observed spectrum (Fig. 4). We checked that ignoring scattering remains a valid assumption for this dust distribution with large grains. By neglecting the dust albedo the visibilities and fluxes are affected by only a few percent (≲3.5%) within the N-band (further details in Paper I).
The temperature structure of the dusty CSE is given by (5)where Tin is the dust temperature at the disc inner radius Rin, i.e., the dust sublimation temperature. To be consistent with our choice of dust composition we require that Tin ≤ 1500 K. The coefficient γ is expected to assume values ≲1.
Finally, because OLBI is sensitive to the projection of the object’s intensity distribution onto the sky, there are two angles related to this projection:
-
the inclination of the disc plane towards the observer i (0° for pole-on view and 90° for equator-on view);
-
the position angle (from North to East) of the maximum elongation of the sky-projected disc PAd. This angle is defined for i ≠ 0°.
Thus, the 10 free parameters (nfree) of the model are: , α, Tin, γ, Rin, i, PAd, A2, nin, and m.
5. Model-fitting with FRACS
Here we use FRACS with the parametrised sgB[e] description defined in the last section in order to interpret the VLTI/MIDI observations of CPD-57° 2874 through a model-fitting procedure.
To ensure spectrally independent observations for the model-fitting we decided to consider one data point every ≃0.5 μm between 8 μm and 13 μm. This step approximately corresponds to twice the used spectral resolution width (Δλ = λ/R). Additionally, to avoid poor visibility calibration owing to the Earth signature of ozone around 9.6 μm we have not included observations in this spectral region in our analysis. Finally, the same spectral sampling (nλ = 10 wavelengths points) was adopted for the VLTI/MIDI visibilities and spectrum. This choice also provides faster calculations because it is not necessary to compute model images at too many wavelengths.
We have performed a χ2 minimisation simultaneously on the VLTI/MIDI visibilities and fluxes using a Levenberg-Marquardt (LM) algorithm (Markwardt 2008). In order to treat the visibilities and fluxes on the same level (similar weights) we have minimised a χ2 like quantity defined as (see further details in Paper I): (6)where and Vj,k are the observed and modelled visibility modulus wavelength index j and for baseline index k, and Fj are the observed and modelled mid-IR fluxes. σV,j,k and σF,j are the estimated errors on the visibilities and fluxes.
The starting parameter values for the fit were determined from physical considerations of the CSE and from the previous results from DS07. Below we consider the reduced χ2 defined by , where nfree = 10. The LM algorithm stops when the relative decrease in is less then 10-3. For the CPD-57° 2874 data, the LM algorithm reaches the minimum () in a few hours (≃2–3 h) on a single CPU.
Figure 2 shows the intensity map of the model corresponding to (best-fit model) for our distance estimate of 1.7 kpc, which also corresponds to the lowest . The visibilities and fluxes for the best-fit model are shown, together with the observations, in Figs. 3 and 4. These plots show that the model well reproduces most observations within their uncertainties for both adopted distances (1.7 and 2.5 kpc). In particular the slightly curved shape of the visibilities is well reproduced by FRACS. The models indicate that this curved shape is probably caused by the combined fact (1) that the intensity maps have different relative contributions from the central source and from the dusty CSE at different wavelengths, (2) that the optical properties of the adopted dust grains are wavelength-dependent even if there is no strong silicate feature seen in the spectrum, and (3) that the angular resolution significantly changes along the observed wavelengths.
The model parameters at and their uncertainties are listed in Table 4. The derived parameters are almost independent of the adopted distance, except of course for those scaling with the distance. The uncertainties of the parameters have been estimated from maps calculated with 21 × 21 points for each pair of free parameters (45 pairs) and centred on the position. All 45 maps are shown in Fig. 5 for d = 1.7 kpc. These maps show that the space presents a well defined , without showing several local minima in the parameter domain explored. Additionally, they provide visual and direct information on the behaviour of the model parameters in the vicinity of , revealing, for instance, potential correlations between certain parameters.
We have estimated the parameter uncertainties in a conservative way by searching for the maximum parameter extension in all maps corresponding to , where Δχ2 = 1 (see contours in Fig. 5). This choice of Δχ2 sets a lower limit confidence region of ≃60% to the parameter uncertainties. This limit results from two extreme assumptions about the data:
-
data points per baseline are completely dependent (correlated):because the same set of stars is used to calibrate all visibilities of agiven baseline, we can consider a limiting case where all thesevisibilities are correlated. This assumption implies that only 10independent visibility observations are available (thiscorresponds to the number of baselines). The flux at each spectralchannel can still be considered to be independent. Thispessimistic assumption leads to the lower limit of ≃60% to the formal confidence level for Δχ2 = 1, corresponding to only 20 independent observations (10 baselines and 10 fluxes),
-
all data points are completely independent (uncorrelated): an upper limit of ≃100% of formal confidence level is obtained if we assume that all data points are independent. Then the uncertainties derived from Δχ2 = 1 are very conservative (overestimated).
Hence the parameter uncertainties given in Table 4 correspond to a confidence level of at least 60%, but most probably they are somewhat overestimated.
In the next section we present a physically motivated discussion of the derived model parameters of CPD-57° 2874.
Fig. 2 Intensity map of CPD-57° 2874 at 10 μm for the best-fit FRACS model obtained for a distance d = 1.7 kpc (see Table 4). The image scale is in log of the specific intensity . |
Fig. 3 VLTI/MIDI visibilities of CPD-57° 2874 (circles) and the best-fit FRACS visibilities obtained from a χ2 minimisation for a distance d = 1.7 kpc (solid curve) and d = 2.5 kpc (dashed curve). The model-fitting was performed simultaneously on the visibilities and spectral flux. The visibilities effectively used for the fit are shown as filled circles together with the corresponding visibility error bars. The good quality of this fit is reflected by a reduced of ≃0.55 for both distances (see details in Table 4). |
Best-fit model parameters and uncertainties derived for CPD-57° 2874 from a χ2 minimisation.
6. Discussion
6.1. Geometrical parameters (PAd, Rin, and i)
Fig. 4 VLTI/MIDI flux of CPD-57° 2874 (thick solid grey curve) and the ± 20% adopted uncertainty (dots). The thin solid and dashed curves are the best-fit model fluxes for assumed distances of d = 1.7 kpc and d = 2.5 kpc, respectively. The wavelengths used for the fit are those from Fig. 3. |
Let us first compare the derived geometrical parameters (PAd, Rin, and i) with those previously obtained by DS07 from elliptical Gaussian models fitted on a sub-set of the VLTI/MIDI data used here.
The geometrical parameter PAd can be directly compared with the major-axis position angle of the ellipse previously determined by DS07 (≃143°−145°). As expected, the two estimates of PAd are identical within their error bars.
Because the bulk of the thermal IR emission comes from the internal regions of the disc, one can expect the inner dust angular radius (Rin/d) to be comparable to (or slightly smaller than) the major-axis half width at half maximum (HWHM) of an elliptical Gaussian. Indeed, the Rin/d derived here agrees with the major-axis HWHM (=0.5FWHM) given by DS07: 4.5 < HWHM (mas) < 8.0.
Domiciano de Souza et al. (2007) estimated a CSE viewing angle i ~ 30°−60° from the minor- to major-axis ratio of the fitted elliptical Gaussian model. This estimate agrees fairly well with the more precise determination of this parameter given here.
The comparison of these parameters shows on one hand that the geometrical parameters obtained with FRACS agree with those obtained from a simpler approach using analytical models. On the other hand, this comparison clearly shows that, at the cost of somewhat higher but similar computing times, FRACS gives us access to physical parameters of CPD-57° 2874 that cannot be extracted from simple geometrical analytical models.
6.2. Continuum emission from the central source ( and α)
In contrast to what may be initially expected, we show below that the central source (star and continuum emission components such as free-free, free-bound) contributes to almost half of the total mid-IR radiation of CPD-57° 2874. The total 10 μm flux of the best model for d = 1.7 kpc is Ftot = 7.9 × 10-13 Wm-2 μm-1 (see Fig. 4). The contribution of the dust CSE alone to the 10 μm flux computed with FRACS is Fd = 4.4 × 10-13 Wm-2 μm-1. It follows therefore that the 10 μm flux from the CSE is ≃56%Ftot and from the central source is ≃44%Ftot. Similar results are obtained for d = 2.5 kpc.
We note that in the particular case of CPD-57° 2874 where i ≃ 60° and where the dust is confined in a relatively narrow disc (opening angle ~7°), these relative flux contributions can be directly obtained from (derived from ; see Table 4). This is valid if there is no absorption of the central regions by the CSE.
Although the uncertainties on the spectral index α are relatively high (≃50%), the derived value (=2.4) suggests an important contribution of free-free continuum radiation from the central regions (Felli & Panagia 1981).
6.3. Temperature structure of the dusty CSE (Tin and γ)
The dust temperature at the inner radius Tin is found to be ≃1500 K (the imposed upper limit for the fit), which is consistent with the definition of Tin itself and with the chosen silicate dust composition. The large upper limit uncertainties in Tin and the fact that the best-fit Tin is 1500 K for d = 2.5 kpc indicate that a Tin slightly higher than 1500 K could still be compatible with the observations. A higher Tin value is consistent with, for instance, different dust compositions. However, we prefer to keep our choice of dust composition and to have a relatively large upper limit uncertainty in Tin, since the present observations do not provide strong constraints on the exact dust composition.
The derived value for the coefficient of the temperature profile γ indicates an almost linear decrease of the disc temperature as a function of r. The steepness of this temperature profile lies between those expected for a non irradiated, adiabatically cooling disc (T ∝ r−4/3) and a reprocessing disc (T ∝ r−3/4) (cf. Porter 2003). This implies that a non negligible part of the reprocessed radiation from the inner parts of the disc escapes without being re-absorbed, so that the disc cools down faster than highly optically thick discs, without of course reaching the limit of a purely adiabatic cooling.
6.4. Parameters related to the density law (nin, A2, and m)
Our results suggest that the observed mid-IR visibilities and fluxes cannot strongly constrain each individual parameter related to the density law: nin, A2, and m. The χ2 maps show a significant correlation between these three parameters, indicating that they are degenerated for the available VLTI/MIDI observations. Even if their uncertainties are significant, the values of these parameters at suggest
-
a low inner density nin corresponding to relatively low CSE optical depth in the mid-IR (≲0.2 along the line of sight around 10 μm);
-
a polar-to-equatorial terminal velocity ratio v∞(0°)/v∞(90°) = 1/(1 + A2) ≃ 50, compatible with the values found in the literature (e.g. Zickgraf 2003);
-
a high value for m, translating into a quite narrow opening angle (≤10°) for the dust disc.
Because nin, A2, and m are not well constrained, we fitted the observations by fixing these parameters to their values in Table 4 in order to investigate their influence on the remaining parameters. We have also fixed the Tin to 1500 K. The fit was performed for d = 1.7 kpc, starting from slightly different values from those in Table 4 . The χ2 and values obtained for the free parameters (, α, γ, Rin, i, PAd) are essentially the same as in Table 4 (differences are only a small fraction of the parameter uncertainties). We have also checked that the uncertainties on the other parameters are not affected by the fact that nin, A2, and m are not well constrained.
6.5. Data analysis from a model with fewer free parameters
Thanks to the data analysis performed here we found that some parameters of the sgB[e] model adopted for CPD-57° 2874 cannot be well constrained from the available VLTI/MIDI data. Of course this is not necessarily a general conclusion because it depends on the nature of the studied target and on the spectro-interferometric data available.
The results from our analysis of CPD-57° 2874 indicate that it is justified to consider a simplified version of the model described in Sect. 4. As shown in the previous section, the parameters related to the density law are the less constrained by the data. Let us thus consider an alternative density law where the number density of dust grains is given by (7)The parameters A1, A2, and m are not present in this simpler density prescription. Only nin and Δθd are necessary to define the density structure. Based on our previous results we have also fixed Tin to 1500 K. The number of free parameters is thus reduced from 10 to 8, namely, , α, γ, Rin, i, PAd, nin, and Δθd.
Following the previous procedure we have performed a χ2 minimisation on the visibilities and fluxes using a LM algorithm. The results are shown in Table 5. Note that the best-fit parameters completely agree (well within the uncertainties) with the previous values obtained with the initial model (Table 4).
Only nin shows an important difference compared to the previously tested model. Contrarily to the previous density law (Eq. (3)), the new one (Eq. (7)) assumes that the density is constant along a given r inside the dust disc. In order to obtain similar fluxes and optical depths as before, nin has to be somewhat smaller than the previous value.
Table 5 also provides the parameter uncertainties estimated with the LM algorithm, which are smaller than those estimated from the maps (given in Table 4). These smaller errors appear because all data points are assumed to be independent and no covariance matrix is used in the error estimations. The estimated errors from LM in Table 5 can then be considered as lower limits, while the errors from the maps can be considered as upper limits to the parameter uncertainties (see discussion in Sect. 5).
7. Conclusions
The dusty CSE of the Galactic sgB[e] CPD-57° 2874 was spatially resolved thanks to mid-IR spectro-interferometric observations performed with the VLTI/MIDI instrument. Several physical parameters and corresponding uncertainties of this star were derived from a χ2 minimisation and from the analysis of χ2 maps. The physical quantities derived include the inner dust radius, relative flux contribution of the central source and of the dusty CSE, dust temperature profile, and disc inclination (refer to Table 4 and Sect. 6 for details).
To our knowledge, this is the first direct determination of physical parameters of the dusty CSE of a B[e] supergiant based on interferometric data and using a model-fitting approach from a χ2 minimisation. This was possible thanks to FRACS, which adopts a parametrised description of the central object and of the dusty CSE combined to a simplified radiative transfer (no scattering). Its main advantage is computation speed (<10 s per monochromatic image with 300 × 300 pixels). Because it is fast, FRACS allows us (1) to explore a large parameter space domain for a true global χ2 minimisation and (2) to more realistically estimate the uncertainties on the best-fit parameters. We would like to recall that contrarily to a model such as FRACS, simple geometrical models do not allow a simple and direct access to physical parameters and uncertainties of the dusty CSE.
Future complementary observations could be included to measure new CSE parameters and/or to reduce the uncertainties and the correlations on some parameters that were not strongly constrained by the VLTI/MIDI observations alone. Consistently with the domain of validity of FRACS, these complementary data should be obtained at wavelengths above the near-IR, where the dust starts to contribute to a significant amount of the stellar flux (UV to visible observations cannot be consistently modelled by FRACS). For example, the closure phase information from
the MATISSE beam combiner (second generation instrument on VLTI, (Lopez et al. 2009; Wolf et al. 2009) could inform on the disc inclination and opening angle, which are parameters that can influence the symmetry of the mid-IR intensity distribution. Also, near-future interferometric observations in the millimetre with ALMA will be more sensitive than MIDI and MATISSE to the thermal emission from the colder regions of the dusty CSE. Thus, these observations will probably better constrain the temperature profile of the dust and also provide direct information on the structure and actual size of the CSE, allowing also for mass estimates. Moreover, this information on the colder dust is very important for the study of the evolutionary history of sgB[e] stars (e.g. mass and angular momentum losses, chemical enrichment and interaction with the close interstellar medium).
Online material
Fig. 5 Reduced χ2 maps from VLTI/MIDI observations of CPD-57° 2874 for all 45 combinations of free parameters of our model (FRACS). The maps are centered on the (=0.54) position and correspond to the distance of 1.7 kpc. Contours are drawn for , with Δχ2 = 0.3,1,3. The parameters are ordered as in Table 4. The map scale is given in logarithmic units. Details on the model-fitting procedure are given in Sect. 5. |
Fig. 5 continued. |
Fig. 5 continued. |
The MIA+EWS software package is available at http://www.mpia-hd.mpg.de/MIDISOFT/ and http://www.strw.leidenuniv.nl/~nevec/MIDI/index.html
Acknowledgments
This research used the SIMBAD and VIZIER databases at the CDS, Strasbourg (France), and NASA’s ADS bibliographic services. M.B.F. acknowledges financial support from the Programme National de Physique Stellaire (PNPS-France) and the Centre National de la Recherche Scientifique (CNRS-France). M.B.F. also acknowledges Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq-Brazil) for the post-doctoral grant. We thank the CNRS for financial support through the collaboration program PICS France-Brazil. We also thank the referee for his useful and constructive comments, which helped us to improve the quality of this work.
References
- Allen, C. W. 1973, in Astrophysical Quantities, University of London (The Athlone Press), 266 [Google Scholar]
- Bjorkman, J. E. 1998, B[e] stars, 233, 189 [NASA ADS] [Google Scholar]
- Carciofi, A. C., Miroshnichenko, A. S., & Bjorkman, J. E. 2010, ApJ, 721, 1079 [NASA ADS] [CrossRef] [Google Scholar]
- Domiciano de Souza, A., Driebe, T., Chesneau, O., et al. 2007, A&A, 464, 81 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Domiciano de Souza, A., Kervella, P., Bendjoya, P., & Niccolini, G. 2008, A&A, 480, L29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Draine, B. T., & Lee, H. M. 1984, ApJ, 285, 89 [Google Scholar]
- Felli, M., & Panagia, N. 1981, A&A, 102, 424 [NASA ADS] [Google Scholar]
- Jaffe, W. J. 2004, in New Frontiers in Stellar Interferometry, ed. W. A. Traub, Bellingham, WA: The International Society for Optical Engineering, Proc. SPIE Conf., 5491, 715 [Google Scholar]
- Lamers, H. J. G. L. M., & Pauldrach, A. W. A. 1991, A&A, 244, L5 [NASA ADS] [Google Scholar]
- Lamers, H. J. G. L. M., Zickgraf, F.-J., de Winter, D., Houziaux, L., & Zorec, J. 1998, A&A, 340, 117 [NASA ADS] [Google Scholar]
- Lamers, H. J. G. L. M., & Cassinelli, J. P. 1999, Introduction to Stellar Winds, ed. J. P. Lamers, H. J. G. L. M., & Cassinelli [Google Scholar]
- Leinert, Ch., van Boekel, R., Waters, L. B. F. M., et al. 2004, A&A, 423, 537 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Lopez, B., Lagarde, S., Wolf, S., et al. 2009, Science with the VLT in the ELT Era, 353 [Google Scholar]
- Magalhães, A. M. 1992, ApJ, 398, 286 [NASA ADS] [CrossRef] [Google Scholar]
- Markwardt, C. B. 2008, in Proc. Astronomical Data Analysis Software and Systems XVIII, Quebec, Canada, ASP Conf. Ser., ed. D. Bohlender, P. Dowler, & D. Durand, 411, 251 [Google Scholar]
- Mathis, J. S., Rumpl, W., & Nordsieck, K. H. 1977, ApJ, 217, 425 [NASA ADS] [CrossRef] [Google Scholar]
- McGregor, P. J., Hyland, A. R., & Hillier, D. J. 1988, ApJ, 324, 1071 [NASA ADS] [CrossRef] [Google Scholar]
- Meynet, G., & Maeder, A. 2006, Stars with the B[e] Phenomenon, 355, 27 [NASA ADS] [Google Scholar]
- Mie, G. 1908, Ann. Phys., 25, 377 [CrossRef] [Google Scholar]
- Millour, F., Chesneau, O., Borges Fernandes, M., et al. 2009, A&A, 507, 317 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Miroshnichenko, A. S. 2007, ApJ, 667, 497 [NASA ADS] [CrossRef] [Google Scholar]
- Niccolini, G., Bendjoya, P., & Domiciano de Souza, A. 2010, A&A, 525, A21 (Paper I) [Google Scholar]
- Panagia, N., & Felli, M. 1975, A&A, 39, 1 [NASA ADS] [Google Scholar]
- Pelupessy, I., Lamers, H. J. G. L. M., & Vink, J. S. 2000, A&A, 359, 695 [NASA ADS] [Google Scholar]
- Stee, P., de Araujo, F. X., Vakili, F., et al. 1995, A&A, 300, 219 [NASA ADS] [Google Scholar]
- Vink, J. S., Davies, B., Harries, T. J., Oudmaijer, R. D., & Walborn, N. R. 2009, A&A, 505, 743 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Wolf, S., Lopez, B., Jaffe, W., et al. 2009, Science with the VLT in the ELT Era, 359 [Google Scholar]
- Zickgraf, F.-J. 2003, A&A, 408, 257 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Zickgraf, F.-J., Wolf, B., Stahl, O., Leitherer, C., & Klare, G. 1985, A&A, 143, 421 [NASA ADS] [Google Scholar]
All Tables
Summary of VLTI/MIDI observations of CPD-57° 2874: data set index, date, Coordinated Universal Time (UTC) of observation, baseline configuration, projected baseline length and position angle.
Observation log and angular diameters of calibrators (values from DS07) used to derive the calibrated N-band visibilities of CPD-57° 2874.
Equivalent widths (EW) of the Na i absorption lines (λλ 5890 Å, 5896 Å), obtained from our FEROS data.
Best-fit model parameters and uncertainties derived for CPD-57° 2874 from a χ2 minimisation.
All Figures
Fig. 1 uv-Plane coverage: projected baselines (length and position angle) for the VLTI/MIDI observations of CPD-57° 2874 (further details are given in Table 1). |
|
In the text |
Fig. 2 Intensity map of CPD-57° 2874 at 10 μm for the best-fit FRACS model obtained for a distance d = 1.7 kpc (see Table 4). The image scale is in log of the specific intensity . |
|
In the text |
Fig. 3 VLTI/MIDI visibilities of CPD-57° 2874 (circles) and the best-fit FRACS visibilities obtained from a χ2 minimisation for a distance d = 1.7 kpc (solid curve) and d = 2.5 kpc (dashed curve). The model-fitting was performed simultaneously on the visibilities and spectral flux. The visibilities effectively used for the fit are shown as filled circles together with the corresponding visibility error bars. The good quality of this fit is reflected by a reduced of ≃0.55 for both distances (see details in Table 4). |
|
In the text |
Fig. 4 VLTI/MIDI flux of CPD-57° 2874 (thick solid grey curve) and the ± 20% adopted uncertainty (dots). The thin solid and dashed curves are the best-fit model fluxes for assumed distances of d = 1.7 kpc and d = 2.5 kpc, respectively. The wavelengths used for the fit are those from Fig. 3. |
|
In the text |
Fig. 5 Reduced χ2 maps from VLTI/MIDI observations of CPD-57° 2874 for all 45 combinations of free parameters of our model (FRACS). The maps are centered on the (=0.54) position and correspond to the distance of 1.7 kpc. Contours are drawn for , with Δχ2 = 0.3,1,3. The parameters are ordered as in Table 4. The map scale is given in logarithmic units. Details on the model-fitting procedure are given in Sect. 5. |
|
In the text |
Fig. 5 continued. |
|
In the text |
Fig. 5 continued. |
|
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.