arXiv:astro-ph/0406316 v1 14 Jun 2004 ioe D228 fi[email protected] 21218; MD timore, rv,Hnll,H 96822 HI Honolulu, Drive, Spain Madrid, 29006 IL YTEGNRU IACA UPR FTHE OF SUPPORT FOUNDATION. FINANCIAL KECK GENEROUS W.M. THE POS- BY MADE SIBLE SPACE WAS OBSERVATORY AND CALI- THE AERONAUTICS OF ADMINISTRATION. NATIONAL UNIVERSITY THE AND THE FORNIA TECHNOLOGY, A IN- OF AS CALIFORNIA STITUTE OPERATED THE AMONG IS PARTNERSHIP WHICH SCIENTIFIC OBSERVATORY, KECK observational W.M. give who (2003) ∗ (2004), Figer Kroupa by & recently Weidner made and are cases promising most uvyn h aayaogtepaeadtedat of dearth masses the with and clusters plane in starburst the incompleteness resolved along Galaxy of the because surveying and stellar observation, pressure maximum radiation and between the opacity, interplay because theory, to complex both the limit escapes quantity of basic upper a Such firm mass. known no AGN. is and galaxies, prod- starburst nearby output mergers, prominent of galactic also of formation are ucts They further they and stars. ends process probes, mass lower formation likely As star presence the of the bursts. their limits in ray upper explosions are the gamma they define energetic as that most seen evidence the is Universe, of there evolution progenitors and early galaxies, the the first in medium. role the the interstellar important of an the on play likely of influence They composition their the and through control environs energetics galaxies they local individual their ingredients, of the and As the evolution when chemical to reionization and born. dynamical of Orion, were age stars as the first such during sites, Universe from early scales, formation distance star and size individual all on phenomena physical lcrncades fi[email protected] address: Electronic AAPEETDHRI EEOTIE TTHE AT OBTAINED WERE HEREIN PRESENTED DATA 2 3 1 ept hi eta oei hs aidtpc,there topics, varied these in role central their Despite astro- of probes and ingredients key are stars Massive rpittpstuigL using 2006 typeset 24, Preprint May version Draft pc eecp cec nttt,30 a atnDie B Drive, Martin San 3700 Institute, Science Telescope Space nttt o srnm,Uiest fHwi,28 Woodla 2680 Hawaii, of University Astronomy, for Institute nttt eEtutr el aei,CI,Srao121, Serrano CSIC, Materia, la de Estructura de Instituto ujc headings: Subject e ihcmoet eaae nvlct by velocity in separated 1806 components with LBV tem of spectroscopy near-infrared high-resolution n ytmms of mass system a and ihms siae ( estimates mass high oecnitn ihteae fnab tr nteLV1806 LBV the in argu We stars V nearby system. of of binary velocity ages massive the the a with is consistent it more yet known, star massive ue ymsiesasi omto n h pct finfall m of most opacity “the 1806 the of LBV between and claims interplay observational test. formation the many in and of stars uncertainties, understanding massive Our by duced stars. of masses ept uhtertcladosrainlpors,ther progress, observational and theoretical much Despite 1. A INTRODUCTION T E tl mltajv 6/22/04 v. emulateapj style X − LSR 0i atclrylmnu bet L object, luminous particularly a is 20 iais pcrsoi tr:eouin—sas early- stars: — evolution stars: — spectroscopic binaries: aaees—sas niiul(B10-0 nrrd st infrared: — (LBV1806-20) individual stars: — parameters oadF Figer F. Donald 3 ms km =36 H OBELNDSETU FLV1806 LBV OF SPECTRUM DOUBLE-LINED THE M ∼ initial 130 M > 200 − ⊙ 1 . > scnitn ihadsac f1. p,almnst f10 of luminosity a kpc, 11.8 of distance a with consistent is 10 M 1 ⊙ 4 rnic Najarro Francisco , ,bsd npr,o t iiaiyt h itlSa.W pre We Star. Pistol the to similarity its on part, in based, ), M rf eso a 4 2006 24, May version Draft ⊙ The . ABSTRACT wn al- ∼ 0k s km 70 vdnefra pe asctffof cutoff mass upper an for evidence htR3 satal asv lse fsas rather M stars, with of star cluster find single massive (1991) supermassive a al. a actually et than Weigelt is instance, such R136 of For that composed Some actually obser- stars. is subsequent object singularity. multiple with the of test, that revealing this vations test passed not crucial have claims the face must limit have. may the this star of beyond a understanding mass significantly a our maximum mass change even single-handedly a Finding would with respectively. star R136, single and cluster Arches 90 aai,Tni aegi18) Likewise, that claim 1983). (2000) al. Tarenghi et & Tanzi, Damineli Panagia, Winkler 1980; & Schmidt-Kaler, Schlosser, (Feitzinger, rsn-a assof masses present-day iia oLmnu leVrals(Bs,ie thas it i.e. characteristics (LBVs), has 1806 Variables source Blue SGR L Luminous this to repeater, that similar ray determine Neugebauer & (1995) gamma Matthews, Kulkarni, soft Kerkwijk, van the Presumably, Star. Pistol near observations. the as the massive of as time nearly is same the FMM362 the at Pis- had Star, the it as tol characteristics that photometric demonstrate and spectroscopic (2000) al. et Geballe yei h ag 9B,laigte opooethat propose to them leading O9-B2, range the in type ie ta.(99)ietf erti otePistol the to 2 near-twin than a less identify FMM362, (1999a) while spectrum. Star, composite mas), al. obviously (14 et an AU show Figer 110 not of do single latter distance is star the projected the a that reveal to speckle former Keck down the their upon spectra; based and single data is star the that show hnasnl trwt M with star single a than o tr nteQitpe lse,a n ftemost the of of mass one initial as an Cluster, with Pis- Quintuplet the known, identify the massive (1998) in al. Star, et tol Figer Center. Galactic the ∼ > − fcus,ec dsoeyo h otmsiestar” massive most the of “discovery each course, Of ukrie l 19)ietf uiossource luminous a identify (1995) al. et Kulkarni oeo h otmsiesnl tr nw r near are known stars single massive most the of Some − 10 ∼ 1 0 hwn hti spsil iaysys- binary a possibly is it that showing 20, fcret hnti ytmi o h most the not is system this then correct, If . 10 2 6 ofP Kudritzki P. Rolf , L 6 sn nw r pe ii othe to limit upper firm known no is e ⊙ siesa”hv aldtesingularity the failed have star” assive L − n aeili ujc otheoretical to subject is material ing ⊙ h mes aito rsuepro- pressure radiation immense the htabnr,o egd ytmis system merged, or binary, a that e n -adsetu niaigaspectral a indicating spectrum K-band a and 0cutr nadto,w n that find we addition, In cluster. 20 o hc oehv lie very claimed have some which for , − ye—sas fundamental stars: — type 20 ∼ ∗ ∼ > 70 3 ars 100 ′ M rmtePso tr and Star, Pistol the from ⊙ M ⊙ and η . a sbnr,with binary, is Car ∼ ∼ ∼ 30 150 6 ∼ . 250-2,000 3 200 sent L M M ⊙ ⊙ M , ⊙ rather , ⊙ nthe in and , − M 20. ⊙ 2 FIGER, NAJARRO, & KUDRITZKI it is an LBV-candidate. Eikenberry et al. (2001, 2004) estimate the properties of the star, finding a luminosity, and thus mass, that is at least comparable to that of the Pistol Star. Eikenberry et al. (2004) obtained high resolution images at Palomar, showing that the star is single down to about 900 AU. Note that the putative binary components of η Car have a maximum separation of 30 AU (Damineli et al. 2000). So, it is possible that LBV 1806−20 is multiple, but it is certainly not a “cluster,” in the traditional sense. Even the most compact cluster in the Galaxy, the Arches cluster, has a half-light radius of ∼40,000 AU (Figer et al. 1999b). In this Letter, we present high-resolution near-infrared spectra of LBV 1806−20 showing that the helium ab- sorption lines are double, suggesting that the star may Fig. 1.— Spectra of the Pistol Star and LBV 1806−20 near be binary. the 1.70 µm He I transition. The feature in the spectrum for LBV 1806−20 is double (marked by arrows), with the blueward 2. OBSERVATIONS AND DATA REDUCTION component having a P Cygni signature. The broad feature on the left at 1.688 µm is attributable to an Fe II transition. The The data were obtained at Keck on 22 June 2003, using two features at 1.691 µm and 1.698 µm are due to Si II. The NIRSPEC (McLean et al. 1998). LBV 1806−20 was ob- weak absorption features near 1.6905 µm and 1.6960 µm are due to telluric absorption. served at high resolution (R∼26,500=λ/∆λFWHM) from 1.68 µm to 2.25 µm in three grating settings, and at low resolution (R∼2,500) from 2.7 µm to 4.2 µm in two grating settings. We obtained two exposures per grat- ing setting, with the telescope nodded along the direc- tion of the slit inbetween exposures. The exposure times were 200 seconds, 100 seconds, and 30 seconds, for the high resolution short-wavelength, long-wavelength, and low resolution pointings, respectively. Quintuplet Star #3 (hereafter “Q3”), which is feature- less in this spectral region (Figer et al. 1998, Figure 1), was observed as a telluric standard (Moneti et al. 1994). Arc lamps containing Ar, Ne, Kr, and Xe, were observed to set the wavelength scale. In addition, a continuum lamp was observed through an etalon filter in order to produce an accurate wavelength scale between arc lamp lines and sky lines (predominantly from OH). A quartz- tungsten-halogen lamp was observed to provide a “flat” Fig. 2.— Spectra of massive stars near the 2.112/2.113 µm He I doublet. The Pistol Star, IRS16NE, and LBV 1806−20 are image which was divided into the background-subtracted LBVs or LBV-candidates (Figer et al. 1998; Tamblyn et al. 1996; target images. van Kerkwijk, Kulkarni, Matthews, & Neugebauer 1995), and the We subtracted background flux, dark current, and remaining stars are B-type supergiants (Figer et al. 1999a; I residual bias by differencing the nod pair frames, and we Genzel et al. 1996). All the spectra show a single He doublet near 2.1125 µm, except that for LBV 1806−20, in which it is obvi- divided the result by the flat field frame in order to com- ously double (marked by arrows). The broad feature near 2.089 µm pensate for non-uniform response. We extracted spectra is common to LBV spectra and corresponds to an Fe II transition. of the objects and applied a wavelength solution based upon the locations of arc lines. Finally, we divided the target spectra by the spectrum for Q3 in order to remove telluric effects. Note that we omitted this last step in the case of the spectra near 1.7 µm because Q3 is too faint at such short wavelengths. The final extracted spectra are shown in Figures 1, 2, 3, and 4. They have been individually shifted along the x-axis so that the spectral features are aligned with their vacuum wavelengths at rest. 3. DISCUSSION 3.1. Spectral Analysis The spectra of LBV 1806−20 and the Pistol Star are virtually indistinguishable, except for the morphology of the Fe II and He I lines. The latter are generated in the photospheres and winds of LBVs and are invariably seen Fig. 3.— Spectra of the Pistol Star and LBV 1806−20 near the − I Brackett−γ line. The two absorption lines near 2.160 µm in the in absorption. In LBV 1806 20, the He lines are obvi- − ± −1 spectrum for LBV 1806 20 (marked by arrows) are due to two ously double, with a velocity separation of 72 3 kms . separate transitions of He I. The short-wavelength component of the 1.70 µm line has LBV 1806−20 3

Fig. 4.— Spectra of the Pistol Star and LBV 1806−20 in the 3- Fig. 5.— Galactic rotation curve toward the line of sight 4 µm region. Both spectra have a strong Brackett−α emission line of LBV 1806−20 from Brand & Blitz (1993). The cross-hatched near 4.05 µm, a Pfund−γ (H I 6-5) emission line near 3.74 µm, and region corresponds to our estimate of the systemic velocity of a set of emission lines in the Humphreys series (H I n-6). Hydro- LBV 1806−20 in the LSR frame. The error band limits are genic transitions of He I can be seen in absorption just shortward ±10 km s−1, where 6 kms−1 represents the measurement error, of the Brackett−α lines. The remaining lines are mostly due to and 4 km s−1 allows for mass ratios between 1:1 and 1:1.5. Fe II.

for that star. a P Cygni signature, indicating that the star responsible Spectra of both stars have similar morphology in the for the bulk of the emission lines is blue-shifted with re- H lines. However, the presence of deeper absorption fea- spect to the other star. Indeed, none of the emission tures in the observed H lines of the LBV 1806−20, as well lines are double, indicating that the red-shifted compo- as their narrower nature, are consistent with underlying nent likely has less luminosity and/or a thinner wind. absorption H lines from the putative companion. The velocity separation between the two minima in the We estimate the line-of-sight velocity for components is 106 kms−1, larger than the separation of LBV 1806−20 by comparing the vacuum wavelength the components near 2.112 µm, as would be expected of the He I line near 2.1125787 µm with its measured given the P Cygni nature of the profile in the blue com- positions at 2.1125065 µm and 2.1130192 µm, as ponent of the 1.700 µm line. The He I absorption com- determined by fitting Gaussian profiles to the two ponent near 2.160 µm appears to be double (marked by absorption components. We adjust the velocity to arrows), with a separation of 99 kms−1; however, we the heliocentric frame by adding +0.2 km s−1 to the believe the lines correspond to two expected He I lines measured velocity to compensate for Earth motion at with velocity seperation of 83 kms−1 near that wave- the time of observations. We add +10 ± 0.4 kms−1 length. Figure 4 does not appear to show any double to this value in order to adjust the velocity to the absorption lines in the spectrum of LBV 1806−20, pre- local standard of rest frame (Dehnen & Binney 1998). sumably due to the relatively low spectral resolution of − Assuming a roughly equal mass binary, we average the the data (70 km s 1 corresponds roughly to 2 pixels in velocities of the two absorption lines, to find a systemic the low resolution mode). The absorption feature near velocity of V =36.4 kms−1; note that a change in this 4.04 µm is likely due to a blend composed of a He I LSR assumption, to a mass ratio of 1:1.5, would only induce doublet (5f3F -4d3D at 4.0377 µm and 5f1F -4d1D at o o a change in our estimate of ±4 kms−1. Our estimate is 4.0409 µm). The He I absorption blend near 3.703 µm similar to those we obtain for three nearby hot stars, for (5d3D-4p3P ) has a FWHM of 183 kms−1, compared to o which V =35.0±5.7 kms−1 (see Figure 2). These a FWHM of 152 kms−1 for the nearby Pfund-γ line at ave,LSR velocities contrast with the value in Eikenberry et al. 3.74 µm. (2004) of 10±20 kms−1. The Fe II lines are broad in the spectra of both the The bulk velocity of the LBV 1806−20 cluster can Pistol star and LBV 1806−20, but they have multiple be interpreted in the context of the Galactic rotation peaks in the latter. The 1.688 µm and 2.089 µm lines curve in Brand & Blitz (1993) in order to estimate its have nearly identical morphology, each with three peaks distance. We find a distance of 4.1 kpc or 11.8 kpc in similar intensity ratios and relative positions in the (see Figure 5). Only the high solution is consistent spectrum of LBV 1806−20. The Gaussian shape of the with a location beyond an NH cloud at a distance of Fe II lines in the spectrum for the Pistol star is likely 3 5.7 kpc (Brand & Blitz 1993). This estimate is consid- caused by recombination closer to photosphere, where erably less than the estimate in Eikenberry et al. (2004) the wind velocity has not reached its constant terminal +1.8 value. In the spectrum for LBV 1806−20, the iron likely of 15.1−1.3 kpc. recombines further out in the wind so that emission is 3.2. − dominated by gas at the terminal velocity, thus causing Mass of LBV 1806 20 the flat top shape of the emission line. This explanation For the far distance estimate (d∼11.8 kpc), and for the variation in the Fe II line shapes in the two stars the photometry in Eikenberry et al. (2004), we find suggests that LBV 1806−20 has slightly higher ionizing MK =−9.5±0.1. It is difficult to convert this value into radiation, consistent with the apparent P-Cygni nature a luminosity estimate because the range of tempera- of the blue component of the He I lines in the spectrum tures, and, therefore the K-band bolometric correction, 4 FIGER, NAJARRO, & KUDRITZKI varies widely for LBVs (Blum, Depoy, & Sellgren 1995). tuplet cluster should not live longer than ∼3 Myr old, However, if we use the BCK for the Pistol Star from whereas the other cluster members are 4 Myr old; how- Figer et al. (1998), then we arrive at a luminosity of ever, in this case, the original natal cloud must have been 6.3 10 L⊙. Using Langer’s stellar evolution models in strongly sheared by the tidal field in the Galactic center Figer et al. (1998), we find an initial mass ∼130 M⊙ in on time scales much less than 2 Myr. A more general the case that the object is a single star. explanation for the age discrepancy in both the Quintu- plet and in the LBV 1806−20 cluster might be that the 3.3. LBV 1806−20 as a Massive Binary most massive stars are binary systems of less massive We suggest that the object is a massive binary, based stars (which can live longer), or that they were binary upon the double lines in its spectrum. For a 1:1 mass systems that recently merged (Figer & Kim 2002). system, the components have initial masses of 65 M⊙, This proposition also resolves the problem of explain- assuming a mass-luminosity relation of M∝Lα, where ing why 200 M⊙ stars are observed in the Quintuplet alpha∼1. Our present analysis is fairly insensitive to the cluster, whereas the upper mass cutoff to the initial mass ∼ mass ratio, a topic that we will revisit in a future study. function in the nearby Arches cluster is 150 M⊙ (Figer ∼ By comparison, the initial masses of the components in 2003). In fact, if single stars with masses of 100 M⊙ η Car are 120 M⊙ and 40 M⊙ (Damineli, Conti, & Lopes exist, then where are the binary systems with such stars? 1997; Pittard & Corcoran 2002). Another massive bi- Certainly, some of the most massive stars exist in binary nary system, WR20a, contains roughly equal mass com- systems, as evidenced by the substantial binary fraction ponents, each having M∼70-80 M⊙ (Rauw et al. 2004; ( 25%) in the WR catalogue (van der Hucht 2001). Per- Bonanos et al. 2004). haps we are witnessing the ancestors of such systems in LBV 1806−20 is similar, in appearance, to a growing these LBVs. − group of objects, i.e. the Pistol Star, FMM362, η Car, The evidence in this paper suggests that LBV 1806 20 and others. Whether these objects are single or multiple, may be binary, but we will continue to explore the they all presumably have at least one evolved component possibility that the double-lined spectra could have an alternate explanation. For instance, perhaps the that is massive (∼>50 M⊙) and some are surrounded by circumstellar material produced by massive eruptions. profile is a composite of an emission feature, gen- erated by a disk or wind, at the systemic veloc- 3.4. The Formation of the Most Massive Stars ity and superposed on a photospheric absorption fea- ture. In addition, we will refine the temperature, lu- Theoretical arguments suggest that stars with ini- minosity, and mass ratio estimates through a more de- tial masses ∼>150 M⊙ have hydrogen-burning (main se- tailed analysis, coupled with wind/atmosphere mod- quence) lifetimes of about 2 Myr (Bond, Arnett, & Carr els Najarro, Figer & Kudritzki (2004). In addition, we 1984), and total lifetimes no more than ∼3 Myr have begun a spectroscopic monitoring program to de- (Meynet et al. 1994). If LBV 1806−20 is single, and it tect changes that might indicate a binary nature for had an initial mass in this range, then it would necessar- LBV 1806−20. ily be much younger than nearby Wolf-Rayet stars in the same cluster that are at least 3 Myr old (Meynet et al. 1994). Eikenberry et al. (2004) suggest that a supernova from a first generation of stars could have triggered the F. N. acknowledges grants AYA2003-02785-E and formation of LBV 1806−20, and perhaps other nearby ESP2002-01627. We thank Keck staff members Randy stars. While this mechanism is often invoked to explain Campbell and Grant Hill. We acknowledge useful dis- mutli-generation star formation, it is not the only expla- cussions with Richard Larson, Stephen Eikenberry, and nation in this case. Note that the LBVs in the Quin- Augusto Damineli.

REFERENCES Blum, R. D., Depoy, D. L., & Sellgren, K. 1995, ApJ, 441, 603 Figer, D. F., Najarro, F., Morris, M., McLean, I. S., Geballe, T. Bonanos, A. Z., et al. 2004, ApJ, submitted, ArXiv Astrophysics, R., Ghez, A. M., & Langer, N. 1998, ApJ, 506, 384 astro-ph/0405338 Geballe, T.R., Figer, D.F., & Najarro, F. 2000, ApJ, 530, 97 Bond, J. R., Arnett, W. D., & Carr, B. J. 1984, ApJ, 280, 825 Genzel, R., Thatte, N., Krabbe, A., Kroker, H., & Tacconi-Garman, Brand, J. & Blitz, L. 1993, A&A, 275, 67 L. E. 1996, ApJ, 472, 153 Damineli, A., Kaufer, A., Wolf, B., Stahl, O., Lopes, D. F., & de Kulkarni, S. R., Matthews, K., Neugebauer, G., Reid, I. N., van Ara´ujo, F. X. 2000, ApJ, 528, L101 Kerkwijk, M. H., & Vasisht, G. 1995, ApJ, 440, L61 Damineli, A., Conti, P. S., & Lopes, D. F. 1997, New Astronomy, McLean et al. 1998, SPIE Vol. 3354, 566 2, 107 Meynet, G., Maeder, A., Schaller, G., Schaerer, D., & Charbonnel, Dehnen, W. & Binney, J. J. 1998, MNRAS, 298, 387 C. 1994, A&A Supp., 103, 97 Eikenberry, S. S., Garske, M. A., Hu, D., Jackson, M. A., Patel, Moneti, A., Glass, I.S. & Moorwood, A.F.M. 1994, MNRAS, 268, S. G., Barry, D. J., Colonno, M. R., & Houck, J. R. 2001, ApJ, 194 563, L133 Najarro, F., Figer, D. F. & Kudritzki, R. P. 2004, in preparation Eikenberry, S. S. et al. 2004, ApJ, submitted, ArXiv Astrophysics Panagia, N., Tanzi, E. G., & Tarenghi, M. 1983, ApJ, 272, 123 e-prints, astro-ph/0404435 Pittard, J. M. & Corcoran, M. F. 2002, A&A, 383, 636 Feitzinger, J. V., Schlosser, W., Schmidt-Kaler, T., & Winkler, C. Rauw, G., et al. 2004, A&A, accepted, ArXiv Astrophysics, 1980, A&A, 84, 50 astro-ph/0404551 Figer, D. F. 2003, IAU Symposium, 212, 487 Tamblyn, P., Reike, G., Hanson, M., Close, L., McCarthy, D., Figer, D. F. & Kim, S. S. 2002, ASP Conf. Ser. 263, 287 Reike, M. 1996, ApJ, 456, 206 Figer, D. F., Kim, S. S., Morris, M., Serabyn, E., Rich, R. M., & van der Hucht, K. A. 2001, New Astronomy Review, 45, 135 McLean, I. S. 1999a, ApJ, 525, 750 van Kerkwijk, M. H., Kulkarni, S. R., Matthews, K., & Neugebauer, Figer, D. F., McLean, I. S., & Morris, M. 1999, ApJ, 514, 202 G. 1995, ApJ, 444, L33 Weidner, C. & Kroupa, P. 2004, MNRAS, 348, 187 LBV 1806−20 5

Weigelt, G. et al. 1991, ApJ, 378, L21