GALPROP WebRun: an internet-based service for calculating galactic cosmic ray propagation and associated photon emissions

GALPROP is a numerical code for calculating the galactic propagation of relativistic charged particles and the diffuse emissions produced during their propagation. The code incorporates as much realistic astrophysical input as possible together with latest theoretical developments and has become a de facto standard in astrophysics of cosmic rays. We present GALPROP WebRun, a service to the scientific community enabling easy use of the freely available GALPROP code via web browsers. In addition, we introduce the latest GALPROP version 54, available through this service.Comment: Accepted for publication in Computer Physics Communications. Version 2 includes improvements suggested by the referee. Metadata completed in version 3 (no changes to the manuscript

Techniques for carrying out radiative transfer calculations for the Martian atmospheric dust

A description is given of the modification of a theory on the reflectance of particulate media so as to apply it to analysis of the infrared spectra obtained by the IRIS instrument on Mariner 9. With the aid of this theory and the optical constants of muscovite mica, quartz, andesite, anorthosite, diopside pyroxenite, and dunite, modeling calculations were made to refine previous estimates of the mineralogical composition of the Martian dust particles. These calculations suggest that a feldspar rich mixture is a very likely composition for the dust particles. The optical constants used for anorthosite and diopside pyroxenite were derived during this program from reflectance measurements. Those for the mica were derived from literature reflectance data. Finally, a computer program was written to invert the measured radiance data so as to obtain the absorption coefficient spectrum which should then be independent of the temperature profile and gaseous component effects

Challenging cosmic ray propagation with antiprotons. Evidence for a "fresh" nuclei component?

Recent measurements of the cosmic ray (CR) antiproton flux have been shown to challenge existing CR propagation models. It was shown that the reacceleration models designed to match secondary to primary nuclei ratios (e.g., boron/carbon) produce too few antiprotons. Matching both the secondary to primary nuclei ratio and the antiproton flux requires artificial breaks in the diffusion coefficient and the primary injection spectrum suggesting the need for other approaches. In the present paper we discuss one possibility to overcome these difficulties. Using the measured antiproton flux AND B/C ratio to fix the diffusion coefficient, we show that the spectra of primary nuclei as measured in the heliosphere may contain a fresh local "unprocessed" component at low energies perhaps associated with the Local Bubble, thus decreasing the measured secondary to primary nuclei ratio. The independent evidence for SN activity in the solar vicinity in the last few Myr supports this idea. The model reproduces antiprotons, B/C ratio, and elemental abundances up to Ni (Z<=28). Calculated isotopic distributions of Be and B are in perfect agreement with CR data. The abundances of three "radioactive clock" isotopes in CR, 10Be, 26Al, 36Cl, are all consistent and indicate a halo size z_h~4 kpc based on the most accurate data taken by the ACE spacecraft.Comment: To be published in The Astrophysical Journal, v.586, 2003 April 1; final version: 19 pages, 24 ps-figures, emulateapj5.sty (modified), natbib.sty, aastex.cls. More details can be found at http://www.gamma.mpe-garching.mpg.de/~aws/aws.htm

Propagation of Light Elements in the Galaxy

The origin and evolution of isotopes of the lightest elements H2, He3, Li, Be, B in the universe is a key problem in such fields as astrophysics of CR, Galactic evolution, non-thermal nucleosynthesis, and cosmological studies. One of the major sources of these species is spallation by CR nuclei in the interstellar medium. On the other hand, it is the B/C ratio in CR and Be10 abundance which are used to fix the propagation parameters and thus the spallation rate. We study the production and Galactic propagation of isotopes of elements Z<6 using the numerical propagation code GALPROP and updated production cross sections.Comment: 4 pages, 6 ps-figures, tsukuba.sty, to appear in the Proc. 28th International Cosmic Ray Conference (Tsukuba, Japan 2003). More details can be found at http://www.gamma.mpe-garching.mpg.de/~aws/aws.htm

Surface circulation in the Great Lakes as observed by LANDSAT-1 August 1972 - December 1973: Southern Lake Michigan

The surface current circulation patterns of southern Lake Michigan were charted for all cardinal and subcardinal wind directions, employing LANDSAT-1 observations of the distribution of natural tracing material borne in the surface waters. These colorants consist chiefly of river discharges composed of suspended sediments, pollutants, and algae; extensive chemical precipitations proved valuable for areas farther from shore. Comparison of the satellite-derived surface current charts with previous theoretical and empirical studies shows good agreement

Propagation of cosmic rays: nuclear physics in cosmic-ray studies

The nuclei fraction in cosmic rays (CR) far exceeds the fraction of other CR species, such as antiprotons, electrons, and positrons. Thus the majority of information obtained from CR studies is based on interpretation of isotopic abundances using CR propagation models where the nuclear data and isotopic production cross sections in p- and alpha-induced reactions are the key elements. This paper presents an introduction to the astrophysics of CR and diffuse gamma rays and discusses some of the puzzles that have emerged recently due to more precise data and improved propagation models. Merging with cosmology and particle physics, astrophysics of CR has become a very dynamic field with a large potential of breakthrough and discoveries in the near future. Exploiting the data collected by the CR experiments to the fullest requires accurate nuclear cross sections.Comment: 6 pages, 13 figures, aip style files. Invited review talk at the Int. Conf. on Nuclear Data-2004 (Santa Fe, Sep. 26 - Oct. 1, 2004). To appear in AIP Conf. Pro

Diffuse MeV Gamma-rays and Galactic 511 keV Line from Decaying WIMP Dark Matter

The origin of both the diffuse high-latitude MeV gamma-ray emission and the 511 keV line flux from the Galactic bulge are uncertain. Previous studies have invoked dark matter physics to independently explain these observations, though as yet none has been able to explain both of these emissions within the well-motivated framework of Weakly-Interacting Massive Particles (WIMPs). Here we use an unstable WIMP dark matter model to show that it is in fact possible to simultaneously reconcile both of these observations, and in the process show a remarkable coincidence: decaying dark matter with MeV mass splittings can explain both observations if positrons and photons are produced with similar branching fractions. We illustrate this idea with an unstable branon, which is a standard WIMP dark matter candidate appearing in brane world models with large extra dimensions. We show that because branons decay via three-body final states, they are additionally unconstrained by searches for Galactic MeV gamma-ray lines. As a result, such unstable long-lifetime dark matter particles provide novel and distinct signatures that can be tested by future observations of MeV gamma-rays.Comment: 19 pages, 4 figure

Diffractive Interaction and Scaling Violation in pp->pi^0 Interaction and GeV Excess in Galactic Diffuse Gamma-Ray Spectrum of EGRET

We present here a new calculation of the gamma-ray spectrum from pp->pi^0 in the Galactic ridge environment. The calculation includes the diffractive pp interaction and incorporates the Feynman scaling violation for the first time. Galactic diffuse gamma-rays come, predominantly, from pi^0->gamma gamma in the sub-GeV to multi-GeV range. Hunter et al. found, however, an excess in the GeV range ("GeV Excess") in the EGRET Galactic diffuse spectrum above the prediction based on experimental pp->pi^0 cross-sections and the Feynman scaling hypothesis. We show, in this work, that the diffractive process makes the gamma-ray spectrum harder than the incident proton spectrum by ~0.05 in power-law index, and, that the scaling violation produces 30-80% more pi^0 than the scaling model for incident proton energies above 100GeV. Combination of the two can explain about a half of the "GeV Excess" with the local cosmic proton (power-law index ~2.7). The excess can be fully explained if the proton spectral index in the Galactic ridge is a little harder (~0.2 in power-law index) than the local spectrum. Given also in the paper is that the diffractive process enhances e^+ over e^- and the scaling violation gives 50-100% higher p-bar yield than without the violation, both in the multi-GeV range.Comment: 35 pages, 11 figures, to appear in Astrophysical Journa

Source population synthesis and the Galactic diffuse gamma-ray emission

Population synthesis is used to study the contribution from undetected sources to the Galactic ridge emission measured by EGRET. Synthesized source counts are compared with the 3rd EGRET catalogue at low and high latitudes. For pulsar-like populations, 5-10% of the emission >100 MeV comes from sources below the EGRET threshold. A steeper luminosity function can increase this to 20% without violating EGRET source statistics. Less luminous populations can produce much higher values without being detected. Since the unresolved source spectrum is different from the interstellar spectrum, it could provide an explanation of the observed MeV and GeV excesses above the predictions, and we give an explicit example of how this could work.Comment: Astrophysics and Space Science, in press. (Proceedings of Conference 'The multi-messenger approach to high-energy gamma-ray sources', Barcelona, 2006). Minor changes for accepted version, updated reference