Abstract
In order to clarify the properties of the secondary clump star HD 226808 (KIC 5307747), we combined four years of data from Kepler space photometry with high-resolution spectroscopy of the High Efficiency and Resolution Mercator Échelle Spectrograph mounted on the Mercator telescope. The fundamental atmospheric parameters, radial velocities, rotation velocities, and elemental abundance for Fe and Li were determined by analyzing line strengths and fitting-line profiles, based on a 1D local thermodynamic equilibrium model atmosphere. Second, we analyzed a photometric light curve obtained by Kepler and we extracted asteroseismic data of this target using Lets Analysis, Use and Report of Asteroseismology, a new seismic tool developed for the study of evolved FGK solar-like stars. We determined the evolutionary status and effective temperature, surface gravity, metallicity, microturbulence, and chemical abundances for Li, Ti, Fe, and Ni for HD 226808, by employing spectroscopy, asteroseismic scaling relations, and evolutionary structure models built in order to match observed data. Our results also show that an accurate synergy between good spectroscopic analysis and asteroseismology can provide a jump toward understanding evolved stars.
Export citation and abstract BibTeX RIS
1. Introduction
Red giants are cool evolved stars with an extended convective envelope, which can, as in main-sequence solar-like stars, stochastically excite modes of oscillation. This stellar evolutionary phase is well studied observationally and theoretically but is still quite intriguing, because after the exhaustion of H in the core, several structural changes occur in the interior of the star during this phase. In fact, the evolution is proceeded by fusion of hydrogen in a shell surrounding an inert core of helium, which continues to increase its mass and contracts while the stellar radius expands. If the stellar core reaches a temperature of ∼100 million K, the fusion of helium through the triple alpha process can start. However, the initial stellar mass will play a defining role in the evolution of the star after helium begins to burn in the core. Stars with masses below ∼1.9M⊙ develop degenerate helium cores that increase their mass while the H is burning in the shell, until the core reaches a critical mass of Mc ∼ 0.47M⊙.
At this point, a helium flash occurs abruptly in the center of the helium core and the degeneracy is gradually lifted. Those stars, due to the degenerate condition of the core, have similar core masses and hence similar luminosities. Thus, they are found in a narrow region of the H–R diagram, forming the so-called "red clump" (RC; e.g., Cannon 1970; Faulkner & Cannon 1973; Girardi et al. 1998; Girardi 1999). On the other hand, stars with masses higher than ∼1.9M⊙ will not develop a degenerate helium core in the post-main-sequence phase. They will not have an explosive liberation of energy like the one seen during the helium flash. Instead, they will have a gradual ignition of helium. Therefore, those stars will occupy a region of the H–R diagram similar to the RC, but at lower luminosities and a wider spread. This group, which includes more massive core helium-burning stars, at ages of about 1 Gyr, is referred as the secondary clump, and its components are the second clump stars (see, e.g., Girardi 1999).
Hawkins et al. (2018) studied in detail the spectroscopic signatures of genuine core helium-burning RC giant stars and how they compare to shell hydrogen-burning red giant branch (RGB) stars. Additionally, spectroscopic studies of genuine RGB or RC stars have been done by Anders et al. (2016) for some APOGEE and CoRoT targets. Red clump stars are natural standard candles (Stanek et al. 1998; Hawkins et al. 2018), while regular RGB stars or more massive secondary RC stars of nearly the same effective temperatures (Teff) are not. In this context, finding and characterizing core helium-burning RC stars is very important to fine-tune stellar evolution and Galactic archeology, and consequently to build a precise cosmic distance ladder (Stanek et al. 1998; Bressan et al. 2013; Bovy et al. 2014; Gontcharov 2017; Hawkins et al. 2018). Distinguishing RC stars from less evolved shell hydrogen-burning RGB stars or more massive secondary RC stars is a critical point and an obstacle to solving some stellar astrophysics problems (Bressan et al. 2013).
Fortunately, continuous photometric observations from space missions such as CoRoT (Baglin et al. 2006) and Kepler (Borucki et al. 2009) provide photometric data of unprecedented quality. This gives us the possibility of better distinguishing the evolutionary status of red giants along the H shell-burning phase before He ignition from those on the He-burning phase after He ignition, as described by Bedding et al. (2011). In particular, using asteroseismic scaling relations calibrated on solar values, it has been possible to determine with good accuracy the stellar mass (M) and radius (R; e.g., Casagrande et al. 2014; Pinsonneault et al. 2014) of hundreds of observed red-giant stars. Much progress has been made during the last decade and asteroseismology provides a way forward to classify RC based on stellar natural oscillation frequencies. However, at present, some RC stars, classified so far as authentic secondary clump giant stars (Mosser et al. 2014), lack a high-resolution spectroscopic study, hence the importance of this analysis.
In this study, we use photometric data from the Kepler space telescope and High Efficiency and Resolution Mercator Échelle Spectrograph (HERMES) ground-based high-resolution spectroscopy to produce a deep analysis of the bright red-giant star HD 226808 (KIC 530774). In particular, this star is one of three brightest objects classified as a secondary clump star and observed by Kepler in long-cadence mode (Mosser et al. 2014). This paper is organized as follows. In Section 2, we present the spectroscopic observations and spectroscopic derived fundamental parameters. In Section 3 we present an asteroseismic analysis. In Section 4 we compare the asteroseismic results with evolutionary models. Our conclusions are presented in Section 5.
2. Spectroscopic Observations
Our target, the star HD 226808 (KIC 530774), with V = 8.67 mag, was observed on 2015 July 3 with the HERMES (Raskin et al. 2011; Raskin 2011) mounted on the 1.2 m Mercator Telescope at the Observatorio del Roque de los Muchachos on La Palma, Canary Islands, Spain. The HERMES spectra covers a wavelength range between 375 and 900 nm with a spectral resolution of R ≃ 85,000. The wavelength reference was obtained from emission spectra of Thorium–Argon–Neon reference frames in close proximity to the individual exposure. The standard steps of the spectral reduction were performed with an instrument-specific pipeline as described by Raskin et al. (2011) and Raskin (2011). The radial velocity (RV) for each individual spectrum was determined from the cross-correlation of the stellar spectrum in the wavelength range between 478 and 653 nm with a standard mask optimized for Arcturus provided by the HERMES pipeline toolbox. For HERMES, the 3σ level of the night-to-night stability for the observing mode described above is ∼300 m s−1, which is used as the classical threshold for RV variations to detect binarity. We corrected individual spectra for the Doppler shift before normalization and combined individual spectra as described in Beck et al. (2016). The total integration time was 1.7 hr, split into four equal parts of 1500 s each. The final spectrum, stacked from individual exposures, shows a signal-to-noise ratio ≃ 150 around the Li iλ 670 nm line.
2.1. Fundamental Parameters from HERMES Spectroscopy
For the fundamental parameters, the analysis of HD 226808 starts by using as a first guess fundamental parameters (effective temperature Teff; surface gravity metallicity [Fe/H] and microturbulence) from the revised Kepler Input Catalog (KIC) by Huber et al. (2014). Next, we used the ARES code (v2) (Sousa et al. 2007, 2015) to measure equivalent widths (EW) of selected spectra absorption Fe i and Fe ii lines, and with the q2 code Ramírez et al. (2014) we determined the fundamental physical parameters based on a line list, as described by Ramírez et al. (2014). Then, we refine the process based on one-dimensional (1D) local thermodynamic equilibrium model atmosphere as described by Beck et al. (2016). We also considered the Arcturus and solar parameters as other reference values, as described by Hinkle & Wallace (2005), Ramírez & Allende Prieto (2011), Meléndez et al. (2012), Monroe et al. (2013), and Ramírez et al. (2014). Final spectroscopic parameters, such as the values of effective temperature , surface gravity , metallicity [Fe/H], and A(Li) of HD 226808 are given in Table 2. Typically, stellar parameter uncertainties were computed as described by Epstein et al. (2010) and Bensby et al. (2014) and are , ±0.08, ±0.04 and ±0.04 for Teff, log g, [Fe/H], and , respectively. For abundance determinations, we used spectral syntheses based on a combination of a 2014 version of the code MOOG (Sneden 1973) with Kurucz atmosphere models (Castelli & Kurucz 2004) and line lists as described in Table 1 (for fundamental parameter determinations, we used a line list as described in Table B1). A low lithium abundance signature is found for this secondary RC giant star. However, on the same evolutionary stage, some lithium-rich stars were found by Bharat Kumar et al. (2018). The values that we obtained for effective temperature and metallicity are comparable with those by Takeda & Tajitsu (2015), who observed this star with the High Dispersion Spectrograph at the Subaru telescope. The resulting value of is slightly higher than the value obtained by Takeda & Tajitsu (2015), but within the error bar. In Figure 1 we show part of the spectrum for HD 226808 in the region around 5086 Å. We also show results for two different synthetic spectra computed with the code MOOG for the same spectroscopic fundamental parameters. For the error analysis in the abundance determination we preserve the same fundamental parameters and slightly vary the abundances along a step size from 0.05 until 0.10 dex. The values are shown in Table 2.
Table 1. Line List Used for HD 226808 (KIC 530774)
Wavelength | Species | χexc | log (gf) | EW |
---|---|---|---|---|
(Å) | (eV) | (Å) | ||
5085.477 | 28.0 | 3.658 | −1.541 | 94.12 |
5085.676 | 26.0 | 4.178 | −2.610 | 9.78 |
5085.837 | 26.0 | 4.495 | −3.881 | 9.27 |
5085.933 | 26.0 | 3.943 | −3.151 | 9.77 |
5086.231 | 606.0 | 0.252 | −0.110 | 24.16 |
5086.251 | 26.0 | 4.988 | −2.624 | 24.16 |
5086.259 | 606.0 | 0.252 | −0.122 | 24.16 |
5086.390 | 606.0 | 0.252 | −0.133 | ⋯ |
5086.765 | 25.0 | 4.435 | −1.324 | 36.07 |
5086.805 | 606.0 | 0.301 | −0.649 | 36.07 |
5086.997 | 10108.0 | 0.040 | −10.286 | 75.40 |
5087.058 | 22.0 | 1.430 | −0.780 | 75.43 |
5087.065 | 606.0 | 0.301 | −0.684 | 75.43 |
Note. For species we are using MOOG standard notation for atomic or molecular identification. For example, 26 represents Fe(26), while the 0 after the decimal indicates neutral and 1 indicates singly ionized. χexc is the line excitation potential in electron volts (eV).
Download table as: ASCIITypeset image
Table 2. Fundamental Parameter Data of HD 226808 (KIC 530774)
R.A. (J2000.0) | |
---|---|
Decl. (J2000.0) | +40°37'161 |
Kp (mag)a | 8.397 |
Spectroscopically Derived Parameters | |
Teff (K) | 5065 |
log g (cgs) | 2.99 ± 0.08 |
Vt (km s−1) | 1.27 ± 0.04 |
[Fe/H]b | 0.08 ± 0.04 dex |
A(Li) | ≤0.58 |
A(C) | 8.35 ± 0.05 dex |
A(Ti) | 4.82 ± 0.06 dex |
A(Fe) | 7.58 ± 0.07 dex |
A(Ni) | 6.58 ± 0.05 dex |
Notes.
aKp is the white light Kepler magnitude taken from the revised KIC (Huber et al. 2014). bAbundances assuming Asplund et al. (2009).Download table as: ASCIITypeset image
3. Analysis of Asteroseismic Data
We have developed a code with the purpose of analyzing stellar oscillations and extracting seismic data from light curves. The code makes use of version 3.6 of the Python programming language and has been calibrated using several stars from the literature, especially the giant stars from Ceillier et al. (2017) and the KIC 4448777 from Di Mauro et al. (2018). The Lets Analysis, Use and Report of Asteroseismology (LAURA; De Moura & De Almeida 2018) searches the Pre-search Data Conditioning Simple Aperture Photometry (PDC-SAP) flux light curves retrieved from the Mikulski Archive for Space Telescopes (MAST) or the Kepler Asteroseismic Science Operations Center (KASOC) databases and concatenates them in order to reduce noise in the data using a high rate of oversample to maintain temporal stability by increasing the resolution of the points. After that, the code proceeds to do a cleanup of excessive high values of frequencies thanks to a box-car type filter of 2 μHz. The time series extracted can then be treated by employing two different tools in order to characterize the periodicity of the star: Lomb–Scargle (LS) periodogram (Scargle 1982) and minimization of phase dispersion (PDM) (Stellingwerf 1978). Figure 2 shows the results obtained for the star HD 226808 (KIC 5307747) using the method of minimization of phase dispersion and showing that the smallest phase occurs for a rotation period of 122 days, in good agreement with the result obtained by the LS periodogram shown in Figure 3. Our method, based on the use of LS and PDM together, has given results consistent with the values found by Ceillier et al. (2017), who applied a different technique to a set of red giants, confirming the potential of our tool to determine the rotational period. The main packages of the code are presented in detail in Appendix A.
Download figure:
Standard image High-resolution imageDownload figure:
Standard image High-resolution imageThe same code has been implemented to perform the Fast Fourier Transform of the time series to analyze the stellar oscillation power spectra and study the excess of power at high frequency. The global approach in the seismic study that involves the adjustment of the entire spectrum has challenged us to compare with other already consolidated codes. From that comparison we finally can conclude that LAURA is in good agreement with other ones using similar methods. For example, in Figure 4 we show a comparison of the background modeling methods between LAURA and DIAMONDS (Corsaro et al. 2015) developed for the analysis of HD 226808. We found that the photometric granulation signatures due to intense convective activities on the surface present only subtle differences and do not affect the seismic analysis (Mathur et al. 2011).
Download figure:
Standard image High-resolution image3.1. Results
The power spectrum of the light curve for HD 226808 shows a clear excess of power, with a typical Gaussian distribution profile, in the range 60–120 μHz due to the radial modes and some additional dipolar and quadrupolar modes, as is shown in Figure 5. Based on our method, we determined the frequency at maximum oscillation power μHz and the so-called large frequency separation between modes with the same harmonic degree μHz. The quoted uncertainties are based on a local minimum optimization strategy called the Levenberg–Marquardt9 method in a specific package, allowing an estimate of the following differences:
where , , and i are the set of measured data, the model calculation, and estimated uncertainty in the data, respectively. Using the universal method for the red-giant oscillation pattern by Mosser et al. (2011, 2017) it was possible to identify pure-pressure modes and the identification of several dipole modes by a grid-search allowed us to estimate the period spacing ΔΠ = 296.85 ± 6.53 s, which places this star on the secondary clump red-giant phase (Mosser et al. 2014). The asteroseismic fundamental parameters found here agree with the results published by Mosser et al. (2012, 2014), who found for this star νmax = 88.2 μHz, and ΔΠ = 296.3 ± 5.9 s, while their large separation, which is Δν = 6.88 ± 0.05 μHz agrees within the errors. Table 3 lists the final set of frequencies for the detected individual modes together with their uncertainties.
Download figure:
Standard image High-resolution imageTable 3. Frequencies and Harmonic Degrees for Observed Oscillation Modes Found Using the LAURA Script
l | νn, l |
---|---|
0 | 66.34916 ± 0.3364 |
0 | 73.23633 ± 0.0292 |
0 | 80.02904 ± 0.2785 |
0 | 93.56723 ± 0.0614 |
1 | 62.81494 ± 0.0205 |
1 | 63.38166 ± 0.0194 |
1 | 68.50571 ± 0.0305 |
1 | 75.96745 ± 0.0583 |
1 | 76.06977 ± 0.0306 |
1 | 80.02891 ± 0.0415 |
1 | 82.29577 ± 0.0455 |
1 | 83.05140 ± 0.0865 |
1 | 83.22456 ± 0.0147 |
1 | 83.35049 ± 0.0611 |
1 | 89.19081 ± 0.0136 |
1 | 89.87559 ± 0.0550 |
1 | 90.06449 ± 0.0399 |
1 | 96.03074 ± 0.0325 |
2 | 72.31541 ± 0.3473 |
2 | 79.10025 ± 0.0461 |
2 | 92.04025 ± 0.1325 |
3 | 67.30144 ± 0.0958 |
3 | 74.47983 ± 0.0431 |
3 | 87.91571 ± 0.0326 |
Note. Uncertainties are measured by Lorentzian (l = 0, 2, 3) and sinc fits (l = 1).
Download table as: ASCIITypeset image
3.2. Fundamental Parameters from Asteroseismic Scaling Laws
The asteroseismic surface gravity, stellar mass, and radius for this star were obtained from the observed Δν and νmax together with the spectroscopic value of following the scaling relations calibrated on solar values as provided by Kjeldsen & Bedding (1995), Kallinger et al. (2010), and Belkacem et al. (2011):
and
where the values for the Sun are Δν⊙ = 134.9 μHz, μHz, and . We obtained , , and .
By comparison, the ratio obtained in the spectroscopic analysis had a value of ∼0.036. These values agree within the errors from asteroseismic radius and mass values computed by Mosser et al. (2012) and Takeda & Tajitsu (2015) and are based on different estimations of maximum power frequency.
4. Comparison with Evolutionary Models
Given the identified pulsation frequencies and the basic atmospheric parameters, we faced the theoretical challenge of interpreting the observed oscillation modes by constructing stellar models that satisfy the spectroscopic and asteroseismic observational constraints. We assumed the effective temperature and surface gravity as calculated in Section 2.1, respectively K, dex and metallicity (see Table 2).
We computed theoretical evolutionary models representative of the present star using the MESA evolutionary code (Paxton et al. 2011), in which we varied the stellar mass and the initial chemical composition in order to match the available fundamental parameters. Figure 6 shows two evolutionary tracks obtained with masses 2.5M⊙ and 2.8M⊙ and fixed initial composition. We used an initial helium abundance of Y0 = 0.27 and initial hydrogen abundance X0 = 0.71 as input to MESA and plotted in an effective temperature–gravity evolutionary diagram. The location of HD 226808 is indicated in Figure 6, which identifies this star as being in the CHeB phase of evolution (see Section 1), with a mass between 2.5 and 2.8M⊙. We point out that the location of evolutionary tracks in the diagram strongly depends on the metallicity of the model, since it influences the opacity and thus the depth of the convective envelope. If a solar metallicity would have been assumed instead, as we can see in the insert in Figure 6, the star would have appeared to be in the RGB phase, burning hydrogen in a shell around the inert He core. As a result, a good determination of the measured metallicity is essential to achieve a correct determination of the evolutionary status of stars located in this region of the HR diagram. The A(Li) also confirm the evolved status of this star. In order to investigate the observed solar-like oscillations, we used the GYRE package (Townsend & Teitler 2013) to compute adiabatic oscillation frequencies with degree for some of the models satisfying the spectroscopic constraints. In Figure 7, the open symbols represent the best output of the frequency model. Theoretical oscillation frequencies have been corrected by a near-surface effects term following the approach proposed by Ball & Gizon (2014), which has been proven to work much better for evolved stars than other prescriptions. The adjustment of the frequencies by a corrective term is a common procedure used for evolved stellar models to overcome the lack of a proper theory for the description of oscillations in the upper surface layers, where frequency behaviors deviate from the adiabatic assumption. Figure 8 shows the difference between modeled and observed frequencies for our best model with and without surface corrections, demonstrating the importance of such corrections. Among our models we selected one that appeared to best fit the observed frequencies of HD 226808, and whose characteristics are given in Table 4. This model is characterized by an age of 0.69 Gyr, a mass M = 2.6 M⊙, and a radius R = 9.74 R⊙. The value of the radius obtained by direct modeling of individual frequencies is now better constrained, and agrees within the error bars obtained using a scaling relation. Furthermore, the value of the mass derived from the evolutionary track resulted in good agreement with respect to the seismic value, obtained using the scaling relation and estimated in Section 3.2. Large discrepancies of up to 50% between masses inferred from both methods were pointed out by Miglio et al. (2012) and Takeda & Tajitsu (2015), who explored large sets of red-clump stars, and attributed to the overlap of the evolutionary tracks of RG stars, with higher masses, and RC stars, with lower masses, at the same clump region of the HR diagram. This difference does not appear in the case of H shell-burning red giants.
Download figure:
Standard image High-resolution imageDownload figure:
Standard image High-resolution imageDownload figure:
Standard image High-resolution imageTable 4. Parameters of the Best Model Using MESA for HD 226808 (KIC 5307747)
Parameters | Best model |
---|---|
M/M⊙ | 2.60 |
R/R⊙ | 9.74 |
A(Li) | 0.58 |
3.695 | |
(cgs) | 2.856 |
Age (Gyr) | 0.69 |
ΔP(s) | 221 |
Download table as: ASCIITypeset image
Several tests considering interferometry (e.g., Huber et al. 2012; White et al. 2013), Hipparcos parallaxes (Silva Aguirre et al. 2012), eclipsing binaries (e.g., Frandsen et al. 2013; Huber 2015; Gaulme et al. 2016), and open clusters (e.g., Miglio et al. 2012, 2016; Stello et al. 2016) have demonstrated that mass estimates from asteroseismic scaling relations are accurate to a few percent for main-sequence stars, while larger discrepancies have been identified for evolved stars. Theoretical studies have suggested corrections to scaling relations by, for example, comparing the large frequency separation (Δν) calculated from individual frequencies with models (e.g., Stello et al. 2009; White et al. 2011; Guggenberger et al. 2016; Sharma et al. 2016) or an extension of the asymptotic relation (Mosser et al. 2013). For the moment, it remains clear that scaling relation corrections should mainly depend on Teff and evolutionary state. However, it is yet unclear how these corrections should be applied in the case of RC stars (e.g., Miglio et al. 2012; Sharma et al. 2016).
5. Conclusion
In this study we perform an asteroseismic analysis of the red giant clump star HD 226808 (KIC 530774) based on its high-precision space photometry and high-resolution ground spectroscopic observations. The spectroscopic analysis was important to understand the metallicity of the star, which, as discussed in Section 4, has an impact on the evolutionary state of the star. Abundances obtained from the fundamental parameter solution, such as, for example, A(Li), produced good agreement between observational and synthesized values. We used the Kepler photometric data as input for seismic diagnostic tools. We characterized the internal structure and evolutionary status of this secondary clump red giant star, and studied its frequencies, especially since this star is one of the brightest star of the Kepler field. We have shown that accurate mass determination of RC stars with such good agreement between methods is possible. The metallicity effect on oscillation modes for this class of stars is still poorly understood. Future developments will be welcome to better understand this point, especially when comparing giants with very close evolutionary states. Finally, the presented method is robust enough to conduct seismic characterization of giants and to obtain measurements of rotational period. Our tools are open source and ready for use.
We are grateful to the "National Council for Scientific and Technological Development" (CNPq, Brazil) for support. J.D.N. is supported by the CNPq Brazil PQ1 grant 310078/2015-6. P.G.B. acknowledges the support of the MINECO under the fellowship program "Juan de la Cierva incorporacion" (IJCI-2015–26034). This work is based on observations obtained with the HERMES spectrograph on the Mercator Telescope, which is supported by the Research Foundation—Flanders (FWO), Belgium, the Research Council of KU Leuven, Belgium, the Fonds National de la Recherche Scientifique (F.R.S.-FNRS), Belgium, the Royal Observatory of Belgium, the Observatoire de Genéve, Switzerland, and the Thüringer Landessternwarte Tautenburg, Germany. Funding for the Kepler mission is provided by NASA.
Software: GYRE (Townsend & Teitler 2013), astropy (The Astropy Collaboration 2013, 2018), matplotlib (Hunter 2007), Numpy (Walt et al. 2011), Scipy (Jones et al. 2001), MOOG (Sneden 1973), ARES v2 (Sousa et al. 2007, 2015).
Appendix A: Master–Laura
A.1. The Algorithm
LAURA is written in Python 3, with five different packages to download, reduce, and analyze raw light curves under the asteroseismology theory. The code uses an object-oriented approach and can be easily used even for non-experts in program language. All the additional python packages that we use are free and easy to install using the pip package. Our code is summarized in two files: LAURA_master.py and LAURA_aux.py. From the LAURA_master.py file, the user can perform all steps to download, reduce, and analyze the light curve.
A.2. Installation
In order to use our code, make sure you have Python 3.x installed with the follow packages: Astropy (Astropy Collaboration et al. 2013), Scipy (Jones et al. 2001), kplr (Foreman-Mackey 2018), peakutils (Lucas Hermann Negri & Vestri 2017), and lmfit (Newville & Otten 2018). Installing the necessary packages in Ubuntu GNU/Linux is done as follows:
sudo apt-get install build-essential |
sudo apt-get install python3 |
sudo apt-get install python3-pip |
sudo pip install numpy |
sudo pip install matplotlib |
sudo pip install scipy |
sudo pip install astropy |
sudo pip install kplr |
sudo pip install peakutils |
sudo pip install lmfit |
Download table as: ASCIITypeset image
export PYTHONPATH=/home/USER_NAME/LAURA-1.0.0/source:$PYTHONPATH |
Download table as: ASCIITypeset image
A.3. Packages
LAURA runs from the box with five principal packages from the acquisition of the light curve to the Echelle diagram as follows:
import LAURAZaux as aux |
aux.TS_PS(ID,Qtin,Qtfin,cleaning,oversample,path="TEMP/") |
aux.FIT_PS(ID) |
aux.LS(ID,kernerforca,Pbeg,Pend,PeriodMod,path='''TEMP/'') |
aux.OBS(ID) |
aux.MOD(ID, largeSep,smallSep,size,Radial_Orders=5) |
aux.ECHELLE(ID, deltav) |
print (''Done!'') |
Download table as: ASCIITypeset image
- 1.ID: Identification number of star.
- 2.Qtin: Inicial quarter for Kepler objects.
- 3.Qtfin: Final quarter for Kepler objects.
- 4.cleaning: Cleaning parameter for systematic errors.
- 5.Oversample: The resolution of the power spectrum.
- 6.Path: The folder in which the .fits files will be download.
When running this package, LAURA will download all the .fits files from all the quarters you select. If you already have the object on your computer, the code will not download it again. A power spectrum will be generated and all files will be stored in a folder with the name of the object.
The FIT_PS package fits a background tendency to the power spectrum generated earlier using only the ID of the object. The first plot will ask you to select the region to fit a Gaussian curve, as shown in Figure 9 following
Download figure:
Standard image High-resolution imageA is the amplitude, μ is the center, and σ is sigma.
Then, you will be asked to select two consecutive Lorentzian fits by clicking the regions of interest. Each click stores a position to use as parameters of the following function:
All three models (one Gaussian and two Lorentzian) are then used to create a background fit for the power spectrum as shown in Figure 4.
The LS package runs a kernel to smooth out the light curve and runs a periodgram either using Lomb–Scargle or PDM. The input parameters are
- 1.ID: Identification number of star.
- 2.Kernel: Strengh of the Gaussian 1D kernel.
- 3.Pin: Initial period of search.
- 4.Pend: Final period of search.
- 5.PeriodMod: Either LS (Lomb–Scargle) or PDM (Phase Dispersion Minimization).
- 6.Path: Where to save the data.
This package will save all statistical results in a file named kic_periodgram.txt.
The OBS package returns all the observed modes from the high-frequency region knowing only the ID of the object. First, the user will be asked to select the high-frequency region from the power spectrum. Then it can select the value of the threshold and the strength of the box kernel. The code then selects all the peaks using a Lorentzian peak selector and saves all the output data to the _data_mods_obs.txt.
The MOD package makes use of the universal pattern (Mosser et al. 2012) to model the possible frequencies of the pure-pressure modes of spherical degrees for a defined range of radial orders. It takes as input:
- 1.ID: Identification number of star.
- 2.LargeSep: large Frequency separation.
- 3.SmallSep: small Frequency separation.
- 4.Size: Size of the region which contains each mode.
- 5.Radial_Orders: how many radial orders to calculate.
The code plots the high-frequency region from the observed data, highlighting the regions l = 0 and l = 1. This package will then save all model modes in a file named_data_mods_model.txt that can be used to compute the Echelle diagram in the future.
For comparison, the example star mentioned here was fully analyzed in 22 minutes: 1 minute to download the light curves, around 20 minutes to process the raw light curves and create the power spectrum, and less than 1 s to automatically output the rotation period and asteroseismic parameters. Once we have the data downloaded and the power spectrum created, the time consumed relies basically on the user interaction with the code, as LAURA is not fully automatic and needs the user to click and select regions in certain steps of the run. Thus, our code cannot, for now, be used automatically to analyze multiple stars without user interaction.
Appendix B: Line List
Table B1 shows the full line list Used for HD 226808 (KIC 530774).
Table B1. Full Line List Used for HD 226808 (KIC 530774)
Wavelength | Species | χexc | log (gf) | EW (Sun) | EW (HD 226808) |
---|---|---|---|---|---|
(Å) | (eV) | (Å) | (Å) | ||
5295.3101 | 26.0 | 4.42 | −1.59 | 29.0 | 50.9 |
5379.5698 | 26.0 | 3.69 | −1.51 | 62.5 | 92.3 |
5386.3301 | 26.0 | 4.15 | −1.67 | 32.6 | 59.5 |
5441.3398 | 26.0 | 4.31 | −1.63 | 32.5 | 58.3 |
5638.2598 | 26.0 | 4.22 | −0.77 | 80.0 | 102.7 |
5679.0229 | 26.0 | 4.65 | −0.75 | 59.6 | 78.7 |
5705.4639 | 26.0 | 4.30 | −1.35 | 38.0 | 62.1 |
5731.7598 | 26.0 | 4.26 | −1.20 | 57.7 | 83.7 |
5778.4531 | 26.0 | 2.59 | −3.44 | 22.1 | 60.7 |
5793.9141 | 26.0 | 4.22 | −1.62 | 34.2 | 61.2 |
5855.0762 | 26.0 | 4.61 | −1.48 | 22.4 | 43.3 |
5905.6699 | 26.0 | 4.65 | −0.69 | 58.6 | 78.0 |
5927.7900 | 26.0 | 4.65 | −0.99 | 42.9 | 63.9 |
5929.6802 | 26.0 | 4.55 | −1.31 | 40.0 | 62.6 |
6003.0098 | 26.0 | 3.88 | −1.06 | 84.0 | 109.6 |
6027.0498 | 26.0 | 4.08 | −1.09 | 64.2 | 91.9 |
6056.0000 | 26.0 | 4.73 | −0.40 | 72.6 | 92.2 |
6079.0098 | 26.0 | 4.65 | −1.02 | 45.6 | 68.2 |
6093.6440 | 26.0 | 4.61 | −1.30 | 30.9 | 55.4 |
6096.6650 | 26.0 | 3.98 | −1.81 | 37.6 | 64.4 |
6151.6182 | 26.0 | 2.17 | −3.28 | 49.8 | 93.0 |
6165.3599 | 26.0 | 4.14 | −1.46 | 44.8 | 72.4 |
6187.9902 | 26.0 | 3.94 | −1.62 | 47.6 | 77.3 |
6240.6460 | 26.0 | 2.22 | −3.29 | 48.2 | 91.9 |
6270.2251 | 26.0 | 2.86 | −2.54 | 52.4 | 92.2 |
6703.5669 | 26.0 | 2.76 | −3.02 | 36.8 | 80.3 |
6705.1021 | 26.0 | 4.61 | −0.98 | 46.4 | 73.5 |
6713.7451 | 26.0 | 4.79 | −1.40 | 21.2 | 39.6 |
6726.6670 | 26.0 | 4.61 | −1.03 | 46.9 | 68.8 |
6793.2588 | 26.0 | 4.08 | −2.33 | 12.8 | 32.4 |
6810.2632 | 26.0 | 4.61 | −0.97 | 50.0 | 73.7 |
6828.5898 | 26.0 | 4.64 | −0.82 | 55.9 | 79.5 |
6842.6899 | 26.0 | 4.64 | −1.22 | 39.1 | 62.4 |
6843.6602 | 26.0 | 4.55 | −0.83 | 60.9 | 87.2 |
6999.8799 | 26.0 | 4.10 | −1.46 | 53.9 | ⋯ |
7022.9502 | 26.0 | 4.19 | −1.15 | 64.5 | ⋯ |
7132.9902 | 26.0 | 4.08 | −1.65 | 43.1 | ⋯ |
4576.3330 | 26.1 | 2.84 | −2.95 | 64.6 | 111.2 |
4620.5132 | 26.1 | 2.83 | −3.21 | 50.4 | 78.8 |
5234.6240 | 26.1 | 3.22 | −2.18 | 82.9 | 103.6 |
5264.8042 | 26.1 | 3.23 | −3.13 | 46.1 | 63.8 |
5414.0718 | 26.1 | 3.22 | −3.58 | 27.3 | 45.2 |
5425.2568 | 26.1 | 3.20 | −3.22 | 41.9 | 60.5 |
6369.4619 | 26.1 | 2.89 | −4.11 | 19.2 | 35.6 |
6432.6758 | 26.1 | 2.89 | −3.57 | 41.3 | 62.7 |
6516.0771 | 26.1 | 2.89 | −3.31 | 54.7 | 51.2 |
Note. For species we are using MOOG standard notation for atomic or molecular identification. For example, 26 represents Fe(26), while the 0 after the decimal indicates neutral and 1 indicates singly ionized. χexc is the line excitation potential in electron volts (eV).
Download table as: ASCIITypeset image
Footnotes
- *
KIC 5307747.
- 9