Implementation of CR Energy SPectrum (CRESP) algorithm in PIERNIK MHD code. II. Propagation of Primary and Secondary nuclei in a magneto-hydrodynamical environment
Fuente:
arXiv
Saved in:
| Main Authors: | , , , , , , |
|---|---|
| Format: | Preprint |
| Published: |
2025
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866917336839094272 |
|---|---|
| author | Baldacchino-Jordan, Antoine Hanasz, Michał Ogrodnik, Mateusz Wóltański, Dominik Gawryszczak, Artur Strong, Andrew W. Girichidis, Philipp |
| author_facet | Baldacchino-Jordan, Antoine Hanasz, Michał Ogrodnik, Mateusz Wóltański, Dominik Gawryszczak, Artur Strong, Andrew W. Girichidis, Philipp |
| contents | We developed a new model for the production and propagation of spectrally resolved primary and secondary Cosmic Ray (CR) nuclei elements within the framework of the Cosmic Ray Energy Spectrum (CRESP) module of the PIERNIK MHD code. We extend the algorithm to several CR nuclei and demonstrate our code's capability to model primary and secondary CR species simultaneously. Primary C, N, and O are accelerated in supernova (SN) remnants. The spallation collisions of the primary nuclei against the thermal ISM protons lead to secondary Li, Be, and B products. All the CR species evolve according to the momentum-dependent Fokker-Planck equations that are dynamically coupled to the MHD system of equations governing the evolution of the ISM. We demonstrate the operation of this system in the gravity stratified box reproducing the Milky Way conditions in the Sun's local environment. We perform a parameter study by investigating the impact of the SN rate, the CR parallel diffusion coefficient $D_\parallel$, and the rigidity-dependent diffusion coefficient power index $δ$. A novel result of our investigation is that the secondary-to-primary flux ratio \BtoC increases with increasing diffusion coefficient, due to the weaker vertical magnetic field resulting from CR buoyancy effects. Moreover, a higher SN rate leads to lower values of \BtoC because of stronger winds and the shorter residence time of primary CR particles in dense disk regions. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2505_15614 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Implementation of CR Energy SPectrum (CRESP) algorithm in PIERNIK MHD code. II. Propagation of Primary and Secondary nuclei in a magneto-hydrodynamical environment Baldacchino-Jordan, Antoine Hanasz, Michał Ogrodnik, Mateusz Wóltański, Dominik Gawryszczak, Artur Strong, Andrew W. Girichidis, Philipp High Energy Astrophysical Phenomena We developed a new model for the production and propagation of spectrally resolved primary and secondary Cosmic Ray (CR) nuclei elements within the framework of the Cosmic Ray Energy Spectrum (CRESP) module of the PIERNIK MHD code. We extend the algorithm to several CR nuclei and demonstrate our code's capability to model primary and secondary CR species simultaneously. Primary C, N, and O are accelerated in supernova (SN) remnants. The spallation collisions of the primary nuclei against the thermal ISM protons lead to secondary Li, Be, and B products. All the CR species evolve according to the momentum-dependent Fokker-Planck equations that are dynamically coupled to the MHD system of equations governing the evolution of the ISM. We demonstrate the operation of this system in the gravity stratified box reproducing the Milky Way conditions in the Sun's local environment. We perform a parameter study by investigating the impact of the SN rate, the CR parallel diffusion coefficient $D_\parallel$, and the rigidity-dependent diffusion coefficient power index $δ$. A novel result of our investigation is that the secondary-to-primary flux ratio \BtoC increases with increasing diffusion coefficient, due to the weaker vertical magnetic field resulting from CR buoyancy effects. Moreover, a higher SN rate leads to lower values of \BtoC because of stronger winds and the shorter residence time of primary CR particles in dense disk regions. |
| title | Implementation of CR Energy SPectrum (CRESP) algorithm in PIERNIK MHD code. II. Propagation of Primary and Secondary nuclei in a magneto-hydrodynamical environment |
| topic | High Energy Astrophysical Phenomena |
| url | https://arxiv.org/abs/2505.15614 |