Hard X-ray astronomical observatories in orbit suffer from a significant amount of background due to radioactivation induced by cosmic-ray protons and/or geomagnetically trapped protons. Within the framework of a full Monte Carlo simulation, we present modeling of in-orbit instrumental background dominated by radioactivation. To reduce the computation time required by straightforward simulations of delayed emissions from activated isotopes, we insert a semi-analytical calculation that converts production probabilities of radioactive isotopes by interaction of the primary protons into decay rates at measurement time of all secondary isotopes. Therefore, our simulation method is separated into three steps: (1) simulation of isotope production, (2) semi-analytical conversion to decay rates, and (3) simulation of decays of the isotopes at measurement time. This method is verified by a simple setup that has a CdTe semiconductor detector, and shows a 100-fold improvement in efficiency over the straightforward simulation. The simulation framework was tested against data measured with a CdTe sensor in the Hard X-ray Imager onboard the Hitomi X-ray Astronomy Satellite, which was put into a low Earth orbit with an altitude of 570 km and an inclination of 31 degrees, and thus experienced a large amount of irradiation from geomagnetically trapped protons during its passages through the South Atlantic Anomaly. The simulation is able to treat full histories of the proton irradiation and multiple measurement windows. The simulation results agree very well with the measured data, showing that the measured background is well described by the combination of proton-induced radioactivation of the CdTe detector itself and thick Bi4Ge3O12 scintillator shields, leakage of cosmic X-ray background and albedo gamma-ray radiation, and emissions from naturally contaminated isotopes in the detector system.

We present the yield and spectra of stable secondary particles (gamma, e^+/-, nu_e, nubar_e, nu_mu, and nubar_mu) of p-p interaction in parameterized formulae to facilitate calculations involving them in astronomical environments. The formulae are derived on the up-to-date p-p interaction model by [Kamae05] which incorporates the logarithmically rising inelastic cross-section, the diffraction dissociation process, and the Feynman scaling violation. To improve fidelity to experimental data in lower energies, two baryon resonance contributions have been added: one representing Delta(1232) and the other multiple resonances around 1600MeV/c^2. The parametrized formulae predict that all secondary particle spectra be harder by about 0.05 in power-law indices than that of the incident proton and their inclusive cross-sections be larger than those predicted by p-p interaction models based on the Feynman scaling.

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.