Research Paper

Chemical Composition of RR Lyn – an Eclipsing Binary System with Am and λ Boo Type Components

Yeuncheol Jeong1, Alexander V. Yushchenko2,http://orcid.org/0000-0002-9325-5840, Dmytry N. Doikov3, Vira F. Gopka4, Volodymyr O. Yushchenko4
Author Information & Copyright
1Faculty of General Education, Sejong University, Seoul 05006, Korea
2Department of Astronomy and Space Science, Sejong University, Seoul 05006, Korea
3Department of Mathematics, Physics and Astronomy, Odessa National Maritime University, Odessa 65001, Ukraine
4Astronomical Observatory, Odessa National University, Odessa 65014, Ukraine
Corresponding Author Tel: +82-10-4093-5325, E-mail: avyushchenko@gmail.com

© The Korean Space Science Society. All rights reserved. This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Apr 14, 2017; Revised: May 15, 2017; Accepted: May 18, 2017

Abstract

High-resolution spectroscopic observations of the eclipsing binary system RR Lyn were made using the 1.8 m telescope at the Bohuynsan Optical Astronomical Observatory in Korea. The spectral resolving power was R = 82,000, with a signal to noise ratio of S/N > 150. We found the effective temperatures and surface gravities of the primary and secondary components to be equal to Teff = 7,920 & 7,210 K and log(g) = 3.80 & 4.16, respectively. The abundances of 34 and 17 different chemical elements were found in the atmospheric components. Correlations between the derived abundances with condensation temperatures and the second ionization potentials of these elements are discussed. The primary component is a typical metallic line star with the abundances of light and iron group elements close to solar values, while elements with atomic numbers Z > 30 are overabundant by 0.5–1.5 dex with respect to solar values. The secondary component is a λ Boo type star. In this type of stars, CNO abundances are close to solar values, while the abundance pattern shows a negative correlation with condensation temperatures.

Keywords: stars; eclipsing binaries (RR Lyn); individual (RR Lyn); abundances; accretion; chemically peculiar

1. INTRODUCTION

The majority of stars are members of binary and multiple systems; however, most investigations of stellar chemical composition do not account for stellar binarity. This is due to the assumption that the components of unevolved binary systems should have similar chemical abundance patterns. It is important to note that investigations of the chemical composition of binary components face a variety of problems that include the necessity to obtain the highest quality spectral observations and extracting the component’s spectra from the composite spectrum of the overall binary system. To understand the atmospheric chemical abundances of binary components, we need to take into account not only the standard stellar evolution theory, but also the effects that the accretion processes of interstellar and circumstellar gases, mainly hydrogen and helium, have on the stellar atmospheres.

It is proposed that accretion from the interstellar and circumstellar environment should strongly influence ob-servations of surface abundances. For example, Fowler et al. (1965) and Proffitt & Michaud (1989) pointed out the importance of mass transfer from the evolved binary companion, and Drobyshevskii (1975) and Cowley (1977) noted the importance of accretion for asteroids and planets. However, the effect of accretion processes on the surfaces of normal stars still requires further investigation. Additionally, other types of accretion processes may exist that are not yet proved or discovered. Yushchenko et al. (2015) looked for evidence that the accretion of interstellar hydrogen and helium had an effect on the atmospheres of ρ Pup (δ Sct type star) and several hundred other B-F type stars. With these effects taken into account, both Yushchenko et al. (1999) and Kang et al. (2013) presented examples of binary systems with very different abundance patterns for the components.

Herein, we present the atmospheric chemical abundances of components for the well-known eclipsing binary system, Lyn. The orbital period of this system is close to 10 days. The most recently published solutions of the radial velocity curve for RR Lyn were published by Bensch et al. (2014) and Tomkin & Fekel (2006) and the light curve solution can be found in Khaliullin et al. (2001). The masses of the binary components are determined to be 1.927 ± 0.008 and 1.507 ± 0.004 solar masses, and the radiuses are equal to 2.57 ± 0.02 and 1.59 ± 0.03 solar radiuses (Tomkin & Fekel 2006). The age of the system was estimated as 1.08 ± 0.15 billion years, and its evolutionary status is classed as before the first mass exchange (Khaliullin et al. 2001).

The chemical compositions of the components of RR Lyn have been previously investigated by Lyubimkov & Rachkovskaya (1995), and it was found that the abundance patterns of the components are different from each other based on 15 and 13 chemical elements present in the atmospheres of the primary and the secondary components, respectively.

Yushchenko et al. (2015) pointed out that this system is one of the five binary systems with a λ Boo type abundance pattern in at least one of the components, based on the criteria established by Lyubimkov & Rachkovskaya (1995). It should be noted that the total number of identified λ Boo type stars is on the order of one hundred, with no clear understanding yet of their nature. Nevertheless, for these stars, Venn & Lambert (1990, 2008) proposed a dust-gas separation mechanism, and Andrievsky (1997) regarded these stars as the descendants of close binary systems. Yushchenko et al. (2015) also suggested increasing the high temperature limit for this group to 14,000 K. Therefore, a new investigation into the abundance patterns of RR Lyn

components can be important for understanding the physical nature of λ Boo type stars. It is expected that the accretion phenomena in a binary system, such as the well-known mass exchange between components, could be responsible for the λ Boo type abundance pattern.

In this study, we aim to refine the chemical composition of the binary component using high quality spectral observations, modern software, and atomic line databases. We present the results of a spectrum synthesis fit of the observed spectrum as well as a comparison between the observed abundances of chemical elements with their second ionization potentials and the condensation temperatures of these elements.

2. SPECTROSCOPIC OBSERVATIONS AND ABUNDANCE ANALYSIS

The high-resolution spectrum of RR Lyn was obtained from the Bohyunsan Echelle Spectrograph (BOES) installed on the 1.8 m telescope at the Bohyunsan National Optical Astronomical Observatory in Korea. The signal to noise ratio (S/N) of the spectrum is higher than 150 at the center of echelle orders, the resolving power was equal to R = 82,000, and the wavelength coverage is from 3,780 Å to 9,500 Å. The observation was obtained on the Julian date 2,453,369.16 and the exposure time was 30 minutes. The extraction of the observed spectrum from the CCD image and the corresponding wavelength calibration was made using the standard IRAF software. Further data processing for the placement of continuum and the identification of spectral lines was performed using the URAN code (Yushchenko 1998).

The separation of the spectral lines for the primary and the secondary components at the moment of observation was near 135 km/s, which allowed both of the component lines to be investigated. To determine the continuum level and identify the various lines, the latest version of Kurucz (1993) SYNTHE code was used to generate the synthetic spectra of both components for the entire observed wavelength region. Line data were taken from the Kurucz (1993) database, Hirata Horaguchi (1995) line list, the VALD database (Piskunov et al. 1995), Morton (2000), Fuhr & Wises (2006), the DREAM database (Biémont et al. 2017), and other sources. Spectral lines of the appropriate quality from the observed spectrum were selected on the basis of comparison with the combined synthetic spectrum of the binary system.

To construct the combined spectrum, the synthetic spectra of both components were shifted to the observed radial velocities, broadened, and co-added. Fig. 1 shows an example of the fit obtained from the observed spectrum combined with the synthetic spectrum of the binary system. More details about the calculation procedure can be found in Yushchenko et al. (2015) and Kang et al. (2012, 2013). It is necessary to note that analyzing composite spectra is often challenging. Systematic errors can be minimized when the difference between the component spectra can be identified. In this study, we attempt to confirm or reject the large observed differences in chemical abundances between the RR Lyn binary components that were first noted by Lyubimkov & Rachkovskaya (1995).

jass-34-2-75_F1
Fig. 1. The RR Lyn example spectrum. The axes are the wavelength in angstroms and the relative flux. The squares in the upper part of the figure denote the observed spectrum; the thin line denotes the combined synthetic spectrum of the binary system. Scaled and shifted synthetic spectra of the components are shown in the middle of the figure. The thin and thick lines are the spectra of the primary and secondary components, respectively. The synthetic spectra are shifted in accordance with the observed radial velocities and smoothed by the instrumental and rotational profiles. The positions of the spectral lines, which were taken into account in calculations, are marked in the lower part of the figure with short and long dashes (faint and strong lines, respectively) for the primary (upper line of marks) and secondary (lower line of marks) components. The identifications are shown for the strongest lines.
Download Original Figure

To determine the effective temperatures, surface gravities, and other parameters of the components, the method de-veloped by Yushchenko et al. (1999) and Gopka et al. (2004) was used. This method uses the iron abundances calculated for spectral lines of iron using the different atmospheric models to determine the stellar parameters. The grids of the atmospheric models with effective temperatures from 7,500 K to 8,100 K (step 10 K), surface gravities log(g) = 3.60 to 4.40 (step 0.02) for the primary component, and from 7,000 K to 7,400 K (step 10 K), surface gravities log(g) = 3.70 - 4.50 (step 0.02) for the secondary component were prepared by interpolation of the Kurucz (1993) atmospheric models. A detailed description of the method can be found in Kang et al. (2012, 2013).

The flux ratio of the binary components at wavelength λ = 5,500 Å was derived to be 3.0. For the primary and secondary components, we derived the following parameters: the effective temperatures Teff = 7,920 ± 50 and 7,210 ± 50 K, surface gravities log(g) = 3.80 ± 0.10 and 4.16 ± 0.10, microturbulent velocities vmicro = 4.60 ± 0.30 and 3.15 ± 0.30 km/s, and projected rotational velocities v sin i = 11.7 ± 1.5 and 6.1 ± 1.5 km/s. The iron abundances [Fe/H] were +0.20 ± 0.10 and -0.34 ± 0.10. The errors indicated are systematic errors associated with the method. Comparison of our values with Lyubimkov & Rachkovskaya (1995), Khaliullin et al. (2001), and Tomkin & Fekel (2006) show that higher systematic errors have been reported. Table 1 shows the parameters of the RR Lyn components published by the above authors alongside our values.

Table 1. Parameters of the primary and secondary components from this study with values obtained from the literature for comparison
jass-34-2-75_T1
Download Excel Table

The set of parameters shown in Table 1 were used for determining the chemical composition of both binary com-ponents. The resulting abundance calculations were then made using the spectrum synthesis method. The latest versions of Kurucz (1993) SYNTHE and Yushchenko (1998) URAN codes were used to make the fit between the observed spectrum and the combined synthetic spectrum of the binary system in semiautomatic mode.

To avoid possible errors in the oscillator strengths, the counterparts of the investigated absorption lines in the solar spectrum were found. The Liege atlas of the solar spectrum (Dellbouille et al. 1973) and the solar photosphere model (Grevesse & Sauval, 1999) were used for this purpose. It should be noted that the use of a comparison star decreases the uncertainties associated with unknown hyperfine structure and isotope shifts for spectral lines. LTE approximation calculations of the synthetic spectra were made.

Table 2 summarizes the abundance calculations. The first two columns of the table are the atomic number, the chemical symbol, and the ionization state for each of the chemical elements. In columns three to six, the number of spectral lines used and the mean atmospheric chemical abundances for both components with respect to the solar photosphere abundances can be found. In brackets after the corresponding relative abundances are the calculated error values. The last four columns of the table show the Lyubimkov & Rachkovskaya (1995) results for comparison. The numbers of used lines and relative abundances for both components are also shown.

Table 2. Chemical composition of RR Lyn (with respect to the Sun)
jass-34-2-75_T2
Download Excel Table

Fig. 2 shows the abundance patterns for both components. The abundances of 34 and 17 different chemical elements were found in the atmospheres of the primary and secondary components, respectively.

jass-34-2-75_F2
Fig. 2. Chemical composition of RR Lyn. The axes are the atomic numbers of chemical elements and the logarithmic abundances with respect to solar values. The horizontal dashed line marks the solar abundances, open circles represent the abundances in the photosphere of the primary companion, and the filled circles represent the secondary companion.
Download Original Figure

3. DISCUSSION

Table 2 and Fig. 2 clearly show that the abundance patterns for the components are significantly different from each other. This confirms the main result of Lyubimkov & Rachkovskaya (1995), which also indicated that there were large differences between the component abundances.

The primary component exhibits a relative underabundance for the CNO elements, and the solar or slightly oversolar abundances for elements with atomic numbers from Z = 11 to Z = 29. The exceptions are calcium and scandium, which are underabundant. The heavier elements, elements with atomic numbers Z > 29, are overabundant in the atmosphere of the primary. These abundance pattern features allow us to conclude that the primary component of RR Lyn is a typical metallic line (Am) star.

In contrast, for the atmosphere of the secondary component, CNO elements are close to solar abundance values. The abun-dances of other elements in the secondary are slightly under-abundant or near to solar values. This type of abundance pattern is typical for λ Boo type stars.

A possible reason for the differences in atmospheric chemical composition for both components of this binary system is that some external influences, for example the accretion of gas and dust from the circumstellar envelope, played a role in modifying what were initially similar atmospheres. The possible types of accretions that could be responsible were briefly pointed out in this manuscript and were also discussed by Yushchenko et al. (2015) in more detail.

If accretion processes occur in the RR Lyn system, it is possible to expect the dependences between the relative chemical abundances, the second ionization potentials and the condensation temperatures of the corresponding elements. Fig. 3 shows the first of these dependencies. Comparison of this figure with figure 6 in Yushchenko et al. (2015) allows us to conclude that the abundance patterns of ρ Pup and the primary component of the RR Lyn binary system have similar properties. This means that the accretion of interstellar and circumstellar gas, mainly in the form of hydrogen and helium, can change the atmospheric chemical composition of the primary component of RR Lyn as well as ρ Pup.

jass-34-2-75_F3
Fig. 3. The relative atmospheric abundances of chemical elements (represented by atomic numbers) for the primary (A) and secondary (B) components of RR Lyn, as a function of second ionization potentials of these elements. The ionization energies of hydrogen and helium, 13.6 and 24.6 eV, respectively, are marked by vertical dashed lines.
Download Original Figure

The main features of this accretion scenario were first discussed by Greenstein (1949) and later developed by Havnes (1971). It was proposed that the charge exchange reactions between hydrogen atoms from the interstellar medium with atoms from the stellar photosphere results in a decrease of the relative abundances of elements with second ionization potentials close to the ionization potential of hydrogen, namely those close to 13.6 eV. The result of possible resonance in the charge exchange reactions is that the escape velocities of atoms with the above highlighted second ionization potentials, and perhaps the production of cosmic rays, would be able to brake the stellar rotational velocities. Böhm-Vitense (2006) found that a similar effect could also be important for helium atoms.

The abundance pattern of the secondary component of the Lyn binary system includes only 17 chemical elements, but the increase in the relative abundance of chemical elements with second ionization potentials between 13.6 and 24.6 eV (the ionization potentials of hydrogen and helium) are clearly detected in the bottom panel of Fig. 3. It is therefore possible that the accretion of hydrogen and helium from outer space can be an important process for both components of the RR Lyn binary system.

Fig. 4 shows a plot of the relative chemical abundances in the atmospheres of RR Lyn components against the condensation temperatures of the corresponding elements. The negative correlation illustrated by this plot is typical for λ Boo type stars (Venn & Lambert 1990, 2008). The correlation can be explained by a dust-gas separation mechanism resulting from accretion in a circumstellar cloud. Elements with higher condensation temperatures condense into dust particles in regions with higher temperatures, which are closer to the stellar surface. After condensation, the dust particles are influenced by solar radiation pressure and thus are effectively removed from the circumstellar environment. This explanation is one of the most plausible to justify the negative correlation observed between the relative abundances and condensation temperatures in the atmospheres of λ Boo type stars.

jass-34-2-75_F4
Figure 4. The relative surface abundances of chemical elements (represented by atomic numbers) in the atmospheres of primary (A) and secondary (B) components of RR Lyn, as a function of condensation temperatures for these elements. The values of condensation temperatures for a solar-photosphere composition gas are taken from Lodders (2003). The explanation of elliptic area in the upper panel is given in the text.
Download Original Figure

Fig. 4 shows that a negative trend is evident for the data presented in the bottom panel of the figure. The correlation coefficient is equal to -0.73 +/- 0.17. This implies that the dust-gas separation mechanism proposed by Venn & Lambert (1990, 2008) can explain the observed distribution for the relative abundances in the atmosphere of the secondary component.

Further discussion is needed for the main component shown in the upper panel of Fig. 4. Elements with condensation temperatures higher than 1,000 K and elements with relative abundances higher than that of strontium (Z = 38) are not usually detected in the spectra of λ Boo type stars. For CNO elements and scandium (second ionization potential = 12.8 eV), the underabundances of these elements can be explained by charge-exchange reactions with hydrogen.

The remaining elements, namely the 15 elements located in the highlighted elliptical area of the upper panel of Fig. 4, also show a strong negative correlation. The correlation coef-ficient is -0.57 +/- 0.17. This evidence gives strength to the conclusion that the dust-gas separation scenario could be the most important influencing process for these elements while the other elements are strongly influenced by other physical effects.

Figs. 3 and 4 allow us to surmise that the different accretion scenarios presented do take place in the RR Lyn binary system. It is worth taking into consideration that the abundance patterns of both components can be the result of a mass transfer process between the two components. The influence of mass transfer on the distribution of the relative abundance in the binary components was not investigated for heavy elements earler.

Finally, we confirm the large chemical composition differences observed for the RR Lyn binary components first published by Lyubimkov & Rachkovskaya (1995). We assert that the atmospheres of these components are strongly influenced by different types of accretion. These are the accretion of hydrogen and helium from interstellar or circumstellar gas (Greenstein 1949;Havnes 1971;Böhm-Vitense 2006), and dust-gas separation (Venn & Lambert 1990, 2008). It is necessary to point out that the net metallicity difference between the two components is significant and that further investigations may be needed to explain this difference in relation to accretion effects. Additionally, the diffusion effects first proposed by Michaud (1970) should be also taken into account during further studies.

More careful analyses of RR Lyn and similar binary systems can help us to understand the physical processes that occur in stellar atmospheres. The theory of accretion effects due to hydrogen and helium from the interstellar or circumstellar environment on the atmospheres of normal stars was proposed by Greenstein (1949) and Havnes (1971), but this scenario had very few observational confirmations until now with the exception of Böhm-Vitense (2006) and Yushchenko et al. (2015). The detailed theory of this effect has not yet been fully developed, and would require new observational results and further research.

ACKNOWLEDGMENTS

This work was supported in part by the Korea Astronomy and Space Science Institute under the R&D program (Project No. 2015-1-320-18) supervised by the Ministry of Science, ICT, and Future Planning. One of the authors (VYu) was given a grant by the Swiss National Science Foundation (SCOPES project No.~IZ73Z0152485).

References

1.

Andrievsky SM, On the possible origin of λ Boo stars, Astron. Astrophys. 321, 838-840 (1997).

2.

Bensch K, Dimitrov W, Żywucka N, Fagas M, Kamiński K, et al., New spectroscopy of multiple stars RR Lyncis and HT Virginis, Inf. Bull. Var. Stars 6121, 1-7 (2014).

3.

Biémont E, Palmeri P, Quinet P, Database on Rare Earths At Mons University (DREAM) [Internet], cited 2017 Apr 14, available from: http://hosting.umons.ac.be/html/agif/databases/dream.html

4.

Böhm-Vitense E, The puzzle of the metallic line stars, Publ. Astron. Soc. Pac. 118, 419-435 (2006).

5.

Cowley CR, An examination of the planetesimal impact hypothesis of the formation of CP stars, Astrophys. Space Sci. 51, 349-362 (1977).

6.

Dellbouille L, Rolland G, Neven L, Photometric Atlas of the Solar Spectrum from λ 3000 to λ 10000 (University of Liège, Liège, 1973).

7.

Drobyshevski EM, Peculiar A-stars and planetary systems, Astrophys. Space Sci. 35, 403-408 (1975).

8.

Fowler W, Burbidge EM, Burbidge GR, Hoyle F, The synthesis and destruction of elements in peculiar stars of types a and B, Astrophys. J. 142, 423-450 (1965).

9.

Fuhr JR, Wiese WL, A critical compilation of atomic transition probabilities for neutral and singly ionized iron, J. Phys. Chem. Ref. Data 35, 1669-1809 (2006).

10.

Gopka VF, Yushchenko AV, Mishenina TV, Kim C, Musaev FA, et al., Atmospheric chemical composition of the halo star HD 221170 from a synthetic-spectrum analysis, Astron. Rep. 48, 577-587 (2004).

11.

Greenstein JL, Analysis of the metallic-line stars. II, Astrophys. J. 109, 121-138 (1949).

12.

Grevesse N, Sauval AJ, The solar abundance of iron and the photospheric model, Astron. Astrophys. 347, 348-354 (1999).

13.

Havnes O, Abundances and acceleration mechanisms of cosmic rays, Nature 229, 548-549 (1971).

14.

Hirata R, Horaguchi T, VizieR Online Data Catalog: Atomic Spectral Line List, SIMBAD Catalog VI/69 (1995). available from: http://vizier.cfa.harvard.edu/viz-bin/Cat?VI/69

15.

Kang YW, Yushchenko A, Hong K, Kim S, Yushchenko, V, Chemical composition of the components of eclipsing binary star ZZ Bootis, Astron. J. 144, A35 (2012).

16.

Kang YW, Yushchenko AV, Hong K, Guinan EF, Gopka VF, Signs of accretion in the abundance patterns of the components of the RS CVn-type eclipsing binary star LX Persei, Astron. J. 145, A167 (2013).

17.

Khaliullin KF, Khaliullina AI, Krylov AV, Precision WBVR photoelectric photometry of the eclipsing system RR Lyncis, Astron. Rep. 45, 888-898 (2001).

18.

Kurucz RL, SYNTHE spectrum synthesis programs and line data, in Kurucz CD-ROM 18 Series (Smithsonian Astrophysical Observatory, Cambridge, 1993).

19.

Lodders K, Solar system abundances and condensation tem-pertures of the elements, Astrophys. J. 591, 1220-1247 (2003).

20.

Lyubimkov LS, Rachkovskaya TM, The Am binary star RR Lyn: chemical composition of the components, Astron. Rep. 39, 56-62 (1995).

21.

Michaud G, Diffusion processes in peculiar A stars, Astrophys. J. 160, 641-658 (1970).

22.

Morton DC, Atomic data for resonance absorption lines. II. Wavelengths longward of the Lyman limit for heavy elements, Astrophys. J. Suppl. Ser. 130, 403-436 (2000).

23.

Piskunov NE, Kupka F, Ryabchikova TA, Weiss WW, Jeffery CS, VALD: the Vienna atomic line data base, Astron. Astrophys. Suppl. Ser. 112, 525-535 (1995).

24.

Proffitt CR, Michaud G, Abundance anomalies in A and B stars and the accretion of nuclear-processed material from supernovae and evolved giants. Astrophys. J. 345, 998-1007 (1989).

25.

Tomkin J, Fekel FC, New precision orbits of bright double-lined spectroscopic binaries. I. RR Lyncis, 12 Bootis, and HR 6169, Astron. J. 131, 2652-2663 (2006).

26.

Venn KA, Lambert DL, The chemical composition of three Lambda Bootis stars, Astrophys. J. 363, 234-244 (1990).

27.

Venn KA, Lambert DL, Could the ultra-metal-poor stars be chemically peculiar and not related to the First stars, Astrophys. J. 677, 572-580 (2008).

28.

Yushchenko AV, URAN: a software system for the analysis of stellar spectra, Proceedings of the 20th Stellar Conference of the Czech and Slovak Astronomical Institutes, ed. Dusek J (the Czech and Slovak Astonomical Institutes, Brno, 1998), 201-203.

29.

Yushchenko AV, Gopka VF, Khokhlova VL, Musaev FA, Bikmaev IF, Atmospheric chemical composition of the "twin" components of equal mass in the CP SB2 system 66 Eri, Astron. Lett. 25, 453-466 (1999).

30.

Yushchenko AV, Gopka VF, Kang YW, Kim C, Lee BC, et al., The chemical composition of ρ Puppis and the signs of accretion in the atmospheres of B-F-type stars, Astron. J. 149, A59 (2015).