Source author record

Julia Becker Tjus

Julia Becker Tjus appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2025arXiv

A Seyfert galaxy as a hidden counterpart to a neutrino-associated blazar

The origin and production mechanisms of high-energy astrophysical neutrinos remain open questions in multimessenger astronomy. Previous studies have hinted at a possible linear correlation between the hard X-ray and high-energy neutrino emission in active galactic nuclei. New \textit{NuSTAR} observations, first presented here, reveal that blazar PKS 1424+240, located within a prominent IceCube neutrino hotspot, is far fainter in hard X-rays than expected from this trend. Motivated by this apparent ambiguity, we identify the nearby Seyfert galaxy NGC 5610, also coincident with the hotspot, whose unabsorbed hard X-ray flux exceeds that of PKS 1424+240 by about an order of magnitude. When the local IceCube neutrino flux is apportioned between the two AGN in proportion to their hard X-ray emission, both align with the previously suggested X-ray-neutrino correlation. This suggests that certain IceCube hotspots may be unresolved blends of multiple AGN, and supports a multimessenger scenario in which high-energy neutrinos and hard X-rays originate from the same hadronic interactions, with the X-ray emission produced through cascade reprocessing.

preprint2022arXiv

CRPropa 3.2 -- an advanced framework for high-energy particle propagation in extragalactic and galactic spaces

The landscape of high- and ultra-high-energy astrophysics has changed in the last decade, largely due to the inflow of data collected by large-scale cosmic-ray, gamma-ray, and neutrino observatories. At the dawn of the multimessenger era, the interpretation of these observations within a consistent framework is important to elucidate the open questions in this field. CRPropa 3.2 is a Monte Carlo code for simulating the propagation of high-energy particles in the Universe. This version represents a major leap forward, significantly expanding the simulation framework and opening up the possibility for many more astrophysical applications. This includes, among others: efficient simulation of high-energy particles in diffusion-dominated domains, self-consistent and fast modelling of electromagnetic cascades with an extended set of channels for photon production, and studies of cosmic-ray diffusion tensors based on updated coherent and turbulent magnetic-field models. Furthermore, several technical updates and improvements are introduced with the new version, such as: enhanced interpolation, targeted emission of sources, and a new propagation algorithm (Boris push). The detailed description of all novel features is accompanied by a discussion and a selected number of example applications.

preprint2022arXiv

Multiwavelength search for the origin of IceCube's neutrinos

The origin of astrophysical high-energy neutrinos detected by the IceCube Neutrino Observatory remains a mystery to be solved. In this paper we search for neutrino source candidates within the $90$% containment area of $70$ track-type neutrino events recorded by the IceCube Neutrino Observatory. By employing the Fermi-LAT 4FGL-DR2, the Swift-XRT 2SXPS and the CRATES catalogs, we identify possible gamma, X-ray and flat-spectrum radio candidate sources of track-type neutrinos. We find that based on the brightness of sources and their spatial correlation with the track-type IceCube neutrinos, the constructed neutrino samples represent special populations of sources taken from the full Fermi-LAT 4FGL-DR2/Swift-XRT 2SXPS/CRATES catalogs with similar significance ($2.1σ$, $1.2σ$, $2σ$ at $4.8~\mathrm{GHz}$, $2.1σ$ at $8.4~\mathrm{GHz}$, respectively, assuming 50% astrophysical signalness). After collecting redshifts and deriving sub-samples of the CRATES catalog complete in the redshift--luminosity plane, we find that the 4.8 GHz ($8.4$~GHz) sub-sample can explain between 4% and 53% ($3$% and $42$%) of the neutrinos (90% C.L.), when the probability to detect a neutrino is proportional to the ($k$-corrected) radio flux. The overfluctuations indicate that a part of the sample is likely to contribute and that more sophisticated schemes in the source catalog selection are necessary to identify the neutrino sources at the $5σ$ level. Our selection serves as a starting point to further select the correct sources.

preprint2020arXiv

A novel analytical model of the magnetic field configuration in the Galactic Center

Context: Cosmic-ray propagation is strongly dependent on the large-scale configuration of the Galactic magnetic field. In particular, the Galactic center (GC) region provides highly interesting cosmic-ray data from gamma-ray maps and it is clear that a large fraction of the cosmic rays detected at Earth originate in this region of the Galaxy. Yet because of confusion from line-of-sight integration, the magnetic field structure in the Galactic center is not well known and no large-scale magnetic field model exists at present. Aims: In this paper, we develop a magnetic field model, derived from observational data on the diffuse gas, nonthermal radio filaments, and molecular clouds. Methods. We derive an analytical description of the magnetic field structure in the central molecular zone by combining observational data with the theoretical modeling of the basic properties of magnetic fields. Results. We provide a first description of the large-scale magnetic field in the Galactic center region. We present first test simulations of cosmic-ray propagation and the impact of the magnetic field structure on the cosmic-ray distribution in the three dimensions. Conclusions: Our magnetic field model is able to describe the main features of polarization maps; it is particularly important to note that they are significantly better than standard global Galactic magnetic field models. It can also be used to model cosmic-ray propagation in the Galactic Center region more accurately.

preprint2020arXiv

Closing in on the origin of Galactic cosmic rays using multimessenger information

In cosmic ray physics extensive progress has been made in recent years, both concerning theory and observation. The vast details in direct, indirect and secondary detections on the one hand provide the basis for a detailed modeling of the signatures via cosmic-ray transport and interaction, paving the way for the identification of Galactic cosmic-ray sources. On the other hand, the large number of constraints from different channels of cosmic-ray observables challenges these models frequently. In this review, we will summarize the state-of-the art of the detection of cosmic rays and their secondaries, followed by a discussion what we can learn from coupling our knowledge of the cosmic-ray observables to the theory of cosmic-ray transport in the Galactic magnetic field. Finally, information from neutral secondaries will be added to draw a multimessenger-picture of the non-thermal sky, in which the hypothesis of supernova remnants as the dominant sources survives best. While this has been known since the 1930s, evidence for this scenario is steadily growing, with the first detection of hadronic signatures at GeV energies detected for three SNRs with Fermi. The existence of SNRs as PeVatrons, however, is not validated yet. The discussion of this and other open questions concerning the level of anisotropy, composition and spectral shape of the cosmic-ray energy spectrum is reviewed. Future perspectives of how to find the smoking cosmic-ray source gun concludes this review.

preprint2020arXiv

On the relative importance of hadronic emission processes along the jet axis of Active Galactic Nuclei

With the coincident detection of a gamma-ray flare and a neutrino from the blazar TXS 0506+056, Active Galactic Nuclei (AGN) have been put into focus as possible sources of the diffuse neutrino flux. We present a space and time-resolved model of the high-energy particle emission of a plasmoid assumed to travel along the axis of an AGN jet at relativistic speed. This was achieved by modifying the publicly available CRPropa (version 3.1+) propagation framework which in our work is capable of being applied to source physics on sub-kpc scales. The propagation of a population of primary protons is modelled in a purely turbulent magnetic field and we take into account interactions of these protons with photons scattered from the accretion disc, synchrotron radiation emitted by ambient relativistic electrons, as well with themselves and with other ambient matter. Our model produces a PeV-neutrino flare caused mainly by photo-hadronic interactions of primaries with the accretion disc field. Secondary high-energy gamma-rays partly attenuate with the ambient photon fields whose combined optical depths achieve their minimal opacity for photons of around 10 TeV. Thus, our model is well capable of producing neutrino flares with a significantly reduced emission of gamma-rays in jets with a hadronic jet component which in the future can be fit to specific AGN flare scenarios.

preprint2016arXiv

The Radio-Gamma Correlation In Starburst Galaxies

We present a systematic study of the non-thermal electron-proton plasma and its emission processes in starburst galaxies in order to explain the correlation between the luminosity in the radio band and the recently observed gamma luminosity. In doing so, a steady state description of the cosmic ray electrons and protons within the spatially homogeneous starburst is considered where continuous momentum losses are included as well as catastrophic losses due to diffusion and advection. The primary source of the relativistic cosmic rays, e.g. supernova remnants, provides a quasi-neutral plasma with a power law spectrum in momentum where we account for rigidity dependent differences between the electron and proton spectrum. We examine the resulting leptonic and hadronic radiation processes by synchrotron radiation, inverse Compton scattering, Bremsstrahlung and hadronic pion production. Finally, the observations of NGC 253, M82, NGC 4945 and NGC 1068 in the radio and gamma-ray band as well as the observed supernova rate are used to constrain a best-fit model. In the case of NGC 253, M82, NGC 4945 our model is able to accurately describe the data, showing that: (i) Supernovae are the dominant particle accelerators for NGC 253, M82 and NGC 4945, but not in the case of NGC 1068. (ii) All considered starburst galaxies are poor proton calorimeters in which for NGC 253 the escape is predominantly driven by the galactic wind, whereas the diffusive escape dominates in NGC 4945 and M82 (at energies > 1TeV). (iii) Secondary electrons from hadronic pion production are important to model the radio flux, but the associated neutrino flux is below the current observation limit.

preprint2015arXiv

Galactic Propagation of Cosmic Rays from Individual Supernova Remnants

It is widely believed that supernova remnants are the best candidate sources for the observed cosmic ray flux up to the knee, i.e. up to ~PeV energies. Indeed, the gamma-ray spectra of some supernova remnants can be well explained by assuming the decay of neutral pions which are created in hadronic interactions. Therefore, fitting the corresponding gamma spectra allows us to derive the spectra of cosmic rays at the source which are locally injected into our Galaxy. Using these spectra as a starting point, we propagate the cosmic rays through the Galaxy using the publicly available GALPROP code. Here, we will present first results on the contribution of those SNRs to the total cosmic ray flux and discuss implications.

preprint2015arXiv

Gamma-ray emitting supernova remnants as the origin of Galactic cosmic rays?

The origin of cosmic rays is one of the long-standing mysteries in physics and astrophysics. Simple arguments suggest that a scenario of supernova remnants (SNRs) in the Milky Way as the dominant sources for the cosmic ray population below the knee could work: in a generic calculation, it can be shown that these objects can provide the energy budget necessary to explain the observed flux of cosmic rays. However, this argument is based on the assumption that all sources behave in the same way, i.e.\ they all have the same energy budget, spectral behavior and maximum energy. In this paper, we investigate if a realistic population of SNRs is capable of producing the cosmic ray flux as it is observed below the knee. We use 21 SNRs that are well-studied from radio wavelengths up to gamma-ray energies. It could be shown previously (Mandelartz & Becker Tjus 2015) that the high-energy bump in the energy spectrum of these 21 sources can be dominated by hadronic emission. Here, gamma-rays are produced via $π^{0}-$decays from cosmic ray interactions in molecular clouds near the supernova remnant, which serves as the cosmic ray accelerator. The cosmic ray spectra show a large variety in their energy budget, spectral behavior and maximum energy. These sources are assumed to be representative for the total class of SNRs, where we assume that about 100 - 200 cosmic ray emitting SNRs should be present today. Finally, we use these source spectra to simulate the cosmic ray transport from individual SNRs in the Galaxy with the GALPROP code for cosmic ray propagation. We find that the cosmic ray budget can be matched well for a diffusion coefficient that is close to $D\propto E^{0.3}$. A stronger dependence on the energy, e.g. $E^{0.5}$, would lead to a spectrum at Earth that is too steep when compared to what is detected and the energy budget cannot be matched, in particular toward high energies.

preprint2014arXiv

Cosmic backgrounds due to the formation of the first generation of supermassive black holes

The statistics of black holes and their masses strongly suggests that their mass distribution has a cutoff towards lower masses near $3 \times 10^{6}$ M$_{\odot}$. This is consistent with a classical formation mechanism from the agglomeration of the first massive stars in the universe. However, when the masses of the stars approach $10^{6}$ M$_{\odot}$, the stars become unstable and collapse, possibly forming the first generation of cosmological black holes. Here we speculate that the claimed detection of an isotropic radio background may constitute evidence of the formation of these first supermassive black holes, since their data are compatible in spectrum and intensity with synchrotron emission from the remnants. The model proposed fulfills all observational conditions for the background, in terms of single-source strength, number of sources, far-infrared and gamma-ray emission. The observed high energy neutrino flux is consistent with our calculations in flux and spectrum. The proposal described in this paper may also explain the early formation and growth of massive bulge-less disk galaxies as derived from the massive, gaseous shell formed during the explosion prior to the formation of a supermassive black hole.

preprint2014arXiv

Impact of Secondary Acceleration in Gamma-Ray Bursts

We discuss the acceleration of secondary muons, pions, and kaons in gamma-ray bursts within the internal shock scenario, and their impact on the neutrino fluxes. We introduce a two-zone model consisting of an acceleration zone (the shocks) and a radiation zone (the plasma downstream the shocks). The acceleration in the shocks, which is an unavoidable consequence of the efficient proton acceleration, requires efficient transport from the radiation back to the acceleration zone. On the other hand, stochastic acceleration in the radiation zone can enhance the secondary spectra of muons and kaons significantly if there is a sufficiently large turbulent region. Overall, it is plausible that neutrino spectra can be enhanced by up to a factor of two at the peak by stochastic acceleration, that an additional spectral peaks appears from shock acceleration of the secondary muons and pions, and that the neutrino production from kaon decays is enhanced. Depending on the GRB parameters, the general conclusions concerning the limits to the internal shock scenario obtained by recent IceCube and ANTARES analyses may be affected by up to a factor of two by secondary acceleration. Most of the changes occur at energies above 10^7 GeV, so the effects for next-generation radio-detection experiments will be more pronounced. In the future, however, if GRBs are detected as high-energy neutrino sources, the detection of one or several pronounced peaks around 10^6 GeV or higher energies could help to derive the basic properties of the magnetic field strength in the GRB.

preprint2014arXiv

Ionisation as indicator for cosmic ray acceleration

Astrospheres and wind bubbles of massive stars are believed to be sources of cosmic rays with energies $E\lesssim 1\,$TeV. These particles are not directly detectable, but their impact on surrounding matter, in particular ionisation of atomic and molecular hydrogen, can lead to observable signatures. A correlation study of both gamma ray emission, induced by proton-proton interactions of cosmic ray protons with kinetic energies $E_\mathrm{p}\ge 280\,$MeV with ambient hydrogen, and ionisation induced by cosmic ray protons of kinetic energies $E_\mathrm{p}< 280\,$MeV can be performed in order to study potential sources of (sub)TeV cosmic rays.

preprint2014arXiv

Prediction of the diffuse neutrino flux from cosmic ray interactions near Supernova Remnants

In this paper, we present high-energy neutrino spectra from 21 Galactic supernova remnants (SNRs), derived from gamma-ray measurements in the GeV-TeV range. We find that only the strongest sources, i.e. G40.5-0.5 in the north and Vela Junior in the south could be detected as single point sources by IceCube or KM3NeT, respectively. For the first time, it is also possible to derive a diffuse signal by applying the observed correlation between gamma-ray emission and radio signal. Radio data from 234 supernova remnants listed in Green's catalog are used to show that the total diffuse neutrino flux is approximately a factor of 2.5 higher compared to the sources that are resolved so far. We show that the signal at above 10 TeV energies can actually become comparable to the diffuse neutrino flux component from interactions in the interstellar medium. Recently, the IceCube collaboration announced the detection of a first diffuse signal of astrophysical high-energy neutrinos. Directional information cannot unambiguously reveal the nature of the sources at this point due to low statistics. A number of events come from close to the Galactic center and one of the main questions is whether at least a part of the signal can be of Galactic nature. In this paper, we show that the diffuse flux from well-resolved SNRs is at least a factor of 20 below the observed flux.

preprint2014arXiv

Theoretical study of ionization profiles of molecular clouds near supernova remnants: Tracing the hadronic origin of GeV gamma radiation

Context: Since a few years, signatures of supernova remnants associated with molecular clouds have been detected in gamma rays. Whether these gamma rays are generated by cosmic ray electrons or by cosmic ray protons is usually not known. The detection of hadronic ionization signatures in spatial coincidence with gamma ray signatures can help to unambiguously identify supernova remnants as sources of cosmic ray protons. Methods: In order to calculate hadronic signatures from cosmic ray-induced ionization for an examination of the origin of the observed gamma rays, the transport equation for cosmic ray protons propagating in a molecular cloud, including the relevant momentum loss processes, is solved analytically and the proton flux at any position in the cloud is determined. Results: Since the solution of the transport equation is obtained for arbitrary source functions, it can be used for a variety of supernova remnants. The corresponding theoretical ionization rate, as a function of the penetration depth, is derived and compared to photoinduced ionization profiles in a case study with four supernova remnants associated with molecular clouds. Three of the remnants show a clear dominance of the hadronically induced ionization rate, while for one remnant, X-ray emission seems to dominate by a factor of 10. Conclusions: This is the first derivation of position-dependent profiles for cosmic ray-induced ionization with an analytic solution for arbitrary cosmic ray source spectra. The cosmic ray-induced ionization has to be compared to X-ray ionization for strong X-ray sources. For sources dominated by cosmic ray-induced ionization (e.g., W49B), the ionization profiles can be used in the future to map the spatial structure of hadronic gamma rays and rotation-vibrational lines induced by cosmic ray protons, helping to identify sources of hadronic cosmic rays.

preprint2013arXiv

Influence of hadronic interaction models and the cosmic ray spectrum on the high energy atmospheric muon and neutrino flux

The recent observations of muon charge ratio up to about 10 TeV and of atmospheric neutrinos up to energies of about 400 TeV has triggered a renewed interest into the high-energy interaction models and cosmic ray primary composition. A reviewed calculation of lepton spectra produced in cosmic-ray induced extensive air showers is carried out with a primary cosmic-ray spectrum that fits the latest direct measurements below the knee. In order to achieve this, we used a full Monte Carlo method to derive the inclusive differential spectra (yields) of muons, muon neutrinos and electron neutrinos at the surface for energies between 80 GeV and hundreds of PeV. The air shower simulator {\sc corsika} 6.990 was used for showering and propagation of the secondary particles through the atmosphere, employing the established high-energy hadronic interaction models {\sc sibyll} 2.1, {\sc qgsjet-01} and {\sc qgsjet-ii 03}. We show that the performance of the interaction models allows makes it possible to predict the spectra within experimental uncertainties, while {\sc sibyll} generally yields a higher flux at the surface than the qgsjet models. The calculation of the flavor and charge ratios has lead to inconsistent results, mainly influenced by the different representations of the K/$π$ ratio within the models. Furthermore, we could quantify systematic uncertainties of atmospheric muon- and neutrino fluxes, associated to the models of the primary cosmic-ray spectrum and the interaction models. For most recent parametrizations of the cosmic-ray primary spectrum, atmospheric muons can be determined with an uncertainty smaller than $^{+15}_{-13}$% of the average flux. Uncertainties of the muon- and electron neutrino fluxes can be calculated within an average error of $^{+32}_{-22}$% and $^{+25}_{-19}$%, respectively.

preprint2013arXiv

Limits on the source properties of FR-I galaxies from high-energy neutrino and gamma observations

Active galactic nuclei (AGN) are believed to be the source of ultra high energy cosmic rays (UHECRs). Particles are assumed to be accelerated in the accretion disk and the plasma jets, produced due to conservation of angular momentum, to the highest energies, where they interact with each other and produce pions, which decay among others in neutrinos. For a known cosmic ray spectral behavior, the main parameters in the calcula- tion of the neutrino flux from proton-proton interactions are the target density nH and the ratio of electrons to protons fe . Using most recent neutrino flux limits from IceCube point source searches, we set limits on the target densities for 33 FR-I galaxies. The densities are shown to be smaller than 30 cm^{-3} to 2000 cm^{-3}, depending on the source and when using a fixed electron to proton ratio of f_{e} = 0.1. This implies that some cosmic ray acceleration sites, espe- cially those close to the core of the AGN, can already be excluded, or else that the ratio of electrons to protons deviates significantly from the commonly used value of 0.1. For Centaurus A (Cen A) and Messier 87 (M 87) we use Fermi observations to model the γ-flux, the neutrino flux and the resulting target density. The detec- tion of these neutrinos will help to find information about acceleration processes in the source.

preprint2013arXiv

The Galactic Wind Haze and its $γ$-spectrum

We study the possibility that the gamma ray emission in the Fermi bubbles observed is produced by cosmic ray electrons with a spectrum similar to Galactic cosmic rays. We argue that the cosmic ray electrons steepen near 1 TeV from $E^{-3}$ to about $E^{-4.2}$, and are partially secondaries derived from the knee-feature of normal cosmic rays. We speculate that the observed feature at $\sim 130$ GeV could essentially be due to inverse Compton emission off a pair-production peak on top of a turn-off in the $γ$ ray spectrum at $\sim 130$ GeV. It suggests that the knee of normal cosmic rays is the same everywhere in the Galaxy. A consequence could be that all supernovae contributing give the same cosmic ray spectrum, with the knee feature given by common stellar properties; in fact, this is consistent with the supernova theory proposed by Bisnovatyi-Kogan (1970), a magneto-rotational mechanism, if massive stars converge to common properties in terms of rotation and magnetic fields just before they explode.