Gamma-ray astronomy

Abstract:

In low-mass X-ray binaries (LMXBs), accretion flows are often associated with either jet outflows or disk winds. Studies of LMXBs with luminosities up to roughly 20% of the Eddington limit indicate that these outflows generally do not co-occur, suggesting that disk winds might inhibit jets. However, previous observations of LMXBs accreting near or above the Eddington limit show that jets and winds can potentially coexist. To investigate this phenomenon, we carried out a comprehensive multiwavelength campaign (using the Very Large Array (VLA), Chandra/High Energy Transmission Grating Spectrometer (HETG), and NICER) on the near-Eddington neutron-star Z-source LMXB GX 13+1. NICER and Chandra/HETG observations tracked GX 13+1 across the entire Z track during high Eddington rates, detecting substantial resonance absorption features originating from the accretion disk wind in all X-ray spectra, which implies a persistent wind presence. Simultaneous VLA observations captured a variable radio jet, with radio emission notably strong during all flaring branch observations—contrary to typical behavior in Z sources—and weaker when the source was on the normal branch. Interestingly, no clear correlation was found between the radio emission and the wind features. Analysis of VLA radio light curves and simultaneous Chandra/HETG spectra demonstrates that an ionized disk wind and jet outflow can indeed coexist in GX 13+1, suggesting that their launching mechanisms are not necessarily linked in this system.

Abstract:

Tidal disruption events (TDEs) can generate non-spherical, relativistic, and optically thick outflows. Simulations show that the radiation we observe is reprocessed by these outflows. According to a unified model suggested by these simulations, the spectral energy distributions (SEDs) of TDEs depend strongly on viewing angle: low [high] optical-to-X-ray ratios (OXRs) correspond to face-on [edge-on] orientations. Post-processing with radiative transfer codes has simulated the emergent spectra but has so far been carried out only in a quasi-1D framework, with three atomic species (H, He, and O). Here, we present 2.5D Monte Carlo radiative transfer simulations which model the emission from a non-spherical outflow, including a more comprehensive set of cosmically abundant species. While the basic trend of OXR increasing with inclination is preserved, the inherently multi-D nature of photon transport through the non-spherical outflow significantly affects the emergent SEDs. Relaxing the quasi-1D approximation allows photons to preferentially escape in (polar) directions of lower optical depth, resulting in a greater variation of bolometric luminosity as a function of inclination. According to our simulations, inclination alone may not fully explain the large dynamic range of observed TDE OXRs. We also find that including metals, other than O, changes the emergent spectra significantly, resulting in stronger absorption and emission lines in the extreme ultraviolet, as well as a greater variation in the OXR as a function of inclination. Whilst our results support previously proposed unified models for TDEs, they also highlight the critical importance of multi-D ionization and radiative transfer.

Abstract:

Supernova remnants (SNRs) are likely sources of hadronic particle acceleration within our galaxy, contributing to the galactic cosmic ray flux. Next-generation instruments, such as the Southern Wide-field Gamma-ray Observatory (SWGO), will be of crucial importance in identifying new candidate SNRs. SWGO will observe two-thirds of the gamma-ray sky, covering the energy range between a few hundreds of GeV and a PeV. In this work, we apply a model of SNR evolution to a catalogue of SNRs in order to predict their gamma-ray spectra, explore the SNR emission phase space, and quantify detection prospects for SWGO. Finally, we validate our model for sources observed with current-generation instruments, fitting it using a Monte-Carlo Markov Chain technique to the observed gamma-ray emission from four SNRs. We anticipate that at least 6, and potentially as many as 11 SNRs will be detected by SWGO within 1 year.

Abstract:

We present numerical simulations of discrete relativistic ejecta from an X-ray binary (XRB) with initial conditions directly informed by observations. XRBs have been observed to launch powerful discrete plasma ejecta during state transitions, which can propagate up to parsec distances. Understanding these ejection events unveils new understanding of jet-launching, jet power, and jet–interstellar medium (ISM) interaction among other implications. Multifrequency quasi-simultaneous radio observations of ejecta from the black hole XRB MAXI J1820+070 produced both size and calorimetry constraints, which we use as initial conditions of a relativistic hydrodynamic simulation. We qualitatively reproduce the observed deceleration of the ejecta in a homogeneous ISM. Our simulations demonstrate that the ejecta must be denser than the ISM, the ISM be significantly low density, and the launch be extremely powerful, in order to propagate to the observed distances. The blob propagates and clears out a high-pressure low-density cavity in its wake, providing an explanation for this pre-existing low-density environment, as well as ‘bubble-like’ environments in the vicinity of XRBs inferred from other studies. As the blob decelerates, we observe the onset of instabilities and a long-lived reverse shock – these mechanisms convert kinetic to internal energy in the blob, responsible for in situ particle acceleration. We transform the outputs of our simulation into pseudo-radio images, incorporating the coverage of the MeerKAT and e-MERLIN telescopes from the original observations with real-sky background. Through this, we maximize the interpretability of the results and provide direct comparison to current data, as well as provide prediction capabilities.