Gamma-ray bursts (GRBs) are the most energetic explosions in the Universe, and their afterglow emission provides an opportunity to probe the physics of relativistic shock waves in an extreme environment. Several key pieces for completing the picture of GRB afterglow physics are still missing, including jet properties, the emission mechanism, and particle acceleration. Here, we present a study of the afterglow emission of GRB 221009A, the most energetic GRB ever observed. Using optical, X-ray, and gamma-ray data up to approximately 2 days after the trigger, we trace the evolution of the multiwavelength spectrum and the physical parameters behind the emission process. The broadband spectrum is consistent with the synchrotron emission emitted by relativistic electrons with its index of p = 2.29 ± 0.02. We identify a break energy at keV and an exponential cutoff at GeV in the observed multiwavelength spectrum. The break energy increases in time from keV at 0.65 days to keV at 1.68 days, favoring a stellar-wind-like profile of the circumburst medium with k = 2.4 ± 0.1 as in ρ(r) ∝ r −k . The high-energy attenuation at around 0.4 to 4 GeV is attributed to the maximum of the particle acceleration in the relativistic shock wave. This study confirms that the synchrotron process can explain the multiwavelength afterglow emission and its evolution.

Understanding the nature and identity of dark matter is a key goal in the physics community. In the case that TeV-scale dark matter particles decay or annihilate into Standard Model particles, very-high-energy (VHE) gamma rays (greater than 100 GeV) will be present in the final state. The Very Energetic Radiation Imaging Telescope Array System (VERITAS) is an imaging atmospheric Cherenkov telescope array that can indirectly detect VHE gamma rays in an energy range of 100 GeV to >30 TeV. Dwarf spheroidal galaxies (dSphs) are ideal candidates in the search for dark matter due to their high dark matter content, high mass-to-light ratios, and their low gamma-ray fluxes from astrophysical processes. This study uses a legacy dataset of 638 h collected on 17 dSphs, built over 11 years with an observing strategy optimized according to the dark matter content of the targets. The study addresses a broad dark matter particle mass range, extending from 200 GeV to 30 PeV. In the absence of a detection, we set the upper limits on the dark matter velocity-weighted annihilation cross section.

A kilonova, the electromagnetic emission produced by compact binary mergers, is formed through a delicate interplay of physical processes, involving r-process nucleosynthesis and interactions between heavy elements and photons through radiative transfer. This complexity makes it difficult to achieve a comprehensive understanding of kilonova spectra. In this study, we aim to enhance our understanding and establish connections between physical parameters and observables through radiative-transfer simulations. Specifically, we investigate how ejecta temperature and element mass influence the resulting kilonova spectrum. For each species, the strength of its line features depends on these parameters, leading to the formation of a distinct region in the parameter space, dubbed the resonance island, where the line signature of that species is notably evident in the kilonova spectrum. We explore its origin and applications. Among explored r-process elements (31 ≤ Z ≤ 92), we find that four species—SrII, YII, BaII, and CeII—exhibit large and strong resonance islands, suggesting their significant contributions to kilonova spectra at specific wavelengths. In addition, we discuss potential challenges and future perspectives in observable heavy elements and their masses in the context of the resonance island.

A kilonova is a short-lived explosive event in the Universe, resulting from the merger of two compact objects. Despite its importance as a primary source of heavy elements through r-process nucleosynthesis, its nature is not well understood due to its rarity. In this work, we introduce a model that determines the density of a radially stratified relativistic ejecta. We apply the model to kilonova ejecta and explore several hypothesized velocity profiles as a function of the merger's ejection time. These velocity profiles result in diverse density profiles of the ejecta, for which we conduct radiative transfer simulations using TARDIS with the solar r-process composition. Consequently, we investigate the impact of the ejecta velocity profile on the resulting evolution of the lightcurve and spectra through the line transitions of heavy elements. The change in the rate at which these elements accumulate in the line-forming region leaves its imprint on the kilonova lightcurve at specific wavelengths, causing the lightcurves to decay at different rates. Furthermore, in several profiles, plateau-like behaviors (slow and/or flat decline) are also observed. In conclusion, this work proposes potential scenarios of the evolution of kilonova due to the ejecta velocity profile.

We perform a detailed analysis of broad pulses in bright gamma-ray bursts (GRBs) to understand the evolution of GRB broad pulses. Using the temporal and spectral properties, we test the high-latitude emission (HLE) scenario in the decaying phase of broad pulses. The HLE originates from the curvature effect of a relativistic spherical jet, where higher-latitude photons are delayed and softer than the observer's line-of-sight emission. The signature of HLE has not yet been identified undisputedly during the prompt emission of GRBs. The HLE theory predicts a specific relation, Fν,Ep ∝ E2p, between the peak energy E p in ν F ν spectra and the spectral flux F ν measured at E p , Fν,Ep. We search for evidence of this relation in 2157 GRBs detected by the Gamma-ray Burst Monitor on board the Fermi Gamma-ray Space Telescope from 2008 to 2017. After imposing unbiased selection criteria in order to minimize contamination in a signal by background and overlaps of pulses, we build a sample of 32 broad pulses in 32 GRBs. We perform a time-resolved spectral analysis on each of these 32 broad pulses and find that the evolution of 18 pulses (56%) is clearly consistent with the HLE relation. For the 18 broad pulses, the exponent δ in the relation of Fν,Ep ∝ E p δ is distributed as a Gaussian function with a median and width of 1.99 and 0.34, respectively. This result provides a constraint on the emission radius of GRBs with the HLE signature.

Dark matter is a key piece of the current cosmological scenario, with weakly interacting massive particles (WIMPs) a leading dark matter candidate. WIMPs have not been detected in their conventional parameter space (100 GeV ≲M χ ≲ 100 TeV), a mass range accessible with current Imaging Atmospheric Cherenkov Telescopes. As ultraheavy dark matter (UHDM; M χ ≳ 100 TeV) has been suggested as an underexplored alternative to the WIMP paradigm, we search for an indirect dark matter annihilation signal in a higher mass range (up to 30 PeV) with the VERITAS γ-ray observatory. With 216 hr of observations of four dwarf spheroidal galaxies, we perform an unbinned likelihood analysis. We find no evidence of a γ-ray signal from UHDM annihilation above the background fluctuation for any individual dwarf galaxy nor for a joint-fit analysis, and consequently constrain the velocity-weighted annihilation cross section of UHDM for dark matter particle masses between 1 TeV and 30 PeV. We additionally set constraints on the allowed radius of a composite UHDM particle.

For decades, searches for electroweak-scale dark matter (DM) have been performed without a definitive detection. This lack of success may hint that DM searches have focused on the wrong mass range. A proposed candidate beyond the canonical parameter space is ultraheavy DM (UHDM). In this work, we consider indirect UHDM annihilation searches for masses between 30 TeV and 30 PeV-extending well beyond the unitarity limit at ~100 TeV-and discuss the basic requirements for DM models in this regime. We explore the feasibility of detecting the annihilation signature, and the expected reach for UHDM with current and future very-high-energy (VHE; >100 GeV) γ-ray observatories. Specifically, we focus on three reference instruments: two Imaging Atmospheric Cherenkov Telescope arrays, modeled on VERITAS and CTA-North, and one extended air shower array, motivated by HAWC. With reasonable assumptions on the instrument response functions and background rate, we find a set of UHDM parameters (mass and cross section) for which a γ-ray signature can be detected by the aforementioned observatories. We further compute the expected upper limits for each experiment. With realistic exposure times, the three instruments can probe DM across a wide mass range. At the lower end, it can still have a point-like cross section, while at higher masses the DM could have a geometric cross section, indicative of compositeness.

We report on the observations of gamma-ray burst (GRB) 190114C by the Fermi Gamma-ray Space Telescope and the Neil Gehrels Swift Observatory. The prompt gamma-ray emission was detected by the Fermi GRB Monitor (GBM), the Fermi Large Area Telescope (LAT), and the Swift Burst Alert Telescope (BAT) and the long-lived afterglow emission was subsequently observed by the GBM, LAT, Swift X-ray Telescope (XRT), and Swift UV Optical Telescope. The early-time observations reveal multiple emission components that evolve independently, with a delayed power-law component that exhibits significant spectral attenuation above 40 MeV in the first few seconds of the burst. This power-law component transitions to a harder spectrum that is consistent with the afterglow emission observed by the XRT at later times. This afterglow component is clearly identifiable in the GBM and BAT light curves as a slowly fading emission component on which the rest of the prompt emission is superimposed. As a result, we are able to observe the transition from internal-shock- to external-shock-dominated emission. We find that the temporal and spectral evolution of the broadband afterglow emission can be well modeled as synchrotron emission from a forward shock propagating into a wind-like circumstellar environment. We estimate the initial bulk Lorentz factor using the observed high-energy spectral cutoff. Considering the onset of the afterglow component, we constrain the deceleration radius at which this forward shock begins to radiate in order to estimate the maximum synchrotron energy as a function of time. We find that even in the LAT energy range, there exist high-energy photons that are in tension with the theoretical maximum energy that can be achieved through synchrotron emission from a shock. These violations of the maximum synchrotron energy are further compounded by the detection of very high-energy (VHE) emission above 300 GeV by MAGIC concurrent with our observations. We conclude that the observations of VHE photons from GRB 190114C necessitates either an additional emission mechanism at very high energies that is hidden in the synchrotron component in the LAT energy range, an acceleration mechanism that imparts energy to the particles at a rate that is faster than the electron synchrotron energy-loss rate, or revisions of the fundamental assumptions used in estimating the maximum photon energy attainable through the synchrotron process.

GRB 131108A is a bright long gamma-ray burst (GRB) detected by the Large Area Telescope and the Gamma-ray Burst Monitor on board the Fermi Gamma-ray Space Telescope. Dedicated temporal and spectral analyses reveal three γ-ray flares dominating above 100 MeV, which are not directly related to the prompt emission in the Gamma-ray Burst Monitor band (10 keV-10 MeV). The high-energy light curve of GRB 131108A (100 MeV-10 GeV) shows an unusual evolution: a steep decay, followed by three flares with an underlying emission, and then a long-lasting decay phase. The detailed analysis of the γ-ray flares finds that the three flares are 6-20 times brighter than the underlying emission and are similar to each other. The fluence of each flare, (1.6 ∼ 2.0) × 10-6 erg cm-2, is comparable to that of emission during the steep decay phase, 1.7 × 10-6 erg cm-2. The total fluence from three γ-ray flares is 5.3 × 10-6 erg cm-2. The three γ-ray flares show properties similar to the usual X-ray flares that are sharp flux increases, occurring in ∼50% of afterglows, in some cases well after the prompt emission. Also, the temporal and spectral indices during the early steep decay phase and the decaying phase of each flare show the consistency with a relation of the curvature effect (\hat{α } = 2 + \hat{β }), which is the first observational evidence of the high-latitude emission in the GeV energy band.

The synchrotron external shock model predicts the evolution of the spectral (β) and temporal (α) indices during the gamma-ray burst (GRB) afterglow for different environmental density profiles, electron spectral indices, electron cooling regimes, and regions of the spectrum. We study the relationship between α and β, the so-called “closure relations” with GRBs detected by Fermi Large Area Telescope (Fermi-LAT) from 2008 August to 2018 August. The spectral and temporal indices for the >100 MeV emission from the Fermi-LAT as determined in the Second Fermi-LAT Gamma-Ray Burst Catalog (2FLGC) are used in this work. We select GRBs whose spectral and temporal indices are well constrained (58 long-duration GRBs and 1 short-duration GRBs) and classify each GRB into the best-matched relation. As a result, we found that a number of GRBs require a very small fraction of the total energy density contained in the magnetic field (∊ B ≲ 10-7). The estimated mean and standard deviation of electron spectral index p are 2.40 and 0.44, respectively. The GRBs satisfying a closure relation of the slow cooling tend to have a softer p value compared to those of the fast cooling. Moreover, the Kolmogorov-Smirnov test of the two p distributions from the fast and slow coolings rejects a hypothesis that the two distributions are drawn from the single reference distribution with a significance of 3.2σ. Lastly, the uniform density medium is preferred over the medium that decreases like the inverse of distance squared for long-duration GRBs.

GRB 160709A is one of the few bright short gamma-ray bursts detected by both the Gamma-ray Burst Monitor and the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. The γ-ray prompt emission of GRB 160709A is adequately fitted by combinations of three distinct components: (i) a nonthermal component described by a power law (PL) with a high-energy exponential cutoff, (ii) a thermal component modeled with a Planck function, and (iii) a second nonthermal component shaped by an additional PL crossing the whole γ-ray spectrum. While the thermal component dominates during ∼0.12 s of the main emission episode of GRB 160709A with an unusually high temperature of ∼340 keV, the nonthermal components dominate in the early and late time. The thermal component is consistent with the photospheric emission resulting in the following parameters: the size of the central engine, R0=3.8-1.8+5.9 × 108 cm, the size of the photosphere, Rph=7.4-1.2+0.8 × 1010 cm, and a bulk Lorentz factor, Γ =728-93+75, assuming a redshift of 1. The slope of the additional PL spectrum stays unchanged throughout the burst duration; however, its flux decreases continuously as a function of time. A standard external shock model has been tested for the additional PL component using the relation between the temporal and spectral indices (the closure relation). Each set of spectral and temporal indices from two energy bands (200 keV-40 MeV and 100 MeV-10 GeV) satisfies a distinct closure relation. From the closure relation test we derived the index for the electron spectral distribution, p = 2.5 ± 0.1. The interaction of the jet with the interstellar environment is preferred over the interaction with the wind medium.