1. Preamble
Massive stars are stars having enough matter to form a non-degenerate iron-nickel core at the end of their lives. Hence, massive stars are born with initial mass (M
) of at least ∼8 M
[
1], and have spectral type O or early-B on the H-burning main sequence, although the spectral type and luminosity class will vary as the star evolves off the main sequence stage. At the end of their lives,
most massive stars implode internally (and simultaneously explode externally) to leave behind neutron stars for M
8–20 M
and black holes for M
20–130 M
and possibly above about 250 M
[
2,
3,
4]. For M
130–250 M
, the pre-supernova star is so luminous that gamma-ray photons in their cores have enough energy to spontaneously create electron-positron pairs that lead to a highly unstable configuration that produces a conflagration supernova, i.e., total dispersal of the entire star via a thermonuclear runaway, resulting in no remnant at all [
5].
The initial mass function describes the relative numbers of stars of different mass that are born in a gravitationally collapsing cloud of gas and dust, and it is approximated by a power law (similar to the turbulence spectrum of the preceding cloudlets that formed the stars) strongly favouring lower masses. Therefore, massive stars are very rare: for every 10 M star there are hundreds of Sun-like stars. However, much more than making up for this low number, an individual massive star outshines its low-mass solar-type cousin by up to a factor of a million. The arms of spiral galaxies are therefore dominated by the light from individual, and clusters of, massive and hence very luminous stars.
However, not for long: after only a few million years, massive stars, that have “lived hard” also “die young”, literally “burning” themselves out quickly as astronomical timescales go. As star formation eventually peters out in any region, this leaves only the lower light level of many more lower-mass stars (with much longer nuclear lifetimes) to shine on for billions of years, each faintly glowing together, as in a globular star-cluster or an elliptical galaxy, where significant star formation has all but ceased long ago. Consequently, while the mass of ordinary matter in the Universe is mostly locked up in low-mass stars, the brilliance of regions of sustained or sporadic star formation, as in most spiral and irregular galaxies, is dominated by the light of massive stellar “queens that radiate in the night”.
The high luminosity of massive stars, especially those of hot O-type stars on the main sequence and their descendants with masses above ∼20 M, adds to yet another of their extreme characteristics, namely the presence of very strong winds of matter accelerated away from the underlying star by radiation pressure. For O-type stars this occurs mostly on ultraviolet (UV) spectral lines especially of iron with the particles making up the wind leaving the stellar photosphere at v km s and being accelerated up to a terminal speed v 2000–3000 km s. A main sequence star of 100 M has a net mass-loss rate that is a billion times greater than that of the Sun, only increasing as it evolves to larger radii, where surface gravity holding back the escaping material, weakens. Thus, massive stars not only have the greatest luminosities among stars but are also the main source of ecological change in the Universe, providing the lion’s share of heavy elements, be it from normal evolution with strong winds, or from core collapse supernovae, or from neutron star/black hole mergers and binary neutron star mergers, all of which come from massive stars.
These massive queens are also among the hottest stars, not only in their cores, but also at their surfaces, since like all stars, they spend most of their lives (about 90%) on the main sequence where they “burn” (i.e., fuse nuclei giving off lots of energy) hydrogen, the most abundant element at 90% by number, into helium. Yet, despite this evolutionary enriching of the interstellar medium, most of the helium and hydrogen still in the Universe today was already formed previously by a similar but much more rapid and intense process of nucleosynthesis in the Big Bang. On the other hand, the majority of the heavier elements are mainly produced in the central nuclear furnaces of massive stars, during their explosions/disintegrations at the end of their lives, or in their afterlives through kilonovae events. After the main sequence H-burning phase, the most massive stars pass through a ten-times shorter He-burning phase as classical Wolf-Rayet stars, with at least ten times stronger winds and, thus, enormous hot spectral emission lines that betray their kind.
2. Enter Puppis
There are numerous regions in the sky where massive stars dominate and are seen visually. In the Milky Way, Orion is the closest example, with much richer regions near the Galactic center and in the Carina Nebula taking a lead role, although less obviously to the naked eye due to their greater distances. The largest collection of massive stars known in the Local Group of galaxies is depicted in
Figure 1.
Among individual bright stars in the sky, the hottest and brightest, both apparently and intrinsically, is the 2nd magnitude, O4If(n)p star
Puppis, also known as Naos (Greek for “Ship”). At a declination of
degrees, it is only visible from southern sites around the globe, best seen in January and adjacent months. Only 7 degrees further south of
Pup lies a multiple system that contains the sky’s brightest Wolf-Rayet star, 2nd magnitude
Velorum, aka Regor.
Pup and
Vel form a bright blue-white pair in the southern Milky Way (
Figure 2).
Pup is a very hot single supergiant star with a mass of 56 M
, recently revised dramatically downwards to potentially be only 25 M
due to a new closer distance derived with revised
HIPPARCOS (HIgh Precision PARallax COllecting Satellite) astrometry [
6].
Pup is so bright that the
Gaia mission is not able to provide reliable astrometry for the star as of the 3rd
Gaia data release, making the
HIPPARCOS measurements the most accurate astrometric measurements available thus far for
Pup. The resulting revised mass would make
Pup an anomaly compared to other similar supergiants, for which “standard” calibrations [
7] give a mass of 58 M
for an O4I star, although such “standard” calibrations indicate that any given single star could be uncertain by to up to 35%–50% in its spectrum-based mass. With an effective surface temperature of ∼40,000 K (as for O4I stars in general [
7]),
Pup is both a runaway (space velocity of ∼60 km s
) and a rapid rotator (
km s
). The latter characteristic is imprinted in its spectral type and makes it belong to the Onfp class of O-type stars [
8] that exhibit a broad, centrally reversed HeII
emission line profile resulting from fast rotation and the so-called “resonance zone effect” [
9,
10].
The fast rotation and runaway status of
Pup are intimately connected within the context of its plausible formation scenarios. One of such scenarios suggests that
Pup might have been the original secondary star in a binary in which the original more massive primary evolved faster and exploded as a supernova, leading via reduced binding gravity to separation of the two stars at relatively high speeds in opposite directions [
11]. The original primary, now an isolated neutron star or more likely a black hole, has never been identified.
Pup would then be the original secondary, spun up by Roche lobe overflow from the rapidly evolving primary just before it went supernova. Alternatively,
Pup could have been spun-up as the product of the merger of at least two massive stars, then gravitationally ejected during dynamical interactions within the Vela R2 R-association [
12,
13].
Being very bright and isolated in the southern sky, Pup has long served as a reference calibrator for O-type stars. However, now as we find more peculiarities, as happens for essentially all stars whenever we look at them more closely, one wonders about its true value as a standard star.
3. Key to Understanding Pup: Variability
The first intensive monitoring effort conducted on
Pup occurred in the mid-seventies and was geared towards investigating the variability of its H
wind emission line, leading to the detection of a periodicity of
d [
14]. This single periodicity was claimed to be associated with stellar rotation along with a speculated global magnetic field that could force trapped hot gas to rotate at the same rate above the star and be observable in H-alpha (see also
Section 4).
Then from 1978 to 1996,
the International Ultraviolet Explorer (IUE) satellite produced a boon in hot-star spectroscopy, including for
Pup. This culminated with a dedicated two-week “Mega”
IUE program in 1995 using all the joint time of the collaborating space agencies, to intensely observe three key luminous hot stars in alternate succession: J Pup (B0.5Ib),
Pup (O4If(n)p), and EZ CMa (WN4b). The end goal was to properly probe the quasi-periodic behaviour of their winds. From that joint effort, the presence of discrete absorption components (DACs) on the blue side of P Cygni spectral lines in the UV for
Pup was confirmed and their characteristics were much better quantified [
15]. Such DACs are now known to be universal in O-star winds, and are best interpreted as the spectroscopic signatures of the presence of large-scale spiral-like structures called Corotating Interaction Regions (CIRs). Some theories then suggested that such large scale wind structures might be triggered by stellar surface perturbations such as bright corotating surface spots [
16]—see
Figure 3. In
Pup, the
IUE Mega data exhibited two periodicities in its DACs: (1) inter-DAC spacing of
h and (2) a second longer period of
d, suspected then to confirm the previously found rotation period of
d. Then, an independent examination of the three years of optical photometry obtained for
Pup by the astrometric satellite
HIPPARCOS yielded a period of
d [
17], which may represent the first harmonic of the 5.1–5.2 d period. Meanwhile, optical spectroscopic monitoring efforts found periods of
h and
h [
18]. The former was claimed to be non-radial pulsations, while the latter coincides with the DAC recurrence time previously found in UV resonance lines.
Another significant discovery was achieved in the late-nineties: the presence of random clumps in an O-star wind—that of
Pup—was revealed for the first time during the last Reticon detector observing-run at the coudé focus of the Canada-France-Hawaii Telescope (CFHT)
-m telescope in 1995 [
19]. The clumps were seen as emission-line subpeaks superposed on the HeII
optical wind line (
Figure 4), interpreted as small-scale wind inhomogeneities propagating radially outward in all parts of the wind. The trajectories of the clumps were found to be entirely compatible with the wind-expansion
-law,
, that had been proposed for all hot-star winds, with
exponent close to unity. This discovery was essentially a follow-up to a similar study made on the >
stronger winds of Wolf-Rayet stars [
20], where the clumps were much more readily seen in repeated spectra with a signal-to-noise ratio of ∼200, although the star itself is hidden by the thick wind. In the case of
Pup, the weakness of even the strongest optical emission lines requires a signal-to-noise ratio of about 1000. Characterizing clumping in hot star winds is crucial as it has been found that unclumped hot star wind models lead to overestimates of the mass-loss rates by a factor of two to ten [
21], which in turn translate into significant uncertainties in the evolutionary tracks of massive stars. Furthermore, for O-type stars, the morphology of clumps affects the morphology of their emission line profiles in the X-ray domain, with radially compressed, pancake-like clumps explaining the observed symmetric X-ray emission line profiles better than spherical clumps [
22].
4. Finally, and Most Recently
The pace picks up in the new millennium. First, the Solar Mass Ejection Imager (SMEI) on board the
Coriolis satellite delivered a mean optical light curve in one filter for
Pup during 2003–2006, that led to a new periodicity of
d [
23]. The monoperiodic signal was tentatively attributed to pulsation.
Then, following the successful launch of the joint Canadian/Austrian/Polish constellation of five operational nanosatellites—the
BRIght Target Explorer (BRITE) mission—we organized a coordinated optical monitoring campaign on
Pup, spanning about
months in 2014–2015. The idea was to use
BRITE to probe light variations coming from the surface of
Pup, and simultaneously monitoring the behaviour of the inner part of its wind by intensively monitoring its HeII
wind emission line from multiple southern ground-based sites involving both amateurs observatories and professional facilities [
24]. The
-d periodicity was detected in the two distinct filters (blue and red) from different
BRITE satellites, although the shape of the phased light curve had changed from one slightly asymmetric peak with SMEI to two, more distinct peaks separated by about
in phase (
Figure 5). The
-d periodicity was also detected in the variations of the HeII
wind emission line. The behaviour revealed by the 2014–2015 observing campaign indicated that this
-d periodicity is unlikely the signature of pulsations, but rather the result of rotational modulation arising from localized corotating features at the surface of the star. More importantly, the correlation observed between the light curve variations and the variations in the wind indicated for the first time that bright surface spots induce CIRs in the wind of
Pup (see
Figure 6 for a cartoon of
Pup’s major CIRs along with a mottled stellar surface).
Although the bright spots are semi-permanent, they might have changed in periodicity on a decade time-scale (
Figure 7, compiling data from various studies [
15,
23,
24,
25,
26]), which is a new twist as to their origin. One speculation for this is that these spots might wander slowly in mean latitude due to differential rotation, going from inter-CIR period of 15 h in 1989 to 19 h in 1995 and remaining constant at
h (i.e., half the
d period) since about 2000. For comparison the Sun’s rotation varies from a period of
d at the equator to almost 38 d at the poles, i.e., with similar ratio as in
Pup, allowing for the fact that we may not have yet seen the full range in
Pup. It is also not yet clear in which direction the period changes in
Pup, i.e., increasing or decreasing from equator to pole. Either could be possible. A more definitive explanation of
Pup’s differential rotation—if this is the correct interpretation—will have to await further constraining data and modeling. Spot modeling via numerical inversion of the 2014–2015
BRITE light curve [
24] indicated that, while the
d period is relatively stable over half a year, there are at least two detectable bright spots at any given time, which change slowly by small amounts in relative intensity and position on timescales of weeks. If
d does in fact represent the equatorial rotation period, then combined with the spectral line broadening
km s
an equatorial radius
R
, one finds the equatorial rotation speed to be close to 500 km s
, i.e., about 80% of break-up, with an inclination of
degrees. Some of these numbers become less extreme for the closer distance from the
HIPPARCOS analysis. Nevertheless, a high rotation speed is in line with its runaway status and history as noted above.
In 2018, a multi-wavelength variability study that compiled X-ray data from the
XMM-Newton and
Swift satellites along with optical monitoring with the
Coriolis and
BRITE satellites found no obvious correlation between the observations taken in these two wavelength regimes [
27]. However, the X-ray observations were sparsely sampled. It then became clearer that there was a need for a longer, more intensive, contemporaneous multi-wavelength monitoring. Hence, in 2018–2019, we co-organized a joint multi-wavelength observing campaign on
Pup with precision optical photometric monitoring during a contiguous 5-month run with
BRITE, in conjunction with the already planned
Chandra X-ray observing campaign (PI Wayne Waldron) for a total of just under ten days spread unevenly over 14 months. This proved to be extremely useful. The
d periodicity was present in both the optical and the X-ray observations [
28]. This may point to the origin of at least some of the X-rays to be the shocks associated with the CIRs and modulated by rotation via absorption and eclipse effects. However, the peak-to-valley amplitude of the
d modulation in X-rays was close to 6%, significantly greater than the 1% amplitude in the optical. Furthermore, although the
d phased X-ray light curve showed similar-shaped “bumps” compared to the nearly simultaneous
BRITE light curve, with the secondary bump in the X-ray observations being much less pronounced, there was a phase shift between the two domains (optical vs. X-ray). The phase shift amounted to
in phase (clearly less than half a rotation) if the two maxima are matched up, or
in phase if the stronger optical peak is made to coincide with the weaker X-ray peak and vice versa (
Figure 8). Theoretically, there is a problem with the larger phase shift, since it would require CIR X-rays to be formed implausibly far out in the wind, given the predictable shape of the CIRs [
16]. The shorter phase shift is supported by similar phase shifts observed earlier in optical emission lines (formed in the wind adjacent to the continuum arising below the wind in the photosphere), that also follow the
d modulation [
24]. However, then the question remains why the more intense X-ray light curve peak does not coincide with the more intense optical light curve peak. A conceivable reason may be that not all CIRs are created equal.
Pup is not the only hot massive star now known to exhibit correlated multi-wavelength periodic variability. Other cases include
Per [
29],
Cep [
30] and
Oph [
31], all runaway O-type stars. However,
Pup is probably the clearest case so far.
Interestingly, while there is no trace of any 5 d period in the BRITE or SMEI data, the Chandra data show a 5 d and 6 d period (although probably only one or the other is real, not both simultaneously) of amplitude only slightly less than that of the d modulation. Although there is currently no clear explanation for this, one of the most exotic speculations is that it is some kind of more or less stable body (dense gas, planet?) in a ∼5-d orbit around Pup, in agreement with the 5 d periodicity seen before in some observations. Unlike the slowly varying CIR period though, such an orbital period, if real, would be much more stable in time. There is a basic problem with this scenario, though: given Pup’s likely history as a runaway after a supernova explosion in a binary, the survivability of a close 5 d orbiting body during the explosion would be questionable.
Returning to the
d period found by
BRITE, it was a great help that different
BRITE satellites each with either a blue or a red optical filter, revealed essentially identical light curves (
Figure 6), with Fourier peaks at
d and its first harmonic and nothing else of significance. Just as important, though, was that the phased light curves showed a level of completely stochastic variability at about the same 1% level peak-to-valley as the
d modulation, not obvious in the SMEI or any previous precision photometry where it is not easy to separate intrinsic from instrumental noise. However, with
BRITE this was possible and it became clear that these random variations in the photosphere were correlated with clumps in the adjacent wind [
24], based on high-quality optical spectra of the HeII
wind emission line obtained simultaneously with the
BRITE photometry. With the advent of even higher-precision photometric monitoring from space, it is becoming obvious that virtually all massive hot stars show similar stochastic variability [
32,
33].
Curiously, one does see stochastic variability in
Pup in the optical; why is this not also clearly seen in X-rays? Could it just be because the X-ray photon count-rate is not high enough to attain a significant level of certainty? On the ond hand, one expects to see X-rays globally from all parts of the wind caused by many individual independent shocks. On the other hand, it has been proposed that the number of shocks (and associated clumps) at any given time might be so great (>
) that their relative variability (<
in amplitude) is drowned in Poisson noise [
34]. In X-rays, the only apparent way forward may be via future X-ray telescopes in space with larger collecting area leading to higher sensitivity to be able to detect any significant stochastic X-ray variability.