A New Infrared Photometric Study of Intrinsic and Extrinsic S-type Stars

We collect all known intrinsic and extrinsic S-type stars to discuss their infrared properties and find their difference in the infrared using photometric data from the Two Micron All Sky Survey, Wide Field Infrared Survey Explorer, and Infrared Astronomical Satellite missions. Then we look for the diagnosis to extract intrinsic S-type stars from the large unclassified sample. We found that, statistically, intrinsic S-type stars have larger infrared excesses than extrinsic S-type stars in the wavelength region of 1–60 μm due to thicker dusty circumstellar envelopes. We also found that only intrinsic S-type stars occupy the reddest color areas in all of the two-color diagrams we presented. Finally, 172 new intrinsic S-type stars are presented in this paper. This makes the number of known intrinsic S-type stars almost double. In addition, some intrinsic and extrinsic S-type stars have power-law distributions in some two-color diagrams with the wavelength longer than 5 μm. The possible reason for this is also discussed.


Introduction
S-type stars have been traditionally considered for a long time as intermediate red giants between M-type stars and carbon stars as in the evolution sequence M-S-C at the asymptotic giant branch (AGB) phase. The spectra of S-type stars are characterized by the absorption bands of ZrO and LaO molecules in the photosphere. However, About 30 yr ago this idea has been changed that some nonvariable S-type stars without the element of Technetium (Tc) may be in the binary system (Jorissen & Mayor 1988, 1992 at the red giant branch (RGB) or the thermally pulsing asymptotic giant branch (TP-AGB) phases (Busso et al. 1992). In fact, Iben & Renzini (1983) hinted the existence of such S-type stars. Now we know that there are actually two categories of S-type stars. They are as follows.
(1) Intrinsic S-type stars that follow the M-S-C evolution sequence and show Tc lines in their optical spectra. They are also called Tc-rich S-type stars or S-type stars with Tc.
(2) Extrinsic S-type stars are in a binary system with a white dwarf (WD) and have enhanced s-process elements at the stellar surface as the result of mass transfer from a former AGB star (now the WD) to a less evolved companion (now the S-type star), but without Tc lines in their optical spectra. They are also called Tc-deficient S-type stars or S-type stars without Tc.
In the intrinsic S-type star the third dredge-up process (the signature of s-process elements) is presently at work while in the extrinsic S-type star s-process elements are transferred to a less evolved companion in a binary system. Among elements present in the atmospheres of intrinsic S-type stars there is the element Tc with a half lifetime of about 2 * 10 5 yr, while in the extrinsic S-type star this element has already decayed. Therefore the presence of Tc is well discriminating the intrinsic S-type star from the extrinsic S-type star (Iben & Renzini 1983;Kipper 1991;Johnson et al. 1993).
Several methods can be used to distinguish between extrinsic S-type stars and intrinsic S-type stars; a direct way is to detect Tc lines (Jorissen 2004). However, this is generally difficult since these lines are very weak and lie in a spectral region crowded with many other lines. Fortunately, some indirect methods have been suggested such as: (1) for extrinsic S-type stars the WD companion in a binary system can be directly detected in the UV region if it is warm enough. This has been done by Johnson et al. (1990Johnson et al. ( , 1993. However, up to now only a small amount of S-type stars have UV observations from the International Ultraviolet Explorer. (2) Variations of the periodic radial velocity can be observed for extrinsic S-type stars in a binary system (Jorissen & Mayor 1988, 1992. However, it is limited by their long orbital periods. (3) The high-excitation emission line of He I (10830 Å) has been detected in the spectra of extrinsic S-type stars, while in the spectra of intrinsic S-type stars this line is lacking ). However, this line almost reaches the cutoff wavelength of the response curve in the current CCD. (4) The different infrared colors of the two categories of S-type stars can serve as another important criterion (Chen & Kwok 1993;Groenewegen 1993;Jorissen et al. 1993). Johnson (1992) well summarized these methods to distinguish between intrinsic S-type stars and extrinsic S-type stars.
The unstable element Tc was first identified for some M, S, and carbon stars by Merrill (1952). Little-Marenin & Little (1979) first made the spectral observation of 90 late-type stars (including 10 S-type and MS-type stars) and concluded that the stars with Tc are in Population I, while the stars without Tc are in Population II. Iben & Renzini (1983) discussed the status of the element Tc in the AGB phase and hinted at the existence of extrinsic S-type stars in the binary system. Then Little et al. (1987) made the spectral observation for about 280 late-type stars (including 24 S-type and MS-type stars) to identify some intrinsic S-type stars and extrinsic S-type stars. They pointed out that single S-type stars all have Tc and that extrinsic S-type stars should be in the binary system.  gave high-resolution spectra for 40 MS-type and S-type stars to search for the presence of Tc. They suggested that intrinsic S-type stars are currently thermally pulsing AGB stars undergoing third dredge-up while extrinsic S-type stars represent the coolest members of the barium star class. Brown et al. (1990) made observations of the radial velocity and the He I (10830 Å) triplet for some MS-type and S-type stars. They suggested that the He I features and radial-velocity variations in extrinsic MS-type and S-type stars are both caused by the presence of WD companions. Groenewegen (1993) studied the discrimination between intrinsic S-type stars and extrinsic S-type stars using their infrared properties from the visual and Infrared Astronomical Satellite (IRAS) observations. Groenewegen found that intrinsic S-type stars have larger infrared excesses from circumstellar envelopes while extrinsic S-type stars are consistent with stellar blackbodies with no or very little infrared excesses. Groenewegen also suggested that extrinsic S-type stars might be 50%-75% of S-type stars in an optically complete sample. Jorissen et al. (1993) investigated the infrared properties for some S-type stars using the Caltech Infrared Catalog (IRC) and the IRAS data. They pointed out that the frequency of intrinsic S-type stars may be less than 50%. Chen & Kwok (1993) studied the circumstellar properties for about 150 S-type stars using the IRAS lowresolution spectra (LRS). They found that the silicate emission feature only appeared for intrinsic S-type stars. Chen et al. (1998) found that an S-type star whose energy distribution can be fitted by a double blackbody curve is likely to be an intrinsic S-type star, whereas candidates of extrinsic S-type stars can be found from sources whose energy distribution can be fitted by a single blackbody curve. Wang & Chen (2002) made the JHK observations for a number of S-type stars. Together with the IRAS photometric and the LRS data they found that about 100 stars are very like extrinsic S-type stars. Lebzelter & Hron (2003) confirmed that only a small fraction of the semiregular variables show Tc lines. However, they also found a significant number of Miras without Tc. Yang et al. (2006) searched the new candidates of intrinsic and extrinsic S-type stars using the Two Micron All Sky Survey (2MASS), IRAS, and Midcourse Space Experiment (MSX) data. Finally they found about 150 candidates of extrinsic S-type stars and 250 candidates of intrinsic S-type stars. Recently infrared properties of about 600 S-type and related stars have been investigated by Guandalini & Busso (2008) using 2MASS, IRAS, Infrared Space Observatory (ISO), and MSX data. They identified about 40 extrinsic S-type stars and about 120 intrinsic S-type stars.
In this paper we attempt to collect all known intrinsic S-type stars and extrinsic S-type stars, and investigate the difference of their infrared properties using 2MASS, Wide Field Infrared Survey Explorer (WISE), and IRAS data. Then the diagnosis to distinguish both kinds of S-type stars is discussed.
It is noted that the WISE data are used first time for the study of these S-type stars in this paper. The WISE mission has completed an all-sky survey in the mid-infrared (Wright et al. 2010). WISE performed observations in four bands, W1 (3.4 μm), W2 (4.6 μm), W3 (12 μm), and W4 (22 μm), and WISE all-sky data were released in 2012 March 14. This survey extended the 2MASS All Sky Survey into the mid-infrared and connected it with far-infrared observations by IRAS.

Sample Selection and Data Processing
We have collected almost all of the known intrinsic S-type stars and extrinsic S-type stars from the literature (the candidates are not included) to study their infrared color properties using 2MASS, WISE, and IRAS data. Finally 151 certain extrinsic S-type stars and 190 certain intrinsic S-type stars are found. The infrared data for extrinsic S-type stars and intrinsic S-type stars are listed in Tables 1 and 2, respectively. In Tables 1 and 2 the contents are: (1) the star number in this paper; (2) the number in the General Catalog of Galactic S-Stars (GCGSS, simply CSS in this paper) from Stepheson (1984) and Stephenson (1990, simply CSS2); (3) the star name; (4) the star position in the epoch of 2000; (5) the 2MASS JHK magnitudes, where if the empty is showing in the observational error position, it means the upper limit value, which is not used for the data analysis below); (6) the WISE magnitudes in four bands, where if the empty is showing in the observational error position, it means the upper limit value, which is not used for the data analysis below; (7) the Galactic extinction coefficient Av derived in this paper; and (8) the identified IRAS Point Source Catalog/Faint Source Catalog counterpart. In addition, in the note of Table 1 the reference of the source origin is indicated for both Tables 1 and 2. Note that in the context only small parts of Tables 1 and 2 are shown. The full versions of Tables 1 and 2 can be found in the supplemental material as the online data.
The cross-identifications between S-type stars and 2MASS/ WISE counterparts are made from Cutri et al. (2012) using the radius of 2″. The cross-identifications of IRAS counterparts are made according to the positional error ellipse of the source, because it has a 95% confidence level (Beichman et al. 1988). We try to extract observational data from the AKARI and Hershel missions, but only very few counterparts are found. Therefore, for statistical purposes those data are not used in this paper.
For the data from 2MASS and WISE, the Galactic extinction correction is needed. We use the Galactic extinction law of Schlegel et al. (1998) for 2MASS data and Yuan et al. (2013) for WISE data, respectively. The Galactic extinction is corrected according to the relation from van Herk (1965), where b is the Galactic latitude and d is the distance in parsecs. We also know another relation, where M K and A K are the absolute magnitude and the Galactic interstellar extinction coefficient in the K band, respectively. In addition, from Schlegel et al. (1998) Feast et al. (1982) pointed out that M K =−8.0 can be adopted for carbon stars, S-type stars, and related stars, thus iterations for d and A K (Av) are made to derive the final A K (Av). Note that the 2MASS and WISE data used below are all corrected for the Galactic extinction according to the method here.

2MASS, WISE, and IRAS Two-color Diagrams for Intrinsic S-type Stars and Extrinsic S-type Stars
In order to reveal the infrared properties of intrinsic S-type stars and extrinsic S-type stars, several two-color diagrams using the 2MASS, WISE, and IRAS photometry data from Tables 1 and 2 are constructed. It is noted that in some previous papers the 2MASS JHK data are often combined with IRAS data, and even the visual data and MSX data, to construct the two-color diagrams and obtain the results. However, many S-type stars are Mira-type variables or semiregular variables. In order to avoid the errors due to observational data used in different time, we use only the data of one mission in each of the two-color diagrams for the analysis in this paper.
Although, as showing in Section 1 of this paper, some authors already used 2MASS and IRAS data to discuss infrared properties of intrinsic S-type stars and extrinsic S-type stars, but samples in this paper are larger to get more convincing results.
In addition, in order to show infrared properties clearly, the blackbody distributions and the power-law distributions are    also presented in these two-color diagrams using the method from Chen et al. (2016).

2MASS Two-color Diagram
Taking 2MASS data from Tables 1 and 2, the (J-H) versus (H-K ) diagram can be shown in Figure 1 for known intrinsic S-type stars and extrinsic S-type stars. From Figure 1 we find that: (1) almost all sources are located around or in the upper left of the blackbody line, indicative of the thermal radiations from the stellar photosphere and/or dust in circumstellar envelopes; and (2) statistically, intrinsic S-type stars have much larger infrared excess than extrinsic S-type stars. In the region of H-K>0.6 and J-H>1.2 no extrinsic S-type stars appear.

WISE Two-color Diagrams
Taking WISE data from Tables 1 and 2 the one of WISE twocolor diagram, the (W2-W3) versus (W1-W2) diagram can be plotted in Figure 2 for known intrinsic S-type stars and extrinsic S-type stars. From Figure 2 it is seen that: (1) the majority of both types of sources are located around or in the bottom right region of the blackbody line, indicative of dust thermal radiations from circumstellar envelopes; (2) the most of extrinsic S-type stars are distributed in the region around W2-W3=0 indicating no infrared excess in this color; (3) it is unexpected that some sources are located in the left region of the power-law line with the high temperate; (4) statistically, intrinsic S-type stars have much larger infrared excesses than   extrinsic S-type stars; and (5) in the region of W2-W3>1.0 all sources are intrinsic S-type stars. Another WISE two-color diagram, the (W3-W4) versus (W2-W3) diagram can be plotted in Figure 3 for known intrinsic S-type stars and extrinsic S-type stars. It can be seen from Figure 3 that: (1) it is unexpected that over half of the sources are located around the power-law line; (2) statistically, intrinsic S-type stars have much larger infrared excesses than extrinsic S-type stars either in the W2-W3 color or the W3-W4 color; and (3) again, in the region of W2-W3>1.0 all of the sources are intrinsic S-type stars.

IRAS Two-color Diagram
From Table 1 it can be seen that 138 (out of 151) extrinsic S-type stars have IRAS counterparts. However, only 16 sources have good quality data in all 12, 25, and 60 μm bands. From Table 2 it is found that 183 (out of 190) intrinsic S-type stars have IRAS counterparts. Fortunately, 78 of them have good quality data in all 12, 25, and 60 μm bands. Taking IRAS data for these two groups of stars and using the equations (Beichman et al. 1988 Figure 4. Here in Equations (3) and (4) the units of IRAS fluxes in 12, 25, and 60 μm are jansky. It can be seen that: (1) it is also unexpected that about half of both kinds of sources are distributed along with the power-law line; (2) statistically, intrinsic S-type stars have much larger infrared excess than extrinsic S-type stars either in the In summary, statistically, in the 1-60 μm wavelength region intrinsic S-type stars in the AGB phase have much larger infrared excesses than extrinsic S-type stars in the RGB/TP-AGB phase due to thick dust circumstellar envelopes for formers, as expected. Another important result is that in the reddest color areas from all of the two-color diagrams above only intrinsic S-type stars can be found. In addition, in the wavelength longer than 5 μm some intrinsic S-stars and extrinsic S-stars have power-law distributions.

Discussion on the Power-law Distribution
It is seen from Figures 2-4 that some intrinsic and extrinsic S-stars have power-law distributions in the wavelength longer than 5 μm.
S-type stars often exhibit substantial mass loss through the dust-driven or pulsation-driven wind in the surrounding envelopes with gas and dust. The dust absorbs the stellar radiation and reradiate in the infrared (see, e.g., Hony et al. 2009). Therefore the possible reasons for power-law distributions for some S-type stars are as follows.
(1) One possibility is that it is caused by the disk-like structure of their dusty envelopes. The two typical S-type stars, π 1 Gruis and RS Cnc, are proved to have the disklike dusty envelopes with bipolar outflows (Hirano et al. 2005;Cruzalebes &Sacuto 2006 andLibert et al. 2010). Perhaps those S-type stars with power-law distributions in this paper have envelope structures similar to π 1 Gruis and RS Cnc.
(2) Another possibility is that it is probably caused by the free-free/bound-free emissions in the stellar wind. Perhaps, besides the neutral gas, these stars may also have the ionized gases in their circumstellar envelopes to produce such emissions. Shetye et al. (2019) indeed detected the ionized material in some S-type stars.

Diagnosis for Finding Intrinsic S-type Stars
In order to further discuss classifications of S-type stars not listed in Tables 1 and 2, the 2MASS, WISE, and IRAS counterparts for all stars in CSS (Stepheson 1984) and CSS2 (Stepheson 1990) are cross-identified, and the Galactic extinctions are obtained for all 2MASS/WISE counterparts here according to the method in Section 2. Note.The source with the symbol " * " appears in Table 3.
(This table is available in machine-readable form.) Figures 1-4 are together replotted in Figure 5. We can see that if we set out a line of Y = −1.50X + 1.90 in the 2MASS diagram of Figure 5, in the upper right region of this line only intrinsic S-type stars are found. Thus we can consider that for unclassified S-type stars this line can be used to extract intrinsic S-type stars. All unclassified S-type stars listed in CSS Note.The source with the symbol " * " appears in Table 3. The source with the symbol "#" appears in Table 4.
(This table is available in machine-readable form.) (Stepheson 1984) and CSS2 (Stepheson 1990) are tested in the 2MASS JHK diagram. Finally, 122 sources in these catalogs are identified as new intrinsic S-type stars, which are listed in Table 3. The contents in Table 3 are: (1) the number in CSS or CSS2; (2) the JHK magnitudes; (3) the derived Galactic extinction coefficient Av; and (4) the de-reddening colors of H-K and J-H.
In the WISE first and second diagrams it is clear that in the region with the W2-W3 color >1.0 only intrinsic S-type stars are found. Thus this is another diagnosis to extract intrinsic S-type stars from the unclassified samples in the CSS and CSS2. Finally, 35 sources in these catalogs are found as new intrinsic S-type stars that are listed in Table 4. The contents of Table 4 are: (1) the number in CSS or CSS2; (2) the magnitudes in W2 and W3; (3) the Av derived; and (4) the de-reddening color of W2-W3. In the table note the symbol " * " indicates the source already listed in Table 3. It is seen that 35 sources are found to be new intrinsic S-type stars and 25 sources are already listed in Table 3.
If we draw a line of Y = −3X + 1.4 in the IRAS diagram of Figure 5, it can be seen that only intrinsic S-type stars are located in the right region of this line. Thus we can consider that for unclassified S-type stars this line can be used to extract intrinsic S-type stars. All unclassified S-type stars listed in CSS (Stepheson 1984) and CSS2 (Stepheson 1990) are tested in this diagram. Finally, 57 sources in these catalogs are identified as new intrinsic S-type stars, which are listed in Table 5. The contents in Table 5 Table 3 and the symbol "#" indicates the source already listed in Table 4. It can be seen that 57 sources are found to be new intrinsic S-type stars, and among them 28 sources are already listed in Tables 3  and 12 sources are already listed in Table 4.
It is seen from Tables 3-5 that 172 sources in total are newly identified intrinsic S-type stars in this paper. It can be seen that the number of previously known intrinsic S-type stars is 190 from Table 2, and the number of new intrinsic S-type stars found in this paper is 172. It means that the number of known intrinsic S-type stars is almost double after this paper.

Summary
We have collected 151 certain extrinsic S-type stars and 190 certain intrinsic S-type stars to study their infrared properties using several two-color diagrams from 2MASS, WISE, and IRAS data in this paper.
The results in this paper reveal that, statistically, intrinsic S-type stars have larger infrared excesses than extrinsic S-type stars in the wavelength region of 1-60 μm due to more thick dust circumstellar envelopes. It is found that in all two-color diagrams, only intrinsic S-type stars occupy the reddest color areas. This result can be used as a diagnostic to extract intrinsic S-type stars. Finally, 172 new intrinsic S-type stars are presented. This makes the number of known intrinsic S-type stars almost double.
In addition, some intrinsic and extrinsic S stars have powerlaw distributions in the two-color diagrams in the wavelength longer than 5 μm. The possible reasons for this are also discussed.