Magnetospheric models and trajectory computations

The calculation of particle trajectories in the Earth's magnetic field has been a subject of interest since the time of Stormer. The fundamental problem is that the trajectory-tracing process involves using mathematical equations that have `no solution in closed form'. This difficulty has forced researchers to use the `brute force' technique of numerical integration of many individual trajectories to ascertain the behavior of trajectory families or groups. As the power of computers has improved over the decades, the numerical integration procedure has grown more tractable and while the problem is still formidable, thousands of trajectories can be computed without the expenditure of excessive resources. As particle trajectories are computed and the characteristics analyzed we can determine the cutoff rigidity of a specific location and viewing direction and direction and deduce the direction in space of various cosmic ray anisotropies. Unfortunately, cutoff rigidities are not simple parameters due to the chaotic behavior of the cosmic-ray trajectories in the cosmic ray penumbral region. As the computational problem becomes more manageable, there is still the issue of the accuracy of the magnetic field models. Over the decades, magnetic field models of increasing complexity have been developed and utilized. The accuracy of trajectory calculations employing contemporary magnetic field models is sufficient that cosmic ray experiments can be designed on the basis of trajectory calculations. However, the Earth's magnetosphere is dynamic and the most widely used magnetospheric models currently available are static. This means that the greatest uncertainly in the application of charged particle trajectories occurs at low energies.