Search

In this paper, we studied the X-ray polarization of the jet in the Vela Pulsar Wind Nebula. After the background subtraction, we found a high degree of polarization in the helical jet region, close to the upper limit allowed by the synchrotron radiation mechanism. This result indicates the magnetic field is highly ordered, even though the jet's morphology is curved and varies rapidly, which suggests that the magnetic collimation process for the jet does not introduce significant turbulence. The polarization angle is perpendicular to the symmetry axis of the nebula, consistent with the results from the magnetohydrodynamic numerical simulations. This finding sheds light on the underlying mechanisms related to the dynamics and evolution of the jet structures. [ABSTRACT FROM AUTHOR]

By applying the Efron–Petrosian method to the fluxes S and distances D of the magnetars listed in the Magnetar Outburst Online Catalogue, we show that the observational data are consistent with the dependence |$S\propto D^{-3/2}$|⁠ , which characterizes the emission from the superluminally moving current sheet in the magnetosphere of a non-aligned neutron star, at substantially higher levels of significance than they are with the dependence |$S\propto D^{-2}$|⁠. This result agrees with that previously obtained by an analysis of the data in the McGill Online Magnetar Catalog and confirms that, contrary to the currently prevalent view, magnetars' X-ray luminosities do not exceed their spin-down luminosities. The X-ray spectra of magnetars, moreover, are congruous with the spectral energy distribution (SED) of a broad-band non-thermal emission mechanism identical to that at play in rotation-powered pulsars: we show that the SED of the caustics that are generated in certain privileged directions by the magnetospheric current sheet single-handedly fits the observed spectra of 4U 0142+61, 1E 1841−045, and XTE J1810−197 over their entire breadths. Magnetars' outbursts and their associated radio bursts are predicted to occur when, as a result of large-scale timing anomalies (such as glitches, quakes, or precession), one of the privileged directions along which the radiation from the current sheet decays more slowly than predicted by the inverse-square law either swings past or oscillates across the line of sight. [ABSTRACT FROM AUTHOR]

Energetic GeV photons expected from the closest and the most energetic Gamma-ray bursts (GRBs) provide a unique opportunity to study the very-high-energy emission as well as the possible correlations with lower energy bands in realistic GRB afterglow models. In the standard GRB afterglow model, the relativistic homogeneous shock is usually considered to be fully adiabatic, however, it could be partially radiative. Based on the external forward-shock scenario in both stellar wind and constant-density medium, we present a radiative-adiabatic analytical model of the synchrotron self-Compton (SSC) and synchrotron processes considering an electron energy distribution with a power-law index of |$1\lt p\lt 2$| and |$2\le p$|⁠. We show that the SSC scenario plays a relevant role in the radiative parameter |$\epsilon$|⁠ , leading to a prolonged evolution during the slow cooling regime. In a particular case, we derive the Fermi /LAT light curves together with the photons with energies |$\ge 100$|  MeV in a sample of nine bursts from the second Fermi /LAT GRB catalogue that exhibited temporal and spectral indices with |$\gtrsim 1.5$| and |$\approx 2$|⁠ , respectively. These events can hardly be described with closure relations of the standard synchrotron afterglow model, and also exhibit energetic photons above the synchrotron limit. We have modelled the multiwavelength observations of our sample to constrain the microphysical parameters, the circumburst density, the bulk Lorentz factor, and the mechanism responsible for explaining the energetic GeV photons. [ABSTRACT FROM AUTHOR]

The variability data for the BL Lacertae object TXS 1902+556 in the optical and |$\gamma$| -ray wavebands were obtained from the 0.76-m Katzman Automatic Imaging Telescope and the Fermi Large Area Telescope (Fermi -LAT), covering periods of 14.4 and 14.7 yr, respectively. The variability properties were systematically analysed, with particular emphasis on the first comprehensive investigation of radiation variation in the optical waveband. Four well-established techniques were employed for this purpose: the Lomb–Scargle periodogram, REDFIT program, Jurkevich method, and discrete correlation function (DCF) approach. The optical waveband exhibits quasi-periodic oscillations (QPO) with a time-scale of |$P_{\rm O}=276.8\pm 6.1$|  d at a significance level |$3.87\sigma$|⁠ , while the |$\gamma$| -ray waveband does not exhibit any significant periodicity. However, it should be noted that the QPO time-scale is consistent with the Sun-gaps in the optical light curve within 2 |$\sigma$| uncertainties. The optical QPO behaviour is most likely attributed to the helical motion of the jet driven by the orbital motion in a supermassive black hole binary system. Moreover, we have provided an explanation for the absence of QPO in the |$\gamma$| -ray light curves. Furthermore, utilizing the DCF method, a weak correlation between the variability in the optical and |$\gamma$| -ray wavebands was observed, suggesting that the emission of TXS 1902+556 may be generated through a combination of synchrotron self-Compton (SSC) and external Compton (EC) processes, or a leptonic–hadronic hybrid process. [ABSTRACT FROM AUTHOR]

Spider pulsars are compact binary systems composed of a millisecond pulsar and a low-mass companion. Their X-ray emission – modulated on the orbital period – is interpreted as synchrotron radiation from high-energy electrons accelerated at the intrabinary shock. We perform global two-dimensional particle-in-cell simulations of the intrabinary shock, assuming that the shock wraps around the companion star. When the pulsar spin axis is nearly aligned with the orbital angular momentum, we find that the magnetic energy of the relativistic pulsar wind – composed of magnetic stripes of alternating field polarity – efficiently converts to particle energy at the intrabinary shock, via shock-driven reconnection. The highest energy particles accelerated by reconnection can stream ahead of the shock and be further accelerated by the upstream motional electric field. In the downstream, further energization is governed by stochastic interactions with the plasmoids/magnetic islands generated by reconnection. We also extend our earlier work by performing simulations that have a larger (and more realistic) companion size and a more strongly magnetized pulsar wind. We confirm that our first-principles synchrotron spectra and light curves are in good agreement with X-ray observations. [ABSTRACT FROM AUTHOR]

We present a comprehensive analysis of the X-ray spectral properties of 198 newly identified active galactic nuclei (AGNs), leveraging archival data from the Chandra X-ray Observatory. All these AGNs exhibit a power-law spectral signature spanning a broad energy range of |$0.5{\!-\!}7.0$| keV, characterized by the photon index (⁠|$\Gamma$|⁠) values ranging from |$0.3^{+0.16}_{-0.14}$| to |$2.54^{+0.14}_{-0.13}$|⁠. Particularly, 76 of these AGNs display discernible levels of intrinsic absorption, after considering the Galactic absorption. The column densities associated with this local absorption (⁠|$n_{\rm H}^{\rm local}$|⁠) are within a range of |${\sim} 10^{19}{\!-\!}10^{22}\ {\rm cm^{-2}}$|⁠. We study the cosmological evolution of AGNs using the variation of |$n_{\rm H}^{\rm local}$| and |$\Gamma$| with their estimated redshift. The intrinsic spectral signature did not reveal any significant cosmological evolution; however, a deficit of hard sources at high redshift is possibly intrinsic. Our sample covers several orders of broad-band intrinsic luminosity (⁠|$L_{\rm B}^{\rm intr}$|⁠) ranging from |$4.59^{+0.41}_{-0.41} \times 10^{42}$| to |$2.4^{+0.12}_{-0.12} \times 10^{46}\, {\rm erg~s}^{-1}$| with peak at 1.84 redshift. We also investigate the hardness–luminosity diagram to further probe the AGNs. We conduct a sanity check by applying our findings to known AGNs, and the results are consistent with our observations. [ABSTRACT FROM AUTHOR]

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}^{-1}$| 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 FROM AUTHOR]

Context. X-shaped radio galaxies (XRGs) are characterised by two pairs of misaligned lobes: active lobes hosting radio jets and the wings. None of the formation mechanisms proposed thus far are able to exhaustively reproduce the diverse features observed among XRGs. Emerging evidence has proposed the existence of sub-populations of XRGs forming via different processes. Aims. The brightest cluster galaxy (BCG) in Abell 3670 (A3670) is a dumbbell system hosting the XRG MRC 2011-298. The morphological and spectral properties of this interesting XRG were first characterised based on Karl G. Jansky Very Large Array (JVLA) data at 1–10 GHz. In the present work, we follow up on MRC 2011-298 with the upgraded Giant Metrewave Radio Telescope (uGMRT) at 120–800 MHz to further constrain its properties and origin. Methods. We carried out a detailed spectral analysis sampling different spatial scales. Integrated radio spectra, spectral index maps, radio colour-colour diagrams, and radiative age maps of both the active lobes and prominent wings were employed to test the origin of the source. Results. We confirm a progressive spectral steepening from the lobes to the wings. The maximum radiative age of the source is ~80 Myr, with the wings being older than the lobes by ≳30 Myr in their outermost regions. Conclusions. The observed properties are in line with an abrupt reorientation of the jets by ~90 deg from the direction of the wings to their present position. This formation mechanism is further supported by the comparison with numerical simulations in the literature, which additionally highlight the role of hydrodynamic processes in the evolution of large wings such as those of MRC 2011-298. It is plausible that the coalescence of supermassive black holes could have triggered the spin-flip of the jets. Moreover, we show that the S-shape of the radio jets is likely driven by precession with a period of P ~ 10 Myr. [ABSTRACT FROM AUTHOR]

Context. Tailed radio galaxies are shaped by ram pressure that is due to the high-velocity motion of their host through the intracluster medium (ICM). Recent works have reported on the increasing complexity of the phenomenology of tailed galaxies, with departures from theoretical ageing models and novel evidence of re-energising mechanisms that are nonetheless unclear. Aims. The nearby (z = 0.0894) galaxy cluster Abell 2142 hosts two tailed galaxies, namely T1 and T2, which exhibit peculiar morphological features. We aim to investigate the properties of T1 and T2 and constrain their spectral evolution, dynamics, and interactions with the ICM. Methods. We combined data from LOw Frequency Array (LOFAR), upgraded Giant Metrewave Radio Telescope (uGMRT), Very Large Array (VLA), and MeerKAT (from 30 MHz to 6.5 GHz) to carry out a detailed spectral analysis of T1 and T2. We analysed the surface brightness profiles, measured integrated and spatially resolved spectral indices. We performed a comparison with single injection ageing models. The Chandra X-ray data were used to search for discontinuities in the ICM properties in the direction of the targets. Results. The spectral properties of T1 at low frequencies can be predicted by ageing models and provide constraints on the 3D dynamics of the host by assuming a constant velocity. However, a more complex scenario is suggested by the sharp transitions along sub-regions of the tail, local surface brightness enhancements, and a spectral shape at high frequencies that is not predicted by any of the models. This scenario may plausibly involve hydrodynamical instabilities and particle mixing. T2 exhibits unusual morphological and surface brightness features, and its spectral behaviour is not predicted by standard models. The effects of two active galactic nucleus (AGN) outburst events during the infall of T2 towards the cluster centre could potentially explain its characteristic properties. [ABSTRACT FROM AUTHOR]

Context. Abell 754 is a rich galaxy cluster at z = 0.0543 and is considered the prototype of a major cluster merger. As many dynamically unrelaxed systems, it hosts diffuse radio emission on megaparsec-scales. Extended synchrotron sources in the intra-cluster medium (ICM) are commonly interpreted as evidence that a fraction of the gravitational energy released during cluster mergers is dissipated into nonthermal components. Aims. Here, we aim to use new MeerKAT UHF- and L-band observations to study nonthermal phenomena in Abell 754. These data are complemented with archival XMM-Newton observations to investigate the resolved spectral properties of both the radio and X-ray cluster emission. Methods. For the first time, we employed the pipeline originally developed to calibrate LOFAR data to MeerKAT observations. This allowed us to perform a direction-dependent calibration and obtain highly sensitive radio images in UHF and L bands that capture the extended emission with unprecedented detail. By using a large XMM-Newton mosaic, we produced thermodynamic maps of the ICM. Results. Our analysis reveals that the radio halo in the cluster center is bounded by the well-known shock in the eastern direction. Furthermore, in the southwest periphery, we discover an extended radio source that we classify as a radio relic that is possibly tracing a shock driven by the squeezed gas compressed by the merger, outflowing in perpendicular directions. The low-luminosity of this relic appears compatible with direct acceleration of thermal pool electrons. We interpret the observed radio and X-ray features in the context of a major cluster merger with a nonzero impact parameter. Conclusions. Abell 754 is a remarkable galaxy cluster showcasing exceptional features associated with the ongoing merger event. The high quality of the new MeerKAT data motivates further work on this system. [ABSTRACT FROM AUTHOR]

Context. Massive stars produce strong stellar winds that consist of continuous outflows of material at speeds of thousands of km s−1. These winds convey large amounts of kinetic power, especially in the case of Wolf–Rayet (WR) stars. When these winds interact with nearby material, they will likely produce shocks. Among other processes, particle acceleration is expected to occur. This is particularly well established in the case of massive binary systems, where the stellar winds collide, allowing these systems to be identified thanks to the detection of synchrotron radio emission, produced by a population of relativistic particles accelerated in the shocks. Aims. Our goal is to investigate the occurrence of particle acceleration among massive stars in their pre-supernova evolution phases. Methods. We observed a subset of five WR stars in the radio domain using the upgraded Giant Metrewave Radio Telescope (uGMRT), located in India. The observations were carried out in bands 4 (550–950 MHz) and 5 (1050–1450 MHz) for all the targets. Results. We detected radio emission for only WR 110 in bands 4 and 5. Its thermal spectrum displays a consistent index of +0.74 down to uGMRT bands. The four other targets were not detected and we derived 3σ upper limits on their flux density. Our upper limits in Band 4 are the first provided for these targets below 1 GHz. None of the targets was identified as a synchrotron radio emitter in these radio bands. Conclusions. If some synchrotron emission is produced in these systems, the non-detection with the uGMRT can be most likely attributed to strong free-free absorption (FFA). This is especially relevant for WR98a which is catalogued as a particle accelerator based on previous measurements at higher radio frequencies. The prominence of FFA constitutes a severe obstacle to identifying particle accelerators in the radio domain. [ABSTRACT FROM AUTHOR]

Context. The evolution of active galactic nuclei (AGNs) is closely connected to their host galaxies and surroundings. Via feedback processes, AGNs can counteract the cooling of the intracluster medium (ICM) and suppress star formation in their host galaxies. Radio observations at low frequencies provide a glimpse into the history of AGN activity. The Virgo cluster is a substantial reservoir of nearby galaxies and provides an ideal laboratory for the study of AGNs as well as their feedback mechanisms. Aims. The aim of our work is to characterise the AGN population within the Virgo cluster down to low radio luminosities, constrain the AGN duty cycle, and investigate environmental feedback in cluster member galaxies. Methods. We analysed 144 MHz and 1.3 GHz radio observations of early-type galaxies from the ACS Virgo Cluster Survey (ACSVCS) taken with LOFAR and MeerKAT. Results. We detect 12 of these galaxies at 144 MHz, five of which show clearly extended radio emission. The radio luminosity shows a strong dependence on the stellar mass of the host galaxy, in agreement with previous results. As a notable outlier, the massive elliptical galaxy NGC 4365 (M* = 2.2 × 1011 M⊙) is not detected as compact source in the LOFAR observations. Instead, it is surrounded by diffuse, low-surface brightness emission, which hints towards a past phase of stronger nuclear activity. Furthermore, we find a cavity in NGC 4472 (=M 49) inflated by the wide-angle tail only visible in the LOFAR data, which implies that the cavity was created by a past outburst. The corresponding cavity power is of the same order of magnitude as the jet power in the present duty cycle of the AGN. [ABSTRACT FROM AUTHOR]

Context. The particle-in-cell approach has proven effective in modeling neutron-star and black-hole magnetospheres from first principles, but global simulations are plagued with an unrealistically small separation between the scales where microphysics operates and the system-size scales due to limited numerical resources. A legitimate concern is whether the scale separation achieved to date is large enough for results to be safely extrapolated to realistic scales. Aims. In this work, our aim is to explore the effect of scaling up physical parameters and to check whether salient features uncovered by pure kinetic models at smaller scales are still valid, with a special emphasis on particle acceleration and high-energy radiation emitted beyond the light cylinder. Methods. To reach this objective, we developed a new hybrid numerical scheme coupling the ideal force-free and the particle-in-cell methods to optimize the numerical cost of global models. We propose a domain decomposition of the magnetosphere based on the magnetic-field topology using the flux function. The force-free model is enforced along open field lines while the particle-in-cell model is restricted to the reconnecting field line region. Results. As a proof of concept, this new hybrid model is applied to simulate a weak millisecond pulsar magnetosphere with realistic scales using high-resolution axisymmetric simulations. Magnetospheric features reported by previous kinetic models are recovered, and strong synchrotron radiation above 100MeV consistent with the Fermi-LAT gamma-ray pulsar population is successfully reproduced. Conclusions.This work further consolidates the shining-reconnecting current sheet scenario as the origin of the gamma-ray emission in pulsars, as well as firmly establishing pulsar magnetospheres as at least teraelectronvolt particle accelerators. [ABSTRACT FROM AUTHOR]

Context. Galaxy cluster mergers are believed to generate large-scale shock waves that are ideal sites for cosmic ray production. In these so-called radio relic shocks, synchrotron radiation is produced mainly as a result of electron acceleration in the presence of intracluster magnetic fields. Aims. We aim to compute radio emission light curves for a sample of galaxy group and cluster mergers simulated in a cosmological context in order to study the dependence of radio luminosity on cluster mass, redshift, and impact parameter. Methods. We used model galaxy clusters from THE THREE HUNDRED project, a sample of 324 simulated high-density regions of radius 15 h−1 Mpc extracted from a cosmological volume, to identify cluster mergers characterised by the two main merging structures, construct their light curves, and follow their evolution throughout the complete simulated cosmic history. Results. We found that the median non-thermal radio relic luminosity light curve produced in galaxy cluster mergers can be described by a skewed Gaussian function abruptly rising after core-passage of the secondary cluster that peaks after ~0.1–0.8 Gyr as a function of M200,1, the mass of the primary, displaying a mass-dependent luminosity output increase of ≲10 to about ≳10–50 times relative to the radio emission measured at core-passage for galaxy groups and clusters, respectively. In general, most merger orbits are fairly radial with a median opening angle of ~20º before the collision. We also found that, independent of the cluster mass, less radial mergers tend to last longer, although the trend is weak. Finally, in agreement with previous works, we found that the peak radio luminosity shows a significant correlation with mass, P1.4 ∝ M200,12.05 demonstrating that this relation holds all the way up from galaxy group scales to the most massive galaxy clusters. Conclusions. We conclude that cluster mass is the primary driver for radio 'gischt' median luminosity, although there are significant variations for a given cluster mass related to the specifics of each merging process. In general, binary mergers are responsible for many of the well-known observed radio relic structures but complex situations involving three or more substructures are also common. Our simulations suggest that the shock-driven, non-thermal radio emission observed on cluster outskirts are the result of massive galaxy cluster mergers at ɀ ≲ 1, peaking at ɀ ~ 0–0.5. [ABSTRACT FROM AUTHOR]

Aims. G321.3–3.9 was first identified as a partial shell at radio frequencies a few decades ago. Although it continued to be observed, no additional studies were undertaken until recently. Methods. In this paper, we present results from a large selection of radio and X-ray data that cover the position of G321.3–3.9. We confirmed G321.3–3.9 as a new supernova remnant (SNR) using data collected by several radio surveys, spanning a frequency range from 200 to 2300 MHz. Stacked eROSITA data from four consecutive all-sky surveys (eRASS:4) provide spectro-imaging information in the energy band 0.2–8.0 keV. Results. G321.3–3.9 has an elliptical shape with major and minor axes of approximately 1.°7 × 1.°1. From CHIPASS and S-PASS data, we calculate a spectral index α = −0.8 ± 0.2, consistent with synchrotron emission from an expanding shell in the radiative phase. The eROSITA data show an X-ray diffuse structure filling almost the entire radio shell. Based on our spectral analysis, we found the temperature to be approximately 0.6 keV and the column absorption density about 1021 cm−2. Comparing this absorption density to optical extinction maps, we estimated the distance to fall within the range of (1.0–1.7) kpc, considering the 1σ uncertainty range. [ABSTRACT FROM AUTHOR]

Context. The electron-cyclotron maser instability (ECMI) is a significant coherent radio emission mechanism widely utilized in various astrophysical radio phenomena. It is well known that the velocity anisotropic distribution of energetic electrons, which leads to an inverted perpendicular population in the vertical direction with ∂fb/∂v⊥ > 0, can provide the free energy necessary for the ECMI. Aims. The initial velocity distribution of energetic electrons leaving the acceleration region is typically isotropic or beam-like. However, as these energetic electrons travel along the magnetic field as fast electron beams (FEBs) in magnetic plasma, various velocity anisotropic distributions can emerge. In this paper, we examine the impact of temperature anisotropy formed by beam electrons traveling along a flare loop on the ECMI. Methods. By neglecting the energy loss of energetic electrons as they traverse the corona and invoking the conservation of energy and magnetic moments, we established the relationship between momentum dispersion and the magnetic field. Utilizing the magnetic field model of the flare loop, we calculated the evolution of momentum dispersion and the growth rates of the ECMI as FEBs precipitate along the flare loop. Results. The results demonstrate that the temperature anisotropy arising as FEBs descend along the flare loop significantly impacts the ECMI. The maximum growth rates of the excited modes exhibit a gradual increase initially and then decline rapidly after reaching a critical height for β0 = 0.2c and 0.15c. The results also show that the growth rates of the O2 mode are one order of magnitude smaller than those of the O1 and X2 modes. This indicates that the harmonic radiation is X-mode polarized. Notably, the temperature anisotropy of FEBs as they precipitate along the flare loop with different magnetic field models or at different heights has similar effects on the ECMI. [ABSTRACT FROM AUTHOR]

We present a complete numerical model of the afterglow of a laterally structured relativistic ejecta from the radio to very high energies (VHE). This includes a self-consistent calculation of the synchrotron radiation, with its maximum frequency, and of synchrotron self-Compton (SSC) scattering that takes the Klein-Nishina regime into account. Attenuation due to pair production is also included. This model is computationally efficient and allows multi-wavelength data fitting. As a validation test, the radiative model was used to fit the broad-band spectrum of GRB 190114C at 90 s up to the TeV range. The full model was then used to fit the afterglow of GW 170817 and predict its VHE emission. We find that the SSC flux at the peak was much dimmer than the upper limit from H.E.S.S. observations. However, we show that either a smaller viewing angle or a higher external density would make similar off-axis events detectable in the future at VHE, even above 100 Mpc with the sensitivity of the Cherenkov telescope array. High external densities are expected in the case of fast mergers, but the existence of a formation channel for these binary neutron stars is still uncertain. We highlight that VHE afterglow detections would help to efficiently probe systems like this. [ABSTRACT FROM AUTHOR]

Motivated by the recent detection of ultralong-period radio transients, we investigate new models of coherent radio emission via low-altitude electron–positron pair production in neutron stars (NSs) beyond rotationally powered curvature radiation deathlines. We find that plastic motion (akin to 'continental drift') and qualitatively similar thermoelectric action by temperature gradients in the crusts of slowly rotating, highly magnetized NSs could impart mild local magnetospheric twists. Regardless of which mechanism drives twists, we find that particle acceleration initiates pair cascades across charge-starved gaps above a mild critical twist. Cascades are initiated via resonant inverse-Compton scattered photons or curvature radiation, and may produce broad-band coherent radio emission. We compute the pair luminosity (maximum allowed radio luminosity) for these two channels, and derive deathlines and 'active zones' in |$P-\dot{P}$| space from a variety of considerations. We find these twist-initiated pair cascades only occur for magnetar-like field strengths |$B \gtrsim 10^{14}$|  G and long periods: |$P_{\rm RICS} \gtrsim 120 \,\, (T/10^{6.5} {\rm K})^{-5} \, {\rm s}$| and |$P_{\rm curv} \gtrsim 150 \,\, ({\rm v_{\rm pl}}/10^{3} {\, \rm cm \, yr^{-1}})^{-7/6} \, {\rm s}$|⁠. Using a simplified geometric model, we find that plastic motion or thermoelectrically driven twists might naturally reproduce the observed luminosities, time-scales, and timing signatures. We further derive 'active zones' in which rotationally powered pair creation occurs via resonantly scattered photons, beyond standard curvature deathlines for pulsars. All cascades are generically accompanied by simultaneous (non-)thermal X-ray/UV counterparts which might be detectable with current instrumentation. [ABSTRACT FROM AUTHOR]

At optical/ultraviolet energies, blazars display an underlying thermal (unpolarized) contribution from the accretion disc, torus, and line-emitting regions, diluting the polarized emission from the jet-component. Optical polarimetry can be used to disentangle the thermal and non-thermal components, and place constraints on the particle populations and acceleration mechanisms responsible for the non-thermal emission. We present the results of a linear optical spectropolarimetric observing campaign of 18 blazars (six BL Lac-type objects and 12 flat-spectrum radio quasars) undertaken with the Southern African Large Telescope between 2016 and 2022. This was done to observe these systems during flaring states, as well as long-term monitoring of PKS1510-089, AP Lib, and PKS 1034-293. The observations traced the frequency dependence of the degree and angle of polarization, as well as changes in the spectral line strengths. We investigated possible correlations between the polarization and other observed characteristics for the sources. While an indication of correlation was found between the frequency dependence and the average level of polarization for some sources, a correlation was not found for the population as a whole. These results highlight that continuous observations and in-depth modelling of polarization and its frequency dependence is required to obtain a more holistic view of TeV blazars. [ABSTRACT FROM AUTHOR]

The recent detection of gamma-ray burst GRB 221009A has attracted attention due to its record brightness and first-ever detection of |$\gtrsim 10$| TeV gamma-rays from a GRB. Despite being the second-nearest GRB ever detected, at a redshift of |$z=0.151$|⁠ , the distance is large enough for severe attenuation of gamma-ray flux at these energies due to |$\gamma \gamma \rightarrow e^\pm$| pair production with the extragalactic background light (EBL). Here, we investigate whether the presence of cosmic voids along the line of sight can significantly impact the detectability of very high energy (VHE, |$\gt $| 100 GeV) gamma-rays from distant sources. Notably, we find that the gamma–gamma opacity for VHE gamma-rays can be reduced by approximately 10 per cent and up to 30 per cent at around 13 TeV, the highest-energy photon detected from GRB 221009A, for intervening cosmic voids along the line of sight with a combined radius of 110 Mpc, typically found from void catalogues, and 250 Mpc, respectively. This reduction is substantially higher for TeV photons compared to GeV photons, attributable to the broader target photon spectrum that TeV photons interact with. This finding implies that VHE photons are more susceptible to variations in the EBL spectrum, especially in regions dominated by cosmic voids. Our study sheds light on the detection of |$\gtrsim 10$| TeV gamma-rays from GRB 221009A in particular, and on the detection of extragalactic VHE sources in general. [ABSTRACT FROM AUTHOR]

Supermassive black hole binaries (SMBHBs) are natural by-products of galaxy mergers and are expected to be powerful multimessenger sources. They can be powered by the accretion of matter and then radiate across the electromagnetic spectrum, much like normal active galactic nuclei (AGNs). Current electromagnetic observatories have a good chance of detecting and identifying these systems in the near future. However, precise observational indicators are needed to distinguish individual AGNs from SMBHBs. In this paper, we propose a novel electromagnetic signature from SMBHBs: non-thermal emission produced by the interaction between the jets ejected by the black holes. We study close SMBHBs, which accrete matter from a circumbinary disc and the mini-discs formed around each hole. Each black hole ejects a magnetically dominated jet in the direction of its spin through the Blandford–Znajek mechanism. We argue that in such a situation, the interaction between the jets can trigger strong magnetic reconnection events, where particles are accelerated and emit non-thermal radiation. Depending on whether the jets are aligned or misaligned, this radiation can have different periodicities. We model the evolution of the particles accelerated during the dual jet interaction and calculate their radiative output, obtaining spectra and providing estimates for the variability time-scales. We finally discuss how this emission compares with that of normal AGNs. [ABSTRACT FROM AUTHOR]

Context. Pulsars are detected over the whole electromagnetic spectrum, from radio wavelengths up to very high energies, in the GeV-TeV range. While the radio emission site for young pulsars is well constrained to occur at altitudes about several percent of the light-cylinder radius and γ-ray emission is believed to be produced in the striped wind, outside the light cylinder, the non-thermal X-ray production site remains unknown. Aims. The aim of this Letter is to localize the non-thermal X-ray emission region based on a multi-wavelength pulse profile fitting for PSR J2229+6114, which stands as a particularly good candidate given its high X-ray brightness. Methods. Based on the geometry deduced from the joint radio and γ-ray pulse profiles, we fixed the magnetic axis inclination angle and the line-of-sight (LoS) inclination angle. However, we left the region of X-ray emission unlocalized, setting it somewhere between the surface and the light cylinder. We localized this region and its extension by fitting the X-ray pulse profile as observed by the NICER, NuSTAR, and RXTE telescopes in the ranges of 2–7 keV, 3–10 keV, and 9.4–22.4 keV, respectively. Results. We constrained the non-thermal X-ray emission to arise from altitudes between 0.2 rL and 0.55 rL where rL is the light-cylinder radius. The magnetic obliquity is approximately α ≈ 45° −50° and the LoS inclination angle is ζ ≈ 32° −48°. Conclusions. This Letter is among the first works to offer a tight constraint on the location of non-thermal X-ray emission from pulsars. We plan to apply this procedure to several other promising candidates to confirm this new result. [ABSTRACT FROM AUTHOR]

Context. Extreme TeV BL Lacs are a class of blazars with unique spectral and temporal features that are not easily reproducible using standard one-zone models based on single shock acceleration. To account for their peculiar properties, we elaborated a two-step acceleration model in which a recollimation shock and the subsequent downstream turbulence energize non-thermal electrons. Aims. We applied the model to a sample of extreme TeV BL Lacs with well-characterized spectral energy distributions. Since we used several sources, we automatized the exploration of the parameter space. This allowed us to derive the parameter distributions and study the correlations among them. Methods. We numerically solved a system of two coupled nonlinear differential equations to obtain the non-thermal particles and turbulence spectra. We calculated the spectral energy distribution via the synchrotron self-Compton emission model. The automatization of the parameter space exploration is possible through a Markov chain Monte Carlo (MCMC) ensemble sampler, in our case emcee. Results. We derived well-defined posterior distributions for the parameters, showing that the model is well constrained by available data and demonstrating the suitability of our method. The cross-correlations among some of the physical parameters are not trivial. Therefore, we conclude that MCMC sampling is a key instrument for characterizing the complexity of our multiparameter phenomenological model. [ABSTRACT FROM AUTHOR]

Context. The recent 230 GHz observations from the Event Horizon Telescope (EHT) collaboration are able to image the innermost structure of the M87 galaxy showing the shadow of the black hole, a photon ring, and a ring-like structure that agrees with thermal synchrotron emission from the accretion disc. However, at lower frequencies, M87 is characterized by a large-scale jet with clear signatures of nonthermal emission. It is necessary to explore the impacts of nonthermal emission on black hole shadow images and extended jets, especially at lower frequencies. Aims. In this study, we aim to compare models with different electron heating prescriptions to one another and to investigate how these prescriptions and nonthermal electron distributions may affect black hole shadow images and the broadband spectrum energy distribution (SED) function. Methods. We performed general relativistic radiative transfer (GRRT) calculations in various two-temperature general relativistic magnetohydrodynamic (GRMHD) models utilizing different black hole spins and different electron heating prescriptions coupled with different electron distribution functions (eDFs). Results. Through a comparison with GRRT images and SEDs, we found that when considering a variable κ eDF, the parameterized prescription of the R − β electron temperature model with Rh = 1 is similar to the model with electron heating in the morphology of images, and the SEDs at a high frequency. This is consistent with previous studies using the thermal eDF. However, the nuance between them could be differentiated through the diffuse extended structure seen in GRRT images, especially at a lower frequency, and the behavior of SEDs at low frequency. The emission from the nearside jet region is enhanced in the case of electron heating provided by magnetic reconnection and it will increase if the contribution from the regions with stronger magnetization is included or if the magnetic energy contribution to κ eDF mainly in the magnetized regions is considered. Compared with the thermal eDF, the peaks of the SEDs shift to a lower frequency when we consider nonthermal eDF. [ABSTRACT FROM AUTHOR]

PSR B1259–63/LS 2883 is a gamma-ray binary system that hosts a pulsar in an eccentric orbit, with a 3.4 yr period, around an O9.5Ve star (LS 2883). At orbital phases close to periastron passages, the system radiates bright and variable non-thermal emission, for which the temporal and spectral properties of this emission are, for now, poorly understood. In this regard, very high-energy (VHE) emission is especially useful to study and constrain radiation processes and particle acceleration in the system. We report on an extensive VHE observation campaign conducted with the High Energy Stereoscopic System, comprised of approximately 100 h of data taken over five months, from tp − 24 days to tp + 127 days around the system's 2021 periastron passage (where tp is the time of periastron). We also present the timing and spectral analyses of the source. The VHE light curve in 2021 is consistent overall with the stacked light curve of all previous observations. Within the light curve, we report a VHE maximum at times coincident with the third X-ray peak first detected in the 2021 X-ray light curve. In the light curve – although sparsely sampled in this time period – we see no VHE enhancement during the second disc crossing. In addition, we see no correspondence to the 2021 GeV flare in the VHE light curve. The VHE spectrum obtained from the analysis of the 2021 dataset is best described by a power law of spectral index Γ = 2.65 ± 0.04stat ± 0.04sys, a value consistent with the spectral index obtained from the analysis of data collected with H.E.S.S. during the previous observations of the source. We report spectral variability with a difference of ΔΓ = 0.56 ± 0.18stat ± 0.10sys at 95% confidence intervals, between sub-periods of the 2021 dataset. We also detail our investigation into X-ray/TeV and GeV/TeV flux correlations in the 2021 periastron passage. We find a linear correlation between contemporaneous flux values of X-ray and TeV datasets, detected mainly after tp + 25 days, suggesting a change in the available energy for non-thermal radiation processes. We detect no significant correlation between GeV and TeV flux points, within the uncertainties of the measurements, from ∼tp − 23 days to ∼tp + 126 days. This suggests that the GeV and TeV emission originate from different electron populations. [ABSTRACT FROM AUTHOR]

We report on new multifrequency Very Large Array (VLA) radio observations and Chandra X-ray observations of a radio-loud quasar with a ∼300 kpc-long jet, PKS 1127−145, during a flaring event detected in γ-rays by the Fermi Large Area Telescope in December 2020. The high angular resolution of the new radio images allows us to disentangle for the first time the kiloparsec-scale inner jet from the core contribution. The inner radio jet, up to 15 kpc from the core, is highly polarized (33 percent) and the magnetic field is parallel to the jet axis. At about 18 arcsecs from the core, the jet slightly bends and we observe a re-brightening of the radio emission and a 90-degree rotation of the magnetic field, likely highlighting the presence of a shock that is compressing the magnetic field to a plane perpendicular to the jet axis, where efficient particle acceleration takes place. At the same position, the X-ray emission fades, suggesting a deceleration of the bulk velocity of the jet after the bend. A change in velocity and collimation of the jet is supported by the widening of the jet profile and the detection of a limb-brightened structure connecting the bending region with the jet termination. The limb-brightened structure might indicate the coexistence of both longitudinal and transverse velocity gradients at the jet bending. There is no evidence of significant brightening of the kiloparsec-scale jet in the radio or X-ray band during the γ-ray flare. The X-ray flux, F2 − 10 keV = (6.24 ± 0.57)×10−12 ergs s−1 cm−2, measured by Chandra from the quasar core is consistent with the flux measured by the X-ray Telescope on board the Neil Gehrels Swift Observatory after the high-energy flare. Our results indicate that the γ-ray flaring region is located within the VLA source core. [ABSTRACT FROM AUTHOR]

Context. In recent years, some polarization measurements of gamma-ray bursts (GRBs) have been reported, and the polarization-angle (PA) rotation in the prompt emission phase has been found in several bursts. The physical mechanism of the PA evolution is still unclear. In this work, we studied the origin of the PA rotation in a toroidal magnetic field. Aims. We aim to provide an explanation for the PA rotation in GRBs and find the physical conditions that lead to the rotation by 90 degrees in the toroidal magnetic-field (MF) model. Moreover, we present some observable polarization properties in the MF model that can be tested in the future. Methods. We calculated the instantaneous polarization degree (PD) from a top-hat jet with different normalized viewing angles (q = θv/θj), jet opening angles (θj), and jet Lorentz factors (Γ) in three wavebands. When the PD changes between positive and negative values, it means that the PA flips by 90 degrees. On these grounds, we can summarize the range of parameters required for these PA flips. Considering these parameter conditions, we can further estimate the observed rate of the GRBs exhibiting such PA rotations. Results. We find that the PA rotation in the toroidal MF is primarily related to three critical factors: the viewing angle, the jet opening angle, and the jet Lorentz factor. Additionally, the PA can experience flips of 90 degrees twice. The conditions for the flips are q ≳ 0.5 (except for q ≃ 1) and yj = (Γθj)2 ≳ 4. However, the two flips in the PA might not be concurrently observable due to the constraint of flux. Taking these conditions into account and assuming a random orientation between the jet axis and the line of sight (LOS), we obtain a theoretical upper limit (without any constraints) for the observed rate of GRBs in the X-ray or γ-ray band displaying the flips in PA as Rch ≲ 80%. We further constrain the observed rate as Rch ∼ 16% according to the maximal post-flip polarized flux level, where the observed rate of single and double flips each account for ∼8%. It should be noted that the observed rates are different in various wavebands. The observed rate of the second PA flip in the optical bands should be higher than that in the X-ray or γ-ray band since the flux in the optical band declines much slower than that in the X-ray or γ-ray band. Moreover, when the LOS is close to the jet edge (q → 1), it is the easiest case in which to observe the 90-degree PA flip due to the relatively high post-flip polarized flux level. The first and second PA flips in a GRB pulse are most likely to occur at the observed times of tobs ∼ [2 − 3]tpeak and ∼[3 − 4]tpeak, respectively, where tpeak is the peak time of the pulse. It is also noted that the PA flip would not happen before the peak time. [ABSTRACT FROM AUTHOR]

Context. The last two decades have witnessed dramatic progress in our understanding of neutron star magnetospheres thanks to force-free and particle-in-cell simulations. However, the associated particle dynamics and its emission mechanisms and locations have not been fully constrained, notably in X-rays. Aims. In this paper, we compute a full atlas of radio, X-ray, and γ-ray pulse profiles, relying on the force-free magnetosphere model. Our goal is to use such a data bank of multi-wavelength profiles to fit a substantial number of radio-loud γ-ray pulsars that have also been detected in non-thermal X-rays to decipher the X-ray radiation mechanism and sites. Using results from the third γ-ray pulsar catalogue (3PC), we investigate the statistical properties of this population. Methods. We assume that radio emission emanates from field lines rooted to the polar caps, at varying height above the surface, close to the surface, at an altitude about 5–10% of the light cylinder radius, rL. The X-ray photons are produced in the separatrix region within the magnetosphere; that is, the current sheet formed by the jump from closed to open magnetic field lines. We allow for substantial variations in emission height. The γ-rays are produced within the current sheet of the striped wind, outside the light cylinder. Results. A comprehensive set of radio, X-ray, and γ-ray light curves was computed. Based on only geometric considerations about magnetic obliquity, line-of-sight inclination, and the radio beam cone opening angle, pulsars can be classified as radio-loud or quiet and as γ-ray-loud or quiet. We found that the 3PC sample is compatible with an isotropic distribution of obliquity and line of sight. Conclusions. The atlases constructed in this work are the fundamental tools with which to explore individual pulsars and fit their multi-wavelength pulse profiles in order to constrain their magnetic topology, the emission sites, and the observer's line of sight. [ABSTRACT FROM AUTHOR]

Context. Blazar flares provide a window onto the extreme physical processes occurring in relativistic outflows. Most numerical codes used for modelling blazar emission during flares use a simplified continuous-loss description of particle cooling due to the inverse Compton (IC) process, neglecting non-continuous (discrete) effects that arise in the Klein–Nishina (KN) regime. The significance of such effects has not yet been explored in detail. Aims. In this study, we investigate the importance of non-continuous Compton cooling losses and their impact on the electron spectrum and spectral energy distribution (SED) of blazars during high flux states (flares), as well as in the low state. Methods. We solve the full transport equation numerically, accounting for large relative jumps in energy by extending our existing blazar flare modelling code EMBLEM. We perform a detailed physical modelling of the brightest γ-ray flare of the archetypal flat-spectrum radio quasar (FSRQ) 3C 279 detected in June 2015. We then compare results obtained using the full cooling term and using the continuous-loss approximation. Results. We show that during flaring states of FSRQs characterised by high Compton dominance, the non-continuous cooling can lead to significant modification of the electron spectrum, introducing a range of distinct features, such as low-energy tails, hardening and/or softening, narrow and extended particle excesses, and shifts in the cooling break position. Such distortion translates to differences in the associated SED of up to ∼50%. This highlights the importance of non-continuous effects and the need to consider them in blazar emission models, particularly applied to extreme γ-ray flares. [ABSTRACT FROM AUTHOR]

Pulsations from the Crab pulsar have been detected by the MAGIC telescopes at energies up to 1.5 TeV, and the pulsed emission from the Vela pulsar was detected by H.E.S.S. reaching tens of TeV. These discoveries, along with the proposed additional emission due to inverse Compton scattering at TeV energies, lead us to consider suitable candidates for detection with current and future extensive air show (EAS) experiments at very-high-energy (VHE; 0.1 |$-$| 100 TeV) ranges. Leveraging energy spectrum data from pulsars as observed by Fermi and Imaging Atmospheric Cherenkov Telescopes (IACTs) and considering the sensitivities of both LHAASO and SWGO, this study evaluates their detectability and estimates the time required for their significant detection. Our results indicate that LHAASO could detect the Crab's pulsed signal within six years, while SWGO might detect Vela's signal within one year. Observations of the most energetic Fermi pulsars with EAS experiments will provide insight into the nature of VHE pulsar emissions, helping us to clarify the primary characteristics of VHE pulsars. [ABSTRACT FROM AUTHOR]

Magnetic reconnection is often considered as the primary particle acceleration mechanism in a magnetized blazar zone environment. The majority of radiation in the reconnection layer comes from plasmoids and their mergers. In particular, plasmoid mergers can produce strong multiwavelength flares and major variations in synchrotron polarization signatures. However, radiative properties of plasmoid mergers have not been well explored due to difficulties in tracking the merging processes. Here we use an image processing method that combines the magnetic vector potential and density to identify isolated and merging plasmoids. We find that this method can clearly distinguish radiation contributions from isolated plasmoids, merging plasmoids, and the primary current sheet of reconnection. This new method enables us to study the radiative properties of plasmoids and mergers statistically. Our results show that isolated plasmoids have similar emissivity regardless of their sizes, and they generally have non-zero polarization degree (PD) due to their quasi-circular shape. Flares due to plasmoid mergers have relative amplitudes that are antiproportional to the size ratio of the plasmoids participating in the mergers. Finally, only mergers between plasmoids of comparable sizes (width ratio ≲5) can lead to significant spectral hardening and polarization angle (PA) variations; the amplitude of the PA variations is between 0 and 180° and has a mean value of 90°. Our analyses on 2D simulations can pave the way for future analyses and machine learning techniques on radiative properties of 3D magnetic reconnection simulations. [ABSTRACT FROM AUTHOR]

The S-shaped swing of the linear polarization position angle (PPA) observed in many pulsars can be interpreted by the rotating vector model (RVM). However, efforts to fit the RVM for a large sample of pulsars observed with the MeerKAT telescope as a part of the Thousand-Pulsar-Array (TPA) programme, only succeeded for about half the cases. High time-resolution studies suggest that the failed cases arise due to the presence of orthogonal polarization modes, or highly disordered distribution of PPA points. One such example is PSR J1645−0317. Recently it has been shown that the RVM can be recovered in this pulsar by using only time samples which are greater than 80 per cent linearly polarized. In this work, we test this novel approach on the brightest 249 pulsars from the TPA sample, of which 177 yield sufficient highly polarized samples to be amenable to our method. Remarkably, only nine of these pulsars (5 per cent) now fail to fit the RVM as opposed to 59 per cent from the original analysis. This result favours the paradigm that the underlying mechanism is coherent curvature radiation. [ABSTRACT FROM AUTHOR]

Models for gamma-ray burst afterglow dynamics and synchrotron spectra are known to exhibit various scale invariances, owing to the scale-free nature of fluid dynamics and the power-law shape of synchrotron spectra. Since GRB 170817A, off-axis jet models including a lateral energy structure in the initial outflow geometry have gained in prominence. Here, we demonstrate how the scale invariance for arbitrary jet structure and dynamical stage can be expressed locally as a function of jet temporal light-curve slope. We provide afterglow flux expressions and demonstrate their use to quickly assess the physical implications of observations. We apply the scaling expressions to the Swift X-ray Telescope sample, which shows a spread in observed fluxes, binned by light-curve slope at time of observation, that increases with increasing light-curve slope. According to the scaling relations, this pattern is inconsistent with a large spread in environment densities if these were the dominant factor determining the variability of light curves. We further show how the late deep Newtonian afterglow stage remains scale-invariant but adds distinct spectral scaling regimes. Finally, we show that for given jet structure a universal curve can be constructed of the centroid offset, image size, and ellipticity (that can be measured using very large baseline interferometry) versus observer angle, in a manner independent of explosion energy and circumburst density. Our results apply to any synchrotron transient characterized by a release of energy in an external medium, including supernova remnants, kilonova afterglows, and soft gamma-repeater flares. [ABSTRACT FROM AUTHOR]

The origin of the radio synchrotron background (RSB) is currently unknown. Its understanding might have profound implications in fundamental physics or might reveal a new class of radio emitters. In this work, we consider the scenario in which the RSB is due to extragalactic radio sources and measure the angular cross-correlation of Low-Frequency Array (LOFAR) images of the diffuse radio sky with matter tracers at different redshifts, provided by galaxy catalogues and cosmic microwave background lensing. We compare these measured cross-correlations to those expected for models of RSB sources. We find that low-redshift populations of discrete sources are excluded by the data, while higher redshift explanations are compatible with available observations. We also conclude that at least 20 per cent of the RSB surface brightness level must originate from populations tracing the large-scale distribution of matter in the Universe, indicating that at least this fraction of the RSB is of extragalactic origin. Future measurements of the correlation between the RSB and tracers of high-redshift sources will be crucial to constraining the source population of the RSB. [ABSTRACT FROM AUTHOR]

A substantial fraction of quasars display broad absorption lines (BALs) in their rest-frame ultraviolet spectra. While the origin of BALs is thought to be related to the accretion disc wind, it remains unclear whether the observed ratio of BAL to non-BAL quasars is a result of orientation. We conducted observations of 48 BAL quasars and the same number of non-BAL quasars at 322 MHz using the Giant Metrewave Radio Telescope. Combined with previous flux measurements ranging from MHz to GHz frequencies, we compared continuum radio spectra between the two quasar groups. These data offer insights into low-frequency radio properties that have been difficult to investigate with previous observations only at GHz frequencies. Our results show that 73 ± 13 per cent of the BAL quasars exhibit steep or peaked spectra, a higher proportion than the 44 ± 14 per cent observed in the non-BAL quasars. In contrast, there are no discernible differences between the two quasar groups in the radio luminosity, peak frequency, and spectral index distributions of sources with steep or peaked spectra and sources with flat or inverted spectra. Generally, as the jet axis and line of sight become closer to parallel, quasars exhibit flat or inverted spectra rather than steep or peaked spectra. Therefore, these results suggest that BAL quasars are more frequently observed farther from the jet axis than non-BAL quasars. However, given that a certain proportion of BAL quasars exhibit flat or inverted spectra, more than the simple orientation scenario is required to elucidate the radio properties of BAL quasars. [ABSTRACT FROM AUTHOR]

This work presents a methodological approach to generate realistic γ-ray light curves of pulsars, resembling reasonably well the observational ones observed by the Fermi -Large Area Telescope (Fermi -LAT) instrument, fitting at the same time their high-energy spectra. The theoretical light curves are obtained from a spectral and geometrical model of the synchro-curvature emission. Despite our model relying on a few effective physical parameters, the synthetic light curves present the same main features observed in the observational γ-ray light-curve zoo, such as the different shapes, variety in the number of peaks, and a diversity of peak widths. The morphological features of the light curves allow us to statistically compare the observed properties. In particular, we find that the proportion of the number of peaks found in our synthetic light curves is in agreement with the observational one provided by the third Fermi -LAT pulsar catalogue. We also found that the detection probability due to beaming is much higher for orthogonal rotators (approaching 100  per cent) than for small inclination angles (less than 20 per cent). The small variation in the synthetic skymaps generated for different pulsars indicates that the geometry dominates over timing and spectral properties in shaping the γ-ray light curves. This means that geometrical parameters such as the inclination angle can be in principle constrained by γ-ray data alone independently of the specific properties of a pulsar. At the same time, we find that γ-ray spectra seen by different observers can slightly differ, opening the door to constraining the viewing angle of a particular pulsar. [ABSTRACT FROM AUTHOR]

The study of the high-energy part (MeV-GeV) of the spectrum of gamma-ray bursts (GRBs) can play a crucial role in investigating the physics of prompt emission, but it is often hampered by low statistics and the paucity of GeV observations. In this work, we analyze the prompt emission spectra of the 22 brightest GRBs which have been simultaneously observed by Fermi/GBM and Fermi/LAT, spanning six orders of magnitude in energy. The high-energy photon spectra can be modeled with a power-law N(E)∝E−β possibly featuring an exponential cutoff. We find that, with the inclusion of the LAT data, the spectral index β is softer than what is typically inferred from the analysis of Fermi/GBM data alone. Under the assumption that the emission is synchrotron, we derived a median value of the index p ∼ 2.79 of the power-law energy distribution of accelerated particles (N(γ)∝γ−p). In nine out of 22 GRB spectra, we find a significant presence of an exponential cutoff at high energy, ranging between 14 and 298 MeV. By interpreting the observed cutoff as a sign of pair-production opacity, we estimate the jet bulk Lorentz factor Γ, finding values in the range 130–330. These values are consistent with those inferred from the afterglow light curve onset time. Finally, by combining the information from the high-energy prompt emission spectrum with the afterglow light curve, we exploited a promising method to derive the distance R from the central engine where the prompt emission occurs. The distances (R > 1013 − 15 cm) inferred for the only two GRBs in our sample that are suitable for the application of this method, which have only lower limits on their cutoff energies, suggest large emitting regions, although they are still compatible with the standard model. Larger samples of GRBs with measured cutoff energies and afterglow deceleration time will allow for more informative values to be derived. These results highlight the importance of including high-energy data, when available, in the study of prompt spectra and their role in addressing the current challenges of the GRB standard model. [ABSTRACT FROM AUTHOR]

Aims. We present the first multiwavelength study of Mrk 501 that contains simultaneous very-high-energy (VHE) γ-ray observations and X-ray polarization measurements from the Imaging X-ray Polarimetry Explorer (IXPE). Methods. We used radio-to-VHE data from a multiwavelength campaign carried out between March 1, 2022, and July 19, 2022 (MJD 59639 to MJD 59779). The observations were performed by MAGIC, Fermi-LAT, NuSTAR, Swift (XRT and UVOT), and several other instruments that cover the optical and radio bands to complement the IXPE pointings. We characterized the dynamics of the broadband emission around the X-ray polarization measurements through its multiband fractional variability and correlations, and compared changes observed in the polarization degree to changes seen in the broadband emission using a multi-zone leptonic scenario. Results. During the IXPE pointings, the VHE state is close to the average behavior, with a 0.2–1 TeV flux of 20%–50% of the emission of the Crab Nebula. Additionally, it shows low variability and a hint of correlation between VHE γ-rays and X-rays. Despite the average VHE activity, an extreme X-ray behavior is measured for the first two IXPE pointings, taken in March 2022 (MJD 59646 to 59648 and MJD 59665 to 59667), with a synchrotron peak frequency > 1 keV. For the third IXPE pointing, in July 2022 (MJD 59769 to 59772), the synchrotron peak shifts toward lower energies and the optical/X-ray polarization degrees drop. All three IXPE epochs show an atypically low Compton dominance in the γ-rays. The X-ray polarization is systematically higher than at lower energies, suggesting an energy stratification of the jet. While during the IXPE epochs the polarization angles in the X-ray, optical, and radio bands align well, we find a clear discrepancy in the optical and radio polarization angles in the middle of the campaign. Such results further support the hypothesis of an energy-stratified jet. We modeled broadband spectra taken simultaneous to the IXPE pointings, assuming a compact zone that dominates in the X-rays and the VHE band, and an extended zone stretching farther downstream in the jet that dominates the emission at lower energies. NuSTAR data allow us to precisely constrain the synchrotron peak and therefore the underlying electron distribution. The change between the different states observed in the three IXPE pointings can be explained by a change in the magnetization and/or the emission region size, which directly connects the shift in the synchrotron peak to lower energies with the drop in the polarization degree. [ABSTRACT FROM AUTHOR]

Fast radio bursts (FRBs) can present a variety of polarization properties and some of them are characterized by narrow spectra. In this work, we study spectral properties from the perspective of intrinsic radiation mechanisms and absorption through the waves propagating in the magnetosphere. The intrinsic radiation mechanisms are considered by invoking quasi-periodic bunch distribution and perturbations on charged bunches moving on curved trajectories. The narrowband emission is likely to reflect some quasi-periodic structure on the bulk of bunches, which may be due to quasi-periodically sparking in a "gap" or quasi-monochromatic Langmuir waves. A sharp spike would appear in the spectrum if the perturbations were to induce a monochromatic oscillation of bunches; however, it is difficult to create a narrow spectrum because the Lorentz factor has large fluctuations, so the spike disappears. Both the bunching mechanism and perturbations scenarios share the same polarization properties, with a uniformly distributed bulk of bunches. We investigated the absorption effects, including Landau damping and curvature self-absorption in the magnetosphere, which are significant at low frequencies. Subluminous O-mode photons cannot escape from the magnetosphere due to the Landau damping, leading to a height-dependent lower frequency cut-off. The spectra can be narrow when the frequency cut-off is close to the characteristic frequency of curvature radiation, however, such conditions cannot always be met. The spectral index is 5/3 at low-frequency bands due to the curvature self-absorption is not as steep as what is seen in observations. The intrinsic radiation mechanisms are more likely to generate the observed narrow spectra of FRBs, rather than the absorption effects. [ABSTRACT FROM AUTHOR]

Recent observations have revealed detailed structures of radio relics across a wide range of frequencies. In this work, we performed three-dimensional magnetohydrodynamical (3D MHD) simulations of merger shocks propagating through a turbulent magnetized intracluster medium. We then employed on-the-fly Lagrangian particles to explore the physical processes behind the origination of radio substructures and their appearance in high and low-frequency observations. We employed two cosmic-ray (CR) electron acceleration models, with a fresh injection of electrons from the thermal pool and the re-acceleration of mildly relativistic electrons. We used the relative surface brightness fluctuations, δSν, to define a "degree of patchiness." First, we found that patchiness is produced if the shock's surface has a distribution of Mach numbers, rather than a single Mach number. Second, radio relics appear patchier if the Mach number distribution consists of a large percentage of low Mach numbers (ℳ ≲ 2.5). Furthermore, as the frequency increases, the patchiness also becomes larger. Nevertheless, if radio relics are patchy at high frequencies (e.g., 18.6 GHz), they necessarily will also be patchy at low frequencies (e.g., 150 MHz). Then, to produce noticeable differences in the patchiness at low and high frequencies, the shock front should have a Mach number spread of σℳ ≳ 0.3 − 0.4. Finally, the extent of the patchiness depends on the Mach number distribution as well as the CR acceleration model. We propose δSν as a potential tool for extracting merger shock properties and information on particle acceleration processes at shocks in radio observations. [ABSTRACT FROM AUTHOR]

Context. Very high-energy (VHE, E > 100 GeV) observations of the blazar Mrk 501 with MAGIC in 2014 provided evidence for an unusual narrow spectral feature at about 3 TeV during an extreme X-ray flaring activity. The one-zone synchrotron-self Compton scenario, widely used in blazar broadband spectral modeling, fails to explain the narrow TeV component. Aims. Motivated by this rare observation, we propose an alternative model for the production of narrow features in the VHE spectra of flaring blazars. These spectral features may result from the decay of neutral pions (π0 bumps) that are in turn produced via interactions of protons (of tens of TeV energy) with energetic photons, whose density increases during hard X-ray flares. Methods. We explored the conditions needed for the emergence of narrow π0 bumps in VHE blazar spectra during X-ray flares reaching synchrotron energies ∼100 keV using time-dependent radiative transfer calculations. We focused on high-synchrotron peaked (HSP) blazars, which comprise the majority of VHE-detected extragalactic sources. Results. We find that synchrotron-dominated flares with peak energies ≳100 keV can be ideal periods for the search of π0 bumps in the VHE spectra of HSP blazars. The flaring region is optically thin to photopion production, its energy content is dominated by the relativistic proton population, and the inferred jet power is highly super-Eddington. Application of the model to the spectral energy distribution of Mrk 501 on MJD 56857.98 shows that the VHE spectrum of the flare is described well by the sum of a synchrotron self-Compton (SSC) component and a distinct π0 bump centered at 3 TeV. Spectral fitting of simulated SSC+π0 spectra for the Cherenkov Telescope Array (CTA) show that a π0 bump could be detected at a 5σ significance level with a 30-min exposure. Conclusions. A harder VHE γ-ray spectrum than the usual SSC prediction or, more occasionally, a distinct narrow bump at VHE energies during hard X-ray flares, can be suggestive of a relativistic hadronic component in blazar jets that otherwise would remain hidden. The production of narrow features or spectral hardenings due to π0 decay in the VHE spectra of blazars is testable with the advent of CTA. [ABSTRACT FROM AUTHOR]

Context. Spiral galaxies M31 and M33 are among the γ-ray sources detected by the Fermi Large Area Telescope (LAT). Aims. We aim to model the broadband non-thermal emission of the central region of M31 (a LAT point source) and of the disk of M33 (a LAT extended source), as part of our continued survey of non-thermal properties of local galaxies that includes the Magellanic Clouds. Methods. We analysed the observed emission from the central region of M31 (R < 5.5 kpc) and the disk-sized emission from M33 (R ∼ 9 kpc). For each galaxy, we self-consistently modelled the broadband spectral energy distribution of the diffuse non-thermal emission based on published radio and γ-ray data. All relevant radiative processes involving relativistic and thermal electrons (synchrotron, Compton scattering, bremsstrahlung, and free–free emission and absorption), along with relativistic protons (π0 decay following interaction with thermal protons), were considered, using exact emissivity formulae. We also used the Fermi-LAT-validated γ-ray emissivities for pulsars. Results. Joint spectral analyses of the emission from the central region of M31 and the extended disk of M33 indicate that the radio emission is composed of both primary and secondary electron synchrotron and thermal bremsstrahlung, whereas the γ-ray emission may be explained as a combination of diffuse pionic, pulsar, and nuclear-BH-related emissions in M31 and plain diffuse pionic emission (with an average proton energy density of 0.5 eV cm−3) in M33. Conclusions. The observed γ-ray emission from M33 appears to be mainly hadronic. This situation is similar to other local galaxies, namely, the Magellanic Clouds. In contrast, we have found suggestions of a more complex situation in the central region of M31, whose emission could be an admixture of pulsar emission and hadronic emission, with the latter possibly originating from both the disk and the vicinity of the nuclear black hole. The alternative modelling of the spectra of M31 and M33 is motivated by the different hydrogen distribution in the two galaxies: The hydrogen deficiency in the central region of M31 partially unveils emissions from the nuclear BH and the pulsar population in the bulge and inner disk. If this were to be the case in M33 as well, these emissions would be outshined by diffuse pionic emission originating within the flat central-peak gas distribution in M33. [ABSTRACT FROM AUTHOR]

We investigate the nature of a Galactic Centre source, G0.17+0.15, lying along the northern extension of the radio arc near l ∼ 0.2°. G0.17+0.15 is an H  ii region located towards the eastern edge of the radio bubble, embedded within the highly polarized Galactic Centre eastern lobe where a number of radio filaments appear to cross through the H  ii region. We report the detection of hydrogen and helium recombination lines with a radial velocity exceeding 140 km s−1 based on Green Bank Telescope and Very Large Array observations. The morphology of G0.17+0.15, aided by kinematics, and spectral index characteristics, suggests the presence of an external pressure dragging and shredding the ionized gas. We argue that this ionized cloud is interacting with a bundle of radio filaments and is entrained by the ram pressure of the radio bubble, which itself is thought to be produced by cosmic ray driven outflows at the Galactic Centre. In this interpretation, the gas streamers on the western side of G0.17+0.15 are stripped, accelerated from 0 to |$\delta v\sim \, 35$| km s−1 over a time-scale roughly 8 × 104 yr, implying that ablating ram pressure is |$\sim 700\, \mathrm{eV\, cm^{-3}}$|⁠ , comparable to the |$\sim 10^3\, \mathrm{eV \, cm^{-3}}$| cosmic ray driven wind pressure in the Galactic Centre region. [ABSTRACT FROM AUTHOR]

In previous studies, the time evolution information is missed in deducing the time-integrated polarizations of gamma-ray burst prompt emission. Here, it is considered and the time-integrated polarizations is investigated through the summation of the time-resolved ones. The statistical properties of the distribution of the time-integrated polarization degree (Π) can be read from the q– Π curve, where q ≡ θ V /θ j . θ V and θ j are the observational and jet half-opening angles, respectively. Hence, only the q– Π curves are studied. In addition to a toroidal magnetic field in the radiation region, an aligned field is also discussed. We found the predicted time-integrated polarization degree is around |$(40\!-\!50)\,\, \mathrm{ per\,cent}$| for High-energy Polarimetry Detector (HPD) onboard POLAR-2 and is roughly |$(30\!-\!40)\,\, \mathrm{ per\,cent}$| for its Low-energy Polarimetry Detector (LPD). Therefore, Π value detected by the HPD will be larger than that of the LPD in statistics and the result of the former estimations will underestimate the value of Π in an ordered field. There are mainly two types of the q –Π curve profiles, corresponding to two ordered magnetic field configurations. [ABSTRACT FROM AUTHOR]

In their catalogue of pulsars' radio spectra, Swainston et al. distinguish between five different forms of these spectra: those that can be fitted with (i) a simple power law, (ii) a broken power law, (iii) a low-frequency turn-over, (iv) a high-frequency turn-over or (v) a double turn-over spectrum. Here, we choose two examples from each of these categories and fit them with the spectral distribution function of the caustics that are generated by the superluminally moving current sheet in the magnetosphere of a non-aligned neutron star. In contrast to the prevailing view that the curved features of pulsars' radio spectra arise from the absorption of the observed radiation in high-density environments, our results imply that these features are intrinsic to the emission mechanism. We find that all observed features of pulsar spectra (including those that are normally fitted with simple or broken power laws) can be described by a single spectral distribution function and regarded as manifestations of a single emission mechanism. From the results of an earlier analysis of the emission from a pulsar's current sheet and the values of the fit parameters for each spectrum, we also determine the physical characteristics of the central neutron star of each considered example and its magnetosphere. [ABSTRACT FROM AUTHOR]

PG 1553 + 113 is one of the few blazars with a convincing quasi-periodic emission in the gamma-ray band. The source is also a very high energy (VHE; >100 GeV) gamma-ray emitter. To better understand its properties and identify the underlying physical processes driving its variability, the MAGIC Collaboration initiated a multiyear, multiwavelength monitoring campaign in 2015 involving the OVRO 40-m and Medicina radio telescopes, REM, KVA, and the MAGIC telescopes, Swift and Fermi satellites, and the WEBT network. The analysis presented in this paper uses data until 2017 and focuses on the characterization of the variability. The gamma-ray data show a (hint of a) periodic signal compatible with literature, but the X-ray and VHE gamma-ray data do not show statistical evidence for a periodic signal. In other bands, the data are compatible with the gamma-ray period, but with a relatively high p -value. The complex connection between the low- and high-energy emission and the non-monochromatic modulation and changes in flux suggests that a simple one-zone model is unable to explain all the variability. Instead, a model including a periodic component along with multiple emission zones is required. [ABSTRACT FROM AUTHOR]

We consider the broad-band spectral energy distribution of the high-energy-peaked blazar Mkn 501 using Swift -XRT/UVOT, NuSTAR and Fermi -LAT observations taken between 2013 and 2022. The spectra were fitted with a one-zone leptonic model using synchrotron and synchrotron self-Compton emission from different particle energy distributions such as a broken power law, log-parabola, as well as distributions expected when the diffusion or the acceleration time-scale are energy-dependent. The jet power estimated for a broken power-law distribution was ∼1047(1044) erg s−1 for a minimum electron energy γmin ∼ 10(103). However, for electron energy distributions with intrinsic curvature (such as the log-parabola form), the jet power is significantly lower at a few times 1042 erg s−1 which is a few per cent of the Eddington luminosity of a 107 M⊙ black hole, suggesting that the jet may be powered by accretion processes. We discuss the implications of these results. [ABSTRACT FROM AUTHOR]

In this work, we provide a detailed analysis of the broad-band temporal and spectral properties of the blazar Ton 599 by using observations from the Fermi Large Area Telescope (LAT) and Swift X-Ray Telescope (XRT)/Ultraviolet–Optical Telescope (UVOT), during its brightest γ-ray flaring. The one-day bin γ-ray light curve exhibits multiple substructures with asymmetric and symmetric profiles. Notably, the γ-ray light curve shows a maximum flux of |$\rm 3.63 \times 10^{-6}\, photon\, cm^{-2}\, s^{-1}$| on MJD 59954.50, which is the highest flux ever observed from this source. The correlation between the γ-ray flux and γ-ray spectral indices suggests a moderate 'harder when brighter' trend. Taking the γ-ray light curve as the reference, a strong correlation is observed with X-ray, optical, and UV energies. Additionally, the γ-rays and optical/UV emission exhibit higher variability compared with X-rays. To understand the parameter variation during the active state of the source, we conducted a statistical broad-band spectral modelling of the source in 10 flux intervals of equal duration. A one-zone leptonic model involving synchrotron, synchrotron-self-Compton, and external Compton processes successfully reproduces the broad-band spectral energy distribution (SED) in each of these flux intervals. We observed that flux variation during the active state is associated mainly with variation in the magnetic field and particle spectral indices. [ABSTRACT FROM AUTHOR]

Blazars, a subclass of active galactic nuclei (AGN), are known to be bright γ-ray sources, frequently exhibiting active (flaring) periods. Blazar PKS 2155-304 is a high synchrotron-peaked BL Lac object located at redshift z  = 0.116. On 2006 July 28, an extremely remarkable outburst of VHE γ-ray emission from this blazar was reported by the H.E.S.S. experiment, with an average flux more than 10 times the low-state level. The variability time-scale of this extraordinary flare was as short as approximately 200 s. In order to guarantee the transparency of the emission region for TeV photons, the fast variability demands an extremely high Doppler factor δD > 50 of the jet within the classical one-zone model, leading to the so-called 'Doppler factor crisis'. Here, we demonstrate that the stochastic dissipation model, which is a multiblob scenario for blazars, can self-consistently explain the giant TeV flares of PKS 2155-304 and the low-state emission before and after the flares, in terms of both multiwavelength spectral and variability characteristics. The required Doppler factor in this model can be as low as 20, which is a reasonable and typical value for blazar jets. The obtained model parameters may shed some light on the physical properties of the relativistic jet. [ABSTRACT FROM AUTHOR]

Previous studies on the fitting of spectral energy distributions (SEDs) often apply the external-Compton process to interpret the high-energy peak of low-synchrotron-peaked (LSP) BL Lac objects (LBLs), despite the lack of strong broad emission lines observed for LBLs. In this work, we collect quasi-simultaneous multiwavelength data of 15 LBLs from the Fermi fourth LAT AGN catalogue (4LAC). We propose an analytical method to assess the necessity of external photon fields in the framework of one-zone scenario. Following derived analytical results, we fit the SEDs of these LBLs with the conventional one-zone leptonic model and study their jet physical properties. Our main results can be summarized as follows. (1) We find that most LBLs cannot be fitted by the one-zone synchrotron-self-Compton (SSC) model. This indicates that external photons play a crucial role in the high-energy emission of LBLs, therefore we suggest that LBLs are masquerading BL Lacs. (2) We suggest that the γ-ray emitting regions of LBLs are located outside the broad-line region and within the dusty torus. (3) By extending the analytical method to all types of LSPs in Fermi -4LAC (using historical data), we find that the high-energy peaks of some flat spectrum radio quasars and blazar candidates of unknown types can be attributed to the SSC emission, implying that the importance of external photons could be minor. We suggest that the variability time-scale may help distinguish the origin of the high-energy peak. [ABSTRACT FROM AUTHOR]