Epoch J2000.0 Equinox J2000.0
|Right ascension||5h 38m 42.39s|
|Declination||−69° 06′ 02.91″|
|Apparent magnitude (V)||12.23|
|Evolutionary stage||Wolf–Rayet star|
|B−V color index||0.03|
|Absolute magnitude (MV)||−8.09|
|Luminosity (visual, LV)||150,000 L☉|
|Surface gravity (log g)||4.0 cgs|
|Temperature||53,000 ± 3000 K|
|Age||0 - 0.8 Myr|
RMC 136a1 (usually abbreviated to R136a1) is a Wolf–Rayet star located at the center of R136, the central condensation of stars of the large NGC 2070 open cluster in the Tarantula Nebula. It lies at a distance of about 50 kiloparsecs (163,000 light-years) in a neighbouring galaxy known as the Large Magellanic Cloud. It has the highest mass and luminosity of any known star, at 315 M☉ and 8.7 million L☉, and is also one of the hottest at around 53,000 K.
In 1960, a group of astronomers working at the Radcliffe Observatory in Pretoria made systematic measurements of the brightness and spectra of bright stars in the Large Magellanic Cloud. Among the objects cataloged was RMC 136, (Radcliffe observatory Magellanic Cloud catalog number 136) the central "star" of the Tarantula Nebula, which the observers concluded was probably a multiple star system. Subsequent observations showed that R136 was located in the middle of a giant region of ionized interstellar hydrogen known as an H II region that was a center of intense star formation in the immediate vicinity of the observed stars.
In 1979, ESO's 3.6 m telescope was used to resolve R136 into three components; R136a, R136b, and R136c. The exact nature of R136a was unclear and a subject of intense discussion. Estimates that the brightness of the central region would require as many as 100 hot O class stars within half a parsec at the centre of the cluster led to speculation that a star 3,000 times the mass of the Sun was the more likely explanation.
The first demonstration that R136a was a star cluster was provided by Weigelt and Beier in 1985. Using the speckle interferometry technique, R136a was shown to be made up of 8 stars within 1 arcsecond at the centre of the cluster, with R136a1 being the brightest.
Final confirmation of the nature of R136a came after the launch of the Hubble Space Telescope. Its Wide Field and Planetary Camera (WFPC) resolved R136a into at least 12 components and showed that R136 contained over 200 highly luminous stars. The more advanced WFPC2 allowed the study of 46 massive luminous stars within half a parsec of R136a and over 3,000 stars within a 4.7 parsec radius.
In the night sky, R136 appears as a 10th magnitude object at the core of the NGC 2070 cluster embedded in the Tarantula Nebula in the Large Magellanic Cloud. It required a 3.6 metre telescope to detect R136a as a component of R136 in 1979, and resolving R136a to detect R136a1 requires a space telescope or sophisticated techniques such as adaptive optics or speckle interferometry.
South of about the 20th parallel south, the LMC is circumpolar, meaning that it can be seen (at least in part) all night every night of the year, weather permitting. In the Northern Hemisphere, it can be visible south of the 20th parallel north. This excludes North America (except southern Mexico), Europe, northern Africa and northern Asia.
The R136a system at the core of R136 is a dense luminous knot of stars containing at least 12 stars, the most prominent being R136a1, R136a2, and R136a3, all of which are extremely luminous and massive WN5h stars. R136a1 is separated from R136a2, the second brightest star in the cluster, by 5,000 AU.
R136 is located approximately 157,000 light-years away from Earth in the Large Magellanic Cloud, positioned on the south-east corner of the galaxy at the centre of the Tarantula Nebula, also known as 30 Doradus. R136 itself is just the central condensation of the much larger NGC 2070 open cluster.
For such a distant star, R136a1 is relatively unobscured by interstellar dust. The reddening causes the visual brightness to be reduced by about 1.8 magnitudes, but only around 0.22 magnitudes in the near infrared.
The distance to R136a1 cannot be determined directly, but is assumed to be at the same distance as the Large Magellanic Cloud at around 50 kiloparsecs.
Although binary systems are very common among the most massive stars, R136a1 appears to be a single star as no evidence of a massive companion has been detected.
X-ray emission was detected from R136 using the Chandra X-ray Observatory. R136a and R136c were both clearly detected, but R136a could not be resolved. Another study separated the R136a1/2 pair from R136a3. R136a1/2 showed relatively soft x-rays not thought to indicate a colliding winds binary.
Rapid Doppler radial velocity variations would be expected from a pair of equal mass stars in a close orbit, but this has not been seen in the R136a1 spectrum. A high orbital inclination, a more distant binary, or a chance alignment of two distant stars cannot be completely ruled out but is thought to be unlikely. Highly unequal binary components are possible, but would not affect the modelling of R136a1's properties.
R136a1 is a high-luminosity WN5h star, placing it on the extreme top left corner of the Hertzsprung-Russell diagram. A Wolf–Rayet star is distinguished by the strong, broad emission lines in its spectrum. This includes ionized nitrogen, helium, carbon, oxygen and occasionally silicon, but with hydrogen lines usually weak or absent. A WN5 star is classified on the basis of ionised helium emission being considerably stronger than the neutral helium lines, and having roughly equal emission strength from NIII, NIV, and NV. The "h" in the spectral type indicates significant hydrogen emission in the spectrum, and hydrogen is calculated to make up 40% of the surface abundance by mass.
WNh stars as a class are massive luminous stars still burning hydrogen at their cores. The emission spectrum is produced in a powerful dense stellar wind, and the enhanced levels of helium and nitrogen arise from convectional mixing of CNO cycle products to the surface.
An evolutionary mass of 265 M☉ is found from near infrared (K-band) spectra using a combination of non-LTE line-blanketed CMFGEN and TLUSTY model atmosphere code. The models were validated against the dynamical masses derived for the WN6h binary NGC 3603-A1. In the worst case of a line of sight alignment or an unexpected wide binary, the stars would each be 150 M☉. R136a1 closely matches the expected properties for an initially rapidly rotating 320 M☉ star with LMC metallicity after shining for about 1.7 million years.
A current mass of 256 M☉ is found in similar analysis using PoWR (Potsdam Wolf Rayet) atmospheric models with optical and ultraviolet spectra and a mass-luminosity relation, assuming a single star.
R136a1 is undergoing extreme mass loss through a stellar wind reaching a velocity of 2,600 ± 150 km/s. This is caused by intense electromagnetic radiation from the very hot photosphere accelerating material away from the surface more strongly than gravity can retain it. Mass loss is largest for high luminosity stars with low surface gravity and enhanced levels of heavy elements in the photosphere. R136a1 loses 5.1 × 10−5 M☉ (3.21 × 1018 kg/s) per year, over a billion times more than the Sun, and is expected to have shed over 50 M☉ since its formation.
At around 8,000,000 L☉, R136a1 is the most luminous star known, radiating more energy in five seconds than the Sun does in a year. If it replaced the Sun in the Solar System, it would outshine the Sun by 94,000 times (MV = −7.6) and would appear from Earth at magnitude -39. Its brightness at a distance of 10 parsecs, the absolute visual magnitude, is -7.6, three magnitudes brighter than Venus ever appears from Earth. Its brightness at the distance of the nearest star to Earth, Proxima Centauri (just over a parsec), would be about the same as the full Moon.
R136a1 supplies ~7% of the ionizing flux of the entire 30 Doradus region, as much as 70 O7 dwarf stars. Along with R136a2, a3, and c, it produces 43-46% of the Lyman continuum radiation of the whole R136 cluster.
Massive stars lie close to the Eddington limit, the luminosity at which the radiation pressure acting outwards at the surface of the star equals the force of the star's gravity pulling it inward. Above the Eddington limit, a star generates so much energy that its outer layers are rapidly thrown off. This effectively restricts stars from shining at higher luminosities for long periods. The classical Eddington luminosity limit is not applicable to stars such as R136a1 that are not in hydrostatic equilibrium, and its calculation is extremely complex for real stars. The empirical Humphrey-Davidson limit has been identified as a luminosity limit for observed stars, but recent models have attempted to calculate useful theoretical Eddington limits applicable to massive stars. R136a1 is currently around 70% of its Eddington luminosity.
R136a1 has a B-V index of about 0.03, which is a typical colour for an F-type star. The "U-V" colour from the HST WFPC2 336 nm and 555 nm filters is −1.28, more indicative of an extremely hot star. This variation of different colour indices relative to a blackbody is the result of interstellar dust causing reddening and extinction. The reddening (EB-V) can be used to estimate the level of visual extinction (AV). EB-V values of 0.29 - 0.37 have been measured, with considerable uncertainty due to contamination from close neighbours such as R136a2 0.1" away, leading to AV around 1.80 and a de-reddened B-V (B-V0) of -0.30.
The effective temperature of a star can be approximated from the colour, but this is not very accurate and spectral fitting to an atmospheric model is necessary to derive the temperature. Temperatures of 53,000 K - 56,000 K are found for R136a1 using different atmospheric models. Older models had produced temperatures around 45,000 K and hence dramatically lower luminosities. The extreme temperature of the star causes its peak radiation to be around 50 nm and nearly 99% of the radiation to be emitted outside the visual range (a bolometric correction around −5).
R136a1 does not have a well-defined visible surface like the Earth or the Sun. The hydrostatic main body of the star is surrounded by a dense atmosphere being accelerated outwards into the stellar wind. An arbitrary point within this wind is defined as the surface for measuring the radius, and different authors may use different definitions. For example, a Rosseland optical depth of 2/3 corresponds approximately to a visible surface while a Rosseland depth of 20 or 100 corresponds more closely to a physical photosphere. Stellar temperatures are typically quoted at the same depth so that the radius and temperature correspond to the luminosity.
R136a1's dimensions are far smaller than the largest stars: red supergiants are several hundred to over a thousand R☉, tens of times larger than R136a1. Despite the large mass and modest dimensions, R136a1 has an average density around 1% of the Sun's; at about 14 kg/m3, it is over 10 times denser than Earth's atmosphere at sea level.
The rotation rate of R136a1 cannot be measured directly since the photosphere is hidden by the dense stellar wind and the photospheric absorption lines used to measure rotational doppler broadening are not present in the spectrum. A NV emission line at 2.1 µm is produced relatively deep in the wind and can be used to estimate rotation. In R136a1 it has a FWHM of about 15Å, indicating a slow or non rotating star, although it could be aligned with its pole facing Earth. R136a2 and a3 are rotating rapidly and the closest evolutionary models for R136a1 match a star still rotating with an equatorial speed of ∼200 km/s after ∼1.75 Myr.
R136a1 is currently fusing hydrogen to helium, predominantly by the CNO cycle due to the high temperatures at the core. Despite the Wolf–Rayet spectral appearance, it is a young star. The emission spectrum is created by a dense stellar wind caused by the extreme luminosity, with the enhanced levels of helium and nitrogen being mixed from the core to the surface by strong convection. It is effectively a main sequence star. Over 90% of the star is convective, with a small non-convective layer at the surface.
Models of star formation by accretion from molecular clouds predict an upper limit to the mass a star can achieve before its radiation prevents further accretion. The most simplistic accretion models at population I metallicities predict a limit as low as 40 M☉, but more complex theories allow masses several times higher. An empirical limit of around 150 M☉ has become widely accepted. R136a1 clearly exceeds all these limits, leading to development of new single star accretion models potentially removing the upper limit, and the potential for massive star formation by stellar mergers.
As a single star formed from accretion, the properties of such a massive star are still uncertain. Synthetic spectra indicate that it would never have a main sequence luminosity class (V), or even a normal O type spectrum. The high luminosity, proximity to the Eddington limit, and strong stellar wind, would be likely to create an If* or WNh spectrum as soon as R136a1 became visible as a star. Helium and nitrogen are rapidly mixed to the surface due to the large convective core and high mass loss, and their presence in the stellar wind creates the characteristic Wolf Rayet emission spectrum. The ZAMS at very high masses curves back to cooler temperatures, and at LMC metallicities the maximum temperature predicted to be around 56,000 K for 150-200 M☉ stars, so R136a1 would have been slightly cooler than some less massive main sequence stars.
During core hydrogen burning, the helium fraction in the core increases and according to the virial theorem the core pressure and temperature will increase. This leads to an increase in luminosity so that R136a1 is somewhat more luminous now than when it first formed. The temperature decreases slightly, but the outer layers of the star have inflated, driving even higher mass loss.
The future development of R136a1 is uncertain, and there are no comparable stars to confirm predictions. The evolution of massive stars depends critically on the amount of mass they can lose, and various models give different results, none of which entirely match observations. It is thought that WNh stars develop into LBVs as hydrogen at the core starts to become depleted. This is an important phase of extreme mass loss which enables a star, at near solar metallicity, to transition to a hydrogen-free Wolf Rayet star. Stars with sufficiently strong mixing from the core to the surface, due to the very large convective core, high metallicity, or additional rotational mixing, may skip the LBV phase and evolve directly from the hydrogen-rich WNh phase to the hydrogen-poor WN phase. Hydrogen fusion lasts for a little over two million years, and the star's mass at the end is expected to be 70-80 M☉. A single star with LMC metallicity, even if it starts out rotating very rapidly, will be braked to near zero rotation by the end of hydrogen burning.
After core helium fusion starts, the remaining hydrogen in the atmosphere is rapidly lost and R136a1 will quickly contract to a hydrogen-free WNE star and the luminosity will decrease. Wolf Rayet stars at this point are mostly helium and they lie on the Zero Age Helium Main Sequence (He-ZAMS), analogous to and parallel to the hydrogen-burning main sequence but at hotter temperatures.
During helium burning, carbon and oxygen will accumulate in the core and heavy mass loss continues. This eventually leads to the development of a WC spectrum, although at LMC metallicity a star is expected to spend most of the helium burning phase with a WN spectrum. Towards the end of helium burning, core temperature increase and mass loss cause an increase in both luminosity and temperature, with the spectral type becoming WO. Several hundred thousand years will be spent fusing helium, but the final stages of heavier element burning take no more than a few thousand years. R136a1 will eventually shrink to a little over 50 M☉, with just half a M☉ of helium left surrounding the core.
Any star which produces a carbon-oxygen (C-O) core more massive than the maximum for a white dwarf (~1.4 M☉) will inevitably suffer core collapse. This usually happens when an iron core has been produced and fusion can no longer produce the energy required to prevent core collapse, although it can happen in other circumstances.
A C-O core between about 64 M☉ and 133 M☉ will become so hot that the gamma radiation will spontaneously produce electron-positron pairs and the sudden loss of energy in the core will cause it to collapse as a pair-instability supernova (PISN), sometimes called a pair-creation supernova (PCSN). A PISN is usually only produced in very low metallicity stars that do not lose sufficient mass to keep their C-O core sizes below 64 M☉. This can also occur at LMC metallicity for very massive stars, but the predicted C-O core size for R136a1 is below 50 M☉ so a PISN is unlikely.
Iron core collapse may produce a supernova explosion, and sometimes a gamma-ray burst (GRB). The type of any supernova explosion will be a type I since the star has no hydrogen, type Ic since it has almost no helium. Particularly massive iron cores may collapse the entire star into a black hole with no visible explosion, or a sub-luminous supernova as radioactive 56Ni falls back onto the black hole. Other models predict that such a large core will produce very large amounts of 56Ni, powering a super-luminous supernova.
A type Ic supernova can produce a GRB if the star is rotating and has an appropriate mass. R136a1 is expected to lose almost all its spin long before core collapse so a GRB is unlikely.
The remnant from a type Ic core collapse supernova is either a neutron star or black hole, depending on the mass of the progenitor core. R136a1 will have a core far above the maximum mass of a neutron star, so a black hole is inevitable.
- Doran, E. I.; Crowther, P. A.; de Koter, A.; Evans, C. J.; McEvoy, C.; Walborn, N. R.; Bastian, N.; Bestenlehner, J. M.; Gräfener, G.; Herrero, A.; Kohler, K.; Maiz Apellaniz, J.; Najarro, F.; Puls, J.; Sana, H.; Schneider, F. R. N.; Taylor, W. D.; van Loon, J. Th.; Vink, J. S. (2013). "The VLT-FLAMES Tarantula Survey - XI. A census of the hot luminous stars and their feedback in 30 Doradus". Astronomy & Astrophysics. 558: A134. Bibcode:2013A&A...558A.134D. arXiv: . doi:10.1051/0004-6361/201321824.
- Hainich, R.; Rühling, U.; Todt, H.; Oskinova, L. M.; Liermann, A.; Gräfener, G.; Foellmi, C.; Schnurr, O.; Hamann, W. -R. (2014). "The Wolf–Rayet stars in the Large Magellanic Cloud". Astronomy & Astrophysics. 565: A27. Bibcode:2014A&A...565A..27H. arXiv: . doi:10.1051/0004-6361/201322696.
- Pietrzyński, G; D. Graczyk; W. Gieren; I. B. Thompson; B. Pilecki; A. Udalski; I. Soszyński; et al. (7 March 2013). "An eclipsing-binary distance to the Large Magellanic Cloud accurate to two per cent". Nature. 495 (7439): 76–79. Bibcode:2013Natur.495...76P. PMID 23467166. arXiv: . doi:10.1038/nature11878.
- Crowther, Paul A.; Caballero-Nieves, S. M.; Bostroem, K. A.; Maíz Apellániz, J.; Schneider, F. R. N.; Walborn, N. R.; Angus, C. R.; Brott, I.; Bonanos, A.; De Koter, A.; De Mink, S. E.; Evans, C. J.; Gräfener, G.; Herrero, A.; Howarth, I. D.; Langer, N.; Lennon, D. J.; Puls, J.; Sana, H.; Vink, J. S. (2016). "The R136 star cluster dissected with Hubble Space Telescope/STIS. I. Far-ultraviolet spectroscopic census and the origin of He II λ1640 in young star clusters". Monthly Notices of the Royal Astronomical Society. 458: 624. Bibcode:2016MNRAS.458..624C. arXiv: . doi:10.1093/mnras/stw273.
- Crowther, P. A.; Schnurr, O.; Hirschi, R.; Yusof, N.; Parker, R. J.; Goodwin, S. P.; Kassim, H. A. (2010). "The R136 star cluster hosts several stars whose individual masses greatly exceed the accepted 150 M⊙ stellar mass limit". Monthly Notices of the Royal Astronomical Society. 408 (2): 731. Bibcode:2010MNRAS.408..731C. arXiv: . doi:10.1111/j.1365-2966.2010.17167.x.
- Feast, M. W.; Thackeray, A. D.; Wesselink, A. J. (1960). "The brightest stars in the Magellanic Clouds". Monthly Notices of the Royal Astronomical Society. 121 (4): 337. Bibcode:1960MNRAS.121..337F. doi:10.1093/mnras/121.4.337.
- Feitzinger, J. V.; Schlosser, W.; Schmidt-Kaler, T; Winkler, C. (April 1980). "The central object R 136 in the gas nebula 30 Doradus - Structure, color, mass and excitation parameter". Astronomy and Astrophysics. 84 (1–2): 50–59. Bibcode:1980A&A....84...50F.
- Ebbets, D. C.; Conti, P. S. (1982). "The optical spectrum of R136a - The central object of the 30 Doradus nebula". The Astrophysical Journal. 263: 108. Bibcode:1982ApJ...263..108E. ISSN 0004-637X. doi:10.1086/160485.
- Weigelt, G.; Baier, G. (1985). "R136a in the 30 Doradus nebula resolved by holographic speckle interferometry". Astronomy and Astrophysics. 150: L18. Bibcode:1985A&A...150L..18W.
- Campbell, Bel; Hunter, Deidre A.; Holtzman, Jon A.; Lauer, Tod R.; Shayer, Edward J.; Code, Arthur; Faber, S. M.; Groth, Edward J.; Light, Robert M.; Lynds, Roger; O'Neil, Earl J., Jr.; Westphal, James A. (1992). "Hubble Space Telescope Planetary Camera images of R136". The Astronomical Journal. 104: 1721. Bibcode:1992AJ....104.1721C. doi:10.1086/116355.
- Hunter, Deidre A.; Shaya, Edward J.; Holtzman, Jon A.; Light, Robert M.; O'Neil, Earl J., Jr.; Lynds, Roger (1995). "The Intermediate Stellar Mass Population in R136 Determined from Hubble Space Telescope Planetary Camera 2 Images". The Astrophysical Journal. 448: 179. Bibcode:1995ApJ...448..179H. doi:10.1086/175950.
- Westerlund, B. E.; Smith, L. F. (1964). "Worlf-Rayet stars in the Large Magellanic Cloud". Monthly Notices of the Royal Astronomical Society. 128 (4): 311. Bibcode:1964MNRAS.128..311W. doi:10.1093/mnras/128.4.311.
- "Large Magellanic Cloud is spectacular from Earth’s Southern Hemisphere". December 26, 2014.
- Massey, P.; Hunter, D. A. (1998). "Star Formation in R136: A Cluster of O3 Stars Revealed by Hubble Space Telescope Spectroscopy". The Astrophysical Journal. 493: 180. Bibcode:1998ApJ...493..180M. doi:10.1086/305126.
- Bestenlehner, J. M.; Vink, J. S.; Gräfener, G.; Najarro, F.; Evans, C. J.; Bastian, N.; Bonanos, A. Z.; Bressert, E.; Crowther, P. A.; Doran, E.; Friedrich, K.; Hénault-Brunet, V.; Herrero, A.; de Koter, A.; Langer, N.; Lennon, D. J.; Maíz Apellániz, J.; Sana, H.; Soszynski, I.; Taylor, W. D. (2011). "The VLT-FLAMES Tarantula Survey". Astronomy & Astrophysics. 530: L14. Bibcode:2011A&A...530L..14B. arXiv: . doi:10.1051/0004-6361/201117043.
- Guerrero, Martín A.; Chu, You‐Hua (2008). "An X‐Ray Survey of Wolf‐Rayet Stars in the Magellanic Clouds. I. TheChandraACIS Data Set". The Astrophysical Journal Supplement Series. 177 (1): 216–237. Bibcode:2008ApJS..177..216G. arXiv: . doi:10.1086/587059.
- Townsley, Leisa K.; Broos, Patrick S.; Feigelson, Eric D.; Garmire, Gordon P.; Getman, Konstantin V. (2006). "AChandraACIS Study of 30 Doradus. II. X-Ray Point Sources in the Massive Star Cluster R136 and Beyond". The Astronomical Journal. 131 (4): 2164–2184. Bibcode:2006AJ....131.2164T. arXiv: . doi:10.1086/500535.
- Smith, Nathan; Conti, Peter S. (2008). "On the Role of the WNH Phase in the Evolution of Very Massive Stars: Enabling the LBV Instability with Feedback". The Astrophysical Journal. 679 (2): 1467. Bibcode:2008ApJ...679.1467S. arXiv: . doi:10.1086/586885.
- Hillier, D. John; Miller, D. L. (1998). "The Treatment of Non‐LTE Line Blanketing in Spherically Expanding Outflows". The Astrophysical Journal. 496 (1): 407–427. Bibcode:1998ApJ...496..407H. ISSN 0004-637X. doi:10.1086/305350.
- Lanz, Thierry; Hubeny, Ivan (2003). "A Grid of Non‐LTE Line‐blanketed Model Atmospheres of O‐Type Stars". The Astrophysical Journal Supplement Series. 146 (2): 417–441. Bibcode:2003ApJS..146..417L. arXiv: . doi:10.1086/374373.
- Hamann, W.-R.; Gräfener, G. (2004). "Grids of model spectra for WN stars, ready for use". Astronomy and Astrophysics. 427 (2): 697–704. Bibcode:2004A&A...427..697H. doi:10.1051/0004-6361:20040506.
- Gräfener, G.; Vink, J. S.; de Koter, A.; Langer, N. (2011). "The Eddington factor as the key to understand the winds of the most massive stars". Astronomy & Astrophysics. 535: A56. Bibcode:2011A&A...535A..56G. arXiv: . doi:10.1051/0004-6361/201116701.
- Schneider, F. R. N.; Langer, N.; De Koter, A.; Brott, I.; Izzard, R. G.; Lau, H. H. B. (2014). "Bonnsai: A Bayesian tool for comparing stars with stellar evolution models". Astronomy & Astrophysics. 570: A66. Bibcode:2014A&A...570A..66S. arXiv: . doi:10.1051/0004-6361/201424286.
- Breysacher, J.; Azzopardi, M.; Testor, G. (1999). "The fourth catalogue of Population I Wolf-Rayet stars in the Large Magellanic Cloud". Astronomy and Astrophysics. Astronomy and Astrophysics Supplement Series. 137 (1): 117–145. Bibcode:1999A&AS..137..117B. doi:10.1051/aas:1999240.
- A. J. van Marle; S. P. Owocki; N. J. Shaviv (2008). "Continuum driven winds from super-Eddington stars. A tale of two limits". AIP Conference Proceedings. AIP Conference Proceedings. 990: 250–253. Bibcode:2008AIPC..990..250V. arXiv: . doi:10.1063/1.2905555.
- Martins, Fabrice (2015). "Empirical Properties of Very Massive Stars". Very Massive Stars in the Local Universe. Very Massive Stars in the Local Universe. Astrophysics and Space Science Library. 412. pp. 9–42. Bibcode:2015ASSL..412....9M. ISBN 978-3-319-09595-0. arXiv: . doi:10.1007/978-3-319-09596-7_2.
- Humphreys, Roberta M.; Davidson, Kris (1994). "The luminous blue variables: Astrophysical geysers". Publications of the Astronomical Society of the Pacific. 106: 1025. Bibcode:1994PASP..106.1025H. doi:10.1086/133478.
- Köhler, K.; Langer, N.; de Koter, A.; de Mink, S. E.; Crowther, P. A.; Evans, C. J.; Gräfener, G.; Sana, H.; Sanyal, D.; Schneider, F. R. N.; Vink, J. S. (2014). "The evolution of rotating very massive stars with LMC composition". Astronomy & Astrophysics. 573: A71. Bibcode:2015A&A...573A..71K. arXiv: . doi:10.1051/0004-6361/201424356.
- Zinnecker, Hans; Yorke, Harold W. (2007). "Toward Understanding Massive Star Formation*". Annual Review of Astronomy and Astrophysics. 45 (1): 481–563. Bibcode:2007ARA&A..45..481Z. arXiv: . doi:10.1146/annurev.astro.44.051905.092549.
- Figer, Donald F. (2005). "An upper limit to the masses of stars". Nature. 434 (7030): 192–194. Bibcode:2005Natur.434..192F. PMID 15758993. arXiv: . doi:10.1038/nature03293.
- Kuiper, Rolf; Klahr, Hubert; Beuther, Henrik; Henning, Thomas (2011). "THREE-DIMENSIONAL SIMULATION OF MASSIVE STAR FORMATION IN THE DISK ACCRETION SCENARIO". The Astrophysical Journal. 732 (1): 20. Bibcode:2011ApJ...732...20K. ISSN 0004-637X. arXiv: . doi:10.1088/0004-637X/732/1/20.
- Oh, Sambaran; Kroupa, Pavel; Oh, Seungkyung (2012). "The emergence of super-canonical stars in R136-type star-burst clusters". Monthly Notices of the Royal Astronomical Society. 1208 (2): 826. Bibcode:2012MNRAS.426.1416B. arXiv: . doi:10.1111/j.1365-2966.2012.21672.x.
- Vink, Jorick S.; Heger, Alexander; Krumholz, Mark R.; Puls, Joachim; Banerjee, S.; Castro, N.; Chen, K. -J.; Chene, A. -N.; Crowther, P. A.; Daminelli, A.; Grafener, G.; Groh, J. H.; Hamann, W. -R.; Heap, S.; Herrero, A.; Kaper, L.; Najarro, F.; Oskinova, L. M.; Roman-Lopes, A.; Rosen, A.; Sander, A.; Shirazi, M.; Sugawara, Y.; Tramper, F.; Vanbeveren, D.; Voss, R.; Wofford, A.; Zhang, Y. (2013). "Very Massive Stars (VMS) in the Local Universe". Proceedings of the International Astronomical Union. 10: 51–79. Bibcode:2015HiA....16...51V. arXiv: [astro-ph.SR]. doi:10.1017/S1743921314004657.
- Langer, N. (2012). "Presupernova Evolution of Massive Single and Binary Stars". Annual Review of Astronomy and Astrophysics. 50 (1): 107–164. Bibcode:2012ARA&A..50..107L. arXiv: . doi:10.1146/annurev-astro-081811-125534.
- Meynet, Georges; Georgy, Cyril; Hirschi, Raphael; Maeder, André; Massey, Phil; Przybilla, Norbert; Nieva, M.-Fernanda (January 2011). "Red Supergiants, Luminous Blue Variables and Wolf-Rayet stars: the single massive star perspective". Société Royale des Sciences de Liège, Bulletin. 80: 266–278. Bibcode:2011BSRSL..80..266M. arXiv: .
- Yusof, Norhasliza; Hirschi, Raphael; Meynet, Georges; Crowther, Paul A.; Ekström, Sylvia; Frischknecht, Urs; Georgy, Cyril; Abu Kassim, Hasan; Schnurr, Olivier (2013). "Evolution and fate of very massive stars". Monthly Notices of the Royal Astronomical Society. 433 (2): 1114. Bibcode:2013MNRAS.433.1114Y. arXiv: . doi:10.1093/mnras/stt794.
- Groh, Jose H.; Meynet, Georges; Georgy, Cyril; Ekström, Sylvia (2013). "Fundamental properties of core-collapse supernova and GRB progenitors: Predicting the look of massive stars before death". Astronomy & Astrophysics. 558: A131. Bibcode:2013A&A...558A.131G. arXiv: . doi:10.1051/0004-6361/201321906.
- O'Connor, Evan; Ott, Christian D. (2011). "BLACK HOLE FORMATION IN FAILING CORE-COLLAPSE SUPERNOVAE". The Astrophysical Journal. 730 (2): 70. Bibcode:2011ApJ...730...70O. ISSN 0004-637X. arXiv: . doi:10.1088/0004-637X/730/2/70.
- Valenti, S.; Pastorello, A.; Cappellaro, E.; Benetti, S.; Mazzali, P. A.; Manteca, J.; Taubenberger, S.; Elias-Rosa, N.; Ferrando, R.; Harutyunyan, A.; Hentunen, V. P.; Nissinen, M.; Pian, E.; Turatto, M.; Zampieri, L.; Smartt, S. J. (2009). "A low-energy core-collapse supernova without a hydrogen envelope". Nature. 459 (7247): 674–677. Bibcode:2009Natur.459..674V. PMID 19494909. arXiv: . doi:10.1038/nature08023.