Search for: All records

Abstract We investigate the process of diffusive shock acceleration of particles with mass number to charge number ratiosA/Q > 1, e.g., partially ionized heavy ions. To this end, we introduce helium- and carbon-like ions at solar abundances into two-dimensional hybrid (kinetic ions–fluid electrons) simulations of nonrelativistic collisionless shocks. This study yields three main results: (1) Heavy ions are preferentially accelerated compared to hydrogen. For typical solar abundances, the energy transferred to accelerated helium ions is comparable to, or even exceeds, that of hydrogen, thereby enhancing the overall shock acceleration efficiency. (2) Accelerated helium ions contribute to magnetic field amplification, which increases the maximum attainable particle energy and steepens the spectra of accelerated particles. (3) The efficient acceleration of helium significantly enhances the production of hadronicγ-rays and neutrinos, likely dominating the one due to hydrogen. These effects should be taken into account, especially when modeling strong space and astrophysical shocks. 

Abstract Superheavy dark matter has been attractive as a candidate of particle dark matter. We propose a “natural” particle model, in which the dark matter serves as the inflaton in natural inflation, while decaying to high-energy particles at energies of 10⁹-10¹³GeV from the prediction of the inflation. A scalar field responsible for diluting the dark matter abundance revives the natural inflation either with or without the recent data from the Atacama Cosmology Telescope (ACT) and baryon acoustic oscillation results from Dark Energy Spectroscopic Instrument.Since the dark matter must be a spin-zero scalar, we carefully study the galactic dark matter 3-body decay into fermions and two body decays into a gluon pair, and point out relevant multi-messenger bounds that constrain these decay modes. Interestingly, the predicted energy scale may coincide with the AMATERASU event and/or the KM3NeT neutrino event, KM3-230213A. We also point out particle models with dark baryon to further alleviateγ-ray bounds. This scenario yields several testable predictions for the UHECR observations, including the highest-energy neutrons that are unaffected by magnetic fields, the tensor-to-scalar ratio, the running of spectral indices,α_s≳ 𝒪(0.001), and the existence of light new colored particles that could be accessible at future collider experiments.Further measurements of high-energy cosmic rays, including their components and detailed directions, may provide insight into not only the origin of the cosmic rays but also inflation. 

Abstract We investigate the external reverse shock (RS) region of relativistic jets as the origin of X-ray afterglows of jetted tidal disruption events (TDEs) that exhibit luminous jets accompanied by fast-declining nonthermal X-ray emissions. We model the dynamics of jet propagating within an external density medium, accounting for continuous energy injection driven by accretion activities. We compute the time-dependent synchrotron and inverse Compton emissions from the RS region. Our analysis demonstrates that the RS scenario can potentially explain the X-ray light curves and spectra of four jetted TDEs, namely, AT 2022cmc, Swift J1644, Swift J2058, and Swift J1112. Notably, the rapid steepening of the late-stage X-ray light curves can be attributed jointly to the jet break and cessation of the central engine as the accretion rate drops below the Eddington limit. Using parameters obtained from X-ray data fitting, we also discuss the prospects forγ-ray and neutrino detection. 

ABSTRACT The light curves of radioactive transients, such as supernovae and kilonovae, are powered by the decay of radioisotopes, which release high-energy leptons through $$\beta ^+$$ and $$\beta ^-$$ decays. These leptons deposit energy into the expanding ejecta. As the ejecta density decreases during expansion, the plasma becomes collisionless, with particle motion governed by electromagnetic forces. In such environments, strong or turbulent magnetic fields are thought to confine particles, though the origin of these fields and the confinement mechanism have remained unclear. Using fully kinetic particle-in-cell (PIC) simulations, we demonstrate that plasma instabilities can naturally confine high-energy leptons. These leptons generate magnetic fields through plasma streaming instabilities, even in the absence of pre-existing fields. The self-generated magnetic fields slow lepton diffusion, enabling confinement, and transferring energy to thermal electrons and ions. Our results naturally explain the positron trapping inferred from late-time observations of thermonuclear and core-collapse supernovae. Furthermore, they suggest potential implications for electron dynamics in the ejecta of kilonovae. We also estimate synchrotron radio luminosities from positrons for Type Ia supernovae and find that such emission could only be detectable with next-generation radio observatories from a Galactic or local-group supernova in an environment without any circumstellar material. 

Abstract We present the results from our extensive hard-to-soft X-ray (NuSTAR, Swift-XRT, XMM-Newton, Chandra) and meter-to-millimeter-wave radio (Giant Metrewave Radio Telescope, Very Large Array, NOEMA) monitoring campaign of the very nearby (d = 6.9 Mpc) Type II supernova (SN) 2023ixf spanning ≈4–165 days post-explosion. This unprecedented data set enables inferences on the explosion’s circumstellar medium (CSM) density and geometry. In particular, we find that the luminous X-ray emission is well modeled by thermal free–free radiation from the forward shock with rapidly decreasing photoelectric absorption with time. The radio spectrum is dominated by synchrotron radiation from the same shock. Similar to the X-rays, the level of free–free absorption affecting the radio spectrum rapidly decreases with time as a consequence of the shock propagation into the dense CSM. While the X-ray and the radio modeling independently support the presence of a dense medium corresponding to an effective mass-loss rate $\dot{M}\approx 10^{-4}{M}_{\odot }{\mathrm{yr}}^{-1}$atR = (0.4–14) × 10¹⁵cm (forv_w = 25 km s⁻¹), our study points at a complex CSM density structure with asymmetries and clumps. The inferred densities are ≈10–100 times those of typical red supergiants, indicating an extreme mass-loss phase of the progenitor in the ≈200 yr preceding core collapse, which leads to the most X-ray luminous Type II SN and the one with the most delayed emergence of radio emission. These results add to the picture of the complex mass-loss history of massive stars on the verge of collapse and demonstrate the need for panchromatic campaigns to fully map their intricate environments. 

Abstract AT 2022cmc is a recently documented tidal disruption event that exhibits a luminous jet, accompanied by fast-declining X-ray and long-lasting radio and millimeter emission. Motivated by the distinct spectral and temporal signatures between the X-ray and radio observations, we propose a multizone model involving relativistic jets with different Lorentz factors. We systematically study the evolution of faster and slower jets in an external density profile, considering the continuous energy injection rate associated with time-dependent accretion rates before and after the mass fallback time. We investigate time-dependent multiwavelength emission from both the forward shock (FS) and reverse shock (RS) regions of the fast and slow jets, in a self-consistent manner. Our analysis demonstrates that the energy injection rate can significantly impact the jet evolution and subsequently influence the lightcurves. We find that the X-ray spectra and lightcurves could be described by electron synchrotron emission from the RS of the faster jet, in which the late-time X-ray upper limits, extending to 400 days after the disruption, could be interpreted as a jet break. Meanwhile, the radio observations can be interpreted as a result of synchrotron emission from the FS region of the slower jet. We also discuss prospects for testing the model with current and future observations. 

Abstract The detection of high-energy neutrino signals from the nearby Seyfert galaxy NGC 1068 provides us with an opportunity to study nonthermal processes near the center of supermassive black holes. Using the IceCube and latest Fermi-LAT data, we present general multimessenger constraints on the energetics of cosmic rays and the size of neutrino emission regions. In the photohadronic scenario, the required cosmic-ray luminosity should be larger than ∼1%−10% of the Eddington luminosity and the emission radius should be ≲15R_Sin low-βplasma and ≲3R_Sin high-βplasma. The leptonic scenario overshoots the NuSTAR or Fermi-LAT data for any emission radii we consider, and the required gamma-ray luminosity is much larger than the Eddington luminosity. The beta-decay scenario also violates not only the energetics requirement but also gamma-ray constraints, especially when the Bethe–Heitler and photomeson production processes are consistently considered. Our results rule out the leptonic and beta-decay scenarios in a nearly model-independent manner and support hadronic mechanisms in magnetically powered coronae if NGC 1068 is a source of high-energy neutrinos. 

Abstract We present a survey of 1D kinetic particle-in-cell simulations of quasi-parallel nonrelativistic shocks to identify the environments favorable for electron acceleration. We explore an unprecedented range of shock speedsv_sh≈ 0.067–0.267c, Alfvén Mach numbers ${\mathcal{M}}_{\mathrm{A}}=5–40$, sonic Mach numbers ${\mathcal{M}}_{\mathrm{s}}=5–160$, as well as the proton-to-electron mass ratiosm_i/m_e= 16–1836. We find that high Alfvén Mach number shocks can channel a large fraction of their kinetic energy into nonthermal particles, self-sustaining magnetic turbulence and acceleration to larger and larger energies. The fraction of injected particles is ≲0.5% for electrons and ≈1% for protons, and the corresponding energy efficiencies are ≲2% and ≈10%, respectively. The extent of the nonthermal tail is sensitive to the Alfvén Mach number; when ${\mathcal{M}}_{\mathrm{A}}\lesssim 10$, the nonthermal electron distribution exhibits minimal growth beyond the average momentum of the downstream thermal protons, independently of the proton-to-electron mass ratio. Acceleration is slow for shocks with low sonic Mach numbers, yet nonthermal electrons still achieve momenta exceeding the downstream thermal proton momentum when the shock Alfvén Mach number is large enough. We provide simulation-based parameterizations of the transition from thermal to nonthermal distribution in the downstream (found at a momentum around $p_{\mathrm{i},\mathrm{e}}/m_{\mathrm{i}}{v}_{\mathrm{sh}}\approx 3\sqrt{m_{\mathrm{i},\mathrm{e}}/m_{\mathrm{i}}}$), as well as the ratio of nonthermal electron to proton number density. The results are applicable to many different environments and are important for modeling shock-powered nonthermal radiation. 

The origin of Galactic cosmic rays (CRs), particularly around the knee region (∼3 PeV), remains a major unanswered question. Recent observations by LHAASO suggest that the knee is shaped mainly by protons, with a transition to heavier elements at higher energies. Microquasars—compact jet-emitting sources—have emerged as possible PeV CR accelerators, especially after detections of ultrahigh-energy gamma rays from these systems. We propose that the observed proton spectrum (hard below a few PeV, steep beyond) arises from the reacceleration of sub-TeV Galactic CRs via shear acceleration in large-scale microquasar jet-cocoon structures. Our model also naturally explains the observed spectrum of energies around a few tens of PeV by summing up heavier nuclei contributions. Additionally, similar reacceleration processes in radio galaxies can contribute to ultrahigh-energy CRs, bridging Galactic and extragalactic origins. Combined with low-energy CRs from supernova remnants and galaxy clusters around the second knee region, this scenario could provide a unified explanation for CRs across the entire energy spectrum.