Study of a red clump giant, KIC 11087027, with high rotation and strong infrared excess - Evidence of tidal interaction for high lithium abundance
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.40.1 days, which translates to rotational velocity Vrot=19.5 1.7 km s-1. Star’s location in the HR-diagram, derived values of = 71 and , 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.
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 () evolution. During transition phase from the degenerate core to the convective core He burning, the () value increases. As a result, very high A(Li) values are mostly seen at relatively lower () values, i.e., immediately after the He-flash. The normal Li RC giants are seen at high values. From this, they hypothesised that the super Li-rich (SLR) giant phase (A(Li)111A(Li)=12+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 2M⊙) giants and none among intermediate (M2M⊙) 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.

3 Analysis and results
3.1 Stellar parameters and Abundances
Stellar parameters (Teff, log , [Fe/H], Vsini, ), 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 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.350.10 dex and from the subordinate line at 6103 Å is A(Li)=3.850.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 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 , [Fe/H], ) and SNR. Abundances of other elements are derived using equivalent widths.
Name | KIC11087027 | reference |
---|---|---|
R.A. | 19:33:14.212 | |
DEC. | +48:41:50.83 | |
V (mag) | APASS | |
AG (mag) | this work | |
Parallax (mas) | Gaia DR3 | |
RUWE | Gaia DR3 | |
Prot (Days) | this work | |
Vrot(km s-1) | 19.5 | this work |
Vsini(km s-1) | this work | |
Inclination angle | 35∘ | this work |
Number of flares | 15 | this wok |
Teff [K] | this work | |
log | this work | |
this work | ||
(km s-1) | this work | |
this work | ||
Mass (M⊙) | this work | |
Radius (R⊙) | this work | |
RVGaia (km s-1) | Gaia DR3 | |
RVLAMOST (km s-1) | Luo+2018 | |
RVHESP,Helio (km s-1) | this work | |
A(Li)NLTE(6103 Å) | this work | |
A(Li)NLTE(6707 Å) | this work | |
this work | ||
this work | ||
this work | ||
this work | ||
this work | ||
this work | ||
this work | ||
this work | ||
this work | ||
this work | ||
this work | ||
this work | ||
this work | ||
this work | ||
this work | ||
this work | ||
this work | ||
this work | ||
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 . 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.950.22 and [Fe/H]= -0.550.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 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 s and 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 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 () and Kepler response function ().
Here and are the lower and upper wavelength limits of the Kepler filter. Bolometric flux for stellar parameters is taken from VOSA 777http://svo2.cab.inta-csic.es and Kepler response function from Kepler website 888https://keplergo.github.io/KeplerScienceWebsite/. Energy of flares is in the range of erg, which makes KIC11087027 as super flaring star.

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 (); dust shell relative thickness y=.
We used a modified minimum chi-square test -test for the best fit.
where fobs, fmod and 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 (), kinematic age () and total mass (Md) are derived following Mallick et al. (2022) and are given in Table 2. The expected mass loss rate () 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.

(M) | Tinner | Td | Vs | td | (M) | Md | |||
---|---|---|---|---|---|---|---|---|---|
(Reimer’s Law) | (K) | (K) | (km s-1) | (cm) | (years) | (DUSTY) | (M⊙) | ||
Inner Shell | 6.6810-11 | ||||||||
500 | 0.09 | 193 | 5.48 | 4.871013 | 2.82 | 5.7810-9 | 8.2610-8 | ||
Outer Shell | |||||||||
300 | 10 | 99 | 1.41 | 2.451014 | 55.09 | 8.29 10-7 | 2.32 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 2-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.380.12 days. Using the estimated star’s radius, R= R⊙ and considering this period as a rotational period (Prot), we found average = 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 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, 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 mas yr-1 and Gaia DR3 are mas yr-1. Significant acceleration () is found in the direction with , 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 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 values (Singh et al., 2021). The estimated value of =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 10 will spin-up the star. Given the Vrot 0.7 km s-1 for the giants at the RGB tip one would expect Vrot 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 D. 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 RD 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