Free Access
Issue
A&A
Volume 566, June 2014
Article Number A88
Number of page(s) 11
Section Stellar atmospheres
DOI https://doi.org/10.1051/0004-6361/201323264
Published online 19 June 2014

© ESO, 2014

1. Introduction

The motivation for our study is to improve our understanding of the circumstellar environment of asymptotic giant branch stars (AGBs) close to the photosphere, to obtain estimates about their fundamental parameters, and to locate them in the Hertzsprung-Russell (HR) diagram. The location of the stars in the HR diagram is very important for calibrating stellar evolutionary models for intermediate-mass stars.

Interferometric techniques at visible and IR wavelengths are important for resolving the stellar disk to better understand the circumstellar environment (Quirrenbach et al. 1993; Perrin et al. 2004). Recent studies with VLTI/AMBER and VLTI/MIDI have provided information about the pulsation and the mass-loss of AGB stars (Ohnaka et al. 2006, 2007; Wittkowski et al. 2007; Chiavassa et al. 2010; Karovicova et al. 2011, 2013) and about the structure of the molecular distribution in AGB stars (Wittkowski et al. 2008, 2011; Martí-Vidal et al. 2011).

Table 1

VLTI/AMBER observations.

Quirrenbach et al. (1993, 2001) studied the TiO band (around 712 nm) in the atmosphere of cold giant stars (spectral type M). Their interferometric observations were made with two filters, one centered on the TiO band, and the other on the continuum close to that band. They observed an increase of the size of the star corresponding to the TiO band with respect to the size in the continuum. After fitting the PHOENIX models (Hauschildt & Baron 1999) to their data, they concluded that the diameter ratio between the TiO band and the continuum agreed with models computed for a mass of 0.5 M; but the evolutionary models predict for these stars masses of about 5 M. This disagreement might be explained by the existence of a transition zone at the base of the stellar wind (Tsuji 2008), which could provide sufficient opacity in the TiO band to make the AGB larger than the size predicted by the PHOENIX model.

Martí-Vidal et al. (2011) observed RS Cap (AGB star of spectral type M6/M7III) with the VLTI/AMBER instrument in the K band. They found that the apparent size of the star increased around 12% in the CO band (2.29–2.47 μm). The fit to the data with MARCS models (Gustafsson et al. 2008) was reasonable, although the lower visibilities in the CO band were not reproduced. These authors added an ad hoc spherical water envelope around the star (Perrin et al. 2004) that made the synthetic visibilities and the observations in the CO bands appear consistent.

Cruzalèbes et al. (2013) observed sixteen red giants and supergiants with VLTI/AMBER over a two-year period. They used MARCS models to fit their data. Their estimates of the angular diameters were moderately dependent on the variation of the model input parameters Teff, log (g), and ξturb. Eight of these sources were studied for the first time but the others had been studied earlier with Long-Baseline Interferometry (LBI), and the angular diameter estimates obtained with both methods were similar.

Cusano et al. (2012) studied five giant stars while investigating planet formation around stars more massive than the Sun. They estimated the uniform disk (UD) and limb-darkened (LD) angular diameters and the effective temperatures of these sources. The measurements of both angular diameters (UD and LD) were consistent within 1.5σ, the differences being smaller than 0.8%. Their estimates were also consistent with the values derived by da Silva et al. (2006).

In this paper we study a sample of cool giant stars with VLTI/AMBER. We locate our targets in the HR diagram and compare our results with those of the red supergiant stars studied in Arroyo-Torres et al. (2013). The remainder of this paper is structured as follows: in Sect. 2, we describe our AMBER observations and the data reduction, in Sect. 3, we present the PHOENIX model used, in Sect. 4, we report and discuss our results. Finally, we present our conclusions in Sect. 5.

2. Observations and data reduction

We observed ϵ Oct (M5 III), β Peg (M2.5 II-III), NU Pav (M6 III), ψ Peg (M3 III), and γ Hya (G8 III) with the ESO Very Large Telescope Interferometer (VLTI), using three of the auxiliary telescopes of 1.8 m diameter, and the AMBER instrument (Astronomical Multi-BEam combineR) with the external fringe tracker FINITO (Petrov et al. 2007). We observed in medium-resolution mode (R ~ 1500) in the K-2.3 μm band (γ Hya has been also observed in band K-2.1). We scheduled our observations as sequences of cal-sci-cal (cal is calibrator and sci is our target), with five scans for each of them. The integration time (DIT) of each frame was 0.2s (for ϵ Oct, β Peg, NU Pav, and ψ Peg) and 1.0s (for γ Hya). The stars ϵ Oct and β Peg were each observed with two different arrays. In Table 1, we show information about our observations and the calibrator used for each target. In Table 2, we show the calibrators used for our observations, selected from the ESO Calibration Selector CalVin, in turn based on the catalog of Lafrasse et al. (2010).

Table 2

Calibration sources.

Raw visibilities and closure phases were obtained from our AMBER data using the latest version of the amdlib data reduction package (Tatulli et al. 2007; Chelli et al. 2009). First, we removed the bad-pixel map and corrected for the flat contribution. Then we calculated the pixel-to-visibility matrix (P2VM) to calibrate our data for the instrumental dispersive effects, and obtained the interferometric observables. Next, we appended all scans of the same source taken consecutively and selected and averaged the resulting visibilities of each frame using appropriate criteria. In our case, the criteria were based on the flux (we selected all frames with flux densities three times higher than the noise) and on the signal-to-noise ratio (S/N). We only used 80% of the remaining frames with the best S/N1.

Using in-house developed IDL (Interactive Data Language) scripts we performed the absolute wavelength calibration by correlating the AMBER flux spectra with a reference spectrum, that of the star BS 4432 (spectral type K4.5 III, similar to our calibrators; Lançon & Wood 2000). A relative flux calibration of the targets was performed by using the instrumental response, estimated from the calibrators and the BS 4432 spectrum. Finally, calibrated visibility spectra were obtained by using the average of two transfer functions obtained from calibrator observations before and after each science target observation: Vscic=Vscim/0.5(𝒯cal1+𝒯cal2)\hbox{$V^c_{\rm sci}=V^m_{\rm sci}/0.5(\mathcal{T}_{\rm cal_1}+\mathcal{T}_{\rm cal_2})$}. In the case of ϵ Oct and NU Pav, we only used one calibrator, because the visibilities of the other were not of sufficient quality. The errors of the calibrated visibilities were estimated by error propagation. For each calibrator, the error of the transfer function \hbox{$\Delta \mathcal{T}$} was calculated as Δ𝒯=Δ𝒯A2+Δ𝒯B2+0.052\hbox{$\Delta \mathcal{T}=\sqrt{\Delta \mathcal{T}^2_{A}+\Delta \mathcal{T}^2_{B}+0.05^2}$}. The first term was obtained as Δ𝒯A=ΔVcalm/Vexp\hbox{$\Delta \mathcal{T}_{A}=\Delta V^{m}_{\rm cal}/V_{\rm exp}$} where ΔVcalm\hbox{$\Delta V^{m}_{\rm cal}$} were the uncertainties of the measured calibrator visibilities and Vexp = 2J1(πθB/λ)/(πθB/λ), with θ the angular diameter adopted for the calibrator. The second term was defined as \hbox{$\Delta \mathcal{T}_{B}=|\mathcal{T}_{\rm cal_1}-\mathcal{T}_{\rm cal_2}|/2$}. The value of 0.05 was a systematic error adopted. This term is relevant when only one calibrator is available (then we make \hbox{$\mathcal{T}_{\rm cal_1}=\mathcal{T}_{\rm cal_2}$} for practical reasons) or when by chance \hbox{$\mathcal{T}_{\rm cal_1}$} and \hbox{$\mathcal{T}_{\rm cal_2}$} are very similar.

3. Modeling: PHOENIX

Our goals are to estimate the angular diameters of our stars and then derive their fundamental parameters. For the purpose of estimating the stellar angular diameter, we compared our observables with those provided by PHOENIX model atmospheres (version 16.03, Hauschildt & Baron 1999, which incorporates a model with a hydrostatic atmosphere, local thermodynamic equilibrium and spherical geometry). We used the grid corresponding to 1 M, that was also used by Arroyo-Torres et al. (2013). This grid includes effective temperatures between 2500 K and 8000 K in steps of 100 K, surface gravities between log (g) = −0.5 and log (g) = 4.0 in steps of 0.5 (in cgs units), and solar metallicity. Within these grids, the model provides intensity profiles for different angles of the star. To compare these models with our data, we need the flux integrated over the stellar disk and the visibility values for the baseline used. To obtain this flux, we tabulated model intensity profiles at 64 viewing angles for wavelengths 21 000 Å to 25 000 Å in steps of 0.01 Å. Then, we averaged the monochromatic intensity profiles to match the spectral channels of the individual observations and computed the data of the model for that match (see a full description of the procedure in Wittkowski et al. 2003).

thumbnail Fig. 1

Left: from top to bottom, observed (black) squared visibility amplitudes, UD diameters predicted from our data (blue) and from the best-fit PHOENIX model (red), closure phases in degrees, and normalized flux of β Peg obtained on 2012 Jun. 25. Right: same as left for data obtained on 2012 Aug. 09.

3.1. Fitting the PHOENIX model to the observations

After obtaining the grids, and using estimates of the distance and the bolometric flux (see Sect. 4.2), we proceeded to fit the PHOENIX model to our data to obtain the fundamental parameters of our stars. This iterative process was explained in Arroyo-Torres et al. (2013). In this case, we selected the initial Teff given by Dyck et al. (1998) and van Belle et al. (2009): Teff = 3330 K for ϵ Oct; Teff = 3890 K for β Peg; Teff = 3248 K for NU Pav; Teff = 3475 K for ψ Peg; and Teff = 5087 K for γ Hya. We adopted an initial surface gravity of log (g) = 0.0 throughout.

After all the iterations, the final values for the PHOENIX model are for ϵ Oct, Teff = 3500 K, log (g) = 0.0; for β Peg, Teff = 3800 K, log (g) = 1.0; for NU Pav, Teff = 3500 K, log (g) = 0.0; for ψ Peg, Teff = 3700 K, log (g) = 1.0; and for γ Hya, Teff = 4800 K, log (g) = 2.0. For all five cases, we used a model with solar metallicity, micro-turbulent velocity of 2 km s-1, and mass of 1 M. We note that the structure of the atmosphere is not very sensitive to variations of the mass (Hauschildt et al. 1999). Certainly, any of those structure variations are below the level of the detectability of our interferometer.

4. Results and discussion

Figure 1 and Figs. 2 to 4 show the visibility amplitudes, the closure phases, the uniform disk diameters for the stars ϵ Oct, β Peg, NU Pav, ψ Peg, and γ Hya, respectively, and also the fluxes for all of them but γ Hya. For ϵ Oct and β Peg, we have observations of two different epochs (2012 June – left panel – and 2012 August – right panel), while for ψ Peg, NU Pav, and γ Hya, we only have observations of one epoch (2012 June, 2012 August, 2013 March, respectively).

The visibility curves are shown in the top panels of Fig. 1 and Figs. 2 to 4. In these figures, we see no visibility decrease at the positions of the CO lines (in contrast with the results shown in Martí-Vidal et al. 2011, for RS Cap, or in Arroyo-Torres et al. 2013, for a sample of RSGs). Similarly, these is no drop in the visibility between 2.3 μm and 2.5 μm either. For our sources (except for β Peg), the synthetic visibilities are consistent with our observations. Since the PHOENIX model uses a hydrostatic atmosphere model and a limb-darkened disk, ϵ Oct, NU Pav, ψ Peg, and γ Hya are compatible with a limb-darkened disk and a hydrostatic atmosphere. On the other hand, β Peg shows a decrease of the visibility in the CO lines greater than the synthetic visibilities, and also a weak systematic trend in the range 2.32.5 μm, like that observed in RS Cap.

The uniform disk diameters calculated from our data as a function of wavelength are shown in the second panel of Fig. 1 and Figs. 2 to 4. In the observations from 2012 August, we see a larger scatter in the data, which prevents us from drawing any firm conclusion about the possible larger diameter in the CO lines (the atmospheric conditions were poor these nights). In the case of ϵ Oct, the data from 2012 Jun. 25 show that the uniform disk diameter is constant across the band. Thus, the CO bandheads do not present a larger size than for the continuum of the star. The synthetic visibilities and our observations are consistent. The β Peg data observed in 2012 Jun. 25 show a size in the CO band similar to that predicted by the PHOENIX model, but this source also presents an additional slope in the data that is not predicted by the PHOENIX model. This slope is perhaps due to a layer of H2O that is not present in the model (as in RS Cap, Martí-Vidal et al. 2011). The increase of the angular diameter in the CO region is about 5.3% with respect to the near-continuum bandpass diameter. The NU Pav data have so much scatter that we cannot discern any change in the angular diameter of the region related to the CO lines. The ψ Peg data show a very small increase of the angular diameter in the CO band, similar to the one synthesized in the PHOENIX model. Finally, the uniform disk diameter for γ Hya is constant in both bands (K-2.1 and K-2.3), with the size at the CO bandheads not larger than in the continuum. In summary, the PHOENIX predictions agree well with the data, showing that the simulated atmospheres are as compact as the observed ones, except for β Peg, with a 5.3% size increase clearly visible in the data but not modeled by PHOENIX.

The closure phases are shown in the third panel of Fig. 1 and Figs. 2 to 4. They show low values (20 deg) across the band, indicating little or no deviations from point symmetry. However, as our measurements lie in the first visibility lobe, we cannot exclude asymmetries on scales smaller than the observed stellar disk.

The fourth panel in Figs. 1 and 2 to 3 shows the normalized flux spectra of our targets. For γ Hya we do not have spectra data because there was an error with the data formatting during the observation. In the normalized spectra, we observe a decreasing flux from 2.25 μm and strong absorption lines of CO. The synthetic spectra of the PHOENIX model agree well with our data, in particular in the CO bandheads. This indicates that the opacities of CO lines in cool giant stars are well reproduced by the PHOENIX model (as found previously by Lançon et al. 2007).

We note that the sizes estimated in the CO bandheads of all the stars in our sample match the extension of the CO layers predicted by the PHOENIX hydrostatic models (with the exception of β Peg, for which an additional layer of water vapor may be necessary). In contrast, much more extended CO layers, which cannot be reproduced by hydrostatic models, have been reported in other AGB and RSG stars (Martí-Vidal et al. 2011; Wittkowski et al. 2008, 2012; Arroyo-Torres et al. 2013).

thumbnail Fig. 5

Average (per baseline) of squared visibility amplitudes taken in the continuum bandpass at 2.152.25 μm for ϵ Oct, β Peg, NU Pav, ψ Peg, and γ Hya (from top to bottom) as a function of spatial frequency. For ϵ Oct and β Peg, the graphics show the data of June and August, and for γ Hya show the data of the both setups. The red lines indicate the best-fit UD models and the blue lines the best-fit PHOENIX models. The dashed lines are the maximum and minimum visibility curves, from which we estimated the errors of the angular diameters.

The stars reported here are semiregular late-type giants (ϵ Oct and NU Pav), an irregular pulsating variable (β Peg), a normal giant (ψ Peg), and a low-amplitude variable (γ Hya). These stars are characterized by short and irregular periods, with a low or very low variability amplitude (ΔV between 0.06 and 0.72; see Table 3). However, the stars reported in Arroyo-Torres et al. (2013) and Martí-Vidal et al. (2011) are semiregular giants with much longer periods and variability amplitudes (ΔV between 2.0 and 2.3).

The larger extension of the CO layers in stars with higher variability may be indicative of the important role of strong pulsations for the mechanical transport of the molecular gas, away from the stellar surface. These pulsations trigger the onset of the stellar wind (by radiative pressure on the condensed dust grains) according to the standard model of strong winds in evolved giant stars.

Table 3

Information about the variability.

4.1. Estimate of the angular diameter

After the best PHOENIX models were fitted to our data, we estimated the angular diameter for each source. To estimate the angular diameter, we used the continuum interferometric data near 2.20 μm, free of contamination from the CO band and well reproduced by a limb-darkened disk.

The angular diameter as obtained from the best-fit to the PHOENIX model corresponds to the size of the outermost layer (0% intensity). To estimate the Rosseland angular diameter (corresponding to the layer where the Rosseland optical depth equals 2/3), we multiplied our value of the angular diameter by the ratio between the outermost layer and the Rosseland layer. This ratio was 0.79 (for ϵ Oct and NU Pav), 0.93 (for β Peg and ψ Peg), and 0.97 (for γ Hya). The resulting Rosseland angular diameters and the angular diameters obtained from the UD model are shown in Table 4. Results for ϵ Oct and β Peg were obtained using the two available epochs (2012 June and 2012 August), and results for γ Hya were estimated using both setups (K-2.1 and K-2.3). The errors include statistical and systematic errors caused by calibration uncertainties. We conservatively estimate these errors from the differences between the visibility curves lying at the maximum and minimum of our data.

The angular diameter of β Peg, ψ Peg, and γ Hya were estimated previously by Richichi et al. (2005). Their angular diameter estimates were θUD = 16.19 ± 0.23 mas and θLD = 16.75 ± 0.24 mas for β Peg, θUD = 6.40 ± 0.60 mas for ψ Peg, θUD = 2.96 ± 0.15 mas for γ Hya. These values are compatible with our estimates shown in Table 4.

Table 4

Summary of estimated angular diameters.

Figure 5 shows averaged continuum visibility data as a function of spatial frequency and the visibility curves corresponding to the extreme value covered by our uncertainty in the angular diameter. The visibility errors shown in this figure were estimated as an average of individual errors, since the errors are dominated by systematic effects. The model fit itself was made using all the individual data points, not those averaged. According to the results of Fig. 5, our PHOENIX estimated angular diameter is compatible with our observations. In β Peg, ψ Peg and γ Hya, we have several points sampling up to the first visibility null. Knowing the position of the null gives higher precision to the estimate of the angular diameter.

Table 5

Fundamental parameters of ϵ Oct, β Peg, NU Pav, ψ Peg, and γ Hya.

4.2. Fundamental parameters

Our angular diameter results allow us to obtain estimates of fundamental parameters of the observed stars, namely, the effective temperature, the radius, and the luminosity. The angular diameters are those calculated in Sect. 4.1 and shown in Table 4. The effective temperature is based on the angular diameter and the bolometric flux; the radius is estimated from the angular diameter and the adopted distance; and the luminosity is derived from the bolometric flux and the distance. We used the distance values from Anderson & Francis (2012).

For the bolometric flux of our targets, we used the BVRIJHKL magnitudes from Kharchenko (2001), Morel & Magnenat (1978), and Cutri et al. (2003). We also used the IRAS flux from IRAS (1988). To convert the magnitudes into fluxes, we used the zero values from Skinner (1996) and the 2MASS (Cohen et al. 2003) system. To deredden the flux values we used the color excess method applied to (VK), because photometric colors with longer spectral baselines provide more accurate results. The color excess EVK is calculated as the difference between the observed color (VK)star of our star and the intrinsic color (VK)0, obtained from Ducati et al. (2001) as a function of the spectral type of the star. We converted the EVK into EBV using the following equations for M-type stars given by Fiorucci & Munari (2003): AKAV=0.12AVEBV=3.69,andEBV=1(1AK/AV)3.69EVK.\begin{eqnarray*} &&\frac{A_K}{A_V}=0.12\\ &&\frac{A_V}{E_{B-V}}=3.69 \: \mathrm{, and}\\ &&E_{B-V}=\frac{1}{(1-A_K/A_V)\ast 3.69}\ast E_{V-K}. \end{eqnarray*}The obtained EVK values are summarized in Table 5. After that, we calculated the absorption in all photometric bands (Aλ) by means of the equation Aλ = (αλ + βλEBV) ∗ EBV, where the values of αλ and βλ are also taken from Fiorucci & Munari (2003). Finally, we corrected all the flux values for interstellar extinction and integrated them to obtain the bolometric flux (we used the Newton-Cotes integration method within the IDL program). Table 5 summarizes the fundamental parameters for our targets.

Figure 6 shows effective temperature vs. spectral type for our targets. We also include in the figures the AGB star RS Cap from Martí-Vidal et al. (2011), and the following sets of RSGs: AH Sco, UY Sct, and Kw Sgr from Arroyo-Torres et al. (2013), VY CMa from Wittkowski et al. (2012), Betelgeuse from Ohnaka et al. (2011), and VX Sgr from Chen et al. (2007) and Chiavassa et al. (2010). To calculate the effective temperature of RS Cap we used the angular diameter from Martí-Vidal et al. (2011) and the bolometric flux from Dyck et al. (1998).

For comparison, Fig. 6 includes the calibrations of the effective temperature scale by Dyck et al. (1998) for cool giants stars and van Belle et al. (2009) for cool giants stars and RSG stars. We also show the effective temperature scale by Levesque et al. (2005) only for RSGs. Given our observational uncertainties, the cool giant stars agree with the calibration of these effective temperature scales.

thumbnail Fig. 6

Effective temperature versus spectral type of our sources in this paper, RS Cap from Martí-Vidal et al. (2011), and RSGs in Arroyo-Torres et al. (2013), VY CMa from Wittkowski et al. (2012), Betelgeuse from Ohnaka et al. (2011), and VX Sgr from Chen et al. (2007) and Chiavassa et al. (2010). We also show the scales of Dyck et al. (1998), Levesque et al. (2005), and van Belle et al. (2009).

4.3. HR-diagram

Figure 7 shows the positions of ϵ Oct, β Peg, NU Pav, ψ Peg, and γ Hya in the Hertzsprung-Russell (HR) diagram, together with the evolutionary tracks from the model STAREVOL (Lagarde et al. 2012). We used the model for which the transport processes in radiative zones are performed by thermohaline mixing and rotation-induced mixing. Thermohaline mixing occurs when material of high mean molecular weight lies on top of material of low mean molecular weight, a situation that is unstable against a blob of material moving downward, causing the mixture. This process is developed along the red giant branch in low-mass stars and on the early-AGB in intermediate-mass stars (more information in Lagarde et al. 2012). We also show the position of RS Cap from Martí-Vidal et al. (2011). We represent evolutionary tracks with masses between 1.0 M and 6 M and solar metallicity. All targets are close to the red limit of these tracks (Hayashi limit). ϵ Oct, ψ Peg, γ Hya, and RS Cap are close to the evolutionary tracks with initial masses between 1 M and 3 M (the age of γ Hya is younger than the others sources). β Peg is consistent with tracks of masses between 1.25 M and 4 M and NU Pav with tracks of masses between 2.5 M and 6 M.

thumbnail Fig. 7

Location of ϵ Oct, β Peg, NU Pav, ψ Peg, and γ Hya in the HR diagram using our determination of the fundamental parameters. We also show the position of RS Cap by Martí-Vidal et al. (2011). The positions of the stars are compared with evolutionary tracks from Lagarde et al. (2012) for masses of 1.0 M, 1.25 M, 1.5 M, 2.0 M, 2.5 M, 3.0 M, 4.0 M, and 6.0 M.

Additionally, we compare the position of our stars with the evolutionary tracks from Ekström et al. (2012). We used the evolutionary tracks in Arroyo-Torres et al. (2013) for comparison with a sample of RSGs (AH Sco, UY Sct, and KW Sgr) because they were a good choice for RSGs. In this work, we show both samples for comparison purposes. We also show the position of RS Cap from Martí-Vidal et al. (2011), VY CMa from Wittkowski et al. (2012), Betelgeuse from Ohnaka et al. (2011), and VX Sgr from Chen et al. (2007) and Chiavassa et al.(2010). With these tracks, we observe that our stars are located considerably to the right of the models, unlike with the Lagarde evolutionary tracks. The positions of ϵ Oct, β Peg, and ψ Peg are close to the evolutionary tracks with initial masses of 5 or 7 M. NU Pav is consistent with tracks of masses of 7 or 9 M, γ Hya with mass of 3 M, and RS Cap with masses of 59 M. In all cases the tracks are shown without/with rotation. The position of the RSGs are close to tracks with masses between 20 M and 40 M. As expected, RSGs are much more massive and luminous than cool giants stars for the same effective temperatures.

The STAREVOL model is complementary to the Ekström model for low-mass stars. The Ekström model is computed only until the helium flash at the RGB tip and does not include thermohaline mixing. Both models use almost the same assumptions and input physics: convection, opacity, mass loss, and nuclear reaction rates. The initial abundances are also similar, although the STAREVOL model considers more species (for more information see Lagarde et al. 2012).

thumbnail Fig. 8

Location of ϵ Oct, β Peg, NU Pav, ψ Peg, and γ Hya in the HR diagram using our determination of the fundamental parameters. We also show the position of RS Cap by Martí-Vidal et al. (2011), and the RSGs studied in Arroyo-Torres et al. (2013), VY CMa from Wittkowski et al. (2012), Betelgeuse from Ohnaka et al. (2011), and VX Sgr from Chen et al. (2007) and Chiavassa et al. (2010). The positions of the stars are compared with evolutionary tracks from Ekström et al. (2012) for masses of 3 M, 4 M, 5 M, 7 M, 9 M, 12 M, 15 M, 20 M, 25 M, 32 M, and 40 M. The solid lines are models without rotation, the dashed lines with rotation.

5. Conclusions

Our spectro-interferometric near-infrared observations of ϵ Oct, NU Pav, ψ Peg, and γ Hya show that synthetic visibilities from hydrostatic atmospheric models are consistent with the observations, concluding that their atmospheres can be modeled with a limb-darkened disk. In ϵ Oct, NU Pav, and γ Hya, the uniform disk diameter is constant across the band, and the CO bandheads present a similar size to that of the continuum. On the other hand, the data of ψ Peg show a low increase in the CO band, similar to the one obtained in the model. According to these results, the atmospheres of ϵ Oct, NU Pav, ψ Peg, and γ Hya are compatible with hydrostatic atmospheres and the role of pulsation does not seem to be important. However, the data from β Peg (at least in the 2012 June epoch) show a layer (possibly of H2O) that is not modeled by PHOENIX, but CO bands similar to those modeled with PHOENIX. The uniform disk diameter of the star at the CO band increases about 5.3% with respect to the continuum (less than the 14% increase of diameter observed in RS Cap).

We used the continuum near 2.2 μm, which is free from molecular band contamination, to estimate the angular diameter of the targets (see Table 4). We also estimated fundamental parameters such as the luminosity, Rosseland radius, and temperature (shown in Table 5).

Finally, we located each of our targets in the HR diagram using the effective temperature and the luminosity calculated from the Rosseland angular diameter, the bolometric flux, and the distance. In the HR diagram, we also showed the evolutionary tracks from Lagarde et al. (2012). The positions of the stars in this HR diagram are close to the Hayashi limit. Their positions are close to evolutionary tracks corresponding to stars of initial masses between 1.0 M and 3 M (ϵ Oct, ψ Peg, γ Hya, and RS Cap), between 1.25 M and 4 M (β Peg), and between 2.5 M and 6 M (γ Hya). We also compared the position of our stars with the evolutionary tracks from Ekström et al. (2012). The STAREVOL model fits the positions of our stars in the HR diagram better than the Ekström model. This is probably because the STAREVOL model is designed for low-mass stars on the red giant branch and for intermediate-mass stars on the early-AGB. It is complementary to the Ekström model for low- and intermediate- mass stars.

Online material

thumbnail Fig. 2

Left: from top to bottom, observed (black) squared visibility amplitudes, UD diameters predicted from our data (blue) and from the best-fit PHOENIX model (red), closure phases in degrees, and normalized flux of ϵ Oct obtained on 2012 Jun. 25. Right: same as left, but for data obtained on 2012 Aug. 02.

thumbnail Fig. 3

Left: from top to bottom, observed (black) squared visibility amplitudes, UD diameters predicted from our data (blue) and from the best-fit PHOENIX model (red), closure phases in degrees, and normalized flux of NU Pav obtained on 2012 Aug. 02. Right: same as left, but for data from ψ Peg obtained on 2012 Jun. 16.

thumbnail Fig. 4

Left: from top to bottom, observed (black) squared visibility amplitudes, UD diameters predicted from our data (blue) and from the best-fit PHOENIX model (red), and closure phases in degrees of γ Hya obtained on 2013 Mar. 16. Right: same as left, but obtained with the MR-K 2.3 μm setting.


1

See AMBER Data Reduction Software User Manual; http://www.jmmc.fr/doc/approved/JMMC-MAN-2720-0001.pdf

Acknowledgments

This research has made use of the AMBER data reduction package of the Jean-Marie Mariotti Center and the SIMBAD database operated at CDS, Strasbourg, France. BAT, J.C.G. and J.M.M. acknowledge support by the Spanish Ministry of Science and Innovation through the grants AYA2009-13036-C02-02 and AYA2012-38491-C02-01.

References

  1. Anderson, E., & Francis, C. 2012, VizieR Online Data Catalog: V/137 [Google Scholar]
  2. Arroyo-Torres, B., Wittkowski, M., Marcaide, J. M., & Hauschildt, P. H. 2013, A&A, 554, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Beichman, C. A., Neugebauer, G., Habing, H. J., Clegg, P. E., & Chester, T. J. 1988, Infrared astronomical satellite (IRAS) catalogs and atlases. Vol. 1: Explanatory supplement, 1 [Google Scholar]
  4. Chelli, A., Utrera, O. H., & Duvert, G. 2009, A&A, 502, 705 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Chen, X., Shen, Z.-Q., & Xu, Y. 2007, Chin. J. Astron. Astrophys., 7, 531 [Google Scholar]
  6. Chiavassa, A., Lacour, S., Millour, F., et al. 2010, A&A, 511, A51 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Cohen, M., Wheaton, Wm. A., & Megeath, S. T. 2003, AJ, 126, 1090 [NASA ADS] [CrossRef] [Google Scholar]
  8. Cruzalèbes, P., Jorissen, A., Rabbia, Y., et al. 2013, MNRAS, 434, 437 [NASA ADS] [CrossRef] [Google Scholar]
  9. Cusano, F., Paladini, C., Richichi, A., et al. 2012, A&A, 539, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, The IRSA 2MASS All-Sky Point Source Catalog, NASA/IPAC Infrared Science Archive [Google Scholar]
  11. da Silva, L., Girardi, L., Pasquini, L., et al. 2006, A&A, 458, 609 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Ducati, J. R., Bevilacqua, C. M., Rembold, S. B., & Ribeiro, D. 2001, ApJ, 558, 309 [NASA ADS] [CrossRef] [Google Scholar]
  13. Dyck, H. M., van Belle, G. T., & Thompson, R. R. 1998, AJ, 116, 981 [NASA ADS] [CrossRef] [Google Scholar]
  14. Ekström, S., Georgy, C., Eggenberger, P., et al. 2012, A&A, 537, A146 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Fiorucci, M., & Munari, U. 2003, A&A, 401, 781 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Hauschildt, P. H., & Baron, E. 1999, J. Comput. Appl. Math, 109, 41 [Google Scholar]
  18. Hauschildt, P. H., Allard, F., Ferguson, J., et al. 1999, ApJ, 525, 871 [NASA ADS] [CrossRef] [Google Scholar]
  19. Karovicova, I., Wittkowski, M., Boboltz, D. A., et al. 2011, A&A, 532, A134 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Karovicova, I., Wittkowski, M., Ohnaka, K., et al. 2013, A&A, 560, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Kharchenko, N. V. 2001, KFNT, 17, 409 [Google Scholar]
  22. Lafrasse, S., Mella, G., Bonneau, D., et al. 2010, VizieR Online Data Catalog: II/300 [Google Scholar]
  23. Lagarde, N., Decressin, T., Charbonnel, C., et al. 2012, A&A, 543, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Lançon, A., & Wood, P. R. 2000, A&AS, 146, 217 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Lançon, A., Hauschildt, P. H., Ladjal, D., & Mouhcine, M. 2007, A&A, 468, 205 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Levesque, E. M., Massey, P., Olsen, K. A. G., et al. 2005, ApJ, 628, 973 [NASA ADS] [CrossRef] [Google Scholar]
  27. Martí-Vidal, I., Marcaide, J. M., Quirrenbach, A., et al. 2011, A&A, 529, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Morel, M., & Magnenat, P. 1978, A&AS, 34, 477 [NASA ADS] [Google Scholar]
  29. Ohnaka, K., Driebe, T., Hofmann, K.-H, et al. 2006, A&A, 445, 1015 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Ohnaka, K., Driebe, T., Weigelt, G., & Wittkowski, M. 2007, A&A, 466, 1099 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Ohnaka, K., Weigelt, G., Millour, F., et al. 2011, A&A, 529, A163 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Perrin, G., Ridgway, S. T., Coudé du Foresto, V., et al. 2004, A&A, 418, 675 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Petrov, R. G., Malbet, F., Weigelt, G., et al. 2007, A&A, 464, 1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Quirrenbach, A., Mozurkewich, D., Armstrong, J. T., et al. 1993, ApJ, 406, 215 [NASA ADS] [CrossRef] [Google Scholar]
  35. Quirrenbach, A., Mozurkewich, D., Armstrong, J. T., et al. 2001, BAAS, 33, 882 [NASA ADS] [Google Scholar]
  36. Richichi, A., Percheron, I., & Khristoforova, M. 2005, A&A, 431, 773 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Samus, N. N., Durlevich, O. V., Kazarovets, E. V., et al. 2007, General Catalogue of Variable Stars (Samus+ 2007−2012), Vizier On-line Data catalog: B/gcvs7gcvs_cat [Google Scholar]
  38. Skinner, C. J. 1996, Space Telescope NICMOS Instrument Science Report, 14 [Google Scholar]
  39. Tatulli, E., Millour, F., Chelli, A., et al. 2007, A&A, 464, 29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Tsuji, T. 2008, A&A, 489, 1271 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. van Belle, G. T., Creech-Eakman, M. J., & Hart, A. 2009, MNRAS, 394, 1925 [NASA ADS] [CrossRef] [Google Scholar]
  42. Watson, C. L. 2006, Society for Astronomical Sciences 25th Annual Symposium on Telescope Science, 25, 47 [Google Scholar]
  43. Wittkowski, M., Aufdenberg, J. P, & Kervella, P. 2003, A&A, 413, 711 [Google Scholar]
  44. Wittkowski, M., Boboltz, D. A., Ohnaka, K., Driebe, T., & Scholz, M. 2007, A&A, 470, 191 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Wittkowski, M., Boboltz, D. A., Driebe, T., et al. 2008, A&A, 479, L21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Wittkowski, M., Boboltz, D. A., Ireland, M., et al. 2011, A&A, 532, L7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Wittkowski, M., Hauschildt, P. H., Arroyo-Torres, B., & Marcaide, J. M. 2012, A&A, 540, L12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

All Tables

Table 1

VLTI/AMBER observations.

Table 2

Calibration sources.

Table 3

Information about the variability.

Table 4

Summary of estimated angular diameters.

Table 5

Fundamental parameters of ϵ Oct, β Peg, NU Pav, ψ Peg, and γ Hya.

All Figures

thumbnail Fig. 1

Left: from top to bottom, observed (black) squared visibility amplitudes, UD diameters predicted from our data (blue) and from the best-fit PHOENIX model (red), closure phases in degrees, and normalized flux of β Peg obtained on 2012 Jun. 25. Right: same as left for data obtained on 2012 Aug. 09.

In the text
thumbnail Fig. 5

Average (per baseline) of squared visibility amplitudes taken in the continuum bandpass at 2.152.25 μm for ϵ Oct, β Peg, NU Pav, ψ Peg, and γ Hya (from top to bottom) as a function of spatial frequency. For ϵ Oct and β Peg, the graphics show the data of June and August, and for γ Hya show the data of the both setups. The red lines indicate the best-fit UD models and the blue lines the best-fit PHOENIX models. The dashed lines are the maximum and minimum visibility curves, from which we estimated the errors of the angular diameters.

In the text
thumbnail Fig. 6

Effective temperature versus spectral type of our sources in this paper, RS Cap from Martí-Vidal et al. (2011), and RSGs in Arroyo-Torres et al. (2013), VY CMa from Wittkowski et al. (2012), Betelgeuse from Ohnaka et al. (2011), and VX Sgr from Chen et al. (2007) and Chiavassa et al. (2010). We also show the scales of Dyck et al. (1998), Levesque et al. (2005), and van Belle et al. (2009).

In the text
thumbnail Fig. 7

Location of ϵ Oct, β Peg, NU Pav, ψ Peg, and γ Hya in the HR diagram using our determination of the fundamental parameters. We also show the position of RS Cap by Martí-Vidal et al. (2011). The positions of the stars are compared with evolutionary tracks from Lagarde et al. (2012) for masses of 1.0 M, 1.25 M, 1.5 M, 2.0 M, 2.5 M, 3.0 M, 4.0 M, and 6.0 M.

In the text
thumbnail Fig. 8

Location of ϵ Oct, β Peg, NU Pav, ψ Peg, and γ Hya in the HR diagram using our determination of the fundamental parameters. We also show the position of RS Cap by Martí-Vidal et al. (2011), and the RSGs studied in Arroyo-Torres et al. (2013), VY CMa from Wittkowski et al. (2012), Betelgeuse from Ohnaka et al. (2011), and VX Sgr from Chen et al. (2007) and Chiavassa et al. (2010). The positions of the stars are compared with evolutionary tracks from Ekström et al. (2012) for masses of 3 M, 4 M, 5 M, 7 M, 9 M, 12 M, 15 M, 20 M, 25 M, 32 M, and 40 M. The solid lines are models without rotation, the dashed lines with rotation.

In the text
thumbnail Fig. 2

Left: from top to bottom, observed (black) squared visibility amplitudes, UD diameters predicted from our data (blue) and from the best-fit PHOENIX model (red), closure phases in degrees, and normalized flux of ϵ Oct obtained on 2012 Jun. 25. Right: same as left, but for data obtained on 2012 Aug. 02.

In the text
thumbnail Fig. 3

Left: from top to bottom, observed (black) squared visibility amplitudes, UD diameters predicted from our data (blue) and from the best-fit PHOENIX model (red), closure phases in degrees, and normalized flux of NU Pav obtained on 2012 Aug. 02. Right: same as left, but for data from ψ Peg obtained on 2012 Jun. 16.

In the text
thumbnail Fig. 4

Left: from top to bottom, observed (black) squared visibility amplitudes, UD diameters predicted from our data (blue) and from the best-fit PHOENIX model (red), and closure phases in degrees of γ Hya obtained on 2013 Mar. 16. Right: same as left, but obtained with the MR-K 2.3 μm setting.

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.