Source author record

Tomonori Totani

Tomonori Totani 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

45works
12topics
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

45 published item(s)

preprint2022arXiv

On the use of CHIME to Detect Long-Duration Radio Transients from Neutron Star Mergers

Short gamma-ray burst (SGRB) GRB 170817A was found to be related to a binary neutron star (BNS) merger. It is uncertain whether all SGRBs are caused by BNS mergers, and also under what conditions a BNS merger can cause a SGRB. As BNS mergers can cause SGRBs, afterglow observations will also provide an alternative measurement of the BNS merger rate independent of gravitational-wave observations. In previous work by Feng et al. (2014), the feasibility of the detection of afterglows was considered using a variety of radio observatories and a simple flux threshold detection algorithm. Here, we consider a more sophisticated detection algorithm for SGRB afterglows, and provide an estimate of the trials factors for a realistic search to obtain an updated estimate of the possibility of observing afterglows with the Canadian Hydrogen Intensity Mapping Experiment (CHIME). We estimate 893 and 312 afterglows per year can be detected using a $3 σ$ confidence level threshold with two jet models, one with half opening angle uniformly distributed in $6^\circ$ to $30^\circ$ and the other uniformly distributed in $3^\circ$ to $8^\circ$ with the median $6^\circ$. We also find 88% and 98%, respectively, of the detectable afterglows for each jet opening distribution are off-axis, which are candidates for orphan afterglows. Our result predicts fewer detectable sources per year than the earlier analysis, but confirms the essential conclusion that using CHIME to search for afterglows will be effective at constraining the astrophysical merger rate.

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

Detectability of radio afterglows from binary neutron star mergers and implications for fast radio bursts

Binary neutron star (BNS) mergers are one of the proposed origins for both repeating and non-repeating fast radio bursts (FRBs), which associates FRBs with gravitational waves and short gamma-ray bursts (GRBs). In this work, we explore detectability of radio afterglows from BNS mergers and compare it to the observed radio limits on FRB afterglows. We calculate the afterglow flux powered by the two components: a relativistic jet and a slower isotropic ejecta, and quantify the detection probability as a function of the source redshift, observing time, and flux sensitivity. The model parameter distributions inferred from short GRB afterglows are adopted, and viewing angle distributions (uniform spherical, gravitational-wave, on-axis biased) are assumed to reflect different searching scenario. Assuming that FRBs are not strongly beamed, we make comparison to FRBs detected with reported radio limits and find the detection probabilities are 1--10% in general, and hence not a strong constraint on the BNS progenitor model considering the small sample number (<10). In particular for some nearby FRBs (e.g. 180916.J0158+65, 190608), we find a high chance of detection (>20% at 10$μ$Jy sensitivity) for the isotropic component that would peak around $\sim$1--10 years after the merger. Therefore a long term radio monitoring of persistent radio emission for these objects is important.

preprint2020arXiv

Emergence of life in an inflationary universe

Abiotic emergence of ordered information stored in the form of RNA is an important unresolved problem concerning the origin of life. A polymer longer than 40--100 nucleotides is necessary to expect a self-replicating activity, but the formation of such a long polymer having a correct nucleotide sequence by random reactions seems statistically unlikely. However, our universe, created by a single inflation event, likely includes more than $10^{100}$ Sun-like stars. If life can emerge at least once in such a large volume, it is not in contradiction with our observations of life on Earth, even if the expected number of abiogenesis events is negligibly small within the observable universe that contains only $10^{22}$ stars. Here, a quantitative relation is derived between the minimum RNA length $l_{\min}$ required to be the first biological polymer, and the universe size necessary to expect the formation of such a long and active RNA by randomly adding monomers. It is then shown that an active RNA can indeed be produced somewhere in an inflationary universe, giving a solution to the abiotic polymerization problem. On the other hand, $l_{\min}$ must be shorter than $\sim$20 nucleotides for the abiogenesis probability close to unity on a terrestrial planet, but a self-replicating activity is not expected for such a short RNA. Therefore, if extraterrestrial organisms of a different origin from those on Earth are discovered in the future, it would imply an unknown mechanism at work to polymerize nucleotides much faster than random statistical processes.

preprint2019arXiv

The Galactic Halo Contribution to the Dispersion Measure of Extragalactic Fast Radio Bursts

A new model of the Milky Way (MW) halo component of the dispersion measure (DM) for extragalactic sources, such as fast radio bursts (FRBs), is presented in light of recent diffuse X-ray observations. In addition to the spherical component of isothermal gas ($kT\sim0.3$ keV) in hydrostatic equilibrium with the Galactic gravitational potential, our model includes a disk-like non-spherical hot gas component to reproduce the directional dependence of the observed X-ray emission measure (EM). The total gas mass ($1.2\times10^{11}\,M_{\odot}$) is dominated by the spherical component, and is consistent with the total baryon mass of the MW expected from the dark matter mass and the cosmic baryon-to-dark-matter ratio. Our model predicts a mean halo DM of $43\:\,{\rm pc\:cm^{-3}}$, with a full range of $30$-$245\:\,{\rm pc\:cm^{-3}}$ over the whole sky. The large scatter seen in the X-ray EM data implies a $\sim0.2$ dex (rms) fluctuation of the MW halo DM. We provide an analytic formula to estimate the MW halo DM of our model along any line of sight, which can be easily used to compute the total MW component of DM toward extragalactic sources, in combination with existing DM models of the warm ionized medium associated with the Galactic disk.

preprint2016arXiv

Testing anthropic reasoning for the cosmological constant with a realistic galaxy formation model

The anthropic principle is one of the possible explanations for the cosmological constant ($Λ$) problem. In previous studies, a dark halo mass threshold comparable with our Galaxy must be assumed in galaxy formation to get a reasonably large probability of finding the observed small value, $P(<$$Λ_{\rm obs})$, though stars are found in much smaller galaxies as well. Here we examine the anthropic argument by using a semi-analytic model of cosmological galaxy formation, which can reproduce many observations such as galaxy luminosity functions. We calculate the probability distribution of $Λ$ by running the model code for a wide range of $Λ$, while other cosmological parameters and model parameters for baryonic processes of galaxy formation are kept constant. Assuming that the prior probability distribution is flat per unit $Λ$, and that the number of observers is proportional to stellar mass, we find $P(<$$Λ_{\rm obs}) = 6.7 \%$ without introducing any galaxy mass threshold. We also investigate the effect of metallicity; we find $P(<$$Λ_{\rm obs}) = 9.0 \%$ if observers exist only in galaxies whose metallicity is higher than the solar abundance. If the number of observers is proportional to metallicity, we find $P(<$$Λ_{\rm obs}) = 9.7 \%$. Since these probabilities are not extremely small, we conclude that the anthropic argument is a viable explanation, if the value of $Λ$ observed in our universe is determined by a probability distribution.

preprint2016arXiv

The Subaru FMOS Galaxy Redshift Survey (FastSound). II. The Emission Line Catalog and Properties of Emission Line Galaxies

We present basic properties of $\sim$3,300 emission line galaxies detected by the FastSound survey, which are mostly H$α$ emitters at $z \sim$ 1.2-1.5 in the total area of about 20 deg$^2$, with the H$α$ flux sensitivity limit of $\sim 1.6 \times 10^{-16} \rm erg \ cm^{-2} s^{-1}$ at 4.5 sigma. This paper presents the catalogs of the FastSound emission lines and galaxies, which will be open to the public in the near future. We also present basic properties of typical FastSound H$α$ emitters, which have H$α$ luminosities of $10^{41.8}$-$10^{43.3}$ erg/s, SFRs of 20--500 $M_\odot$/yr, and stellar masses of $10^{10.0}$--$10^{11.3}$ $M_\odot$. The 3D distribution maps for the four fields of CFHTLS W1--4 are presented, clearly showing large scale clustering of galaxies at the scale of $\sim$ 100--600 comoving Mpc. Based on 1,105 galaxies with detections of multiple emission lines, we estimate that contamination of non-H$α$ lines is about 4% in the single-line emission galaxies, which are mostly [OIII]$λ$5007. This contamination fraction is also confirmed by the stacked spectrum of all the FastSound spectra, in which H$α$, [NII]$λλ$6548,6583, [SII]$λλ$6717, 6731, and [OI]$λλ$6300,6364 are seen.

preprint2016arXiv

The Subaru FMOS galaxy redshift survey (FastSound). IV. New constraint on gravity theory from redshift space distortions at $z\sim 1.4$

We measure the redshift-space correlation function from a spectroscopic sample of 2783 emission line galaxies from the FastSound survey. The survey, which uses the Subaru Telescope and covers the redshift ranges of $1.19<z<1.55$, is the first cosmological study at such high redshifts. We detect clear anisotropy due to redshift-space distortions (RSD) both in the correlation function as a function of separations parallel and perpendicular to the line of sight and its quadrupole moment. RSD has been extensively used to test general relativity on cosmological scales at $z<1$. Adopting a LCDM cosmology with the fixed expansion history and no velocity dispersion $σ_{\rm v}=0$, and using the RSD measurements on scales above 8Mpc/h, we obtain the first constraint on the growth rate at the redshift, $f(z)σ_8(z)=0.482\pm 0.116$ at $z\sim 1.4$ after marginalizing over the galaxy bias parameter $b(z)σ_8(z)$. This corresponds to $4.2σ$ detection of RSD. Our constraint is consistent with the prediction of general relativity $fσ_8\sim 0.392$ within the $1-σ$ confidence level. When we allow $σ_{\rm v}$ to vary and marginalize it over, the growth rate constraint becomes $fσ_8=0.494^{+0.126}_{-0.120}$. We also demonstrate that by combining with the low-z constraints on $fσ_8$, high-z galaxy surveys like the FastSound can be useful to distinguish modified gravity models without relying on CMB anisotropy experiments.

preprint2015arXiv

High Precision Analyses of Lyman alpha Damping Wing of Gamma-Ray Bursts in the Reionization Era: On the Controversial Results from GRB 130606A at z = 5.91

The unprecedentedly bright afterglow of Swift GRB 130606A at z = 5.91 gave us a unique opportunity to probe the reionization era by high precision analyses of the redward damping wing of Ly alpha absorption, but the reported constraints on the neutral hydrogen fraction (f_HI) in intergalactic medium (IGM) derived from spectra taken by different telescopes are in contradiction. Here we examine the origin of this discrepancy by analyzing the spectrum taken by VLT with our own analysis code previously used to fit the Subaru spectrum. Though the VLT team reported no evidence for IGM HI using the VLT spectrum, we confirmed our previous result of preferring non-zero IGM HI (the best-fit f_HI ~ 0.06, when IGM HI extends to the GRB redshift). The fit residuals of the VLT spectrum by the model without IGM HI show the same systematic trend as the Subaru spectrum. We consider that the likely origin of the discrepancy between the two teams is the difference of the wavelength ranges adopted in the fittings; our wavelength range is wider than that of the VLT team, and also we avoided the shortest wavelength range of deep Ly alpha absorption (lambda_obs < 8426 A), because this region is dominated by HI in the host galaxy and the systematic uncertainty about host HI velocity distribution is large. We also study the sensitivity of these results to the adopted Ly alpha cross section formulae, ranging from the classical Lorentzian function to the most recent one taking into account fully quantum mechanical scattering. It is found that the preference for non-zero IGM HI is robust against the choice of the cross section formulae, but it is quantitatively not negligible and hence one should be careful in future analyses.

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

The Subaru FMOS Galaxy Redshift Survey (FastSound). I. Overview of the Survey Targeting on H$α$ Emitters at $z \sim 1.4$

FastSound is a galaxy redshift survey using the near-infrared Fiber Multi-Object Spectrograph (FMOS) mounted on the Subaru Telescope, targeting H$α$ emitters at $z \sim 1.18$--$1.54$ down to the sensitivity limit of H$α$ flux $\sim 2 \times 10^{-16} \ \rm erg \ cm^{-2} s^{-1}$. The primary goal of the survey is to detect redshift space distortions (RSD), to test General Relativity by measuring the growth rate of large scale structure and to constrain modified gravity models for the origin of the accelerated expansion of the universe. The target galaxies were selected based on photometric redshifts and H$α$ flux estimates calculated by fitting spectral energy distribution (SED) models to the five optical magnitudes of the Canada France Hawaii Telescope Legacy Survey (CFHTLS) Wide catalog. The survey started in March 2012, and all the observations were completed in July 2014. In total, we achieved $121$ pointings of FMOS (each pointing has a $30$ arcmin diameter circular footprint) covering $20.6$ deg$^2$ by tiling the four fields of the CFHTLS Wide in a hexagonal pattern. Emission lines were detected from $\sim 4,000$ star forming galaxies by an automatic line detection algorithm applied to 2D spectral images. This is the first in a series of papers based on FastSound data, and we describe the details of the survey design, target selection, observations, data reduction, and emission line detections.

preprint2015arXiv

The Subaru FMOS Galaxy Redshift Survey (FastSound). III. The mass-metallicity relation and the fundamental metallicity relation at $z\sim1.4$

We present the results from a large near-infrared spectroscopic survey with Subaru/FMOS (\textit{FastSound}) consisting of $\sim$ 4,000 galaxies at $z\sim1.4$ with significant H$α$ detection. We measure the gas-phase metallicity from the [N~{\sc ii}]$λ$6583/H$α$ emission line ratio of the composite spectra in various stellar mass and star-formation rate bins. The resulting mass-metallicity relation generally agrees with previous studies obtained in a similar redshift range to that of our sample. No clear dependence of the mass-metallicity relation with star-formation rate is found. Our result at $z\sim1.4$ is roughly in agreement with the fundamental metallicity relation at $z\sim0.1$ with fiber aperture corrected star-formation rate. We detect significant [S~{\sc ii}]$λλ$6716,6731 emission lines from the composite spectra. The electron density estimated from the [S~{\sc ii}]$λλ$6716,6731 line ratio ranges from 10 -- 500 cm$^{-3}$, which generally agrees with that of local galaxies. On the other hand, the distribution of our sample on [N~{\sc ii}]$λ$6583/H$α$ vs. [S~{\sc ii}]$λλ$6716,6731/H$α$ is different from that found locally. We estimate the nitrogen-to-oxygen abundance ratio (N/O) from the N2S2 index, and find that the N/O in galaxies at $z\sim1.4$ is significantly higher than the local values at a fixed metallicity and stellar mass. The metallicity at $z\sim1.4$ recalculated with this N/O enhancement taken into account decreases by 0.1 -- 0.2 dex. The resulting metallicity is lower than the local fundamental metallicity relation.

preprint2014arXiv

FIELD: An automated emission-line detection software for Subaru/FMOS near-infrared spectroscopy

We describe the development of automated emission line detection software for the Fiber Multi-Object Spectrograph (FMOS), which is a near-infrared spectrograph fed by $400$ fibers from the $0.2$ deg$^2$ prime focus field of view of the Subaru Telescope. The software, FIELD (FMOS software for Image-based Emission Line Detection), is developed and tested mainly for the FastSound survey, which is targeting H$α$ emitting galaxies at $z \sim 1.3$ to measure the redshift space distortion as a test of general relativity beyond $z \sim 1$. The basic algorithm is to calculate the line signal-to-noise ratio ($S/N$) along the wavelength direction, given by a 2-D convolution of the spectral image and a detection kernel representing a typical emission line profile. A unique feature of FMOS is its use of OH airglow suppression masks, requiring the use of flat-field images to suppress noise around the mask regions. Bad pixels on the detectors and pixels affected by cosmic-rays are efficiently removed by using the information obtained from the FMOS analysis pipeline. We limit the range of acceptable line-shape parameters for the detected candidates to further improve the reliability of line detection. The final performance of line detection is tested using a subset of the FastSound data; the false detection rate of spurious objects is examined by using inverted frames obtained by exchanging object and sky frames. The false detection rate is $< 1$\% at $S/N > 5$, allowing an efficient and objective emission line search for FMOS data at the line flux level of $\gtrsim 1.0 \times 10^{-16}$[erg/cm$^2$/s].

preprint2014arXiv

Galaxy luminosity function and its cosmological evolution: Testing a new feedback model depending on galaxy-scale dust opacity

We present a new version of a semi-analytic model of cosmological galaxy formation, incorporating a star formation law with a feedback depending on the galaxy-scale mean dust opacity and metallicity, motivated by recent observations of star formation in nearby galaxies and theoretical considerations. This new model is used to investigate the effect of such a feedback on shaping the galaxy luminosity function and its evolution. Star formation activity is significantly suppressed in dwarf galaxies by the new feedback effect, and the faint-end slope of local luminosity functions can be reproduced with a reasonable strength of supernova feedback, which is in contrast to the previous models that require a rather extreme strength of supernova feedback. Our model can also reproduce the early appearance of massive galaxies manifested in the bright-end of high redshift K-band luminosity functions. Though some of the previous models also succeeded in reproducing this, they assumed a star formation law depending on the galaxy-scale dynamical time, which is not supported by observations. We argue that the feedback depending on dust opacity (or metal column density) is essential, rather than that simply depending on gas column density, to get these results.

preprint2014arXiv

On the Systematic Errors of Cosmological-Scale Gravity Tests using Redshift Space Distortion: Non-linear Effects and the Halo Bias

Redshift space distortion (RSD) observed in galaxy redshift surveys is a powerful tool to test gravity theories on cosmological scales, but the systematic uncertainties must carefully be examined for future surveys with large statistics. Here we employ various analytic models of RSD and estimate the systematic errors on measurements of the structure growth-rate parameter, $fσ_8$, induced by non-linear effects and the halo bias with respect to the dark matter distribution, by using halo catalogues from 40 realisations of $3.4 \times 10^8$ comoving $h^{-3}$Mpc$^3$ cosmological N-body simulations. We consider hypothetical redshift surveys at redshifts z=0.5, 1.35 and 2, and different minimum halo mass thresholds in the range of $5.0 \times 10^{11}$ -- $2.0 \times 10^{13} h^{-1} M_\odot$. We find that the systematic error of $fσ_8$ is greatly reduced to ~5 per cent level, when a recently proposed analytical formula of RSD that takes into account the higher-order coupling between the density and velocity fields is adopted, with a scale-dependent parametric bias model. Dependence of the systematic error on the halo mass, the redshift, and the maximum wavenumber used in the analysis is discussed. We also find that the Wilson-Hilferty transformation is useful to improve the accuracy of likelihood analysis when only a small number of modes are available in power spectrum measurements.

preprint2014arXiv

Probing Intergalactic Neutral Hydrogen by the Lyman Alpha Red Damping Wing of Gamma-Ray Burst 130606A Afterglow Spectrum at z = 5.913

The unprecedentedly bright optical afterglow of GRB 130606A located by Swift at a redshift close to the reionization era (z = 5.913) provides a new opportunity to probe the ionization status of intergalactic medium (IGM). Here we present an analysis of the red Ly alpha damping wing of the afterglow spectrum taken by Subaru/FOCAS during 10.4-13.2 hr after the burst. We find that the minimal model including only the baseline power-law and HI absorption in the host galaxy does not give a good fit, leaving residuals showing concave curvature in 8400-8900 A with an amplitude of about 0.6% of the flux. Such a curvature in the short wavelength range cannot be explained either by extinction at the host with standard extinction curves, intrinsic curvature of afterglow spectra, or by the known systematic uncertainties in the observed spectrum. The red damping wing by intervening HI gas outside the host can reduce the residual by about 3 sigma statistical significance. We find that a damped Ly alpha system is not favored as the origin of this intervening HI absorption, from the observed Ly beta and metal absorption features. Therefore absorption by diffuse IGM remains as a plausible explanation. A fit by a simple uniform IGM model requires HI neutral fraction of f_HI ~ 0.1-0.5 depending on the distance to the GRB host, implying high f_HI IGM associated with the observed dark Gunn-Peterson (GP) troughs. This gives a new evidence that the reionization is not yet complete at z ~ 6.

preprint2014arXiv

Unveiling the Origin of Fast Radio Bursts by Optical Follow Up Observations

We discuss how we can detect and identify counterparts of fast radio bursts (FRBs) in future optical follow up observations of FRBs if real-time alert of FRBs becomes available. We consider kilonovae as candidates of FRB optical counterparts, as expected in the case that FRBs originate from mergers of double neutron star binaries. Although theoretical predictions on luminosities of kilonovae are still highly uncertain, recent models suggest that kilonovae can be detected at redshifts up to z $\sim$ 0.3 within the range of the uncertainties. We expect $\sim$ 1--5 unrelated supernovae (SNe) down to a similar variability magnitude in 5 days interval within the typical error radius of a FRB. We show that, however, a kilonova can be distinguished from these SNe by its rapid decay and/or color evolution, making it possible to verify the existence of a kilonova associated with a FRB. We also discuss the case that SNe Ia are FRB optical counterparts, as it might be if FRBs originate from double white dwarf binaries. Verification of this scenario is also possible, since the chance probability of finding a SNe Ia having consistent explosion time with that of a FRB within the FRB error region is small (typically $\lesssim$ 0.01).

preprint2013arXiv

A Study of Selection Methods for H alpha Emitting Galaxies at z~1.3 for the Subaru/FMOS Galaxy Redshift Survey for Cosmology (FastSound)

The efficient selection of high-redshift emission galaxies is important for future large galaxy redshift surveys for cosmology. Here we describe the target selection methods for the FastSound project, a redshift survey for H alpha emitting galaxies at z=1.2-1.5 using Subaru/FMOS to measure the linear growth rate fσ8 via Redshift Space Distortion (RSD) and constrain the theory of gravity. To select ~400 target galaxies in the 0.2 deg^2 FMOS field-of-view from photometric data of CFHTLS-Wide (u*g'r'i'z'), we test several different methods based on color-color diagrams or photometric redshift estimates from spectral energy distribution (SED) fitting. We also test the improvement in selection efficiency that can be achieved by adding near-infrared data from the UKIDSS DXS (J). The success rates of H alpha detection with FMOS averaged over two observed fields using these methods are 11.3% (color-color, optical), 13.6% (color-color, optical+NIR), 17.3% (photo-z, optical), and 15.1% (photo-z, optical+NIR). Selection from photometric redshifts tends to give a better efficiency than color-based methods, although there is no significant improvement by adding J band data within the statistical scatter. We also investigate the main limiting factors for the success rate, by using the sample of the HiZELS H alpha emitters that were selected by narrow-band imaging. Although the number density of total H alpha emitters having higher H alpha fluxes than the FMOS sensitivity is comparable with the FMOS fiber density, the limited accuracy of photometric redshift and H alpha flux estimations have comparable effects on the success rate of <~20% obtained from SED fitting.

preprint2013arXiv

Constraints on f(R) theory and Galileons from the latest data of galaxy redshift surveys

The growth rate of matter density perturbations has been measured from redshift-space distortion (RSD) in the galaxy power spectrum. We constrain the model parameter space for representative modified gravity models to explain the dark energy problem, by using the recent data of f_m(z)sigma_8(z) at the redshifts z = 0.06--0.8 measured by WiggleZ, SDSS LRG, BOSS, and 6dFGRS. We first test the Hu-Sawicki's f(R) dark energy model, and find that only the parameter region close to the standard Lambda Cold Dark Matter (Lambda-CDM) model is allowed (lambda > 12 and 5 for n = 1.5 and 2, respectively, at 95% CL). We then investigate the covariant Galileon model and show that the parameter space consistent with the background expansion history is excluded by the RSD data at more than 10 sigma because of the too large growth rate predicted by the theory. Finally, we consider the extended Galileon scenario, and we find that, in contrast to the covariant Galileon, there is a model parameter space for a tracker solution that is consistent with the RSD data within a 2 sigma level.

preprint2013arXiv

Cosmological Fast Radio Bursts from Binary Neutron Star Mergers

Fast radio bursts (FRBs) at cosmological distances have recently been discovered, whose duration is about milliseconds. We argue that the observed short duration is difficult to explain by giant flares of soft gamma-ray repeaters, though their event rate and energetics are consistent with FRBs. Here we discuss binary neutron star (NS-NS) mergers as a possible origin of FRBs. The FRB rate is within the plausible range of NS-NS merger rate and its cosmological evolution, while a large fraction of NS-NS mergers must produce observable FRBs. A likely radiation mechanism is coherent radio emission like radio pulsars, by magnetic braking when magnetic fields of neutron stars are synchronized to binary rotation at the time of coalescence. Magnetic fields of the standard strength (~ 10^{12-13} G) can explain the observed FRB fluxes, if the conversion efficiency from magnetic braking energy loss to radio emission is similar to that of isolated radio pulsars. Corresponding gamma-ray emission is difficult to detect by current or past gamma-ray burst satellites. Since FRBs tell us the exact time of mergers, a correlated search would significantly improve the effective sensitivity of gravitational wave detectors.

preprint2013arXiv

Supernova Neutrino Light Curves and Spectra for Various Progenitor Stars: From Core Collapse to Proto-neutron Star Cooling

We present a new series of supernova neutrino light curves and spectra calculated by numerical simulations for a variety of progenitor stellar masses (13-50Msolar) and metallicities (Z = 0.02 and 0.004), which would be useful for a broad range of supernova neutrino studies, e.g., simulations of future neutrino burst detection by underground detectors, or theoretical predictions for the relic supernova neutrino background. To follow the evolution from the onset of collapse to 20 s after the core bounce, we combine the results of neutrino-radiation hydrodynamic simulations for the early phase and quasi-static evolutionary calculations of neutrino diffusion for the late phase, with different values of shock revival time as a parameter that should depend on the still unknown explosion mechanism. We here describe the calculation methods and basic results including the dependence on progenitor models and the shock revival time. The neutrino data are publicly available electronically.

preprint2012arXiv

Infrared Spectral Energy Distribution of Galaxies in the AKARI All Sky Survey: Correlations with Galaxy Properties, and Their Physical Origin

We have studied the properties of more than 1600 low-redshift galaxies by utilizing high-quality infrared flux measurements of the AKARI All-Sky Survey and physical quantities based on optical and 21-cm observations. Our goal is to understand the physics determining the infrared spectral energy distribution (SED). The ratio of the total infrared luminosity L_TIR, to the star-formation rate (SFR) is tightly correlated by a power-law to specific SFR (SSFR), and L_TIR is a good SFR indicator only for galaxies with the largest SSFR. We discovered a tight linear correlation for normal galaxies between the radiation field strength of dust heating, estimated by infrared SED fits (U_h), and that of galactic-scale infrared emission (U_TIR ~ L_TIR/R^2), where R is the optical size of a galaxy. The dispersion of U_h along this relation is 0.3 dex, corresponding to 13% dispersion in the dust temperature. This scaling and the U_h/U_TIR ratio can be explained physically by a thin layer of heating sources embedded in a thicker, optically-thick dust screen. The data also indicate that the heated fraction of the total dust mass is anti-correlated to the dust column density, supporting this interpretation. In the large U_TIR limit, the data of circumnuclear starbursts indicate the existence of an upper limit on U_h, corresponding to the maximum SFR per gas mass of ~ 10 Gyr^{-1}. We find that the number of galaxies sharply drops when they become optically thin against dust-heating radiation, suggesting that a feedback process to galaxy formation (likely by the photoelectric heating) is working when dust-heating radiation is not self-shielded on a galactic scale. Implications are discussed for the M_HI-size relation, the Kennicutt-Schmidt relation, and galaxy formation in the cosmological context.

preprint2012arXiv

Suzaku Observation of Abell 1555 and Abell 1558: Searching for Non-thermal Emission from Large Scale Structure Formation

We report X-ray observations of two galaxy clusters Abell 1555 and Abell 1558 with Suzaku, which are included in a large scale filamentary structure and a supercluster, to search for non-thermal emission driven by shocks produced in structure formation. These two clusters are detected by Suzaku/XIS for the first time in the X-ray band of 0.5-7 keV. No significant flux is detected by HXD in the energy band of 13-40 keV, and upper limits are reported. From the analysis of the XIS data, we find that the spectrum of A1555 is fit by a thermal plus power-law model, significantly better than a single-temperature pure thermal spectrum. If this power-law component is due to inverse-Compton scattering, the fraction of total baryon energy imparted to non-thermal electrons is consistent with the typical value inferred from the observation of other clusters. However, other scenarios (e.g., under lying AGNs, multi-temperature thermal models) cannot be excluded and further investigation of this system is desired. Basic physical properties of A1555 (e.g., total mass) are also reported.

preprint2011arXiv

Supernovae in the Subaru Deep Field: the rate and delay-time distribution of type Ia supernovae out to redshift 2

The type Ia supernova (SN Ia) rate, when compared to the cosmic star formation history (SFH), can be used to derive the delay-time distribution (DTD) of SNe Ia, which can distinguish among progenitor models. We present the results of a SN survey in the Subaru Deep Field (SDF). Over a period of 3 years, we have observed the SDF on 4 independent epochs with Suprime-Cam on the Subaru 8.2-m telescope, with 2 nights of exposure per epoch, in the R, i', and z' bands. We have discovered 150 SNe out to redshift z~2. Using 11 photometric bands from the observer-frame far-ultraviolet to the near-infrared, we derive photometric redshifts for the SN host galaxies (for 24 we also have spectroscopic redshifts). This information is combined with the SN photometry to determine the type and redshift distribution of the SN sample. Our final sample includes 28 SNe Ia in the range 1.0<z<1.5 and 10 in the range 1.5<z<2.0. Our SN Ia rate measurements are consistent with those derived from the HST/GOODS sample, but the overall uncertainty of our 1.5<z<2.0 measurement is a factor of 2 smaller, of 35-50%. We find that the SN Ia rate evolution levels off at 1.0<z<2.0, but shows no sign of declining. Combining our SN Ia rate measurements and those from the literature, and comparing to a wide range of possible SFHs, the best-fit DTD is a power law of the form Psi(t) ~ t^beta, with beta = -1.1 +/- 0.1 (statistical) +/- 0.17 (systematic). This result is consistent with other recent DTD measurements at various redshifts and environments, and is in agreement with a generic prediction of the double-degenerate progenitor scenario for SNe Ia. By combining the contribution to iron production from core-collapse SNe, based on the wide range of SFHs, with that from SNe Ia, calculated with the best-fit DTD, we predict that the mean present-day cosmic iron abundance is in the range Z(Fe) = (0.09-0.37) Z(Fe,Sun) (abridged).

preprint2011arXiv

Suzaku X-ray Follow-up Observations of Seven Unassociated Fermi-LAT Gamma-ray Sources at High Galactic Latitudes

We report on our second-year campaign of X-ray follow-up observations of unidentified Fermi-LAT γ-ray sources at high Galactic latitudes (|b|>10 degree) using the X-ray Imaging Spectrometer onboard the Suzaku X-ray Observatory. In this second year of the project, seven new targets were selected from the First Fermi-LAT Catalog, and studied with 20-40 ks effective Suzaku exposures. We detected an X-ray point source coincident with the position of the recently discovered millisecond pulsar PSR J2302+4442 within the 95% confidence error circle of 1FGL J2302.8+4443. The X-ray spectrum of the detected counterpart was well fit by a blackbody model with temperature of kT ~0.3 keV, consistent with an origin of the observed X-ray photons from the surface of a rotating magnetized neutron star. For four other targets which were also recently identified with a normal pulsar (1FGL J0106.7+4853) and millisecond pulsars (1FGL J1312.6+0048, J1902.0-5110, and J2043.2+1709), only upper limits in the 0.5-10 keV band were obtained at the flux levels of ~10^{-14} erg cm^{-2} s^{-1}. A weak X-ray source was found in the field of 1FGL J1739.4+8717, but its association with the variable γ-ray emitter could not be confirmed with the available Suzaku data alone. For the remaining Fermi-LAT object 1FGL J1743.8-7620 no X-ray source was detected within the LAT 95% error ellipse. We briefly discuss the general properties of the observed high Galactic-latitude Fermi-LAT objects by comparing their multiwavelength properties with those of known blazars and millisecond pulsars.

preprint2011arXiv

Unraveling the Nature of Unidentified High Galactic Latitude Fermi/LAT Gamma-ray Sources with Suzaku

We report on the results of deep X-ray follow-up observations of four unidentified Fermi/LAT gamma-ray sources at high Galactic latitudes using Suzaku. The studied objects were detected with high significance during the first 3 months of Fermi/LAT operation, and subsequently better localized in the Fermi/LAT 1 year catalog (1FGL). Possible associations with pulsars and active galaxies have subsequently been discussed, and our observations provide an important contribution to this debate. In particular, an X-ray point source was found within the 95% confidence error circle of 1FGL J1231.1-1410. X-ray spectrum is well-fitted by a blackbody with an additional power-law. This supports the recently claimed identification of this source with a millisecond pulsar (MSP) PSR J1231-1411. Concerning 1FGL J1311.7-3429, two X-ray sources were found within the LAT error circle. Even though the X-ray spectral and variability properties were accessed, their nature and relationship with the gamma-ray source remain uncertain. We found several weak X-ray sources in the field of 1FGL J1333.2+5056, one coinciding with CLASS J1333+5057. We argue the available data are consistent with the association between these two objects. Finally, we have detected an X-ray source in the vicinity of 1FGL J2017.3+0603. This object was recently suggested to be associated with a newly discovered MSP PSR J2017+0603, because of the spatial-coincidence and the gamma-ray pulse detection. We have only detected the X-ray counterpart of the CLASS J2017+0603, while we determined an X-ray flux upper limit at the pulsar position. All in all, our studies indicate while a significant fraction of unidentified high Galactic latitude gamma-ray sources is related to the pulsar and blazar phenomena, associations with other classes of astrophysical objects are still valid options.

preprint2010arXiv

Contribution from Star-Forming Galaxies to the Cosmic Gamma-Ray Background Radiation

We present a new theoretical calculation of the contribution to the extragalactic gamma-ray background radiation (EGRB) from star-forming galaxies, based on a state-of-the-art model of hierarchical galaxy formation that is in quantitative agreement with a variety of observations of local and high-redshift galaxies. Gamma-ray luminosity ($L_γ$) and spectrum of galaxies are related to star formation rate ($ψ$), gas mass ($M_gas$), and star formation mode (quiescent or starburst) of model galaxies using latest observed data of nearby galaxies. We try the two limiting cases about gamma-ray production: the escape limit ($L_γ\propto ψM_gas$) and the calorimetric limit ($L_γ\propto ψ$), and our standard model predicts 7 and 4% contribution from star-forming galaxies to the total EGRB flux (including bright resolved sources) recently reported by the Fermi Gamma-Ray Space Telescope. Systematic uncertainties do not allow us to determine the EGRB flux better than by a factor of ~2. The predicted number of nearby galaxies detectable by Fermi is consistent with the observation. Intergalactic absorption by pair-production attenuates the EGRB flux only by a modest factor of ~1.3 at the highest Fermi energy band, and the reprocessed cascade emission does not significantly alter EGRB at lower photon energies. The sum of the known contributions from AGNs and star-forming galaxies can explain a large part of EGRB, with a remarkable agreement between the predicted model spectrum and observation.

preprint2010arXiv

Luminosity Distribution of Gamma-Ray Burst Host Galaxies at redshift z=1 in Cosmological Smoothed Particle Hydrodynamic Simulations: Implications for the Metallicity Dependence of GRBs

We study the relationship between the metallicity of gamma-ray burst (GRB) progenitors and the probability distribution function (PDF) of GRB host galaxies as a function of luminosity using cosmological hydrodynamic simulations of galaxy formation. We impose a maximum limit to the gas metallicity in which GRBs can occur, and examine how the predicted luminosity PDF of GRB host galaxies changes in the simulation. We perform the Kolmogorov-Smirnov test, and show that the result from our simulation agrees with the observed luminosity PDF of core-collapse supernovae (SNe) host galaxies when we assume that the core-collapse SNe trace star formation. When we assume that GRBs occur only in a low-metallicity environment with $Z\lesssim 0.1 \Zsun$, GRBs occur in lower luminosity galaxies, and the simulated luminosity PDF becomes quantitatively consistent with the observed luminosity PDF. The observational bias against the host galaxies of optically dark GRBs owing to dust extinction may be another reason for the lower luminosities of GRB host galaxies, but the observed luminosity PDF of GRB host galaxies cannot be reproduced solely by the dust bias in our simulation.

preprint2010arXiv

Lyman Alpha Emitters at z=7 in the Subaru/XMM-Newton Deep Survey Field: Photometric Candidates and Luminosity Function

We conducted a deep narrowband NB973 (FWHM = 200 A centered at 9755 A) survey of z=7 Lyman alpha emitters (LAEs) in the Subaru/XMM-Newton Deep Survey Field, using the fully depleted CCDs newly installed on the Subaru Telescope Suprime-Cam, which is twice more sensitive to z=7 Lyman alpha at ~ 1 micron than the previous CCDs. Reaching the depth 0.5 magnitude deeper than our previous survey in the Subaru Deep Field that led to the discovery of a z=6.96 LAE, we detected three probable z=7 LAE candidates. Even if all the candidates are real, the Lyman alpha luminosity function (LF) at z=7 shows a significant deficit from the LF at z=5.7 determined by previous surveys. The LAE number and Lyman alpha luminosity densities at z=7 is ~ 7.7-54% and ~5.5-39% of those at z=5.7 to the Lyman alpha line luminosity limit of L(Ly-alpha) >~ 9.2 x 10^{42} erg s^{-1}. This could be due to evolution of the LAE population at these epochs as a recent galaxy evolution model predicts that the LAE modestly evolves from z=5.7 to 7. However, even after correcting for this effect of galaxy evolution on the decrease in LAE number density, the z=7 Lyman alpha LF still shows a deficit from z=5.7 LF. This might reflect the attenuation of Lyman alpha emission by neutral hydrogen remaining at the epoch of reionization and suggests that reionization of the universe might not be complete yet at z=7. If we attribute the density deficit to reionization, the intergalactic medium (IGM) transmission for Lyman alpha photons at z=7 would be 0.4 <= T_{Ly-alpha}^{IGM} <= 1, supporting the possible higher neutral fraction at the earlier epochs at z > 6 suggested by the previous surveys of z=5.7-7 LAEs, z ~ 6 quasars and z > 6 gamma-ray bursts.

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

Lyman Alpha Emitters in Hierarchical Galaxy Formation II. UV Continuum Luminosity Function and Equivalent Width Distribution

We present theoretical predictions of UV continuum luminosity function (UV LF) and Lya equivalent width (EW) distribution of Lyman alpha emitters (LAEs) in the framework of the hierarchical clustering model of galaxy formation. The model parameters about LAEs were determined by fitting to the observed Lya LF at z=5.7 in our previous study, and the fit indicates that extinction of Lya photons by dust is significantly less effective than that of UV continuum photons, implying clumpy dust distribution in interstellar medium. We then compare the predictions about UV LFs and EW distributions with a variety of observations at z ~ 3-6, allowing no more free parameters and paying careful attention to the selection conditions of LAEs in each survey. We find that the predicted UV LFs and EW distributions are in nice agreement with observed data, and especially, our model naturally reproduces the existence of large EW LAEs (> 240 A) without introducing Pop III stars or top-heavy initial mass function. We show that both the stellar population (young age and low metallicity) and extinction by clumpy dust are the keys to reproduce large EW LAEs. The evidence of EW enhancement by clumpy dust is further strengthened by the quantitative agreement between our model and recent observations about a positive correlation between EW and extinction. The observed trend that brighter LAEs in UV continuum tend to have smaller mean EW is also reproduced, and the clumpy dust is playing an important role again for this trend. We suggested in our previous study that the transmission of intergalactic medium for Lya emission rapidly decreases from z ~ 6 to 7 by the fitting to Lya LFs, and this evidence is quantitatively strengthened by the comparison with the UV LF and EW distribution at z ~ 6.6.

preprint2009arXiv

Size Evolution of Early-Type Galaxies and Massive Compact Objects as the Dark Matter

The dramatic size evolution of early-type galaxies from z ~ 2 to 0 poses a new challenge in the theory of galaxy formation, which may not be explained by the standard picture. It is shown here that the size evolution can be explained if the non-baryonic cold dark matter is composed of compact objects having a mass scale of ~10^5 M_sun. This form of dark matter is consistent with or only weakly constrained by the currently available observations. The kinetic energy of the dark compact objects is transferred to stars by dynamical friction, and stars around the effective radius are pushed out to larger radii, resulting in a pure size evolution. This scenario has several good properties to explain the observations, including the ubiquitous nature of size evolution and faster disappearance of higher density galaxies.

preprint2009arXiv

Subaru FOCAS Spectroscopic Observations of High-Redshift Supernovae

We present spectra of high-redshift supernovae (SNe) that were taken with the Subaru low resolution optical spectrograph, FOCAS. These SNe were found in SN surveys with Suprime-Cam on Subaru, the CFH12k camera on the Canada-France-Hawaii Telescope (CFHT), and the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST). These SN surveys specifically targeted z>1 Type Ia supernovae (SNe Ia). From the spectra of 39 candidates, we obtain redshifts for 32 candidates and spectroscopically identify 7 active candidates as probable SNe Ia, including one at z=1.35, which is the most distant SN Ia to be spectroscopically confirmed with a ground-based telescope. An additional 4 candidates are identified as likely SNe Ia from the spectrophotometric properties of their host galaxies. Seven candidates are not SNe Ia, either being SNe of another type or active galactic nuclei. When SNe Ia are observed within a week of maximum light, we find that we can spectroscopically identify most of them up to z=1.1. Beyond this redshift, very few candidates were spectroscopically identified as SNe Ia. The current generation of super red-sensitive, fringe-free CCDs will push this redshift limit higher.

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.

preprint2008arXiv

Delay Time Distribution Measurement of Type Ia Supernovae by the Subaru/XMM-Newton Deep Survey and Implications for the Progenitor

The delay time distribution (DTD) of type Ia supernovae (SNe Ia) from star formation is an important clue to reveal the still unknown progenitor system of SNe Ia. Here we report on a measurement of the SN Ia DTD in a delay time range of t_Ia = 0.1-8.0 Gyr by using the faint variable objects detected in the Subaru/XMM-Newton Deep Survey (SXDS) down to i' ~ 25.5. We select 65 SN candidates showing significant spatial offset from nuclei of the host galaxies having old stellar population at z ~ 0.4-1.2, out of more than 1,000 SXDS variable objects. Although spectroscopic type classification is not available for these, we quantitatively demonstrate that more than ~80% of these should be SNe Ia. The DTD is derived using the stellar age estimates of the old galaxies based on 9 band photometries from optical to mid-infrared wavelength. Combined with the observed SN Ia rate in elliptical galaxies at the local universe, the DTD in t_Ia ~ 0.1-10 Gyr is well described by a featureless power-law as f_D(t_Ia) \propto t_Ia^{-1}. The derived DTD is in excellent agreement with the generic prediction of the double-degenerate scenario, giving a strong support to this scenario. In the single-degenerate (SD) scenario, although predictions by simple analytic formulations have broad DTD shapes that are similar to the observation, DTD shapes calculated by more detailed binary population synthesis tend to have strong peaks at characteristic time scales, which do not fit the observation. This result thus indicates either that the SD channel is not the major contributor to SNe Ia in old stellar population, or that improvement of binary population synthesis theory is required. Various sources of systematic uncertainties are examined and tested, but our main conclusions are not affected significantly.

preprint2008arXiv

The Reionization and Galaxy Evolution Probed by z=7 Lyman Alpha Emitters

We made a narrowband NB973 (bandwidth of 200A at 9755A) imaging of the Subaru Deep Field (SDF) and found two z=7 Lyman alpha emitter (LAE) candidates down to NB973=24.9. Carrying out deep follow-up spectroscopy, we identified one of them as a real z=6.96 LAE. This has shown that galaxy formation was in progress just 750 Myr after the Big Bang. Meanwhile, the Lyman alpha line luminosity function of LAE is known to decline from z=5.7 to 6.6 in the SDF. L* at z=6.6 is 40-60% of that at z=5.7. We also confirm that the number density of z=7 LAE is only 17% of the density at z=6.6 comparing the latest SDF LAE samples. This series of significant decreases in LAE density with increasing redshift can be the result of galaxy evolution during these epochs. However, using the UV continuum luminosity functions of LAEs, those of Lyman break galaxies, and a LAE evolution model based on the hierarchical clustering, we find that galaxy evolution alone cannot explain all the decrease in density. This extra density deficit can be interpreted as the attenuation of the Lyman alpha photons from LAEs due to a rapid evolution of neutral hydrogen fraction during the ongoing cosmic reionization at z~6.6-7.

preprint2006arXiv

A RIAF Interpretation for the Past Higher Activity of the Galactic Center Black Hole and the 511 keV Annihilation Emission

There are several lines of evidence that the super-massive black hole at the Galactic center had higher activities in the past than directly observed at present. Here I show that these lines of evidence can quantitatively and consistently be explained if the mean accretion rate during the past ~10^7 yrs has been ~10^{3-4} times higher than the current rate, by the picture of radiatively inefficient accretion flow (RIAF) and associated outflow that has been successfully applied to Sgr A*. I argue that this increased rate and its duration are theoretically reasonable in the Galactic center environment, while the accretion rate suddenly dropped about 300 years ago most likely because of the shell passage of the supernova remnant Sgr A East. The chance probability of witnessing Sgr A* in such a low state is not extremely small (~0.5%). The outflow energetics is sufficient to keep the hot (~8 keV) diffuse gas observed in the Galactic center region. Then, I show that a significant amount of positrons should have been created around the event horizon during the higher activity phase, and injected into interstellar medium by the outflow. The predicted positron production rate and propagation distance are close to those required to explain the observed 511 keV annihilation line emission from the Galactic bulge, giving a natural explanation for the large bulge-to-disk ratio of the emission. The expected injection energy into interstellar medium is about 1 MeV, which is also favorable as an explanation of the 511 keV line emission.

preprint2006arXiv

Gamma-Ray Luminosity Function of Blazars and the Cosmic Gamma-Ray Background: Evidence for the Luminosity Dependent Density Evolution

We present a comprehensive study for the gamma-ray luminosity function (GLF) of blazars and their contribution to the extragalactic diffuse gamma-ray background (EGRB). Radio and gamma-ray luminosity correlation is introduced to take into account the radio detectability which is important for the blazar identification. Previous studies considered only pure luminosity evolution (PLE) or pure density evolution, but we introduce the luminosity dependent density evolution (LDDE) model, which is favored from the evolution of X-ray luminosity function (XLF) of AGNs. The model parameters are constrained by likelihood analyses about the observed redshift and gamma-ray flux distributions of the EGRET blazars. We find that the LDDE model gives a better fit to the observed distributions than the PLE model, indicating that the LDDE model is also appropriate for gamma-ray blazars, and that the jet activity is universally correlated with the accretion history of AGNs. The normalization between the GLF and XLF is consistent with the unified picture of AGNs, when the beaming and a reasonable duty cycle of jet activity are taken into account. We then find that only 25--50% of the EGRB can be explained by unresolved blazars with the best-fit LDDE parameters. Unresolved blazars can account for all the EGRB only with a steeper index of the faint-end slope of the GLF, which is marginally consistent with the EGRET data but inconsistent with that of the XLF. Therefore unresolved AGNs cannot be the dominant source of the EGRB, unless there is a new population of gamma-ray emitting AGNs that evolves differently from the XLF of AGNs. Predictions for the GLAST mission are made, and we find that the best-fit LDDE model predicts about 3000 blazars in the entire sky, which is considerably fewer than a previous estimate.

preprint2006arXiv

Implications for the Cosmic Reionization from the Optical Afterglow Spectrum of the Gamma-Ray Burst 050904 at z = 6.3

The gamma-ray burst (GRB) 050904 at z = 6.3 provides the first opportunity of probing the intergalactic medium (IGM) by GRBs at the epoch of the reionization. Here we present a spectral modeling analysis of the optical afterglow spectrum taken by the Subaru Telescope, aiming to constrain the reionization history. The spectrum shows a clear damping wing at wavelengths redward of the Lyman break, and the wing shape can be fit either by a damped Ly alpha system with a column density of log (N_HI/cm^{-2}) ~ 21.6 at a redshift close to the detected metal absorption lines (z_metal = 6.295), or by almost neutral IGM extending to a slightly higher redshift of z_{IGM,u} ~ 6.36. In the latter case, the difference from z_metal may be explained by acceleration of metal absorbing shells by the activities of the GRB or its progenitor. However, we exclude this possibility by using the light transmission feature around the Ly beta resonance, leading to a firm upper limit of z_{IGM,u} < 6.314. We then show an evidence that the IGM was largely ionized already at z=6.3, with the best-fit neutral fraction of IGM, x_HI = N_HI/N_H = 0.00, and upper limits of x_HI < 0.17 and 0.60 at 68 and 95% C.L., respectively. This is the first direct and quantitative upper limit on x_HI at z > 6. Various systematic uncertainties are examined, but none of them appears large enough to change this conclusion. To get further information on the reionization, it is important to increase the sample size of z >~ 6 GRBs, in order to find GRBs with low column densities (log N_HI <~ 20) within their host galaxies, and for statistical studies of Ly alpha line emission from host galaxies.

preprint1997arXiv

Cosmological Gamma-Ray Bursts and Evolution of Galaxies

Evolution of the rate density of cosmological gamma-ray bursts (GRBs) is calculated and compared to the BATSE brightness distribution in the context of binary neutron-star mergers as the source of GRBs, taking account of the realistic star formation history in the universe and evolution of compact binary systems. We tried two models of the evolution of cosmic star formation rate (SFR): one is based on recent observations of SFRs at high redshifts, while the other is based on a galaxy evolution model of stellar population synthesis that reproduces the present-day colors of galaxies. It is shown that the binary merger scenario of GRBs naturally results in the comoving rate-density evolution of \propto (1+z)^{2-2.5} up to z ~ 1, that has been suggested independently from the compatibility between the number-brightness distribution and duration-brightness correlation. If the cosmic SFR has its peak at z ~ 1--2 as suggested by recent observations, the effective power-index of GRB photon spectrum, α>~ 1.5$ is favored, that is softer than the recent observational determination of α= 1.1 \pm 0.3. However, high redshift starbursts (z >~ 5) in elliptical galaxies, that have not yet been detected, can alleviate this discrepancy. The redshift of GRB970508 is likely about 2, just below the upper limit that is recently determined, and the absorption system at z = 0.835 seems not to be the site of the GRB.

preprint1996arXiv

Resonant Spin-Flavor Conversion of Supernova Neutrinos and Deformation of the Electron Antineutrino Spectrum

The neutrino spin-flavor conversion of \barν_e and ν_μwhich is induced by the interaction of the Majorana neutrino magnetic moment and magnetic fields in the collapse-driven supernova is investigated in detail. We calculate the conversion probability by using the latest precollapse models of Woosley and Weaver (1995), and also those of Nomono and Hashimoto (1988), changing the stellar mass and metallicity in order to estimate the effect of the astrophysical uncertainties. Contour maps of the conversion probability are given for all the models as a function of neutrino mass squared difference and the neutrino magnetic moment times magnetic fields. It is shown that in the solar metallicity models some observational effects are expected with Δm^2 = 10^{-5}--10^{-1} [eV^2] and μ_ν>~ 10^{-12} (10^9 G / B_0) [μ_B], where B_0 is the strength of the magnetic fields at the surface of the iron core. We also find that although the dependence on the stellar models or stellar mass is not so large, the metallicity of precollapse stars has considerable effects on this conversion. Such effects may be seen in a supernova in the Large or Small Magellanic Clouds, and should be taken into account when one considers an upper bound on μ_νfrom the SN1987A data.

preprint1995arXiv

Spectrum of the relic neutrino background from past supernovae and cosmological models

It is greatly expected that the relic neutrino background from past supernovae is detected by Superkamiokande (SK) which is now under construction. We calculate the spectrum and the event rate at SK systematically by using the results of simulations of a supernova explosion and reasonable supernova rates. We also investigate the effect of a cosmological constant, $Λ$, on the spectrum, since some recent cosmological observations strongly suggest the existence of $Λ$. We find following results. 1) The spectrum has a peak at about 3 MeV, which is much lower than that of previous estimate ($6 \sim$ 10 MeV). 2) The event rate at SK in the range from 10 MeV to 50 MeV, where the relic neutrinos from past supernovae is dominant, is about $25 \, {h_{50}}^2 \, \left( \frac{R_{SN}}{0.1 {\rm yr^{-1}}} \right) \left(\frac{n_G \, h_{50}^{-3}}{0.02{\rm Mpc}^{-3}}\right)$ events per year, where $R_{SN}$ is the supernova rate in a galaxy, $n_G$ is the number density of galaxies, and $h_{50} = H_0 /$ (50km/s/Mpc), where $H_0$ is the Hubble constant. 3) The event rate is almost insensitive to $Λ$. The flux increases in the low energy side ($< 10$ MeV) with increasing $Λ$, but decreases in the high energy side (10 MeV $<$) in models in which the integrated number of supernovae in one galaxy is fixed.

preprint1995arXiv

Spectrum of the Supernova Relic Neutrino Background and Evolution of Galaxies

The spectrum of the supernova relic neutrino background (SRN) from collapse-driven supernovae ever occurred in the universe is calculated by using a realistic, time-dependent supernova rate derived from a standard model of galaxy evolution based on the population synthesis method. The SRN spectrum we show here is the most realistic at present, because the largest uncertainty in previous theoretical predictions has come from unrealistic assumptions of the supernova rate so far made. The SRN is one of the targets of the Superkamiokande (SK) detector which will be constructed in a year and the SRN, if at all detected, would provide a new tool to probe the history of supernova explosions in the universe. The expected event rate at the SK is therefore calculated in this paper. Our major results include: (1) the supernova rate is much higher in the early phase of evolution of galaxies and there appears a hump in the SRN spectrum in the low-energy region of $\ltilde 5$ MeV, (2) the SRN flux depends on the Hubble constant ($H_0$) in a way approximately proportional to $H_0^2$ and only weakly on the density parameter of the universe ($Ω_0$) and a cosmological constant ($λ_0$), and (3) the plausible event rate at the SK is 1.2 yr$^{-1}$ in the observable energy range. Such a low event rate is due mainly to a quite low supernova rate at present which is averaged over the morphological types of galaxies. The most optimistic rate in our model is found to be 4.7 yr$^{-1}$, and if more events are detected, we will have to reconsider our current understanding of collapse-driven supernovae and evolution of galaxies.