Research Paper

The Chemical Composition of HD47536: A Planetary Host Halo Giant with Possible λ Bootis Features and Signs of Interstellar Matter Accretion

Alexander Yushchenko1,2, Dmytry Doikov3, Sergei Andrievsky2,4,5, Yeuncheol Jeong6,, Volodymyr Yushchenko2,7, Pakakaew Rittipruk8, Valery Kovtyukh2,5, Aizat Demessinova9, Vira Gopka2, Alexander Raikov10,11, Kyung Sook Jeong1,2
Author Information & Copyright
1Astrocamp Contents Research Institute, Goyang 10329, Korea
2Astronomical Observatory, Odessa National University, Odessa 65014, Ukraine
3Department of Natural and Technical Sciences, Odessa National Maritime University, Odessa 65029, Ukraine
4Observatoire de Paris-Meudon, CNRS, Universite Paris Diderot, Meudon Cedex 92125, France
5Institut für Astronomie und Astrophysik, Kepler Center for Astro and Particle Physics, Universität Tübingen, Tübingen D-72076, Germany
6Department of History, Sejong University, Seoul 05006, Korea
7Main Astronomical Observatory of National Academy of Sciences of Ukraine, Kyiv 03143, Ukraine
8National Astronomical Research Institute of Thailand, Chiangmai 50180, Thailand
9Physico-Technical Department, Al Farabi Kazakh National University, Almaty 050040, Kazakhstan
10Pulkovo Observatory, Russian Academy of Sciences, Saint Petersburg 196140, Russia
11Special Astrophysical Observatory, Russian Academy of Sciences, Nizhny Arkhyz 369167, Russia
12Institute of Liberal Education, Incheon National University, Incheon 22012, Korea
Corresponding Author Tel: +82-2-3408-3134, E-mail:

© 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 ( which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Nov 15, 2022; Revised: Dec 13, 2022; Accepted: Dec 13, 2022


We investigated the chemical composition of the planetary host halo star HD47536 via high-resolution spectral observations recorded using a 1.5 meter Cerro Tololo Inter-American Observatory (CTIO) telescope (Chile). Furthermore, we determined the abundances of 38 chemical elements. Both light and heavy elements were overabundant compared to the iron group elements. The abundance pattern of HD47536 was similar to that of halo-type stars, with an enrichment of heavy elements. We analyzed the relationships between the relative abundances of chemical elements and their second ionization potentials and condensation temperatures. We demonstrated that the interplay of charge-exchange reactions owing to the accretion of interstellar matter and the gas-dust separation mechanism can influence the initial abundances and can be used to qualitatively explain the abundance patterns in the atmosphere of HD47536.

Keywords: stars: individual (HD47536); stars: abundances; stars: circumstellar matter; stars: atmospheres; stars: population II; physical data and processes: accretion


The present study investigated the chemical composition of HD47536 in detail and analyzed specific chemical abundance features observed in similar stars. HD47536 is a bright (V = 5.25) halo or an intermediate-population K-type giant. Different investigators have determined that the iron abundance in HD47536 ranges from [Fe/H] = –0.54 ± 0.10 (Sadakane et al. 2005) to [Fe/H] = –0.68 (da Silva et al. 2006).

It hosts one or two planets, of which the first was discovered by Setiawan et al. (2003). The mass of this planet can be more than 4.0 ± 0.4 times that of Jupiter, the semimajor axis of its orbit is 1.12 ± 0.005 astronomical units, and the eccentricity of its orbit is 0.3 ± 0.1 (Soto et al. 2015). Notably, the radius of the star is 21.36 ± 1.47 R (van Belle & von Braun 2009). This means that the shortest distance between the planet and the host star can be as small as 0.78 astronomical units, or near eight radii of the host star. Soto et al. (2015) examined the existence of a second planet.

Importantly, only 42 halo or thick disk stars are known as planetary hosts (Perottoni et al. 2021), and HD47536 is one of the brightest.

Numerous studies have compared the chemical abundances of stars with and without planets. The survey of 1111 FGK stars by Adibekyan et al. (2012) revealed that for stars with giant planets up to 0.5 dex, the abundances of N, Mg, Al, Si, Ca, and iron peak elements were higher than for stars without planets or those with Neptunian and super-Earth planets.

The chemical composition of HD47536 has been previously analyzed. Metallicity (Sadakane et al. 2005) and the abundance of a few chemical elements, specifically the upper limits of lithium and beryllium, were determined by Gálvez-Ortiz et al. (2011). The abundance of oxygen was determined by Ecuvillon et al. (2006); sodium, magnesium, and aluminum by Beirão et al. (2005); sulfur by Caffau et al. (2005); and Si, Ca, Sc, Ti, V, Cr, Mn, Co, and Ni by Gilli et al. (2006). The star was also included in the Luck (2015) Survey.

The preliminary abundance distribution in HD47536 was determined by Yushchenko et al. (2013), who derived the abundances of 33 chemical elements using the model atmosphere method. These results were used as the basis for calculating the individual atmospheric model using the ATLAS12 code.

The second and third sections of this paper describe the observations, data reduction, and model atmospheric parameters of HD47536. In the fourth section, we report the abundances of 38 chemical elements in the atmosphere of HD47536 derived using the differential spectrum synthesis method.

In the last section of the paper, we discuss the obtained results and outline possible ways to develop a theoretical explanation regarding the observed correlations of the abundances of chemical elements in the atmosphere of HD47536 with condensation temperatures and second ionization potentials of the elements.


Four spectra of HD47536 were obtained using the Chiron spectrograph installed on a 1.5 meter SMART Cerro Tololo Inter-American Observator (CTIO) telescope (Chile) in 2012. These spectra were observed at heliocentric Julian dates 2,444,897.834, 2,444,897.835, 2,444,900.781, and 2,444,900.746 with exposure times of 50, 50, 300, and 600 s, respectively. The 3-pixel sampling spectral resolving power was R = 30,000, the wavelength range was 4,105–8,170 angstroms, and the signal-to-noise ratio S/N at the centers of echelle orders in the red part of the wavelength region exceeded 100 and 250, respectively, for the first two and the last two spectra.

The initial reduction of spectra was performed using the standard Image Reduction and Analysis Facility (IRAF) package, whereas the final processing of spectra, including continuum placement, co-adding of the spectra, identification of spectral lines, and equivalent width measurements, was performed using the latest version of the URAN code (Yushchenko 1997).

We found no significant differences between the observed spectra; therefore, all measurements were performed using the average observed spectrum of the star. The four observed spectra were shifted in accordance with radial velocity corrections and were co-added with weights proportional to exposure times.

An example of the spectrum in the vicinity of the thorium line 5,989.045 Å is shown in Fig. 1, which shows the averaged profile over the four observed spectra compared to the synthetic spectrum. The synthetic spectrum of HD47536 was calculated for the entire observed wavelength region using the Kurucz (1993) SYNTHE code. This was considered when placing the continuum and identifying spectral lines.

Fig. 1. The coadded observed spectrum of HD47536 (points) and its approximation by synthetic spectra (solid lines). The axes show the wavelength in angstroms and the relative fluxes. The positions of spectral lines taken into account in the calculations are marked in the bottom part of the figure. Identifications are given for some of the strongest lines. The position of Th II λ 5,989.045 Å line is marked by a vertical dotted line. Note that three synthetic spectra are shown in the vicinity of thorium line. These synthetic spectra are calculated for the best thorium abundance derived from this line (log N(Th) = +0.02), and the abundances which are lower and higher by 0.5 dex with respect to the best value (log N(Th) = –0.48 and +0.52 respectively). Two boxes in the central part of the plot compare the coadded observed spectrum (upper box) and three individual observed spectra (lower box).
Download Original Figure


The atmospheric parameters derived by Yushchenko et al. (2013) (Teff = 4,400K, log g = 1.80, vmicro = 1.5 km/s) and the abundances of chemical elements were used as input information for calculating the individual atmosphere model with the ATLAS12 code (Kurucz 1993). To test the accuracy of this model, we calculated the correlation coefficients of iron abundance determined using the Kurucz’s (1993) WIDTH9 code for 90 lines of neutral iron as a function of log (Eq. width/λ) and the excitation potential of the lower level for these lines. They were determined to be equal to 0.00 ± 0.11 and 0.07 ± 0.11, respectively. The value of the microturbulent velocity used for these calculations was slightly changed: vmicro = 1.44 km/s.

The abundance of iron derived from 90 lines of neutral iron and 18 lines of ionized iron were determined to be equal to log N(Fe) = 6.99 ± 0.11 and log N(Fe) = 7.01 ± 0.15, respectively. The oscillator strengths for lines of neutral and ionized iron were taken primarily from Fuhr & Wiese (2006) –73 lines, for 25 lines Kurucz (1993) values were used, and for 10 lines, we calculated the solar oscillator strength.

The above correlations and abundances of iron derived from neutral and ionized iron lines implied that the ATLAS12 model could be used to determine the abundances of all chemical elements using the spectrum synthesis method. To estimate the influence of possible uncertainties in atmospheric parameters on the derived abundances, we calculated two additional atmospheric models with slightly shifted parameters, specifically the effective temperature increased by 100K and the surface gravity increased by 0.4 dex. All calculations were made for these three atmosphere models as well as for the first model with the microturbulent velocity increased by 0.1 km/s. The spectral synthesis method was extremely sensitive to the observed spectral line profiles. Therefore, special attention was paid while determining the broadening parameters. First, we found that the spectral resolving power at the time of observation was higher than that described in the previous section. Based on terrestrial line profiles, we estimated the spectral resolution to be R = 48,000. Subsequently, the additional line broadening in the spectrum of HD47536 was approximated by the macroturbulent velocity vmacro = 7 km/s. The broadening of the observed spectrum by the rotational profile was neglected.

Notably, the difference between the standard resolving power of the spectrograph and the value obtained from our observed spectra can be explained by better assessment at the moments of our observations, primarily by different methods of deriving the resolving power in our investigation. The value of 3-pixel sampling spectral resolution is strongly dependent on the charge-coupled device used; however, the value of the instrumental profile of the spectral optical system used can be different.


The atmospheric parameters described in the previous section were used to calculate the synthetic spectrum of HD47536 over the entire observation region. A comparison of this spectrum with the observed spectra of HD47536 made it possible to reliably identify the spectral lines and select the best lines for determining the chemical composition. The abundances of all chemical elements (except iron) in the HD47536 atmosphere were derived using the spectrum synthesis method. We used Kurucz’s (1993) SYNTHE code and our URAN (Yushchenko 1997) software to perform calculations in semiautomatic mode.

We managed to find the counterparts of most lines in the solar spectrum. Furthermore, the abundances of the solar lines were derived using the spectrum synthesis method. We used the Liége solar atlas (Delbouille et al. 1973) and the Grevesse & Sauval (1999) solar atmosphere models. The adopted values of microturbulent and macroturbulent velocities for this model were 0.8 km/s and 1.8 km/s, respectively. The continuum level in the solar spectrum was adjusted according to Ardeberg & Virderfors (1979) and Rutten & van der Zalm (1984). The same codes were used to build a synthetic spectrum, identify lines, and determine the abundances.

Table 1 presents the average relative and absolute values of abundances in the atmosphere of HD47536 for 38 chemical elements. The first two columns indicate atomic numbers and ion identification. The third column lists the number of lines used. The two subsequent groups of columns present the mean abundances in the atmosphere of HD47536 and their errors for the four sets of atmospheric parameters, specifically for the best values and for effective temperature, surface gravity, and microturbulent velocity increased by 100K, 0.4 dex, and 0.1 km/s, respectively. The first and second groups of columns depict relative and absolute abundances, respectively. The relative abundances (*–⊙) and the absolute abundances (log N) in the atmosphere of HD47536 are the first columns of each group.

Table 1. Abundances of chemical elements in the atmosphere of HD47536. Mean values (Continued on the next page)
Download Excel Table

Fig. 2 shows the relative abundances of chemical elements in the atmosphere of HD47536 as a function of the atomic number.

Fig. 2. Chemical composition of HD47536. The axes show atomic number and logarithmic abundance with respect to solar values. Small, medium, and large circles in this and the next two figures designate the abundances calculated using 1–2, 3–4, and 5 or more lines, respectively. The lithium relative abundance (Δlog N(Li) = –1.68) is much lower than the minimum value of Fig. 2, the direction to this point is shown in the lower left corner of the plot.
Download Original Figure


From Table 1 and Fig. 2, we can conclude that there is a well-expressed regularity in the [El/Fe] distribution. The light elements (from carbon to silicon) appear to be slightly overabundant, the iron-peak elements are deficient (being nearly at the same level of iron abundance), and the scenario for heavier species is rather controversial. Some elements are more abundant than iron (e.g., Rb, Mo, Ru, Ba, Sm, Eu, and others), while others are equally abundant (e.g., Y, Zr, and Ce).

In Fig. 3, we depict the abundances in HD47536 superimposed on the scaled abundances of the solar system r-process elements. The deviations from the distribution of the solar r-process elements were fitted by the s-processes for Ba, Nd, and heavier lanthanides. The abundances of light lanthanides, La (Z = 57), Ce (58), and Pr (59), decreased with respect to the s-process enrichment of the solar system.

Fig. 3. The upper panel shows the comparison of the surface abundances in HD47536 (circles) with the solar system r-process abundance distribution (Simmerer et al. 2004) scaled at the observed Eu abundance (line). The bottom panel shows the deviations of the observed abundances in HD47536 from the scaled solar system r-abundances (circles). The line is the deviations of solar photosphere abundances from solar r-process abundance distribution. The maximums of this curve are expected for the elements with the highest relative s-process contributions.
Download Original Figure

A possible reason for this decrease can be attributed to charge-exchange reactions. We describe this in more detail in this section. The abundances of elements with atomic numbers of 37 ≤ Z ≤ 42 can be explained by the s-process enrichment only for Sr (38), Zr (40), and Mo (42). The abundances of Rb (37), Y (39), and Nb (41) were not strongly influenced by the s-process.

We investigated whether it is possible to explain the different behaviors of the two groups of chemical elements with odd and even atomic numbers and whether it is sufficient to have only information regarding nuclear processes from previous stellar generations. To make the situation more understandable, it is necessary to examine the relationship between the abundance of the studied chemical elements and some of their physical properties. Among these are properties such as the second ionization potential (SIP) and condensation temperature.

These two characteristics may be important when considering the possible observational manifestation of such chemodynamical processes, such as interaction between atoms in the stellar atmosphere and interaction of the stellar atmosphere with circumstellar or interstellar media. Figs. 4 and 5 show the derived elemental abundances in our program stars depending on the SIP and condensation temperature of the studied chemical elements. The relationships obtained are discussed below.

Fig. 4. The left panel is the plot of relative abundances of chemical elements in HD47536 as a function of the SIP of these elements. Vertical dashed lines mark the positions of ionizations energies of hydrogen and helium. The values of lithium’s SIP (75.6 eV) and relative abundance (Δlog N(Li) = –1.68) define the position right and down from the upper panel. The direction to lithium point is shown in the lower right corner of this panel. Solid circles show 10 chemical elements with a second ionization potentials less than 20 eV. The correlation coefficient of the relative abundance of these 10 elements as a function of the SIP is –0.17 ± 0.31. The abundances and correlation coefficient calculated for these 10 elements will be used in our next paper, that is why we show it here by other symbols. For 17 chemical elements with second ionization potentials in the range from 12.5 to 20 eV the correlation coefficient of relative abundances against the SIPs is +0.19 ± 0.23. For 9 lanthanides with SIPs in the range from 10.55 to 12.09 eV the similar correlation coefficient is +0.51 ± 0.25. The right panel is the enlarged left part of the left panel. SIP, second ionization potential.
Download Original Figure
Fig. 5. Plot of the relative abundance of chemical elements in HD47536 as a function of their condensation temperatures. The condensation temperature values for the solar photosphere composition gas are taken from Lodders (2003). Filled circles show 10 chemical elements with a second ionization potentials less than 20 eV. Correlation coefficient of the relative abundances of these ten elements as a function of the condensation temperature is +0.25 ± 0.30. The lithium’s relative abundance (Δlog N(Li) = –1.68) place this element out of the plot. The direction to lithium point is shown against its condensation temperature (Tcond = 1,135K). The right panel is the enlarged right part of the left panel.
Download Original Figure
5.1 Second Ionization Potential (SIP)

Greenstein (1949) was the first to reveal that in some metallic-line stars, the abundances of elements such as Na, Mg, Al, Si, Sc, Ti, V, and Zr are reduced. He supposed that these elements can persist in the second ionization stage, and therefore, the corresponding lines are not seen in the visual part of the spectra of A–F stars. This mechanism is possible in a non-equilibrium medium if the SIPs are close to the hydrogen ionization potential.

In such a case, the overionization of hydrogen may cause an excess of the above-mentioned elements in the second ionization stage and the respective reduction of their abundances derived from the lines of neutral and first ionization species. This can be achieved either due to the UV line and continuous emission of hydrogen or due to the charge-exchange reactions between ionized hydrogen and the mentioned atoms in the first ionization stage.

Charge-exchange reactions are possible in the case of resonance between the hydrogen ionization potential and SIP of certain chemical elements. Böhm-Vitense (2006) and Yushchenko et al. (2015) observed the same effect for helium, specifically the decrease of relative abundances for chemical elements with SIPs close to or greater than SIP = 24.6 eV.

The value of helium’s first ionization potential is 24.6 eV. Havnes (1971) and Havnes & Conti (1971) suggested that the chemical peculiarity of magnetic stars can be caused by the selective capture of atoms from the ionized interstellar medium. If such a mechanism is implemented to explain some chemical anomalies in HD47536, then we need to suppose that this star was a mid- to late-type dwarf at the main sequence with a rather strong magnetic field.

The charged atoms with a smaller charge-to-mass ratio easily penetrate the main sequence star magnetosphere and reach the stellar atmosphere, whereas the lighter species are reflected in the zones of locally increased magnetic field. In principle, such a mechanism can also work in the circumstellar envelope, possibly with a higher efficiency, because the expected gas concentration in the circumstellar envelope should be much higher than that in the interstellar medium. SIPs for Sr and Zr are close to the ionization potential of hydrogen, while rare Earth’s have smaller SIPs (as a rule, lower than 12 eV), which may explain their increased abundance in the HD47536 star (Fig. 4). However, it should be noted that Borra & Landstreet (1980) did not succeed in detecting strong magnetic fields in metallic stars. North (1984) also doubts the validity of the findings of Havnes (1971) and Havnes & Conti (1971).

Importantly, Greenstein (1949) discovered the observation phenomenon, specifically the deficiency of relative abundances of chemical elements with SIPs between 12 and 16 eV and proposed an explanation for this phenomenon. His physical scenario did not require a magnetic field. Havnes (1971) and Havnes & Conti (1971) tried to develop Greenstein’s scenario for stars with strong magnetic fields; however, their hypothesis faced major problems and needed improvement.

Böhm-Vitense (2006) comprehensively discussed several possible scenarios regarding the formation of metallic-line stars. She came to a conclusion that only the charge-exchange mechanism is able to explain observed abundance peculiarities, particularly the deficiency of the chemical elements with SIPs that are close to 13.6 eV (hydrogen ionization potential) and 24.6 eV (helium first ionization potential).

5.2 Condensation Temperature

Here, we focus on Fig. 5, which shows the relative abundances of 38 chemical elements as a function of their condensation temperature. We note that elements with condensation temperatures below 1,200–1,400K apparently show progressively decreasing abundance with respect to iron. The higher the condensation temperature, the lower the abundance. The lowest abundances in this group of elements were iron-peak elements (refractory siderophyles: Cr, Mn, Fe, and Ni; see Fig. 5).

After this minimum, elements with a higher condensation temperature showed an increased abundance. Among these elements, there are many whose abundances were derived using one or two lines and cannot be considered very reliable. A similar excess was observed for refractory lithophile elements such as Mg and Si, whose abundances were derived using a statistically significant number of lines. Notably, their abundances were higher than those of the iron peak elements, although their condensation temperatures were almost the same as that of iron. It is also interesting to note that the four elements: Na, K, Zn, and Rb, which are classified as volatile lithophile elements, showed a fairly close abundance in HD47536.

First, we discuss the behavior of elements with condensation temperatures below 1,400K, including refractory siderophile elements. This behavior resembles that observed in λ Bootype stars. The chemical peculiarities of λ Boo-type stars are due to a specific mechanism of contamination of their atmospheres by material from the circumstellar envelope, which has previously experienced dust and gas separation. In a certain temperature regime, the gaseous component around the central star consists of volatile fractions with low condensation temperatures and refractory elements with high condensation temperatures. These refractory elements can form dust particles even in a high-temperature environment.

Gravitational attraction acts with equal efficiency on the gas and dust components; however, the radiative pressure of the central star more efficiently removes dust particles farther from the star. Over time, this process reduces the fraction of refractory elements in the surroundings nearest to the star. While the gas component returns to the stellar atmosphere, some elements enclosed in the dust grains are lost to the stellar atmosphere. Thus, as a result, the stellar atmosphere acquires a specific chemical composition. It is not uncommon that dust grains containing refractory elements may coagulate at a greater distance from the star, forming planetesimals and subsequently, planet-like bodies. Recall that the star in our program HD47536 has at least one planet. Thus, mid- and late-type main sequence stars with planets may be deficient in refractory elements.

This scenario was first qualitatively proposed by Venn & Lambert (1990) and later investigated by Charbonneau (1991, 1993), Stuerenburg (1993), Andrievsky (1997), Andrievsky & Paunzen (2000), Andrievsky (2006), and Venn & Lambert (2008). λ Boo-type stars are the main sequence of A–F dwarfs, whose atmospheres are rather stable against convection. This convection may erase any signs of chemical peculiarity at the stellar surface layers, because a nonmodified stellar gas from the stellar interior is brought to the surface and mixed with a modified gas.

According to Jofré et al. (2015), HD47536 has a mass of approximately 0.9 solar masses; hence, on the main sequence it was a G, or possibly a late F dwarf (if some mass was lost during expansion). Therefore, HD47536 may be a borderline case between λ Boo- and non-λ Boo-type stars. Moreover, even for G main-sequence stars, a mechanism of dust-gas separation in the circumstellar medium is possible. Meléndez et al. (2009) reported that according to the LTE analysis of high-quality spectral material obtained for solar twin stars, our Sun has a small deficiency of refractory elements compared to the average solar twin abundances, which may be an indication of the loss of refractory elements from the proto-Sun owing to terrestrial-like planet formation from this lost material. This conclusion was based on the dependence of the relative abundances on the condensation temperature.

A similar conclusion was drawn by Ramírez et al. (2009) for a large sample of solar twins. These published results suggest that (if the described scenario is correct) the Sun has a somewhat unusual composition of refractory elements, as such a feature is not observed on average for the solar twins. It is possible to estimate the distance from the central star at which the temperature is close to the condensation temperature of siderophile refractory elements, such as Cr, Mn, Fe, and Ni (Tcond ≈ 1,300K). For this, we use the dilution factor W = 0.5 · (1 – (1 – (r*/r)2)0.5), where r* is the radius of the central (main sequence) star. For the effective temperature of the star, that is 6,000K, we assume that the required temperature of 1,300K takes place at approximately R equals to ten stellar radii, where the formation of the terrestrial-like planets can take place. As pointed out earlier, the planet of the HD47536 system approached the main star at a distance of approximately eight radii from the red giant star.

The abundance distribution of the high condensation temperature elements (> 1,300K) resembles that observed in the solar system chondrites (see Wood et al. 2019; Fig. 4). Therefore, the relative iron abundances of the lithophile refractory elements (Mg, Al, Si, Sc, Sr, Zr, Ni, and Ba, as well as all the studied rare earth elements and thorium) had positive values (see Fig. 5). We can conclude that the observed peculiarity of the chemical element distribution in the K-giant halo star HD47536 does not contradict the suggestion that it was formed under certain conditions connected with the dust-gas separation in the pristine envelope surrounding the central proto-star. It is also possible to conclude that the processes of star formation are pertinent to those that have taken place in the primitive solar system.

Additional evidence that HD47536 has signs of remnant circumstellar shells is its detectable IRAS flux. Fig. 6 shows the fit of the broadband photometry and IRAS fluxes for our program star using the theoretical fluxes calculated by Kurucz (1993).

Fig. 6. The fit of wideband photometry and Infrared Astronomical Satellite (IRAS) fluxes at the wavelengths 12, 25, and 60 microns for HD47536 by Kurucz model flux (Kurucz 1993) with metallicity [Fe/H] = –0.5 interpolated for Teff = 4,400K, log g = 1.8. The values of photometrical magnitudes are taken from Set of Identifications, Measurements, and Bibliography for Astronomical Data (SIMBAD) database, the IRAS fluxes - from IRAS point sources catalogue (PSC) and faint sources catalogue (FSC). IRAS PSC and IRAS FSC were distributed by Astronomical Data Center [NASA/Goddard Space Flight Center (GSFC)] CD-ROMs. The software for comparison of Kurucz’a fluxes with observed photometrical values was written by one of the authors of this paper (Alexander Yushchenko). The observed IRAS flux exceeds calculated model flux distribution approximately by 60 percents at the wavelength 60 microns.
Download Original Figure
5.3 Interplay of Several Scenarios

As revealed by Greenstein (1949), the deficiency of relative abundances is observed for elements with SIP between 11.8 and 16 eV; therefore, it was concluded that resonance can be observed at the SIP=13.61.8+2.4 eV. Yushchenko et al. (2015) showed that the plot of relative abundances vs. SIPs for ρ Pup is a sawtooth graphic with two cutoffs near the values of SIP = 12 eV and SIP = 24 eV (Fig. 6 in Yushchenko et al. 2015). The sawtooth shape was also found for the main component of RR Lyn, as illustrated in Fig. 3 by Jeong et al. (2017), and for barium star ζ Cyg (Fig. 1 in Jeong et al. 2018).

For the abundance pattern of halo giant HD47536, the sawtooth picture was observed only for elements with SIPs of less than approximately 20 eV. The elements with higher SIP values in the sawtooth shape of relative abundances vs. SIP plots are not clearly visible. This may be explained by the lower density of the halo interstellar medium.

The chemical elements with SIPs between the two cutoffs, as well as the elements with SIPs less than the first cutoff and greater than the second cutoff, show positive correlations of relative abundances vs. SIPs in the above-mentioned three stars. Yushchenko et al. (2015) and Yushchenko et al. (2021) found similar correlations in the published observations of several hundred B–F type stars and disk red giants for chemical elements with SIPs between the two cutoffs and greater than the second cutoff. It was also shown that these correlations are dependent on the effective temperature and disappear for stars with convective energy transfer in their atmospheres. Note that not all relative abundances of chemical elements between the cutoffs followed the described sawtooth shape in the relative abundances vs. SIPs plots.

A portion of the elements showed higher relative abundance. In the case of ρ Pup, these are Mo (Z = 42), Ru (44), Cd (48), and elements with atomic numbers Z ≥ 71, namely Lu, Hf, Os, Ir, Pt, Au, and Pb. In the atmosphere of the primary component of RR Lyn, Lu, in the atmosphere of barium star ζ Cyg, an increase in relative abundance was observed for Mg (12), Nb (41), Mo (42), Ru (44), Pd (46), Hf (72), and Os (76). It is natural to expect that several physical scenarios occur in the atmosphere of certain stars, and the relative abundances of the above-mentioned chemical elements are influenced not only by charge-exchange reactions but also by other scenarios.

Let us return to the chemical elements with 37 ≤ Z ≤ 42 in the atmosphere of HD47536 and attempt to understand the low deviations of Rb (37), Y (39), and Nb (41) from the scaled solar system r-process izotopes fit (see Fig. 3). The SIPs of these elements were 27.29, 12.22, and 14.00 eV, respectively. This allows Rb, Y, and Nb to be influenced by charge exchange reactions, and we can expect a decrease in the relative abundance of these chemical elements with respect to s-process predictions. This is illustrated in Fig. 3. Rb has condensation Tcond = 798K, and its abundance can also be decreased owing to the dust-gas separation scenario. Sr (38) and Mo (42) had SIPs 11.03 of 16.16 eV, respectively. These values are outside the resonance region in the relative abundance vs. SIPs plots. Therefore, the relative abundances of these two chemical elements should not be strongly decreased by charge-exchange reactions, and s-process enrichment can be detected, as shown in Fig. 3.

Zr (40) has an SIP of 13.1 eV and should be influenced by charge-exchange reactions; however, Fig. 3 shows significant enrichment by the s-process. At present, we are not able to explain the abundance of this element; perhaps, additional physical effects, such as example dust-gas separation, radiative diffusion in stellar atmospheres, or other scenarios should be taken into account.

Similar considerations can be made for lanthanides. Fig. 3 shows that the first three lanthanides, La, Ce, and Pr are underabundant with respect to the solar system s-process curve. The SIPs of these three chemical elements are 11.05, 10.85, and 10.55 eV, respectively, while the SIPs of other six lanthanides investigated in the atmosphere of HD47536, specifically Nd, Sm, Eu, Gd, Dy, and Er, are from 10.72 eV for Nd to 11.07 eV for Sm and higher values for the other four elements. The highest SIPs among the investigated lanthanides were 12.09 eV for Gd.

We can conclude that the three lanthanides with the lowest SIPs exhibit negative deviations from the solar system s-process approximation. For the other six lanthanides with higher SIPs, the fit of their abundances by the solar system s-process distribution results in close to zero or positive (for Gd) deviations.

Fig. 4 shows the positive correlations between the relative abundances of lanthanides and SIPs. This allows us to suppose that the abundances of lanthanides varied owing to charge-exchange reactions, and we should detect a decrease in relative abundance for lanthanides, with the lowest SIP values. The observed relative abundances show negative deviations from the solar system s-process fit for the first three lanthanides with the lowest SIPs, close to zero deviations for the other lanthanides, and the element with the highest SIP – for Gd.

In summary, the solar system s-process isotope abundance distribution fits the observed relative abundances in the atmosphere of HD47536. Negative deviations are observed for Rb (37), Y (39), Nb (41), La (57), Ce (58), and Pr (59). The abundances of these six chemical elements can be influenced by charge exchange reactions, and we should expect a decrease in the relative abundances of these elements over time.

Note that the SIP for Th is 11.9 eV, and the abundance of this radioactive element can be changed not only because of the natural radioactive decay but also because of chargeexchange reactions. The SIP for Th was in the first highest region of the sawtooth picture observed in the plot of relative abundances against SIPs. Therefore, we can expect that the influence of charge-exchange reactions on thorium abundance is relatively weak with respect to the abundance of other chemical elements; however, this conclusion is qualitative and requires additional investigation.


The chemical composition of the planetary host halo star HD47536 was investigated. The abundances of 38 chemical elements were determined using the spectrum synthesis method. It was shown that the abundance pattern of this star is influenced by r- and s-process enrichment, chargeexchange reactions, and dust-gas separation. The interplay of these three scenarios allows us to qualitatively explain the atmospheric abundances of this star.

Our results obtained for HD47536 and other stars indicate that the plot of relative abundances of chemical elements vs. SIPs of these elements can be a sawtooth picture with two cutoffs, with SIPs near 12 or 13 eV and near 24 eV. The relative abundances of chemical elements with SIPs in three intervals, namely SIP < 12 eV, 12.5 eV < SIP < 20 eV, and SIP > 24.6 eV show the positive correlation with respect to SIP.

In summary, one can conclude that both the dust-gas separation mechanism and charge-exchange reactions discussed above can be involved in the interaction between the star’s atmosphere and circumstellar/interstellar medium; this interaction produces a specific abundance pattern, which is observed in our program star HD47536 as well as in other stars.

Note that the dependence of relative abundances on SIPs was discovered by Greenstein (1949) as soon as the use of high dispersion spectra and quantitative analysis of stellar abundances allowed for the first stellar abundance patterns. However, the explanation for this dependence remains qualitative.

It should be mentioned that Jeong et al. (2018) found the possibility of age dependence for correlations of relative abundances vs. SIPs in barium stars. This means that the abundance of stable (non-radioactive) chemical elements in stellar atmospheres varies with time. If this possibility is confirmed for other types of stars, it would be difficult to use Th/Eu or other radioactive to stable element abundance ratios to determine the age of the star. We found both Th and Eu abundances in the atmosphere of halo giant HD47536. We also found a clear dependence of lanthanide abundance on the SIP (see Fig. 4). The SIP of thorium is 11.9 eV, which is similar to that of lanthanides. It seems natural to expect that Eu and Th atoms in the atmosphere of HD47536 can be influenced by charge-exchange reactions owing to accretion of interstellar matter, and the abundance ratio of these elements can change with time owing to the natural radioactive decay of Th atoms. Therefore, we did not attempt to estimate the age of HD47536 by using the Th/Eu abundance ratio.

A detailed analysis of the correlations between the relative abundances against second ionization potentials and condensation temperatures found in this study is the addition of previous determinations of these values for different types of stars and will be used for constructing the quantitative theory of charge-exchange reactions in stellar atmospheres.


This study used the SIMBAD database, operated at CDS, Strasbourg, France. SA is thankful to the administration of ESO for financial support during his stay in Garching within a program of Ukrainian astronomers at risk and to Dr. F. Primas and N. Silva for excellent organization of their short-term visits. VK is grateful to the Vector-Stiftung at Stuttgart, Germany, for support within the program “2022 – Immediate help for Ukrainian refugee scientists” under grant P2022-0064. AD was funded by the Science Committee of the Ministry of Science and Higher Education of the Republic of Kazakhstan (grant No. AP14972694).



Adibekyan VZ, Sousa SG, Santos NC, Delgado Mena E, González Hernández JI, et al., Chemical abundances of 1111 FGK stars from the HARPS GTO planet search program, Astron. Astrophys. 545, A32 (2012).


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


Andrievsky S, The sodium abundance in λ Bootis stars, Astron. Astrophys. 449, 345-347 (2006).


Andrievsky S, Paunzen E, Towards the solution of the λ Bootis problem, Mon. Not. R. Astron. Soc. 313, 547-552 (2000).


Ardeberg A, Virdefors B, Solar line blocking for lambda lambda 4006-6860, Astron. Astrophys. Suppl. Ser. 36, 317-321 (1979).


Beirão P, Santos NC, Israelian G, Mayor M, Abundances of Na, Mg and Al in stars with giant planets, Astron. Astrophys. 438, 251- 256 (2005).


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


Borra EF, Landstreet JD, The magnetic fields of the AP stars, Astrophys. J. Suppl. Ser. 42, 421-445 (1980).


Caffau E, Bonifacio P, Faraggiana R, François P, Gratton RG, et al., Sulphur abundance in Galactic stars, Astron. Astrophys. 441, 533-548 (2005).


Charbonneau P, A simple accretion/diffusion model for lambda Bootis stars, Astrophys. J. Lett. 372, L33 (1991).


Charbonneau P, Particle transport and the lambda Bootis phenomenon. I. The diffusion/mass-loss model revisited, Astrophys. J. 405, 720 (1993).


da Silva L, Girardi L, Pasquini L, Setiawan J, von der Lühe O, et al., Basic physical parameters of a selected sample of evolved stars, Astron. Astrophys. 458, 609-623 (2006).


Delbouille L, Rolland G, Neven L, Atlas Photometrique du Spectre Solaire de λ 3000 a λ 10000 (Universite de Liege, Institut d’Astrophysique, Liege, 1973).


Ecuvillon A, Israelian G, Santos NC, Shchukina NG, Mayor M, et al., Oxygen abundances in planet-harbouring stars: comparison of different abundance indicators, Astron. Astrophys. 445, 633-645 (2006).


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


Gálvez-Ortiz MC, Delgado-Mena E, González Hernández JI, Israelian G, Santos NC, et al., Beryllium abundances in stars with planets:extending the sample, Astron. Astrophys. 530, A66 (2011).


Gilli G, Israelian G, Ecuvillon A, Santos NC, Mayor M, Abundances of refractory elements in the atmospheres of stars with extrasolar planets, Astron. Astrophys. 449, 723-736 (2006).


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


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


Havnes O, Magnetic stars as generators of cosmic rays, Astron. Astrophys. 13, 52-57 (1971).


Havnes O, Conti PS, Magnetic accretion processes in peculiar A stars, Astron. Astrophys. 14, 1-11 (1971).


Jeong Y, Yushchenko AV, Doikov DN, The interaction between accretion from the interstellar medium and accretion from the evolved binary component in barium stars, J. Astron. Space Sci. 35, 75-82 (2018).


Jeong Y, Yushchenko AV, Doikov DN, Gopka VF, Yushchenko VO, Chemical composition of RR Lyn - an eclipsing binary system with Am and λ Boo type components, J. Astron. Space Sci. 34, 75-82 (2017).


Jofré E, Petrucci R, Saffe C, Saker L, Artur de la Villarmois E, et al., Stellar parameters and chemical abundances of 223 evolved stars with and without planet, Astron. Astrophys. 574, A50 (2015).


Kurucz RL, SYNTHE Spectrum Synthesis Programs and Line Data (Smithsonian Astrophysical Observatory, Cambridge, 1993).


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


Luck RE, Abundances in the local region. I. G and K giants, Astron. J. 150, 88 (2015).


Meléndez J, Asplund M, Gustafsson B, Yong D, The peculiar solar composition and its possible relation to planet formation, Astrophys. J. 704, L66-L70 (2009).


North P, The rotation of AP stars, Astron. Astrophys. 141, 328-340 (1984).


Perottoni HD, Amarante JAS, Limberg G, Rocha-Pinto HJ, Rossi S, et al., Searching for extragalactic exoplanetary systems: the curious case of BD+20 2457, Astrophys. J. 913, L3 (2021).


Ramírez I, Meléndez J, Asplund M, Accurate abundance patterns of solar twins and analogs: does the anomalous solar chemical composition come from planet formation?, Astron. Astrophys. 508, L17-L20 (2009).


Rutten RJ, van der Zalm EBJ, Revision of solar equivalent widths, Fe I oscillator strengths and the solar iron abundance, Astron. Astrophys. Suppl. Ser. 55, 143-161 (1984).


Sadakane K, Ohnishi T, Ohkubo M, Takeda Y, Metallicities in four planet-harbouring K-type giants: HD 47536, HD 59686, HD 137759, and HD 219449, Publ. Astron. Soc. Jpn. 57, 127-133 (2005).


Setiawan J, Hatzes AP, von der Lühe O, Pasquini L, Naef D, et al., Evidence of a sub-stellar companion around HD 47536, Astron. Astrophys. 398, L19-L23 (2003).


Simmerer J, Sneden C, Cowan JJ, Collier J, Woolf VM, et al., The rise of the s-process in the galaxy, Astrophys. J. 617, 1091- 1114 (2004).


Soto MG, Jenkins JS, Jones MI, RAFT - I. Discovery of new planetary candidates and updated orbits from archival FEROS spectra, Mon. Not. R. Astron. Soc. 451, 3131-3144 (2015).


Stuerenburg S, Abundance analysis of lambda Bootis stars, Astron. Astrophys. 277, 139-154 (1993).


van Belle GT, von Braun K, Directly determined linear radii and effective temperatures of exoplanet host stars, Astrophys. J. 694, 1085-1098 (2009).


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


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


Wood BJ, Smithe DJ, Harrison T, The condensation temperatures of the elements: a reappraisal, Am. Mineral. 104, 844-856 (2019).


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, Brno, Czech Republic, 5-7 Nov 1997.


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, 59 (2015).


Yushchenko AV, Kim S, Jeong Y, Demessinova A, Yushchenko V, et al., The possible signs of hydrogen and helium accretion from interstellar medium on the atmospheres of F-K giants in the local region of the galaxy, J. Astron. Space Sci. 38, 175- 183 (2021).


Yushchenko AV, Rittipruk P, Yushchenko VA, Kang YW, The planetary host red giant HD47536 - chemical composition and signs of accretion, Odessa Astron. Publ. 26, 131-136 (2013).