We present a multi-messenger model for the prompt emission from GRB 221009A within the internal shock scenario. We consider the time-dependent evolution of the outflow with its impact on the observed light curve from multiple collisions, as well as the self-consistent generation of the electromagnetic spectrum in synchrotron and inverse Compton-dominated scenarios. Our lepto-hadronic model includes UHE protons potentially accelerated in the outflow, and their feedback on spectral energy distribution and on the neutrino emission. We find that we can roughly reproduce the observed light curves with an engine with varying ejection velocity of ultrarelativistic material, which has an intermediate quiescent period of about 200 s and a variability timescale of ~1 s. We consider baryonic loadings of 3 and 30 that are compatible with the hypothesis that the highest-energetic LHAASO photons might come from UHECR interactions with the extragalactic background light, and the paradigm that energetic GRBs may power the UHECR flux. For these values and the high dissipation radii considered, we find consistency with the nonobservation of neutrinos and no significant signatures on the electromagnetic spectrum. Inverse Compton-dominated scenarios from the prompt emission are demonstrated to lead to about an order of magnitude higher fluxes in the HE range; this enhancement is testable via its spectral impact in the Fermi-GBM and LAT ranges.

Blazars are a rare class of active galactic nuclei (AGNs) with relativistic jets pointing towards the observer. Jets are thought to be launched as Poynting-flux dominated outflows that accelerate to relativistic speeds at the expense of the available magnetic energy. In this work, we consider electron-proton jets and assume that particles are energized via magnetic reconnection in parts of the jet where the magnetization is still high (σ ≥ 1). The magnetization and bulk Lorentz factor Γ are related to the available jet energy per baryon as μ = Γ(1 + σ). We adopt an observationally motivated relation between Γ and the mass accretion rate into the black hole $\dot{m}$, which also controls the luminosity of external radiation fields. We numerically compute the photon and neutrino jet emission as a function of μ and σ. We find that the blazar SED is produced by synchrotron and inverse Compton radiation of accelerated electrons, while the emission of hadronic-related processes is subdominant except for the highest magnetization considered. We show that low-luminosity blazars (Lγ ≲ 1045 erg s-1) are associated with less powerful, slower jets with higher magnetizations in the jet dissipation region. Their broad-band photon spectra resemble those of BL Lac objects, and the expected neutrino luminosity is $L_{\nu +\bar{\nu }}\sim (0.3-1)\, L_{\gamma }$. High-luminosity blazars (Lγ ≫ 1045 erg s-1) are associated with more powerful, faster jets with lower magnetizations. Their broad-band photon spectra resemble those of flat spectrum radio quasars, and they are expected to be dim neutrino sources with $L_{\nu +\bar{\nu }}\ll L_{\gamma }$.

For over a century, the identification of high-energy cosmic ray (CR) sources remains an open question. For Galactic CRs with energy up to 1015 eV, supernova remnants (SNRs) have traditionally been thought the main candidate source. However, recent TeV γ-ray observations have questioned the SNR paradigm. Propagating CRs are deflected by the Galactic magnetic field, hence, γ-rays and neutrinos produced via inelastic hadronic interactions are the only means for unveiling the CR sources. In this work, we study the γ-ray and neutrino emission produced by CRs accelerated inside Galactic jets of stellar-mass black holes in X-ray binaries (BHXBs). We calculate the intrinsic neutrino emission of two prototypical BHXBs , Cygnus X-1 and GX 339-4, for which we have high-quality, quasi-simultaneous multiwavelength spectra. Based on these prototypical sources, we discuss the likelihood of the 35 known Galactic BHXBs to be efficient CR accelerators. Moreover, we estimate the potential contribution to the CR spectrum of a viable population of BHXBs that reside in the Galactic plane. When these BHXBs go into outburst, they may accelerate particles up to hundreds of TeV that contribute to the diffuse γ-ray and neutrino spectra while propagating in the Galactic medium. Using HERMES, an open-source code that calculates the hadronic processes along the line of sight, we discuss the contribution of BHXBs to the diffuse γ-ray and neutrino fluxes, and compare these to their intrinsic γ-ray and neutrino emissions. Finally, we discuss the contribution of BHXBs to the observed spectrum of Galactic CRs.

Gamma-ray flares of blazars may be accompanied by high-energy neutrinos due to interactions of high-energy cosmic rays in the jet with photons, as suggested by the detection of the high-energy neutrino IceCube-170922A during a major gamma-ray flare from blazar TXS 0506+056 at the ~3σ significance level. In this work, we present a statistical study of gamma-ray emission from blazars to constrain the contribution of gamma-ray flares to their neutrino output. We construct weekly binned light curves for 145 gamma-ray bright blazars in the Fermi Large Area Telescope Monitored Source List adding TXS 0506+056. We derive the fraction of time spent in the flaring state (flare duty cycle) and the fraction of energy released during each flare from the light curves with a Bayesian blocks algorithm. We find that blazars with lower flare duty cycles and energy fractions are more numerous among our sample. We identify a significant difference in flare duty cycles between blazar subclasses at a significance level of 5%. Then using a general scaling relation for the neutrino and gamma-ray luminosities, ${L}_{\nu }\propto {({L}_{\gamma })}^{\gamma }$ with a weighting exponent of γ = 1.0-2.0, normalized to the quiescent gamma-ray or X-ray flux of each blazar, we evaluate the neutrino energy flux of each gamma-ray flare. The gamma-ray flare distribution indicates that blazar neutrino emission may be dominated by flares for γ ≳ 1.5. The neutrino energy fluxes for 1 week and 10 yr bins are compared with the decl.-dependent IceCube sensitivity to constrain the standard neutrino emission models for gamma-ray flares. Finally, we present the upper-limit contribution of blazar gamma-ray flares to the isotropic diffuse neutrino flux.

Relativistic hadronic plasmas can become, under certain conditions, supercritical, abruptly and efficiently releasing the energy stored in protons through photon outbursts. Past studies have tried to relate the features of such hadronic supercriticalities (HSCs) to the phenomenology of gamma-ray burst (GRB) prompt emission. In this work we investigate, for the first time, HSC in adiabatically expanding sources. We examine the conditions required to trigger HSC, study the role of expansion velocity, and discuss our results in relation to GRB prompt emission. We find multipulse light curves from slowly expanding regions (≲ 0.01c) that are a manifestation of the natural HSC quasi-periodicity, while single-pulse light curves with a fast rise and slow decay are found for higher velocities. The formation of the photon spectrum is governed by an in-source electromagnetic cascade. The peak photon energy is approximately $1 \cdot \frac{\Gamma }{100} \frac{1+z}{3}$ MeV for maximum proton energies $(1-10) \cdot \frac{\Gamma }{100} \frac{1+z}{3}$ PeV, while the peak γ-ray luminosities are in the range $(10^{49}-10^{52}) \cdot (\frac{\Gamma }{100})^4$ erg s-1. HSC bursts peaking in the MeV energy band are also copious neutrino emitters with peak energies $\sim 10 \cdot \frac{\Gamma }{100} \frac{1+z}{3}$ TeV and an all-flavour neutrino fluence $\sim 10~{{\ \rm per\ cent}}$ of the γ-ray one. The hypothesis that long-duration GRBs are powered by HSCs could be applied therefore only to the most luminous GRBs observed assuming bulk Lorentz factors Γ ≤ 100.

For a subpopulation of energetic gamma-ray bursts (GRBs), a moderate baryonic loading may suffice to power ultra-high-energy cosmic rays (UHECRs). Motivated by this, we study the radiative signatures of cosmic-ray protons in the prompt phase of energetic GRBs. Our framework is the internal shock model with multicollision descriptions of the relativistic ejecta (with different emission regions along the jet), plus time-dependent calculations of photon and neutrino spectra. Our GRB prototypes are motivated by Fermi-Large Area Telescope-detected GRBs (including GRB 221009A) for which further, owing to the large energy flux, neutrino nonobservation of single events may pose a strong limit on the baryonic loading. We study the feedback of protons on electromagnetic spectra in synchrotron- and inverse Compton-dominated scenarios to identify the multiwavelength signatures, to constrain the maximally allowed baryonic loading, and to point out the differences between hadronic and inverse Compton signatures. We find that hadronic signatures appear as correlated flux increases in the optical-UV to soft X-ray and GeV-TeV gamma-ray ranges in the synchrotron scenarios, whereas they are difficult to identify in inverse Compton-dominated scenarios. We demonstrate that baryonic loadings around 10, which satisfy the UHECR energetic requirements, do not distort the predicted photon spectra in the Fermi Gamma-Ray Burst Monitor range and are consistent with constraints from neutrino data if the collision radii are large enough (i.e., the time variability is not too short). It therefore seems plausible that under the condition of large dissipation radii a population of energetic GRBs can be the origin of the UHECRs.

We present the first results for the dust-scattering rings of GRB 221009A, coined as the gamma-ray bursts (GRBs) of the century, as observed by the Neil Gehrels Swift Observatory. We perform analysis of both time resolved observations and stacked data. The former approach enable us to study the expansion of the most prominent rings, associate their origin with the prompt X-ray emission of the GRB and determine the location of the dust layers. The stacked radial profiles increase the signal-to-noise ratio of the data and allows detection of fainter and overlapping peaks in the angular profile. We find a total of 16 dust concentrations (with hints of even more) that span about 15kpc in depth and could be responsible for the highly structured X-ray angular profiles. By comparing the relative scattered fluxes of the five most prominent rings we show that the layer with the largest amount of dust is located at about 0.44 kpc away from us. We finally compare the location of the dust layers with results from experiments that study the 3D structure of our Galaxy via extinction or CO radio observations, and highlight the complementarity of dust X-ray tomography to these approaches.

In this study, we present a method to estimate posterior distributions for standard accretion torque model parameters and binary orbital parameters for X-ray binaries using a nested sampling algorithm for Bayesian parameter estimation. We study the spin evolution of two Be X-ray binary systems in the Magellanic Clouds, RX J0520.5-6932 and RX J0209-7427, during major outbursts, in which they surpassed the Eddington limit. Moreover, we apply our method to the recently discovered Swift J0243.6+6124, the only known Galactic pulsating ultra-luminous X-ray source. This is an excellent candidate for studying the disc evolution at super-Eddington accretion rates, because its luminosity spans several orders of magnitude during its outburst, with a maximum LX that exceeded the Eddington limit by a factor of ~10. Our method, when applied to RX J0520.5-6932 and RX J0209-7427, is able to identify the more favourable torque model for each system, while yielding meaningful ranges for the NS and orbital parameters. Our analysis for Swift J0243.6+6124 illustrates that, contrary to the standard torque model predictions, the magnetospheric radius (Rm) and the Alfvén radius (RA) are not proportional to each other when surpassing the Eddington limit. Reported distance estimates of this source range between 5 and 7 kpc. Smaller distances require non-typical neutron star properties (i.e. mass and radius) and possibly lower radiative efficiency of the accretion column.

The blazar PKS 0735+178 is possibly associated with multiple neutrino events observed by the IceCube, Baikal, Baksan, and KM3NeT neutrino telescopes while it was flaring in the γ-ray, X-ray, ultraviolet, and optical bands. We present a detailed study of this peculiar blazar to investigate the temporal and spectral changes in the multiwavelength emission when the neutrino events were observed. The analysis of Swift-XRT snapshots reveal a flux variability of more than a factor 2 in about 5 × 103 s during the observation on 2021 December 17. In the γ-ray band, the source was in its historical highest flux level at the time of the arrival of the neutrinos. The observational comparison between PKS 0735+178 and other neutrino source candidates, such as TXS 0506+056, PKS 1424+240, and GB6 J1542+6129, shows that all these sources share similar spectral energy distributions, very high radio and γ-ray powers, and parsec scale jet properties. Moreover, we present strong supporting evidence for PKS 0735+178 to be, like all the others, a masquerading BL Lac. We perform comprehensive modelling of the multiwavelength emission from PKS 0735+178 within one-zone lepto-hadronic models considering both internal and external photon fields and estimate the expected accompanying neutrino flux. The most optimistic scenario invokes a jet with luminosity close to the Eddington value and the interactions of ~ PeV protons with an external UV photon field. This scenario predicts ~0.067 muon and anti-muon neutrinos over the observed 3-week flare. Our results are consistent with the detection of one very high-energy neutrino like IceCube-211208A.

Centaurus A is one of the closest radio galaxies to Earth. Its proximity allowed us to extensively study its active galactic nucleus but the core emission mechanism remains elusive because of local strong dust and gas obscuration. The capability of polarimetry to shave-off contaminating emission has been exploited without success in the near-infrared by previous studies but the very recent measurement of the 2-8 keV polarization by the Imaging X-ray Polarimetry Explorer (IXPE) brought the question back to the fore. To determine what is the prevalent photon generation mechanism to the multiwavelength emission from the core of Centaurus A, we retrieved from the archives the panchromatic polarization measurements of the central compact component. We built the total and polarized flux spectral energy distributions of the core and demonstrated that synchrotron self-Compton models nicely fit the polarized flux from the radio to the X-ray band. The linear polarization of the synchrotron continuum is perpendicular to the jet radio axis from the optical to the radio band, and parallel to it at higher energies. The observed smooth rotation of the polarization angle in the ultraviolet band is attributed to synchrotron emission from regions that are getting closer to the particle acceleration site, where the orientation of the jet's magnetic fields become perpendicular to the jet axis. This phenomenon support the shock acceleration mechanism for particle acceleration in Centaurus A, in line with IXPE observations of several high-synchrotron peak blazars.

Magnetic reconnection is often invoked as a source of high-energy particles, and in relativistic astrophysical systems it is regarded as a prime candidate for powering fast and bright flares. We present a novel analytical model-supported and benchmarked with large-scale three-dimensional kinetic particle-in-cell simulations in electron-positron plasmas-that elucidates the physics governing the generation of power-law energy spectra in relativistic reconnection. Particles with Lorentz factor γ ≳ 3σ (here, σ is the magnetization) gain most of their energy in the inflow region, while meandering between the two sides of the reconnection layer. Their acceleration time is ${t}_{\mathrm{acc}}\sim \gamma \,{\eta }_{\mathrm{rec}}^{-1}{\omega }_{{\rm{c}}}^{-1}\simeq 20\,\gamma \,{\omega }_{{\rm{c}}}^{-1}$ , where η rec ≃ 0.06 is the inflow speed in units of the speed of light and ω c = eB 0/mc is the gyrofrequency in the upstream magnetic field. They leave the region of active energization after t esc, when they get captured by one of the outflowing flux ropes of reconnected plasma. We directly measure t esc in our simulations and find that t esc ~ t acc for σ ≳ few. This leads to a universal (i.e., σ-independent) power-law spectrum ${{dN}}_{\mathrm{free}}/d\gamma \propto {\gamma }^{-1}$ for the particles undergoing active acceleration, and ${dN}/d\gamma \propto {\gamma }^{-2}$ for the overall particle population. Our results help to shed light on the ubiquitous presence of power-law particle and photon spectra in astrophysical nonthermal sources.

The origin of extremely fast variability is one of the long-standing questions in the gamma-ray astronomy of blazars. While many models explain the slower, lower energy variability, they cannot easily account for such fast flares reaching hour-to-minute timescales. Magnetic reconnection, a process where magnetic energy is converted to the acceleration of relativistic particles in the reconnection layer, is a candidate solution to this problem. In this work, we employ state-of-the-art particle-in-cell simulations in a statistical comparison with observations of a flaring episode of a well-known blazar, Mrk 421, at a very high energy (VHE, E > 100 GeV). We tested the predictions of our model by generating simulated VHE light curves that we compared quantitatively with methods that we have developed for a precise evaluation of theoretical and observed data. With our analysis, we can constrain the parameter space of the model, such as the magnetic field strength of the unreconnected plasma, viewing angle and the reconnection layer orientation in the blazar jet. Our analysis favours parameter spaces with magnetic field strength 0.1 G, rather large viewing angles (6 − 8°), and misaligned layer angles, offering a strong candidate explanation for the Doppler crisis often observed in the jets of high synchrotron peaking blazars. Full Tables B.1-B.10 are available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (ftp://130.79.128.5) or via https://cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/678/A140