Free Access
Issue
A&A
Volume 615, July 2018
Article Number A120
Number of page(s) 4
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/201832717
Published online 24 July 2018

© ESO 2018

1 Introduction

There has recently come to light a new class of eruptive events whose outbursts can be attributed to a stellar merger. These events generally display high luminosity at maximum (Bond et al. 2003) (Mbol~10,Mv~<9$M_{\textrm{bol}}\,{\sim}\,{-}{10}, M_v\,{{\,\stackrel{{<}}{\sim}\,}}\,{-}9$) as evidenced by the detection of such variables in M31 (Bond et al. 2003; Williams et al. 2015) and in other galaxies (Kasliwal 2012; Smith et al. 2016). While the nature of these eruptions was initially unclear (with nova eruptions, planet-swallowing stars, and very late thermal pulses having been suggested) the “best-bet” scenario – based on the behaviour of V1309 Sco – is the merger of two stars.

V1309 Sco is the “Rosetta Stone” of stellar mergers, in that its pre-eruptive behaviour and the subsequent eruption are not only consistent with, but conclusively point to, a stellar merger event in a contact binary. The progenitor had an orbital period of 1.4 days that decreased up to the 2008 eruption (Tylenda et al. 2011). The eruptive event that marked the merger began in 2008 March, reaching peak luminosity ~3 × 104 L (Tylenda et al. 2011). Infrared (IR) observations (McCollum et al. 2014) show that it had an IR excess at least a year before the merger. Pejcha (2014) has shown that V1309 Sco experienced mass loss through the outer Lagrangian point, which eventually obscured the binary. Zhu et al. (2013) note that the presence of a significant amount of dust around V1309 Sco suggests that the ejecta following a stellar merger is an ideal environment for dust formation and growth; they showed that ~5.2 × 10−4 M of silicate and iron dust was produced in the merger.

KIC 9832227 was identified by Molnar et al. (2015) as a potential merger candidate, and we present here a near-infrared (NIR) spectrum of this candidate merger object.

2 KIC 9832227

KIC 9832227 is both an eclipsing and a contact binary, with an orbital period of 0.458 days and amplitude ~0.2 mag in the optical. Molnar et al. (2017, hereafter M17) found that the orbital period as determined from Kepler data was significantly shorter than in previous years, and that the period derivative was rapidly decreasing, a sure signature of an impending merger; they estimate the time of merger to be 2022.2 ± 0.7. M17 determine that the primary (secondary) star has mass 1.4 M (0.32 M), radius 2.6 R (0.78 R), Teff = 5800 K (5920 K) and logg = 4.19 (4.10). There is also some evidence for the presence of a third component in the system, with orbital period 590 ± 8 days and M sini = 0, 11M. A detailed discussion of the merger potential of KIC 9832227 is given by M17.

3 Observations

3.1 IRTF observation

KIC 9832227 was observed with the SpeX spectrograph (Rayner et al. 2003) on the 3 m NASA Infra-Red Telescope Facility (IRTF), Hawaii, on 2017 July 7.54 UT (MJD 57941.54). SpeX was used in the cross-dispersed mode using a 0.5″ × 15″ slit resulting in a spectral coverage 0.77–2.50 μm, at resolving power R = λδλ = 1 200. The total integration time was 717 s. The A0V star HD194354 was used as the telluric standard. The data were reduced and calibrated using the Spextool software (Cushing et al. 2005), and corrections for telluric absorption were performed using the IDL tool Xtellcor (Vacca et al. 2003). The observed spectrum, which was dereddened by E(BV) = 0.03 (M17), is shown in Fig. 1.

thumbnail Fig. 1

Left panel: dereddened IRTF spectrum of KIC 9832227, together with the fit with b = 1; the case b = 0 is also shown for comparison. See text for details. Right panel: minimisation parameter S as a functionof b.

3.2 Complementary observations

KIC 9832227 is included in the 2MASS (Skrutskie et al. 2006) and WISE (Wright et al. 2010) IR surveys; the source was not detected in other surveys, such as IRAS, AKARI, and Herschel PACS. The 2MASS and WISE fluxes are given in Table 1 and shown in Fig. 2; KIC 9832227 was not detected in WISE Band 4 (22 μm). The IR spectral energy distribution (SED) is photospheric out to at least 12 μm, indicating that dust – if present in significant quantities, as in V1309 Sco – has temperature ~<425${{\,\stackrel{{<}}{\sim}\,}}425$ K.

Table 1

Infrared fluxes from 2MASS and WISE.

thumbnail Fig. 2

IR SED of KIC 9832227 from 2MASS (red) and WISE (blue) surveys. Black curve is a 5930 K black body, corresponding to the stellar components of KIC 9832227.

3.3 Orbital phase of KIC 9832227 at the time of the IRTF observation

KIC 9832227 is an eclipsing binary: in order to interpret the IR spectrum, we need to know the orbital phase and stellar velocity separation at the time of the IRTF observation. M17 do not give the full information necessary to deduce these quantities, so this process was non-trivial. We neglected the effect of the third body, which affects the orbital timings by only ±0.002 cycles (Fig. 10 in M17), and the proposed three-body solutions suggested by M17 because this body was not detected in the spectrum.

We adopted the parameters of the exponential fit given in Fig. 12 of M17 and the linear orbital ephemeris from the Kepler Eclipsing Binary Catalogue1 (see Kirk et al. 2016). We further imposed a normalisation which forces the first point in the year 2004 to be at zero OC. The resulting orbital phases correctly reproduce the OC, instantaneous orbital period, and dP∕dt as plotted in Figs. 12, 13 and 14 in M17.

The UTC time of mid-observation of our spectrum (2 457 942.04199) was converted into the BJD(TDB) timescale using the IDL routines from Eastman et al. (2010), giving a time of 2 457 942.04135. This corresponds to an orbital phase of 0.010, so our spectrum was obtained around a time of inferior conjunction. The velocity separation of the two stars at this time was negligible compared to their spectral line broadening, so the lines of the two stars are superposed in the spectrum.

4 Theoretical spectra and SEDs

Synthetic spectra and SEDs were computed using the WITA6 program for the SAM12 model one-dimensional (1D) atmospheres (Pavlenko 2003) with Te(in K)/ logg/[Fe/H] = 5920/4.10/0.0 and 5800/4.19/0.0 for the hotter (A) and cooler (B) components, respectively, as determined by M17. Significant molecular features are not seen in the observed spectrum, therefore we accounted for only atomic lines from the VALD3 database (Ryabchikova et al. 2015). Our procedure for producing the theoretical spectra is described elsewhere (see Pavlenko 2003; Pavlenko et al. 2008). The synthetic spectra were convolved with a Gaussian profile having the instrumental broadening corresponding to the spectral resolution R, and rotational profiles corresponding to v sin i = 85 km s−1 (A) and 150 km s−1 (B; M17). To account for rotational broadening of the spectra we used formulae given by Gray (1976).

5 Results

5.1 The overall spectrum

We computed theoretical SEDs for the A and B components, and created a grid of combined spectra Ft = b × FA + (1 − b) × FB, where 0 < b < 1. We computed the best fits to the dereddened spectrum for the cases b = 0.0, 0.2, 0.4, 0.6, 0.8, 1.0; the procedure to determine the best fit used a χ2 minimising algorithm, details of which are given elsewhere (see Pavlenko et al. 2008). The minimisation parameter S= (FλtFλobs)2$S = \sum(F_{\lambda}^{\textrm{t}}-F_{\lambda}^{\textrm{obs}})^2$ indicates a minimum S at b = 1 (see Fig. 1). This is consistent with the fact that the observed spectrum was obtained at conjunction, when the hotter component dominates the combined spectrum (although a more compelling reason for not seeing the second component is the very similar Teff). Furthermore, as the eclipses are very shallow, component A dominates the spectrum at all orbital phases.

5.2 Absorption lines in KIC 9832227’s spectrum

Generally speaking, both the intermediate resolution of the observed spectra (R = 1200) and the rapid rotation of both components complicate the use of atomic lines for abundance analysis. However, we see that practically all notable lines in the NIR spectrum are well fitted with close-to-solar abundances; the likely errors do not exceed 0.3 dex. In the right panel of Fig. 3 we show a few well-fitted Si I lines in the spectral range 10 500–11 000 Å. We find fits of the same quality to metallic lines in other spectral regions.

We draw particular attention to the modelling of the moderately strong hydrogen lines, to which our model provides a satisfactory fit. In the left panel of Fig. 3 we show the hydrogen lines computed for the case of lower log g = 2.0. In the observed spectrum the upper members of the Brackett series are much weaker; this can be taken as independent evidence of the higher log g, as determined by M17. Indeed, our fits of the theoretical spectra computed for log g = 4.19 and log g = 4.10 to the observed spectra are far superior.

We also note the presence of a broad feature at 10 700 Å (identified by “??” in Fig. 3). Given the width of the feature, a possible identification is C I (3 P −2S) at 10 686.01, 10 688.27, 10 694.17 Å; these lines have relatively high excitation potentials (~ 7.5 eV).

thumbnail Fig. 3

Left panel: best fit to the observed spectrum with b = 1 in the region of the Ca II triplet and He I; the case of b = 0 is also shown, as is the hydrogen line spectrum computed for a model atmosphere with log g = 2.0. Right panel:fit to the spectrum containing Si I lines, marked by vertical arrows; the He I line is also shown. We note the broad feature at 10 700 Å, marked by “??”; see text.

5.3 The Ca II triplet: evidence of a chromosphere

Absorption in the Ca II triplet (8 498.023, 8 542.091, 8 662.14 Å) – by comparison with other metallic lines – is a notable feature of KIC 9832227’s spectrum (see Fig. 3). This figure shows the relevant portion of Fig. 1 on a larger scale, to compare the observed Ca II lines with those computed for the case of a classical model atmosphere. The observed lines are notably weaker by comparison with their computed counterparts. The discrepancies between computed and observed Ca II line strengths may be explained by several factors, such as NLTE effects or a Ca deficit. In the former case, NLTE effects in other lines should also be apparent, but this is not the case. Also H I and metallic lines are well fitted for the case of solar abundances.

The components of the KIC 9832227 system consist of rapidly rotating, tightly bound stars in a contact binary system and we would not expect the atmospheres of such stars so be adequately described by the classical model. In particular, rapid stellar rotation is likely to lead to strong chromospheric (and related) activity, which could certainly account for the weakness of the Ca triplet.

We have simulated a chromosphere by imposing a temperature minimum of 5000 K, and a temperature gradient G = dT∕d(logm) (m being the column density) in the model atmosphere. We fit the Ca II line profiles by choosing the value of G that minimises the fitting parameter S (details of our approach to this calculation are given in Pavlenko 1998). A model with a weak chromosphere (G = −200, cf. the solar value G = −800, with T in K, m in g cm−2) provides a good fit to the data in the region of the CaII triplet (see Fig. 4).

We also note here an alternative interpretation: our classical model atmospheres were computed for Rosseland mean optical depth in the range τRoss = 10−5–100 but such models lead to Ca II lines that are too strong. Interestingly, a “stripped” model atmosphere for the A component that lacks atmospheric layers in the τRoss range [10−5 –3 × 10−2] provides an equally good fit to the Ca II triplet (see Fig. 4).

However, while both approaches provide equally good fits to the data, for reasons to be discussed below, we prefer the chromospheric interpretation.

thumbnail Fig. 4

Top panel: structure of the “classical” and “chromospheric” model atmosphere of component A. Middle panel: best fit to the observed Ca II triplet, with classical and “chromospheric” model atmospheres. Bottom panel: fit to the strongest component of the Ca II triplet shown on a larger scale. In all three panels, the “stripped model atmosphere”refers to the model briefly discussed at the end of Sect. 5.3.

5.4 He I at 10 830 Å

A strong absorption feature is present at 10 830 Å (see Fig. 3); we attribute this to He I 10 830 Å (3S −3P0), which is also observed in absorption in the spectra of metal-poor subdwarfs (Smith et al. 2012) and red giants (Smith et al. 2014). This line is noticeably narrower than the rotationally broadened photospheric lines, and most likely forms outside the stellar photospheres, possibly in the common envelope.

In the case of the Sun, the line is also observed in absorption; Avrett et al. (1994) suggested that coronal radiation penetrates into the upper solar chromosphere, causing sufficient helium ionisation to populate the lower level of the He I 10 830 Å transition, thus resulting in optically thin absorption against the photospheric continuum. This, together with the presence of the C I lines at ~ 10 690Å, lends support to the chromospheric explanation offered in Sect. 5.3.

6 Conclusions

We have carried out an analysis of the NIR spectrum of KIC 9832227. We computed synthetic spectra for two classical 1D model atmospheresand fitted the dereddened spectrum. While we acknowledge that the computed model atmospheres may be overly simplistic for the case of the complex KIC 9832227 system, we nonetheless obtain a satisfactory reproduction of the observed SEDs. This may be taken as evidence that, at the time of observation, the 1D+LTE approach is valid.

Within the constraints of our simplified model, we conclude that our ephemeris analysis and spectral synthesis show that, at the time of observation, the contribution of the hotter component prevails in the combined spectrum. The strength of the H I lines is well reproduced by our models, providing independent confirmation of the high log g (4.10, component A; 4.19, component B; M17) at least in the atmosphere of component A; model atmosphereswith lower logg result in Brackett and Paschen lines that are too srong. A comparatively narrow feature at 10 830 Å is identified as the He I (3 S −3P) line, and is formed in the outer part of the common envelope. A satisfactory fit to the Ca II IR triplet is obtained by introducing a chromosphere into the model atmosphere; this is likely a consequence of rapid stellar rotation. Despite the intermediate resolution of observed spectra and the fast rotation of both components, we confirm that component A has near-solar abundances.

Acknowledgements

We warmly thank David Sand for making possible the IRTF observations, and for help with the data reduction. The Infrared Telescope Facility is operated by the University of Hawaii under contract NNH14CK55B with the National Aeronautics and Space administration. We thank the referee for their helpful comments.

References

  1. Avrett, E. H., Fontenla, J. M., & Loeser, R. 1994, in Infrared Solar Physics, eds. D. M. Rabin, J. T. Jefferies, & C. Lindsey (Dordrecht: Kluwer Academic Publishers), IAU Symp., 154, 35 [Google Scholar]
  2. Bond, H. E., Henden, A., Levay, Z. G., et al. 2003, Nature, 422, 405 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  3. Cushing M. C., Rayner J. T., & Vacca W. D. 2005, ApJ, 623, 1115 [NASA ADS] [CrossRef] [Google Scholar]
  4. Eastman, J., Siverd, R., & Gaudi, B. S. 2010, PASP, 122, 935 [NASA ADS] [CrossRef] [Google Scholar]
  5. Gray, D. F. 1976, The Observation and Analysis of Stellar Photospheres (New York: Wiley Interscience) [Google Scholar]
  6. Kasliwal, M. M. 2012, PASA, 29, 482 [NASA ADS] [CrossRef] [Google Scholar]
  7. Kirk, B., Conroy, K. E., Prša, A., et al. 2016, AJ, 151, 68 [NASA ADS] [CrossRef] [Google Scholar]
  8. McCollum, B., Laine, S., Väisänen, P., et al. 2014, AJ, 147, 11 [NASA ADS] [CrossRef] [Google Scholar]
  9. Molnar, L. A., Van Noord, D. M., Steenwyk, S. D., Spedden, C. J., & Kinemuchi, K. 2015, BAAS, 225, 415.05 [Google Scholar]
  10. Molnar, L. A., Van Noord, D. M., Kinemuchi, K., et al. 2017, ApJ, 840, 1 [NASA ADS] [CrossRef] [Google Scholar]
  11. Nicholls, C. P., Melis, C., Soszyński, I., et al. 2013, MNRAS, 431, L33 [NASA ADS] [CrossRef] [Google Scholar]
  12. Pavlenko, Ya. V. 1998, Astron. Rep., 42, 501 [Google Scholar]
  13. Pavlenko, Ya. V. 2003, Astron. Rep., 47, 59 [NASA ADS] [CrossRef] [Google Scholar]
  14. Pavlenko, Ya. V., Evans, A., Kerr, T., et al. 2008, A&A, 485, 541 [CrossRef] [EDP Sciences] [Google Scholar]
  15. Pejcha, O. 2014, ApJ, 788, 22 [NASA ADS] [CrossRef] [Google Scholar]
  16. Rayner J. T., Toomey D. W., Onaka P. M., et al. 2003, PASP, 115, 389 [NASA ADS] [CrossRef] [Google Scholar]
  17. Ryabchikova, T., Piskunov, N., Kurucz, R. L., et al. 2015, Phys. Scr., 90, 054005 [NASA ADS] [CrossRef] [Google Scholar]
  18. Smith, G. H., Dupree, A. K., & Strader, J. 2012, PASP, 124, 1252 [CrossRef] [Google Scholar]
  19. Smith, G. H., Dupree, A. K., & Strader, J. 2014, PASP, 126, 901 [Google Scholar]
  20. Smith, N., Andrews, J. E., Van Dyk, S. D., et al. 2016, MNRAS, 458, 950 [NASA ADS] [CrossRef] [Google Scholar]
  21. Strutskie, M. F., Cutri, R. M., Stiening, R. et al. 2006, AJ, 131, 1163 [NASA ADS] [CrossRef] [Google Scholar]
  22. Tylenda, R., Hajduk, M., Kamiński, T., et al. 2011, A&A, 528, A114 [Google Scholar]
  23. Vacca, W. D., Cushing, M. C., & Rayner, J. T. 2003, PASP, 115, 389 [NASA ADS] [CrossRef] [Google Scholar]
  24. Williams, S. C., Darnley, M. J., Bode, M. F., & Steele, I. A. 2015, ApJ, 805, L18 [NASA ADS] [CrossRef] [Google Scholar]
  25. Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868 [NASA ADS] [CrossRef] [Google Scholar]
  26. Zhu, C., Lü, G., & Wang, Z. 2013, ApJ, 777, 23 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Infrared fluxes from 2MASS and WISE.

All Figures

thumbnail Fig. 1

Left panel: dereddened IRTF spectrum of KIC 9832227, together with the fit with b = 1; the case b = 0 is also shown for comparison. See text for details. Right panel: minimisation parameter S as a functionof b.

In the text
thumbnail Fig. 2

IR SED of KIC 9832227 from 2MASS (red) and WISE (blue) surveys. Black curve is a 5930 K black body, corresponding to the stellar components of KIC 9832227.

In the text
thumbnail Fig. 3

Left panel: best fit to the observed spectrum with b = 1 in the region of the Ca II triplet and He I; the case of b = 0 is also shown, as is the hydrogen line spectrum computed for a model atmosphere with log g = 2.0. Right panel:fit to the spectrum containing Si I lines, marked by vertical arrows; the He I line is also shown. We note the broad feature at 10 700 Å, marked by “??”; see text.

In the text
thumbnail Fig. 4

Top panel: structure of the “classical” and “chromospheric” model atmosphere of component A. Middle panel: best fit to the observed Ca II triplet, with classical and “chromospheric” model atmospheres. Bottom panel: fit to the strongest component of the Ca II triplet shown on a larger scale. In all three panels, the “stripped model atmosphere”refers to the model briefly discussed at the end of Sect. 5.3.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.