Cosmic Ray Production Research Articles

SOFIA/HAWC+ Far-infrared Polarimetric Large-area CMZ Exploration Survey. IV. Relative Magnetic Field Orientation throughout the CMZ

Abstract The nature of the magnetic field structure throughout the Galactic Center (GC) has long been of interest. The recent Far-InfraREd Polarimetric Large-Area Central Molecular Zone (CMZ) Exploration (FIREPLACE) Survey reveals preliminary connections between the seemingly distinct vertical and horizontal magnetic field distributions previously observed in the GC. We use the statistical techniques of the Histogram of Relative Orientation and the Projected Rayleigh Statistic to assess whether the CMZ magnetic field preferentially aligns with the structure of the CMZ molecular clouds or the morphology of the nonthermal emission of the GC nonthermal filament (NTF) population. We find that there is a range of magnetic field orientations throughout the population of CMZ molecular clouds, ranging from parallel to perpendicular orientation. We posit these orientations depend on the prevalence of gravitational shear in the GC, in contrast with what is observed in Galactic Disk star-forming regions. We also compare the magnetic field orientation from dust polarimetry with individual prominent NTFs, finding a preferred perpendicular relative orientation. This perpendicular orientation indicates that the vertical field component found in the FIREPLACE observations is not spatially confined to the NTFs, providing evidence for a more pervasive vertical field in the GC. From dynamical arguments, we estimate an upper limit on the magnetic field strength for this vertical field, finding B ≤ 4 mG. A field close to this upper limit would indicate that the NTFs are not local enhancements of a weaker background field and that the locations of the NTFs depend on proximity to sites of cosmic-ray production.

Mass and wind luminosity of young Galactic open clusters in Gaia DR2

Context. Star clusters constitute a significant part of the stellar population in our Galaxy. The feedback processes they exert on the interstellar medium impact multiple physical processes from the chemical to the dynamical evolution of the Galaxy. In addition, young and massive stellar clusters might act as efficient particle accelerators and contribute to the production of cosmic rays. Aims. We aim at evaluating the wind luminosity driven by the young (< 30 Myr) Galactic open stellar clusters observed by the Gaia space mission. This is crucial for determining the energy channeled into accelerated particles. Methods. To do this, we developed a method relying on the number, magnitude, and line-of-sight extinction of the stars observed per cluster. Assuming that the stellar mass function follows a Kroupa mass distribution and accounting for the maximum stellar mass allowed by the age and mass of the parent cluster, we conservatively estimated the mass and wind luminosity of 387 local clusters within the second data release of Gaia. Results. We compared the results of our computation with recent estimates of young cluster masses. With respect to these, our sample is three times more abundant, particularly above a few thousand solar masses. This is of the utmost relevance for predicting the gamma-ray emission resulting from the interaction of accelerated particles. The cluster wind luminosity distribution we obtained extends up to 3 × 1038 erg s−1. This is a promising feature in terms of potential particle acceleration scenarios.

Floor of Cosmogenic Neutrino Fluxes above 1017 eV

The search for neutrinos with energies greater than 1017 eV is being actively pursued. Although normalization of the dominant neutrino flux is highly uncertain, a floor level is guaranteed by the interactions of extragalactic cosmic rays with milky Way gas. We estimate that this floor level gives an energy flux of E2ϕν≃10−13−0.5+0.5GeVcm−2sr−1s−1 at 1018 eV, where uncertainties arise from the modeling of the gas distribution and the experimental determination of the mass composition of ultra–high-energy cosmic rays on Earth. Based on a minimal model of cosmic-ray production to explain the mass-discriminated energy spectra observed on Earth above 5 × 1018 eV, we also present generic estimates of the neutrino fluxes expected from extragalactic production that generally exceed the aforementioned guaranteed floor. The prospects for detecting neutrinos above 1018 eV remain challenging, however, unless proton acceleration to the highest energies is at play in a subdominant population of cosmic-ray sources or new physical phenomena are at work.

Origin of 60Fe nuclei in cosmic rays: the contribution of local OB associations

ABSTRACT The presence of live 60Fe nuclei (lifetime of 3.8 Myr) in cosmic rays detected by the Advanced Composition Explorer/Cosmic Ray Isotope Spectrometer instrument suggests a nearby nucleosynthesis source. 60Fe is primarily produced in core-collapse supernovae, and we aim to clarify whether the detected 60Fe nuclei can be associated with a particular local supernova. We consider 25 OB associations and subgroups located within 1 kpc of the Solar system based on recent Gaia census. A model is developed that combines stellar population synthesis within these OB associations, cosmic ray acceleration within associated superbubbles, and cosmic ray transport to the Solar system. The most critical model parameter impacting 60Fe cosmic ray production is the explodability criterion, which determines if a massive star ends its life as a supernova. Our study points to the Scorpius–Centaurus (Sco–Cen) OB association as the most probable origin of the observed 60Fe nuclei, particularly suggesting they were accelerated in the Sco–Cen superbubble by a young supernova aged ≤500 kyr with a progenitor mass of approximately 13–20 M⊙. A less likely source is the supernova at the origin of the Geminga pulsar 342 kyr ago, if the progenitor originated in the Orion OB1 association. The contribution of local OB associations to the cosmic ray density of stable 56Fe is estimated to be around 20 per cent, with some sensitivity to cosmic ray acceleration efficiency and diffusion coefficient. These findings shed light on the origins of cosmic ray nuclei, connecting them to nucleosynthesis events within our local cosmic neighbourhood.

The role of supernovae inside AGN jets in UHECR acceleration

Context. Jets of active galactic nuclei are potential accelerators of ultra high-energy cosmic rays. Supernovae can occur inside these jets and contribute to cosmic ray acceleration, particularly of heavy nuclei, but that contribution has been hardly investigated so far. Aims. We carried out a first dedicated exploration of the role of supernovae inside extragalactic jets in the production of ultra high-energy cosmic rays. Methods. We characterized the energy budget of supernova-jet interactions, and the maximum possible energies of the particles accelerated in those events, likely dominated by heavy nuclei. This allowed us to assess whether these interactions can be potential acceleration sites of ultra high-energy cosmic rays, or at least of their seeds. For that, we estimated the cosmic ray luminosity for different galaxy types, and compared the injection rate of cosmic ray seeds into the jet with that due to galactic cosmic ray entrainment. Results. Since the supernova is fueled for a long time by the luminosity of the jet, the energy of a supernova-jet interaction can be several orders of magnitude greater than that of an isolated supernova. Thus, despite the low rate of supernovae expected to occur in the jet, they could still provide more seeds for accelerating ultra high-energy particles than cosmic ray entrainment from the host galaxy. Moreover, these interactions can create sufficiently efficient accelerators to be a source of cosmic rays with energies ≳10 EeV. Conclusions. Supernova-jet interactions can contribute significantly to the production of ultra high-energy cosmic rays, either directly by accelerating these particles themselves or indirectly by providing pre-accelerated seeds.

Evidence for Thermal X-Ray Emission from the Synchrotron-dominated Shocks in Tycho’s Supernova Remnant

Young supernova remnant (SNR) shocks are believed to be the main sites of galactic cosmic-ray production, showing X-ray synchrotron-dominated spectra in the vicinity of their shock. While a faint thermal signature left by the shocked interstellar medium (ISM) should also be found in the spectra, proofs for such an emission in Tycho’s SNR have been lacking. We perform an extended statistical analysis of the X-ray spectra of five regions behind the blast wave of Tycho’s SNR using Chandra archival data. We use Bayesian inference to perform extended parameter space exploration and sample the posterior distributions of a variety of models of interest. According to Bayes factors, spectra of all five regions of analysis are best described by composite three-component models taking nonthermal emission, ejecta emission, and shocked ISM emission into account. The shocked ISM stands out the most in the northern limb of the SNR. We find for the shocked ISM a mean electron temperature keV for all regions and a mean ionization timescale cm−3 s resulting in a mean ambient density cm−3 around the remnant. We performed an extended analysis of the northern limb and show that the measured synchrotron cutoff energy is not well constrained in the presence of a shocked ISM component. Such results cannot currently be further investigated by analyzing emission lines in the 0.5–1 keV range, because of the low Chandra spectral resolution in this band. We show with simulated spectra that Athena X-ray Integral Field Unit future performances will be crucial to address this point.

Prospects for ultra-high-energy particle acceleration at relativistic shocks

ABSTRACT We study the acceleration of charged particles by ultra-relativistic shocks using test-particle Monte Carlo simulations. Two field configurations are considered: (i) shocks with uniform upstream magnetic field in the plane of the shock, and (ii) shocks in which the upstream magnetic field has a cylindrical geometry. Particles are assumed to diffuse in angle due to frequent non-resonant scattering on small-scale fields. The steady-state distribution of particles’ Lorentz factors is shown to approximately satisfy dN/dγ ∝ γ−2.2 provided the particle motion is scattering dominated on at least one side of the shock. For scattering dominated transport, the acceleration rate scales as tacc ∝ t1/2, though recovers Bohm scaling tacc ∝ t if particles become magnetized on one side of the shock. For uniform field configurations, a limiting energy is reached when particles are magnetized on both sides of the shock. For the cylindrical field configuration, this limit does not apply, and particles of one sign of charge will experience a curvature drift that redirects particles upstream. For the non-resonant scattering model considered, these particles preferentially escape only when they reach the confinement limit determined by the finite system size, and the distribution approaches the escapeless limit dN/dγ ∝ γ−1. The cylindrical field configuration resembles that expected for jets launched by the Blandford & Znajek mechanism, the luminous jets of active galactic nuclei and gamma-ray bursts thus provide favourable sites for the production of ultra-high-energy cosmic rays.

Interpretations of the cosmic ray secondary-to-primary ratios measured by DAMPE

Precise measurements of the boron-to-carbon and boron-to-oxygen ratios by DAMPE show clear hardenings around $100$ GeV/n, which provide important implications on the production, propagation, and interaction of Galactic cosmic rays. In this work we investigate a number of models proposed in literature in light of the DAMPE findings. These models can roughly be classified into two classes, driven by propagation effects or by source ones. Among these models discussed, we find that the re-acceleration of cosmic rays, during their propagation, by random magnetohydrodynamic waves may not reproduce sufficient hardenings of B/C and B/O, and an additional spectral break of the diffusion coefficient is required. The other models can properly explain the hardenings of the ratios. However, depending on simplifications assumed, the models differ in their quality in reproducing the data in a wide energy range. The models with significant re-acceleration effect will under-predict low-energy antiprotons but over-predict low-energy positrons, and the models with secondary production at sources over-predict high-energy antiprotons. For all models high-energy positron excess exists.

Modeling CO Line Profiles in Shocks of W28 and IC 443

Molecular emission arising from the interactions of supernova remnant (SNR) shock waves and molecular clouds provide a tool for studying the dispersion and compression that might kick-start star formation as well as understanding cosmic-ray production. Purely rotational CO emission created by magnetohydrodynamic shock in the SNR–molecular cloud interaction is an effective shock tracer, particularly for slow-moving, continuous shocks into cold inner clumps of the molecular cloud. In this work, we present a new theoretical radiative transfer framework for predicting the line profile of CO with the Paris–Durham 1D shock model. We generated line profile predictions for CO emission produced by slow, magnetized C shocks into gas of density ∼104 cm−3 with shock speeds of 35 and 50 km s−1. The numerical framework to reproduce the CO line profile utilizes the large velocity gradient (LVG) approximation and the omission of optically thick plane-parallel slabs. With this framework, we generated predictions for various CO spectroscopic observations up to J = 16 in SNRs W28 and IC 443, obtained with SOFIA, IRAM-30 m, APEX, and KPNO. We found that CO line profile prediction offers constraints on the shock velocity and pre-shock density independent of the absolute line brightness and requires fewer CO lines than diagnostics using a rotational excitation diagram.

Cosmic ray production in superbubbles

ABSTRACT We compute the production of cosmic rays (CRs) in the dynamical superbubble (SB) produced by a cluster of massive stars. Stellar winds, supernova remnants, and turbulence are found to accelerate particles so efficiently that the non-linear feedback of the particles must be taken into account in order to ensure that the energy balance is not violated. High-energy particles do not scatter efficiently on the turbulence and escape quickly after each supernova explosion, which makes both their intensity inside the bubble and injection in the interstellar medium intermittent. On the other hand, the stochastic acceleration of low-energy particles hardens the spectra at GeV energies. Because CRs damp the turbulence cascade, this hardening is less pronounced when non-linearities are taken into account. Nevertheless, spectra with hard components extending up to 1–10 GeV and normalized to an energy density of 1–100 eV cm−3 are found to be typical signatures of CRs produced in SBs. Efficient shock reacceleration within compact clusters is further shown to produce hard, slightly concave spectra, while the presence of a magnetized shell is shown to enhance the confinement of CRs in the bubble and therefore the collective plasma effects acting on them. We eventually estimate the overall contribution of SBs to the Galactic CR content and show typical gamma-ray spectra expected from hadronic interactions in SB shells. In both cases, a qualitative agreement with observations is obtained.

Cosmic Ray production in SNRs of Pulsar Wind Nebulae type at different ages

The overall observations of plerions from radio to the very high energy gamma-rays could provide information about the evolution of PWN from the young Crab-like to the older stages. The extended MeV-TeV emission from Geminga middle-aged pulsar of ~ 3.4 × 105year age in SHALON, Milagro, HAWC observations, and Fermi-LAT detection could arise from the PWN associated with the Geminga SNR. The TeVγ-ray emission from 3C 58 PWN whose estimated age varies from ~ 800 to (5 − 7) × 103years was first detected in year 2011 by SHALON. The overall spectral energy distribution and information about the extension of PWN from radio to GeV-TeV energies from Fermi-LAT and SHALON observations can contribute to particle transport models and also to the understanding of the mechanisms of PWN expanding, which is, in turn, can shed light on the age of 3C 58 and the history of progenitor SN explosion.

A Simulation Study of Ultra-relativistic Jets. II. Structures and Dynamics of FR-II Jets

Abstract We study the structures of ultra-relativistic jets injected into the intracluster medium and the associated flow dynamics, such as shocks, velocity shear, and turbulence, through three-dimensional relativistic hydrodynamic (RHD) simulations. To that end, we have developed a high-order accurate RHD code, equipped with a weighted essentially non-oscillatory scheme and a realistic equation of state. Using the code, we explore a set of jet models with the parameters relevant to FR-II radio galaxies. We confirm that the overall jet morphology is primarily determined by the jet power, and the jet-to-background density and pressure ratios play secondary roles. Jets with higher powers propagate faster, resulting in more elongated structures, while those with lower powers produce more extended cocoons. Shear interfaces in the jet are dynamically unstable, and hence, chaotic structures with shocks and turbulence develop. We find that the fraction of the jet-injected energy dissipated through shocks and turbulence is greater in less powerful jets, although the actual amount of the dissipated energy is larger in more powerful jets. In lower power jets, the backflow is dominant in the energy dissipation owing to the broad cocoon filled with shocks and turbulence. In higher power jets, by contrast, both the backflow and jet-spine flow are important for the energy dissipation. Our results imply that different mechanisms, such as diffusive shock acceleration, shear acceleration, and stochastic turbulent acceleration, may be involved in the production of ultra-high energy cosmic rays in FR-II radio galaxies.

Millisecond pulsars modify the radio-star-formation-rate correlation in quiescent galaxies

The observed correlation between the far-infrared and radio luminosities of galaxies illustrates the close connection between star formation and cosmic-ray production. Intriguingly, recent gamma-ray observations indicate that recycled/millisecond pulsars (MSPs), which do not trace recent star formation, may also efficiently accelerate cosmic-ray electrons. We study the contribution of MSPs to the galactic non-thermal radio emission, finding that they can dominate the emission from massive quiescent galaxies. This model can explain recent LOFAR observations that found a peculiar radio excess in galaxies with high stellar masses and low star-formation rates. We show that MSP-based models provide a significantly improved fit to LOFAR data. We discuss the implications for the radio-FIR correlation, the observation of radio excesses in nearby galaxies, and local electron and positron observations.

Connecting turbulent velocities and magnetic fields in galaxy cluster simulations with active galactic nuclei jets

ABSTRACT The study of velocity fields of the hot gas in galaxy clusters can help to unravel details of microphysics on small scales and to decipher the nature of feedback by active galactic nuclei (AGN). Likewise, magnetic fields as traced by Faraday rotation measurements (RMs) inform about their impact on gas dynamics as well as on cosmic ray production and transport. We investigate the inherent relationship between large-scale gas kinematics and magnetic fields through non-radiative magnetohydrodynamical simulations of the creation, evolution, and disruption of AGN jet-inflated lobes in an isolated Perseus-like galaxy cluster, with and without pre-existing turbulence. In particular, we connect cluster velocity measurements with mock RM maps to highlight their underlying physical connection, which opens up the possibility of comparing turbulence levels in two different observables. For single-jet outbursts, we find only a local impact on the velocity field, i.e. the associated increase in velocity dispersion is not volume-filling. Furthermore, in a setup with pre-existing turbulence, this increase in velocity dispersion is largely hidden. We use mock X-ray observations to show that at arcmin resolution, the velocity dispersion is therefore dominated by existing large-scale turbulence and is only minimally altered by the presence of a jet. For the velocity structure of central gas uplifted by buoyantly rising lobes, we find fast, coherent outflows with low velocity dispersion. Our results highlight that projected velocity distributions show complex structures, which pose challenges for the interpretation of observations.

Nonthermal element abundances at astrophysical shocks

The nonthermal source abundances of elements play a crucial role in the understanding of cosmic ray phenomena from a few GeV up to several tens of EeV . We present a first systematic approach to describe the change of the abundances from the thermal to the nonthermal state via diffusive shock acceleration by a temporally evolving shock. We consider hereby not only ionization states of elements contained in the ambient gas, which we allow to be time dependent due to shock heating, but also elements condensed on solid, charged dust grains which can be injected into the acceleration process as well. Our generic parametrized model is then applied to the case of particle acceleration by supernova remnants in various ISM phases, for which we use state-of-the-art computation packages to calculate the ionization states of all elements. The resulting predictions for low energy cosmic ray (LECR) source abundances are compared with the data obtained by various experiments. We obtain excellent agreement for shocks in warm ionized ISM environments, which include HII regions, if dust grains are injected into the diffusive shock acceleration process with a much higher efficiency than ions. Less dependence of the fit quality is found on the mass-to-charge ratio of ions. For neutral environments, assuming that there are shocks in the weakly ionized component, and for the hot ionized medium we obtain generally inferior fits, but except for the cold neutral medium we do not exclude them as subdominant sites of Galactic cosmic ray production. The key challenge is found to be putting the LECR abundance of pure gas phase elements like neon and the (semi-)volatile elements phosphorus, sulfur and chlorine into the right balance with silicon, calcium and elements of the iron group. We present a brief outlook to the potential consequences of our results for the understanding of the composition around the second knee or the cosmic ray spectrum, or for the viability of explaining ultra-high energy cosmic rays with a dominant contribution by radio galaxies.

Modeling of Cosmic-Ray Production and Transport and Estimation of Gamma-Ray and Neutrino Emissions in Starburst Galaxies

Abstract Starburst galaxies (SBGs) with copious massive stars and supernova (SN) explosions are the sites of active cosmic-ray production. Based on the predictions of nonlinear diffusive shock acceleration theory, we model the cosmic-ray proton (CRP) production by both pre-SN stellar winds (SWs) and supernova remnants (SNRs) from core-collapse SNe inside the starburst nucleus. Adopting different models for the transport of CRPs, we estimate the γ-ray and neutrino emissions due to pp collisions from nearby SBGs such as M82, NGC 253, and Arp 220. We find that with the current γ-ray observations by Fermi-LAT, Veritas, and H.E.S.S., it would be difficult to constrain CRP production and transport models. Yet, the observations are better reproduced with (1) the combination of a single power-law (PL) momentum distribution for SNR-produced CRPs and a diffusion model in which the CRP diffusion is mediated by the strong Kolmogorov-type turbulence of δ B / B ∼ 1 , and (2) the combination of a double-PL model for SNR-produced CRPs and a diffusion model in which the scattering of CRPs is controlled mostly by self-excited waves rather than the preexisting turbulence. The contribution of SW-produced CRPs could be substantial in Arp 220, where the star formation rate is higher and the slope of the initial mass function would be flatter. We suggest that M82 and NGC 253 might be detectable as point sources of high-energy neutrinos in the upcoming KM3NET and IceCube-Gen2, when optimistic models are applied. Future observations of neutrinos, as well as γ-rays, would provide constraints for the production and diffusion of CRPs in SBGs.

Constraints on millicharged particles from cosmic-ray production

We study cosmic-ray-atmosphere collisions as a permanent production source of exotic millicharged particles (MCPs) for all terrestrial experiments. [MCPs are also known as charged massive particles (CHAMPs).] Based on data from Super-K, this allows us to derive new limits on MCPs that are competitive with, or improve, the currently leading bounds from accelerator-based searches for masses up to 1.5 GeV. In models where a subdominant component of dark matter (DM) is fractionally charged, these constraints probe parts of the parameter space that is inaccessible for conventional direct-detection DM experiments, independently of assumptions about the DM abundance.

Cosmogenic activation of silicon

The production of $^{3}$H, $^{7}$Be, and $^{22}$Na by interactions of cosmic-ray particles with silicon can produce radioactive backgrounds in detectors used to search for rare events. Through controlled irradiation of silicon CCDs and wafers with a neutron beam that mimics the cosmic-ray neutron spectrum, followed by direct counting, we determined that the production rate from cosmic-ray neutrons at sea level is ($112 \pm 24$) atoms/(kg day) for $^{3}$H, ($8.1 \pm 1.9 $) atoms/(kg day) for $^{7}$Be, and ($43.0 \pm 7.1 $) atoms/(kg day) for $^{22}$Na. Complementing these results with the current best estimates of activation cross sections for cosmic-ray particles other than neutrons, we obtain a total sea-level cosmic-ray production rate of ($124 \pm 24$) atoms/(kg day) for $^{3}$H, ($9.4 \pm 2.0 $) atoms/(kg day) for $^{7}$Be, and ($49.6 \pm 7.3 $) atoms/(kg day) for $^{22}$Na. These measurements will help constrain background estimates and determine the maximum time that silicon-based detectors can remain unshielded during detector fabrication before cosmogenic backgrounds impact the sensitivity of next-generation rare-event searches.

Cosmic-ray antinuclei as messengers of new physics: status and outlook for the new decade

The precise measurement of cosmic-ray antinuclei serves as an important means for identifying the nature of dark matter and other new astrophysical phenomena, and could be used with other cosmic-ray species to understand cosmic-ray production and propagation in the Galaxy. For instance, low-energy antideuterons would provide a “smoking gun” signature of dark matter annihilation or decay, essentially free of astrophysical background. Studies in recent years have emphasized that models for cosmic-ray antideuterons must be considered together with the abundant cosmic antiprotons and any potential observation of antihelium. Therefore, a second dedicated Antideuteron Workshop was organized at UCLA in March 2019, bringing together a community of theorists and experimentalists to review the status of current observations of cosmic-ray antinuclei, the theoretical work towards understanding these signatures, and the potential of upcoming measurements to illuminate ongoing controversies. This review aims to synthesize this recent work and present implications for the upcoming decade of antinuclei observations and searches. This includes discussion of a possible dark matter signature in the AMS-02 antiproton spectrum, the most recent limits from BESS Polar-II on the cosmic antideuteron flux, and reports of candidate antihelium events by AMS-02; recent collider and cosmic-ray measurements relevant for antinuclei production models; the state of cosmic-ray transport models in light of AMS-02 and Voyager data; and the prospects for upcoming experiments, such as GAPS. This provides a roadmap for progress on cosmic antinuclei signatures of dark matter in the coming years.

Supermassive Black Holes as Possible Sources of Ultrahigh-energy Cosmic Rays

Abstract The production and acceleration mechanisms of ultrahigh-energy cosmic rays (UHECRs) of energy >1020 eV, clearly beyond the GZK cutoff limit, remain unclear, which points to the exotic nature of the phenomena. Recent observations of extragalactic neutrinos may indicate that the source of UHECRs is an extragalactic supermassive black hole (SMBH). We demonstrate that ultraefficient energy extraction from a rotating SMBH driven by the magnetic Penrose process (MPP) could indeed fit the bill. We envision ionization of neutral particles, such as neutron beta decay, skirting close to the black hole horizon that energizes protons to over 1020 eV for an SMBH of mass 109 M ⊙ and magnetic field 104 G. Applied to the Galactic center SMBH, we have a proton energy of order ≈1015.6 eV that coincides with the knee of the cosmic-ray spectra. We show that large γ z factors of high-energy particles along the escaping directions occur only in the presence of an induced charge of the black hole, which is known as the Wald charge in the case of a uniform magnetic field. It is remarkable that the process requires neither an extended acceleration zone nor fine-tuning of accreting-matter parameters. Further, this leads to certain verifiable constraints on the SMBH’s mass and magnetic field strength as the source of UHECRs. This clearly makes the ultraefficient regime of the MPP one of the most promising mechanisms for fueling the UHECR powerhouse.