Cut-off rigidity model

Cut-Off Rigidity model (COR model) is a model for simulation of cosmic rays propagation in the Earth's magnetosphere. The model provides an evaluation of cut-off rigidities for stations at Earth surface or detectors measuring on orbits inside the magnetosphere.

Cut-off rigidities are evaluated for multiple incoming directions. A registered user could run simulations to evaluate cut-off rigidities for 577 incoming directions including the vertical one. To use the model, register here.

Complete simulation for one selected point with a 0.01GV step takes several hours. After submitting a request for simulation, your request will be processed and published when it will be done.

The catalogue of all to date realized simulations is in the part Simulations.

COR model uses a combined internal (IGRF) and external (Tsyganenko 96, Tsyganenko 05) geomagnetic fields. COR with Tsyganenko 96 model is available for years 1968 - 2016. COR with Tsyganenko 05 model is available for years 1995 - 2016. For more information about the model, results visit the section Publications.

COR model is also available for period 0 - 1968. Used geomagnetic models are described in article [Kudela K., Bobik P., 2004].

For educational purpose COR provides visualization of single cosmic ray trajectory in section New Trajectory.

For a collection of significant results, see this link. For updates, see this link.

COR model code

Code to solve Lorentz force equation by Runge-Kutta scheme of 6th order in Tsyganenko 04 model for simulating Earth's external magnetic field and IGRF (gen. 9 - 13) for simulating Earth's internal magnetic field is publicly available at Github repository Trajectories_IGRF_T04_C.

Example of model results

Effective cutoff rigidity sky map for Athens neutron monitor station with latitude 37.97 deg. N and longitude 23.78 deg E. The cutoff rigidity sky map was evaluated for 7. June 1998, with the Tsyganenko 05 model of the external geomagnetic field.

Example of single trajectory visualization

Example of the trajectory of a particle with energy 2.859 GeV (rigidity 3.68 GV) coming to a position with latitude 50 and longitude 0 coming from a direction from the zenith angle 10 and the azimuthal angle 70. The used backtracking model is described in [Bobik P et. al., 2005].

Authors and their affiliation at the time of contribution:

COR model:
Pavol Bobík, Ján Villim, Martin Vaško
Institute of Experimental Physics, Slovak Academy of Sciences; Technical University of Košice

Simulations management system and web application:
Daniel Gecášek, Michal Vrábel and Ján Genči
Technical University of Košice

Continoous integration and testing:
Dominik Leskoviansky, Pavol Kolcun
Technical University of Košice

Jupyter integration:
Jakub Pavelica, Anastasiia Zaytseva
Technical University of Košice

Visualisation module:
Martin Vaško and Ivan Bernat
Technical University of Košice