High Energy Astrophysical Phenomena (astro-ph.HE)

Abstract: Type II supernovae represent the most common stellar explosions in the Universe, for which the final stage evolution of their hydrogen-rich massive progenitors towards core-collapse explosion are elusive. The recent explosion of SN 2023ixf in a very nearby galaxy, Messier 101, provides a rare opportunity to explore this longstanding issue. With the timely high-cadence flash spectra taken within 1-5 days after the explosion, we can put stringent constraints on the properties of the surrounding circumstellar material around this supernova. Based on the rapid fading of the narrow emission lines and luminosity/profile of $\rm H\alpha$ emission at very early times, we estimate that the progenitor of SN 2023ixf lost material at a mass-loss rate $\dot{\rm M} \approx 6 \times 10^{-4}\, \rm M_{\odot}\,yr^{-1}$ over the last 2-3 years before explosion. This close-by material, moving at a velocity $v_{\rm w} \approx 55\rm \, km\,s^{-1}$, accumulates a compact CSM shell at the radius smaller than $7 \times 10^{14}$ cm from the progenitor. Given the high mass-loss rate and relatively large wind velocity presented here, together with the pre-explosion observations made about two decades ago, the progenitor of SN 2023ixf could be a short-lived yellow hypergiant that evolved from a red supergiant shortly before the explosion.

Abstract: We perform a new general-relativistic viscous-radiation hydrodynamics simulation for supernova-like explosion associated with stellar core collapse of rotating massive stars to a system of a black hole and a massive torus paying particular attention to large-mass progenitor stars with the zero-age main-sequence mass of $M_\mathrm{ZAMS}=$20, 35, and 45$M_\odot$ of Ref.~\cite{Aguilera-Dena2020oct}. Assuming that a black hole is formed in a short timescale after the onset of the stellar collapse, the new simulations are started from initial data of a spinning black hole and infalling matter that self-consistently satisfy the constraint equations of general relativity. It is found that with a reasonable size of the viscous parameter, the supernova-like explosion is driven by the viscous heating effect in the torus around the black hole irrespective of the progenitor mass. The typical explosion energy and ejecta mass for the large-mass cases ($M_\mathrm{ZAMS}=35$ and $45M_\odot$) are $\sim 10^{52}$ erg and $\sim 5M_\odot$, respectively, with $^{56}$Ni mass larger than $0.15M_\odot$. These are consistent with the observational data of stripped-envelope and high-energy supernovae such as broad-lined type Ic supernovae. This indicates that rotating stellar collapses of massive stars to a black hole surrounded by a massive torus can be a central engine for high-energy supernovae. By artificially varying the angular velocity of the initial data, we explore the dependence of the explosion energy and ejecta mass on the initial angular momentum and find that the large explosion energy $\sim 10^{52}$ erg and large $^{56}$Ni mass $\geq 0.15M_\odot$ are possible only when a large-mass compact torus with mass $\gtrsim 1M_\odot$ is formed.

Abstract: Axion-like particles (ALPs) could mix with photons in the presence of astrophysical magnetic fields, and result in oscillations in the high energy $\gamma$-ray spectra observed by experiments. In this work, we investigate the ALP-photon oscillation effect through the blazar Mrk 421 spectra of 15 periods observed by Major Atmospheric Gamma Imaging Cherenkov Telescopes (MAGIC) and Fermi Large Area Telescope (Fermi-LAT). Compared with previous studies, we generate the mock data under the ALP hypothesis and apply the ${\rm CL_s}$ method to set constraints on the ALP parameters. This method is widely employed in high energy experiments and could avoid the possibility of excluding some parameter regions due to the fluctuation. We find that the ALP-photon coupling $g_{a\gamma}$ is constrained to be smaller than $\sim 2\times10^{-11}$ GeV$^{-1}$ for ALP mass ranging from $10^{-9}$ eV to $10^{-7}$ eV at a 95\% confidence level. The constraints obtained with the method based on the TS distribution under the null hypothesis, which is adopted in many previous astrophysical ALP studies, are also shown. Our results demonstrate that the joint constraints of all the periods from both methods are consistent. However, the latter method fails to provide constraints for some observation periods, whereas the ${\rm CL_s}$ method remains effective in such cases.

Abstract: Jets of active galactic nuclei are potential accelerators of ultra high-energy cosmic rays. Supernovae can occur inside these jets and contribute to cosmic ray acceleration, particularly of heavy nuclei, but that contribution has been hardly investigated so far. We carried out a first dedicated exploration of the role of supernovae inside extragalactic jets in the production of ultra high-energy cosmic rays. We characterized the energy budget of supernova-jet interactions, and the maximum possible energies of the particles accelerated in those events, likely dominated by heavy nuclei. This allowed us to assess whether these interactions can be potential acceleration sites of ultra high-energy cosmic rays, or at least of their seeds. For that, we estimated the cosmic ray luminosity for different galaxy types, and compared the injection rate of cosmic ray seeds into the jet with that due to galactic cosmic ray entrainment. Since the supernova is fueled for a long time by the luminosity of the jet, the energy of a supernova-jet interaction can be several orders of magnitude greater than that of an isolated supernova. Thus, despite the low rate of supernovae expected to occur in the jet, they could still provide more seeds for accelerating ultra high-energy particles than cosmic ray entrainment from the host galaxy. Moreover, these interactions can create sufficiently efficient accelerators to be a source of cosmic rays with energies $\gtrsim 10$~EeV. Supernova-jet interactions can contribute significantly to the production of ultra high-energy cosmic rays, either directly by accelerating these particles themselves or indirectly by providing pre-accelerated seeds.

Abstract: Based on the expected population of core collapse supernova remnants and the huge number of detected pulsars in the Galaxy, still representing only a fraction of the real population, pulsar wind nebulae are likely to constitute one of the largest classes of {extended} Galactic sources in many energy bands. For simple evolutionary reasons, the majority of the population is made of evolved systems, whose detection and identification are complicated by their reduced luminosity, the possible lack of X-ray emission (that fades progressively away with the age of the pulsar), and by their modified morphology with respect to young systems. Nevertheless they have gained renewed attention in recent years, following the detection of misaligned X-ray tails protruding from an increasing number of nebulae created by fast moving pulsars, and of extended TeV halos surrounding aged systems. Both these features are clear signs of an efficient escape of particles, with energy close to the maximum acceleration limit of the pulsar. Here we discuss the properties of those evolved systems and what we have understood about the process of particle escape, and the formation of observed features.

Abstract: Recent radio observations and coincident neutrino detections suggest that some tidal disruption events (TDEs) exhibit late-time activities, relative to the optical emission peak, and these may be due to delayed outflows launched from the central supermassive black hole. We investigate the possibility that jets launched with a time delay of days to months, interact with a debris that may expand outwards. We discuss the effects of the time delay and expansion velocity on the outcomes of jet breakout and collimation. We find that a jet with an isotropic-equivalent luminosity of $\lesssim 5 \times 10^{45}\,{\rm erg/s}$ is likely to be choked for a delay time of $\sim 3$ months. We also study the observational signatures of such delayed choked jets. The jet-debris interaction preceding the breakout would lead to particle acceleration and the resulting synchrotron emission can be detected by current and near-future radio, optical and X-ray telescopes, and the expanding jet-driven debris could explain late-time radio emission. We discuss high-energy neutrino production in delayed choked jets, and the time delay can significantly alleviate the difficulty of the hidden jet scenario in explaining neutrino coincidences.

Abstract: Pulse profiles of accreting millisecond pulsars can be used to determine neutron star (NS) parameters, such as their masses and radii, and therefore provide constraints on the equation of state of cold dense matter. Information obtained by the Imaging X-ray Polarimetry Explorer (IXPE) can be used to decipher pulsar inclination and magnetic obliquity, providing ever tighter constraints on other parameters. In this paper, we develop a new emission model for accretion-powered millisecond pulsars based on thermal Comptonization in an accretion shock above the NS surface. The shock structure was approximated by an isothermal plane-parallel slab and the Stokes parameters of the emergent radiation were computed as a function of the zenith angle and energy for different values of the electron temperature, the Thomson optical depth of the slab, and the temperature of the seed blackbody photons. We show that our Compton scattering model leads to a significantly lower polarization degree of the emitted radiation compared to the previously used Thomson scattering model. We computed a large grid of shock models, which can be combined with pulse profile modeling techniques both with and without polarization included. In this work, we used the relativistic rotating vector model for the oblate NS in order to produce the observed Stokes parameters as a function of the pulsar phase. Furthermore, we simulated the data to be produced by IXPE and obtained constraints on model parameters using nested sampling. The developed methods can also be used in the analysis of the data from future satellites, such as the enhanced X-ray Timing and Polarimetry mission.

Abstract: Measurements of neutron star masses, radii, and tidal deformability have direct connections to nuclear physics via the equation of state (EoS), which for the cold, catalyzed matter in neutron star cores is commonly represented as the pressure as a function of energy density. Microscopic models with exotic degrees of freedom display nontrivial structure in the speed of sound ($c_s$) in the form of first-order phase transitions and bumps, oscillations, and plateaus in the case of crossovers and higher-order phase transitions. We present a procedure based on Gaussian processes to generate an ensemble of EoSs that include nontrivial features. Using a Bayesian analysis incorporating measurements from X-ray sources, gravitational wave observations, and perturbative QCD results, we show that these features are compatible with current constraints. We investigate the possibility of a global maximum in $c_s$ that occurs within the densities realized in neutron stars -- implying a softening of the EoS and possibly an exotic phase in the core of massive stars -- and find that such a global maximum is consistent with, but not required by, current constraints.