Free Access
Issue
A&A
Volume 600, April 2017
Article Number L2
Number of page(s) 5
Section Letters
DOI https://doi.org/10.1051/0004-6361/201730438
Published online 20 March 2017

© ESO, 2017

1. Introduction

Red giant branch (RGB) stars have evolved from the main sequence and are powered by hydrogen burning in a thin shell surrounding their helium core. This evolutionary phase precedes the asymptotic giant branch (AGB) and is characterised during the evolved states by more expanded and deformed outer layers. RGB stars are bright candles in galaxies, and the accurate determination of their fundamental parameters and chemical composition is essential for tracing the morphology and the evolution of the Galaxy, for probing distant stellar populations, and for characterising globular clusters. Their masses are typically lower than ~2.0 M (Salaris et al. 2002) with 4000 ≲ Teff ≲ 5100 K depending on metallicity, 1.5 ≲ log g ≲ 3.5, and 3 ≲ R ≲ 70 R (van Belle et al. 1999; Baines et al. 2010).

Like all late-type stars, red giant atmospheres are made complex by convective motions and turbulent flows. Convection contributes significantly to the transportation of energy from the stellar interior to the outer layer, and in the photosphere, it manifests itself as a granulation pattern characterised by dark intergranular lanes of downflowing cooler plasma and bright areas (the granules) where hot plasma rises (see the review of Nordlund et al. 2009). The granules cause an inhomogeneous stellar surface that changes with time. The granulation potentially acts as an intrinsic noise in stellar parameters, radial velocity, and chemical abundance determinations. In addition to this, a magnetic field may be present, as detected in several RGB stars by Aurière et al. (2015), and this may contribute to a bias in the stellar parameter determination. The characterisation of the dynamics and morphology of RGB stars is important to quantify the effect of the granulation and magnetic fields, and thanks to its high angular resolution, interferometry is the ideal tool for this purpose.

In this Letter, we present the detection of an interferometric signal at high spatial frequencies for three RGB stars using the MIRC instrument mounted at the CHARA interferometer. We analyse the possible causes of this signal.

Table 1

Log of observations.

Table 2

Parameters of the RGB stars.

thumbnail Fig. 1

Top and central panels: squared visibilities with the limb-darkening fit residuals (top and central panels) at the wavelength band 1.6004 ± 0.0036 μm. The continuous line in the top panel corresponds to the limb-darkening fit whose parameters are reported in Table 3. The black horizontal dashed lines in the central panel correspond to the value of 1σ, and the red line shows the value of 3σ. Bottom panel: closure phase data points of β Oph (Table 2) for all the wavelengths. The red data correspond to closure phase departures larger than 3σ (see Fig. 2). The horizontal dashed lines in the bottom panel display the zero or ± 180° values.

2. Interferometric observations with MIRC at CHARA

We collected observations of three stars (Table 1) using the Michigan Infrared Combiner (MIRC) on the Georgia State University Center for High Angular Resolution Astronomy (CHARA). The CHARA array is located on Mount Wilson, CA, and consists of six 1 m telescopes for a total of 15 baselines ranging in length from 34 m to 331 m, resulting in an angular resolution of ~0.5 mas in the H band (ten Brummelaar et al. 2005). The detailed parameters of the RGB stars are reported in Table 2.

thumbnail Fig. 2

Closure phase departures from zero or ±180° for all the observed stars and for all the wavelengths. For every data point, the lowest value between |data-0°|, |data-180°| , and | data+180°| is computed and normalised by the corresponding observed error, σ. β Oph is shown in the left panel, ϵ Cyg in the central, and HR 7633 in the right panel. The horizontal black line corresponds to the value of 1σ and the red line to 3σ.

The MIRC (Monnier et al. 2004, 2012) is a six-beam combiner operating in the H band (1.51.8 μm) at low spectral resolution (R = 30). Each observing block consists of observations of a calibrator, a target, and when possible, a second calibrator. Time spent collecting data on the target ranged within 10 to 30 min, excluding background and other calibration frames, depending on observing conditions. We used the latest version of the MIRC reduction pipeline (as of October 2016) and as previously described in Monnier et al. (2007) and Monnier et al. (2012). The pipeline uses Fourier transforms to compute squared visibilities, which are then averaged and corrected for biases. We determined the bispectrum using the phases and amplitudes of three baselines in a closed triangle. We calibrated photometric amplitudes using a beam splitter following spatial filtering (Che et al. 2010)1. Because we do not expect significant brightness variation over short time periods for these targets, we combined the two nights of observations for each star into single files, accounting for systematic error by applying multiplicative and additive errors as described in Monnier et al. (2012). At low visibilities (10-3), the signal-to-noise ratio of the data decreases because of cross-talk. We therefore remain cautious when interpreting data at such low visibilities.

The observations were collected in eight different wavelength bands: 1.7379 ± 0.0031, 1.7055 ± 0.0033, 1.6711 ± 0.0035, 1.6361 ± 0.0035, 1.6004 ± 0.0036, 1.5642 ± 0.0037, 1.5273 ± 0.0035, and 1.4833 ± 0.0033 μm. In the following, we use the wavelength band centred at 1.6004 ± 0.0036 μm for the apparent diameter determination with visibility curves because it corresponds to the H continuous opacity minimum (consequently closest to the continuum forming region). For all the closure phase plots, we used the full set of wavelength bands.

Table 3

Apparent diameters of the observed stars at 1.6004 ± 0.0036 μm.

3. Discussion

We fitted the data, based only on the squared visibilities, with an power-law limb-darkened disk model whose parametric values are reported in Table 3. Figure 1 displays the example of β Oph, while the other stars are reported in Figs. A.1 and A.2. The limb-darkening fit is very good (Table 3) with larger residuals for HR 7633. We report the first measure of the radius for HR 7633, and while the radius of ϵ Cyg is in good agreement with Reffert et al. (2015), the radius of β Oph is slightly smaller.

The observed closure phases display values different from zero or ± 180° for all the observed stars. To determine the amplitude of these deviations, we selected the lowest value between |data |, |data 180°| , and |data+180°| for each data point and then normalised it by the corresponding observed error, σ. We plot the data departures in Fig. 2, which shows several points diverging from the centrosymmetric case for values higher than 3σ: the closure phase departures are smaller for β Oph (33 points over 920, 3.6%, higher than 3σ), intermediate for ϵ Cyg (70 points over 1056, 6.6%, higher than 3σ), and larger for HR 7633 (88 points over 1384, 6.3%, higher than 3σ). For β Oph and ϵ Cyg, the spatial frequencies spanned extend to the fourth lobe, while for HR 7633, they only reach (partially) the third lobe. The contribution of the small-scale structures increases with the frequency, and that HR 7633 displays departures already on the second lobe indicates that this star is most likely the most asymmetric of the three.

Moreover, the closure phase departures of the three stars seem to be correlated with log g of the stars (Table 2), the latter is an indicator of the evolutionary phase: the largest deviations are for HR 7633, which has log g = 1.70. This denotes that the size of the granules become more significant with respect to the disk size as the surface gravity decreases, and therefore brightness fluctuations become more important. This idea is supported by previous work showing even larger departures from centrosymmetry for very evolved stars such as AGBs (Wittkowski et al. 2016; Chiavassa et al. 2010c; Ragland et al. 2006) and red supergiant stars (Montargès et al. 2016; Chiavassa et al. 2010b). However, we also note that we detected the largest deviations for the faintest star (HR 7633), and this may indicate that we underestimated the errors.

We now present a tentative explanation of these closure phase departures. A more complete analysis will be presented in a forthcoming paper.

Stellar surface asymmetries in the brightness distribution can be either due to convection-related and/or activity-related structures, to a companion, or to a clumpy dust envelope around the stars. In the following, we analyse the different possibilities.

A first hypothesis concerns convection-related surface structures affecting the interferometric observables. The expected convection turnover time in such a star is between hours to days or weeks, depending on the stellar fundamental parameters. Chiavassa et al. (2010a, 2014) showed that stellar granulation manifests itself as surface asymmetries in the brightness distribution, and more precisely, in the closure phase signal.

A second hypothesis is the stellar magnetic activity. Chiavassa et al. (2014) and Ligi et al. (2015) showed that star spots caused by the stellar magnetic field affect the closure phase signal in a similar way as the granulation. To determine its impact, we estimated the indicator SMW based on the historic stellar activity data of chromospheric line emission. It measures the strength of the chromospheric emission core of the H and K lines of Ca II (Middelkoop 1982). For the three stars, we converted SMW into (Table 4) value using the calibration provided by Mittag et al. (2013) because is independent of the stellar type. In particular, we adopted their Eq. (13) and calibration for giant stars (their Eq. (23)), which is valid for stars with 0.76 <BV< 1.18. It should be noted that the BV colour of HR 7633 lies beyond the range of validity of the calibration, which implies that the activity indicator may not be reliable. The values for β Oph and ϵ Cyg fall into or very close to the stellar inactivity region (Figs. 3 and 4 in Wright 2004). However, the star spots may still exist even if the activity is low, and their signature can only be distinguished by performing more interferometric observations coupled with spectroscopic observations.

The source HR 7633 appears to be an active star with . In addition to this, the stellar diameter of 178.11 R of HR 7633 (Table 3) as well as its low surface gravity (Table 2) seem to indicate that this star approaches the red giant tip towards the AGB phase. This evolutionary step is characterised by prominent stellar granulation accompanied by non-negligible dynamics (Ludwig & Kučinskas 2012).

Concerning any possible companions, only ϵ Cyg is a known double system (Poveda et al. 1994) where the primary star (ϵ Cyg) is 25 000 brighter than its companion (dwarf M4 star) with mV = 2.45 (mH = 0.20, Ducati 2002) and mV = 11.98 (mH = 9.582), respectively. In addition to this, the two stars lie 78.1′′ from each other. For these reasons, we estimate that the secondary star has a negligible effect on the observed data.

No dust production is detected around RGBs for β Oph and HR 7633 (McDonald et al. 2011; and McDonald, priv. comm.). The same is expected for ϵ Cyg. Dust is a sign of strong mass-loss; in evolved AGB stars, for example, the mass loss can be 1000 greater than in RGB stars.

Stellar convection-related and/or the magnetic-related surface activity are the most plausible explanation. To firmly establish the cause, it is necessary to monitor these stars with new interferometric observations, possibly coupled with spectropolarimetric data in the visible in order to measure their magnetic activity strength.

4. Conclusions

We presented observations of three RGB stars using the MIRC instrument at the CHARA interferometer. We showed that for all stars, the limb-darkening fit is very good with larger residuals for HR 7633. We measured the apparent diameters of HD 197989 (ϵ Cyg) = 4.61 ± 0.02 mas, HD 189276 (HR 7633) = 2.95 ± 0.01 mas, and HD 161096 (β Oph) = 4.43 ± 0.01 mas.

Moreover, the closure phases denote departure points from the centrosymmetric case (closure phases not equal to 0 or ± 180°) with values greater than 3σ. The departures seem to be qualitatively correlated with log g of the observed stars. HR 7633, with the lowest log g of the sample, shows the highest deviations: the more the star evolves, the more significant the size of the granules becomes with respect to the disk size.

We explored the possible causes of the break in symmetry in the brightness distribution and found that a possible explanation of this interferometric signal is the granulation and/or the stellar magnetic activity at its surface. However, it is not possible to confirm this at the moment, and a more complete analysis will be presented in a forthcoming paper.

Table 4

Photometric colours and stellar activity indicator of the observed RGB stars.


1

This research has made use of the Jean-Marie Mariotti Center SearchCal service Bonneau et al. (2006), available at http://www.jmmc.fr/searchcal, co-developped by FIZEAU and LAOG/IPAG, and of CDS Astronomical Databases SIMBAD and VIZIER, available at http://cdsweb.u-strasbg.fr/

2

Estimate based on PHOENIX models (e.g. Allard et al. 1997,and perso.ens-lyon.fr/france.allard/)

Acknowledgments

This work is based upon observations obtained with the Georgia State University Center for High Angular Resolution Astronomy Array at Mount Wilson Observatory. The CHARA Array is supported by the National Science Foundation under Grants Nos. AST-1211929 and AST-1411654. Institutional support has been provided from the GSU College of Arts and Sciences and the GSU Office of the Vice President for Research and Economic Development. A.C., M.M., P.K., G.P. acknowledge financial support from “Programme National de Physique Stellaire” (PNPS) of CNRS/INSU, France. R.N. acknowledges the Fizeau exchange visitors program in optical interferometry – WP14 OPTICON/FP7 (20132016, grant number 312430). R.C. acknowledges the funding provided by The Danish National Research Foundation (Grant DNRF106).

References

  1. Allard, F., Hauschildt, P. H., Alexander, D. R., & Starrfield, S. 1997, ARA&A, 35, 137 [NASA ADS] [CrossRef] [Google Scholar]
  2. Aurière, M., Konstantinova-Antova, R., Charbonnel, C., et al. 2015, A&A, 574, A90 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Baines, E. K., Döllinger, M. P., Cusano, F., et al. 2010, ApJ, 710, 1365 [NASA ADS] [CrossRef] [Google Scholar]
  4. Bonneau, D., Clausse, J.-M., Delfosse, X., et al. 2006, A&A, 456, 789 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Boyajian, T. S., McAlister, H. A., van Belle, G., et al. 2012, ApJ, 746, 101 [NASA ADS] [CrossRef] [Google Scholar]
  6. Che, X., Monnier, J. D., & Webster, S. 2010, in Optical and Infrared Interferometry II, Proc. SPIE, 7734, 77342V [NASA ADS] [CrossRef] [Google Scholar]
  7. Chiavassa, A., Collet, R., Casagrande, L., & Asplund, M. 2010a, A&A, 524, A93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Chiavassa, A., Haubois, X., Young, J. S., et al. 2010b, A&A, 515, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Chiavassa, A., Lacour, S., Millour, F., et al. 2010c, A&A, 511, A51 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Chiavassa, A., Ligi, R., Magic, Z., et al. 2014, A&A, 567, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Ducati, J. R. 2002, VizieR Online Data Catalog: 2/237 [Google Scholar]
  12. Duncan, D. K., Vaughan, A. H., Wilson, O. C., et al. 1991, ApJS, 76, 383 [NASA ADS] [CrossRef] [Google Scholar]
  13. Gray, R. O., Corbally, C. J., Garrison, R. F., McFadden, M. T., & Robinson, P. E. 2003, AJ, 126, 2048 [NASA ADS] [CrossRef] [Google Scholar]
  14. Lafrasse, S., Mella, G., Bonneau, D., et al. 2010, in Optical and Infrared Interferometry II, Proc. SPIE, 7734, 77344E–77344E–11 [CrossRef] [Google Scholar]
  15. Laney, C. D., Joner, M. D., & Pietrzyński, G. 2012, MNRAS, 419, 1637 [NASA ADS] [CrossRef] [Google Scholar]
  16. Ligi, R., Mourard, D., Lagrange, A.-M., Perraut, K., & Chiavassa, A. 2015, A&A, 574, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Ludwig, H.-G., & Kučinskas, A. 2012, A&A, 547, A118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. McDonald, I., Boyer, M. L., van Loon, J. T., et al. 2011, ApJS, 193, 23 [NASA ADS] [CrossRef] [Google Scholar]
  19. Middelkoop, F. 1982, A&A, 107, 31 [NASA ADS] [Google Scholar]
  20. Mittag, M., Schmitt, J. H. M. M., & Schröder, K.-P. 2013, A&A, 549, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Monnier, J. D., Berger, J.-P., Millan-Gabet, R., & tenBrummelaar, T. A. 2004, in New Frontiers in Stellar Interferometry, ed. W. A. Traub, Proc. SPIE, 5491, 1370 [Google Scholar]
  22. Monnier, J. D., Zhao, M., Pedretti, E., et al. 2007, Science, 317, 342 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  23. Monnier, J. D., Che, X., Zhao, M., et al. 2012, ApJ, 761, L3 [NASA ADS] [CrossRef] [Google Scholar]
  24. Montargès, M., Kervella, P., Perrin, G., et al. 2016, A&A, 588, A130 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Nordlund, Å., Stein, R. F., & Asplund, M. 2009, Liv. Rev. Sol. Phys., 6, 2 [Google Scholar]
  26. Oja, T. 1993, A&AS, 100, 591 [NASA ADS] [Google Scholar]
  27. Poveda, A., Herrera, M. A., Allen, C., Cordero, G., & Lavalley, C. 1994, Rev. Mex. Astron. Astrofis., 28, 43 [NASA ADS] [Google Scholar]
  28. Ragland, S., Traub, W. A., Berger, J.-P., et al. 2006, ApJ, 652, 650 [NASA ADS] [CrossRef] [Google Scholar]
  29. Reffert, S., Bergmann, C., Quirrenbach, A., Trifonov, T., & Künstler, A. 2015, A&A, 574, A116 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Salaris, M., Cassisi, S., & Weiss, A. 2002, PASP, 114, 375 [NASA ADS] [CrossRef] [Google Scholar]
  31. ten Brummelaar, T. A., McAlister, H. A., Ridgway, S. T., et al. 2005, ApJ, 628, 453 [NASA ADS] [CrossRef] [Google Scholar]
  32. van Belle, G. T., Lane, B. F., Thompson, R. R., et al. 1999, AJ, 117, 521 [NASA ADS] [CrossRef] [Google Scholar]
  33. van Leeuwen, F. 2007, A&A, 474, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Wittkowski, M., Chiavassa, A., Freytag, B., et al. 2016, A&A, 587, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Wright, J. T. 2004, AJ, 128, 1273 [NASA ADS] [CrossRef] [Google Scholar]
  36. Zhao, M., Gies, D., Monnier, J. D., et al. 2008, ApJ, 684, L95 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Additional figures for ϵ Cyg and HR 7633

thumbnail Fig. A.1

Same as in Fig. 1, but for the star ϵ Cyg.

thumbnail Fig. A.2

Same as in Fig. 1, but for the star HR 7633.

All Tables

Table 1

Log of observations.

Table 2

Parameters of the RGB stars.

Table 3

Apparent diameters of the observed stars at 1.6004 ± 0.0036 μm.

Table 4

Photometric colours and stellar activity indicator of the observed RGB stars.

All Figures

thumbnail Fig. 1

Top and central panels: squared visibilities with the limb-darkening fit residuals (top and central panels) at the wavelength band 1.6004 ± 0.0036 μm. The continuous line in the top panel corresponds to the limb-darkening fit whose parameters are reported in Table 3. The black horizontal dashed lines in the central panel correspond to the value of 1σ, and the red line shows the value of 3σ. Bottom panel: closure phase data points of β Oph (Table 2) for all the wavelengths. The red data correspond to closure phase departures larger than 3σ (see Fig. 2). The horizontal dashed lines in the bottom panel display the zero or ± 180° values.

In the text
thumbnail Fig. 2

Closure phase departures from zero or ±180° for all the observed stars and for all the wavelengths. For every data point, the lowest value between |data-0°|, |data-180°| , and | data+180°| is computed and normalised by the corresponding observed error, σ. β Oph is shown in the left panel, ϵ Cyg in the central, and HR 7633 in the right panel. The horizontal black line corresponds to the value of 1σ and the red line to 3σ.

In the text
thumbnail Fig. A.1

Same as in Fig. 1, but for the star ϵ Cyg.

In the text
thumbnail Fig. A.2

Same as in Fig. 1, but for the star HR 7633.

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.