We present a toy model for radio emission in high-mass X-ray binaries (HMXBs) with strongly magnetized neutron stars (NSs) where a wind-collision region is formed by the NS outflow and the stellar wind of the massive companion. Radio emission is expected from the synchrotron radiation of shock-accelerated electrons and the free-free emission of the stellar wind. We found that the predicted relation between the GHz luminosity (LR) and the accretion X-ray luminosity (LX) can be written as $L_\mathrm{ R} \propto L_\mathrm{ X}^b$ for most parameters. No correlation with X-rays is expected (b = 0) when the thermal emission of the stellar wind dominates in radio. We typically find a steep correlation (b = 12/7) for sub-Eddington X-ray luminosities and a more shallow one [b = 2(p - 1)/7] for super-Eddington X-ray luminosities, where p is the power-law index of accelerated electrons. The maximum predicted radio luminosity is independent of the NS properties, while it depends on the stellar wind momentum, binary separation distance, and the minimum electron Lorentz factor. Using a Bayesian approach, we modelled the radio observations of Swift J0243.6+6124 that cover a wide range of mass accretion rates. Our results support a shock origin for the radio detections at sub-Eddington X-ray luminosities. However, no physically meaningful parameters could be found for the super-Eddington phase of the outburst, suggesting a different origin. Future observations with more sensitive instruments might reveal a large number of HMXBs with strongly magnetized NSs in radio, allowing the determination of the slope in the LR-LX relation, and putting the wind-collision scenario into test.

The broad high-energy spectral component in blazars is usually attributed to various inverse Compton scattering processes in the relativistic jet, but has not been clearly identified in most cases due to degeneracies in physical models. AP Librae, a low-synchrotron-peaking BL Lac object (LBL) detected in 2015 by H.E.S.S. at very high energies (VHE; >0.5 TeV), has an extremely broad high-energy spectrum, covering ~9 decades in energy. Standard synchrotron self-Compton models generally fail to reproduce the VHE emission, which has led to the suggestion that it might arise not from the blazar core, but on kiloparsec scales from inverse Compton (IC) scattering of cosmic microwave background (CMB) photons by a still-relativistic jet (IC/CMB). IC/CMB models for the TeV emission of AP Librae in prior works have implied a high level of infrared emission from the kiloparsec-scale jet. With newly obtained Hubble Space Telescope (HST) imaging, we obtain a deep upper limit on the kiloparsec-scale jet emission at 1.6 μm, well below the expected level. High-resolution Atacama Large Millimeter/submillimeter Array imaging in bands 3-9 reveals a residual dust-disk signature after core subtraction, with a clearly thermal spectrum, and an extent (~500 pc) that matches with a nonjet residual emission seen after point-spread function subtraction in our 1.6 μm HST imaging. We find that the unusually broad GeV and VHE emission in AP Librae can be reproduced through the combined IC scattering of photons from the CMB and the dust disk, respectively, by electrons in both the blazar core and subkiloparsec jet.

Blazars are the most extreme subclass of active galactic nuclei with relativistic jets emerging from a super-massive black hole and forming a small angle with respect to our line of sight. Blazars are also known to be related to flaring activity as they exhibit large flux variations over a wide range of frequency and on multiple timescales, ranging from a few minutes to several months. The detection of a high-energy neutrino from the flaring blazar TXS 0506+056 and the subsequent discovery of a neutrino excess from the same direction have naturally strengthened the hypothesis that blazars are cosmic neutrino sources. While neutrino production during gamma-ray flares has been widely discussed, the neutrino yield of X-ray flares has received less attention. Motivated by a theoretical scenario where high energy neutrinos are produced by energetic protons interacting with their own X-ray synchrotron radiation, we make neutrino predictions over a sample of a sample of X-ray blazars. This sample consists of all blazars observed with the X-ray Telescope (XRT) on board Swift more than 50 times from November 2004 to November 2020. The statistical identification of a flaring state is done using the Bayesian Block algorithm to the 1 keV XRT light curves of frequently observed blazars. We categorize flaring states into classes based on their variation from the time-average value of the data points. During each flaring state, we compute the expected muon plus anti-muon neutrino events as well as the total signal for each source using the point-source effective area of Icecube for different operational seasons. We find that the median of the total neutrino number (in logarithm) from flares with duration $<30$ d is $\mathcal{N}^{(\rm tot)}_{\nu_{\mu}+\bar{\nu}_{\mu}} \sim 0.02$.

Blazar hadronic models have been developed in the past decades as an alternative to leptonic ones. In hadronic models the gamma-ray emission is associated with synchrotron emission by protons, and/or secondary leptons produced in proton-photon interactions. Together with photons, hadronic emission models predict the emission of neutrinos that are therefore the smoking gun for acceleration of relativistic hadrons in blazar jets. The simulation of proton-photon interactions and all associated radiative processes is a complex numerical task, and different approaches to the problem have been adopted in the literature. So far, no systematic comparison between the different codes has been performed, preventing a clear understanding of the underlying uncertainties in the numerical simulations. To fill this gap, we have undertaken the first comprehensive comparison of blazar hadronic codes, and the results from this effort will be presented in this contribution.

Several models have been suggested to explain the fast gamma-ray variability observed in blazars, but its origin is still debated. One scenario is magnetic reconnection, a process that can efficiently convert magnetic energy to energy of relativistic particles accelerated in the reconnection layer. In our study, we compare results from state-of-the-art particle-in-cell simulations with observations of blazars at Very High Energy (VHE, E > 100 GeV) gamma-rays. Our goal is to test our model predictions on fast gamma-ray variability with data and to constrain the parameter space of the model, such as the magnetic field strength of the unreconnected plasma and the reconnection layer orientation in the blazar jet. For this first comparison, we used the remarkably well-sampled VHE gamma-ray light curve of Mrk 421 observed with the MAGIC and VERITAS telescopes in 2013.The simulated VHE light curves were generated using the observable parameters of Mrk 421, such as the jet power, bulk Lorentz factor, and the jet viewing angle, and sampled as real data. Our results pave the way for future model-to-data comparison with next-generation Cherenkov telescopes, which will help further constrain the different variability models.

Blazar flares have been suggested as ideal candidates for enhanced neutrino production. While the neutrino signal of γ-ray flares has been widely discussed, the neutrino yield of X-ray flares has received less attention. Here, we compute the predicted neutrino signal from X-ray flares detected in 66 blazars observed more than 50 times with the X-ray Telescope (XRT) on board the Neil Gehrels Swift Observatory. We consider a scenario where X-ray flares are powered by synchrotron radiation of relativistic protons, and neutrinos are produced through photomeson interactions between protons with their own synchrotron X-ray photons. Using the 1 keV X-ray light curves for flare identification, the 0.5-10 keV fluence of each flare as a proxy for the all-flavour neutrino fluence, and the IceCube point-source effective area for different detector configurations, we calculate the number of muon and antimuon neutrinos above 100 TeV expected for IceCube from each flaring source. The bulk of the neutrino events from the sample originates from flares with durations ~1-10 d. Accounting for the X-ray flare duty cycle of the sources in the sample, which ranges between ~2 and 24 per cent, we compute an average yearly neutrino rate for each source. The median of the distribution (in logarithm) is ~0.03 yr-1, with Mkn 421 having the highest predicted rate 1.2 ± 0.3 yr-1, followed by 3C 273 (0.33 ± 0.03 yr-1) and PG 1553+113 (0.25 ± 0.02 yr-1). Next-generation neutrino detectors together with regular X-ray monitoring of blazars could constrain the duty cycle of hadronic X-ray flares.

Eight years after the first detection of high-energy astrophysical neutrinos by IceCube, we are still almost clueless as regards to their origin, although the case for blazars being neutrino sources is getting stronger. After the first significant association at the $3\!-\!3.5\, \sigma$ level in time and space with IceCube neutrinos, i.e. the blazar TXS 0506+056 at z = 0.3365, some of us have in fact selected a unique sample of 47 blazars, out of which ~16 could be associated with individual neutrino track events detected by IceCube. Building upon our recent spectroscopy work on these objects, here we characterize them to determine their real nature and check if they are different from the rest of the blazar population. For the first time we also present a systematic study of the frequency of masquerading BL Lacs, i.e. flat-spectrum radio quasars with their broad lines swamped by non-thermal jet emission, in a γ-ray- and IceCube-selected sample, finding a fraction >24 per cent and possibly as high as 80 per cent. In terms of their broad-band properties, our sources appear to be indistinguishable from the rest of the blazar population. We also discuss two theoretical scenarios for neutrino emission, one in which neutrinos are produced in interactions of protons with jet photons and one in which the target photons are from the broad-line region. Both scenarios can equally account for the neutrino-blazar correlation observed by some of us. Future observations with neutrino telescopes and X-ray satellites will test them out.

Blazars are a class of active galactic nuclei that host relativistic jets oriented close to the observer's line of sight. Blazars have very complex variability properties. Flares, namely flux variations around the mean value with a well-defined shape and duration, are one of the identifying properties of the blazar phenomenon. Blazars are known to exhibit multiwavelength flares, but also "orphan" flares, namely flux changes that appear only in a specific energy range. Various models, sometimes at odds with each other, have been proposed to explain specific flares even for a single source, and cannot be synthesized into a coherent picture. In this paper, we propose a unified model for explaining orphan and multiwavelength flares from blazars in a common framework. We assume that the blazar emission during a flare consists of two components: (i) a quasistable component that arises from the superposition of numerous but comparatively weak dissipation zones along the jet, forming the background (low-state) emission of the blazar, and (ii) a transient component, which is responsible for the sudden enhancement of the blazar flux, forming at a random distance along the jet by a strong energy dissipation event. Whether a multiwavelength or orphan flare is emitted depends on the distance from the base of the jet where the dissipation occurs. Generally speaking, if the dissipation occurs at a small/large distance from the supermassive black hole, the inverse Compton/synchrotron radiation dominates and an orphan γ -ray/optical flare tends to appear. On the other hand, we may expect a multiwavelength flare if the dissipation occurs at an intermediate distance. We show that the model can successfully describe the spectral energy distribution of different flares from the flat spectrum radio quasar 3C 279 and the BL Lacertae object PKS 2155-304.

The origin of Petaelectronvolt (PeV) astrophysical neutrinos is fundamental to our understanding of the high-energy Universe. Apart from the technical challenges of operating detectors deep below ice, oceans, and lakes, the phenomenological challenges are even greater than those of gravitational waves; the sources are unknown, hard to predict, and we lack clear signatures. Neutrino astronomy therefore represents the greatest challenge faced by the astronomy and physics communities thus far. The possible neutrino sources range from accretion disks and tidal disruption events, to relativistic jets and galaxy clusters with blazar TXS~0506+056 the most compelling association thus far. Since that association, immense effort has been put into proving or disproving that jets are indeed neutrino emitters, but to no avail. By generating simulated neutrino counterpart samples, we explore the potential of detecting a significant correlation of neutrinos with jets from active galactic nuclei. We find that, given the existing challenges, even our best experiments could not have produced a $>3\sigma$ result. Larger programs over the next few years will be able to detect a significant correlation only if the brightest radio sources, rather than all jetted active galactic nuclei, are neutrino emitters. We discuss the necessary strategies required to steer future efforts into successful experiments.

X-Ray pulsars are systems powered by accretion, the majority of which is found in Be X-ray binaries (BeXRBs). The study of Giant outbursts (L$_{X}$ > 10$^{38}$ erg s$ ^{-1}$) in such systems becomes very relevant in the advent of the recent discoveries of pulsating ultra-luminous X-ray sources (PULXs) with an apparent isotropic luminosity above the Eddington limit for a typical neutron star (NS) demonstrating that stable accretion onto NSs is possible at super - Eddington rates. Given that PULXs host magnetized NSs, several attempts have been made to estimate the magnetic field of the NS using standard torque models. At the same time theoretical studies have demonstrated that it is required to adjust these models due to changes in the accretion disc structure when exceeding the Eddington limit. Motivated by these findings we studied torque models during Giant outbursts of BeXRBs monitored by Fermi/GBM and Swift/BAT. We developed a code to estimate posterior distributions for the parameters of standard accretion models and binary orbital parameters using a nested sampling algorithm for Bayesian Parameter Estimation. Most notably we applied our method to the recently discovered Swift J0243.6+6124 (i.e., the only known Galactic PULX) and we illustrate that the standard torque models need adjustment to explain the observed spin evolution of the NS. Finally, we discuss the implications the newest GAIA distances have on the NS equation of state.

The power-law X-ray spectrum of black-hole X-ray Binaries (BHXRBs) in the hard and hard-intermediate spectral states most probably arises from the interaction of the hot corona with the cold accretion disk. Soft photons from the disk are Comptonized in the corona and some Comptonized photons return to the disk, where they are reprocessed back into soft photons. This system is inherently non-linear as the soft radiation that cools the electrons of the corona is partially produced by the same particle population. We write simple time-dependent equations that describe energy conservation in the system corona-accretion disk and we show that these are of the Lotka-Volterra type. Solving them for constant mass accretion rate, we find that the energy densities of the hot electrons and the Comptonized photons exhibit oscillations that die out after a few cycles. A crucial requirement for the emergence of transient oscillations is that most of the accretion power goes into the corona and the disk is cold. The characteristic frequencies of these oscillations, which depend on the mass accretion rate, the corona size, and the reflection fraction, lie in the range of 1-60 Hz for typical parameter values. We also demonstrate that the natural frequencies persist even in the case of variable accretion rates. We make, therefore, a tentative claim that type-C Quasi Periodic Oscillations in the hard state of BHXRBs could, in principle, arise from the interaction of the hot Comptonizing corona with the much colder accretion disk.

A correlation between high-energy cosmic neutrinos and active galaxies with relativistic jets, known as blazars, has been remarked by several authors in the past few years. It has also been suggested that blazar flares (i.e. significant changes in flux) are ideal candidates for enhanced neutrino emission. So far, several works have discussed the connection between the X-ray and γ-ray flares with the neutrino signal, while many blazars that are positionally associated with high energy neutrinos detected by IceCube exhibit strong parsec-scale radio emission. In this work, we discuss the connection between the X-ray and radio frequencies during periods of enhanced neutrino emission. In our model, injection of relativistic hadrons into a compact region of the jet produces an X-ray and high-energy neutrino flare. At the same time, relativistic electron and positrons are produced through photo-hadronic interactions and photon-photon pair production (secondary pairs). We investigate the evolution of the electromagnetic cascade by following the secondary pairs upon their injection into an expanding spherical blob. Initially, the radio part of the cascade spectrum is self-absorbed, but the source becomes transparent to GHz frequencies as the blob expands, producing a delayed "radio flare". We investigate how the time lags between the X-ray and the radio flares are formed and what is the effect of physical parameters, such as the expansion velocity V$ _{exp}$ and the magnetic field strength B.

In the most powerful astrophysical sources, reconnection and turbulence operate in the "relativistic" regime, where the magnetic field energy exceeds even the rest mass energy of the plasma. Here, reconnection and turbulence can lead to fast dissipation rates and efficient particle acceleration, thus being prime candidates for powering the observed fast and bright flares of high-energy non-thermal emission. With fully-kinetic particle-in-cell (PIC) simulations and analytical theory, we investigate the physics of relativistic reconnection and turbulence, and demonstrate that they can be the "engines" behind: (1) high-energy flares in blazar jets; and (2) the hard-state spectra of black hole X-ray binaries and Active Galactic Nuclei.

Context. Black-hole X-ray binaries (BHXRBs) in the hard and hard-intermediate spectral (and temporal) states exhibit in their power spectra characteristic frequencies called type-C quasi-periodic oscillations (QPOs). Various models that can explain them with various degrees of success have been proposed, but a definitive answer is still missing. Aims: The hot Comptonizing corona interacting with the cold accretion disk, both of which are central in understanding BHXRBs, is essentially a dynamical system. Our aim is to investigate if the radiative coupling between the two components can produce QPOs. Methods: We write and solve the time-dependent equations that describe energy conservation in the system corona - accretion disk. We examine both constant and variable mass accretion rates. By necessity, in this first investigation we use a simple model, but it contains all the essential ingredients. Results: For a constant mass accretion rate and certain justifiable conditions, the dynamic corona - disk system exhibits oscillations, which die out after a few cycles. The characteristic frequencies of these oscillations are similar to the ones observed in the power spectra of BHXRBs. For most parameters, the natural frequencies persist even in the case of variable accretion rates. Conclusions: We argue that type-C QPOs in BHXRBs could, in principle, arise from the interaction of the hot Comptonizing corona with the much colder accretion disk. If this picture is correct, it has immediate implications for other systems that contain the above constituents, such as active galactic nuclei.

Network meta-analysis (NMA) is widely used in evidence synthesis to estimate the effects of several competing interventions for a given clinical condition. One of the challenges is that it is not possible in disconnected networks. Component network meta-analysis (CNMA) allows technically 'reconnecting' a disconnected network with multicomponent interventions. The additive CNMA model assumes that the effect of any multicomponent intervention is the additive sum of its components. This assumption can be relaxed by adding interaction component terms, which improves the goodness of fit but decreases the network connectivity. Model selection aims at finding the model with a reasonable balance between the goodness of fit and connectivity (selected CNMA model). We aim to introduce a forward model selection strategy for CNMA models and to investigate the performance of CNMA models for connected and disconnected networks. We applied the methods to a real Cochrane review dataset and simulated data with additive, mildly, or strongly violated intervention effects. We started with connected networks, and we artificially constructed disconnected networks. We compared the results of the additive and the selected CNMAs from each connected and disconnected network with the NMA using the mean squared error and coverage probability. CNMA models provide good performance for connected networks and can be an alternative to standard NMA if additivity holds. On the contrary, model selection does not perform well for disconnected networks, and we recommend conducting separate analyses of subnetworks.

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 (HSC) 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 multi-pulse light curves from slowly expanding regions ($u_{\rm exp}\lesssim 0.01 c)$ 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 $\sim 1$ MeV ($\sim 1$ GeV) for maximum proton energies $\sim 1-10$ PeV ($1-10$ EeV) assuming a jet Lorentz factor 100. Peak $\gamma$-ray luminosities are in the range $10^{49}-10^{52}$ erg s$^{-1}$, with the MeV-peaked spectra being $\sim 100-300$ times more luminous than their GeV-peaked analogues. HSC bursts peaking in the MeV are also copious $\sim 10$ TeV neutrino emitters, with an all-flavour fluence $\sim 10\%$ of the $\gamma$-ray one. The hypothesis that typical long-duration GRBs are powered by HSC could be tested in the near future with more sensitive neutrino telescopes like IceCube-Gen2.

PKS 0735+178 is a bright radio and $\gamma$-ray blazar that is possibly associated with multiple neutrino events observed by the IceCube, Baikal, Baksan, and KM3NeT neutrino telescopes. The source was found to undergo a major flaring activity in $\gamma$-ray, X-ray, ultraviolet (UV) and optical bands. We present a long-term detailed study of this peculiar blazar to investigate the temporal and spectral changes in the multi-wavelength 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\times10^3$ seconds during the observation on December 17, 2021. In the $\gamma$-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 $\gamma$-ray powers, and parsec scale jet properties. Moreover, PKS 0735+178, like all the others, is 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 $\sim$ PeV protons with an external UV photon field. This scenario predicts $\sim 0.067$ muon and antimuon neutrinos over the observed 3-week flare. Our results are consistent with the detection of one very-high-energy neutrino like IceCube-211208A.