Tuneyoshi Kamae

Recent observations found that electrons are accelerated to ∼10 GeV and emit synchrotron hard X-rays in two magnetic white dwarfs (WDs), also known as cataclysmic variables (CVs). In nova outbursts of WDs, multi-GeV gamma-rays were detected, implying that protons are accelerated to 100 GeV or higher. In recent optical surveys, the WD density is found to be higher near the Sun than in the Galactic disk by a factor ∼2.5. The cosmic rays (CRs) produced by local CVs and novae will accumulate in the local bubble for 10⁶–10⁷ yr. On these findings, we search for CRs from historic CVs and novae in the observed CR spectra. We model the CR spectra at the heliopause as sums of Galactic and local components based on observational data as much as possible. The initial Galactic CR electron and proton spectra are deduced from the gamma-ray emissivity, the local electron spectrum from the hard X-ray spectra at the CVs, and the local proton spectrum from gamma-ray spectra at novae. These spectral shapes are then expressed in a simple set of polynomial functions of CR energy and regressively fitted until the high-energy (>100 GeV) CR spectra near Earth and the Voyager-1 spectra at the heliopause are reproduced. We then extend the modeling to nuclear CR spectra and find that one spectral shape fits all local nuclear CRs, and that the apparent hardening of the nuclear CR spectra is caused by the roll-down of local nuclear spectra around 100–200 GeV. All local CR spectra populate a limited energy band below 100–200 GeV and enhance gamma-ray emissivity below ∼10 GeV. Such an enhancement is observed in the inner Galaxy, suggesting the CR fluxes from CVs and novae are substantially higher there.

The Polarized Gamma-ray Observer, PoGOLite, is a balloon experiment with the capability of detecting 10% polarization from a 200 mCrab celestial object between the energy-range 25–80 keV in one 6 hour flight. Polarization measurements in soft gamma-rays are expected to provide a powerful probe into high-energy emission mechanisms in/around neutron stars, black holes, supernova remnants, active-galactic nuclei etc. The “pathfinder” flight was performed in July 2013 for 14 days from Sweden to Russia. The polarization is measured using Compton scattering and photoelectric absorption in an array of 61 well-type phoswich detector cells (PDCs) for the pathfinder instrument. The PDCs are surrounded by 30 BGO crystals which form a side anti-coincidence shield (SAS) and passive polyethylene neutron shield. There is a neutron detector consisting of LiCaAlF6 (LiCAF) scintillator covered with BGOs to measure the background contribution of atmospheric neutrons. The data acquisition system treats 92 PMT signals from 61 PDCs + 30 SASs + 1 neutron detector, and it is developed based on SpaceWire spacecraft communication network. Most of the signal processing is done by digital circuits in Field Programmable Gate Arrays (FPGAs). This enables the reduction of the mass, the space and the power consumption. The performance was calibrated before the launch.

We report on the gamma-ray observations of giant molecular clouds Orion A and B with the Large Area Telescope (LAT) on-board the Fermi Gamma-ray Space Telescope. The gamma-ray emission in the energy band between \sim100 MeV and \sim100 GeV is predicted to trace the gas mass distribution in the clouds through nuclear interactions between the Galactic cosmic rays (CRs) and interstellar gas. The gamma-ray production cross-section for the nuclear interaction is known to \sim10% precision which makes the LAT a powerful tool to measure the gas mass column density distribution of molecular clouds for a known CR intensity. We present here such distributions for Orion A and B, and correlate them with those of the velocity integrated CO intensity (WCO) at a 1{\deg} \times1{\deg} pixel level. The correlation is found to be linear over a WCO range of ~10 fold when divided in 3 regions, suggesting penetration of nuclear CRs to most of the cloud volumes. The Wco-to-mass conversion factor, Xco, is found to be \sim2.3\times10^20 cm-2(K km s-1)-1 for the high-longitude part of Orion A (l > 212{\deg}), \sim1.7 times higher than \sim1.3 \times 10^20 found for the rest of Orion A and B. We interpret the apparent high Xco in the high-longitude region of Orion A in the light of recent works proposing a non-linear relation between H2 and CO densities in the diffuse molecular gas. Wco decreases faster than the H2 column density in the region making the gas "darker" to Wco.

We present a model of cosmic ray injection into the Galactic space based on recent gamma-ray observations of supernova remnants (SNRs) and pulsar wind nebulae (PWNe) by the Fermi Large Area Telescope and imaging atmospheric Cherenkov telescopes. Steady-state (SS) injection of nuclear particles and electrons (e^-) from the Galactic ensemble of SNRs, and electrons and positrons (e^+) from the Galactic ensemble of PWNe are assumed, with their injection spectra inferred under guidance of gamma-ray observations and recent development of evolution and emission models. The ensembles of SNRs and PWNe are assumed to share the same spatial distributions. Assessment of possible secondary CR contribution from dense molecular clouds interacting with SNRs is also given. CR propagation in the interstellar space is handled by GALPROP. Different underlying source distribution models and Galaxy halo sizes are employed to estimate the systematic uncertainty of the model. We show that this observation-based model reproduces the positron fraction e^+/(e^- + e^+) and antiproton-to-proton ratio reported by Pamela and other previous missions reasonably well, without calling for any speculative sources. A discrepancy remains, however, between the total e^- + e^+ spectrum measured by Fermi and our model below 20GeV, for which the potential causes are discussed. Important quantities for Galactic CRs including their energy injection, average lifetime, and mean gas density along their typical propagation path are also estimated.

We report observation of the supernova remnant (SNR) IC 443 (G189.1+3.0) with the Fermi Gamma-ray Space Telescope Large Area Telescope (LAT) in the energy band between 200 MeV and 50 GeV. IC 443 is a shell-type SNR with mixed morphology located off the outer Galactic plane where high-energy emission has been detected in the X-ray, GeV and TeV gamma-ray bands. Past observations suggest IC 443 has been interacting with surrounding interstellar matter. Proximity between dense shocked molecular clouds and GeV–TeV gamma-ray emission regions detected by EGRET, MAGIC, and VERITAS suggests an interpretation that cosmic-ray (CR) particles are accelerated by the SNR. With the high gamma-ray statistics and broad energy coverage provided by the LAT, we accurately characterize the gamma-ray emission produced by the CRs accelerated at IC 443. The emission region is extended in the energy band with θ₆₈ = 0.°27 ± 0.°01(stat) ± 0.°03(sys) for an assumed two-dimensional Gaussian profile and overlaps almost completely with the extended source region of VERITAS. Its centroid is displaced significantly from the known pulsar wind nebula (PWN) which suggests the PWN is not the major contributor in the present energy band. The observed spectrum changes its power-law slope continuously and continues smoothly to the MAGIC and VERITAS data points. The combined gamma-ray spectrum (200 MeV <E< 2 TeV) is reproduced well by decays of neutral pions produced by a broken power-law proton spectrum with a break around 70 GeV.

We present a non-iterative method to deconvolve the spatial response function or the point spread function (PSF) from images taken with the Suzaku X-ray Imaging Spectrometer (XIS). The method is optimized for analyses of extended sources with high photon statistics. Suzaku has four XIS detectors each with its own X-ray CCD and X-Ray Telescope (XRT) and has been providing unique opportunities in spatially-resolved spectroscopic analyses of extended objects. The detectors, however, suffer from broad and position-dependent PSFs with their typical half-power density (HPD) of about 110''. In the authors' view, this shortcoming has been preventing the high collecting area and high spectral resolution of Suzaku to be fully exploited. The present method is intended to recover spatial resolution to ~15'' over a dynamic range around 1:100 in the brightness without assuming any source model. Our deconvolution proceeds in two steps: An XIS image is multiplied with the inverse response matrix calculated from its PSF after rebinning CCD pixels to larger-size tiles (typically 6''x 6''); The inverted image is then adaptively smoothed to obtain the final deconvolved image. The PSF is modeled on a ray-tracing program and an observed point-source image. The deconvolution method has been applied to images of Centaurus A, PSR B1509-58 and RCW 89 taken by one XIS (XIS-1). The results have been compared with images obtained with Chandra to conclude that the spatial resolution has been recovered to ~20'' down to regions where surface brightness is about 1:50 of the brightest tile in the image. We believe the spatial resolution and the dynamic range can be improved in the future with higher fidelity PSF modeling and higher precision pointing information.

We describe a new balloon-borne instrument (PoGOLite) capable of detecting 10% polarisation from 200mCrab point-like sources between 25 and 80keV in one 6 hour flight. Polarisation measurements in the soft gamma-ray band are expected to provide a powerful probe into high-energy emission mechanisms as well as the distribution of magnetic fields, radiation fields and interstellar matter. At present, only exploratory polarisation measurements have been carried out in the soft gamma-ray band. Reduction of the large background produced by cosmic-ray particles has been the biggest challenge. PoGOLite uses Compton scattering and photo-absorption in an array of 217 well-type phoswich detector cells made of plastic and BGO scintillators surrounded by a BGO anticoincidence shield and a thick polyethylene neutron shield. The narrow FOV (1.25msr) obtained with well-type phoswich detector technology and the use of thick background shields enhance the detected S/N ratio. Event selections based on recorded phototube waveforms and Compton kinematics reduce the background to that expected for a 40-100mCrab source between 25 and 50keV. A 6 hour observation on the Crab will differentiate between the Polar Cap/Slot Gap, Outer Gap, and Caustic models with greater than 5 sigma; and also cleanly identify the Compton reflection component in the Cygnus X-1 hard state. The first flight is planned for 2010 and long-duration flights from Sweden to Northern Canada are foreseen thereafter.

We present the angular distribution of gamma rays produced by proton-proton interactions in parameterized formulae to facilitate calculations in astrophysical environments. The parameterization is derived from Monte Carlo simulations of the up-to-date proton-proton interaction model by Kamae et al. (2005) and its extension by Kamae et al. (2006). This model includes the logarithmically rising inelastic cross section, the diffraction dissociation process and Feynman scaling violation. The extension adds two baryon resonance contributions: one representing the Delta(1232) and the other representing multiple resonances around 1600 MeV/c^2. We demonstrate the use of the formulae by calculating the predicted gamma-ray spectrum for two different cases: the first is a pencil beam of protons following a power law and the second is a fanned proton jet with a Gaussian intensity profile impinging on the surrounding material. In both cases we find that the predicted gamma-ray spectrum to be dependent on the viewing angle.

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.