Asymmetries on red giant branch surfaces from CHARA/MIRC optical interferometry

TL;DR

This study uses optical interferometry to measure the surface asymmetries of RGB stars, suggesting convection or magnetic activity as potential causes, and highlights the need for further multi-modal observations.

Contribution

First interferometric measurements of surface asymmetries on RGB stars, linking asymmetries to stellar surface activity and emphasizing the importance of combined observational approaches.

Findings

Detected surface asymmetries in all three RGB stars.

Asymmetries tend to be larger in stars with lower surface gravity.

Possible link between asymmetries and convection or magnetic activity.

Abstract

Context. Red giant branch (RGB) stars are very bright objects in galaxies and are often used as standard candles. Interferometry is the ideal tool to characterize the dynamics and morphology of their atmospheres. Aims. We aim at precisely characterising the surface dynamics of a sample of RGB stars. Methods. We obtained interferometric observations for three RGB stars with the MIRC instrument mounted at the CHARA interfer- ometer. We looked for asymmetries on the stellar surfaces using limb-darkening models. Results. We measured the apparent diameters of HD197989 (Epsilon Cyg) = 4.61+-0.02 mas, HD189276 (HR7633) = 2.95+-0.01 mas, and HD161096 (Beta Oph) = 4.43+-0.01 mas. We detected departures from the centrosymmetric case for all three stars with the tendency of a greater effect for lower logg of the sample. We explored the causes of this signal and conclude that a possible explanation to the interferometric signal is the convection-related and/or the magnetic-related surface activity. However, it is necessary to monitor these stars with new observations, possibly coupled with spectroscopy, in order to firmly establish the cause.

Tables

Table 1. Table 1: Log of observations

Date (UT)	Target	$N_{block}$	TTO	Calibrators
			(min)
2016 Jul 6	HD161096 ($\beta$ Oph)	2	20	$ϵ$ Ser, 72 Oph,$\lambda$ Aql
2016 Jul 6	HD189276 (HR 7633)	1	10	$\iota$ Cyg, HR 8185
2016 Jul 6	HD197989 ($ϵ$ Cyg)	1	15	$\gamma$ Lyr
2016 Jul 7	$\beta$ Oph	1	25	$ϵ$ Ser, 72 Oph
2016 Jul 7	HR 7633	1	20	$\iota$ Cyg,
2016 Jul 7	$ϵ$ Cyg	1	30	17 Cyg, $\sigma$ Cyg

Table 2. Table 2: Parameters of the RGB stars.

Stars	Spectral	H mag	$M$	[Fe/H]	$T_{\mathrm{eff}}$	log$g$	R
	typea𝑎aa𝑎a$a$Gray et al. (2003)		$\left[M_{\odot }\right]$		[K]	[cgs]	$\left[R_{\odot }\right]$
$ϵ$ Cyg	K0III-IV	0.200b𝑏bb𝑏b$b$Ducati (2002)	1.84$\pm$0.31d𝑑dd𝑑d$d$Reffert et al. (2015)	-0.11$\pm$0.10d𝑑dd𝑑d$d$Reffert et al. (2015)	4778$\pm$49d𝑑dd𝑑d$d$Reffert et al. (2015)	2.62$\pm$0.10d𝑑dd𝑑d$d$Reffert et al. (2015)	11.08$\pm$0.25d𝑑dd𝑑d$d$Reffert et al. (2015)
HR7633	K4.5IIIa	0.919b𝑏bb𝑏b$b$Ducati (2002)	–	–	4050e𝑒ee𝑒e$e$Lafrasse et al. (2010)	1.70e𝑒ee𝑒e$e$Lafrasse et al. (2010)	–
$\beta$ Oph	K2III	0.354c𝑐cc𝑐c$c$Laney et al. (2012)	1.63$\pm$0.18d𝑑dd𝑑d$d$Reffert et al. (2015)	0.13$\pm$0.10d𝑑dd𝑑d$d$Reffert et al. (2015)	4520$\pm$44d𝑑dd𝑑d$d$Reffert et al. (2015)	2.42$\pm$0.07d𝑑dd𝑑d$d$Reffert et al. (2015)	13.13$\pm$0.32d𝑑dd𝑑d$d$Reffert et al. (2015)

Table 3. Table 3: Apparent diameters of the observed stars at 1.6004 ± plus-or-minus \pm 0.0036 μ 𝜇 \mu m

Star	LD power-law	${\theta }_{LD}$	${\theta }_{LD}$	$\overline{{\chi }^2}$
	exponent	[mas]	$\left[R_{\odot }\right]$
$ϵ$ Cyg	0.25$\pm$0.01	4.61$\pm$0.02	22.09$\pm$0.15	2.45
HR7633	0.14$\pm$0.01	2.95$\pm$0.01	178.11$\pm$11.11	1.32
$\beta$ Oph	0.25$\pm$0.01	4.43$\pm$0.01	23.89$\pm$0.16	1.01

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

Star	B-V	Parallaxc𝑐cc𝑐c$c$van Leeuwen (2007)	SMWd𝑑dd𝑑d$d$Duncan et al. (1991)	$\log R_{\mathrm{HK}}^{\prime }$
		[mas]
$ϵ$ Cyg	1.04a𝑎aa𝑎a$a$Ducati (2002)	44.86$\pm$0.12	0.10	-4.910
HR7633	1.55a𝑎aa𝑎a$a$Ducati (2002)	3.56$\pm$0.21	0.29	-4.132
$\beta$ Oph	1.18b𝑏bb𝑏b$b$Oja (1993)	39.85$\pm$0.17	0.11	-4.738

Full text

*Aims. *We aim at precisely characterising the surface dynamics of a sample of RGB stars.

*Methods. *We obtained interferometric observations for three RGB stars with the MIRC instrument mounted at the CHARA interferometer. We looked for asymmetries on the stellar surfaces using limb-darkening models.

*Results. *We measured the apparent diameters of HD197989 ($\epsilon$ Cyg) = 4.61$\pm$0.02 mas, HD189276 (HR7633) = 2.95$\pm$0.01 mas, and HD161096 ($\beta$ Oph) = 4.43$\pm$0.01 mas. We detected departures from the centrosymmetric case for all three stars with the tendency of a greater effect for lower log$g$ of the sample. We explored the causes of this signal and conclude that a possible explanation to the interferometric signal is the convection-related and/or the magnetic-related surface activity. However, it is necessary to monitor these stars with new observations, possibly coupled with spectroscopy, in order to firmly establish the cause.

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 $\sim$2.0 $M_{\odot}$ (Salaris et al., 2002) with $4\,000\lesssim T_{\rm{eff}}\lesssim 5\,100$ K depending on metallicity, $1.5\lesssim\log{g}\lesssim 3.5$, and $3\lesssim R_{\star}\lesssim 70$ $R_{\odot}$ (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.

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 $\sim$0.5 mas in the H band (ten Brummelaar et al., 2005). The detailed parameters of the RGB stars are reported in Table 2.

The MIRC (Monnier et al., 2004, 2012) is a six-beam combiner operating in the H band (1.5-1.8 $\mu$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 minutes, 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)333This 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/. 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 ($\lesssim 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$\pm$0.0031, 1.7055$\pm$0.0033, 1.6711$\pm$0.0035, 1.6361$\pm$0.0035, 1.6004$\pm$0.0036, 1.5642$\pm$0.0037, 1.5273$\pm$0.0035, and 1.4833$\pm$0.0033 $\mu$m. In the following, we use the wavelength band centred at 1.6004$\pm$0.0036 $\mu$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.

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 $\beta$ Oph, while the other stars are reported in Figs. 3 and 4. The limb-darkening fit is very good (Table 3) with larger residuals for HR7633. We report the first measure of the radius for HR7633, and while the radius of $\epsilon$ Cyg is in good agreement with Reffert et al. (2015), the radius of $\beta$ Oph is slightly smaller.

The observed closure phases display values different from zero or $\pm 180^{\circ}$ for all the observed stars. To determine the amplitude of these deviations, we selected the lowest value between $\mid data-0^{\circ}\mid,\mid data-180^{\circ}\mid,$ and $\mid data+180^{\circ}\mid$ for each data point and then normalised it by the corresponding observed error, $\sigma$. We plot the data departures in Fig. 2, which shows several points diverging from the centrosymmetric case for values higher than 3 $\sigma$: the closure phase departures are smaller for $\beta$ Oph (33 points over 920, 3.6$\%$, higher than 3$\sigma$), intermediate for $\epsilon$ Cyg (70 points over 1056, 6.6$\%$, higher than 3$\sigma$), and larger for HR7633 (88 points over 1384, 6.3$\%$, higher than 3$\sigma$). For $\beta$ Oph and $\epsilon$ Cyg, the spatial frequencies spanned extend to the fourth lobe, while for HR7633, they only reach (partially) the third lobe. The contribution of the small-scale structures increases with the frequency, and that HR7633 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 HR7633, 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 (HR7633), 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 $\log R^{\prime}_{\mathrm{HK}}$ (Table 4) value using the calibration provided by Mittag et al. (2013) because $\log R^{\prime}_{\mathrm{HK}}$ 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 ¡ B-V ¡ 1.18. It should be noted that the B-V colour of HR7633 lies beyond the range of validity of the calibration, which implies that the activity indicator may not be reliable. The values for $\beta$ Oph and $\epsilon$ 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 HR7633 appears to be an active star with $\log R^{\prime}_{\mathrm{HK}}=-4.132$. In addition to this, the stellar diameter of 178.11 R⊙ of HR7633 (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 $\epsilon$ Cyg is a known double system (Poveda et al., 1994) where the primary star ($\epsilon$ Cyg) is $\approx 25000$ brighter than its companion (dwarf M4 star) with $m_{V}$=2.45 ($m_{H}$=0.20, Ducati, 2002) and $m_{V}$=11.98 ($m_{H}$=9.58555estimate based on PHOENIX models (e.g. Allard et al., 1997, and perso.ens-lyon.fr/france.allard/)), 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 $\beta$ Oph and HR7633 (McDonald et al., 2011, and McDonald private communication). The same is expected for $\epsilon$ 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 HR7633. We measured the apparent diameters of HD197989 ($\epsilon$ Cyg) = 4.61$\pm$0.02 mas, HD189276 (HR7633) = 2.95$\pm$0.01 mas, and HD161096 ($\beta$ Oph) = 4.43$\pm$0.01 mas.

Moreover, the closure phases denote departure points from the centrosymmetric case (closure phases not equal to 0 or $\pm 180^{\circ}$) with values greater than 3$\sigma$. The departures seem to be qualitatively correlated with $\log{g}$ of the observed stars. HR7633, 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.

Acknowledgements.

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

AC, MM, PK, GP acknowledge financial support from ”Programme National de Physique Stellaire” (PNPS) of CNRS/INSU, France.

RN acknowledges the Fizeau exchange visitors program in optical interferometry - WP14 OPTICON/FP7 (2013-2016, grant number 312430).

RC acknowledges the funding provided by The Danish National Research Foundation (Grant DNRF106).

Appendix A Additional figures for $\epsilon$ Cyg and HR7633