Source author record

Keitaro Takahashi

Keitaro Takahashi 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

66works
11topics
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

66 published item(s)

preprint2026arXiv

Realistic assessment of a single gravitational wave source localization taking into account precise pulsar distances with pulsar timing arrays

Pulsar timing arrays (PTAs) are anticipated to detect continuous gravitational waves (GWs) from individual supermassive black hole binaries (SMBHBs) in the near future. To identify the host galaxy of a GW source, PTAs require significantly improved angular resolution beyond the typical range of 100-1000 square degrees achieved by recent continuous GW searches. In this study, we investigate how precise pulsar distance measurements can enhance the localization of a single GW source. Accurate distance information, comparable to or better than the GW wavelength (typically 1~pc) can refine GW source localization. In the near future, with the advent of Square Kilometre Array (SKA), such high-precision distance measurements will be feasible for a few nearby pulsars. We focus on the relatively nearby pulsars J0437-4715 (156 pc) and J0030+0451 (331 pc), incorporating their actual distance uncertainties based on current VLBI measurements and the anticipated precision of the SKA-era. By simulating 87 pulsars with the GW signal and Gaussian white noise in the timing residuals, we assess the impact of the pulsar distance information on GW source localization. Our results show that without precise distance information, localization remains insufficient to identify host galaxies under 10 ns noise. However, incorporating SKA-era distance precision for nearby pulsars J0437-4715 and J0030+0451 can reduce localization uncertainties to the required level of $10^{-3}$ $\rm deg^{2}$. Localization accuracy strongly depends on the geometric configuration of pulsars with well-measured distances and improves notably near and between such pulsars. The improvement of the localization will greatly aid in identifying the host galaxy of a GW source and constructing an SMBHB catalog. It will further enable follow-up electromagnetic observations to investigate the SMBHB in greater detail.

preprint2025arXiv

Effects of Coronal Mass Ejection on PSR J1022+1001 and Possible Mode Change of PSR J2145-0750 in the InPTA DR2

The Indian Pulsar Timing Array (InPTA) has recently published its second data release (DR2), comprising the timing analysis of seven years of data on 27 millisecond pulsars (MSPs), observed simultaneously in the 300-500 MHz (band 3) and 1260-1460 MHz (band 5), using the upgraded Giant Metrewave Radio Telescope (uGMRT). The low-frequency data, particularly in band 3, is highly sensitive to propagation effects such as dispersion measure (DM) fluctuations, which can be imprints of some astrophysical phenomena (scientific outliers). Here, we analyze the two outliers of possible astrophysical origin coming from the band 3 DM time series of two pulsars: PSR J1022+1001, with an ecliptic latitude of -0.06 degree, and PSR J2145-0750, one of the brightest MSPs, with multi-component profile morphology. Our study reveals compelling evidence for a coronal mass ejection (CME) event traced in the data of PSR J1022+1001, and reports evidence for a potential mode-changing event in PSR J2145-0750. By contrasting these two cases, we show that DM fluctuations due to CME interacions and intrinsic mode-changing events produce distinct observational signatures, enabling a physically informed classification of scientific outliers in PTA datasets. Extending the analyses presented here to the full sample of InPTA-DR2 pulsars is expected to reveal additional CME events, and possible mode-changing events. Such detections will not only improve our understanding of solar and pulsar magnetospheric plasma interactions but will also enable more accurate modelling of DM variations, leading to improved pulsar timing solutions, which are crucial for high-precision Pulsar Timing Array (PTA) science.

preprint2022arXiv

A Parkes "Murriyang" Search for Pulsars and Transients in the Large Magellanic Cloud

We report on the latest results of a Parkes multibeam survey for pulsars and dispersed radio bursts in the Large Magellanic Cloud (LMC). We conducted both periodicity and single-pulse searches at a much larger range of trial dispersion measures (DMs) than previously searched. We detected 229 single pulses with signal-to-noise ratio ($\rm{S/N}) > 7$ that were classified by the deep learning network FETCH as being real (with $> 90$% likelihood), of which 9 were from the known giant-pulse-emitting pulsar PSR B0540$-$69. Two possibly repeating sources were detected with DMs suggesting that they lie within the LMC, but these require confirmation. Only 3 of the 220 unknown pulses had S/N greater than 8, and the DM distribution for these detected pulses follows an exponential falloff with increasing DM and does not show any excess of signals at DM values expected for the LMC. These features suggest that the detected pulses are not likely to be real, although they are visually compelling. We also report the discovery of a new pulsar (PSR J0556$-$67) in our periodicity search. This pulsar has a spin period of 791 ms, a DM of 71 cm$^{-3}$ pc, an estimated 1400 MHz flux density of $\sim 0.12$ mJy, and shows no evidence of binary motion. Future observations may be able to confirm whether any of the weak but promising candidates in our single pulse and periodicity searches, including our two possible repeaters, are real or not.

preprint2022arXiv

Low-frequency wideband timing of InPTA pulsars observed with the uGMRT

High-precision measurements of the pulsar dispersion measure (DM) are possible using telescopes with low-frequency wideband receivers. We present an initial study of the application of the wideband timing technique, which can simultaneously measure the pulsar times of arrival (ToAs) and DMs, for a set of five pulsars observed with the upgraded Giant Metrewave Radio Telescope (uGMRT) as part of the Indian Pulsar Timing Array (InPTA) campaign. We have used the observations with the 300-500 MHz band of the uGMRT for this purpose. We obtain high precision in DM measurements with precisions of the order 10^{-6}cm^{-3}pc. The ToAs obtained have sub-μs precision and the root-mean-square of the post-fit ToA residuals are in the sub-μs range. We find that the uncertainties in the DMs and ToAs obtained with this wideband technique, applied to low-frequency data, are consistent with the results obtained with traditional pulsar timing techniques and comparable to high-frequency results from other PTAs. This work opens up an interesting possibility of using low-frequency wideband observations for precision pulsar timing and gravitational wave detection with similar precision as high-frequency observations used conventionally.

preprint2022arXiv

On the Potential of Faraday Tomography to Identify Shock Structures in Supernova Remnants

Knowledge about the magnetic fields in supernova remnants (SNRs) is of paramount importance for constraining Galactic cosmic ray acceleration models. It could also indirectly provide information on the interstellar magnetic fields. In this paper, we predict the Faraday dispersion functions (FDFs) of SNRs for the first time. For this study, we use the results of three dimensional (3D) ideal magnetohydrodynamic (MHD) simulations of SNRs expanding into a weak, regular magnetic field. We present the intrinsic FDFs of the shocked region of SNRs for different viewing angles. We find that the FDFs are generally Faraday complex, which implies that conventional rotation measure study is not sufficient to obtain the information on the magnetic fields in the shocked region and Faraday tomography is necessary. We also show that the FDF allows to derive the physical-depth distribution of polarization intensity when the line of sight is parallel to the initial magnetic field orientation. Furthermore, we demonstrate that the location of contact discontinuity can be identified from the radial profile of the width of the FDF with the accuracy of 0.1-0.2 pc.

preprint2022arXiv

Testing the Non-circularity of the Spacetime around Sagittarius A* with Orbiting Pulsars

A disformal Kerr black hole solution is a rotating black hole solution in a modified gravity theory which breaks the circular condition of spacetime differently from the case of the Kerr spacetime. In this paper, assuming that Sagittarius A* (Sgr A*) is a disformal Kerr black hole, we examine the potential to test the spacetime geometry with a hypothetical pulsar whose orbital elements are similar to those of the S2/S0-2 star. By numerically solving the equations of motion for the pulsar and photons emitted from it, we calculate the apparent position of the pulsar and the time of arrival (TOA) of the emitted pulse signals. Our analysis shows that the magnitude of the difference in the TOAs reaches the order of $10\>{\rm ms}$ if the deviation from the Kerr spacetime is significant. The time difference is mainly caused by the non-circularity of the spacetime at the $1.5$ post-Newtonian order. The accuracy of the TOA measurement by a future radio telescope named the Square Kilometer Array (SKA) is between about $0.1\>{\rm ms}$ and $10\>{\rm ms}$ for a normal pulsar. Thus, we expect that the SKA can distinguish the disformal Kerr black hole from the Kerr black hole through the non-circularity of the spacetime around Sgr A*.

preprint2022arXiv

Wavelets and sparsity for Faraday tomography

Faraday tomography through broadband polarimetry can provide crucial information on magnetized astronomical objects, such as quasars, galaxies, or galaxy clusters. However, the limited wavelength coverage of the instruments requires that we solve an ill-posed inverse problem when we want to obtain the Faraday dispersion function (FDF), a tomographic distribution of the magnetoionic media along the line of sight. This paper explores the use of wavelet transforms and the sparsity of the transformed FDFs in the form of wavelet shrinkage (WS) for finding better solutions to the inverse problem. We recently proposed the Constraining and Restoring iterative Algorithm for Faraday Tomography (CRAFT; Cooray et al. 2021), a new flexible algorithm that showed significant improvements over the popular methods such as Rotation Measure Synthesis. In this work, we introduce CRAFT+WS, a new version of CRAFT incorporating the ideas of wavelets and sparsity. CRAFT+WS exhibit significant improvements over the original CRAFT when tested for a complex FDF of realistic Galactic model. Reconstructions of FDFs demonstrate super-resolution in Faraday depth, uncovering previously unseen Faraday complexities in observations. The proposed approach will be necessary for effective cosmic magnetism studies using the Square Kilometre Array and its precursors.

preprint2020arXiv

Artificial neural networks for selection of pulsar candidates from the radio continuum surveys

Pulsar search with time-domain observation is very computationally expensive and data volume will be enormous with the next generation telescopes such as the Square Kilometre Array. We apply artificial neural networks (ANNs), a machine learning method, for efficient selection of pulsar candidates from radio continuum surveys, which are much cheaper than time-domain observation. With observed quantities such as radio fluxes, sky position and compactness as inputs, our ANNs output the "score" that indicates the degree of likeliness of an object to be a pulsar. We demonstrate ANNs based on existing survey data by the TIFR GMRT Sky Survey (TGSS) and the NRAO VLA Sky Survey (NVSS) and test their performance. Precision, which is the ratio of the number of pulsars classified correctly as pulsars to that of any objects classified as pulsars, is about 96$\%$. Finally, we apply the trained ANNs to unidentified radio sources and our fiducial ANN with five inputs (the galactic longitude and latitude, the TGSS and NVSS fluxes and compactness) generates 2,436 pulsar candidates from 456,866 unidentified radio sources. These candidates need to be confirmed if they are truly pulsars by time-domain observations. More information such as polarization will narrow the candidates down further.

preprint2020arXiv

Detectability of 21cm signal during the Epoch of Reionization with 21cm-Lyman-$α$ emitter cross-correlation. III. Model dependence

Detecting $\rm H_I$ 21cm line in the intergalactic medium (IGM) during the Epoch of Reionization (EoR) suffers from foreground contamination such as Galactic synchrotron and extragalactic radio sources. Cross-correlation between the 21cm line and Lyman-$α$ emitter (LAE) galaxies is a powerful tool to identify the 21cm signal since the 21cm line emission has correlation with LAEs while the LAEs are statistically independent of the foregrounds. So far, the detectability of 21cm-LAE cross-power spectrum has been investigated with simple LAE models where the observed Ly$α$ luminosity is proportional to the dark matter halo mass. However, the previous models were inconsistent with the latest observational data of LAEs obtained with Subaru/Hyper Suprime-Cam (HSC). Here, we revisit the detectability of 21cm-LAE cross-power spectrum adopting a state-of-the-art LAE model consistent with all Subaru/HSC observations such as the Ly$α$ luminosity function, LAE angular auto-correlation, and the LAE fractions in the continuum selected galaxies. We find that resultant cross-power spectrum with the updated LAE model is reduced at small scales ($k\sim 1\ \rm Mpc^{-1}$) compared to the simple models, while the amplitudes at large scales ($k \lesssim 0.2 \ \rm Mpc^{-1}$) are not affected so much. We conclude that the large-scale signal would be detectable with Square Kilometre Array (SKA) and HSC LAE cross-correlation but detecting the small scale signal would require an extended HSC LAE survey with an area of $\sim 75\ \rm deg^2$ or 3000 hrs observation time of 21cm line with SKA.

preprint2020arXiv

Faraday dispersion function of disk galaxies with axisymmetric global magnetic fields I

Faraday tomography is a novel method to probe 3-dimensional structure of magnetic fields of polarized radio sources. In this paper, we investigate Faraday dispersion function (FDF) of disk galaxies extending a simple analytic model of galactic magnetic fields developed in Ideguchi et al. (2017). The model consists of axisymmetric coherent fields and turbulent fields and we consider the effects of inclination, relative amplitude of coherent and turbulent magnetic fields and pitch angle of coherent fields. Our simple model makes it easy to obtain physical interpretation of FDFs and helps understanding observational results. We find FDFs have two peaks when galaxies are observed with non zero inclination. The gap and relative height of two peaks are dependent on the inclination angle and pitch angle but are not affected by relative amplitude of coherent and turbulent magnetic fields so much. These findings give us an important caution that two peaks in observed FDFs do not necessarily imply the presence of two separate radio sources within a beam.

preprint2020arXiv

Probing Small Scale Primordial Power Spectrum with 21cm Line Global Signal

We argue that the global signal of neutral hydrogen 21cm line can be a powerful probe of primordial power spectrum on small scales. Since the amplitude of small scale primordial fluctuations is important to determine the early structure formation and the timing when the sources of Lyman $α$ photons are produced, they in turn affect the neutral hydrogen 21cm line signal. We show that the information of the position of the absorption trough can severely constrain the small scale amplitude of primordial fluctuations once astrophysical parameters relevant to the 21cm line signal are fixed. We also discuss how the uncertainties of astrophysical parameters affect the constraints.

preprint2016arXiv

21cm-line bispectrum as method to probe Cosmic Dawn and Epoch of Reionization

Redshifted 21cm signal is a promising tool to investigate the state of intergalactic medium (IGM) in the Cosmic Dawn (CD) and Epoch of Reionization(EoR). In our previous work (Shimabukuro et al 2015), we studied the variance and skewness to give a clear interpretation of 21cm power spectrum and found that skewness is a good indicator of the epoch when X-ray heating becomes effective. Thus, the non-Gaussian feature of the spatial distribution of the 21cm signal is expected to be useful to investigate the astrophysical effects in the CD and EoR. In this paper, in order to investigate such a non-Gaussian feature in more detail, we focus on the bispectrum of the 21cm signal. It is expected that the 21cm brightness temperature bispectrum is produced by non-gaussianity due to the various astrophysical effects such as Wouthysen-Field (WF) effect, X-ray heating and reionization. We study the various properties of 21cm bispectrum such as scale dependence, shape dependence and redshift evolution. And also we study the contribution from each component of 21cm bispectrum. We find that the contribution from each component has characteristic scale-dependent feature, and it is expected that we could obtain more detailed information on the IGM in the CD and EoR by using the 21cm bispectrum in the future experiments, combined with the power spectrum and skewness.

preprint2016arXiv

Cosmology with the Square Kilometre Array by SKA-Japan

In the past several decades, the standard cosmological model has been established and its parameters have been measured to a high precision, while there are still many of the fundamental questions in cosmology; such as the physics in the very early Universe, the origin of the cosmic acceleration and the nature of the dark matter. The future world's largest radio telescope, Square Kilometre Array (SKA), will be able to open the new frontier of cosmology and will be one of the most powerful tools for cosmology in the next decade. The cosmological surveys conducted by the SKA would have the potential not only to answer these fundamental questions but also deliver the precision cosmology. In this article we briefly review the role of the SKA from the view point of the modern cosmology. The cosmology science led by the SKA-Japan Consortium (SKA-JP) Cosmology Science Working Group is also discussed.

preprint2016arXiv

Performance test of RM CLEAN and its evaluation with chi-square value

RM CLEAN is a standard method to reconstruct the distribution of cosmic magnetic fields and polarized sources along the line of sight (LOS) from observed polarization spectrum. This method is similar to the CLEAN algorithm for aperture synthesis radio telescope images but it is rather unclear in what cases RM CLEAN works well. In this paper, we evaluate the performance of RM CLEAN by simulating spectro-polarimetric observations of two compact sources located in the same LOS, varying the relative initial polarization angle and Faraday depth systematically. Especially, we focus on if the two polarized sources can be resolved in the Faraday depth space and how well the source parameters can be estimated. We confirm the previous studies that two sources cannot be resolved when they are closely located in the Faraday depth space for specific values of the relative initial polarization angle. Further, we calculate the chi-square value for the fit between the mock data of polarization spectrum and the one from RM CLEAN. Then we find that the chi-squared value is not always significantly large even when RM CLEAN gives wrong results.

preprint2016arXiv

Probing higher-order primordial non-Gaussianity with galaxy surveys

With a radio continuum galaxy survey by the Square Kilometre Array (SKA), a photometric galaxy survey by Euclid and their combination, we forecast future constraints on primordial non-Gaussianity. We focus on the potential impact of local-type higher-order nonlinear parameters on the parameter estimation and particularly the confirmation of the inflationary consistency inequality. Nonstandard inflationary models, such as multifield models, introduce the scale-dependent stochastic clustering of galaxies on large scales, which is a unique probe of mechanism for generating primordial density fluctuations. Our Fisher matrix analysis indicates that a deep and wide survey provided by SKA is more advantageous to constrain $τ_{\rm NL}$, while Euclid has a strong constraining power for $f_{\rm NL}$ due to the redshift information, suggesting that the joint analysis between them is quite essential to break the degeneracy between the nonlinear parameters. The combination of the full SKA and Euclid will achieve the precision level needed to confirm the consistency inequality even for $f_{\rm NL}\approx 1.5$ and $τ_{\rm NL}\approx 17$, though it is still hard for a single survey to confirm it when $f_{\rm NL}\lesssim 2.7$.

preprint2016arXiv

Rapid Merger of Binary Primordial Black Holes: an Implication for GW150914

We propose a new scenario for the evolution of a binary of primordial black holes (PBHs). We consider a dynamical friction by ambient dark matter, scattering of dark matter particles with a highly eccentric orbit besides the standard two-body relaxation process to refill the loss cone, and interaction between the binary and a circumbinary disk, assuming that PBHs do not constitute the bulk of dark matter. Binary PBHs lose the energy and angular momentum by these processes, which could be sufficiently efficient for a typical configuration. Such a binary coalesces due to the gravitational wave emission in a time scale much shorter than the age of the universe. We estimate the density parameter of the resultant gravitational wave background. Astrophysical implication concerning the formation of intermediate-mass to supermassive black holes is also discussed.

preprint2016arXiv

SKA-Japan Pulsar Science with the Square Kilometre Array

The Square Kilometre Array will revolutionize pulsar studies with its wide field-of-view, wide-band observation and high sensitivity, increasing the number of observable pulsars by more than an order of magnitude. Pulsars are of interest not only for the study of neutron stars themselves but for their usage as tools for probing fundamental physics such as general relativity, gravitational waves and nuclear interaction. In this article, we summarize the activity and interests of SKA-Japan Pulsar Science Working Group, focusing on an investigation of modified gravity theory with the supermassive black hole in the Galactic Centre, gravitational-wave detection from cosmic strings and binary supermassive black holes, a study of the physical state of plasma close to pulsars using giant radio pulses and determination of magnetic field structure of Galaxy with pulsar pairs.

preprint2016arXiv

Studying topological structure in the epoch of reionization with 3D-Minkowski functionals of 21cm line fluctuations

The brightness temperature of the redshifted 21cm line brings rich information on the Inter Galactic Medium (IGM) from the Cosmic Dawn and Epoch of Reionization (EoR). While the power spectrum is a useful tool to statistically investigate the 21cm signal, the 21cm brightness temperature field is highly non-Gaussian, and the power spectrum is inadequate to characterize the non-Gaussianity. The Minkowski Functionals (MFs) are promising tools to extract non-gaussian features of the 21cm signal and give topological information such as morphology of ionized bubbles. In this work, we study the 21cm line signal in detail with MFs. To promote understanding of basic features of the 21cm signal, we calculate the MFs of not only the hydrogen neutral fraction but the matter density and spin temperature, which contribute to the brightness temperature fluctuations. We find that the structure of the brightness temperature depends mainly on the ionized fraction and the spin temperature at late and early stages of the EoR, respectively. Further, we investigate the redshift evolution of the MFs at $7 < z < 20$. Then, we find that, after the onset of reionization, the MFs reflect mainly the ionized bubble property. In addition, the MFs are sensitive to model parameters which are related to the topology of ionized bubbles and we consider the possibility of constraining the parameters by the future 21cm signal observations.

preprint2015arXiv

Euclid & SKA Synergies

Over the past few years two of the largest and highest fidelity experiments conceived have been approved for construction: Euclid is an ESA M-Class mission that will map three-quarters of the extra galactic sky with Hubble Space Telescope resolution optical and NIR imaging, and NIR spectroscopy, its scientific aims (amongst others) are to create a map of the dark Universe and to determine the nature of dark energy. The Square Kilometre Array (SKA) has similar scientific aims (and others) using radio wavelength observations. The two experiments are synergistic in several respects, both through the scientific objectives and through the control of systematic effects. SKA Phase-1 and Euclid will be commissioned on similar timescales offering an exciting opportunity to exploit synergies between these facilities.

preprint2015arXiv

Magnetic field spectrum at cosmological recombination revisited

If vector type perturbations are present in the primordial plasma before recombination, the generation of magnetic fields is known to be inevitable through the Harrison mechanism. In the context of the standard cosmological perturbation theory, non-linear couplings of first-order scalar perturbations create second-order vector perturbations, which generate magnetic fields. Here we reinvestigate the generation of magnetic fields at second-order in cosmological perturbations on the basis of our previous study, and extend it by newly taking into account the time evolution of purely second-order vector perturbations with a newly developed second-order Boltzmann code. We confirm that the amplitude of magnetic fields from the product-terms of the first-order scalar modes is consistent with the result in our previous study. However, we find, both numerically and analytically, that the magnetic fields from the purely second-order vector perturbations partially cancel out the magnetic fields from one of the product-terms of the first-order scalar modes, in the tight coupling regime in the radiation dominated era. Therefore, the amplitude of the magnetic fields on small scales, $k \gtrsim 10~h{\rm Mpc}^{-1}$, is smaller than the previous estimates. The amplitude of the generated magnetic fields at cosmological recombination is about $B_{\rm rec} =5.0\times 10^{-24}~{\rm Gauss}$ on $k = 5.0 \times 10^{-1}~h{\rm Mpc}^{-1}$. Finally, we discuss the reason of the discrepancies that exist in estimates of the amplitude of magnetic fields among other authors.

preprint2015arXiv

Magnetic Field Tomography in Nearby Galaxies with the Square Kilometre Array

Magnetic fields play an important role in shaping the structure and evolution of the interstellar medium (ISM) of galaxies, but the details of this relationship remain unclear. With SKA1, the 3D structure of galactic magnetic fields and its connection to star formation will be revealed. A highly sensitive probe of the internal structure of the magnetoionized ISM is the partial depolarization of synchrotron radiation from inside the volume. Different configurations of magnetic field and ionized gas within the resolution element of the telescope lead to frequency-dependent changes in the observed degree of polarization. The results of spectro-polarimetric observations are tied to physical structure in the ISM through comparison with detailed modeling, supplemented with the use of new analysis techniques that are being actively developed and studied within the community such as Rotation Measure Synthesis. The SKA will enable this field to come into its own and begin the study of the detailed structure of the magnetized ISM in a sample of nearby galaxies, thanks to its extraordinary wideband capabilities coupled with the combination of excellent surface brightness sensitivity and angular resolution.

preprint2015arXiv

Overview of Complementarity and Synergy with Other Wavelengths in Cosmology in the SKA era

We give an overview of complementarity and synergy in cosmology between the Square Kilometre Array and future survey projects in other wavelengths. In the SKA era, precision cosmology will be limited by systematic errors and cosmic variance, rather than statistical errors. However, combining and/or cross-correlating multi-wavelength data, from the SKA to the cosmic microwave background, optical/infrared and X-ray, substantially reduce these limiting factors. In this chapter, we summarize future survey projects and show highlights of complementarity and synergy, which can be very powerful to probe major cosmological problems such as dark energy, modified gravity and primordial non-Gaussianity.

preprint2015arXiv

Sensitivity for 21cm Bispectrum from Epoch of Reionization

The 21cm line brightness temperature brings rich information about Epoch of Reionizaton (EoR) and high-$z$ universe (Cosmic Dawn and Dark Age). While the power spectrum is a useful tool to investigate the EoR signal statistically, higher-order statistics such as bispectrum are also valuable because the EoR signal is expected to be highly non-Gaussian. In this paper, we develop a formalism to calculate the bispectrum contributed from the thermal noise taking array configularion of telescopes into account, by extending a formalism for the power spectrum \cite{2006ApJ...653..815M}. We apply our formalism to the ongoing and future telescopes such as expanded Murchison Widefield Array (MWA), LOw Frequency ARray (LOFAR) and Square Kilometre Array (SKA). We find that expanded MWA does not have enough sensitivity to detect the bispectrum signal. On the other hand, LOFAR has better sensitivity and will be able to detect the peaks of the bispectrum as a function of redshift at large scales with comoving wavenumber $k \lesssim 0.03 {\rm Mpc}^{-1}$. The SKA has enough sensitivity to detect the bispectrum for much smaller scales $k \lesssim 0.3 {\rm Mpc}^{-1}$ and redshift $z \lesssim 20$

preprint2015arXiv

Studying 21cm power spectrum with one-point statistics

The redshifted 21cm line signal from neutral hydrogens is a promising tool to probe the cosmic dawn and the epoch of reionization (EoR). Ongoing and future low-frequency radio experiments are expected to detect its fluctuations, especially through the power spectrum. In this paper, we give a physical interpretation of the time evolution of the power spectrum of the 21cm brightness temperature fluctuations, which can be decomposed into dark matter density, spin temperature and neutral fraction of hydrogen fluctuations. From the one-point statistics of the fluctuations, such as variance and skewness, we find that the peaks and dips in the time evolution are deeply related to X-ray heating of the intergalactic gas, which controls the spin temperature. We suggest the skewness of the brightness temperature distribution is a key observable to identify the onset of X-ray heating.

preprint2015arXiv

X-raying extended emission and rapid decay of short gamma-ray bursts

Extended emission is a mystery in short gamma-ray bursts (SGRBs). By making time resolved spectral analyses of brightest nine events observed by ${\it Swift}$ XRT, we obviously classify the early X-ray emission of SGRBs into two types. One is the extended emission with exponentially rapid decay, which shows significant spectral softening during hundreds seconds since the SGRB trigger and is also detected by ${\it Swift}$-BAT. The other is a dim afterglow only showing power-law decay over $10^4$ s. The correlations between the temporal decay and spectral indices of the extended emissions are inconsistent with the $α$-$β$ correlation expected for the high-latitude curvature emission from a uniform jet. The observed too-rapid decay suggests the emission from a photosphere or a patchy surface, and manifests the stopping central engine via such as magnetic reconnection at the black hole.

preprint2014arXiv

CMB ISW-lensing bispectrum from cosmic strings

We study the effect of weak lensing by cosmic (super-)strings on the higher-order statistics of the cosmic microwave background (CMB). A cosmic string segment is expected to cause weak lensing as well as an integrated Sachs-Wolfe (ISW) effect, the so-called Gott-Kaiser-Stebbins (GKS) effect, to the CMB temperature fluctuation, which are thus naturally cross-correlated. We point out that, in the presence of such a correlation, yet another kind of the post-recombination CMB temperature bispectra, the ISW-lensing bispectra, will arise in the form of products of the auto- and cross-power spectra. We first present an analytic method to calculate the autocorrelation of the temperature fluctuations induced by the strings, and the cross-correlation between the temperature fluctuation and the lensing potential both due to the string network. In our formulation, the evolution of the string network is assumed to be characterized by the simple analytic model, the velocity-dependent one scale model, and the intercommutation probability is properly incorporated in orderto characterize the possible superstringy nature. Furthermore, the obtained power spectra are dominated by the Poisson-distributed string segments, whose correlations are assumed to satisfy the simple relations. We then estimate the signal-to-noise ratios of the string-induced ISW-lensing bispectra and discuss the detectability of such CMB signals from the cosmic string network. It is found that in the case of the smaller string tension, $Gμ\ll 10^{-7}$\,, the ISW-lensing bispectrum induced by a cosmic string network can constrain the string-model parameters even more tightly than the purely GKS-induced bispectrum in the ongoing and future CMB observations on small scales.

preprint2014arXiv

Constraining primordial non-Gaussianity via a multitracer technique with surveys by Euclid and Square Kilometre Array

We forecast future constraints on local-type primordial non-Gaussianity parameter $f_{\rm NL}$ with a photometric galaxy survey by Euclid, a continuum galaxy survey by Square Kilometre Array (SKA), and their combination. We derive a general expression for the covariance matrix of the power spectrum estimates of multiple tracers to show how the so-called multitracer technique improves constraints on $f_{\rm NL}$. In particular we clarify the role of the overlap fraction of multiple tracers and the division method of the tracers. Our Fisher matrix analysis indicates that stringent constraints of $σ(f_{\rm NL})\lesssim 1$ can be obtained even with a single survey, assuming five mass bins. When Euclid and SKA phase 1 (2) are combined, constraints on $f_{\rm NL}$ are improved to $σ(f_{\rm NL})= 0.61~(0.50)$.

preprint2014arXiv

Exploring the Intergalactic Magnetic Field by Means of Faraday Tomography

Unveiling the intergalactic magnetic field (IGMF) in filaments of galaxies is a very important and challenging subject in modern astronomy. In order to probe the IGMF from rotation measures (RMs) of extragalactic radio sources, we need to separate RMs due to other origins such as the source, intervening galaxies, and our Galaxy. In this paper, we discuss observational strategies for the separation by means of Faraday tomography (Faraday RM Synthesis). We consider an observation of a single radio source such as a radio galaxy or a quasar viewed through the Galaxy and the cosmic web. We then compare the observation with another observation of a neighbor source with a small angular separation. Our simulations with simple models of the sources suggest that it would be not easy to detect the RM due to the IGMF of order ~ 1 rad/m/m, an expected value for the IGMF through a single filament. Contrary to it, we find that the RM of at least ~10 rad/m/m could be detected with the SKA or its pathfinders/precursors, if we achieve selections of ideal sources. These results would be improved if we incorporate decomposition techniques such as RMCLEAN and QU-fitting. We discuss feasibility of the strategies for cases with complex Galactic emissions as well as with effects of observational noise and radio frequency interferences.

preprint2014arXiv

Faraday dispersion functions of galaxies

The Faraday dispersion function (FDF), which can be derived from an observed polarization spec- trum by Faraday rotation measure synthesis, is a profile of polarized emissions as a function of Faraday depth. We study intrinsic FDFs along sight lines through face-on, Milky-Way-like galaxies by means of a sophisticated galactic model incorporating 3D MHD turbulence, and investigate how much the FDF contains information intrinsically. Since the FDF reflects distributions of thermal and cosmic- ray electrons as well as magnetic fields, it has been expected that the FDF could be a new probe to examine internal structures of galaxies. We, however, find that an intrinsic FDF along a sight line through a galaxy is very complicated, depending significantly on actual configurations of turbulence. We perform 800 realizations of turbulence, and find no universal shape of the FDF even if we fix the global parameters of the model. We calculate the probability distribution functions of the standard deviation, skewness, and kurtosis of FDFs and compare them for models with different global pa- rameters. Our models predict that the presence of vertical magnetic fields and large scale-height of cosmic-ray electrons tend to make the standard deviation relatively large. Contrastingly, differences in skewness and kurtosis are relatively less significant.

preprint2014arXiv

Properties of Intrinsic Polarization Angle Ambiguities in Faraday Tomography

Faraday tomography is a powerful method to diagnose polarizations and Faraday rotations along the line of sight. The quality of Faraday tomography is, however, limited by several conditions. Recently, it is reported that Faraday tomography indicates false signals in some specific situations. In this paper, we systematically investigate the condition of the appearance of false signals in Faraday tomography. We study this by pseudo-observing two sources within a beam, and change in the intrinsic polarization angles, rotation measures, intensities, and frequency coverage. We find that false signals arise when rotation measure between the sources is less than 1.5 times the full width at half maximum of the rotation measure spread function. False signals also depend on the intensity ratio between the sources and are reduced for large ratio. On the other hand, the appearance of false signals does not depend on frequency coverage, meaning that the uncertainty should be correctly understood and taken into consideration even with future wide-band observations such as Square Kilometer Array (SKA).

preprint2014arXiv

Short Gamma Ray Burst Formation Rate from BATSE data using E_p-L_p correlation and the minimum gravitational wave event rate of coalescing compact binary

Using 72 Short Gamma Ray Bursts (SGRBs) with well determined spectral data observed by BATSE, we determine their redshift and the luminosity by applying $E_p$--$L_p$ correlation for SGRBs found by \cite{tsutsui13}. For 53 SGRBs with the observed flux brighter than $4 \times 10^{-6}~{\rm erg~cm^{-2}s^{-1}}$, the cumulative redshift distribution up to $z=1$ agrees well with that of 22 {\it Swift}~SGRBs. This suggests that the redshift determination by the $E_p$--$L_p$ correlation for SGRBs works well. The minimum event rate at $z=0$ is estimated as $ρ_{SGRB}(0) = 6.3_{-3.9}^{+3.1} \times 10^{-10}~{\rm events~Mpc^{-3}yr^{-1}}$ so that the minimum beaming angle is $0.6^\circ-7.8^\circ$ assuming the merging rate of $10^{-7}-4\times 10^{-6}~{\rm events~Mpc^{-3}yr^{-1}}$ suggested from the binary pulsar data. Interestingly, this angle is consistent with that for SGRB130603B of $\sim 4^\circ-8^\circ$\citep{fong13b}. On the other hand, if we assume the beaming angle of $\sim 6^\circ$ suggested from four SGRBs with the observed value of beaming angle, the minimum event rate including off-axis SGRBs is estimated as $ρ_{SGRB,all}^{min}(0)=1.15_{-0.71}^{+0.57}\times 10^{-7}~{\rm events~Mpc^{-3}yr^{-1}}$. If SGRBs are induced by coalescence of binary neutron stars (NSs) and/or black holes (BHs), this event rate leads to the minimum gravitational-wave detection rate of $\rm 3.9_{-2.4}^{+1.9} (152_{-94}^{+75})~events~y^{-1}$ for NS-NS (NS-BH) binary, respectively, by a worldwide network with KAGRA, advanced-LIGO, advanced-Virgo, and GEO.

preprint2013arXiv

Constraints on the multi-lognormal magnetic fields from the observations of the cosmic microwave background and the matter power spectrum

Primordial magnetic fields (PMFs), which were generated in the early Universe before recombination, affect the motion of plasma and then the cosmic microwave background and the matter power spectrum. We consider constraints on PMFs with a characteristic correlation length from the observations of the anisotropies of the cosmic microwave background and the matter power spectrum. The spectrum of PMFs is modeled with multi-lognormal distributions, rather than power-law distribution, and we derive constraints on the strength $|\mathbf{B}_k|$ at each wave number $k$ along with the standard cosmological parameters in the flat Universe and the foreground sources. We obtain upper bounds on the field strengths at $k=10^{-1}, 10^{-2}, 10^{-4}$, and $10^{-5}$ Mpc$^{-1}$ as 4.7 nG, 2.1 nG, 5.3 nG, and 10.9 nG ($2σ$ C.L.) respectively, while the field strength at $k=10^{-3} $Mpc$^{-1}$ turns out to have a finite value as $|\mathbf{B}_{k = 10^{-3}}| = 6.2 \pm 1.3 $ nG ($1σ$ C.L.). This finite value is attributed to the finite values of BB mode data at $\ell > 300$ obtained from the QUAD experiment. If we do not include the BB mode data, we obtain only the upper bound on $B_{k=10^{-3}}$.

preprint2013arXiv

Fisher Analysis on Wide-Band Polarimetry for Probing the Intergalactic Magnetic Field

We investigate the capability of ongoing radio telescopes for probing Faraday rotation measure (RM) due to the intergalactic magnetic field (IGMF) in the large-scale structure of the universe which is expected to be of order $O(1) {\rm rad/m^2}$. We consider polarization observations of a compact radio source such as quasars behind a diffuse source such as the Galaxy, and calculate Stokes parameters $Q$ and $U$ assuming a simple model of the Faraday dispersion functions with Gaussian shape. Then, we perform the Fisher analysis to estimate the expected errors in the model parameters from QU-fitting of polarization intensity, accounting for sensitivities and frequency bands of Australian Square Kilometer Array Pathfinder, Low Frequency Array, and the Giant Meterwave Radio Telescope. Finally, we examine the condition on the source intensities which are required to detect the IGMF. Our analysis indicates that the QU-fitting is promising for forthcoming wideband polarimetry to explore RM due to the IGMF in filaments of galaxies.

preprint2013arXiv

Forecast constraints on cosmic strings from future CMB, pulsar timing and gravitational wave direct detection experiments

We study future observational constraints on cosmic string parameters from various types of next-generation experiments: direct detection of gravitational waves (GWs), pulsar timing array, and the cosmic microwave background (CMB). We consider both GW burst and stochastic GW background searches by ground- and space-based interferometers as well as GW background detection in pulsar timing experiments. We also consider cosmic string contributions to the CMB temperature and polarization anisotropies. These different types of observations offer independent probes of cosmic strings and may enable us to investigate cosmic string properties if the signature is detected. In this paper, we evaluate the power of future experiments to constrain cosmic string parameters, such as the string tension Gmu, the initial loop size alpha, and the reconnection probability p, by performing Fisher information matrix calculations. We find that combining the information from the different types of observations breaks parameter degeneracies and provides more stringent constraints on the parameters. We also find future space-borne interferometers independently provide a highly precise determination of the parameters.

preprint2013arXiv

Lower Bounds on Magnetic Fields in Intergalactic Voids from Long-Term GeV-TeV Light Curves of the Blazar Mrk 421

Lower bounds are derived on the amplitude B of intergalactic magnetic fields (IGMFs) in the region between Galaxy and the blazar Mrk 421, from constraints on the delayed GeV pair-echo flux that are emitted by secondary e^-e^+ produced in γγinteractions between primary TeV gamma-rays and the cosmic infrared background. The distribution of galaxies mapped by the Sloan Digital Sky Survey shows that this region is dominated by a large intergalactic void. We utilize data from long-term, simultaneous GeV-TeV observations by the Fermi Large Area Telescope and the ARGO-YBJ experiment extending over 850 days. For an assumed value of B, we evaluate the daily GeV pair-echo flux expected from the TeV data, select the dates where this exceeds the Fermi 2-σsensitivity, compute the probability that this flux is excluded by the Fermi data for each date, and then combine the probabilities using the inverse normal method. Consequently, we exclude B < 10^{-20.5} G for a field coherence length of 1 kpc at 4- σlevel, as long as plasma instabilities are unimportant for cooling of the pair beam. This is much more significant than the 2-σbounds we obtained previously from observations of Mrk 501, by virtue of more extensive data from the ARGO-YBJ, as well as improved statistical analysis. Compared with most other studies of IGMF bounds, the evidence we present here for a non-zero IGMF is more robust as it does not rely on unproven assumptions on the primary TeV emission during unobserved periods.

preprint2013arXiv

Possible existence of Ep-Lp and Ep-Eiso correlations for Short Gamma-Ray Bursts with a factor 5 to 100 dimmer than those for Long Gamma-Ray Bursts

We analyzed correlations among the rest frame spectral peak energy $E_{\rm p}$, the observed frame 64ms peak isotropic luminosity $L_{\rm p}$ and the isotropic energy $E_{\rm iso}$ for 13 Short Gamma Ray Burst (SGRB) candidates having the measured redshift $z$, $T_{\rm 90}^{\rm obs}/(1+z)<2$ sec and well determined spectral parameters. A SGRB candidate is regarded as a misguided SGRB if it is located in the 3-$σ_{\rm int}$ dispersion region from the best-fit function of the $E_{\rm p}$--$E_{\rm iso}$ correlation for Long GRBs (LGRBs) while the others are regarded as secure SGRBs possibly from compact star mergers. Using 8 secure SGRBs out of 13 SGRB candidates, we tested whether $E_{\rm p}$--$E_{\rm iso}$ and $E_{\rm p}$--$L_{\rm p}$ correlations exist for SGRBs. We found that $E_{\rm p}$--$E_{\rm iso}$ correlation for SGRBs($E_{\rm iso} =10^{51.42 \pm 0.15}{\rm erg s^{-1}} ({E_{\rm p}}/{\rm 774.5 keV})^{1.58 \pm 0.28}$) seems to exist with the correlation coefficeint $r=0.91$ and chance probability $p=1.5\times10^{-3}$. We found also that the $E_{\rm p}$--$L_{\rm p}$ correlation for SGRBs($L_{\rm p} = 10^{52.29 \pm 0.066}{\rm erg s^{-1}} ({E_{\rm p}}/{\rm 774.5 keV})^{1.59 \pm 0.11}$) is tighter than $E_{\rm p}$--$E_{\rm iso}$ correlation since $r=0.98$ and $p=1.5\times10^{-5}$. Both correlations for SGRBs are dimmer than those of LGRBs for the same $E_{\rm p}$ by factors $\sim$100 ($E_{\rm p}$--$E_{\rm iso}$) and $\sim$ 5($E_{\rm p}$--$L_{\rm p}$). Applying the tighter $E_{\rm p}$--$L_{\rm p}$ correlation for SGRBs to 71 bright BATSE SGRBs, we found that pseudo redshift $z$ ranges from 0.097 to 2.258 with the mean $<z>$ of 1.05. The redshifts of SGRBs apparently cluster at lower redshift than those of LGRBs ($<z>\sim 2.2 $), which supports the merger scenario of SGRBs.

preprint2013arXiv

Type-I cosmic string network

We study the network of Type-I cosmic strings using the field-theoretic numerical simulations in the Abelian-Higgs model. For Type-I strings, the gauge field plays an important role, and thus we find that the correlation length of the strings is strongly dependent upon the parameter β, the ratio between the self-coupling constant of the scalar field and the gauge coupling constant, namely, β=λ/2e^2. In particular, if we take the cosmic expansion into account, the network becomes densest in the comoving box for a specific value of βfor β<1.

preprint2012arXiv

Femto-lensing due to a Cosmic String

We consider the femto-lensing due to a cosmic string. If a cosmic string with the deficit angle $Δ\sim 100$ [femto-arcsec] $\sim10^{-18}$ [rad] exists around the line of sight to a gamma-ray burst, we may observe characteristic interference patterns caused by gravitational lensing in the energy spectrum of the gamma-ray burst. This "femto-lensing" event was first proposed as a tool to probe small mass primordial black holes. In this paper, we propose use of the femto-lensing to probe cosmic strings with extremely small tension. Observability conditions and the event rate are discussed. Differences between the cases of a point mass and a cosmic string are presented.

preprint2012arXiv

Forecast constraints on cosmic string parameters from gravitational wave direct detection experiments

Gravitational waves (GWs) are one of the key signatures of cosmic strings. If GWs from cosmic strings are detected in future experiments, not only their existence can be confirmed but also their properties might be probed. In this paper, we study the determination of cosmic string parameters through direct detection of GW signatures in future ground-based GW experiments. We consider two types of GWs, bursts and the stochastic GW background, which provide us with different information about cosmic string properties. Performing the Fisher matrix calculation on the cosmic string parameters, such as parameters governing the string tension $Gμ$ and initial loop size $α$ and the reconnection probability $p$, we find that the two different types of GW can break degeneracies in some of these parameters and provide better constraints than those from each measurement.

preprint2012arXiv

Gamma-Ray Bursts are precise distance indicators similar to Type Ia Supernovae?

We estimate the distance modulus to long gamma-ray bursts (LGRBs) using the Type I Fundamental Plane, a correlation between the spectral peak energy $E_{\rm p}$, the peak luminosity $L_{\rm p}$, and the luminosity time $T_{\rm L}$ ($\equiv E_{\rm iso}/L_{\rm p}$ where $E_{\rm iso}$ is isotropic energy) for small Absolute Deviation from Constant Luminosity(ADCL). The Type I Fundamental Plane of LGRBs is calibrated using 8 LGRBs with redshift $z<1.4$. To avoid any assumption on the cosmological model, we use the distance modulus of 557 Type Ia supernovae (SNeIa) from the Union 2 sample. This calibrated Type I Fundamental Plane is used to measure the distance moduli to 9 high-redshift LGRBs with the mean error $\bar σ_μ=0.31$, which is comparable with that of SNe Ia $\bar σ_μ=0.26$ where $μ$ stands for the distance modulus. The Type I Fundamental Plane is so tight that our distance moduli have very small uncertainties. From those distance moduli, we obtained the constraint $Ω_{\rm M}=0.22\pm0.04$ for flat $Λ$CDM universe. Adding 9 LGRBs distance moduli ($z>1.4$) to 557 SNeIa distance moduli ($z<1.4$) significantly improves the constraint for non-flat $Λ$CDM universe from ($Ω_{\rm M}, Ω_{\rm Λ}$)=($0.29\pm0.10$, $0.76\pm0.13$) for SNeIa only to ($Ω_{\rm M}, Ω_{\rm Λ}$)=($0.23\pm0.06$, $0.68\pm0.08$) for SNeIa and 9 LGRBs.

preprint2012arXiv

Identifying Subclasses of Long Gamma-Ray Bursts with Cumulative Light Curve Morphology of Prompt Emissions

We argue a new classification scheme of long gamma-ray bursts (LGRBs) using the morphology of the cumulative light curve of the prompt emission. We parametrize the morphology by the absolute deviation from their constant luminosity ($ADCL$) and derive the value for 36 LGRBs which have spectropic redshifts, spectral parameters determined by the Band model, 1-second peak fluxes, fluences, and 64-msec resolution light curves whose peak counts are 10 times larger than background fluctuations. Then we devide the sample according to the value of ADCL into two groups ($ADCL < 0.17$ and $ADCL > 0.17$) and, for each group, derive the spectral peak energy $E_{\rm p}$ - peak luminosity $L_{\rm p}$ correlation and the Fundamental Plane of LGRBs, which is a correlation between the spectral peak energy $E_{\rm p}$, the luminosity time $T_{\rm L}$ ($\equiv E_{\rm iso}/L_{\rm p}$ where $E_{\rm iso}$ is isotropic energy) and the peak luminosity $L_{\rm p}$. We find that both of the correlations for both groups are statistically more significant compared with ones derived from all samples. The Fundamental Planes with small and large ADCL are given by $L_{\rm p}=10^{52.53\pm 0.01}(E_{\rm p}/10^{2.71}{\rm keV})^{1.84\pm 0.03} (T_{\rm L}/10^{0.86}{\rm sec})^{0.29\pm0.08}$ with $χ^2_ν=10.93/14$ and $L_{\rm p}=10^{52.98\pm0.08}(E_{\rm p}/10^{2.71}{\rm keV})^{1.82\pm 0.09} (T_{\rm L}/10^{0.86}{\rm sec})^{0.85\pm 0.27}$ with $χ^2_ν=7.58/8$, respectively. This fact implies the existence of subclasses of LGRBs characterized by the value of $ADCL$. Also there is a hint for the existence of the intermediate-$ADCL$ class which deviates from both fundamental planes. Both relations are so tight that our result provides a new accurate distance measurement scheme up to the high redshift universe.

preprint2012arXiv

Possible Origins of Dispersion of the Peak Energy--Brightness Correlations of Gamma-Ray Bursts

We collect and reanalyze about 200 GRB data of prompt-emission with known redshift observed until the end of 2009, and select 101 GRBs which were well observed to have good spectral parameters to determine the spectral peak energy ($E_p$), 1-second peak luminosity ($L_p$) and isotropic energy ($E_{\rm iso}$). Using our newly-constructed database with 101 GRBs, we first revise the $E_p$--$L_p$ and $E_p$--$E_{\rm iso}$ correlations. The correlation coefficients of the revised correlations are 0.889 for 99 degree of freedom for the $E_p$--$L_p$ correlation and 0.867 for 96 degree of freedom for the $E_p$--$E_{\rm iso}$ correlation. These values correspond to the chance probability of $2.18 \times 10^{-35}$ and $4.27 \times 10^{-31}$, respectively. It is a very important issue whether these tight correlations are intrinsic property of GRBs or caused by some selection effect of observations. In this paper, we examine how the truncation of the detector sensitivity affects the correlations, and we conclude they are surely intrinsic properties of GRBs. Next we investigate origins of the dispersion of the correlations by studying their brightness and redshift dependence. Here the brightness (flux or fluence) dependence would be regarded as an estimator of the bias due to the detector threshold. We find a weak fluence-dependence in the $E_p$--$E_{\rm iso}$ correlations and a redshift dependence in the $E_p$--$L_p$ correlation both with 2 $σ$ statistical level. These two effects may contribute to the dispersion of the correlations which is larger than the statistical uncertainty. We discuss a possible reason of these dependence and give a future prospect to improve the correlations.

preprint2012arXiv

Stochastic Gravitational Wave Background originating from Halo Mergers

The stochastic gravitational wave background (GWB) from halo mergers is investigated by a quasi-analytic method. The method we employ consists of two steps. The first step is to construct a merger tree by using the Extended Press-Schechter formalism or the Sheth & Tormen formalism, with Monte-Carlo realizations. This merger tree provides evolution of halo masses. From $N$-body simulation of two-halo mergers, we can estimate the amount of gravitational wave emission induced by the individual merger process. Therefore the second step is to combine this gravitaional wave emission to the merger tree and obtain the amplitude of GWB. We find $Ω_{GW}\sim 10^{-19}$ for $f\sim 10^{-17}-10^{-16}$ Hz, where $Ω_{GW}$ is the energy density of the GWB. It turns out that most of the contribution on the GWB comes from halos with masses below $10^{15} M_\odot$ and mergers at low redshift, i.e., $0<z<0.8$.

preprint2012arXiv

Weak lensing of CMB by cosmic (super-)strings

We study the effect of weak lensing by cosmic (super-)strings on the anisotropies of cosmic microwave background (CMB). In developing a method to calculate the lensing convergence field due to strings, and thereby temperature and polarization angular power spectra of CMB, we clarify how the nature of strings, characterized by the intercommuting probability, can influence these CMB anisotropies. Assuming that the power spectrum is dominated by Poisson-distributed string segments, we find that the convergence spectrum peaks at low multipoles and is mostly contributed from strings located at relatively low redshifts. As the intercommuting probability decreases, the spectra of the convergence and hence the lensed temperature and polarizations are gained because the number density of string segments becomes larger. An observationally important feature of the string-induced CMB polarizations is that the string-lensed spectra decay more slowly on small scales compared with primordial scalar perturbations from standard inflation.

preprint2011arXiv

CMB Bispectrum from Primordial Scalar, Vector and Tensor non-Gaussianities

We present an all-sky formalism for the Cosmic Microwave Background (CMB) bispectrum induced by the primordial non-Gaussianities not only in scalar but also in vector and tensor fluctuations. We find that the bispectrum can be formed in an explicitly rationally invariant way by taking into account the angular and polarization dependences of the vector and tensor modes. To demonstrate this and present how to use our formalism, we consider a specific example of the correlation between two scalars and a graviton as the source of non-Gaussianity. As a result, we show that the CMB reduced bispectrum of the intensity anisotropies is evaluated as a function of the multipole and the coupling constant between two scalars and a graviton denoted by $g_{tss}$; $|b_{\ell \ell \ell}| \sim \ell^{-4} \times 8 \times 10^{-18} |g_{tss}|$. By estimating the signal-to-noise ratio, we find that the constraint as $|g_{tss}| < 6$ will be expected from the PLANCK experiment.

preprint2011arXiv

Computation approach for CMB bispectrum from primordial magnetic fields

We present a detailed calculation of our previous short paper [M. Shiraishi, D. Nitta, S. Yokoyama, K. Ichiki, and K. Takahashi, Phys. Rev. D 82, 121302 (2010).] in which we have investigated a constraint on the magnetic field strength through comic microwave background temperature bispectrum of vector modes induced from primordial magnetic fields. By taking into account full angular dependence of the bispectrum with spin spherical harmonics and Wigner symbols, we explicitly show that the cosmic microwave background bispectrum induced from the statistical-isotropic primordial vector fluctuations can be also described as an angle-averaged form in the rotationally invariant way. We also study the cases with different spectral indices of the power spectrum of the primordial magnetic fields.

preprint2011arXiv

Constraint on the primordial vector mode and its magnetic field generation from seven-year Wilkinson Microwave Anisotropy Probe Observations

A primordial vector mode and its associated magnetic field generation are investigated. Firstly, we put an observational constraint on the amount of the primordial vector mode from the seven-year WMAP data. The constraint is found as $r_v \lesssim -\frac{r}{40}+0.012$, where $r_v$ and $r$ are the amounts of vector and tensor perturbation amplitudes with respect to the scalar one, respectively. Secondly, we calculate the spectrum of magnetic fields inevitably created from the primordial vector mode, given the constraint on $r_v$. It is found that the maximum amount of magnetic fields generated from the vector mode is given by $B\lesssim 10^{-22}{\rm G} (\frac{r_v}{0.012})^{1/2} (\frac{k}{0.002})^{(n_v+1)/2}$ with $n_v$ being a spectral index of the vector mode. We find a non-trivial cancellation of the magnetic field generation in the radiation dominated era, which creates a characteristic cut off in the magnetic field spectrum around $k\approx 1.0$ Mpc$^{-1}$.

preprint2011arXiv

Cosmic microwave background bispectrum of tensor passive modes induced from primordial magnetic fields

If the seed magnetic fields exist in the early Universe, tensor components of their anisotropic stresses are not compensated prior to neutrino decoupling and the tensor metric perturbations generated from them survive passively. Consequently, due to the decay of these metric perturbations after recombination, the so-called integrated Sachs-Wolfe effect, the large-scale fluctuations of CMB radiation are significantly boosted. This kind of CMB anisotropy is called the "tensor passive mode." Because these fluctuations deviate largely from the Gaussian statistics due to the quadratic dependence on the strength of the Gaussian magnetic field, not only the power spectrum but also the higher-order correlations have reasonable signals. With these motives, we compute the CMB bispectrum induced by this mode. When the magnetic spectrum obeys a nearly scale-invariant shape, we obtain an estimation of a typical value of the normalized reduced bispectrum as $\ell_1(\ell_1 + 1)\ell_3(\ell_3+1)|b_{\ell_1\ell_2\ell_3}| \sim (130-6) \times 10^{-16} (B_{1 \rm Mpc} / 4.7 {\rm nG})^6$ depending on the energy scale of the magnetic field production from $10^{14}$GeV to $10^3$GeV. Here, $B_{1 {\rm Mpc}}$ is the strength of the primordial magnetic field smoothed on $1 {\rm Mpc}$. From the above estimation and the current observational constraint on the primordial non-Gaussianity, we get a rough constraint on the magnetic field strength as $B_{1 {\rm Mpc}} < 2.6 - 4.4 {\rm nG}$.

preprint2011arXiv

Effects of a primordial magnetic field with log-normal distribution on the cosmic microwave background

We study the effect of primordial magnetic fields (PMFs) on the anisotropies of the cosmic microwave background (CMB). We assume the spectrum of PMFs is described by log-normal distribution which has a characteristic scale, rather than power-law spectrum. This scale is expected to reflect the generation mechanisms and our analysis is complementary to previous studies with power-law spectrum. We calculate power spectra of energy density and Lorentz force of the log-normal PMFs, and then calculate CMB temperature and polarization angular power spectra from scalar, vector, and tensor modes of perturbations generated from such PMFs. By comparing these spectra with WMAP7, QUaD, CBI, Boomerang, and ACBAR data sets, we find that the current CMB data set places the strongest constraint at $k\simeq 10^{-2.5}$ Mpc$^{-1}$ with the upper limit $B\lesssim 3$ nG.

preprint2011arXiv

Lower Bounds on Intergalactic Magnetic Fields from Simultaneously Observed GeV-TeV Light Curves of the Blazar Mrk 501

We derive lower bounds on intergalactic magnetic fields (IGMFs) from upper limits on the pair echo emission from the blazar Mrk 501, that is, delayed GeV emission from secondary $e^{-}e^{+}$ pairs produced via interactions of primary TeV gamma rays with the cosmic infrared background. Utilizing only simultaneous GeV-TeV light curves observed by VERITAS, MAGIC and {\it Fermi}-LAT during a multiwavelength campaign in 2009 that included a TeV flare, bounds are deduced on the IGMF strength of $B \gtrsim 10^{-20} {\rm G}$ at 90% confidence level for a field coherence length of 1 kpc. Since our analysis is based firmly on the observational data alone and nearly free of assumptions concerning the primary TeV flux in unobserved periods or spectral bands, our evaluation of the pair echo flux is conservative and the evidence for a non-zero IGMF is more robust compared to previous studies.

preprint2011arXiv

Primordial magnetic fields generated by the non-adiabatic fluctuations at pre-recombination era

In the pre-recombination era, cosmological density fluctuations can naturally generate magnetic fields through Thomson scatterings. In previous studies, only the magnetic field generation from the initially-adiabatic fluctuations has been considered. Here we investigate the generation of cosmological magnetic fields sourced by the primordial non-adiabatic fluctuations based on the cosmological perturbation theory, using the tight-coupling approximations between photon and baryon fluids. It is found that the magnetic fields from the non-adiabatic fluctuations can arise at the first-order expansion of the tight coupling approximation. This result is in contrast to the case of adiabatic initial fluctuations, where the magnetic fields can be generated only at the second-order. In a general case where the primordial density perturbations contain small non-adiabatic fluctuations on the top of the dominant adiabatic ones, we show that the leading source of magnetic fields is given by the second-order coupling of the adiabatic and non-adiabatic fluctuations. We calculate the power spectrum of the generated magnetic fields when the non-adiabatic fluctuations have a blue power spectrum, which has been suggested by recent cosmological observations.

preprint2011arXiv

The Improved Ep-TL-Lp Diagram and a Robust Regression Method

The accuracy and reliability of gamma-ray bursts (GRBs) as distance indicators are strongly restricted by their systematic errors which are larger than statistical errors. These systematic errors might come from either intrinsic variations of GRBs, or systematic errors in observations. In this paper, we consider the possible origins of systematic errors in the following observables, (i) the spectral peak energies (Ep) estimated by Cut-off power law (CPL) function, (ii) the peak luminosities (Lp) estimated by 1 second in observer time. Removing or correcting them, we reveal the true intrinsic variation of the Ep-TL-Lp relation of GRBs. Here TL is the third parameter of GRBs defined as TL ~ Eiso / Lp. Not only the time resolution of Lp is converted from observer time to GRB rest frame time, the time resolution with the largest likelihood is sought for. After removing obvious origin of systematic errors in observation mentioned above, there seems to be still remain some outliers. For this reason, we take account another origin of the systematic error as below, (iii) the contamination of short GRBs or other populations. To estimate the best fit parameters of the Ep-TL-Lp relations from data including outliers, we develop a new method which combine robust regression and an outlier identification technique. Using our new method for 18 GRBs with σEp/Ep < 0.1, we detect 6 outliers and find the Ep-TL-Lp relation become the tightest around 3 second.

preprint2010arXiv

Analytic formulae of the CMB bispectra generated from non-Gaussianity in the tensor and vector perturbations

We present a complete set of formulae for calculating the bispectra of CMB temperature and polarization anisotropies generated from non-Gaussianity in the vector and tensor mode perturbations. In the all-sky analysis it is found that the bispectrum formulae for the tensor and vector-mode non-Gaussianity formally take complicated forms compared to the scalar mode one because the photon transfer functions in the tensor and vector modes depend on the azimuthal angle between the direction of the wave number vector of the photon's perturbation and that of the line of sight. We demonstrate that flat-sky approximations remove this difficulty because this kind of azimuthal angle dependence apparently vanishes in the flat-sky limit. Through the flat-sky analysis, we also find that the vector or tensor bispectrum of $B$-mode polarization vanishes in the squeezed limit, unless the cosmological parity is violated at the nonlinear level.

preprint2010arXiv

Analytical model for CMB temperature angular power spectrum from cosmic (super-)strings

We present a new analytical method to calculate the small angle CMB temperature angular power spectrum due to cosmic (super-)string segments. In particular, using our method, we clarify the dependence on the intercommuting probability $P$. We find that the power spectrum is dominated by Poisson-distributed string segments. The power spectrum for a general value of $P$ has a plateau on large angular scales and shows a power-law decrease on small angular scales. The resulting spectrum in the case of conventional cosmic strings is in very good agreement with the numerical result obtained by Fraisse et al.. Then we estimate the upper bound on the dimensionless tension of the string for various values of $P$ by assuming that the fraction of the CMB power spectrum due to cosmic (super-)strings is less than ten percents at various angular scales up to $\ell=2000$. We find that the amplitude of the spectrum increases as the intercommuting probability. As a consequence, strings with smaller intercommuting probabilities are found to be more tightly constrained.

preprint2010arXiv

Can we probe the Lorentz factor of gamma-ray bursts from GeV-TeV spectra integrated over internal shocks?

We revisit the high-energy spectral cutoff originating from the electron-positron pair creation in the prompt phase of gamma-ray bursts (GRBs) with numerical and analytical calculations. We show that the conventional exponential and/or broken power law cutoff should be drastically modified to a shallower broken power-law in practical observations that integrate emissions from different internal shocks. Since the steepening is tiny for observations, this "smearing" effect can generally reduce the previous estimates of the Lorentz factor of the GRB outflows. We apply our formulation to GRB 080916C, recently detected by the Large Area Telescope detector on the Fermi satellite, and find that the minimum Lorentz factor can be ~600 (or even smaller values), which is below but consistent with the previous result of ~900. Observing the steepening energy (so-called "pair-break energy") is crucial to diagnose the Lorentz factor and/or the emission site in the future observations, especially current and future Cherenkov telescopes such as MAGIC, VERITAS, and CTA.

preprint2010arXiv

Cosmic microwave background bispectrum of vector modes induced from primordial magnetic fields

We calculate the cosmic microwave background bispectrum of vector modes induced from primordial magnetic fields. We take into account the full angular dependence of the bispectrum and discuss the amplitude and also the shape of the bispectrum in $\ell$ space. In the squeezed limit as $\ell_1 = \ell_2 \gg \ell_3$, we estimate a typical values of the normalized reduced bispectrum as $\ell_1(\ell_1 + 1)\ell_3(\ell_3+1)|b_{\ell_1\ell_2\ell_3}| \sim 2 \times 10^{-19}$, for the strength of the primordial magnetic field smoothed on $1 {\rm Mpc}$ scale $B_{1 \rm Mpc} = 4.7 \rm nG$ assuming a nearly scale-invariant spectrum of magnetic fields. We find that a constraint on the magnetic field strength will be roughly estimated as $B_{1 \rm Mpc} < 10 {\rm nG}$ if PLANCK places a limit on the nonlinearity parameter of local-type configuration as $|f^{\rm local}_{\rm NL}| < 5$.

preprint2010arXiv

Generation of Seed Magnetic Fields in Primordial Supernova Remnants

Origin of the magnetic field ubiquitous in the Universe is studied based on the Biermann mechanism, which is expected to work in the non-barotropic region. We perform a series of two-dimensional MHD simulations of the first generation supernova remnant (SNR) expanding in the inhomogeneous interstellar matter (ISM) and study the Biermann mechanism working in the interior of the SNR. Especially, we pay attention to the relaxation process of electron and ion temperatures via the Coulomb interaction. In the early SNR in which the electron temperature is much lower than the ion temperature, the Biermann mechanism is ineffective, since the gradient of electron pressure is small. Magnetic fields begin to be generated just behind the shock front when the electron temperature is sufficiently relaxed. Assuming the explosion energy of 10^52 erg, the total magnetic energy generated reaches about 10^26 erg and does not depend strongly on the parameters of either SNR or ISM. Analytic expression to estimate the magnetic total energy is presented and it is shown this agrees well with our numerical results. Finally we evaluate the expected amplitude of magnetic fields in protogalaxies as ~10^-19 G, which is sufficient for seed fields of the observed galactic magnetic fields.

preprint2010arXiv

Intrisic Dispersion of Correlations among Ep, Lp, and Eiso of Gamma Ray Bursts depends on the quality of Data Set

We reconsider correlations among the spectral peak energy ($E_p$), 1-second peak luminosity ($L_p$) and isotropic energy (\Eiso), using the database constructed by \citet{yonetoku10} which consists of 109 Gamma-Ray Bursts (GRBs) whose redshifts are known and $E_p$, $L_p$ and \Eiso are well determined. We divide the events into two groups by their data quality. One (gold data set) consists of GRBs with peak energies determined by the Band model with four free parameters. On the other hand, GRBs in the other group (bronze data set) have relatively poor energy spectra so that their peak energies were determined by the Band model with fixed spectral index (i.e. three free parameters) or by the Cut-off power law (CPL) model with three free parameters. Using only the gold data set we found the intrinsic dispersion in $\log L_p$ ($=σ_{\rm int}$) is 0.13 and 0.22 for \tsutsui correlation ($T_L \equiv E_{\rm iso}/L_p$) and \yonetoku correlation, respectively. We also find that GRBs in the bronze data set have systematically larger $E_p$ than expected by the correlations constructed with the gold data set. This means that the intrinsic dispersion of correlations among $E_p$, $L_p$, and \Eiso of GRBs depends on the quality of data set. At present, using \tsutsui correlation with gold data set, we would be able to determine the luminosity distance with $\sim 16%$ error, which might be useful to determine the nature of the dark energy at high redshift $z > 3$.

preprint2010arXiv

Probing Early Cosmic Magnetic Fields through Pair Echos from High-Redshift GRBs

We discuss the expected properties of pair echo emission from gamma-ray bursts (GRBs) at high redshifts ($z \gtrsim 5$), their detectability, and the consequent implications for probing intergalactic magnetic fields (IGMFs) at early epochs. Pair echos comprise inverse Compton emission by secondary electron-positron pairs produced via interactions between primary gamma-rays from the GRB and low-energy photons of the diffuse intergalactic radiation, arriving with a time delay that depends on the nature of the intervening IGMFs. At sufficiently high $z$, the IGMFs are unlikely to have been significantly contaminated by astrophysical outflows, and the relevant intergalactic radiation may be dominated by the well-understood cosmic microwave background (CMB). Pair echoes from luminous GRBs at $z \sim 5-10$ may be detectable by future facilities such as the Cherenkov Telescope Array or the Advanced Gamma-ray Imaging System, as long as the GRB primary emission extends to multi-TeV energies, the comoving IGMFs at these redshifts are $B \sim 10^{-16}-10^{-15}~{\rm Gauss}$, and the non-CMB component of the diffuse intergalactic radiation is relatively low. Observations of pair echos from high-$z$ GRBs can provide a unique, in-situ probe of weak IGMFs during the epochs of early structure formation and cosmic reionization.

preprint2010arXiv

Skewness in CMB temperature fluctuations from curved cosmic (super-)strings

We compute the one-point probability distribution function of small-angle cosmic microwave background temperature fluctuations due to curved cosmic (super-)strings with a simple model of string network by performing Monte Carlo simulations. Taking into account of the correlation between the curvature and the velocity of string segments, there appear non-Gaussian features, specifically non-Gaussian tails and a skewness, in the one-point pdf. The obtained sample skewness for the conventional field-theoretic cosmic strings is $g_1\approx -0.14$, which is consistent with the result reported by Fraisse et al.. We also discuss the dependence of the pdf on the intercommuting probability. We find that the standard deviation of the Gaussian part increases and non-Gaussian features are suppressed as the intercommuting probability decreases. For sufficiently small intercommuting probability, the skewness is given by $\lesssim$ (a\ few) $\times 10^{-2}$.

preprint2010arXiv

Systematic Study of Gravitational Waves from Galaxy Merger

A systematic study of gravitational waves from galaxy mergers, through N-body simulations, was performed. In particular, we investigated the relative importance of galaxy components (disk, bulge and halo) and effects of initial relative velocity, relative angular momentum and mass ratio of the galaxies. We found that the features of light curve of gravitational waves, such as peak width and luminosity, are reliably simulated with particle numbers larger than ~10^4. Dominant contribution to gravitational wave emission came from the halo component, while peak luminosity amounted to 10^31 erg/sec for the collision of two halos with mass 3.8 x10^12Msun/h. We also found that the initial relative velocity in the direction of the initial separation did not significantly affect gravitational wave emission, while the initial relative angular momentum broadened the peak width and suppressed the luminosity. Mass dependence of the peak luminosity was also investigated, and we obtained evidence that the luminosity is proportional to the cubic mass when the scaling relation is satisfied. This behavior was considered by a simple analysis.

preprint2009arXiv

Non-Gaussianity in Cosmic Microwave Background Temperature Fluctuations from Cosmic (Super-)Strings

We compute analytically the small-scale temperature fluctuations of the cosmic microwave background from cosmic (super-)strings and study the dependence on the string intercommuting probability $P$. We develop an analytical model which describes the evolution of a string network and calculate the numbers of string segments and kinks in a horizon volume. Then we derive the probability distribution function (pdf) which takes account of finite angular resolution of observation. The resultant pdf consists of a Gaussian part due to frequent scatterings by long string segments and a non-Gaussian tail due to close encounters with kinks. The dispersion of the Gaussian part is reasonably consistent with that obtained by numerical simulations by Fraisse et al.. On the other hand, the non-Gaussian tail contains two phenomenological parameters which are determined by comparison with the numerical results for P=1. Extrapolating the pdf to the cases with $P<1$, we predict that the non-Gaussian feature is suppressed for small $P$.

preprint2009arXiv

Possible Effects of Pair Echoes on Gamma-Ray Burst Afterglow Emission

High-energy emission from gamma-ray bursts (GRBs) is widely expected but had been sparsely observed until recently when the Fermi satellite was launched. If >TeV gamma rays are produced in GRBs and can escape from the emission region, they are attenuated by the cosmic infrared background photons, leading to regeneration of GeV-TeV secondary photons via inverse-Compton scattering. This secondary emission can last for a longer time than the duration of GRBs, and it is called a pair echo. We investigate how this pair echo emission affects spectra and light curves of high energy afterglows, considering not only prompt emission but also afterglow as the primary emission. Detection of pair echoes is possible as long as the intergalactic magnetic field (IGMF) in voids is weak. We find (1) that the pair echo from the primary afterglow emission can affect the observed high-energy emission in the afterglow phase after the jet break, and (2) that the pair echo from the primary prompt emission can also be relevant, but only when significant energy is emitted in the TeV range, typically E_{gamma, >0.1 TeV} > (Y/(1+Y)) epsilon_e E_k. Even non-detections of the pair echoes could place interesting constraints on the strength of IGMF. The more favorable targets to detect pair echoes may be the "naked" GRBs without conventional afterglow emission, although energetic naked GRBs would be rare. If the IGMF is weak enough, it is predicted that the GeV emission extends to >30-300 s.

preprint2008arXiv

Probing Intergalactic Magnetic Fields in the GLAST Era through Pair Echo Emission from TeV Blazars

More than a dozen blazars are known to be emitters of multi-TeV gamma rays, often with strong and rapid flaring activity. By interacting with photons of the cosmic microwave and infrared backgrounds, these gamma rays inevitably produce electron-positron pairs, which in turn radiate secondary inverse Compton gamma rays in the GeV-TeV range with a characteristic time delay that depends on the properties of the intergalactic magnetic field (IGMF). For sufficiently weak IGMF, such "pair echo" emission may be detectable by the Gamma-ray Large Area Space Telescope (GLAST), providing valuable information on the IGMF. We perform detailed calculations of the time-dependent spectra of pair echos from flaring TeV blazars such as Mrk 501 and PKS 2155-304, taking proper account of the echo geometry and other crucial effects. In some cases, the presence of a weak but non-zero IGMF may enhance the detectability of echos. We discuss the quantitative constraints that can be imposed on the IGMF from GLAST observations, including the case of non-detections.

preprint2005arXiv

Gravity is controlled by cosmological constant

We discuss a Randall-Sundrum-type two D-braneworld model in which D-branes possess different values of the tensions from those of the charges, and derive an effective gravitational equation on the branes. As a consequence, the Einstein-Maxwell theory is realized together with the non-zero cosmological constant. Here an interesting point is that the effective gravitational constant is proportional to the cosmological constant. If the distance between two D-branes is appropriately tuned, the cosmological constant can have a consistent value with the current observations. From this result we see that, in our model, the presence of the cosmological constant is naturally explained by the presence of the effective gravitational coupling of the Maxwell field on the D-brane.

preprint2004arXiv

Gravity in Randall-Sundrum two D-brane model

We analyse Randall-Sundrum two D-brane model by linear perturbation and then consider the linearised gravity on the D-brane. The qualitative contribution from the Kaluza-Klein modes of gauge fields to the coupling to the gravity on the brane will be addressed. As a consequence, the gauge fields localised on the brane are shown not to contribute to the gravity on the brane at large distances. Although the coupling between gauge fields and gravity appears in the next order, the ordinary coupling cannot be realised.