Search

Extensive air showers produce high-energy muons that can be utilized to probe hadronic interaction models in cosmic ray interactions. Most muons originate from pion and kaon decays, called $\textit{conventional}$ muons, while a smaller fraction, referred to as $\textit{prompt}$ muons, arises from the decay of heavier, short-lived hadrons. The $\texttt{EHISTORY}$ option of the air shower simulation tool $\texttt{CORSIKA7}$ is used in this work to investigate the prompt and conventional muon flux in the energy range of 100 TeV to 100 PeV, utilizing the newly developed open-source python software $\texttt{PANAMA}$. Identifying the muon parent particles allows for scaling the contribution of prompt particles, which can be leveraged by future experimental analyses to measure the normalization of the prompt muon flux. Obtained prompt muon spectra from $\texttt{CORSIKA7}$ are compared to $\texttt{MCEq}$ results. The relevance to large-volume neutrino detectors, such as IceCube and KM3NeT, and the connection to hadronic interaction models is discussed., Comment: 19 pages, 13 figures, software freely available at https://github.com/The-Ludwig/PANAMA

The sources of cosmic rays between the knee and the ankle are still debated. The Galactic wind and its termination shock have been proposed to contribute to this transition between Galactic and extragalactic origin, but another possibility is large-scale shock structures from local sources in the Milky Way. In this paper, we investigate CR transport in a time-dependent landscape of shocks in the Galactic halo. These shocks could result from local outbursts, e.g. starforming regions and superbubbles. CRs re-accelerated at such shocks can reach energies above the knee. Since the shocks are closer to the Galaxy than a termination shock and CRs escape downstream, they can propagate back more easily. With such outbursts happening frequently, shocks will interact. This interaction could adjust the CR spectrum, particularly for the particles that are able to be accelerated at two shocks simultaneously. The transport and acceleration of CRs at the shock is modeled by Stochastic Differential Equations (SDEs) within the public CR propagation framework CRPropa. We developed extensions for time-dependent wind profiles and for the first time connected the code to hydrodynamic simulations, which were run with the public Athena++ code. We find that, depending on the concrete realization of the diffusion tensor, a significant fraction of CRs can make it back to the Galaxy. These could contribute to the observed spectrum around and above the CR knee ($E \gtrsim 10\,\mathrm{PeV}$). In contrast to simplified models, a simple power-law does not describe the energy spectra well. Instead, for single shocks, we find a flat spectrum ($E^{-2}$) at low energies, which steepens gradually until it reaches an exponential decline. When shocks collide, the energy spectra transiently become harder than $E^{-2}$ at high energies.