220 likes | 238 Views
GALPROP simulates gamma ray and cosmic ray skies with initial conditions for comparison with astronomical data. GaDGET optimizes GALPROP models using Fermi data, refining cosmic ray propagation processes for precise results. Study Milky Way parameters and extragalactic gamma ray spectrum.
E N D
Simulating the Gamma Ray Sky Andrew McLeod SASS August 12, 2009
Astrophysics Measurable Quantities • Particle Flux from a given region of the sky • Particle Energies (binned)
Astrophysics EGRET gamma ray sky, > 100 MeV “Modeling of the Galactic diffuse continuum gamma-ray emission”, Igor V. Moskalenko Source: http://galprop.stanford.edu/web_galprop/galprop_manual/manual2.html
GALPROP • Simulates the gamma ray and cosmic ray sky given a set of initial conditions and physical parameters • Allows a priori predictions to be compared to astronomical data
GALPROP Method Generates a field of cosmic rays given a cosmic ray source distribution (SNR, pulsars) by solving the transport equation:
GALPROP Method Computes the interaction of these cosmic rays with the interstellar medium and radiation field • Gas model based on 21-cm (atomic H) and CO (H2) surveys • 3-D or 2-D options; 2-D models have radial symmetry and model gas in 17 concentric rings
“Propagation of cosmic rays: nuclear physics in cosmic ray studies”, Igor V. Moskalenko Source: http://galprop.stanford.edu/web_galprop/galprop_manual/manual2.html GALPROP
GALPROP Method The contribution of these processes (and others) to the gamma ray spectrum are calculated as a function of radial distance
GALPROP Output • Cosmic ray spectra • Bremsstrahlung gamma ray spectrum • Neutral Pion Decay gamma ray spectrum • Inverse-Compton Scattering gamma ray spectrum All output in FITS format, as a function of radius (or gas ring) and energy
GALPROP Predicted gamma ray sky from Bremsstrahlung (gas ring 13)
GALPROP Predicted gamma ray sky from Neutral Pion Decay (gas ring 13)
GaDGET • Optimizes GALPROP models to best-fit observed gamma ray sky (Fermi data) • Fit-weights are computed for each component’s energy bins • Indicates how GALPROP models can be improved
GaDGET Method • Sums over gamma ray contribution of inner, local, and outer galaxy (relative to earth, ~8.5 kpc) • Assigns re-normalizing fit-weights to each region in incremental energy bins
GaDGET For resulting fit-weight adjusted sky-maps • Produces residual sky-maps • Calculates likelihood of model being correct
GaDGET Residual Map, normalized to 1σ
Optimization GALPROP parameters can be varied • Galactic Dimensions • Cosmic Ray Injection Spectra • Source Distribution • Diffusion Coefficient ~ 40 dimensional parameter space
Optimization Propagation processes can be altered • Reacceleration • Convection
Optimization With a full year of Fermi data, models can now be refined with an unprecedented level of precision
Potential Applications • Indirect determination of Milky Way parameters • Better understand the processes by which cosmic rays propagate • Study extragalactic gamma ray spectrum
Works Cited • Moskalenko, Igor. “Modeling of the Galactic diffuse continuum gamma-ray emission” 6th INTEGRAL Workshop, Moscow, Russia. 2006. • Moskalenko, Igor. “Propagation of Cosmic Rays and Diffuse Galactic Gamma Rays” Nuclear Data for Science and Technology, Santa Fe, New Mexico. 2004. • Strong, Andrew. “GALPROP: a Cosmic-ray propagation and Gamma-ray code” Tools for SUSY, Annecy, France. 2006.