Source author record

Yoshiyuki Inoue

Yoshiyuki Inoue appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

43works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

43 published item(s)

preprint2026arXiv

Coherence of Supermassive Black Hole Binary Demographics with the nHz Stochastic Gravitational Wave Background

We present a refined estimation of the stochastic gravitational wave background (SGWB) based on observed dual active galactic nuclei (AGNs) together with AGN X-ray luminosity functions, in light of recent Pulsar Timing Array detections of an nHz SGWB. We identify a characteristic luminosity dependence in dual AGN fractions by compiling recent observational datasets, providing crucial constraints on supermassive black hole binary (SMBHB) populations. Our AGN-based model reproduces the current SGWB measurements within PTA observational uncertainties of $2 - 4 σ$ uncertainties, demonstrating consistency between electromagnetic and gravitational wave observations. These findings establish SMBHBs as the dominant source of the nHz gravitational wave signal, providing crucial insights into their demographics and evolution.

preprint2026arXiv

Proton Acceleration by Collisionless Shocks in Supermassive Black Hole Coronae: Implications for High-Energy Neutrinos

Recent observations by the IceCube Neutrino Observatory have revealed a significant excess of high-energy neutrinos from nearby Seyfert galaxies, such as NGC~1068, without a corresponding flux of high-energy gamma-rays. This suggests that neutrinos are produced via hadronic interactions in a region opaque to gamma-rays, likely a hot corona surrounding the central supermassive black hole. However, the mechanism responsible for accelerating the parent protons to the required energies ($\sim 100$ TeV) remains an open question. In this study, we investigate diffusive shock acceleration (DSA) in active galactic nucleus (AGN) coronae using a suite of one-dimensional Particle-in-cell (PIC) simulations spanning a broad range of plasma parameters. We find that DSA is a robust and efficient mechanism for proton acceleration, consistently channeling approximately 10\% of the shock's kinetic energy into non-thermal ions, even for shocks with sonic Mach number as low as $ M_s \approx 2$. In contrast, the efficiency of electron acceleration is highly variable and less efficient ($<1\%$) in our parameter survey. These findings provide strong, first-principles support for the hadronic models of neutrino production in AGN and offer quantitative constraints that can explain the observed gamma-ray deficit.

preprint2022arXiv

Advancing the Landscape of Multimessenger Science in the Next Decade

The last decade has brought about a profound transformation in multimessenger science. Ten years ago, facilities had been built or were under construction that would eventually discover the nature of objects in our universe could be detected through multiple messengers. Nonetheless, multimessenger science was hardly more than a dream. The rewards for our foresight were finally realized through IceCube's discovery of the diffuse astrophysical neutrino flux, the first observation of gravitational waves by LIGO, and the first joint detections in gravitational waves and photons and in neutrinos and photons. Today we live in the dawn of the multimessenger era. The successes of the multimessenger campaigns of the last decade have pushed multimessenger science to the forefront of priority science areas in both the particle physics and the astrophysics communities. Multimessenger science provides new methods of testing fundamental theories about the nature of matter and energy, particularly in conditions that are not reproducible on Earth. This white paper will present the science and facilities that will provide opportunities for the particle physics community renew its commitment and maintain its leadership in multimessenger science.

preprint2022arXiv

ALMA detection of parsec-scale blobs at the head of kiloparsec-scale jet in the nearby Seyfert galaxy NGC 1068

We present Atacama Large Millimeter/submillimeter Array observations at $\approx100$ GHz with $0.05$ arcsec (3 pc) resolution of the kiloparsec-scale jet seen in the nearby Seyfert galaxy NGC 1068, and we report the presence of parsec-scale blobs at the head of the jet. The combination of the detected radio flux ($\approx0.8$ mJy), spectral index ($\approx0.5$), and the blob size ($\approx10$ pc) suggests a strong magnetic field of $B\approx240\,μ$G. Such a strong magnetic field most likely implies magnetic field amplification by streaming cosmic rays. The estimated cosmic-ray power by the jet may exceed the limit set by the star formation activity in this galaxy. This result suggests that even modest-power jets can increase the galactic cosmic-ray content while propagating through the galactic bulge.

preprint2022arXiv

Event reconstruction of Compton telescopes using a multi-task neural network

We have developed a neural network model to perform event reconstruction of Compton telescopes. This model reconstructs events that consist of three or more interactions in a detector. It is essential for Compton telescopes to determine the time order of the gamma-ray interactions and whether the incident photon deposits all energy in a detector or it escapes from the detector. Our model simultaneously predicts these two essential factors using a multi-task neural network with three hidden layers of fully connected nodes. For verification, we have conducted numerical experiments using Monte Carlo simulation, assuming a large-area Compton telescope using liquid argon to measure gamma rays with energies up to $3.0\,\mathrm{MeV}$. The reconstruction model shows excellent performance of event reconstruction for multiple scattering events that consist of up to eight hits. The accuracies of hit order prediction are around $60\%$ while those of escape flags are higher than $70\%$ for up to eight-hit events of $4π$ isotropic photons. Compared with two other algorithms, a classical model and a physics-based probabilistic one, the present neural network method shows high performance in estimation accuracy particularly when the number of scattering is small, 3 or 4. Since simulation data easily optimize the network model, the model can be flexibly applied to a wide variety of Compton telescopes.

preprint2022arXiv

High energy emission component, population, and contribution to the extragalactic gamma-ray background of gamma-ray emitting radio galaxies

In this study, we systematically studied the X-ray to GeV gamma-ray spectra of 61 {\it Fermi} Large Area Telescope (LAT) detected radio galaxies. We found an anticorrelation between peak frequency and peak luminosity in the high-energy spectral component of radio galaxies, similar to blazars. With this sample, we also constructed a gamma-ray luminosity function (GLF) of gamma-ray-loud radio galaxies. We found that blazar-like GLF shapes can reproduce their redshift and luminosity distribution, but the log$N$-log$S$ relation prefers models with more low-$z$ radio galaxies. This indicates many low-$z$ gamma-ray-loud radio galaxies. By utilizing our latest GLF, the contribution of radio galaxies to the extragalactic gamma-ray background is found to be 1--10\%. We further investigated the nature of gamma-ray-loud radio galaxies. Compared to radio or X-ray flux-limited radio galaxy samples, the gamma-ray selected sample tends to lack high radio power galaxies like FR-II radio galaxies. We also found that only $\sim$10\% of radio galaxies are GeV gamma-ray loud. Radio galaxies may contribute to the cosmic MeV gamma-ray background comparable to blazars if gamma-ray-quiet radio galaxies have X-ray to gamma-ray spectra like Cen A, with a small gamma-ray to X-ray flux ratio.

preprint2022arXiv

Interstellar Gas Heating by Primordial Black Holes

Interstellar gas heating is a powerful cosmology-independent observable for exploring the parameter space of primordial black holes (PBHs) formed in the early Universe that could constitute part of the dark matter (DM). We provide a detailed analysis of the various aspects for this observable, such as PBH emission mechanisms. Using observational data from the Leo T dwarf galaxy, we constrain the PBH abundance over a broad mass-range, $M_{\rm PBH} \sim \mathcal{O}(1) M_{\odot}-10^7 M_{\odot}$, relevant for the recently detected gravitational wave signals from intermediate-mass BHs. We also consider PBH gas heating of systems with bulk relative velocity with respect to the DM, such as Galactic clouds.

preprint2022arXiv

Reconstruction of multiple Compton scattering events in MeV gamma-ray Compton telescopes towards GRAMS: the physics-based probabilistic model

Aimed at progress in mega-electron volt (MeV) gamma-ray astronomy, which has not yet been well-explored, Compton telescope missions with a variety of detector concepts have been proposed so far. One of the key techniques for these future missions is an event reconstruction algorithm that is able to determine the scattering orders of multiple Compton scattering events and to identify events in which gamma rays escape from the detectors before they deposit all of their energies. We revisit previous event reconstruction methods and propose a modified algorithm based on a probabilistic method. First, we present a general formalism of the probabilistic model of Compton scattering describing physical interactions inside the detector and measurement processes. Then, we also introduce several approximations in the calculation of the probability functions for efficient computation. For validation, the developed algorithm has been applied to simulation data of a Compton telescope using a liquid argon time projection chamber, which is a new type of Compton telescope proposed for the GRAMS project. We have confirmed that it works successfully for up to 8-hit events, including correction of incoming gamma-ray energies for escape events. The proposed algorithm can be used for next-generation MeV gamma-ray missions featured by large-volume detectors, e.g., GRAMS.

preprint2022arXiv

Spatially resolved study of the SS 433/W50 west region with Chandra: X-ray structure and spectral variation of non-thermal emission

The X-ray binary SS 433, embedded in the W50 nebula (or supernova remnant W50), shows bipolar jets that are ejected with mildly relativistic velocities, and extend toward the east and west out to scales of tens of parsecs. Previous X-ray observations revealed twin lobes along the jet precession axis that contain compact bright knots dominated by synchrotron radiation, which provide evidence of electron acceleration in this system. Particle acceleration in this system is substantiated by the recently detected gamma rays with energies up to at least 25 TeV. To further elucidate the origin of the knots and particle acceleration sites in SS 433/W50, we report here on detailed, spatially resolved X-ray spectroscopy of its western lobe with Chandra. We detect synchrotron emission along the jet precession axis, as well as optically thin thermal emission that is more spatially extended. Between the two previously known knots, w1 and w2, we discover another synchrotron knot, which we call w1.5. We find no significant synchrotron emission between SS 433 and the innermost X-ray knot (w1), suggesting that electrons only begin to be accelerated at w1. The X-ray spectra become gradually steeper from w1 to w2, and then rapidly so immediately outside of w2. Comparing with a model taking into account electron transport and cooling along the jet, this result indicates that the magnetic field in w2 is substantially enhanced, which also explains its brightness. We discuss possible origins of the enhanced magnetic field of w2 as well as scenarios to explain the other two knots.

preprint2021arXiv

Impacts of Jets and Winds From Primordial Black Holes

Primordial black holes (PBHs) formed in the early Universe constitute an attractive candidate for dark matter. Within the gaseous environment of the interstellar medium, PBHs with accretion disks naturally launch outflows such as winds and jets. PBHs with significant spin can sustain powerful relativistic jets and generate associated cocoons. Jets and winds can efficiently deposit their kinetic energies and heat the surrounding gas through shocks. Focusing on the Leo T dwarf galaxy, we demonstrate that these considerations can provide novel tests of PBHs over a significant $\sim 10^{-2} M_{\odot} - 10^6 M_{\odot}$ mass range, including the parameter space associated with gravitational wave observations by the LIGO and VIRGO Collaborations. Observing the morphology of emission could allow to distinguish between jet and wind contributions, and hence indirectly detect spinning PBHs.

preprint2021arXiv

Origin of Galactic Spurs: New Insight from Radio/X-ray All-sky Maps

In this study, we analyze giant Galactic spurs seen in both radio and X-ray all-sky maps to reveal their origins. We discuss two types of giant spurs: one is the brightest diffuse emission near the map's center, which is likely to be related to Fermi bubbles (NPSs/SPSs, north/south polar spurs, respectively), and the other is weaker spurs that coincide positionally with local spiral arms in our Galaxy (LAS, local arm spur). Our analysis finds that the X-ray emissions, not only from the NPS but from the SPS are closer to the Galactic center by ~5 deg compared with the corresponding radio emission. Furthermore, larger offsets of 10-20 deg are observed in the LASs; however, they are attributed to different physical origins. Moreover, the temperature of the X-ray emission is kT ~ 0.2 keV for the LAS, which is systematically lower than those of the NPS and SPS (kT ~ 0.3 keV) but consistent with the typical temperature of Galactic halo gas. We argue that the radio/X-ray offset and the slightly higher temperature of the NPS/SPS X-ray gas are due to the shock compression/heating of halo gas during a significant Galactic explosion in the past, whereas the enhanced X-ray emission from the LAS may be due to the weak condensation of halo gas in the arm potential or star formation activity without shock heating.

preprint2020arXiv

Blazar Radio and Optical Survey (BROS): A catalog of blazar candidates showing flat radio spectrum and their optical identification in Pan-STARRS1 Surveys

Utilizing the latest and the most sensitive radio and optical catalogs, we completed a new blazar candidate catalog, Blazar Radio and Optical Survey (BROS), which includes 88,211 sources located at declination $δ> -40^{\circ}$ and outside the galactic plane ($|b| > 10^{\circ}$). We list compact flat-spectrum radio sources of $α> -0.6$ ($α$ is defined as $F_ν \propto ν^α$ ) from 0.15~GHz TGSS and 1.4~GHz NVSS catalogs. We further identify optical counterparts of the selected sources by cross-matching with Pan-STARRS1 photometric data. Color-color and color-magnitude plots for the selected BROS sources clearly show two distinct populations, An "quasar-like" population consisting of both flat-spectrum radio quasars and BL Lac type objects. On the other hand, an "elliptical-like" population of mostly BL Lac-type objects is buried in the elliptical galaxy. We emphasize that the latter population is missed by previous catalogs but newly emerged in the present BROS catalog, due to the lower radio flux threshold of our selection. Model calculations show that the "elliptical-like" population consists of elliptical galaxies located at redshift z $\lesssim$ 0.5, which is also supported by the logN-logS distribution of the power-law index of $1.49 \pm 0.05$. This BROS catalog is useful for identifying the electromagnetic counterparts of ultra-high-energy cosmic rays and PeV neutrinos recently detected by IceCube, as well as nearby BL Lac objects detectable by future high sensitivity TeV telescopes, such as Cherenkov Telescope Array.

preprint2020arXiv

Cosmological Evolution of Flat-Spectrum Radio Quasars based on the {\it Swift}/BAT 105-month Catalog and Their Contribution to the Cosmic MeV Gamma-ray Background Radiation

We present a new X-ray luminosity function of flat-spectrum radio quasars (FSRQs) utilizing the latest {\it Swift}/BAT 105-month X-ray source catalog. Contrary to previous studies of FSRQs in the X-ray band, using the luminosity-dependent density evolution model, we find that FSRQs show evolutionary peaks at $z\sim1-2$ depending on luminosities. Our result is rather consistent with the evolution of FSRQs seen in the radio and GeV bands, although the number density is a factor of 5--10 smaller. We further explore the contribution of FSRQs to the cosmic MeV gamma-ray background radiation. We find that FSRQs can explain only $\sim3$\% of the observed MeV gamma-ray background fluxes around 1 MeV, indicating other populations are required. Future MeV gamma-ray observations will be keys for understanding the origin of the MeV gamma-ray background radiation.

preprint2020arXiv

Follow-up Observations for IceCube-170922A: Detection of Rapid Near-Infrared Variability and Intensive Monitoring of TXS 0506+056

We present our follow-up observations to search for an electromagnetic counterpart of the IceCube high-energy neutrino, IceCube-170922A. Monitoring observations of a likely counterpart, TXS 0506+056, are also described. First, we quickly took optical and near-infrared images of 7 flat-spectrum radio sources within the IceCube error region right after the neutrino detection and found a rapid flux decline of TXS 0506+056 in Kanata/HONIR J-band data. Motivated by this discovery, intensive follow-up observations of TXS 0506+056 are continuously done, including our monitoring imaging observations, spectroscopic observations, and polarimetric observations in optical and near-infrared wavelengths. TXS 0506+056 shows a large amplitude (~1.0 mag) variability in a time scale of several days or longer, although no significant variability is detected in a time scale of a day or shorter. TXS 0506+056 also shows a bluer-when-brighter trend in optical and near-infrared wavelengths. Structure functions of variabilities are examined and indicate that TXS 0506+056 is not a special blazar in terms of optical variability. Polarization measurement results of TXS 0506+056 are also discussed.

preprint2020arXiv

Multi-wavelength Emission from Galactic Jets: the Case of the Microquasar SS433

SS433 is a Galactic microquasar with powerful jets, where very-high-energy particles are produced. We study particle acceleration in the jets of SS433 in the light of the recent multi-wavelength data from radio to TeV gamma ray. We first present a general framework for the particle acceleration, cooling, and transport in relativistic jets. We then apply this to two X-ray knots in the jets of SS433, focusing on leptonic emission. Our detailed treatment of particle transport and evolution produces substantially different predictions from previous papers. For both regions, our model can account for the multi-wavelength data except for the GeV data. This suggests that GeV emission is mostly from different regions and/or mechanisms. We find that the acceleration process should be efficient, which could be realized by diffusive shock acceleration close to the Bohm limit. Provided that protons are accelerated at the same efficiency as electrons, our results imply that SS433 is a PeVatron, i.e., can accelerate protons beyond a PeV. Future hard X-ray and MeV gamma-ray observations can critically test our models by detecting the spectral turnover or cutoff.

preprint2020arXiv

On the Origin of High Energy Neutrinos from NGC 1068: The Role of Non-Thermal Coronal Activity

NGC 1068, a nearby type-2 Seyfert galaxy, is reported as the hottest neutrino spot in the 10-year survey data of IceCube. Although there are several different possibilities for the generation of high-energy neutrinos in astrophysical sources, feasible scenarios allowing such emission in NGC 1068 have not yet been firmly defined. We show that the flux level of GeV and neutrino emission observed from NGC 1068 implies that the neutrino emission can be produced only in the vicinity of the supermassive black hole in the center of the galaxy. The coronal parameters, such as magnetic field strength and corona size, making this emission possible are consistent with the spectral excess registered in the millimeter range. The suggested model and relevant physical parameters are similar to those revealed for several nearby Seyferts. Due to the internal gamma-ray attenuation, the suggested scenario cannot be verified by observations of NGC 1068 in the GeV and TeV gamma-ray energy bands. However, the optical depth is expected to become negligible for MeV gamma rays, thus future observations in this band will be able to prove our model.

preprint2020arXiv

Snowmass 2021 Letter of Interest: The GRAMS Project: MeV Gamma-Ray Observations and Antimatter-Based Dark Matter Searches

The Gamma-Ray and AntiMatter Survey (GRAMS) project is a proposed next-generation balloon/satellite mission targeting both MeV gamma-ray observations and antimatter-based dark matter searches. A cost-effective, large-scale Liquid Argon Time Projection Chamber (LArTPC) detector technology will allow GRAMS to have a significantly improved sensitivity to MeV gamma rays while extensively probing dark matter parameter space via antimatter measurements.

preprint2016arXiv

Baryon Loading Efficiency and Particle Acceleration Efficiency of Relativistic Jets: Cases For Low Luminosity BL Lacs

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 $ε_{\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 $η_{\rm jet}\sim1.5$ and the mass loading efficiency of $ξ_{\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 $η_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.

preprint2016arXiv

Binary Black Hole Merger Rates Inferred from Luminosity Function of Ultra-Luminous X-ray Sources

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} λ^{-0.6} \exp{(-0.30λ)} \ {\rm Gpc^{-3}\ yr^{-1}}$, where $t_{\rm ULX}$ is the typical duration of the ULX phase and $λ$ 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.

preprint2016arXiv

High-frequency excess in the radio continuum spectrum of the type-1 Seyfert galaxy NGC 985

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 $α=-1.10 \pm 0.03$) changes at ~10 GHz into an inverted spectrum at higher frequencies ($α=+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.

preprint2016arXiv

Lower Bound on the Cosmic TeV Gamma-ray Background Radiation

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.

preprint2016arXiv

The Quiescent Intracluster Medium in the Core of the Perseus Cluster

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.

preprint2015arXiv

Metal Enrichment in the Fermi Bubbles as a Probe of Their Origin

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.

preprint2015arXiv

Prospect for Future MeV Gamma-ray Active Galactic Nuclei Population Studies

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.

preprint2014arXiv

Cosmic Gamma-ray Background Radiation

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.

preprint2014arXiv

Diffuse Neutrino Intensity from the Inner Jets of Active Galactic Nuclei: Impacts of External Photon Fields and the Blazar Sequence

We study high-energy neutrino production in inner jets of radio-loud active galactic nuclei (AGN), taking into account effects of external photon fields and the blazar sequence. We show that the resulting diffuse neutrino intensity is dominated by quasar-hosted blazars, in particular, flat spectrum radio quasars, and that PeV-EeV neutrino production due to photohadronic interactions with broadline and dust radiation is unavoidable if the AGN inner jets are ultrahigh-energy cosmic-ray (UHECR) sources. Their neutrino spectrum has a cutoff feature around PeV energies since target photons are due to Ly$α$ emission. Because of infrared photons provided by the dust torus, neutrino spectra above PeV energies are too hard to be consistent with the IceCube data unless the proton spectral index is steeper than 2.5, or the maximum proton energy is $\lesssim100$ PeV. Thus, the simple model has difficulty in explaining the IceCube data. For the cumulative neutrino intensity from blazars to exceed $\sim{10}^{-8}~{\rm GeV}~{\rm cm}^{-2}~{\rm s}^{-1}~{\rm sr}^{-1}$, their local cosmic-ray energy generation rate would be $\sim10-100$ times larger than the local UHECR emissivity, but is comparable to the averaged gamma-ray blazar emissivity. Interestingly, future detectors such as the Askaryan Radio Array can detect $\sim0.1-1$ EeV neutrinos even in more conservative cases, allowing us to indirectly test the hypothesis that UHECRs are produced in the inner jets. We find that the diffuse neutrino intensity from radio-loud AGN is dominated by blazars with gamma-ray luminosity of $\gtrsim10^{48}~{\rm erg}~{\rm s}^{-1}$, and the arrival directions of their $\sim1-100$ PeV neutrinos correlate with the luminous blazars detected by Fermi.

preprint2014arXiv

Photopion Production in Black-Hole Jets and Flat-Spectrum Radio Quasars as PeV Neutrino Sources

The IceCube collaboration has reported neutrinos with energies between ~30 TeV and a few PeV that are significantly enhanced over the cosmic-ray induced atmospheric background. Viable high-energy neutrino sources must contain very high-energy and ultra-high energy cosmic rays while efficiently making PeV neutrinos. Gamma-ray Bursts (GRBs) and blazars have been considered as candidate cosmic-ray accelerators. GRBs, including low-luminosity GRBs, can be efficient PeV neutrino emitters for low bulk Lorentz factor outflows, although the photopion production efficiency needs to be tuned to simultaneously explain ultra-high-energy cosmic rays. Photopion production efficiency of cosmic-rays accelerated in the inner jets of flat spectrum radio quasars (FSRQs) is ~1-10% due to interactions with photons of the broad-line region (BLR), whereas BL Lac objects are not effective PeV neutrino sources due to the lack of external radiation fields. Photopion threshold effects with BLR photons suppress neutrino production below ~1 PeV, which imply that neutrinos from other sources would dominate over the diffuse neutrino intensity at sub-PeV energies. Reduction of the >> PeV neutrino flux can be expected when curving cosmic-ray proton distributions are employed. Considering a log-parabolic function to describe the cosmic-ray distribution, we discuss possible implications for particle acceleration in black-hole jets. Our results encourage a search for IceCube PeV neutrino events associated with gamma-ray loud FSRQs using Fermi-LAT data. In our model, as found in our previous work, the neutrino flux is suppressed below 1 PeV, which can be tested with increased IceCube exposure.

preprint2014arXiv

The ASTRO-H X-ray Astronomy Satellite

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.

preprint2014arXiv

Unveiling the Nature of Coronae in Active Galactic Nuclei through Sub-mm Observations

The heating mechanism of a corona above an accretion disk in active galactic nuclei (AGNs) is still unknown. One possible mechanism is magnetic reconnection heating requiring energy equipartition between magnetic energy and gas energy in the disk. Here, we investigate the expected observed properties in radio band from such a magnetized corona. A magnetized corona can generate synchrotron radiation since a huge amount of electrons exists. Although most of radiation would be absorbed by synchrotron self-absorption, high-frequency end of synchrotron emission can escape from a corona and appears at the sub-mm range. If only thermal electrons exist in a corona, the expected flux from nearby Seyferts is below the Atacama Large Millimeter/ submillimeter Array (ALMA) sensitivity. However, if non-thermal electrons coexist in a corona, ALMA can measure the non-thermal tail of the synchrotron radiation from a corona. Such non-thermal population is naturally expected to exist if the corona is heated by magnetic reconnections. Future ALMA observations will directly probe the coronal magnetic field strength and the existence of non-thermal electrons in coronae of AGNs.

preprint2014arXiv

Variable optical polarization during high state in gamma-ray loud narrow line Seyfert 1 galaxy 1H 0323+342

We present results of optical polarimetric and multi-band photometric observations for gamma-ray loud narrow-line Seyfert 1 galaxy 1H 0323+342. This object has been monitored by 1.5 m Kanata telescope since 2012 September but following a gamma-ray flux enhancement detected by Fermi-LAT on MJD 56483 (2013 July 10) dense follow-up was performed by ten 0.5-2.0 m telescopes in Japan over one week. The 2-year R_C-band light curve showed clear brightening corresponding to the gamma-ray flux increase and then decayed gradually. The high state as a whole lasted for ~20 days, during which we clearly detected optical polarization from this object. The polarization degree (PD) of the source increased from 0-1% in quiescence to ~3% at maximum and then declined to the quiescent level, with the duration of the enhancement of less than 10 days. The moderate PD around the peak allowed us to precisely measure the daily polarization angle (PA). As a result, we found that the daily PAs were almost constant and aligned to the jet axis, suggesting that the magnetic field direction at the emission region is transverse to the jet. This implies either a presence of helical/toroidal magnetic field or transverse magnetic field compressed by shock(s). We also found small-amplitude intra-night variability during the 2-hour continuous exposure on a single night. We discuss these findings based on the turbulent multi-zone model recently advocated by Marscher (2014). Optical to ultraviolet (UV) spectrum showed a rising shape in the higher frequency and the UV magnitude measured by Swift/UVOT was steady even during the flaring state, suggesting that thermal emission from accretion disk is dominant in that band.

preprint2013arXiv

Probing the Cosmic X-ray and MeV Gamma-ray Background Radiation through the Anisotropy

While the cosmic soft X-ray background is very likely to originate from individual Seyfert galaxies, the origin of the cosmic hard X-ray and MeV gamma-ray background is not fully understood. It is expected that Seyferts including Compton thick population may explain the cosmic hard X-ray background. At MeV energy range, Seyferts having non-thermal electrons in coronae above accretion disks or MeV blazars may explain the background radiation. We propose that future measurements of the angular power spectra of anisotropy of the cosmic X-ray and MeV gamma-ray backgrounds will be key to deciphering these backgrounds and the evolution of active galactic nuclei (AGNs). As AGNs trace the cosmic large-scale structure, spatial clustering of AGNs exists. We show that e-ROSITA will clearly detect the correlation signal of unresolved Seyferts at 0.5-2 keV and 2-10 keV bands and will be able to measure the bias parameter of AGNs at both bands. Once the future hard X-ray all sky satellites achieve the sensitivity better than 10^{-12} erg/cm^2/s at 10-30 keV or 30-50 keV - although this is beyond the sensitivities of current hard X-ray all sky monitors - angular power spectra will allow us to independently investigate the fraction of Compton-thick AGNs in all Seyferts. We also find that the expected angular power spectra of Seyferts and blazars in the MeV range are different by about an order of magnitude, where the Poisson term, so-called shot noise, is dominant. Current and future MeV instruments will clearly disentangle the origin of the MeV gamma-ray background through the angular power spectrum.

preprint2013arXiv

Prospects for future very high-energy gamma-ray sky survey: impact of secondary gamma rays

Very high-energy gamma-ray measurements of distant blazars can be well explained by secondary gamma rays emitted by cascades induced by ultra-high-energy cosmic rays. The secondary gamma rays will enable one to detect a large number of blazars with future ground based gamma-ray telescopes such as Cherenkov Telescope Array (CTA). We show that the secondary emission process will allow CTA to detect 100, 130, 150, 87, and 8 blazars above 30 GeV, 100 GeV, 300 GeV, 1 TeV, and 10 TeV, respectively, up to $z\sim8$ assuming the intergalactic magnetic field (IGMF) strength $B=10^{-17}$ G and an unbiased all sky survey with 0.5 hr exposure at each Field of View, where total observing time is $\sim540$ hr. These numbers will be 79, 96, 110, 63, and 6 up to $z\sim5$ in the case of $B=10^{-15}$ G. This large statistics of sources will be a clear evidence of the secondary gamma-ray scenarios and a new key to studying the IGMF statistically. We also find that a wider and shallower survey is favored to detect more and higher redshift sources even if we take into account secondary gamma rays.

preprint2013arXiv

Revisiting the Cosmic Star Formation History: Caution on the Uncertainties in Dust Correction and Star Formation Rate Conversion

The cosmic star formation rate density (CSFRD) has been observationally investigated out to redshift z~10. However, most of theoretical models for galaxy formation underpredict the CSFRD at z>1. Since the theoretical models reproduce the observed luminosity functions (LFs), luminosity densities (LDs), and stellar mass density at each redshift, this inconsistency does not simply imply that theoretical models should incorporate some missing unknown physical processes in galaxy formation. Here, we examine the cause of this inconsistency in UV wavelengths by using a mock catalog of galaxies generated by a semi-analytic model of galaxy formation. We find that this inconsistency is due to two observational uncertainties: dust obscuration correction and conversion from UV luminosity to star formation rate (SFR). The methods for correction of obscuration and SFR conversion used in observational studies result in the overestimation of CSFRD by ~ 0.1-0.3 dex and ~ 0.1-0.2 dex, respectively, compared to the results obtained directly from our mock catalog. We present new empirical calibrations for dust attenuation and conversion from observed UV LFs and LDs into CSFRD.

preprint2013arXiv

Upper Bound on the First Star Formation History

Our understanding of the nature of the extragalactic background light (EBL) has improved with the recent development of gamma-ray observation techniques. An open subject in the context of the EBL is the reionization epoch, which is an important probe of the formation history of first stars, the so-called Population III (Pop III) stars. Although the mechanisms for the formation of Pop III stars are rather well understood on theoretical grounds, their formation history is still veiled in mystery because of their faintness. To shed light into this matter, we study jointly the gamma-ray opacity of distant objects and the reionization constraints from studies of intergalactic gas. By combining these studies, we obtain a sensitive upper bound on the Pop III star formation rate density as $\dotρ_{*}(z)<0.01[(1+z)/{(1+7.0)}]^{3.4}({f_{\rm esc}}/{0.2})^{-1}({C}/{3.0})\ {\rm M}_{\odot} {\rm yr}^{-1}\ {\rm Mpc}^{-3}$ at $z\ge7$, where $f_{\rm esc}$ and $C$ are the escape fraction of ionizing photons from galaxies and the clumping factor of the intergalactic hydrogen gas. This limit is a $\sim10$ times tighter constraint compared with previous studies that take into account gamma-ray opacity constraints only. Even if we do not include the current gamma-ray constraints, the results do not change. This is because the detected gamma-ray sources are still at $z\le4.35$ where the reionization has already finished.

preprint2012arXiv

The ASTRO-H X-ray Observatory

The joint JAXA/NASA ASTRO-H mission is the sixth in a series of highly successful X-ray missions initiated by the Institute of Space and Astronautical Science (ISAS). ASTRO-H will investigate the physics of the high-energy universe via a suite of four instruments, covering a very wide energy range, from 0.3 keV to 600 keV. These instruments include a high-resolution, high-throughput spectrometer sensitive over 0.3-2 keV with high spectral resolution of Delta E < 7 eV, enabled by a micro-calorimeter array located in the focal plane of thin-foil X-ray optics; hard X-ray imaging spectrometers covering 5-80 keV, located in the focal plane of multilayer-coated, focusing hard X-ray mirrors; a wide-field imaging spectrometer sensitive over 0.4-12 keV, with an X-ray CCD camera in the focal plane of a soft X-ray telescope; and a non-focusing Compton-camera type soft gamma-ray detector, sensitive in the 40-600 keV band. The simultaneous broad bandpass, coupled with high spectral resolution, will enable the pursuit of a wide variety of important science themes.

preprint2012arXiv

Upper Limit on the Cosmological Gamma-ray Background

We show that the current extragalactic gamma-ray background (EGB) measurement below 100 GeV sets an upper limit on EGB itself at very high energy (VHE) above 100 GeV. The limit is conservative for the electromagnetic cascade emission from VHE EGB interacting with the cosmic microwave-to-optical background radiation not to exceed the current EGB measurement. The cascade component fits the measured VHE EGB spectrum rather well. However, once we add the contribution from known source classes, the Fermi VHE EGB observation exceeds or even violates the limit, which is approximated as E^2dN/dE < 4.5x10^-5 (E/100 GeV)^-0.7 MeV/cm^2/s/sr. The upper limit above 100 GeV is useful in the future to probe the EGB origin and the new physics like axion-like particles and Lorentz-invariance violation.

preprint2011arXiv

Contribution of the Gamma-ray Loud Radio Galaxies Core Emissions to the Cosmic MeV and GeV Gamma-Ray Background Radiation

The Fermi gamma-ray satellite has recently detected gamma-ray emissions from radio galaxy cores. From these samples, we first examine the correlation between the luminosities at 5 GHz, L_{5GHz}, and at 0.1-10 GeV, L_{gamma}, of these gamma-ray loud radio galaxies. We find that the correlation is significant with L_{gamma} \propto L_{5GHz}^{1.16} based on a partial correlation analysis. Using this correlation and the radio luminosity function (RLF) of radio galaxies, we further explore the contribution of gamma-ray loud radio galaxies to the unresolved extragalactic gamma-ray background (EGRB). The gamma-ray luminosity function is obtained by normalizing the RLF to reproduce the source count distribution of the Fermi gamma-ray loud radio galaxies. We find that gamma-ray loud radio galaxies will explain ~25% of the unresolved Fermi EGRB flux above 100 MeV and will also make a significant contribution to the EGRB in the 1-30 MeV energy band. Since blazars explain 22% of the EGRB above 100 MeV, radio loud active galactic nuclei (AGNs) population explains ~47% of the unresolved EGRB. We further make an interpretation on the origin of the EGRB. The observed EGRB spectrum at 0.2-100 GeV does not show an absorption signature by the extragalactic background light. Thus, the dominant population of the origin of EGRB at very high energy (>30 GeV) might be nearby gamma-ray emitting sources or sources with very hard gamma-ray spectrum.

preprint2011arXiv

Extragalactic Gamma-ray Background Radiation from Beamed and Unbeamed Active Galactic Nuclei

The origin of the extragalactic gamma-ray background (EGRB) radiation has been a mystery in astrophysics for a long time. Recently the Fermi gamma-ray satellite (Fermi) has revealed that ~22% of the unresolved EGRB would be explained by blazars, which are one population of beamed active galactic nuclei (AGNs). The remaining ~78% of the unresolved EGRB is still unknown. We estimate the contribution of gamma-ray loud radio galaxies, which are misaligned radio loud AGNs recently detected by Fermi, to EGRB using the radio luminosity function of radio-loud AGNs with the correlation between the radio and gamma-ray luminosities. We find that ~25% of the unresolved EGRB will be explained by gamma-ray loud radio galaxy population. We also discuss further about the origin of EGRB by comparing the Fermi EGRB data with our studies on various AGN populations' contribution to EGRB, which are radio quiet AGNs, blazars, and gamma-ray loud radio galaxies.

preprint2010arXiv

High Energy Gamma-ray Absorption and Cascade Emission in Nearby Starburst Galaxies

High energy gamma-ray emission from two nearby bright starburst galaxies, M82 and NGC 253, have recently been detected by Fermi, H.E.S.S., and VERITAS. Since starburst galaxies have a high star formation rate and plenty of dust in the central starburst region, infrared emissions are strong there. Gamma-ray photons are absorbed by the interstellar radiation field photons via electron and positron pair creation. The generated electron and positron pairs up scatter the interstellar photons to very high energy gamma-ray photons via cascade emission through inverse Compton scattering. In this paper, we evaluate the contribution of this cascade emission to the gamma-ray spectra of M82 and NGC 253. Although it would be difficult to see direct gamma- ray evidence of cosmic-rays with an energy > 10 TeV due to the gamma-ray attenuation, the resulting cascade emission would be indirect evidence. By including the cascade component, we find that the total flux above 1 TeV increases ~18% and ~45% compared with the absorbed flux assuming the maximum kinetic proton energy as 45.3 TeV and 512 TeV, respectively. Future gamma-ray observatories such as CTA would be able to see the indirect evidence of cosmic-ray with an energy > 10 TeV by comparing with theoretical emission models including this cascade effect.

preprint2010arXiv

Prospects for Very High Energy Blazar Survey by the Next Generation Cherenkov Telescopes

The prospects for future blazar surveys by next-generation very high energy (VHE) gamma-ray telescopes, such as Advanced Gamma-ray Imaging System (AGIS) and Cherenkov Telescope Array (CTA), are investigated using the latest model of blazar luminosity function and its evolution which is in good agreement with the flux and redshift distribution of observed blazars as well as the extragalactic gamma-ray background. We extend and improve the template of spectral energy distributions (SEDs) based on the blazar SED sequence paradigm, to make it reliable also in the VHE bands (above 100 GeV) by comparing with the existing VHE blazar data. Assuming the planned CTA sensitivities, a blind survey using a total survey time of ~100 hrs could detect ~3 VHE blazars, with larger expected numbers for wider/shallower surveys. We also discuss following-up of Fermi blazars. Detectability of VHE blazars in the plane of Fermi flux and redshift is presented, which would be useful for future survey planning. Prospects and strategies are discussed to constrain the extragalactic background light (EBL) by using the absorption feature of brightest blazar spectra, as well as cut-offs in the redshift distribution. We will be able to get useful constraints on EBL by VHE blazars at different redshifts ranging 0.3-1 TeV corresponding to z=0.10-0.36.

preprint2010arXiv

Searching for the most distant blazars with the Fermi Gamma-ray Space Telescope

We investigate the prospects for discovering blazars at very high-redshifts (z>3-6) with the Fermi Gamma-Ray Space Telescope (Fermi), employing a model for the evolving gamma-ray luminosity function (GLF) of the blazar population. Our previous GLF model is used as a basis, which features luminosity-dependent density evolution implied from X-ray data on active galactic nuclei, as well as the blazar sequence paradigm for their spectral energy distribution, and which is consistent with EGRET and current Fermi observations of blazars.Here we augment the high-redshift evolution of this model by utilizing the luminosity function of quasars from the Sloan Digital Sky Survey (SDSS), which is well-constrained up to z~5. We find that Fermi may discover a few blazars up to z~6 in the entire sky during its 5-year survey. We further discuss how such high-redshift blazar candidates may be efficiently selected in future Fermi data.

preprint2009arXiv

The Blazar Sequence and the Cosmic Gamma-Ray Background Radiation in the Fermi Era

(Abridged) We present a new model of the blazar gamma-ray luminosity function (GLF) and the spectrum of the extragalactic gamma-ray background (EGRB), which is consistent with the observed distributions of EGRET blazars. The unified sequence of blazar spectral energy distribution (SED) is taken into account to make a non-trivial prediction for the EGRB spectrum and more realistic comparison with the data than previous studies. We then try to explain the EGRB data by the two AGN populations: one is blazars, and the other is non-blazar AGNs that are responsible for the EGRB in the MeV band. We find that ~80% of the EGRB photon flux at > 100 MeV can be explained by the sum of the two populations, while ~45% can be accounted for only by blazars. The predicted EGRB spectrum is in agreement with a wide range of the observed data from X-ray to GeV, within the systematic uncertainties in the EGRB determination by EGRET. These results indicate that AGNs including blazars are the primary source of EGRB. Blazars are dominant in EGRB at higher energy bands of >100 MeV, while non-blazar AGNs dominate at <100 MeV. Almost all of the EGRB flux from blazars will be resolved into discrete sources by the Fermi Gamma-ray Space Telescope, while that from non-blazar AGNs will largely remain unresolved. Therefore, comparison between the integrated source counts and diffuse EGRB flux as a function of photon energy will give a simple and clear test of our model. Various quantitative predictions for Fermi observations are also made. Especially, our model predicts 600-1200 blazars in all sky down to 2e-9 photons/cm^-2/s (>100 MeV), which is considerably smaller than most of previous studies.

preprint2009arXiv

The Cosmological Evolution of Blazars and the Extragalactic Gamma-Ray Background in the Fermi Era

The latest determination of the extragalactic gamma-ray background (EGRB) radiation by Fermi is compared with the theoretical prediction of the blazar component by Inoue & Totani (2009; hereafter IT09). The Fermi EGRB spectrum is in excellent agreement with IT09, indicating that blazars are the dominant component of the EGRB, and contributions from any other sources (e.g., dark matter annihilations) are minor. It also indicates that the blazar SED (spectral energy distribution) sequence taken into account in IT09 is a valid description of mean blazar SEDs. The possible contribution of MeV blazars to the EGRB in the MeV band is also discussed. In five total years of observations, we predict that Fermi will detect ~1200 blazars all sky down to the corresponding sensitivity limit. We also address the detectability of the highest-redshift blazars. Updating our model with regard to high-redshift evolution based on SDSS quasar data, we show that Fermi may find some blazars up to z~6 during the five-year survey. Such blazars could provide a new probe of early star and galaxy formation through GeV spectral attenuation signatures induced by high-redshift UV background radiation.