Massive Stellar Progenitor Research Articles

The Jets and the Neutron Star Kick Velocity of the Supernova Remnant CTB 1

We identify jet-shaped morphology in the core-collapse supernova remnant (SNR) CTB 1 that includes two opposite structural features. We identify these as the imprints of a pair of jets that were among the last jets to explode the massive stellar progenitor of CTB 1. We find the projected angle between the jets’ axis and the direction of the pulsar velocity, which is the neutron star natal kick, to be α = 78°. We tentatively identify possible signatures of a second pair of opposite jets along a different direction. If this identification holds, SNR CTB 1 has a point-symmetric structure. The morphology and large angle of the jets’ axis to kick velocity are expected in the jittering jets explosion mechanism of core-collapse supernovae.

Density Profiles of Collapsed Rotating Massive Stars Favor Long Gamma-Ray Bursts

Long-duration gamma-ray bursts (lGRBs) originate in relativistic collimated outflows—jets—that drill their way out of collapsing massive stars. Accurately modeling this process requires realistic stellar profiles for the jets to propagate through and break out of. Most previous studies have used simple power laws or pre-collapse models for massive stars. However, the relevant stellar profile for lGRB models is in fact that of a star after its core has collapsed to form a compact object. To self-consistently compute such a stellar profile, we use the open-source code GR1D to simulate the core-collapse process for a suite of low-metallicity rotating massive stellar progenitors that have undergone chemically homogeneous evolution. Our models span a range of zero-age main-sequence (ZAMS) masses: M ZAMS = 13, 18, 21, 25, 35, 40, and 45M ☉. All of these models, at the onset of core-collapse, feature steep density profiles, ρ ∝ r −α , with α ≈ 2.5, which would result in jets that are inconsistent with lGRB observables. We follow the collapses of four of the seven models until they form black holes (BHs) and the other three models until they form proto-neutron stars (PNSs). We find, across all models, that the density profile outside the newly formed BH or PNS is well represented by a flatter power law with α ≈ 1.35–1.55. Such flat density profiles are conducive to the successful formation and breakout of BH-powered jets and are, in fact, required to reproduce observable properties of lGRBs. Future models of lGRBs should be initialized with shallower post-collapse stellar profiles, like those presented here, instead of the much steeper pre-collapse profiles that are typically used.

Late-time Evolution and Modeling of the Off-axis Gamma-Ray Burst Candidate FIRST J141918.9+394036

Abstract We present new radio and optical data, including very-long-baseline interferometry, as well as archival data analysis, for the luminous, decades-long radio transient FIRST J141918.9+394036. The radio data reveal a synchrotron self-absorption peak around 0.3 GHz and a radius of around 1.3 mas (0.5 pc) 26 yr post-discovery, indicating a blastwave energy ∼5 × 1050 erg. The optical spectrum shows a broad [O iii]λ4959,5007 emission line that may indicate collisional excitation in the host galaxy, but its association with the transient cannot be ruled out. The properties of the host galaxy are suggestive of a massive stellar progenitor that formed at low metallicity. Based on the radio light curve, blastwave velocity, energetics, nature of the host galaxy and transient rates, we find that the properties of J1419+3940 are most consistent with long gamma-ray burst (LGRB) afterglows. Other classes of (optically discovered) stellar explosions as well as neutron star mergers are disfavored, and invoking any exotic scenario may not be necessary. It is therefore likely that J1419+3940 is an off-axis LGRB afterglow (as suggested by Law et al. and Marcote et al.), and under this premise the inverse beaming fraction is found to be f b − 1 ≃ 280 − 200 + 700 , corresponding to an average jet half-opening angle < θ j > ≃ 5 − 2 + 4 degrees (68% confidence), consistent with previous estimates. From the volumetric rate we predict that surveys with the Very Large Array, Australian Square Kilometre Array Pathfinder, and MeerKAT will find a handful of J1419+3940-like events over the coming years.

Studying the environment of AT 2018cow with MUSE

ABSTRACT AT 2018cow was the nearest and best-studied example of a new breed of extragalactic, luminous, and rapidly evolving transient. Both the progenitor systems and explosion mechanisms of these rapid transients remain a mystery – the energetics, spectral signatures, and time-scales make them challenging to interpret in established classes of supernovae and tidal disruption events. The rich, multiwavelength data set of AT 2018cow has still left several interpretations viable to explain the nature of this event. In this paper, we analyze integral-field spectroscopic data of the host galaxy, CGCG 137-068, to compare environmental constraints with leading progenitor models. We find the explosion site of AT 2018cow to be very typical of core-collapse supernovae (known to form from stars with MZAMS ∼ 8−25 M⊙), and infer a young stellar population age at the explosion site of few × 10 Myr, at slightly sub-solar metallicity. When comparing to expectations for exotic intermediate-mass black hole (IMBH) tidal disruption events, we find no evidence for a potential host system of the IMBH. In particular, there are no abrupt changes in metallicity or kinematics in the vicinity of the explosion site, arguing against the presence of a distinct host system. The proximity of AT 2018cow to strong star formation in the host galaxy makes us favour a massive stellar progenitor for this event.

Prospects for multi-messenger extended emission from core-collapse supernovae in the Local Universe

Multi-messenger emissions from SN1987A and GW170817/GRB170817A suggest a Universe rife with multi-messenger transients associated with black holes and neutron stars. For LIGO-Virgo, soon to be joined by KAGRA, these observations promise unprecedented opportunities to probe the central engines of core-collapse supernovae (CC-SNe) and gamma-ray bursts. Compared to neutron stars, central engines powered by black hole-disk or torus systems may be of particular interest to multi-messenger observations by the relatively large energy reservoir $E_J$ of angular momentum, up to 29\% of total mass in the Kerr metric. These central engines are expected from relatively massive stellar progenitors and compact binary coalescence involving a neutron star. We review prospects of multi-messenger emission by catalytic conversion of $E_J$ by a non-axisymmetric disk or torus. Observational support for this radiation process is found in a recent identification of ${\cal E}\simeq (3.5\pm1)\%M_\odot c^2$ in Extended Emission to GW170817 at a significance of 4.2\,$\sigma$ concurrent with GRB170817A. A prospect on similar emissions from nearby CC-SNe justifies the need for all-sky blind searches of long duration bursts by heterogeneous computing.

SN 2016coi (ASASSN-16fp): An Energetic H-stripped Core-collapse Supernova from a Massive Stellar Progenitor with Large Mass Loss

Abstract We present comprehensive observations and analysis of the energetic H-stripped SN 2016coi (a.k.a. ASASSN-16fp), spanning the γ-ray through optical and radio wavelengths, acquired within the first hours to ∼420 days post explosion. Our observational campaign confirms the identification of He in the supernova (SN) ejecta, which we interpret to be caused by a larger mixing of Ni into the outer ejecta layers. By modeling the broad bolometric light curve, we derive a large ejecta-mass-to-kinetic-energy ratio (M ej ∼ 4–7 M ⊙, E k ∼ (7–8) × 1051 erg). The small [Ca ii] λλ7291,7324 to [O i] λλ6300,6364 ratio (∼0.2) observed in our late-time optical spectra is suggestive of a large progenitor core mass at the time of collapse. We find that SN 2016coi is a luminous source of X-rays (L X > 1039 erg s−1 in the first ∼100 days post explosion) and radio emission (L 8.5 GHz ∼ 7 × 1027 erg s−1 Hz−1 at peak). These values are in line with those of relativistic SNe (2009bb, 2012ap). However, for SN 2016coi, we infer substantial pre-explosion progenitor mass loss with a rate ∼ (1–2) × and a sub-relativistic shock velocity v sh ∼ 0.15c, which is in stark contrast with relativistic SNe and similar to normal SNe. Finally, we find no evidence for a SN-associated shock breakout γ-ray pulse with energy E γ > 2 × 1046 erg. While we cannot exclude the presence of a companion in a binary system, taken together, our findings are consistent with a massive single-star progenitor that experienced large mass loss in the years leading up to core collapse, but was unable to achieve complete stripping of its outer layers before explosion.

An external-shock model for gamma-ray burst afterglow 130427A

The complex multiwavelength emission of GRB afterglow 130427A (monitored in the radio up to 10 days, in the optical and X-ray until 50 days, and at GeV energies until 1 day) can be accounted for by a hybrid reverse-forward shock synchrotron model, with inverse-Compton emerging only above a few GeV. The high ratio of the early optical to late radio flux requires that the ambient medium is a wind and that the forward-shock synchrotron spectrum peaks in the optical at about 10 ks. The latter has two consequences: the wind must be very tenuous and the optical emission before 10 ks must arise from the reverse-shock, as suggested also by the bright optical flash that Raptor has monitored during the prompt emission phase (<100 s). The VLA radio emission is from the reverse-shock, the Swift X-ray emission is mostly from the forward-shock, but the both shocks give comparable contributions to the Fermi GeV emission. The weak wind implies a large blast-wave radius (8 t_{day}^{1/2} pc), which requires a very tenuous circumstellar medium, suggesting that the massive stellar progenitor of GRB 130427A resided in a super-bubble.

Modifying two-body relaxation in N-body systems by gas accretion

We consider the effects that accretion from the interstellar medium onto the particles of an N-body system has on the rate of two-body relaxation. To this end, we derive an accretion-modified relaxation time by adapting Spitzer's two-component model to include the damping effects of accretion. We consider several different mass-dependencies and efficiency factors for the accretion rate, as well as different mass ratios for the two components of the model. The net effect of accretion is to accelerate mass segregation by increasing the average mass $\bar{m}$, since the relaxation time is inversely proportional to $\bar{m}$. Under the assumption that the accretion rate increases with the accretor mass, there are two additional effects that accelerate mass segregation. First, accretion acts to increase the range of any initial mass spectrum, quickly driving the heaviest members to even higher masses. Second, accretion acts to reduce the velocities of the accretors due to conservation of momentum, and it is the heaviest members that are affected the most. Using our two-component model, we quantify these effects as a function of the accretion rate, the total cluster mass, and the component masses. We conclude by discussing the implications of our results for the dynamical evolution of primordial globular clusters, primarily in the context of black holes formed from the most massive stellar progenitors.

LOCAL ABSORPTION OF HIGH-ENERGY EMISSION FROM GAMMA-RAY BURSTS

High-energy photons emitted from gamma-ray bursts (GRBs) are subject to pair-production interactions with lower energy photons, leading to an effective optical depth. In this Letter, we estimate the opacity resulting from photon fields located at various distances from long GRB sites: that of the binary companion to the massive stellar progenitor, that of the star-forming molecular cloud containing the GRB, and the total photon field of the host galaxy. The first two photon fields are found to be transparent for most reasonable sets of assumptions about these systems. In the case of galactic radiation fields, we have performed several numerical simulations to calculate the expected opacities for different line-of-sight geometries through the host galaxy, and include a full accounting of the infrared radiation produced by the absorption and re-radiation of starlight by dust. The optical depth for GeV gamma-rays, due to direct starlight is less than unity for all host galaxies. At higher energies, $>$10 TeV, a spectral cutoff can occur due to the rapidly increasing number of mid- to far-IR intra-galactic photons reradiated by dust. Photons in the extragalactic background light therefore remain the only relevant source of photon-photon opacity for ongoing GRB observations with Fermi LAT, and potential future detections with ground-based gamma-ray telescopes.

A ultraluminous X-ray source associated with a cloud collision in M 99

The Sc galaxy M 99 in the Virgo Cluster has been strongly affected by tidal interactions and recent close encounters, responsible for an asymmetric spiral pattern and a high star formation rate. Our XMM–Newton study shows that the inner disc is dominated by hot plasma at kT≈ 0.30 keV, with a total X-ray luminosity of ≈1041 erg s−1 in the 0.3–12 keV band. At the outskirts of the galaxy, away from the main star-forming regions, there is an ultraluminous X-ray source (ULX) with an X-ray luminosity of ≈2 × 1040 erg s−1 and a hard spectrum well fitted by a power law of photon index Γ≈ 1.7. This source is close to the location where a massive H i cloud appears to be falling on to the M 99 disc at a relative speed of >100 km s−1. We suggest that there may be a direct physical link between fast cloud collisions and the formation of bright ULXs, which may be powered by accreting black holes with masses ∼100 M⊙. External collisions may trigger large-scale dynamical collapses of protoclusters, leading to the formation of very massive (≳200 M⊙) stellar progenitors; we argue that such stars may later collapse into massive black holes if their metal abundance is sufficiently low.

400 new Planetary Nebulae in the Galactic Bulge

About 500 MASH(Macquarie/AAO/Strasbourg H $\alpha$ planetary nebula project) PNe were discovered towards the galactic bulge (Parker et al. 2006). For 405 MASH PNe with a diameter $ 35 arcsec, we obtain surface brightness, diameter and dynamical age of the nebulae. From line intensity ratios of 133 GBPN observed with the 6dF device, we obtain their density and ionized mass, and their central star temperatures. We discover 15 bipolar and/or nitrogen-rich GBPNe having probably massive stellar progenitors (Peyaud, 2005).

Neutrino tomography of gamma ray bursts and massive stellar collapses

Neutrinos at energies above TeV can serve as probes of the stellar progenitor and jet dynamics of gamma ray bursts arising from stellar core collapses. They can also probe collapses which do not lead to gamma rays, which may be much more numerous. We calculate detailed neutrino spectra from shock accelerated protons in jets just below the outer stellar envelope, before their emergence. We present neutrino flux estimates from such preburst jets for two different massive stellar progenitor models. These should be distinguishable by IceCube, and we discuss the implications.

Soft X-ray emission lines in the afterglow spectrum of GRB 011211: A detailed XMM-Newton analysis

We report on an XMM-Newtonobservation of the X-ray afterglow of the Gamma Ray Burst GRB 011211, originally detected by Beppo-SAX on 11th December 2001. The early afterglow spectrum obtained by XMM-Newton, observed 11 hours after the initial burst, appeared to reveal decaying H-like K emission lines of Mg, Si, S, Ar and Ca, arising in enriched material with an outflow velocity of order 0.1c (Reeves et al. 2002). This was attributed to matter ejected from a massive stellar progenitor occurring shortly before the burst itself. Here, we present a detailed re-analysis of the XMM-NewtonEPIC observations of GRB 011211. In particular, we show that the detection of the soft X-ray line emission appears robust, regardless of detector background, calibration, spectral binning, or the spectral model that is assumed. We demonstrate that thermal emission, from an optically thin plasma, is the most plausible model that can account for the soft X-ray emission, which appears to be the case for at least two burst afterglow spectra observed by XMM-Newton. The X-ray spectrum of GRB 011211 appears to evolve with time after the first 10 ks of the XMM-Newtonobservation as the Si and S emission lines are only detected during the first 10 ks of observation. The observations suggest that thermal emission is present during the early afterglow spectrum, whilst a power-law component dominates the latter stages. Finally we estimate the mass of the ejected material in GRB 011211 to be of the order 4-20 solar masses.

A Radio Flare from GRB 020405: Evidence for a Uniform Medium around a Massive Stellar Progenitor

We present observations of GRB 020405 starting 1.2 days after the burst that reveal a rapidly fading radio flare. Based on its temporal and spectral properties, we interpret the flare as emission from the reverse shock. This scenario rules out a circumburst medium with a radial density profile ρ ∝ r-2 expected around a mass-losing massive star since in that case, the reverse-shock emission decays on the timescale of the burst duration t ~ 102 s. Using published optical and X-ray data, along with the data presented here, we further show that a self-consistent model requires collimated ejecta with an opening angle, θj ~ 6° (tj ≈ 0.95 days). As a consequence of the early jet break, the late-time (t > 10 days) emission measured with the Hubble Space Telescope significantly deviates from an extrapolation of the early, ground-based data. This, along with an unusually red spectrum, Fν ∝ ν-3.9, strengthens the case for a supernova that exploded at about the same time as GRB 020405, thus pointing to a massive stellar progenitor for this burst. This is the first clear association of a massive progenitor with a uniform medium, indicating that a ρ ∝ r-2 profile is not a required signature and, in fact, may not be present on the length scales probed by the afterglow in the majority of bursts.

Collapsar Uncorking and Jet Eruption in Gamma‐Ray Bursts

We show that the collapsar model of gamma-ray bursts results in a series of successive shocks and rarefaction waves propagating in the cork of stellar material being pushed ahead of the jet, as it emerges from the massive stellar progenitor. Our results are derived from analytical calculations assuming a hot, ultrarelativistic one-dimensional flow with an initial Lorentz factor Γj ~ 100. The shocks result in a series of characteristic, increasingly shorter and harder thermal X-ray pulses, as well as a nonthermal γ-ray pulse, which precede the usual nonthermal MeV γ-rays. We consider jets escaping from both compact (CO or He dwarf) and blue supergiant stellar progenitors. The period, fluence, and hardness of the pulses serves as a diagnostic for the size and density of the outer envelope of the progenitor star.

A Submillimeter and Radio Survey of Gamma‐Ray Burst Host Galaxies: A Glimpse into the Future of Star Formation Studies

We present the first comprehensive search for submillimeter and radio emission from the host galaxies of 20 well-localized gamma-ray bursts (GRBs). With the exception of a single source, all observations were undertaken months to years after the GRB explosions to ensure negligible contamination from the afterglows. We detect the host galaxy of GRB 000418 in both the submillimeter and radio, and the host galaxy of GRB 000210 in only the submillimeter. These observations, in conjunction with the previous detections of the host galaxies of GRB 980703 and GRB 010222, indicate that about 20% of GRB host galaxies are ultraluminous (L > 1012 L☉) and have star formation rates of about 500 M☉ yr-1. As an ensemble, the nondetected hosts have a star formation rate of about 100 M☉ yr-1 (5 σ) based on their radio emission. This, in conjunction with an average luminosity for the entire sample that is approximately 20% fainter than the local starburst galaxy Arp 220, indicates that GRB hosts probe a more representative population of star-forming galaxies than those uncovered in blank submillimeter and radio surveys. The detected and ensemble star formation rates exceed the values determined from various optical estimators by an order of magnitude, indicating significant dust obscuration. In the same vein, the ratio of bolometric dust luminosity to UV luminosity for the hosts detected in the submillimeter and radio bands ranges from about ~30 to 500 and follows the known trend of increasing obscuration with increasing bolometric luminosity. We also show that the GRB host sample as a whole, and the submillimeter- and radio-detected hosts individually, have significantly bluer R-K colors as compared with galaxies selected in the submillimeter and radio in the same redshift range. This possibly indicates that the stellar populations in the GRB hosts are on average younger, supporting the massive stellar progenitor scenario for GRBs, but it is also possible that GRB hosts are on average less dusty. For the nondetected GRB hosts, the difference in R-K color may also be a manifestation of their more representative bolometric luminosities relative to the highly luminous submillimeter- and radio-selected galaxies. Beyond the specific results presented in this paper, the submillimeter and radio observations serve as an observational proof-of-concept in anticipation of the upcoming launch of the Swift GRB mission and SIRTF. These new facilities will possibly bring GRB host galaxies into the forefront of star formation studies.

Association of supernovae with recent star formation regions in late type galaxies

A statistical study has been carried out to constrain the various models for the progenitors of supernovae of types Ia, Ib/c, and II. Star formation regions in late type galaxies are traced out by H II regions, imaged through photographic and CCD observations, and possible supernova association with these regions is based on the ratio of the angular separation of each supernova from its nearest H II region to the angular extent of the H II region in the direction of the supernova. The specific problems of supernova classification and positional uncertainties, as well as the probability of chance superposition, are also considered. The results suggest that type Ia supernovae do not arise from massive, short- lived stellar populations. Type Ib/c and type II supernovae, however, are very likely to be associated with H II regions and therefore with massive stellar progenitors. The single Wolf-Rayet star progenitor model for the type Ib/c supernovae is not supported. It is not clear from the results whether a difference between the association of type Ib/c supernovae and the association of type II supernovae with H II regions exists, but this probably results from the effects of small-number statistics for the type Ib/c supernovae. Not all type Ib/c and type II supernovae were found to be associated with detected H II regions; these supernovae either may have had runaway O star progenitors or were associated with H II regions below detectability.

The Galactic Pulsar Population and Neutron Star Birth

The radio pulsars in the Galaxy are found predominantly in the disk, with a scale height of several hundred parsecs. After allowing for pulsar velocities, the data are consistent with the hypothesis that single pulsars form from massive stellar progenitors. The number of active single pulsars in the Galaxy is ∼ 1.5 × 105, and their birthrate is 1 per ∼ 60 yrs. There is some evidence that many single pulsars, particularly those with high magnetic fields, are born spinning slowly, with initial periods ∼ 0.5–1s. This could imply an origin through binary “recycling” followed by orbit disruption, or might suggest that the pre-supernova stellar core efficiently loses angular momentum to the envelope through magnetic coupling. The birthrate of binary radio pulsars, particularly of the millisecond variety, seems to be much larger than previous estimates, and might suggest that these systems do not originate in low mass X-ray binary systems.