Magnetic reconnection is invoked as an efficient particle accelerator in a variety of astrophysical sources of non-thermal high-energy radiation. With large-scale two-dimensional particle-in-cell simulations of relativistic reconnection (i.e. with magnetization σ ≫ 1) in pair plasmas, we study the long-term evolution of the power-law slope and high-energy cut-off of the spectrum of accelerated particles. We find that the high-energy spectral cut-off does not saturate at γcut ∼ 4σ, as claimed by earlier studies, but it steadily grows with time as long as the reconnection process stays active. At late times, the cut-off scales approximately as γ _cut∝ √{t}, regardless of the flow magnetization and initial temperature. We show that the particles dominating the high-energy spectral cut-off reside in plasmoids, and in particular in a strongly magnetized ring around the plasmoid core. The growth of their energy is driven by the increase in the local field strength, coupled with the conservation of the first adiabatic invariant. We also find that the power-law slope of the spectrum (p = -d log N/d log γ) evolves with time. For σ ≳ 10, the spectrum is hard at early times (p ≲ 2), but it tends to asymptote to p ∼ 2; the steepening of the power-law slope allows the spectral cut-off to extend to higher and higher energies, without violating the fixed energy budget of the system. Our results demonstrate that relativistic reconnection is a viable candidate for accelerating the high-energy particles emitting in relativistic astrophysical sources.

Frequency-dependent brightness fluctuations of radio sources, the so-called extreme scattering events (ESEs), have been observed over the last three decades. They are caused by Galactic plasma structures whose geometry and origin are still poorly understood. In this paper, we construct axisymmentric two-dimensional (2D) column density profiles for the plasma lens and explore the resulting ESEs for both point-like and extended sources. A quantity that becomes relevant is the impact parameter b, namely the distance of the observer's path from the lens' symmetry axis. We demonstrate its effects on the shape of ESE light curves and use it for a phenomenological classification of ESEs into four main types. Three of them are unique outcomes of the 2D model and do not show a characteristic U-shaped dip in the light curve, which has been traditionally used as an identification means of ESEs. We apply our model to five well-studied ESEs and show that elongated plasma tubes or quasi-spherical clouds are favoured over plasma sheets for four of them, while the remaining one is compatible with both lens geometries.

Blazars are a sub-category of radio-loud active galactic nuclei with relativistic jets pointing towards to the observer. They are well-known for their non-thermal variable emission, which practically extends over the whole electromagnetic spectrum. Despite the plethora of multi-wavelength observations, the issue about the origin of the γ-ray and radio emission in blazar jets remains unsettled. Here, we construct a parametric leptonic model for studying the connection between the γ-ray and radio emission in both steady-state and flaring states of blazars. Assuming that relativistic electrons are injected continuously at a fixed distance from the black hole, we numerically study the evolution of their population as it propagates to larger distances while losing energy due to expansion and radiative cooling. In this framework, γ-ray photons are naturally produced at small distances (e.g., 10‑3 pc) when the electrons are still very energetic, whereas the radio emission is produced at larger distances (e.g., 1 pc), after the electrons have cooled and the emitting region has become optically thin to synchrotron self-absorption due to expansion. We present preliminary results of our numerical investigation for the steady-state jet emission and the predicted time lags between γ-rays and radio during flares.

We consider implications of high-energy neutrino emission from blazar flares, including the recent event IceCube-170922A and the 2014-2015 neutrino flare that could originate from TXS 0506+056. First, we discuss their contribution to the diffuse neutrino intensity taking into account various observational constraints. Blazars are likely to be subdominant in the diffuse neutrino intensity at sub-PeV energies, and we show that blazar flares like those of TXS 0506+056 could make ≲1%-10% of the total neutrino intensity. We also argue that the neutrino output of blazars can be dominated by the flares in the standard leptonic scenario for their γ-ray emission, and energetic flares may still be detected with a rate of ≲ 1 {yr}}-1. Second, we consider multi-messenger constraints on the source modeling. We show that luminous neutrino flares should be accompanied by luminous broadband cascade emission, emerging also in X-rays and γ-rays. This implies that not only γ-ray telescopes like Fermi but also X-ray sky monitors such as Swift and MAXI are critical to test the canonical picture based on the single-zone modeling. We also suggest a two-zone model that can naturally satisfy the X-ray constraints while explaining the flaring neutrinos via either photomeson or hadronuclear processes.

Detection of the IceCube-170922A neutrino coincident with the flaring blazar TXS 0506+056, the first and only ∼3σ high-energy neutrino source association to date, offers a potential breakthrough in our understanding of high-energy cosmic particles and blazar physics. We present a comprehensive analysis of TXS 0506+056 during its flaring state, using newly collected Swift, NuSTAR, and X-shooter data with Fermi observations and numerical models to constrain the blazar’s particle acceleration processes and multimessenger (electromagnetic (EM) and high-energy neutrino) emissions. Accounting properly for EM cascades in the emission region, we find a physically consistent picture only within a hybrid leptonic scenario, with γ-rays produced by external inverse-Compton processes and high-energy neutrinos via a radiatively subdominant hadronic component. We derive robust constraints on the blazar’s neutrino and cosmic-ray emissions and demonstrate that, because of cascade effects, the 0.1-100 keV emissions of TXS 0506+056 serve as a better probe of its hadronic acceleration and high-energy neutrino production processes than its GeV-TeV emissions. If the IceCube neutrino association holds, physical conditions in the TXS 0506+056 jet must be close to optimal for high-energy neutrino production, and are not favorable for ultrahigh-energy cosmic-ray acceleration. Alternatively, the challenges we identify in generating a significant rate of IceCube neutrino detections from TXS 0506+056 may disfavor single-zone models, in which γ-rays and high-energy neutrinos are produced in a single emission region. In concert with continued operations of the high-energy neutrino observatories, we advocate regular X-ray monitoring of TXS 0506+056 and other blazars in order to test single-zone blazar emission models, clarify the nature and extent of their hadronic acceleration processes, and carry out the most sensitive possible search for additional multimessenger sources.

A unique and often overlooked property of a source loaded with relativistic protons is that it can become supercritical, i.e. it can undergo an abrupt transition from a radiatively inefficient to a radiatively efficient state once its proton energy density exceeds a certain threshold. In this paper, we investigate the temporal variability of hadronic systems in this hardly explored regime. We show that there exists a range of proton densities that prevent the system from reaching a steady state, but drive it instead in a quasi-periodic mode. The escaping radiation then exhibits limit cycles, even if all physical parameters are held constant in time. We extend our analysis to cases where the proton injection rate varies with time and explore the variability patterns of escaping radiation as the system moves in and out from the supercritical regime. We examine the relevance of our results to the variability of the prompt gamma-ray burst emission and show that, at least on a phenomenological level, some interesting analogies exist.

We report an optical/UV jet and counterjet in M84, previously unreported in archival Hubble Space Telescope imaging. With archival VLA, ALMA, and Chandra imaging, we examine the first well-sampled spectral energy distribution of the inner jet of M84, where we find that multiple co-spatial spectral components are required. In particular, the ALMA data reveal that the radio spectrum of all four knots in the jet turns over at approximately 100 GHz, which requires a second component for the bright optical/UV emission. Further, the optical/UV has a soft spectrum and is inconsistent with the relatively flat X-ray spectrum, which indicates a third component at higher energies. Using archival VLA imaging, we have measured the proper motion of the innermost knots at 0.9 ± 0.6 and 1.1 ± 0.4c, which when combined with the low jet-to-counterjet flux ratio yields an orientation angle for the system of {74}-18+9°. In the radio, we find high fractional polarization of the inner jet of up to 30% while in the optical no polarization is detected (<8%). We investigate different scenarios for explaining the particular multicomponent spectral energy distribution (SED) of the knots. Inverse Compton models are ruled out due to the extreme departure from equipartition and the unrealistically high total jet power required. The multicomponent SED can be naturally explained within a leptohadronic scenario, but at the cost of very high power in relativistic protons. A two-component synchrotron model remains a viable explanation, but more theoretical work is needed to explain the origin and properties of the electron populations.

Plasmoids, overdense blobs of plasma containing magnetic fields and high-energy particles, are a self-consistent outcome of the reconnection process in the relativistic regime. Recent two-dimensional particle-in-cell (PIC) simulations have shown that plasmoids can undergo a variety of processes (e.g. mergers, bulk acceleration, growth, and advection) within the reconnection layer. We developed a Monte Carlo code, benchmarked with the recent PIC simulations, to examine the effects of these processes on the steady-state size and momentum distributions of the plasmoid chain. The differential plasmoid size distribution is shown to be a power law, ranging from a few plasma skin depths to ∼0.1 of the reconnection layer's length. The power-law slope is shown to be linearly dependent upon the ratio of the plasmoid acceleration and growth rates, which slightly decreases with increasing plasma magnetization. We perform a detailed comparison of our results with those of recent PIC simulations and briefly discuss the astrophysical implications of our findings through the representative case of flaring events from blazar jets.

We report on the X-ray and optical properties of two high-mass X-ray binary systems located in the Large Magellanic Cloud (LMC). Based on the obtained optical spectra, we classify the massive companion as a supergiant star in both systems. Timing analysis of the X-ray events collected by XMM-Newton revealed the presence of coherent pulsations (spin period ∼2013 s) for XMMU J053108.3-690923 and fast flaring behaviour for XMMU J053320.8-684122. The X-ray spectra of both systems can be modelled sufficiently well by an absorbed power law, yielding hard spectra and high intrinsic absorption from the environment of the systems. Due to their combined X-ray and optical properties, we classify both systems as SgXRBs: the 19th confirmed X-ray pulsar and a probable supergiant fast X-ray transient in the LMC, the second such candidate outside our Galaxy.

PSR J2032+4127 is a young and rapidly rotating pulsar on a highly eccentric orbit around the high-mass Be star MT91 213. X-ray monitoring of the binary system over an ∼4000 d period with Swift has revealed an increase of the X-ray luminosity which we attribute to the synchrotron emission of the shocked pulsar wind. We use Swift X-ray observations to infer a clumpy stellar wind with r-2 density profile and constrain the Lorentz factor of the pulsar wind to 105 < γw < 106. We investigate the effects of an axisymmetric stellar wind with polar gradient on the X-ray emission. Comparison of the X-ray light curve hundreds of days before and after the periastron can be used to explore the polar structure of the wind.