This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Study of a red clump giant, KIC 11087027, with high rotation and strong infrared excess - Evidence of tidal interaction for high lithium abundance

Raghubar Singh National Astronomical Observatories, Chinese Academy of Sciences, 20A Datun Road, Chaoyang District, Beijing, China (100101) Anohita Mallick Indian Institute of Astrophysics, 560034, 100ft road Koramangala, Bangalore, India Pondicherry University, R. V. Nagara, Kala Pet, 605014, Puducherry, India Bacham E. Reddy Indian Institute of Astrophysics, 560034, 100ft road Koramangala, Bangalore, India Jeewan C. Pandey Aryabhatta Research Institute of Observational Sciences (ARIES), Manora Peak, Nainital 263001, India Gang Zhao National Astronomical Observatories, Chinese Academy of Sciences, 20A Datun Road, Chaoyang District, Beijing, China (100101)
Abstract

This paper presents results from Kepler photometric light curves and high-resolution spectroscopic study of a super Li-rich giant KIC11087027. Using the light curve analysis, we measured the star’s rotational period Prot=30.4±\pm0.1 days, which translates to rotational velocity Vrot=19.5 ±\pm 1.7 km s-1. Star’s location in the HR-diagram, derived values of C12/13C{}^{12}C/^{13}C = 7±\pm1 and [C/N]=0.95±0.2[C/N]=-0.95\pm 0.2, and the inferred asteroseismic parameters from secondary calibration based on spectra suggest star is a low-mass red clump giant in the He-core burning phase. Using Gaia data, we found evidence of variation in radial velocity and proper motion, indicative of presence of an unresolved binary. The large Vrot is probably a result of tidal synchronization combined with the after-effects of He-flash, in which the size of the star is reduced significantly. The simultaneous presence of features like high rotation, very high Li abundance, strong dust shell, and strong flares in a single star is relatively uncommon, suggesting that the star experiencing tidal synchronization has recently undergone He-flash. The results pose a question whether the binary interaction, hence the high rotation, is a prerequisite for dredging-up of the high amounts of Li from the interior to the photosphere during or immediately after the He-flash event.

low mass stars — rotation — stellar flares — Li-rich giants—mass loss—stellar activity
journal: ApJ

1 Introduction

A small group of red giants show very high Li abundances contrary to the general understanding that the Li gets destroyed in stars. This has been an anomaly ever since the discovery of the Li-rich giant by Wallerstein & Sneden (1982). Recent studies demonstrated that high lithium abundance among red clump (RC) giants is common (Kumar et al., 2020) and showed that the He-flash event, immediately preceding the RC, holds the key for Li enhancement (Kumar et al., 2011; Casey et al., 2019; Deepak & Reddy, 2019; Singh et al., 2019; Yan et al., 2021). The average Li abundance of RC giants, post He-flash, is about a factor of 40 more than their counterparts on the upper red giant branch (RGB), closer to the RGB tip (Kumar et al., 2020). In addition, the study by Singh et al. (2021) showed that A(Li) evolves with the g-mode period spacing (ΔΠ1\Delta\Pi_{1}) evolution. During transition phase from the degenerate core to the convective core He burning, the (ΔΠ1\Delta\Pi_{1}) value increases. As a result, very high A(Li) values are mostly seen at relatively lower (ΔΠ1\Delta\Pi_{1}) values, i.e., immediately after the He-flash. The normal Li RC giants are seen at high ΔΠ1\Delta\Pi_{1} values. From this, they hypothesised that the super Li-rich (SLR) giant phase (A(Li)111A(Li)=12+log(N(Li)(N(H))\rm log(\frac{N(Li)}{(N(H)})>>3.2 dex) is a short-lived phenomenon, and giants with very high Li abundance have undergone Li enhancement very recently, suggesting the He-flash as the source of Li enhancement. Further, the study by Mallick et al. (2023) showed occurrence of Li-rich giants only among low-mass (M\leq 2M) giants and none among intermediate (M>>2M) mass RC giants providing an indirect evidence that Li enhancement occurs during the He-flash as the He-flash expected to happen in only low-mass giants (Miller Bertolami et al., 2020). We also note that there are Li abundance studies among cluster giants. Some of the super Li-rich giants in the clusters found to be on the RGB (Sanna et al., 2020; Nagarajan et al., 2023; Tsantaki et al., 2023). If these giants are indeed ascending the RGB first time then there may be multiple ways for Li-enrichment in giants.

However, it is not understood what physical mechanism is responsible for the production and transportation of Li to the photosphere. Recent attempts to explain the physical process include internal gravity waves generated during He-flash (e.g. Schwab, 2020; Jermyn & Fuller, 2022) and thermohaline mixing (Gao et al., 2022). There are proposals in literature to drive the convection and mixing of A(Li) with the upper atmospheres, like the differential rotation due to faster core rotation (Fekel & Balachandran, 1993; Simon & Drake, 1989). Stars may get faster rotations due to tidal interaction in binary systems driving the desired mixing process(Denissenkov & Herwig, 2004; Casey et al., 2019). Any external events like mergers or tidal interactions with companions may spin-up the star, generating a strong magnetic field and creating circumstellar envelopes. One would expect all or most of these features if the Li enhancement occurred recently due to external events like mergers.

In this paper, we analysed high-resolution spectra of KIC11087027 for the first time. The star is a super Li-rich giant with features such as high rotation, infrared excess (IR-excess) and high chromospheric activity, indicating some recent binary interaction.

2 Observations and Data Reduction

KIC11087027 was found in the sample of 12500 stars common in LAMOST (Zhao et al., 2006; Cui et al., 2012; Zhao et al., 2012) and Kepler (Borucki et al., 2010) field studied by Singh et al. (2019). This star was excluded in Singh et al. (2019) because of the absence of oscillation modes in Kepler Power Density Spectra (PDS). The star has peculiar features like the presence of rotational modulation, flares, and IR-excess. To supplement the photometric data, we obtained high-resolution (R=60,000) optical spectra using the 2.0-m Hanle Chandra Telescope (HCT) equipped with Hanle Echelle Spectrograph (HESP)222https://www.iiap.res.in/hesp/ for abundance analysis. We took three frames, each with a 40-minute exposure. The spectra cover the wavelength range from 3700 to 9300 Å. We obtained calibration image frames - bias, flat, and Th-Ar lamp for wavelength calibration. We also obtained high-rotation hot star spectra for removing telluric lines. We followed the standard data reduction procedure for Echelle spectral reduction using IRAF 333https://iraf-community.github.io/. The resultant spectral image has a signal-to-noise ratio (SNR) of 86. The extracted spectrum was wavelength calibrated and continuum fitted for further analysis and derivation of stellar atmospheric parameters and elemental abundances.

Refer to caption
Figure 1: (a): Location of KIC11087027 (solid square), in the Hertzsprung-Russell diagram along with MESA-MIST (Paxton et al., 2011; Dotter, 2016) evolutionary tracks for [Fe/H]=-0.55 dex with different masses. Panels (b),(c), and (d) are the spectral synthesis of Li resonance lines at 6707 Å  Li subordinate line at 6103 Å  and CN red line band near 8004 Å, respectively.

3 Analysis and results

3.1 Stellar parameters and Abundances

Stellar parameters (Teff, log gg, [Fe/H], Vsini, ξt\xi_{\rm t}), see Table  1, are derived from high-resolution spectra. We used FASMA444https://github.com/MariaTsantaki/fasma-synthesis code (Tsantaki et al., 2018) for spectrum synthesis to derive stellar parameters. Star’s luminosity is derived using the Gaia Gmag (Gaia Collaboration et al., 2016, 2018) with appropriate bolometric correction (Andrae et al., 2018) and the parallax given in Gaia. Extinction coefficient in G band, AG, is calculated by dustapproxdustapprox code (Fouesneau et al., 2022) 555https://mfouesneau.github.io/dustapprox. Observed Heliocentric radial velocity is derived by cross-correlating the continuum-fitted spectra with the template spectra of Arcturus.

Elemental abundances and carbon isotopic ratio are derived by matching the predicted spectra with the observed ones. Model spectra are produced using the local thermodynamic equilibrium (LTE) stellar atmospheric models (Castelli & Kurucz, 2004) and the 1-D radiative transfer code MOOG (Sneden, 1973). The atomic and molecular line list is taken from Linemake666https://github.com/vmplacco/linemake code (Placco et al., 2021). LTE Li abundance measured from the Li resonance line at 6707.78 Å  is A(Li)=4.35±\pm0.10 dex and from the subordinate line at 6103 Å  is A(Li)=3.85±\pm0.1 dex. The Li abundance is corrected for NLTE effects by following correction values provided by Lind et al. (2009) (Table 1). Carbon abundance is derived from atomic Carbon lines at 5086 and 5384.3 Å  whereas the Nitrogen abundance is derived from C1412N\rm{}^{12}C^{14}N molecular line at 6486.4 Å. The carbon isotopic ratio is measured from 13C14N molecular line in the red band of the spectra at 8004.72 Å. Errors in Li abundance and carbon isotopic ratio are derived as the quadratic sum of uncertainties in stellar parameters (Teff, log gg, [Fe/H], ξt\xi_{\rm t}) and SNR. Abundances of other elements are derived using equivalent widths.

Table 1: Stellar parameters of KIC11087027. Elemental abundance w.r.t. to the Sun from Asplund et al. (2009).
Name KIC11087027 reference
R.A. 19:33:14.212
DEC. +48:41:50.83
V (mag) 11.23±0.0911.23\pm 0.09 APASS
AG (mag) 0.110.11 this work
Parallax (mas) 0.86±0.050.86\pm 0.05 Gaia DR3
RUWE 4.384.38 Gaia DR3
Prot (Days) 30.38±0.1230.38\pm 0.12 this work
Vrot(km s-1) 19.5 this work
Vsini(km s-1) 10.1±0.210.1\pm 0.2 this work
Inclination angle 35 this work
Number of flares 15 this wok
Teff [K] 4920±604920\pm 60 this work
log gg 2.64±0.052.64\pm 0.05 this work
[Fe/H]\rm[Fe/H] 0.54±0.05-0.54\pm 0.05 this work
ξt\xi_{t}(km s-1) 0.95±0.050.95\pm 0.05 this work
log(L/L)\log(L/L_{\odot}) 1.77±0.051.77\pm 0.05 this work
Mass (M) 1.014±0.21.014\pm 0.2 this work
Radius (R) 11.69±111.69\pm 1 this work
RVGaia (km s-1) 26.70±0.3826.70\pm 0.38 Gaia DR3
RVLAMOST (km s-1) 22.71±4.1122.71\pm 4.11 Luo+2018
RVHESP,Helio (km s-1) 19.98±0.2619.98\pm 0.26 this work
A(Li)NLTE(6103 Å) 3.95±0.13.95\pm 0.1 this work
A(Li)NLTE(6707 Å) 3.85±0.13.85\pm 0.1 this work
[CI/Fe]\rm[CI/Fe] 0.18±0.1-0.18\pm 0.1 this work
[NI/Fe]\rm[NI/Fe] 0.77±0.20.77\pm 0.2 this work
[NaI/Fe]\rm[NaI/Fe] 0.01±0.150.01\pm 0.15 this work
[MgI/Fe]\rm[MgI/Fe] 0.03±0.050.03\pm 0.05 this work
[AlI/Fe]\rm[AlI/Fe] 0.01±0.250.01\pm 0.25 this work
[SiI/Fe]\rm[SiI/Fe] 0.02±0.3-0.02\pm 0.3 this work
[CaI/Fe]\rm[CaI/Fe] 0.03±0.3-0.03\pm 0.3 this work
[ScII/Fe]\rm[ScII/Fe] 0.09±0.1-0.09\pm 0.1 this work
[TiII/Fe]\rm[TiII/Fe] 0.0±0.10.0\pm 0.1 this work
[VI/Fe]\rm[VI/Fe] 0.13±0.220.13\pm 0.22 this work
[CrI/Fe]\rm[CrI/Fe] 0.02±0.3-0.02\pm 0.3 this work
[CoI/Fe]\rm[CoI/Fe] 0.28±0.18-0.28\pm 0.18 this work
[NiI/Fe]\rm[NiI/Fe] 0.34±0.21-0.34\pm 0.21 this work
[ZnI/Fe]\rm[ZnI/Fe] 0.23±0.20.23\pm 0.2 this work
[SrI/Fe]\rm[SrI/Fe] 0.51±0.2-0.51\pm 0.2 this work
[BaII/Fe]\rm[BaII/Fe] 0.04±0.060.04\pm 0.06 this work
[NdII/Fe]\rm[NdII/Fe] 0.22±0.10.22\pm 0.1 this work
[EuII/Fe]\rm[EuII/Fe] 0.15±0.080.15\pm 0.08 this work
C12/13C\rm{}^{12}C/^{13}C 7±17\pm 1 this work

3.2 Evolutionary Phase

The evolutionary phase of a star is determined using the star’s location in the HR diagram combined with evolutionary tracks, and using abundances ratios [C/N] and C12/13C\rm{}^{12}C/^{13}C. As shown in Fig 1, the star’s luminosity and Teff places it in a region in the HR diagram that overlaps with the red clump region. Similarly, observed [C/N]= -0.95±\pm0.22 and [Fe/H]= -0.55±\pm0.05, locates KIC11087027 in the red clump region in [C/N] vs [Fe/H] diagram (see, Fig 5 top panel in Hawkins et al., 2018). The very low value of C12/13C=7±1\rm{}^{12}C/^{13}C=7\pm 1 suggests that the star is highly evolved and, at least, beyond the luminosity bump. Though we have Kepler photometric data for 18 quarters, the power density spectrum (PDS) of KIC11087027 shown in Figure 2(b) do not show Gaussian power excess with oscillation modes. Inhibition of oscillation modes is caused by enhanced magnetic activity (Chaplin et al., 2011; Gaulme et al., 2014). As a result, we resorted to adopt asteroseismic parameters using secondary calibrations by (Wang et al., 2023) who have derived ΔΠ1=187\rm\Delta\Pi_{1}=187 s and Δν=4.4μHz\rm\Delta\nu=4.4~{}\mu Hz using spectral indices from LAMOST spectra. The values suggest KIC11087027 is a red clump giant in He-core burning phase (Bedding et al., 2011).

3.3 Detection of stellar flare and measuring flare energy

Stellar flares are manifestations of magnetic energy. Flares are common features among main sequence stars but uncommon among red giants due to reduced magnetic field. Sometimes, stellar superflares may arise in binary companions or due to recent activity with stellar companion-like planets (Cuntz et al., 2000). We have investigated long cadence (29.4 minutes) white light photometric data from the Kepler mission obtained over four years. The visual inspection of the light curve showed many fast rises and exponential decay-like profiles in amplitude, suggesting the presence of flares. We used FLATW,RM\rm FLATW^{,}RM code (Vida & Roettenbacher, 2018) to detect stellar flares and measure flare parameters including equivalent duration and amplitude. We discovered 15 white light flares (WLFs) in KIC11087027 of the duration of a few hours each; see Figure 2 (b). WLF energy is derived from the equivalent duration (ED) of flare and quiescent luminosity. Quiescent luminosity in the Kepler wavelength band is measured by integrating the convolution of bolometric flux (FλF_{\lambda}) and Kepler response function (SKp\rm S_{Kp}).

Eflare=4πR2×EDλ1λ2FλSKp𝑑λE_{flare}=4\pi R^{2}\times ED\int_{\lambda_{1}}^{\lambda_{2}}F_{\lambda}*S_{Kp}d\lambda

Here λ1\lambda_{1} and λ2\lambda_{2} are the lower and upper wavelength limits of the Kepler filter. Bolometric flux FλF_{\lambda} for stellar parameters is taken from VOSA 777http://svo2.cab.inta-csic.es and Kepler response function SKp\rm S_{Kp} from Kepler website 888https://keplergo.github.io/KeplerScienceWebsite/. Energy of flares is in the range of 1035103710^{35}-10^{37} erg, which makes KIC11087027 as super flaring star.

Refer to caption
Figure 2: Kepler light curve analysis of KIC11087027: a) black is raw and red is a median smoothed light curve. b) comparison of the power spectral density of KIC11087027 (black) with a similar mass star KIC12645107 (red) with no activity. c) 7th quarter Kepler LC of KIC11087027 and three flares are marked. (d) The detrended light curve in the flare region and fitted flare model. Measured flare parameters are given in the figure. Panel (e): The Lomb-Scargle periodogram of the entire 18-quarter light curve. The vertical red dashed line is measured period. (f): The rotation period measured in different quarters (Q2—Q16) gives different values.

3.4 IR-excess and chromospheric activity

We looked for infrared excess to understand whether the star had undergone any mass loss due to external merger events, He-flash event or binary interactions. We collected observed fluxes from near UV to far IR wavebands. The observed spectral energy distribution (SED) of KIC11087027 is shown in Fig 3(b). The observed SED is compared with the model SED of the star’s derived stellar atmospheric model taken from Kurucz grid of flux models and is constructed using VOSA999http://svo2.cab.inta-csic.es/theory/vosa/. Star shows significant IR-excess in near IR and far IR wavelengths, indicating an extended circumstellar environment due to episodic mass loss. Star also shows a double peak Hα emission profile as shown in Figure 3(a), which is asymmetric and variable, indicating the presence of chromospheric activity and mass motion in stellar atmosphere (Dupree et al., 1984).

3.5 SED modelling

The 1-D radiative code DUSTY (Ivezic et al., 1999) combined with the Kurucz model fluxes is used to derive dust parameters for a given set of atmospheric parameters and the assumed dust grain distribution. We have assumed a spherically symmetric dust shell with an oxygen-rich environment of warm silicate particles (Draine & Lee, 1984) and standard MRN grain size distribution from Mathis et al. (1977). A model SED is computed and compared with the observed SED (normalized to Ks-band fluxes) by varying dust parameters: inner dust temperature (Tin), optical depth (τν\tau_{\nu}); dust shell relative thickness y=(RoutRin)\rm\left(\frac{R_{out}}{R_{in}}\right).

We used a modified minimum chi-square test χ2\chi^{2}-test for the best fit.

χ2=(fobsfmodσobs)2/(Np1)\rm\chi^{2}=\sum\left(\frac{f_{obs}-f{mod}}{\sigma_{obs}}\right)^{2}/(N-p-1)

where fobs, fmod and σobs\sigma_{\rm obs} are observed flux, model flux and error in observed flux respectively. (N-p-1) represents the degrees of freedom, where N is the number of observed data points, and p is the number of independent parameters. Observed data is best fitted with two component dust shells: an inner hot and an outer cooler dust shell. The two-component dust shell best-fit model using DUSTY is shown in Fig.3(b).

The dust shell inner radii and expansion velocities have been obtained using scaling relations. Other parameters of dust shell: mass loss rates (M˙\rm\dot{M}), kinematic age (τd\tau_{\rm d}) and total mass (Md) are derived following Mallick et al. (2022) and are given in Table 2. The expected mass loss rate (M˙R\dot{M}_{\rm R}) is calculated using the modified Reimers (1975) formula for derived stellar parameters (see Table 2).

Stellar flares on a star’s surface may lead to distinct inner warm and outer cool dust shells, coupled with episodic mass loss (Osten & Wolk, 2015). Intense stellar radiation can lead to sublimation of dust grains closer to the star, creating an inner warm dust shell. The episodic events, shaped by flare variability may contribute to episodic mass loss from the two dust shells.

Refer to caption
Figure 3: (a): Hα emission in KIC11087027. Three spectra are taken at different times and resolutions: LAMOST LRS (R=1800), MRS and HESP (R=7500, 60000). At higher-resolution, the two Hα emission peaks become visible. (b): SED fitting to observed fluxes.
Table 2: Dust shell parameters of KIC11087027.
M˙R\dot{M}_{R} (M.yr1{}_{\odot}.yr^{-1}) Tinner τν\tau_{\nu} Td Vs RinR_{\rm in} td M˙\dot{M} (Myr1{}_{\odot}\rm yr^{-1}) Md
(Reimer’s Law) (K) (K) (km s-1) (cm) (years) (DUSTY) (M)
Inner Shell 6.68×\times10-11
500 0.09 193 5.48 4.87×\times1013 2.82 5.78×\times10-9 8.26×\times10-8
Outer Shell
300 10 99 1.41 2.45×\times1014 55.09 8.29×\times 10-7 2.32×\times 10-4

3.6 Rotational velocity

We estimated the star’s rotational velocity, Vrot, from spectra and photometric light curve data. The observed spectral width results from blending of other lines in the vicinity, instrumental broadening, and stellar phenomena like macroturbulence and rotation. We used two well-defined Fe I lines at 6703 and 6705 Å for extracting the rotation part of the line width by computing a set of profiles for given line abundance values, stellar parameters and instrumental broadening as measured from the Th-Ar calibration lines in the vicinity. The macroturbulence velocity, Vmac=3.9 km s-1, is calculated using calibration relation involving Teff and log g (Hekker & Meléndez, 2007). Using χ\chi2-test we found the observed profiles best fit for Vsini=10.1 km s-1. We also made use of Kepler photometric data to determine the period of the star. The data showed modulation in the light curve, probably due to large stellar spots that co-rotate with the stellar surface. We performed periodogram analysis by using the Lomb-Scargle method (Scargle, 1982) in the Kepler data spanning over four years. Periods ranging from 26 to 33 days during different quarters, as shown in the bottom panel of Fig.2(f), are found. The average period is found to be 30.38±\pm0.12 days. Using the estimated star’s radius, R=11.69±111.69\pm 1 R and considering this period as a rotational period (Prot), we found average VrotV_{\rm rot}=19.5±1.719.5\pm 1.7 km s-1, which is quite large for an evolved star like KIC11087027. The quarter-to-quarter variation, with a median absolute deviation of 2.1, in the rotation is probably caused by latitudinal differential rotation (Suto et al., 2022).

3.7 Binarity: Tidal synchronisation

We did not find any resolved binary of KIC11087027 within 10′′{}^{{}^{\prime\prime}} in the Gaia survey, and also, no signature of eclipse is present in the Kepler light curve. However, the presence of any unresolved companion can not be ruled out. The star has been observed at four different epochs in APOGEE, LAMOST, and HCT-HESP over a span of 7 years, and we found a small difference in the radial velocity (see Table 1). Gaia astrometry has been used for detecting unresolved companions of stars within 1′′{}^{{}^{\prime\prime}}, which can affect astrometric measurements. Gaia renormalized unit weight error (RUWE) value is greater than 1.4 for multiple system (Ziegler et al., 2020) and RUWE of KIC11087027 is 4.38. This inflated value of RUWE indicates a close unresolved binary or multiple system, as a wide binary does not inflate RUWE. This is further supported by the proper motion anomaly (PMa) found for KIC11087027. Proper motion values of KIC11087027 measured by Tycho-2 in the updated UCAC4 catalog (Zacharias et al., 2013) are μT=(2.1±0.7,5.0±0.5)\mu_{\rm T}=(-2.1\pm 0.7,-5.0\pm 0.5) mas yr-1 and Gaia DR3 are μG=(1.547±0.067,7.289±0.059)\mu_{\rm G}=(-1.547\pm 0.067,-7.289\pm 0.059) mas yr-1. Significant acceleration (σ=4.55\sigma=4.55) is found in the δ\delta direction with μTμG=(0.55±0.70,2.29±0.50)\mu_{\rm T}-\mu_{\rm G}=(-0.55\pm 0.70,2.29\pm 0.50), a signature that the star may be an unresolved binary or has multiple systems. The absence of oscillation modes (see Figure 2b) and the presence of strong magnetic activity in low mass fast rotating evolved star (Gaulme et al., 2014) suggests that the star KIC11087027 is probably in a tidally synchronised system. Other evidence is the star’s high projected rotational velocity (Vsini=10 km s-1) for its colour, (B-V)=1.04 (De Medeiros et al., 2002). The observations suggest KIC11087027 is tidally locked with an unresolved close binary companion with an orbital period of 30.4 days.

4 Discussion

Red giants with very high Li and strong IR-excess are rare because both properties seem transient and evolve differently. Though we do not have much information about IR-excess, the large surveys showed no evidence of IR-excess among giants (Bharat Kumar et al., 2015), and only a couple of Li-rich giants and SLRs are known to have IR-excess (Mallick et al., 2022). The general lack of IR excess among Li-rich giants indicates either no ejection of mass during the event that caused Li enrichment or the mass loss occurred, but the resultant IR-excess probably diluted faster than the Li depletion. In that sense, the star KIC11087027 is unusual, with very high Li and large IR-excess. This suggests Li-enrichment occurred recently, assuming the high Li and IR-excess result from a single event that is responsible for high Li among RC giants. Further, it is demonstrated that ΔΠ1\Delta\Pi_{1} is a proxy to the time evolution of a degenerate core into a convective core i.e., older RC giants with normal Li have relatively much higher ΔΠ1\Delta\Pi_{1} values (Singh et al., 2021). The estimated value of ΔΠ1\Delta\Pi_{1}=187 s from spectra Wang et al. (2023) indicates that KIC 11087027 is a young RC and Li-enrichment occurred very recently. The star also shows very high stellar rotation (Vrot=19.5 km s-1), which causes the star to have high chromospheric activity with flares. What process might have created all the above signatures in a single star is unclear. If these properties are a consequence of a single event like He-flash and evolve at different time scales, one would expect stars with all these properties to be uncommon. Here, we also note, the study by Sneden et al. (2022) which shows a subtle correlation between the strength of chromospheric line at He I 10830Å  and the high A(Li) in giants. They also show higher binary fraction among Li-rich giants compared to Li-poor implying some kind of binary interaction for Li production and chromospheric activity.

In the case of KIC 11087027, we discuss two plausible scenarios. One is an in-situ process in single-star evolution, and the other is the stars’ binary evolution, where orbital synchronisation causes high rotation aiding high Li abundance. It is highly likely that He-flash holds the key to high Li anomaly among red clump giants (Casey et al., 2019; Deepak & Reddy, 2019; Singh et al., 2019, 2021; Mallick et al., 2023). Of course, we need to understand the nature of some of the very high Li-rich giants reported in clusters. These giants appear to be on the RGB and well before the RGB tip (Sanna et al., 2020; Nagarajan et al., 2023) i.e before the He-flash event. It is important to analyse the cluster Li-rich giants using asteroseismology to find whether these are ascending the RGB. We will address this in the later study.

As per the RC Li-rich giants, the question remains about the physical process that mixes internally produced Li with the upper layers. Is He-flash the sole process? And does external aid, such as tidal locking, also play a role? Supporting the high levels of Li among red clump giants Schwab (2020) constructed models in which the powerful gravity waves generated during the He-flash provide a conducive environment for large-scale convection. For this, they required diffusion coefficient Dmix=1011cm2s-1, which is a few orders more than expected in stars ascending the RGB.

Let’s assume Li gets produced inside the star through the Cameron-Fowler mechanism (Cameron & Fowler, 1971), and the He-flash aids large-scale mixing. Though the models predict no observable feature like IR-excess due to the He-flash at the centre, it is not unreasonable to suggest mass loss occurs due to the flash, and the severity of mass loss depends on the intensity of the He-flash. The resultant IR-excess is diluted faster since it is a one-time event. This may be why red clump giants with IR-excess are rare. Importantly, all the SLRs or Li-rich giants don’t need to have IR-excess as they may evolve at different time-scales. The red clump giants with IR-excess are found more likely to be super Li-rich or Li-rich (Mallick et al., 2022) suggesting that IR-excess dilutes faster than Li depletion.

Concerning the high rotational velocity of KIC11087027, it is expected that the He-flash and the consequent reduction in star size by a factor of \sim10 will spin-up the star. Given the Vrot \sim0.7 km s-1 for the giants at the RGB tip one would expect Vrot \sim 7 km s-1 for a young RC giant. This expected value due to a sudden drop in the size of the giant from the RGB tip to the RC is much smaller than the observed value of 19.5 km s-1. This implies high rotation in KIC11087027 and other Li-rich RC stars found in the literature probably arose due to external factors like mergers, or tidal locking. A primary giant at or near the RGB tip, which is tidally locked with its binary star, may have a larger rotational velocity than a single star. In the case of KIC11087027, observed Vrot=19.5 km s-1 may be the combination of tidal locking and the after-effect of the He-flash.

If we assume the star in tidal locking with close binary has slightly more Vrot say 2 km s-1 rather the expected value of 0.7 km s-1, the resulting post He-flash Vrot value matches with the observed high velocity. The large Vrot will increase the diffusion coefficient to 1011 cm2 s-1 (Casey et al., 2019) as DmixΩ2{}_{\rm mix}\propto\Omega^{2}. Models show that a significant increase in the value of Dmix may result in the observed lithium abundance in the photosphere(Denissenkov & VandenBerg, 2003).

5 Conclusion

In this study, we reported results from Kepler photometric light curve and high-resolution spectral data. Based on variation in radial velocity and Gaia astrometric data, we suggest that KIC 11087027 has an unresolved companion. The observed high Vrot=19.5 km s-1 may be a result of a combination of tidal synchronization in a binary system and the effect of He-flash. The high rotation probably results in a high diffusion coefficient hence, the high Li abundance. The presence of several transient features like very high Li abundance, multiple dust shells, high rotation and strong flares provide evidence that either the star has undergone a merger-induced He-flash (Mallick et al., 2022) or the tidally locked primary star has undergone the He-flash very recently. Also, to have a comprehensive understanding of origin of Li enhanced giants it is imperative to understand the nature of Li-rich giants in clusters which are reported to be ascending the RGB.

6 Acknowledgement

We thank anonymous referees for useful comments which has improved our manuscript. This study is supported by the National Natural Science Foundation of China under grant No. 11988101, and the National Key R&\&D Program of China No.2019YFA0405500. Guoshoujing Telescope (the Large Sky Area Multi-Object Fiber Spectroscopic Telescope LAMOST) is a National Major Scientific Project built by the Chinese Academy of Sciences. LAMOST is operated and managed by the National Astronomical Observatories, Chinese Academy of Sciences. This work presents results from the European Space Agency (ESA) space mission Gaia. Gaia data are being processed by the Gaia Data Processing and Analysis Consortium (DPAC). Funding for the DPAC is provided by national institutions, in particular the institutions participating in the Gaia MultiLateral Agreement (MLA). The Gaia mission website is https://www.cosmos.esa.int/gaia. We gratefully acknowledge the entire team of Kepler space telescope which is funded by NASA’s Science Mission.

References

  • Andrae et al. (2018) Andrae, R., Fouesneau, M., Creevey, O., et al. 2018, A&A, 616, A8
  • Asplund et al. (2009) Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481
  • Bedding et al. (2011) Bedding, T. R., Mosser, B., Huber, D., et al. 2011, Nature, 471, 608
  • Bharat Kumar et al. (2015) Bharat Kumar, Y., Reddy, B. E., Muthumariappan, C., & Zhao, G. 2015, A&A, 577, A10
  • Borucki et al. (2010) Borucki, W. J., Koch, D., Basri, G., et al. 2010, Science, 327, 977
  • Cameron & Fowler (1971) Cameron, A. G. W., & Fowler, W. A. 1971, ApJ, 164, 111
  • Casey et al. (2019) Casey, A. R., Ho, A. Y. Q., Ness, M., et al. 2019, ApJ, 880, 125
  • Castelli & Kurucz (2004) Castelli, F., & Kurucz, R. L. 2004, ArXiv Astrophysics e-prints, astro-ph/0405087
  • Chaplin et al. (2011) Chaplin, W. J., Bedding, T. R., Bonanno, A., et al. 2011, ApJ, 732, L5
  • Cui et al. (2012) Cui, X.-Q., Zhao, Y.-H., Chu, Y.-Q., et al. 2012, Research in Astronomy and Astrophysics, 12, 1197
  • Cuntz et al. (2000) Cuntz, M., Saar, S. H., & Musielak, Z. E. 2000, ApJ, 533, L151
  • De Medeiros et al. (2002) De Medeiros, J. R., Da Silva, J. R. P., & Maia, M. R. G. 2002, ApJ, 578, 943
  • Deepak & Reddy (2019) Deepak, & Reddy, B. E. 2019, MNRAS, 484, 2000
  • Denissenkov & Herwig (2004) Denissenkov, P. A., & Herwig, F. 2004, ApJ, 612, 1081
  • Denissenkov & VandenBerg (2003) Denissenkov, P. A., & VandenBerg, D. A. 2003, ApJ, 593, 509
  • Dotter (2016) Dotter, A. 2016, ApJS, 222, 8
  • Draine & Lee (1984) Draine, B. T., & Lee, H. M. 1984, ApJ, 285, 89
  • Dupree et al. (1984) Dupree, A. K., Hartmann, L., & Avrett, E. H. 1984, ApJ, 281, L37
  • Fekel & Balachandran (1993) Fekel, F. C., & Balachandran, S. 1993, ApJ, 403, 708
  • Fouesneau et al. (2022) Fouesneau, M., Andrae, R., Sordo, R., & Dharmawardena, T. 2022, dustapprox, v0.1, , . https://github.com/mfouesneau/dustapprox
  • Gaia Collaboration et al. (2018) Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2018, ArXiv e-prints, arXiv:1804.09365
  • Gaia Collaboration et al. (2016) Gaia Collaboration, Prusti, T., de Bruijne, J. H. J., et al. 2016, A&A, 595, A1
  • Gao et al. (2022) Gao, J., Zhu, C., Yu, J., et al. 2022, A&A, 668, A126
  • Gaulme et al. (2014) Gaulme, P., Jackiewicz, J., Appourchaux, T., & Mosser, B. 2014, ApJ, 785, 5
  • Hawkins et al. (2018) Hawkins, K., Ting, Y.-S., & Walter-Rix, H. 2018, ApJ, 853, 20
  • Hekker & Meléndez (2007) Hekker, S., & Meléndez, J. 2007, A&A, 475, 1003
  • Ivezic et al. (1999) Ivezic, Z., Nenkova, M., & Elitzur, M. 1999, arXiv e-prints, astro
  • Jermyn & Fuller (2022) Jermyn, A. S., & Fuller, J. 2022, arXiv e-prints, arXiv:2206.13479
  • Kumar et al. (2020) Kumar, Y. B., Reddy, B. E., Campbell, S. W., et al. 2020, Nature Astronomy, 4, 1059
  • Kumar et al. (2011) Kumar, Y. B., Reddy, B. E., & Lambert, D. L. 2011, ApJ, 730, L12
  • Lind et al. (2009) Lind, K., Asplund, M., & Barklem, P. S. 2009, A&A, 503, 541
  • Mallick et al. (2022) Mallick, A., Reddy, B. E., & Muthumariappan, C. 2022, MNRAS, 511, 3741
  • Mallick et al. (2023) Mallick, A., Singh, R., & Reddy, B. E. 2023, ApJ, 944, L5
  • Mathis et al. (1977) Mathis, J. S., Rumpl, W., & Nordsieck, K. H. 1977, ApJ, 217, 425
  • Miller Bertolami et al. (2020) Miller Bertolami, M. M., Battich, T., Córsico, A. H., Christensen-Dalsgaard, J., & Althaus, L. G. 2020, Nature Astronomy, 4, 67
  • Nagarajan et al. (2023) Nagarajan, N., Sneden, C., Afsar, M., & Pilachowski, C. 2023, in American Astronomical Society Meeting Abstracts, Vol. 241, American Astronomical Society Meeting Abstracts, 367.01
  • Osten & Wolk (2015) Osten, R. A., & Wolk, S. J. 2015, ApJ, 809, 79
  • Paxton et al. (2011) Paxton, B., Bildsten, L., Dotter, A., et al. 2011, ApJS, 192, 3
  • Placco et al. (2021) Placco, V. M., Sneden, C., Roederer, I. U., et al. 2021, Research Notes of the American Astronomical Society, 5, 92
  • Reimers (1975) Reimers, D. 1975, Memoires of the Societe Royale des Sciences de Liege, 8, 369
  • Sanna et al. (2020) Sanna, N., Franciosini, E., Pancino, E., et al. 2020, A&A, 639, L2
  • Scargle (1982) Scargle, J. D. 1982, ApJ, 263, 835
  • Schwab (2020) Schwab, J. 2020, ApJ, 901, L18
  • Simon & Drake (1989) Simon, T., & Drake, S. A. 1989, ApJ, 346, 303
  • Singh et al. (2021) Singh, R., Reddy, B. E., Campbell, S. W., Kumar, Y. B., & Vrard, M. 2021, ApJ, 913, L4
  • Singh et al. (2019) Singh, R., Reddy, B. E., & Kumar, Y. B. 2019, MNRAS, 482, 3822
  • Sneden et al. (2022) Sneden, C., Afşar, M., Bozkurt, Z., et al. 2022, ApJ, 940, 12
  • Sneden (1973) Sneden, C. A. 1973, PhD thesis, THE UNIVERSITY OF TEXAS AT AUSTIN.
  • Suto et al. (2022) Suto, Y., Sasaki, S., Nakagawa, Y., & Benomar, O. 2022, PASJ, 74, 857
  • Tsantaki et al. (2018) Tsantaki, M., Andreasen, D. T., Teixeira, G. D. C., et al. 2018, MNRAS, 473, 5066
  • Tsantaki et al. (2023) Tsantaki, M., Delgado-Mena, E., Bossini, D., et al. 2023, A&A, 674, A157
  • Vida & Roettenbacher (2018) Vida, K., & Roettenbacher, R. M. 2018, A&A, 616, A163
  • Wallerstein & Sneden (1982) Wallerstein, G., & Sneden, C. 1982, ApJ, 255, 577
  • Wang et al. (2023) Wang, C., Huang, Y., Zhou, Y., & Zhang, H. 2023, A&A, 675, A26
  • Yan et al. (2021) Yan, H.-L., Zhou, Y.-T., Zhang, X., et al. 2021, Nature Astronomy, 5, 86
  • Zacharias et al. (2013) Zacharias, N., Finch, C. T., Girard, T. M., et al. 2013, AJ, 145, 44
  • Zhao et al. (2006) Zhao, G., Chen, Y.-Q., Shi, J.-R., et al. 2006, Chinese J. Astron. Astrophys., 6, 265
  • Zhao et al. (2012) Zhao, G., Zhao, Y.-H., Chu, Y.-Q., Jing, Y.-P., & Deng, L.-C. 2012, Research in Astronomy and Astrophysics, 12, 723
  • Ziegler et al. (2020) Ziegler, C., Tokovinin, A., Briceño, C., et al. 2020, AJ, 159, 19