Current Cherenkov telescopes have identified a population of ultra-high-frequency peaked BL Lac objects (UHBLs), also known as extreme blazars, that exhibit exceptionally hard TeV spectra, including ...1ES 0229+200, 1ES 0347−121, RGB J0710+591, 1ES 1101−232, and 1ES 1218+304. Although one-zone synchrotron-self-Compton (SSC) models have been generally successful in interpreting the high-energy emission observed in other BL Lac objects, they are problematic for UHBLs, necessitating very large Doppler factors and/or extremely high minimum Lorentz factors of the emitting leptonic population. In this context, we have investigated alternative scenarios where hadronic emission processes are important, using a newly developed (lepto-)hadronic numerical code to systematically explore the physical parameters of the emission region that reproduces the observed spectra while avoiding the extreme values encountered in pure SSC models. Assuming a fixed Doppler factor δ = 30, two principal parameter regimes are identified, where the high-energy emission is due to: (1) proton-synchrotron radiation, with magnetic fields B ∼ 1–100 G and maximum proton energies E
p; max ≲ 1019 eV; and (2) synchrotron emission from p–γ-induced cascades as well as SSC emission from primary leptons, with B ∼ 0.1–1 G and E
p; max ≲ 1017 eV. This can be realized with plausible, sub-Eddington values for the total (kinetic plus magnetic) power of the emitting plasma, in contrast to hadronic interpretations for other blazar classes that often warrant highly super-Eddington values.
ABSTRACT
While active galactic nuclei with relativistic jets have long been prime candidates for the origin of extragalactic cosmic rays and neutrinos, the BL Lac object TXS 0506+056 is the first ...astrophysical source observed to be associated with some confidence (∼3σ) with a high-energy neutrino, IceCube-170922A, detected by the IceCube Observatory. The source was found to be active in high-energy γ-rays with Fermi-LAT and in very-high-energy γ-rays with the MAGIC telescopes. To consistently explain the observed neutrino and multiwavelength electromagnetic emission of TXS 0506+056, we investigate in detail single-zone models of leptohadronic emission, assuming cospatial acceleration of electrons and protons in the jet, and synchrotron photons from the electrons as targets for photohadronic neutrino production. The parameter space concerning the physical conditions of the emission region and particle populations is extensively explored for scenarios where the γ-rays are dominated by either (1) proton synchrotron emission or (2) synchrotron-self-Compton emission, with a subdominant but non-negligible contribution from photohadronic cascades in both cases. We find that the latter can be compatible with the neutrino observations, while the former is strongly disfavoured due to the insufficient neutrino production rate.
Context. The wealth of recent data from Imaging Air Cherenkov telescopes (IACTs), ultra-high energy cosmic-ray experiments and neutrino telescopes have fuelled a renewed interest in hadronic emission ...models for γ-loud blazars. Aims. We explore physically plausible solutions for a lepto-hadronic interpretation of the stationary emission from high-frequency peaked BL Lac objects (HBLs). The modelled spectral energy distributions are then searched for specific signatures at very high energies that could help to distinguish the hadronic origin of the emission from a standard leptonic scenario. Methods. By introducing a few basic constraints on parameters of the model, such as assuming the co-acceleration of electrons and protons, we significantly reduced the number of free parameters. We then systematically explored the parameter space of the size of the emission region and its magnetic field for two bright γ-loud HBLs, PKS 2155-304 and Mrk 421. For all solutions close to equipartition between the energy densities of protons and of the magnetic field, and with acceptable jet power and light-crossing timescales, we inspected the spectral hardening in the multi-TeV domain from proton-photon induced cascades and muon-synchrotron emission inside the source. Very-high-energy spectra simulated with the available instrument functions from the future Cherenkov Telescope Array (CTA) were evaluated for detectable features as a function of exposure time, source redshift, and flux level. Results. A range of hadronic scenarios are found to provide satisfactory solutions for the broad band emission of the sources under study. The TeV spectrum can be dominated either by proton-synchrotron emission or by muon-synchrotron emission. The solutions for HBLs cover a parameter space that is distinct from the one found for the most extreme BL Lac objects in an earlier study. Over a large range of model parameters, the spectral hardening due to internal synchrotron-pair cascades, the “cascade bump”, should be detectable for acceptable exposure times with the future CTA for a few nearby and bright HBLs.
ABSTRACT
Various attempts have been made in the literature at describing the origin and the physical mechanisms behind flaring events in blazars with radiative emission models, but detailed ...properties of multiwavelength (MWL) light curves still remain difficult to reproduce. We have developed a versatile radiative code, based on a time-dependent treatment of particle acceleration, escape, and radiative cooling, allowing us to test different scenarios to connect the continuous low-state emission self-consistently with that during flaring states. We consider flares as weak perturbations of the quiescent state and apply this description to the 2010 February MWL flare of Mrk 421, the brightest very high energy (VHE) flare ever detected from this archetypal blazar, focusing on interpretations with a minimum number of free parameters. A general criterion is obtained, which disfavours a one-zone model connecting low and high state under our assumptions. A two-zone model combining physically connected acceleration and emission regions yields a satisfactory interpretation of the available time-dependent MWL light curves and spectra of Mrk 421, although certain details remain difficult to reproduce. The two-zone scenario finally proposed for the complex quiescent and flaring VHE emitting region involves both Fermi-I and Fermi-II acceleration mechanisms, respectively, at the origin of the quiescent and flaring emission.
ABSTRACT
The processes operating in blazar jets are still an open question. Modelling the radiation emanating from an extended part of the jet allows one to capture these processes on all scales. ...Kinetic codes solving the Fokker–Planck equation along the jet flow are well suited to this task, as they can efficiently derive the radiation and particle spectra without the need for computationally demanding plasma physical simulations. Here, we present a new extended hadro-leptonic jet code – ExHaLe-jet– which considers simultaneously the processes of relativistic protons and electrons. Within a pre-set geometry and bulk flow, the particle evolution is derived self-consistently. Highly relativistic secondary electrons (and positrons) are created through γ–γ pair production, Bethe–Heitler pair production, and pion/muon decay. These secondaries are entrained in the jet flow decreasing the ratio of protons to electrons with distance from the jet base. For particle–photon interactions, we consider all internal and many external photon fields, such as the accretion disc, broad-line region, and the dusty torus. The external fields turn out to be the most important source for particle–photon interactions governing the resulting photon and neutrino spectra. In this paper, we present the code and an initial parameter study, while in follow-up works we present extensions of the code and more specific applications.
Context.
Standing and moving shocks in relativistic astrophysical jets are very promising sites for particle acceleration to large Lorentz factors and for the emission from the radio up to the
γ
-ray ...band. They are thought to be responsible for at least part of the observed variability in radio-loud active galactic nuclei.
Aims.
We aim to simulate the interactions of moving shock waves with standing recollimation shocks in structured and magnetized relativistic jets and to characterize the profiles of connected flares in the radio light curve.
Methods.
Using the relativistic magneto-hydrodynamic code
MPI-AMRVAC
and a radiative transfer code in post-processing, we explore the influence of the magnetic-field configuration and transverse stratification of an over-pressured jet on its morphology, on the moving shock dynamics, and on the emitted radio light curve. First, we investigate different large-scale magnetic fields with their effects on the standing shocks and on the stratified jet morphology. Secondly, we study the interaction of a moving shock wave with the standing shocks. We calculated the synthetic synchrotron maps and radio light curves and analyze the variability at two frequencies 1 and 15.3 GHz and for several observation angles. Finally, we compare the characteristics of our simulated light curves with radio flares observed from the blazar 3C 273 with the Owens Valley Radio Observatory and Very Long Baseline Array in the MOJAVE survey between 2008 and 2019.
Results.
We find that in a structured over-pressured relativistic jet, the presence of the large-scale magnetic field structure changes the properties of the standing shock waves and leads to an opening in the jet. The interaction between waves from inner and outer jet components can produce strong standing shocks. When crossing such standing shocks, moving shock waves accompanying overdensities injected in the base of the jet cause very luminous radio flares. The observation of the temporal structure of these flares under different viewing angles probes the jet at different optical depths. At 1 GHz and for small angles, the self-absorption caused by the moving shock wave becomes more important and leads to a drop in the observed flux after it interacts with the brightest standing knot. A weak asymmetry is seen in the shape of the simulated flares, resulting from the remnant emission of the shocked standing shocks. The characteristics of the simulated flares and the correlation of peaks in the light curve with the crossing of moving and standing shocks favor this scenario as an explanation of the observed radio flares of 3C 273.
Context.
Among the blazar class, extreme blazars have exceptionally hard intrinsic X-ray/TeV spectra, and extreme peak energies in their spectral energy distribution (SED). Observational evidence ...suggests that the non-thermal emission from extreme blazars is typically non-variable. All these unique features present a challenging case for blazar emission models, especially regarding those sources with hard TeV spectra.
Aims.
We aim to explore the X-ray and GeV observational features of a variety of extreme blazars, including extreme-TeV, extreme-synchrotron (extreme-Syn), and regular high-frequency-peaked BL Lac objects (HBLs). Furthermore, we aim to test the applicability of various blazar emission models that could explain the very hard TeV spectra.
Methods.
We conducted a detailed spectral analysis of X-ray data collected with
AstroSat
and
Swift
-XRT, along with quasi-simultaneous
γ
-ray data from
Fermi
-LAT, for five sources: 1ES 0120+340, RGB J0710+591, 1ES 1101−232, 1ES 1741+196, and 1ES 2322−409. We took three approaches to modelling the SEDs: (1) a steady-state one-zone synchrotron-self-Compton (SSC) code, (2) another leptonic scenario of co-accelerated electrons and protons on multiple shocks applied to the extreme-TeV sources only (
e
–
p
co-acceleration scenario), and (3) a one-zone hadro-leptonic (O
NE
H
A
L
E
) code. The latter code is used twice to explain the
γ
-ray emission process: proton synchrotron and synchrotron emission of secondary pairs.
Results.
Our X-ray analysis provides well-constrained estimates of the synchrotron peak energies for both 1ES0120+340 and 1ES1741+196. These findings categorise these latter objects as extreme-synchrotron sources, as they consistently exhibit peak energies above 1 keV in different flux states. The multi-epoch X-ray and GeV data reveal spectral and flux variabilities in RGB J0710+591 and 1ES 1741+196, even on timescales of days to weeks. As anticipated, the one-zone SSC model adequately reproduces the SEDs of regular HBLs but encounters difficulties in explaining the hardest TeV emission. Hadronic models offer a reasonable fit to the hard TeV spectrum, though with the trade-off of requiring extreme jet powers. On the other hand, the lepto-hadronic scenario faces additional challenges in fitting the GeV spectra of extreme-TeV sources. Finally, the
e
–
p
co-acceleration scenario naturally accounts for the observed hard electron distributions and effectively matches the hardest TeV spectrum of RGB J0710+591 and 1ES 1101−232.
Context. The one-zone synchrotron-self-Compton (SSC) model aims to describe the spectral energy distribution (SED) of BL Lac objects via synchrotron emission by a non-thermal population of electrons ...and positrons in a single homogeneous emission region, partially upscattered to γ-rays by the particles themselves. Aims. The model is usually considered as degenerate, given that the number of free parameters is higher than the number of observables. It is thus common to model the SED by choosing a single set of values for the SSC model parameters that provide a good description of the data, without studying the entire parameter space. We present here a new numerical algorithm that permits us to find the complete set of solutions, using the information coming from the detection in the GeV and TeV energy bands. Methods. The algorithm is composed of three separate steps: we first prepared a grid of simulated SEDs and extracted from each SED the values of the observables; we then parametrized each observable as a function of the SSC parameters; we finally solved the system for a given set of observables. We iteratively solved the system to take into account uncertainties in the values of the observables, producing a family of solutions. Results. We present a first application of our algorithm to the typical high-frequency-peaked BL Lac object 1RXS J101015.9 - 311909, provide constraints on the SSC parameters, and discuss the result in terms of our understanding of the blazar emitting region.
Context.
Sagittarius A*, the supermassive black hole at the center of our Galaxy, exhibits episodic near-infrared flares. The recent monitoring of three such events with the GRAVITY instrument has ...shown that some flares are associated with orbital motions in the close environment of the black hole. The GRAVITY data analysis indicates a super-Keplerian azimuthal velocity, while (sub-) Keplerian velocity is expected for the hot flow surrounding the black hole.
Aims.
We develop a semi-analytic model of the Sagittarius A* flares based on an ejected large plasmoid, inspired by recent particle-in-cell global simulations of black hole magnetospheres. We model the infrared astrometric and photometric signatures associated with this model.
Methods.
We considered a spherical macroscopic hot plasma region that we call a large plasmoid. This structure was ejected along a conical orbit in the vicinity of the black hole. This plasmoid was assumed to be formed by successive mergers of smaller plasmoids produced through magnetic reconnection that we did not model. Nonthermal electrons were injected into the plasmoid. We computed the evolution of the electron-distribution function under the influence of synchrotron cooling. We solved the radiative transfer problem associated with this scenario and transported the radiation along null geodesics of the Schwarzschild space time. We also took the quiescent radiation of the accretion flow into account, on top of which the flare evolves.
Results.
For the first time, we successfully account for the astrometric and flux variations of the GRAVITY data with a flare model that incorporates an explicit modeling of the emission mechanism. The prediction of our model and recent data agree well. In particular, the azimuthal velocity of the plasmoid is set by the magnetic field line to which it belongs, which is anchored in the inner parts of the accretion flow, hence the super-Keplerian motion. The astrometric track is also shifted with respect to the center of mass due to the quiescent radiation, in agreement with the difference measured with the GRAVITY data.
Conclusions.
These results support the hypothesis that magnetic reconnection in a black hole magnetosphere is a viable model for the infrared flares of Sagittarius A*.
Spreading of conservative solutes in groundwater due to aquifer heterogeneity is quantified by the macrodispersivity, which was found to be scale dependent. It increases with travel distance, ...stabilizing eventually at a constant value. However, the question of its asymptotic behavior at very large scale is still a matter of debate. It was surmised in the literature that macrodispersivity scales up following a unique scaling law. Attempts to define such a law were made by fitting a regression line in the log‐log representation of an ensemble of macrodispersivities from multiple experiments. The functional relationships differ among the authors, based on the choice of data. Our study revisits the data basis, used for inferring unique scaling, through a detailed analysis of literature marcodispersivities. In addition, values were collected from the most recent tracer tests reported in the literature. We specified a system of criteria for reliability and reevaluated the reliability of the reported values. The final collection of reliable estimates of macrodispersivity does not support a unique scaling law relationship. On the contrary, our results indicate, that the field data can be explained as a collection of macrodispersivities of aquifers with varying degree of heterogeneity where each exhibits its own constant asymptotic value. Our investigation concludes that transport, and particularly the macrodispersivity, is formation‐specific, and that modeling of transport cannot be relegated to a unique scaling law. Instead, transport requires characterization of aquifer properties, e.g., spatial distribution of hydraulic conductivity, and the use of adequate models.
Key Points:
Collection of macrodispersivity from literature and evaluation of reliability
Demonstration of inapplicability of an unique scaling law for macrodispersivity
Emphasis on the need for aquifer characterization for transport modeling