The Fermi bubbles were possibly created by large injections of energy into the Galactic Center (GC), either by an active galactic nucleus (AGN) or by nuclear starburst more than ~10 Myr ago. However, the origin of the diffuse gamma-ray emission associated with Loop I, a radio continuum loop spanning across 100 deg on the sky, is still being debated. The northern-most part of Loop I, known as the North Polar Spur (NPS), is the brightest arm and is even clearly visible in the ROSAT X-ray sky map. In this paper, we present a comprehensive review on the X-ray observations of the Fermi bubbles and their possible association with the NPS and Loop I structures. Using uniform analysis of archival Suzaku and Swift data, we show that X-ray plasma with kT ~ 0.3 keV and low metal abundance (Z ~ 0.2 Z_solar) is ubiquitous in both the bubbles and Loop I and is naturally interpreted as weakly shock-heated Galactic halo gas. However, the observed asymmetry of the X-ray-emitting gas above and below the GC has still not been resolved; it cannot be fully explained by the inclination of the axis of the Fermi bubbles to the Galactic disk normal. We argue that the NPS and Loop I may be asymmetric remnants of a large explosion that occurred before the event that created the Fermi bubbles, and that the soft gamma-ray emission from Loop I may be due to either pi^0 decay of accelerated protons or electron bremsstrahlung.

We present results from the Hitomi X-ray observation of a young composite-type supernova remnant (SNR) G21.5$-$0.9, whose emission is dominated by the pulsar wind nebula (PWN) contribution. The X-ray spectra in the 0.8-80 keV range obtained with the Soft X-ray Spectrometer (SXS), Soft X-ray Imager (SXI) and Hard X-ray Imager (HXI) show a significant break in the continuum as previously found with the NuSTAR observation. After taking into account all known emissions from the SNR other than the PWN itself, we find that the Hitomi spectra can be fitted with a broken power law with photon indices of $\Gamma_1=1.74\pm0.02$ and $\Gamma_2=2.14\pm0.01$ below and above the break at $7.1\pm0.3$ keV, which is significantly lower than the NuSTAR result ($\sim9.0$ keV). The spectral break cannot be reproduced by time-dependent particle injection one-zone spectral energy distribution models, which strongly indicates that a more complex emission model is needed, as suggested by recent theoretical models. We also search for narrow emission or absorption lines with the SXS, and perform a timing analysis of PSR J1833$-$1034 with the HXI and SGD. No significant pulsation is found from the pulsar. However, unexpectedly, narrow absorption line features are detected in the SXS data at 4.2345 keV and 9.296 keV with a significance of 3.65 $\sigma$. While the origin of these features is not understood, their mere detection opens up a new field of research and was only possible with the high resolution, sensitivity and ability to measure extended sources provided by an X-ray microcalorimeter.

The present paper investigates the temperature structure of the X-ray emitting plasma in the core of the Perseus cluster using the 1.8--20.0 keV data obtained with the Soft X-ray Spectrometer (SXS) onboard the Hitomi Observatory. A series of four observations were carried out, with a total effective exposure time of 338 ks and covering a central region $\sim7'$ in diameter. The SXS was operated with an energy resolution of $\sim$5 eV (full width at half maximum) at 5.9 keV. Not only fine structures of K-shell lines in He-like ions but also transitions from higher principal quantum numbers are clearly resolved from Si through Fe. This enables us to perform temperature diagnostics using the line ratios of Si, S, Ar, Ca, and Fe, and to provide the first direct measurement of the excitation temperature and ionization temperature in the Perseus cluster. The observed spectrum is roughly reproduced by a single temperature thermal plasma model in collisional ionization equilibrium, but detailed line ratio diagnostics reveal slight deviations from this approximation. In particular, the data exhibit an apparent trend of increasing ionization temperature with increasing atomic mass, as well as small differences between the ionization and excitation temperatures for Fe, the only element for which both temperatures can be measured. The best-fit two-temperature models suggest a combination of 3 and 5 keV gas, which is consistent with the idea that the observed small deviations from a single temperature approximation are due to the effects of projection of the known radial temperature gradient in the cluster core along the line of sight. Comparison with the Chandra/ACIS and the XMM-Newton/RGS results on the other hand suggests that additional lower-temperature components are present in the ICM but not detectable by Hitomi SXS given its 1.8--20 keV energy band.

The Hitomi SXS spectrum of the Perseus cluster, with $\sim$5 eV resolution in the 2-9 keV band, offers an unprecedented benchmark of the atomic modeling and database for hot collisional plasmas. It reveals both successes and challenges of the current atomic codes. The latest versions of AtomDB/APEC (3.0.8), SPEX (3.03.00), and CHIANTI (8.0) all provide reasonable fits to the broad-band spectrum, and are in close agreement on best-fit temperature, emission measure, and abundances of a few elements such as Ni. For the Fe abundance, the APEC and SPEX measurements differ by 16%, which is 17 times higher than the statistical uncertainty. This is mostly attributed to the differences in adopted collisional excitation and dielectronic recombination rates of the strongest emission lines. We further investigate and compare the sensitivity of the derived physical parameters to the astrophysical source modeling and instrumental effects. The Hitomi results show that an accurate atomic code is as important as the astrophysical modeling and instrumental calibration aspects. Substantial updates of atomic databases and targeted laboratory measurements are needed to get the current codes ready for the data from the next Hitomi-level mission.

We present Hitomi observations of N132D, a young, X-ray bright, O-rich core-collapse supernova remnant in the Large Magellanic Cloud (LMC). Despite a very short observation of only 3.7 ks, the Soft X-ray Spectrometer (SXS) easily detects the line complexes of highly ionized S K and Fe K with 16-17 counts in each. The Fe feature is measured for the first time at high spectral resolution. Based on the plausible assumption that the Fe K emission is dominated by He-like ions, we find that the material responsible for this Fe emission is highly redshifted at ~800 km/s compared to the local LMC interstellar medium (ISM), with a 90% credible interval of 50-1500 km/s if a weakly informative prior is placed on possible line broadening. This indicates (1) that the Fe emission arises from the supernova ejecta, and (2) that these ejecta are highly asymmetric, since no blue-shifted component is found. The S K velocity is consistent with the local LMC ISM, and is likely from swept-up ISM material. These results are consistent with spatial mapping that shows the He-like Fe concentrated in the interior of the remnant and the S tracing the outer shell. The results also show that even with a very small number of counts, direct velocity measurements from Doppler-shifted lines detected in extended objects like supernova remnants are now possible. Thanks to the very low SXS background of ~1 event per spectral resolution element per 100 ks, such results are obtainable during short pointed or slew observations with similar instruments. This highlights the power of high-spectral-resolution imaging observations, and demonstrates the new window that has been opened with Hitomi and will be greatly widened with future missions such as the X-ray Astronomy Recovery Mission (XARM) and Athena.

We report a Hitomi observation of IGR J16318-4848, a high-mass X-ray binary system with an extremely strong absorption of N_H~10^{24} cm^{-2}. Previous X-ray studies revealed that its spectrum is dominated by strong fluorescence lines of Fe as well as continuum emission. For physical and geometrical insight into the nature of the reprocessing material, we utilize the high spectroscopic resolving power of the X-ray microcalorimeter (the soft X-ray spectrometer; SXS) and the wide-band sensitivity by the soft and hard X-ray imager (SXI and HXI) aboard Hitomi. Even though photon counts are limited due to unintended off-axis pointing, the SXS spectrum resolves Fe K{\alpha_1} and K{\alpha_2} lines and puts strong constraints on the line centroid and width. The line width corresponds to the velocity of 160^{+300}_{-70} km s^{-1}. This represents the most accurate, and smallest, width measurement of this line made so far from any X-ray binary, much less than the Doppler broadening and shift expected from speeds which are characteristic of similar systems. Combined with the K-shell edge energy measured by the SXI and HXI spectra, the ionization state of Fe is estimated to be in the range of Fe I--IV. Considering the estimated ionization parameter and the distance between the X-ray source and the absorber, the density and thickness of the materials are estimated. The extraordinarily strong absorption and the absence of a Compton shoulder component is confirmed. These characteristics suggest reprocessing materials which are distributed in a narrow solid angle or scattering primarily with warm free electrons or neutral hydrogen.

The origin of the narrow Fe-K{\alpha} fluorescence line at 6.4 keV from active galactic nuclei has long been under debate; some of the possible sites are the outer accretion disk, the broad line region, a molecular torus, or interstellar/intracluster media. In February-March 2016, we performed the first X-ray microcalorimeter spectroscopy with the Soft X-ray Spectrometer (SXS) onboard the Hitomi satellite of the Fanaroff-Riley type I radio galaxy NGC 1275 at the center of the Perseus cluster of galaxies. With the high energy resolution of ~5 eV at 6 keV achieved by Hitomi/SXS, we detected the Fe-K{\alpha} line with ~5.4 {\sigma} significance. The velocity width is constrained to be 500-1600 km s$^{-1}$ (FWHM for Gaussian models) at 90% confidence. The SXS also constrains the continuum level from the NGC 1275 nucleus up to ~20 keV, giving an equivalent width ~20 eV of the 6.4 keV line. Because the velocity width is narrower than that of broad H{\alpha} line of ~2750 km s$^{-1}$, we can exclude a large contribution to the line flux from the accretion disk and the broad line region. Furthermore, we performed pixel map analyses on the Hitomi/SXS data and image analyses on the Chandra archival data, and revealed that the Fe-K{\alpha} line comes from a region within ~1.6 kpc from the NGC 1275 core, where an active galactic nucleus emission dominates, rather than that from intracluster media. Therefore, we suggest that the source of the Fe-K{\alpha} line from NGC 1275 is likely a low-covering fraction molecular torus or a rotating molecular disk which probably extends from a pc to hundreds pc scale in the active galactic nucleus system.

Extending the earlier measurements reported in Hitomi collaboration (2016, Nature, 535, 117), we examine the atmospheric gas motions within the central 100~kpc of the Perseus cluster using observations obtained with the Hitomi satellite. After correcting for the point spread function of the telescope and using optically thin emission lines, we find that the line-of-sight velocity dispersion of the hot gas is remarkably low and mostly uniform. The velocity dispersion reaches maxima of approximately 200~km~s$^{-1}$ toward the central active galactic nucleus (AGN) and toward the AGN inflated north-western `ghost' bubble. Elsewhere within the observed region, the velocity dispersion appears constant around 100~km~s$^{-1}$. We also detect a velocity gradient with a 100~km~s$^{-1}$ amplitude across the cluster core, consistent with large-scale sloshing of the core gas. If the observed gas motions are isotropic, the kinetic pressure support is less than 10\% of the thermal pressure support in the cluster core. The well-resolved optically thin emission lines have Gaussian shapes, indicating that the turbulent driving scale is likely below 100~kpc, which is consistent with the size of the AGN jet inflated bubbles. We also report the first measurement of the ion temperature in the intracluster medium, which we find to be consistent with the electron temperature. In addition, we present a new measurement of the redshift to the brightest cluster galaxy NGC~1275.

Thanks to its high spectral resolution (~5 eV at 6 keV), the Soft X-ray Spectrometer (SXS) on board Hitomi enables us to measure the detailed structure of spatially resolved emission lines from highly ionized ions in galaxy clusters for the first time. In this series of papers, using the SXS we have measured the velocities of gas motions, metallicities and the multi-temperature structure of the gas in the core of the Perseus cluster. Here, we show that when inferring physical properties from line emissivities in systems like Perseus, the resonant scattering (RS) effect should be taken into account. In the Hitomi waveband, RS mostly affects the FeXXV He$\alpha$ line ($w$) - the strongest line in the spectrum. The flux measured by Hitomi in this line is suppressed by a factor ~1.3 in the inner ~30 kpc, compared to predictions for an optically thin plasma; the suppression decreases with the distance from the center. The $w$ line also appears slightly broader than other lines from the same ion. The observed distortions of the $w$ line flux, shape and distance dependence are all consistent with the expected effect of the resonant scattering in the Perseus core. By measuring the ratio of fluxes in optically thick ($w$) and thin (FeXXV forbidden, He$\beta$, Ly$\alpha$) lines, and comparing these ratios with predictions from Monte Carlo radiative transfer simulations, the velocities of gas motions have been obtained. The results are consistent with the direct measurements of gas velocities from line broadening described elsewhere in this series, although the systematic and statistical uncertainties remain significant. Further improvements in the predictions of line emissivities in plasma models, and deeper observations with future X-ray missions will enable RS measurements to provide powerful constraints on the amplitude and anisotropy of clusters gas motions.

We set a new upper limit on the abundance of primordial black holes (PBH) based on existing X-ray data. PBH interactions with interstellar medium should result in significant fluxes of X-ray photons, which would contribute to the observed number density of compact X-ray objects in galaxies. The data constrain PBH number density in the mass range from a few $M_\odot$ to $2\times 10^7 M_\odot$. PBH density needed to account for the origin of black holes detected by LIGO is marginally allowed.

One of the major unsolved questions on the understanding of the AGN population is the origin of the dichotomy between radio-quiet and radio-loud quasars. The most promising explanation is provided by the spin paradigm, which suggests radio-loud quasars have higher black hole spin. However, the measurement of black hole spin remains extremely challenging. We here aim at comparing the mean radiative efficiencies of carefully matched samples of radio-loud and radio-quiet SDSS quasars at 0.3<z<0.8. We use the [OIII] luminosity as an indirect average tracer of the ionizing continuum in the extreme-UV regime where differences in the SED due to black hole spin are most pronounced. We find that the radio-loud sample shows an enhancement in [OIII] line strength by a factor of at least 1.5 compared to a radio-quiet sample matched in redshift, black hole mass and optical continuum luminosity or accretion rate. We argue that this enhancement is caused by differences in the SED, suggesting higher average bolometric luminosities at fixed accretion rate in the radio-loud population. This suggests that the radio-loud quasar population has on average systematically larger radiative efficiencies and therefore higher black hole spin than the radio-quiet population, providing observational support for the black hole spin paradigm.

The Crab nebula originated from a core-collapse supernova (SN) explosion observed in 1054 A.D. When viewed as a supernova remnant (SNR), it has an anomalously low observed ejecta mass and kinetic energy for an Fe-core collapse SN. Intensive searches were made for a massive shell that solves this discrepancy, but none has been detected. An alternative idea is that the SN1054 is an electron-capture (EC) explosion with a lower explosion energy by an order of magnitude than Fe-core collapse SNe. In the X-rays, imaging searches were performed for the plasma emission from the shell in the Crab outskirts to set a stringent upper limit to the X-ray emitting mass. However, the extreme brightness of the source hampers access to its vicinity. We thus employed spectroscopic technique using the X-ray micro-calorimeter onboard the Hitomi satellite. By exploiting its superb energy resolution, we set an upper limit for emission or absorption features from yet undetected thermal plasma in the 2-12 keV range. We also re-evaluated the existing Chandra and XMM-Newton data. By assembling these results, a new upper limit was obtained for the X-ray plasma mass of <~ 1Mo for a wide range of assumed shell radius, size, and plasma temperature both in and out of the collisional equilibrium. To compare with the observation, we further performed hydrodynamic simulations of the Crab SNR for two SN models (Fe-core versus EC) under two SN environments (uniform ISM versus progenitor wind). We found that the observed mass limit can be compatible with both SN models if the SN environment has a low density of <~ 0.03 cm-3 (Fe core) or <~ 0.1 cm-3 (EC) for the uniform density, or a progenitor wind density somewhat less than that provided by a mass loss rate of 10-5 Mo yr-1 at 20 km s-1 for the wind environment.

The Cherenkov Telescope Array (CTA) is about to enter construction phase and one of its main key science projects is to perform an unbiased survey in search of extragalactic sources. We make use of both the latest blazar gamma--ray luminosity function and spectral energy distribution to derive the expected number of detectable sources for both the planned Northern and Southern arrays of the CTA observatory. We find that a shallow, wide survey of about 0.5 hour per field of view would lead to the highest number of blazar detections. Furthermore, we investigate the effect of axion-like particles and secondary gamma rays from propagating cosmic rays on the source count distribution, since these processes predict different spectral shape from standard extragalactic background light attenuation. We can generally expect more distant objects in the secondary gamma-ray scenario, while axion-like particles do not significantly alter the expected distribution. Yet, we find that, these results strongly depend on the assumed magnetic field strength during the propagation. We also provide source count predictions for the High Altitude Water Cherenkov observatory (HAWC), the Large High Altitude Air Shower Observatory (LHAASO) and a novel proposal of a hybrid detector.

We study the disk-jet connection in supermassive black holes by investigating the properties of their optical and radio emissions utilizing the SDSS-DR7 and the NVSS catalogs. Our sample contains 7017 radio-loud quasars with detection both at 1.4~GHz and SDSS optical spectrum. Using this radio-loud quasar sample, we investigate the correlation among the jet power ($P_{\rm jet}$), the bolometric disk luminosity ($L_{\rm disk}$), and the black hole mass ($M_{\rm BH}$) in the standard accretion disk regime. We find that the jet powers correlate with the bolometric disk luminosities as $\log P_{\rm jet} = (0.96\pm0.012)\log L_{\rm disk} + (0.79 \pm 0.55)$. This suggests that the jet production efficiency of $\eta_{\rm jet}\simeq1.1_{-0.76}^{+2.6}\times10^{-2}$ assuming the disk radiative efficiency of $0.1$ implying low black hole spin parameters and/or low magnetic flux for radio-loud quasars. But it can be also due to dependence of the efficiency on geometrical thickness of the accretion flow which is expected to be small for quasars accreting at the disk Eddington ratios $0.01 \lesssim \lambda \lesssim 0.3$. This low jet production efficiency does not significantly increase even if we set the disk radiative efficiency of 0.3. We also investigate the fundamental plane in our samples among $P_{\rm jet}$, $L_{\rm disk}$, and $M_{\rm BH}$. We could not find a statistically significant fundamental plane for radio-loud quasars in the standard accretion regime.

High-resolution X-ray spectroscopy with Hitomi was expected to resolve the origin of the faint unidentified E=3.5 keV emission line reported in several low-resolution studies of various massive systems, such as galaxies and clusters, including the Perseus cluster. We have analyzed the Hitomi first-light observation of the Perseus cluster. The emission line expected for Perseus based on the XMM-Newton signal from the large cluster sample under the dark matter decay scenario is too faint to be detectable in the Hitomi data. However, the previously reported 3.5 keV flux from Perseus was anomalously high compared to the sample-based prediction. We find no unidentified line at the reported high flux level. Taking into account the XMM measurement uncertainties for this region, the inconsistency with Hitomi is at a 99% significance for a broad dark-matter line and at 99.7% for a narrow line from the gas. We do not find anomalously high fluxes of the nearby faint K line or the Ar satellite line that were proposed as explanations for the earlier 3.5 keV detections. We do find a hint of a broad excess near the energies of high-n transitions of Sxvi (E=3.44 keV rest-frame) -- a possible signature of charge exchange in the molecular nebula and another proposed explanation for the unidentified line. While its energy is consistent with XMM pn detections, it is unlikely to explain the MOS signal. A confirmation of this interesting feature has to wait for a more sensitive observation with a future calorimeter experiment.

In the extragalactic sky, microquasars and ultra-luminous X-ray sources (ULXs) are known as energetic compact objects locating at off-nucleus positions in galaxies. Some of these objects are associated with expanding bubbles with a velocity of 80-250 ${\rm km~s^{-1}}$. We investigate the shock acceleration of particles in those expanding nebulae. The nebulae having fast expansion velocity $\gtrsim120~{\rm km~s^{-1}}$ are able to accelerate cosmic rays up to $\sim100$ TeV. If 10% of the shock kinetic energy goes into particle acceleration, powerful nebulae such as the microquasar S26 in NGC 7793 would emit gamma rays up to several tens TeV with a photon index of $\sim2$. These nebulae will be good targets for future Cherenkov Telescope Array observations given its sensitivity and angular resolution. They would also contribute to $\sim7$% of the unresolved cosmic gamma-ray background radiation at $\ge0.1~{\rm GeV}$. In contrast, particle acceleration in slowly expanding nebulae $\lesssim120~{\rm km~s^{-1}}$ would be less efficient due to ion-neutral collisions and result in softer spectra at $\gtrsim10$ GeV.

Clusters of galaxies are the most massive gravitationally-bound objects in the Universe and are still forming. They are thus important probes of cosmological parameters and a host of astrophysical processes. Knowledge of the dynamics of the pervasive hot gas, which dominates in mass over stars in a cluster, is a crucial missing ingredient. It can enable new insights into mechanical energy injection by the central supermassive black hole and the use of hydrostatic equilibrium for the determination of cluster masses. X-rays from the core of the Perseus cluster are emitted by the 50 million K diffuse hot plasma filling its gravitational potential well. The Active Galactic Nucleus of the central galaxy NGC1275 is pumping jetted energy into the surrounding intracluster medium, creating buoyant bubbles filled with relativistic plasma. These likely induce motions in the intracluster medium and heat the inner gas preventing runaway radiative cooling; a process known as Active Galactic Nucleus Feedback. Here we report on Hitomi X-ray observations of the Perseus cluster core, which reveal a remarkably quiescent atmosphere where the gas has a line-of-sight velocity dispersion of 164+/-10 km/s in a region 30-60 kpc from the central nucleus. A gradient in the line-of-sight velocity of 150+/-70 km/s is found across the 60 kpc image of the cluster core. Turbulent pressure support in the gas is 4% or less of the thermodynamic pressure, with large scale shear at most doubling that estimate. We infer that total cluster masses determined from hydrostatic equilibrium in the central regions need little correction for turbulent pressure.

The Advanced Laser Interferometer Gravitational-Wave Observatory (aLIGO) has detected direct signals of gravitational waves (GWs) from GW150914. The event was a merger of binary black holes whose masses are $36^{+5}_{-4}M_{\odot}$ and $29^{+4}_{-4}M_{\odot}$. Such binary systems are expected to be directly evolved from stellar binary systems or formed by dynamical interactions of black holes in dense stellar environments. Here we derive the binary black hole merger rate based on the nearby ultra-luminous X-ray source (ULX) luminosity function (LF) under the assumption that binary black holes evolve through X-ray emitting phases. We obtain the binary black hole merger rate as $5.8 ({t}_{\rm ULX}/{0.1 \ \rm Myr})^{-1} \lambda^{-0.6} \exp{(-0.30\lambda)} \ {\rm Gpc^{-3}\ yr^{-1}}$, where $t_{\rm ULX}$ is the typical duration of the ULX phase and $\lambda$ is the Eddington ratio in luminosity. This is coincident with the event rate inferred from the detection of GW150914 as well as the predictions based on binary population synthesis models. Although we are currently unable to constrain the Eddington ratio of ULXs in luminosity due to the uncertainties of our models and measured binary black hole merger event rates, further X-ray and GW data will allow us to narrow down the range of the Eddington ratios of ULXs. We also find the cumulative merger rate for the mass range of $5M_\odot\le M_{\rm BH}\le100M_\odot$ inferred from the ULX LF is consistent with that estimated by the aLIGO collaboration considering various astrophysical conditions such as the mass function of black holes.

Relativistic jets launched by supermassive black holes, so-called as active galactic nuclei (AGNs), are known as the most energetic particle accelerators in the universe. However, the baryon loading efficiency onto the jets from the accretion flows and their particle acceleration efficiencies have been veiled in mystery. With the latest data sets, we perform multi-wavelength spectral analysis of quiescent spectra of 13 TeV gamma-ray detected high-frequency-peaked BL Lacs (HBLs) following one-zone static synchrotron-self-Compton (SSC) model. We determine the minimum, cooling break, and maximum electron Lorentz factors following the diffusive shock acceleration (DSA) theory. We find that HBLs have $P_B/P_e\sim6.3\times10^{-3}$ and the radiative efficiency $\epsilon_{\rm rad,jet}\sim6.7\times10^{-4}$ where $P_B$ and $P_e$ is the Poynting and electron power, respectively. By assuming 10 leptons per one proton, the jet power relates to the black hole mass as $P_{\rm jet}/L_{\rm Edd}\sim0.18$ where $P_{\rm jet}$ and $L_{\rm Edd}$ is the jet power and the Eddington luminosity, respectively. Under our model assumptions, we further find that HBLs have the jet production efficiency of $\eta_{\rm jet}\sim1.5$ and the mass loading efficiency of $\xi_{\rm jet}\gtrsim5\times10^{-2}$. We also investigate the particle acceleration efficiency in the blazar zone by including the most recent Swift/BAT data. Our samples ubiquitously have the particle acceleration efficiency of $\eta_g\sim10^{4.5}$, which is inefficient to accelerate particles up to the ultra-high-energy-cosmic-ray (UHECR) regime. This implies that the UHECR acceleration sites should be other than the blazar zones of quiescent low power AGN jets, if one assumes the one-zone SSC model based on the DSA theory.

The Seyfert galaxy NGC 985 is known to show a high-frequency excess in its radio continuum spectrum in a milli-Jansky level on the basis of previous observations at 1.4--15 GHz; a steep spectrum at low frequencies (a spectral index of $\alpha=-1.10 \pm 0.03$) changes at ~10 GHz into an inverted spectrum at higher frequencies ($\alpha=+0.86 \pm 0.09$). We conduct new observations at 15--43 GHz using the Very Large Array and at 100 GHz using the Nobeyama Millimeter Array. As a result, the high-frequency excess continuing at even higher radio frequencies up to 43 GHz has been confirmed. The non-detection at 100 GHz was not so strong constraint, and therefore the spectral behavior above 43 GHz remains unclear. The astrometric position of the high-frequency excess component coincides with the optical position of the Seyfert nucleus and the low-frequency radio position to an accuracy of 0.1 arcsec, corresponding to ~80 pc; the radio source size is constrained to be <0.02 arcsec, corresponding to <16 pc. We discuss the physical origin of the observed high-frequency excess component. Dust emission at the Rayleigh-Jeans regime, free--free emission from X-ray radiating high-temperature plasma, free--free emission from the ensemble of broad-line region clouds, or thermal synchrotron from hot accretion flow cannot be responsible for the observed radio flux. Compact jets under synchrotron self-absorption may be unlikely in terms of observed time scales. Alternatively, we cannot rule out the hypotheses of synchrotron jets free--free absorbed by a circumnuclear photo-ionized region, and self-absorbed nonthermal synchrotron from disk corona, as the origin of the high-frequency excess component.

The Fermi gamma-ray space telescope has revolutionized our understanding of the cosmic gamma-ray background radiation in the GeV band. However, investigation on the cosmic TeV gamma-ray background radiation still remains sparse. Here, we report the lower bound on the cosmic TeV gamma-ray background spectrum placed by the cumulative flux of individual detected extragalactic TeV sources including blazars, radio galaxies, and starburst galaxies. The current limit on the cosmic TeV gamma-ray background above 0.1 TeV is obtained as $2.8\times10^{-8} (E/100~{\rm GeV})^{-0.55} \exp(-E/2100~{\rm GeV})~{\rm [GeV/cm^2/s/sr]} < E^2dN/dE < 1.1\times10^{-7} (E/100~{\rm GeV})^{-0.49}~{\rm [GeV/cm^2/s/sr]}$, where the upper bound is set by requirement that the cascade flux from the cosmic TeV gamma-ray background radiation can not exceed the measured cosmic GeV gamma-ray background spectrum (Inoue & Ioka, 2012). Two nearby blazars, Mrk 421 and Mrk 501, explain ~70% of the cumulative background flux at 0.8-4 TeV, while extreme blazars start to dominate at higher energies. We also provide the cumulative background flux from each population, i.e. blazars, radio galaxies, and starburst galaxies which will be the minimum requirement for their contribution to the cosmic TeV gamma-ray background radiation.

While the X-ray, GeV gamma-ray, and TeV gamma-ray skies have been extensively studied, the MeV gamma-ray sky is not well investigated after the Imaging Compton Telescope (COMPTEL) scanned the sky about two decades ago. In this paper, we investigate prospects for active galactic nuclei population studies with future MeV gamma-ray missions using recent spectral models and luminosity functions of Seyfert and flat spectrum radio quasars (FSRQs). Both of them are plausible candidates as the origins of the cosmic MeV gamma-ray background. If the cosmic MeV gamma-ray background radiation is dominated by non-thermal emission from Seyferts, the sensitivity of 10^-12 erg cm^-2 s^-1 is required to detect several hundred Seyferts in the entire sky. If FSRQs make up the cosmic MeV gamma-ray background, the sensitivity of ~4 x 10^-12 erg cm^-2 s^-1 is required to detect several hundred FSRQs following the recent FSRQ X-ray luminosity function. However, based on the latest FSRQ gamma-ray luminosity function, with which FSRQs can explain up to ~30% of the MeV background, we can expect several hundred FSRQs even with the sensitivity of 10^-11 erg cm^-2 s^-1 which is almost the same as the sensitivity goal of the next generation MeV telescopes.

The Fermi bubbles are gigantic gamma-ray structures in our Galaxy. The physical origin of the bubbles is still under debate. The leading scenarios can be divided into two categories. One is the nuclear star forming activity similar to extragalactic starburst galaxies and the other is the past active galactic nucleus (AGN) like activity of the Galactic center supermassive black hole. In this paper, we propose that metal abundance measurements will provide an important clue to probe their origin. Based on a simple spherically symmetric bubble model, we find that the generated metallicity and abundance pattern of the bubbles' gas strongly depend on assumed star formation or AGN activities. Star formation scenarios predict higher metallicities and abundance ratios of [O/Fe] and [Ne/Fe] than AGN scenarios do because of supernovae ejecta. Furthermore, the resultant abundance depends on the gamma-ray emission process because different mass injection histories are required for the different gamma-ray emission processes due to the acceleration and cooling time scales of non-thermal particles. Future X-ray missions such as ASTRO-H and Athena will give a clue to probe the origin of the bubbles through abundance measurements with their high energy resolution instruments.

The cosmic gamma-ray background radiation is one of the most fundamental observables in the gamma-ray band. Although the origin of the cosmic gamma-ray background radiation has been a mystery for a long time, the Fermi gamma-ray space telescope has recently measured it at 0.1-820 GeV and revealed that the cosmic GeV gamma-ray background is composed of blazars, radio galaxies, and star-forming galaxies. However, Fermi still leaves the following questions. Those are dark matter contribution, origins of the cosmic MeV gamma-ray background, and the connection to the IceCube TeV-PeV neutrino events. In this proceeding, I will review the current understandings of the cosmic gamma-ray background and discuss future prospects of cosmic gamma-ray background radiation studies. I also briefly review the current status of cosmic infrared/optical background radiation studies.

The joint JAXA/NASA ASTRO-H mission is the sixth in a series of highly successful X-ray missions developed by the Institute of Space and Astronautical Science (ISAS), with a planned launch in 2015. The ASTRO-H mission is equipped with a suite of sensitive instruments with the highest energy resolution ever achieved at E > 3 keV and a wide energy range spanning four decades in energy from soft X-rays to gamma-rays. The simultaneous broad band pass, coupled with the high spectral resolution of Delta E < 7 eV of the micro-calorimeter, will enable a wide variety of important science themes to be pursued. ASTRO-H is expected to provide breakthrough results in scientific areas as diverse as the large-scale structure of the Universe and its evolution, the behavior of matter in the gravitational strong field regime, the physical conditions in sites of cosmic-ray acceleration, and the distribution of dark matter in galaxy clusters at different redshifts.