Arxiv:1201.0843V1 [Astro-Ph.SR] 4 Jan 2012 H Ipro Aaaewihbln Othe in Stars
Total Page:16
File Type:pdf, Size:1020Kb
Draft version January 5, 2012 A Preprint typeset using LTEX style emulateapj v. 5/2/11 ASTEROSEISMOLOGY OF THE NEARBY SN-II PROGENITOR: RIGEL PART I. THE MOST HIGH PRECISION PHOTOMETRY AND RADIAL VELOCITY MONITORING Ehsan Moravveji1,2, Edward F. Guinan2, Matt Shultz3, Michael H. Williamson4, and Andres Moya5 (Dated: last updated: January 5, 2012) Draft version January 5, 2012 ABSTRACT Rigel (β Ori, B8 Ia) is a nearby blue supergiant displaying α Cyg type variability, and is one of the nearest type-II supernova progenitors. As such it is an excellent test bed to study the internal structure of pre core-collapse stars. In this study, for the first time, we present 28 days of high precision MOST photometry and over 6 years of spectroscopic monitoring. We report nineteen significant pulsation modes of SNR&4.6 from radial velocities, with variability time scales ranging from 1.21 to 74.7 days, which are associated with high order low degree gravity modes. While the radial velocity variations show a degree of correlation with the flux changes, there is no clear interplay between the equivalent widths of different metallic and Hα lines. 1. INTRODUCTION days. They were later definitely identified as gravity Blue supergiants (BSGs) such as Rigel are post-main mode oscillations by a Non-LTE spectroscopic analysis of Lefever et al. (2007) who compared the current position sequence (MS) massive stars (M & 15M⊙) evolving quickly across the Hertzsprung-Russel diagram. They of their stars with the instability strip of g-mode domi- are amongst the most intensively studied objects in con- nated pulsators. Beyond the Milky Way, Bresolin et al. temporary astrophysics (see Maeder & Meynet 2011, for (2004) detected a handful of variable BSGs in NGC 300, a recent review). BSG stars are intriguing objects be- and found the periods of two of these stars to be 72.5 and cause they end their lives as Type II supernovae (SN-II) 96 days. Recently, Aerts et al. (2010b) correlated the which serve as catalysts for star formation and primary sudden amplitude decline in the spectroscopically pecu- contributors to the heating and chemical enrichment liar CoRoT target HD 50064 (B1-6 Ia) to radial strange mode variability with a 37 day period; they inferred a of their host galaxies. These stars also offer a poten- ∼ tially valuable addition to the distance ladder by way of mass of 45M⊙ for this luminous mass-losing BSG. De- the BSG wind-momentum-luminosity (Kudritzki et al. spite several decades of observational efforts on ground 1999) and flux-weighted gravity-luminosity relationships based photometry and spectroscopy of α Cyg type pul- (Kudritzki et al. 2003) which give very consistent dis- sating stars, it is not certain what portion of their vari- tance determinations to the Local Group of galaxies ability is periodic, nor how far they deviate from strict (Kudritzki et al. 2008; U et al. 2009). periodicity. Fascinating features of BSGs include their microvari- The nature of the aforementioned variability in high ability in flux (e.g. Sterken 1977), complex changes in mass post-MS variables is still poorly understood. The radial velocity (RV), as well as variations in the equiv- photometric variability of HD 163899 (B2 Ib/II) found alent widths (EW) and line profiles of Balmer (espe- by MOST observations revealed a total of 48 pressure cially Hα) and metallic lines (Kaufer et al. 1996, 1997; (p) and gravity (g) modes. Saio et al. (2006) showed Richardson et al. 2011). Kaufer and coworkers moni- that these modes can be simultaneously excited by the tored the optical spectra of six BA supergiants (includ- κ-mechanism and reach the surface if they arrive at the ing Rigel and α Cyg). They analyzed the radial ve- base of the intermediate convection zone (ICZ) with an appropriate phase. Gautschy (2009) searched for the arXiv:1201.0843v1 [astro-ph.SR] 4 Jan 2012 locity measures of these stars and associated their ob- served variability with non-radial pulsations since the origin of long-period variabilities in the prototype α traveling features in the dynamical spectra could not Cyg. Interestingly his Figure 5 shows a gap where no be reconciled with the rotational periods of the stars. instabilities are predicted for evolutionary tracks with Waelkens et al. (1998) discovered 32 pulsating BSGs in 3.95 . log Teff . 4.15. Rigel lies in this gap. In con- the Hipparcos database which belong to the α Cyg class trast to this result, the study of Saio (2011) predicts that of variable stars, with periods ranging from 1.5 to 24 Rigel should be unstable against non-radial convective g− modes. This contradiction could arise from different [email protected] physical ingredients (such as rotational and overshoot 1 Department of Physics, Institute for Advanced Studies in mixing) and various numerical techniques in pulsation Basic Sciences (IASBS), Zanjan 45137-66731, Iran codes. This is not surprising since evolved massive stars 2 Department of Astronomy, Villanova University, 800 Lan- like Rigel are demanding to model. caster Ave, Villanova, PA 19085, USA The radial strange mode is proposed as another mech- 3 Royal Military College of Canada, PO Box 17000, Station Forces, Kingston, ON K7K 4B4, Canada anism to induce instability and interplay with mass 4 Center of Excellence in Information Systems, Tennessee loss efficiencies in these massive stars with log(L/M) & State University, Nashville, USA 4 log(L⊙/M⊙) (see e.g. Dziembowski & Slawinska 2005; 5 Departamento de Astrof´ısica, Centro de Astrobiolog´ıa (INTA-CSIC), PO BOX 78, 28691 Villanueva de la Ca˜nada, Saio 2011; Aerts et al. 2010b, and references therein). Madrid, Spain This requirement is also marginally fulfilled by Rigel. 2 Moravveji et al. Godart et al. (2009) investigated the destructive impact of core overshooting and mass loss during the MS life- TABLE 1 Updates to Physical Parameters of Rigel time on the extent of the ICZ, and showed that models with wider ICZs are more likely to destabilize stellar os- Parameter Value Reference cillations. These studies show that asteroseismology of ± dHip [pc] 264 24 van Leeuwen (2007) slowly pulsating blue supergiants can reveal a wealth of ± Teff [K] 12 100 150 Przybilla et al. (2010) information about the internal structure of SN-II progen- +0.07 log(L/L⊙) 5.08−0.10 This work itors. This can be complementary to the understanding Ys 0.32 ± 0.04 Przybilla et al. (2010) of the physical properties of pulsars which originate from v sin i [km s−1] 25 ± 3 Przybilla et al. (2010) the cores of massive stars (Heger et al. 2005). Sim´on-D´ıaz et al. (2010) This is the first paper in a series of investigations aim- Inclination 60◦ . i . 90◦ This work ± ing at probing the details of the internal structure of θLD [mas] 2.75 0.01 Aufdenberg et al. (2008) R/R⊙ 78.9 ± 7.4 This work BSGs through asteroseismic study, and investigating the − ± Mbol [mag] 7.92 0.28 This work −1 −7 origin of their short- and long-period variability. Rigel M˙ [M⊙ yr ] 1-2×10 Chesneau et al. (2010) was selected as an ideal test subject because of its appar- Bd [G] . 25 Shultz et al. (2011) ± ent brightness, proximity and current evolutionary sta- Pshortest [d] 1.2191 0.0001 This work ± tus. The latest measurements of the physical parame- Plongest [d] 74.74 0.28 This work ters of Rigel collected from the literature are summarized in Section 2. The space-based MOST photometry and (Kaufer et al. 1996; Morrison et al. 2008) and other ground-based spectroscopy are presented in Section 3. metal-line equivalent widths due to non-radial pulsations The results of the multimode pulsation frequencies are (Kaufer et al. 1997) and mass loss (Chesneau et al. 2010) finally discussed in Section 4. The interpretation of the is already published. pulsation frequencies will appear in Moravveji et al. (in Rigel shows variations in the Hα spectral feature with preparation). some outburst events being recorded by Israelian et al. 2. ABOUT RIGEL (1997) and Chesneau (VLTI campaign, in preparation). The morphology of this line is studied by Morrison et al. th Rigel (β Ori; HD 34085; B8 Ia; V∼0.12 mag) is the 6 (2008). Moreover, Chesneau et al. (2010) detected an brightest star in the night sky and the most luminous star extended rotating Hα region from optical interferometry in the solar neighborhood. It is a member of a multiple with VEGA/CHARA (Mourard et al. 2009) and specu- star system where its companion, Rigel B, is a spectro- ′′ late that Rigel is observed with its rotation axis along scopic binary about 9.5 distant (Sanford 1942). From the North-South direction at significantly high inclina- the revised Hipparcos parallax the distance to Rigel is tion angle. They estimated the current mass loss rate ± dHip = 264 24 pc (van Leeuwen 2007) which is smaller ˙ − × −7 −1 ± is M = 1 2 10 M⊙ yr . Assuming a high incli- than 360 40 pc adopted by Hoffleit & Jaschek (1982) nation angle, the rotation velocity is well below 20% of who assumed that Rigel originates from the τ Ori R1 the estimated critical breakup rotation rate of ∼185 km complex. s−1 (Eq. 2.19 in Maeder 2009), so spherical symmetry Fundamental parameters of Rigel which impose valu- can be safely assumed for its geometry. The current best able constraints on the equilibrium model of the star values of the physical parameters of Rigel are presented are already measured. They are effective temperature in Table 1. Additional parameters for Rigel are given by Teff = 12100 ± 150 K, surface gravity log g =1.75 ± 0.10, +0.07 Przybilla et al.