Stellar Collapse Research Articles

On the synthesis of heavy nuclei in protomagnetar outflows and implications for ultra-high energy cosmic rays

ABSTRACT It has been suggested that strongly magnetized and rapidly rotating protoneutron stars (PNSs) may produce long duration gamma-ray bursts (GRBs) originating from stellar core collapse. We explore the steady-state properties and heavy element nucleosynthesis in neutrino-driven winds from such PNSs whose magnetic axis is generally misaligned with the axis of rotation. We consider a wide variety of central engine properties such as surface dipole field strength, initial rotation period, and magnetic obliquity to show that heavy element nuclei can be synthesized in the radially expanding wind. This process is facilitated provided the outflow is Poynting-flux dominated such that its low entropy and fast expansion time-scale enables heavy nuclei to form in a more efficient manner as compared to the equivalent thermal GRB outflows. We also examine the acceleration and survival of these heavy nuclei and show that they can reach sufficiently high energies ≳ 1020 eV within the same physical regions that are also responsible for powering gamma-ray emission, primarily through magnetic dissipation processes. Although these magnetized outflows generally fail to achieve the production of elements heavier than lanthanides for our explored electron fraction range 0.4–0.6, we show that they are more than capable of synthesizing nuclei near and beyond iron peak elements.

The Effect of Supernova Convection On Neutron Star and Black Hole Masses

Our understanding of the convective-engine paradigm driving core-collapse supernovae has been used for two decades to predict the remnant mass distribution from stellar collapse. These predictions improve as our understanding of this engine increases. In this paper, we review our current understanding of convection (in particular, the growth rate of convection) in stellar collapse and study its effect on the remnant mass distribution. We show how the depth of the mass gap between neutron stars and black holes can help probe this convective growth. We include a study of the effects of stochasticity in both the stellar structure and the convective seeds caused by stellar burning. We study the role of rotation and its effect on the pair-instability mass gap. Under the paradigm limiting stellar rotation to those stars in tight binaries, we determine the effect of rotation on the remnant mass distribution.

Energy conserving and well-balanced discontinuous Galerkin methods for the Euler–Poisson equations in spherical symmetry

ABSTRACT This paper presents high-order Runge–Kutta (RK) discontinuous Galerkin methods for the Euler–Poisson equations in spherical symmetry. The scheme can preserve a general polytropic equilibrium state and achieve total energy conservation up to machine precision with carefully designed spatial and temporal discretizations. To achieve the well-balanced property, the numerical solutions are decomposed into equilibrium and fluctuation components that are treated differently in the source term approximation. One non-trivial challenge encountered in the procedure is the complexity of the equilibrium state, which is governed by the Lane–Emden equation. For total energy conservation, we present second- and third-order RK time discretization, where different source term approximations are introduced in each stage of the RK method to ensure the conservation of total energy. A carefully designed slope limiter for spherical symmetry is also introduced to eliminate oscillations near discontinuities while maintaining the well-balanced and total-energy-conserving properties. Extensive numerical examples – including a toy model of stellar core collapse with a phenomenological equation of state that results in core bounce and shock formation – are provided to demonstrate the desired properties of the proposed methods, including the well-balanced property, high-order accuracy, shock-capturing capability, and total energy conservation.

Search for intermediate-mass black hole binaries in the third observing run of Advanced LIGO and Advanced Virgo

Intermediate-mass black holes (IMBHs) span the approximate mass range 100−105 M⊙, between black holes (BHs) that formed by stellar collapse and the supermassive BHs at the centers of galaxies. Mergers of IMBH binaries are the most energetic gravitational-wave sources accessible by the terrestrial detector network. Searches of the first two observing runs of Advanced LIGO and Advanced Virgo did not yield any significant IMBH binary signals. In the third observing run (O3), the increased network sensitivity enabled the detection of GW190521, a signal consistent with a binary merger of mass ∼150 M⊙ providing direct evidence of IMBH formation. Here, we report on a dedicated search of O3 data for further IMBH binary mergers, combining both modeled (matched filter) and model-independent search methods. We find some marginal candidates, but none are sufficiently significant to indicate detection of further IMBH mergers. We quantify the sensitivity of the individual search methods and of the combined search using a suite of IMBH binary signals obtained via numerical relativity, including the effects of spins misaligned with the binary orbital axis, and present the resulting upper limits on astrophysical merger rates. Our most stringent limit is for equal mass and aligned spin BH binary of total mass 200 M⊙ and effective aligned spin 0.8 at 0.056 Gpc−3 yr−1 (90% confidence), a factor of 3.5 more constraining than previous LIGO-Virgo limits. We also update the estimated rate of mergers similar to GW190521 to 0.08 Gpc−3 yr−1.

Supernova Fallback as Origin of Neutron Star Spins and Spin-kick Alignment

Abstract Natal kicks and spins are characteristic properties of neutron stars (NSs) and black holes (BHs). Both offer valuable clues to dynamical processes during stellar core collapse and explosion. Moreover, they influence the evolution of stellar multiple systems and the gravitational-wave signals from their inspiral and merger. Observational evidence of a possibly generic spin-kick alignment has been interpreted as an indication that NS spins are either induced with the NS kicks or inherited from the progenitor rotation, which thus might play a dynamically important role during stellar collapse. Current three-dimensional supernova simulations suggest that NS kicks are transferred in the first seconds of the explosion, mainly by anisotropic mass ejection and, on a secondary level, anisotropic neutrino emission. By contrast, the NS spins are only determined minutes to hours later by the angular momentum associated with the fallback of matter that does not become gravitationally unbound in the supernova. Here, we propose a novel scenario to explain spin-kick alignment as a consequence of tangential vortex flows in the fallback matter that is accreted mostly from the direction of the NS’s motion. For this effect the initial NS kick is crucial, because it produces a growing offset of the NS away from the explosion center, thus promoting one-sided accretion. In this new scenario conclusions based on traditional concepts are reversed. For example, pre-kick NS spins are not required, and rapid progenitor core rotation can hamper spin-kick alignment. We also discuss implications for natal BH kicks and the possibility of tossing the BH’s spin axis during its formation.

Triangulating black hole forming stellar collapses through neutrinos

In the event of a black hole (BH) forming stellar collapse, the neutrino signal should terminate abruptly at the moment of BH formation, after a phase of steady accretion. Since neutrinos are expected to reach Earth hours before the electromagnetic signal, the combined detection of the neutrino burst through multiple neutrino telescopes could allow to promptly determine the angular location of a nearby stellar collapse in the sky with high precision. In this paper, we contrast the triangulation pointing procedure that relies on the rise time of the neutrino curve, often considered in the literature, to the one that takes advantage of the cutoff of the neutrino curve at the moment of BH formation. By forecasting the neutrino signal expected in the IceCube Neutrino Observatory, Hyper-Kamiokande and DUNE, we devise a strategy to optimize the identification of the rise and cutoff time of the neutrino curve. We show that the triangulation method developed by employing the end tail of the neutrino curve allows to achieve at least one order of magnitude improvement in the pointing precision for a galactic burst, while being insensitive to the neutrino mixing scenario. The triangulation pointing method based on the cutoff of the neutrino curve will also guarantee a better performance for BH forming collapses occurring beyond our own Galaxy.

Detailed Study of Stellar Evolution and Stellar Gravitational Collapse Leading to The Formation of Black Holes

In this paper we talk about the formation mechanism of the stellar gravitational singularities. To begin with we talk about the formation and evolution of stars and understand the Hertzsprung Russell diagram that teaches us how to classify stars. With the help of the diagram, we categorize the stars based of their physical parameters such as color, temperature, and mass. We then talk about the death cycle of different mass stars and what comes after when their fuels have been exhausted. Stars below the Chandrashekhar limit form a white dwarf at the end of their lives, while stars above the limit form a neutron star of a singularity. Further to find which of the heavy mass stars forms a singularity we look at the Tolman-Oppenheimer-Volkoff limit that states stars with mass above the limit will form singularities. The types of singularities formed depend on the solution of the general theory of relativity given by Schwarzschild, Kerr, Kerr-Newman, and Reissner-Nordström. Each of the theory aspect of the four solutions has been described to give a better understanding of the structure of the singularity formed. The paper also explains theories such as Wormholes and time travel in brief to try and explain what can replace the singularity.

Structure of two-and three-alpha systems in cold neutron matter

We present stability and structure of two-and three-alpha systems embedded in dilute cold neutron matter. By solving a few-alpha Schrödinger equation with quasiparticle properties, i.e.,effective mass and induced two-and three-alpha interactions, which are evaluated in terms of Fermi polarons,itis shown that 8Be and the Hoyle state become bound at densities of about10−4 and 10−3 of the saturation density, respectively. It is also seen that, under cold neutron matter environment, both systems become smaller than the corresponding systems in vacuum. Our results would affect astrophysical models for stellar collapse and neutron star mergers, as well as relevant reaction rates for nucleosynthesis.

Probes of the progenitors, engines and physics behind stellar collapse

Super-novas (SNs) are one of the most powerful explosions in the universe and astronomers have invoked the collapse of a stellar core down to a neutron star as a potential power source behind these cosmic blasts. The current paradigm behind core-collapse SN relies on convection in the region just above the newly formed neutron star. This engine was driven and confirmed by observations. We review this observational evidence, and the potential for further observational constraints in this paper.

Sensitivity of multi-PMT optical modules in Antarctic ice to supernova neutrinos of MeV energy

New optical sensors with a segmented photosensitive area are being developed for the next generation of neutrino telescopes at the South Pole. In addition to increasing sensitivity to high-energy astrophysical neutrinos, we show that this will also lead to a significant improvement in sensitivity to MeV neutrinos, such as those produced in core-collapse supernovae (CCSN). These low-energy neutrinos can provide a detailed picture of the events after stellar core collapse, testing our understanding of these violent explosions. We present studies on the event-based detection of MeV neutrinos with a segmented sensor and, for the first time, the potential of a corresponding detector in the deep ice at the South Pole for the detection of extra-galactic CCSN. We find that exploiting temporal coincidences between signals in different photocathode segments, a 27 mathrm {M}_{odot } progenitor mass CCSN can be detected up to a distance of 341 kpc with a false detection rate of {0.01},hbox {year}^{-1} with a detector consisting of 10,000 sensors. Increasing the number of sensors to 20,000 and reducing the optical background by a factor of 70 expands the range such that a CCSN detection rate of 0.1 per year is achieved, while keeping the false detection rate at {0.01},{hbox {year}^{-1}}.

Scalar gravitational wave signals from core collapse in massive scalar–tensor gravity with triple-scalar interactions

The spontaneous scalarization during the stellar core collapse in the massive scalar–tensor theories of gravity introduces extra polarizations (on top of the plus and cross modes) in gravitational waves, whose amplitudes are determined by several model parameters. Observations of such scalarization-induced gravitational waveforms therefore offer valuable probes into these theories of gravity. Considering a triple-scalar interactions in such theories, we find that the self-coupling effects suppress the magnitude of the scalarization and thus reduce the amplitude of the associated gravitational wave signals. In addition, the self-interacting effects in the gravitational waveform are shown to be negligible due to the dispersion throughout the astrophysically distant propagation. As a consequence, the gravitational waves observed on the Earth feature the characteristic inverse-chirp pattern. Although not with the on-going ground-based detectors, we illustrate that the scalarization-induced gravitational waves may be detectable at a signal-to-noise ratio level of with future detectors, such as Einstein Telescope and Cosmic Explorer.

Mapping the likelihood of GW190521 with diverse mass and spin priors

We map the likelihood of GW190521, the heaviest detected binary black hole (BBH) merger, by sampling under different mass and spin priors designed to be uninformative. We find that a source-frame total mass of $\sim$$150 M_{\odot}$ is consistently supported, but posteriors in mass ratio and spin depend critically on the choice of priors. We confirm that the likelihood has a multi-modal structure with peaks in regions of mass ratio representing very different astrophysical scenarios. The unequal-mass region ($m_2 / m_1 < 0.3$) has an average likelihood $\sim$$e^6$ times larger than the equal-mass region and a maximum likelihood $\sim$$e^2$ larger. Using ensembles of samples across priors, we examine the implications of qualitatively different BBH sources that fit the data. We find that the equal-mass solution has poorly constrained spins and at least one black hole mass that is difficult to form via stellar collapse due to (pulsational) pair instability. The unequal-mass solution can avoid this mass gap entirely but requires a negative effective spin and a precessing primary. Both of these scenarios are more easily produced by dynamical formation channels than field binary co-evolution. Drawing representative samples from each region of the likelihood map, we find a sensitive comoving volume-time $\mathcal{O}(10)$ times larger in the mass gap region than the gap-avoiding region. Accounting for this distance effect, the likelihood still reverses the advantage to favor the gap-avoiding scenario by a factor of $\mathcal{O}(100)$ before including mass and spin priors. Posterior samplers can be driven away from this high-likelihood region by common prior choices meant to be uninformative, making GW190521 parameter inference sensitive to the choice of mass and spin priors. This may be a generic issue for similarly heavy events given current detector sensitivity and waveform degeneracies.

Dynamical Formation of Scalarized Black Holes and Neutron Stars through Stellar Core Collapse.

In a certain class of scalar-Gauss-Bonnet gravity, the black holes and the neutron stars can undergo spontaneous scalarization-a strong gravity phase transition triggered by a tachyonic instability due to the nonminimal coupling between the scalar field and the spacetime curvature. Studies of this phenomenon have, so far, been restricted mainly to the study of the tachyonic instability and stationary scalarized black holes and neutron stars. To date, no realistic physical mechanism for the formation of isolated scalarized black holes and neutron stars has been proposed. We study, for the first time, the spherically symmetric fully nonlinear stellar core collapse to a black hole and a neutron star in scalar-Gauss-Bonnet theories allowing for a spontaneous scalarization. We show that the core collapse can produce scalarized black holes and scalarized neutron stars starting with a nonscalarized progenitor star. The possible paths to reach the end (non)scalarized state are quite rich leading to interesting possibilities for observational manifestations.

Detection of LIGO-Virgo binary black holes in the pair-instability mass gap

By probing the population of binary black hole (BBH) mergers detected by LIGO-Virgo, we can infer properties about the underlying black hole formation channels. A mechanism known as pair-instability (PI) supernova is expected to prevent the formation of black holes from stellar collapse with mass greater than $\ensuremath{\sim}40--65\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$ and less than $\ensuremath{\sim}120\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$. Any BBH merger detected by LIGO-Virgo with a component black hole in this so-called PI mass gap likely originated from an alternative formation channel. Here, we firmly establish GW190521 as an outlier to the stellar-mass BBH population if the PI mass gap begins at or below $65\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$. In addition, for a PI lower boundary of $40--50\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$, we find it unlikely that the remaining distribution of detected BBH events, excluding GW190521, is consistent with the stellar-mass BBH population.

A transient radio source consistent with a merger-triggered core collapse supernova.

A core collapse supernova occurs when exothermic fusion ceases in the core of a massive star, which is typically caused by exhaustion of nuclear fuel. Theory predicts that fusion could be interrupted earlier by merging of the star with a compact binary companion. We report a luminous radio transient, VT J121001+495647, found in the Very Large Array Sky Survey. The radio emission is consistent with supernova ejecta colliding with a dense shell of material, potentially ejected by binary interaction in the centuries before explosion. We associate the supernova with an archival x-ray transient, which implies that a relativistic jet was launched during the explosion. The combination of an early relativistic jet and late-time dense interaction is consistent with expectations for a merger-driven explosion.

Press–Schechter primordial black hole mass functions and their observational constraints

ABSTRACT We present a modification of the Press–Schechter (PS) formalism to derive general mass functions for primordial black holes (PBHs), considering their formation as being associated with the amplitude of linear energy density fluctuations. To accommodate a wide range of physical relations between the linear and non-linear conditions for collapse, we introduce an additional parameter to the PS mechanism, and that the collapse occurs at either a given cosmic time, or as fluctuations enter the horizon. We study the case where fluctuations obey Gaussian statistics and follow a primordial power spectrum of broken power-law form with a blue spectral index for small scales. We use the observed abundance of supermassive black holes (SMBH) to constrain the extended mass functions taking into account dynamical friction. We further constrain the modified PS by developing a method for converting existing constraints on the PBH mass fraction, derived assuming monochromatic mass distributions for PBHs, into constraints applicable for extended PBH mass functions. We find that when considering well-established monochromatic constraints, there are regions in parameter space where all the dark matter can be made of PBHs. Of special interest is the region for the characteristic mass of the distribution ${\sim}10^2\, \mathrm{M}_\odot$, for a wide range of blue spectral indices in the scenario where PBHs form as they enter the horizon, where the linear threshold for collapse is of the order of the typical overdensities, as this is close to the black hole masses detected by LIGO, which are difficult to explain by stellar collapse.

High Mass but Low Spin: An Exclusion Region to Rule Out Hierarchical Black Hole Mergers as a Mechanism to Populate the Pair-instability Mass Gap

Abstract The occurrence of pair-instability supernovae is predicted to prevent the formation of black holes with masses ≳50 M ⊙. Recent gravitational-wave detections in this mass range require an explanation beyond that of standard stellar collapse. Current modeling strategies include the hierarchical assembly of previous generations of black hole mergers as well as other mechanisms of astrophysical nature (lowered nuclear-reaction rates, envelope retention, stellar mergers, accretion, dredge-up episodes). In this paper, we point out the occurrence of an exclusion region that cannot be easily populated by hierarchical black hole mergers. A future gravitational-wave detection of a black hole with mass ≳50 M ⊙ and spin ≲0.2 will indicate that the pair-instability mass gap is polluted in some other way. Such a putative outlier can be explained using hierarchical mergers only with considerable fine-tuning of both mass ratio and spins of the preceding black hole merger—an assumption that can then be cross-checked against the bulk of the gravitational-wave catalog.

The neutrino gravitational memory from a core collapse supernova: phenomenology and physics potential

General Relativity predicts that the passage of matter or radiation from an asymmetrically-emitting source should cause a permanent change in the local space-time metric. This phenomenon, called the gravitational memory effect, has never been observed, however supernova neutrinos have long been considered a promising avenue for its detection in the future. With the advent of deci-Hertz gravitational wave interferometers, observing the supernova neutrino memory will be possible, with important implications for multimessenger astronomy and for tests of gravity. In this work, we develop a phenomenological (analytical) toy model for the supernova neutrino memory effect, which is overall consistent with the results of numerical simulations. This description is then generalized to several case studies of interest. We find that, for a galactic supernova, the dimensionless strain, h(t), is of order ∼ 10-22 - 10-21, and develops over a typical time scale that varies between ∼ 0.1 - 10 s, depending on the time-evolution of the anisotropy of the neutrino emission. The characteristic strain, h c(f), has a maximum at a frequency f max ∼ \U0001d4aa(10-1) - \U0001d4aa(1) Hz. The detailed features of the time- and frequency-structure of the memory strain will inform us of the matter dynamics near the collapsed core, and allow to distinguish between different stellar collapse scenarios. Next generation gravitational wave detectors like DECIGO and BBO will be sensitive to the neutrino memory effect for supernovae at typical galactic distances and beyond; with Ultimate DECIGO exceeding a detectability distance of 10 Mpc.

Revealing x-ray and gamma ray temporal and spectral similarities in the GRB 190829A afterglow.

Gamma-ray bursts (GRBs), which are bright flashes of gamma rays from extragalactic sources followed by fading afterglow emission, are associated with stellar core collapse events. We report the detection of very-high-energy (VHE) gamma rays from the afterglow of GRB 190829A, between 4 and 56 hours after the trigger, using the High Energy Stereoscopic System (H.E.S.S.). The low luminosity and redshift of GRB 190829A reduce both internal and external absorption, allowing determination of its intrinsic energy spectrum. Between energies of 0.18 and 3.3 tera-electron volts, this spectrum is described by a power law with photon index of 2.07 ± 0.09, similar to the x-ray spectrum. The x-ray and VHE gamma-ray light curves also show similar decay profiles. These similar characteristics in the x-ray and gamma-ray bands challenge GRB afterglow emission scenarios.

Determining the Structure of Rotating Massive Stellar Cores with Gravitational Waves

Abstract The gravitational wave (GW) signal resulting from stellar core collapse encodes a wealth of information about the physical parameters of the progenitor star and the resulting core-collapse supernova (CCSN). We present a novel approach to constrain CCSN progenitor properties at collapse using two of the most detectable parts of the GW signal: the core-bounce signal and evolution of the dominant frequency mode from the protoneutron star. We focus on the period after core bounce but before explosion and investigate the predictive power of GWs from rotating CCSNe to constrain properties of the progenitor star. We analyze 34 2D and four 3D neutrino-radiation-hydrodynamic simulations of stellar core collapse in progenitors of varied initial mass and rotation rate. Extending previous work, we verify the compactness of the progenitor at collapse to correlate with the early ramp-up slope, and in rotating cases, also with the core angular momentum. Combining this information with the bounce signal, we present a new analysis method to constrain the pre-collapse core compactness of the progenitor. Because these GW features occur less than a second after core bounce, this analysis could allow astronomers to predict electromagnetic properties of a resulting CCSN even before shock breakout.