VY Canis Majoris

VY Canis Majoris (abbreviated to VY CMa) is an extreme oxygen-rich (O-rich) red hypergiant (RHG) or red supergiant (RSG) and pulsating variable star located at 1.2 kpc away from Earth in the constellation of Canis Major. It is one of the largest known stars by radius, and is also one of the most luminous and massive red supergiants, as well as one of the most luminous stars in the Milky Way.

VY CMa is a single star with a large infrared (IR) excess, making it one of the brightest objects in the sky at wavelengths of between 5 and 20 microns (µm) and indicating a dust shell or disk heated by the star. It is about $3,820$ times the mass of the Sun. It is also surrounded by a complex asymmetric circumstellar envelope (CSE) caused by mass loss from the star itself. It produces strong molecular maser emission and was one of the first radio masers discovered. VY CMa is embedded within the large molecular cloud Sharpless 310 (Sh2-310), one of largest star-forming H region with a diameter of 480 arcminutes (') or 681 ly.

The diameter of VY CMa is around 1,420–2,000 times that of the Sun, which is close to the Hayashi limit and corresponds to a volume nearly 3–10 billion times bigger than the Sun. A hypothetical object travelling at the speed of light would take 6–9 hours to travel around the star's circumference, compared to 14.5 seconds for the Sun. If placed at the center of the Solar System, VY CMa's surface would extend beyond the orbit of Jupiter, although there is still considerable variation in estimates of the radius, with some making it larger than the orbit of Saturn.

Observational history
The first known recorded observation of VY Canis Majoris is in the star catalogue of the French astronomer Jérôme Lalande, on 7 March 1801, which lists it as a 7th magnitude star. Further 19th- and 20th-century studies of its apparent magnitude suggested that the star had been fading since 1850.

Since 1847, VY Canis Majoris has been described as a crimson star. During the 19th century, observers measured at least six discrete components, suggesting that it might be a multiple star. These discrete components are now known to be bright areas in the surrounding nebula. Visual observations in 1957 and high-resolution imaging in 1998 showed that there are no companion stars. VY CMa was also discovered to be a strong source of OH (1612 MHz), (22235.08 MHz), and  (43122 MHz) masers emission, which is typical of an OH/IR star. Many molecules, such as, , , , , , , and , have also been detected.

The variation in VY CMa's brightness was first described in 1931 when it was listed (in German) as a long period variable with a photographic magnitude range of 9.5 to 11.5. It was given the variable star designation VY Canis Majoris in 1939, the 43th variable star of the constellation Canis Major.

Surroundings
VY Canis Majoris is surrounded by an extensive and dense asymmetric red reflection nebula with a total ejected mass of and a temperature of 800 K based on a DUSTY model atmosphere that has been formed by material expelled from its central star. The diameter of the inner shell is found to be 0".12, corresponding to 140 AU at a distance of 1.2 kpc, whereas that of the outer one is 10", corresponding to 12000 AU. This nebula is so bright that it was discovered in 1917 with an 18 cm telescope, and also contains condensations that were once regarded as companion stars. It has been extensively studied with the aid of the Hubble Space Telescope (HST), showing that the nebula has a complex structure that includes filaments and arcs, which were caused by past eruptions; this structure is similar to that of the nebula surrounding the post-red supergiant (Post-RSG) or yellow hypergiant (YHG) IRC +10420. The similarity has led astronomers to propose that VY CMa will evolve blueward on the Hertzsprung–Russell diagram (HR diagram) to become a yellow hypergiant, then a Luminous blue variable (LBV), and finally a Wolf–Rayet star (WR star).

Observations combining Hubble Space Telescope data with that obtained from the Keck Telescope, located in Hawaii, it was possible to make a three-dimensional reconstruction of the material surrounding VY CMa. This reconstruction showed that VY CMa's mass loss is much more complex than expected for any red supergiant or hypergiant. It became clear that the bows and nodules appeared at different times; the jets are randomly oriented, which leads us to believe that they derive from explosions occurred in active regions of the stellar photosphere. The spectroscopic data showed that the jets move away from the star at different speeds, indicating that they occurred at different times and originate from different areas of the surface of VY CMa. Multiple asymmetric mass loss events and also the ejection of the outermost material are deduced to have occurred within the last 500–1,000 years, while that of a knot near the star would be less than 100 years. The mass loss is due to strong convection in the tenuous outer layers of the star, associated with magnetic fields. This is similar to the sunspots and coronal ejections of the Sun but on a much larger scale.

Distance
In 1976, Charles J. Lada and Mark J. Reid published observations of the bright-rimmed molecular cloud Sh2-310, which is 15' east of VY Canis Majoris. At the edge of the cloud bordered by the bright rim, an abrupt decrease in the CO emission and an increase in brightness of the emission were observed, indicating possible destruction of molecular material and enhanced heating at the cloud-rim interface, respectively. Lada and Reid assumed the distance of Sh2-310 is approximately equal to that of the stars, which are members of the open cluster NGC 2362, that ionize the rim. NGC 2362 has a distance of $1,170$ (kpc) or roughly $17$ (ly) as determined from its color-magnitude diagram. VY CMa is projected onto the tip of the cloud rim, suggesting its association with Sh2-310. In addition to that, the velocity of Sh2-310 is very close to the velocity of the star. This further indicates the association of the star with Sh2-310, and consequently with NGC 2362, which means VY Canis Majoris is at the same distance. A more recent measurement of the distance to NGC 2362 gives 1.2 kpc or about 3,910 ly.

Stellar distances can be calculated by measuring parallaxes as Earth orbits around the Sun. However, VY Canis Majoris has a small parallax due to its great distance, and standard visual observations have a margin of error too large for a hypergiant star with an extended CSE to be useful, for example, the Hipparcos Catalogue (1997) gives a parallax of $1,420$ (mas), yelding a distance of $2,000$ or about $270,000$. The parallax of VY CMa can be precisely measured with high accuracy from the observation of masers using a long baseline interferometry. In 2008, observations of masers using VERA interferometry from the National Astronomical Observatory of Japan give a parallax of $3,490$, corresponding to a distance of $0.6$ or about $17$. In 2012, observations of masers using Very-long-baseline interferometry (VLBI) from Very Long Baseline Array (VLBA) independently derived a parallax of $1.5 kiloparsecs$, corresponding to a distance of $4,890 light-years$ or about $1.78 milliarcseconds$. These new distance estimates imply that either Sh2-310 is less remote than usually estimated or that VY CMa is a foreground object.

The Gaia mission should provide visual parallaxes with sufficient accuracy to constrain the distance to VY CMa, but the data release 2 value of $561.8 pc$ is not meaningful.

Variability
VY Canis Majoris is a variable star that varies in brightness from apparent magnitude of +6.5 to +9.6. in the General Catalogue of Variable Stars (GCVS) it is classed a semiregular variable of sub-type SRc, indicating a cool supergiant, with an estimated pulsational period of 956 days, although it is classed as a type LC slow irregular variable star in the American Association of Variable Star Observers (AAVSO) Variable Star Index. Other periods of 1,600 and 2,200 days have been derived.

VY CMa is sometimes considered as the prototype for a class of heavily mass-losing OH/IR supergiants, distinct from the more common asymptotic giant branch OH/IR stars.

Spectrum
The spectrum of VY Canis Majoris is that of a high luminosity M class star. The hydrogen lines, however, have P Cygni profiles fit for Luminous Blue Variables. The spectrum is dominated by TiO bands whose strengths suggest a classification of M5. The H-alpha (Hα) line is not visible yet and there are unusual emission lines of neutral elements such as sodium and calcium. The luminosity class as determined from different spectral features varies from bright giant (II) to bright supergiant (Ia), with a compromise being given as M5eIbp. Early attempts at classification were confused by the interpretation of surrounding nebulosity as companion stars.

The derived spectral class varies depending on the features examined. The spectral features also vary noticeably over time. It is considered to be unambiguously cooler, and thus redder, than M2, and is usually classified between M3 and M5. Classes as extreme as M2.5 and M5 have been given. The luminosity class is likewise confused and often given only as I, partly because luminosity classes are poorly defined in the red and infrared portions of the spectrum. One study though, gives a luminosity class of Ia+ which means a hypergiant or extremely luminous supergiant.

Properties
A very large and luminous star, VY CMa is among the most extreme stars in the Milky Way and has an effective temperature below 4,000 K. It occupies the upper-right hand corner of the HR diagram although its exact luminosity and temperature are uncertain. Most of the properties of the star depend directly on its distance.

Luminosity
The bolometric luminosity (Lbol) of VY CMa can be calculated from Spectral energy distribution (SED) or bolometric flux, which can be determined from photometry in several visible and infrared bands. Earlier calculations of the luminosity using an assumed distance of 1.5 kpc gave luminosities between 200,000 and 500,000 solar luminosities, very close to the empirical upper limit of luminosity for the cool hypergiants (e.g. Humphreys-Davidson limit). One study gave nearly at a distance of 2.1 kpc. In 2006 a luminosity of was calculated by integrating the total fluxes over the entire nebula, since most of the radiation coming from the star is reprocessed by the dust in the surrounding cloud. More recent estimates of the luminosity derive values below based on distances below 1.2 kpc. Most of the output of VY CMa is emitted as infrared radiation, with a maximum emission at 5–10 μm, which is in part caused by reprocessing of the radiation by the circumstellar nebula. Many older luminosity estimates are consistent with current ones if they are re-scaled to the distance of 1.2 kpc. Despite being one of the most luminous stars in the Milky Way, much of the visible light of VY CMa is absorbed by the circumstellar envelope, so it cannot be seen with the naked eye and needs a telescope to be observed. It would be a naked eye star if no light was absorbed.

Mass
Since VY CMa has no companion star, its mass cannot be measured directly through gravitational interactions. Comparison of the effective temperature and bolometric luminosity of VY CMa compared to evolutionary tracks for massive stars suggest an initial mass of $1,832.34 ly$ for a rotating star with a current mass of, or for a non-rotating star with a current mass of  , and an age of 8.2 million years (MYr). Older studies have found much higher initial masses (thus also higher current masses) or a progenitor mass of based on old luminosity estimates.

Mass loss
VY CMa has a strong stellar wind and is losing a lot of material due to its high luminosity and relatively low surface gravity. It has an average mass loss rate of $0.88 mas$ per year, among the highest known and unusually high even for a red supergiant, as evidenced by the extensive circumstellar envelope. It is consequently one of most important stars for understanding the high-mass loss episodes near the end of massive star evolution. The mass loss rate probably exceeded $1.14 kpc$/yr during the most violent mass loss events.

Temperature
The effective temperature of VY CMa is uncertain but as the spectrum of VY CMa varies, then also the temperature may vary. Early estimates of the temperature of VY CMa assumed values below 3,000 K based on a spectral class of M5. In 2006, Dr. Philip Massey used a MARCS model atmosphere to directly determine the temperature of VY CMa and calculated that the effective temperature of VY CMa is as high as $3,720 ly$, corresponding to a spectral class of M2.5, although VY CMa is usually considered as an M4-M5 star. Adopting the spectral class of M4-M5 with the temperature scale proposed by Emily Levesque gives a range of between 3,450 and 3,535 K.

Size
The calculation of the radius of VY CMa is complicated by the extensive circumstellar envelope of the star and also VY CMa is a pulsating star, so its size changes with time. Earlier direct measurements of the radius at infrared (K-band = 2.2 µm) wavelength gave an angular diameter of $0.83 mas$, corresponding to radii above at an assumed distance of 1.5 kpc, considerably larger than expected for any red supergiant or hypergiant. However, this is probably larger than the actual size of the underlying star and the angular diameter estimate appears exceedingly large due to interference by the circumstellar envelope. In 2006–2007 radii of 1,800 - 2,100 solar radius have been dervied by Humphreys from the estimated luminosity of and temperatures of 3,200 - 3,535 K.

On 6 and 7 March 2011, VY CMa was observed by Wittkowski et al. using VLTI/AMBER interferometry with the Very Large Telescope (VLT) in a range of 2.20–2.25 μm. The size of VY CMa was calculated using the Rosseland Radius, the location at which the optical depth is $1.2 kpc$, with two modern distances of $3,910 ly$ and $-5.92 mas$. The angular diameter of VY CMa was directly and precisely measured at $25$ based on a PHOENIX atmosphere model, which corresponds to a radius of $6$ at a distance of $$. The high spectral resolution of these observations allowed to minimize the effect of contamination by the circumstellar layers. An effective temperature of $3,650 K$ was then derived from the radius and a luminosity of $18.7 mas$ which is based on distance and measured flux of $2/3$, resulting to a spectral class of M4. In 2013, Alcolea et al. referred to VY Canis Majoris as the largest star discovered so far, referring to the Wittkowski radius, which is considered the highest radius among well-characterised stars. Several other red supergiants or hypergiants are possibly larger, but they have less accurate radius estimates.

Most radius estimates of the VY CMa are considered as the size for the optical photosphere while the size of the star for the radio photosphere is calculated to be twice that of the size of the star for the optical photosphere.

Despite the mass and very large size, VY CMa has an average density of 5.33 to 8.38 mg/m3 (0.00000533 to 0.00000838 kg/m3), it is over 100,000 times less dense than Earth's atmosphere at sea level (1.2 kg/m3).

Largest star
VY Canis Majoris has been known to be an extreme object since the middle of the 20th century, although its true nature was uncertain. In the late 20th century, it was accepted that it was a post-main sequence red supergiant. Its angular diameter had been measured and found to be significantly different depending on the observed wavelength. The first meaningful estimates of its properties showed a very large star, well in excess of.

In contrast to prevailing opinion, a 2006 study, based on an assumed effective temperature of 3,650 K and distance of $1.14$, derived a luminosity of, suggesting an initial mass of and radius of. On this basis they considered VY CMa and another notable extreme cool hypergiant star, NML Cygni, as normal early-type red supergiants. They assert that earlier very high luminosities of $1.2 kpc$ and very large radii of  (or even  ) were based on effective temperatures below 3,000 K that were unreasonably low.

Almost immediately another paper published a size estimate of and concluded that VY CMa is a true hypergiant. This was based on the most recent temperature calibrations, giving an effective temperature of 3,450 - 3,535 K, and a luminosity of based on SED integration and a distance of $11.3 mas$.

Since then, the size of VY CMa has been calculated more accurately to be somewhat lower, for example, while larger sizes have been published for a number of other galactic red supergiants such as Westerlund 1-26. Despite this, VY Canis Majoris is still often described as the largest known star, sometimes with caveats to account for the highly uncertain sizes of all these stars. Another recent estimate published a size of, similar to the Humphreys radius.

Evolution
VY Canis Majoris is a highly evolved star with an age less than 10 MYr although some older authors argued that it would rather either be a very young protostar or a massive pre-main sequence star with an age of only 1 MYr and also a circumstellar disk around the star. It has probably evolved from an O9 main sequence star with a radius of. The star has evolved rapidly because of its high mass. The time spent to the red hypergiant phase is estimated to be between 100,000 and 500,000 years, and thus VY CMa had left its main sequence more than a million years ago.

The future evolution of VY CMa is uncertain, but like the most cool supergiants, the star will certainly explode as a supernova. It has begun to fuse helium into carbon, whereas a main sequence star fuses hydrogen into helium. In 2009, Nathan Smith of the University of California, Ken Hinkle and Nils Ride of the Lund Observatory have used the PHOENIX near-infrared spectrograph on Gemini South to study the geometry and kinematics of the active circumstellar envelopes of two large supergiant stars: Betelgeuse and VY Canis Majoris. These two stars are losing huge amounts of mass in their "last gasps" and could explode as a supernova at any time now although the stars like VY CMa (>) are too massive to produce type II-P supernovae at red supergiant stage and will probably evolve to higher temperature, likely producing a different type of supernova explosion. It is also expected that VY CMa will explode in the next 100,000 years. VY CMa has been discovered to be very unstable, having a prodigious mass loss rate. In this case, the CO emission is coincident with the bright KI shell in its asymmetric nebula. VY CMa will produce a moderately luminous and long-lasting type IIn supernova (SN IIn), or even a hypernova or superluminous supernova (SLSN) comparable to SN 1988Z (or possibly a type Ib supernova), but it is unlikely that the supernova would be extremely luminous like SN 2006tf or SN 2006gy. The explosion can be associated with gamma-ray bursts (GRB) and it will also produce a shock wave of a speed of a few thousand kilometers-per-second that could slam into the surrounding envelope of material and detectable interaction with the surrounding envelope for about a decade after the explosion. This demonstrates that LBVs, such as Eta Carinae, are not the only progenitors of SNe IIn or HN, but it stressed the requirement that these types of supernovae increased episodic mass loss just prior to the explosion. For a star as massive as VY CMa, the remnant would be likely a black hole rather than a neutron star.