Source author record

Y. Inoue

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

72works
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

72 published item(s)

preprint2023arXiv

Model-based cross-correlation search for gravitational waves from the low-mass X-ray binary Scorpius X-1 in LIGO O3 data

We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to balance sensitivity with computing cost. The search covered a range of gravitational-wave frequencies from 25Hz to 1600Hz, as well as ranges in orbital speed, frequency and phase determined from observational constraints. No significant detection candidates were found, and upper limits were set as a function of frequency. The most stringent limits, between 100Hz and 200Hz, correspond to an amplitude h0 of about 1e-25 when marginalized isotropically over the unknown inclination angle of the neutron star's rotation axis, or less than 4e-26 assuming the optimal orientation. The sensitivity of this search is now probing amplitudes predicted by models of torque balance equilibrium. For the usual conservative model assuming accretion at the surface of the neutron star, our isotropically-marginalized upper limits are close to the predicted amplitude from about 70Hz to 100Hz; the limits assuming the neutron star spin is aligned with the most likely orbital angular momentum are below the conservative torque balance predictions from 40Hz to 200Hz. Assuming a broader range of accretion models, our direct limits on gravitational-wave amplitude delve into the relevant parameter space over a wide range of frequencies, to 500Hz or more.

preprint2022arXiv

All-sky search for gravitational wave emission from scalar boson clouds around spinning black holes in LIGO O3 data

This paper describes the first all-sky search for long-duration, quasi-monochromatic gravitational-wave signals emitted by ultralight scalar boson clouds around spinning black holes using data from the third observing run of Advanced LIGO. We analyze the frequency range from 20~Hz to 610~Hz, over a small frequency derivative range around zero, and use multiple frequency resolutions to be robust towards possible signal frequency wanderings. Outliers from this search are followed up using two different methods, one more suitable for nearly monochromatic signals, and the other more robust towards frequency fluctuations. We do not find any evidence for such signals and set upper limits on the signal strain amplitude, the most stringent being $\approx10^{-25}$ at around 130~Hz. We interpret these upper limits as both an "exclusion region" in the boson mass/black hole mass plane and the maximum detectable distance for a given boson mass, based on an assumption of the age of the black hole/boson cloud system.

preprint2022arXiv

Constraints from LIGO O3 data on gravitational-wave emission due to r-modes in the glitching pulsar PSR J0537-6910

We present a search for continuous gravitational-wave emission due to r-modes in the pulsar PSR J0537-6910 using data from the LIGO-Virgo Collaboration observing run O3. PSR J0537-6910 is a young energetic X-ray pulsar and is the most frequent glitcher known. The inter-glitch braking index of the pulsar suggests that gravitational-wave emission due to r-mode oscillations may play an important role in the spin evolution of this pulsar. Theoretical models confirm this possibility and predict emission at a level that can be probed by ground-based detectors. In order to explore this scenario, we search for r-mode emission in the epochs between glitches by using a contemporaneous timing ephemeris obtained from NICER data. We do not detect any signals in the theoretically expected band of 86-97 Hz, and report upper limits on the amplitude of the gravitational waves. Our results improve on previous amplitude upper limits from r-modes in J0537-6910 by a factor of up to 3 and place stringent constraints on theoretical models for r-mode driven spin-down in PSR J0537-6910, especially for higher frequencies at which our results reach below the spin-down limit defined by energy conservation.

preprint2022arXiv

First joint observation by the underground gravitational-wave detector, KAGRA, with GEO600

We report the results of the first joint observation of the KAGRA detector with GEO600. KAGRA is a cryogenic and underground gravitational-wave detector consisting of a laser interferometer with three-kilometer arms, and located in Kamioka, Gifu, Japan. GEO600 is a British--German laser interferometer with 600 m arms, and located near Hannover, Germany. GEO600 and KAGRA performed a joint observing run from April 7 to 20, 2020. We present the results of the joint analysis of the GEO--KAGRA data for transient gravitational-wave signals, including the coalescence of neutron-star binaries and generic unmodeled transients. We also perform dedicated searches for binary coalescence signals and generic transients associated with gamma-ray burst events observed during the joint run. No gravitational-wave events were identified. We evaluate the minimum detectable amplitude for various types of transient signals and the spacetime volume for which the network is sensitive to binary neutron-star coalescences. We also place lower limits on the distances to the gamma-ray bursts analysed based on the non-detection of an associated gravitational-wave signal for several signal models, including binary coalescences. These analyses demonstrate the feasibility and utility of KAGRA as a member of the global gravitational-wave detector network.

preprint2022arXiv

Narrowband searches for continuous and long-duration transient gravitational waves from known pulsars in the LIGO-Virgo third observing run

Isolated neutron stars that are asymmetric with respect to their spin axis are possible sources of detectable continuous gravitational waves. This paper presents a fully-coherent search for such signals from eighteen pulsars in data from LIGO and Virgo's third observing run (O3). For known pulsars, efficient and sensitive matched-filter searches can be carried out if one assumes the gravitational radiation is phase-locked to the electromagnetic emission. In the search presented here, we relax this assumption and allow the frequency and frequency time-derivative of the gravitational waves to vary in a small range around those inferred from electromagnetic observations. We find no evidence for continuous gravitational waves, and set upper limits on the strain amplitude for each target. These limits are more constraining for seven of the targets than the spin-down limit defined by ascribing all rotational energy loss to gravitational radiation. In an additional search we look in O3 data for long-duration (hours-months) transient gravitational waves in the aftermath of pulsar glitches for six targets with a total of nine glitches. We report two marginal outliers from this search, but find no clear evidence for such emission either. The resulting duration-dependent strain upper limits do not surpass indirect energy constraints for any of these targets.

preprint2022arXiv

Search for anisotropic gravitational-wave backgrounds using data from Advanced LIGO and Advanced Virgo's first three observing runs

We report results from searches for anisotropic stochastic gravitational-wave backgrounds using data from the first three observing runs of the Advanced LIGO and Advanced Virgo detectors. For the first time, we include Virgo data in our analysis and run our search with a new efficient pipeline called {\tt PyStoch} on data folded over one sidereal day. We use gravitational-wave radiometry (broadband and narrow band) to produce sky maps of stochastic gravitational-wave backgrounds and to search for gravitational waves from point sources. A spherical harmonic decomposition method is employed to look for gravitational-wave emission from spatially-extended sources. Neither technique found evidence of gravitational-wave signals. Hence we derive 95\% confidence-level upper limit sky maps on the gravitational-wave energy flux from broadband point sources, ranging from $F_{α, Θ} < {\rm (0.013 - 7.6)} \times 10^{-8} {\rm erg \, cm^{-2} \, s^{-1} \, Hz^{-1}},$ and on the (normalized) gravitational-wave energy density spectrum from extended sources, ranging from $Ω_{α, Θ} < {\rm (0.57 - 9.3)} \times 10^{-9} \, {\rm sr^{-1}}$, depending on direction ($Θ$) and spectral index ($α$). These limits improve upon previous limits by factors of $2.9 - 3.5$. We also set 95\% confidence level upper limits on the frequency-dependent strain amplitudes of quasimonochromatic gravitational waves coming from three interesting targets, Scorpius X-1, SN 1987A and the Galactic Center, with best upper limits range from $h_0 < {\rm (1.7-2.1)} \times 10^{-25},$ a factor of $\geq 2.0$ improvement compared to previous stochastic radiometer searches.

preprint2022arXiv

Search for continuous gravitational wave emission from the Milky Way center in O3 LIGO--Virgo data

We present a directed search for continuous gravitational wave (CW) signals emitted by spinning neutron stars located in the inner parsecs of the Galactic Center (GC). Compelling evidence for the presence of a numerous population of neutron stars has been reported in the literature, turning this region into a very interesting place to look for CWs. In this search, data from the full O3 LIGO--Virgo run in the detector frequency band $[10,2000]\rm~Hz$ have been used. No significant detection was found and 95$\%$ confidence level upper limits on the signal strain amplitude were computed, over the full search band, with the deepest limit of about $7.6\times 10^{-26}$ at $\simeq 142\rm~Hz$. These results are significantly more constraining than those reported in previous searches. We use these limits to put constraints on the fiducial neutron star ellipticity and r-mode amplitude. These limits can be also translated into constraints in the black hole mass -- boson mass plane for a hypothetical population of boson clouds around spinning black holes located in the GC.

preprint2022arXiv

Search for continuous gravitational waves from 20 accreting millisecond X-ray pulsars in O3 LIGO data

Results are presented of searches for continuous gravitational waves from 20 accreting millisecond X-ray pulsars with accurately measured spin frequencies and orbital parameters, using data from the third observing run of the Advanced LIGO and Advanced Virgo detectors. The search algorithm uses a hidden Markov model, where the transition probabilities allow the frequency to wander according to an unbiased random walk, while the $\mathcal{J}$-statistic maximum-likelihood matched filter tracks the binary orbital phase. Three narrow sub-bands are searched for each target, centered on harmonics of the measured spin frequency. The search yields 16 candidates, consistent with a false alarm probability of 30% per sub-band and target searched. These candidates, along with one candidate from an additional target-of-opportunity search done for SAX J1808.4$-$3658, which was in outburst during one month of the observing run, cannot be confidently associated with a known noise source. Additional follow-up does not provide convincing evidence that any are a true astrophysical signal. When all candidates are assumed non-astrophysical, upper limits are set on the maximum wave strain detectable at 95% confidence, $h_0^{95\%}$. The strictest constraint is $h_0^{95\%} = 4.7\times 10^{-26}$ from IGR J17062$-$6143. Constraints on the detectable wave strain from each target lead to constraints on neutron star ellipticity and $r$-mode amplitude, the strictest of which are $ε^{95\%} = 3.1\times 10^{-7}$ and $α^{95\%} = 1.8\times 10^{-5}$ respectively. This analysis is the most comprehensive and sensitive search of continuous gravitational waves from accreting millisecond X-ray pulsars to date.

preprint2022arXiv

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in the Second and Third LIGO-Virgo Observing Runs

We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the $l=m=2$ mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the $l=2, m=1,2$ modes with a frequency of both once and twice the rotation frequency (dual harmonic). No evidence of GWs was found so we present 95\% credible upper limits on the strain amplitudes $h_0$ for the single harmonic search along with limits on the pulsars' mass quadrupole moments $Q_{22}$ and ellipticities $\varepsilon$. Of the pulsars studied, 23 have strain amplitudes that are lower than the limits calculated from their electromagnetically measured spin-down rates. These pulsars include the millisecond pulsars J0437\textminus4715 and J0711\textminus6830 which have spin-down ratios of 0.87 and 0.57 respectively. For nine pulsars, their spin-down limits have been surpassed for the first time. For the Crab and Vela pulsars our limits are factors of $\sim 100$ and $\sim 20$ more constraining than their spin-down limits, respectively. For the dual harmonic searches, new limits are placed on the strain amplitudes $C_{21}$ and $C_{22}$. For 23 pulsars we also present limits on the emission amplitude assuming dipole radiation as predicted by Brans-Dicke theory.

preprint2021arXiv

All-sky search for long-duration gravitational-wave bursts in the third Advanced LIGO and Advanced Virgo run

After the detection of gravitational waves from compact binary coalescences, the search for transient gravitational-wave signals with less well-defined waveforms for which matched filtering is not well-suited is one of the frontiers for gravitational-wave astronomy. Broadly classified into "short" $ \lesssim 1~$\,s and "long" $ \gtrsim 1~$\,s duration signals, these signals are expected from a variety of astrophysical processes, including non-axisymmetric deformations in magnetars or eccentric binary black hole coalescences. In this work, we present a search for long-duration gravitational-wave transients from Advanced LIGO and Advanced Virgo's third observing run from April 2019 to March 2020. For this search, we use minimal assumptions for the sky location, event time, waveform morphology, and duration of the source. The search covers the range of $2~\text{--}~ 500$~s in duration and a frequency band of $24 - 2048$ Hz. We find no significant triggers within this parameter space; we report sensitivity limits on the signal strength of gravitational waves characterized by the root-sum-square amplitude $h_{\mathrm{rss}}$ as a function of waveform morphology. These $h_{\mathrm{rss}}$ limits improve upon the results from the second observing run by an average factor of 1.8.

preprint2021arXiv

All-sky search for short gravitational-wave bursts in the third Advanced LIGO and Advanced Virgo run

This paper presents the results of a search for generic short-duration gravitational-wave transients in data from the third observing run of Advanced LIGO and Advanced Virgo. Transients with durations of milliseconds to a few seconds in the 24--4096 Hz frequency band are targeted by the search, with no assumptions made regarding the incoming signal direction, polarization or morphology. Gravitational waves from compact binary coalescences that have been identified by other targeted analyses are detected, but no statistically significant evidence for other gravitational wave bursts is found. Sensitivities to a variety of signals are presented. These include updated upper limits on the source rate-density as a function of the characteristic frequency of the signal, which are roughly an order of magnitude better than previous upper limits. This search is sensitive to sources radiating as little as $\sim$10$^{-10} M_{\odot} c^2$ in gravitational waves at $\sim$70 Hz from a distance of 10~kpc, with 50\% detection efficiency at a false alarm rate of one per century. The sensitivity of this search to two plausible astrophysical sources is estimated: neutron star f-modes, which may be excited by pulsar glitches, as well as selected core-collapse supernova models.

preprint2021arXiv

All-sky, all-frequency directional search for persistent gravitational-waves from Advanced LIGO's and Advanced Virgo's first three observing runs

We present the first results from an all-sky all-frequency (ASAF) search for an anisotropic stochastic gravitational-wave background using the data from the first three observing runs of the Advanced LIGO and Advanced Virgo detectors. Upper limit maps on broadband anisotropies of a persistent stochastic background were published for all observing runs of the LIGO-Virgo detectors. However, a broadband analysis is likely to miss narrowband signals as the signal-to-noise ratio of a narrowband signal can be significantly reduced when combined with detector output from other frequencies. Data folding and the computationally efficient analysis pipeline, {\tt PyStoch}, enable us to perform the radiometer map-making at every frequency bin. We perform the search at 3072 {\tt{HEALPix}} equal area pixels uniformly tiling the sky and in every frequency bin of width $1/32$~Hz in the range $20-1726$~Hz, except for bins that are likely to contain instrumental artefacts and hence are notched. We do not find any statistically significant evidence for the existence of narrowband gravitational-wave signals in the analyzed frequency bins. Therefore, we place $95\%$ confidence upper limits on the gravitational-wave strain for each pixel-frequency pair, the limits are in the range $(0.030 - 9.6) \times10^{-24}$. In addition, we outline a method to identify candidate pixel-frequency pairs that could be followed up by a more sensitive (and potentially computationally expensive) search, e.g., a matched-filtering-based analysis, to look for fainter nearly monochromatic coherent signals. The ASAF analysis is inherently independent of models describing any spectral or spatial distribution of power. We demonstrate that the ASAF results can be appropriately combined over frequencies and sky directions to successfully recover the broadband directional and isotropic results.

preprint2021arXiv

Overview of KAGRA: Calibration, detector characterization, physical environmental monitors, and the geophysics interferometer

KAGRA is a newly built gravitational wave observatory, a laser interferometer with a 3 km arm length, located in Kamioka, Gifu, Japan. In this series of articles, we present an overview of the baseline KAGRA, for which we finished installing the designed configuration in 2019. This article describes the method of calibration (CAL) used for reconstructing gravitational wave signals from the detector outputs, as well as the characterization of the detector (DET). We also review the physical environmental monitors (PEM) system and the geophysics interferometer (GIF). Both are used for characterizing and evaluating the data quality of the gravitational wave channel. They play important roles in utilizing the detector output for gravitational wave searches. These characterization investigations will be even more important in the near future, once gravitational wave detection has been achieved, and in using KAGRA in the gravitational wave astronomy era.

preprint2021arXiv

Search for Gravitational Waves Associated with Gamma-Ray Bursts Detected by Fermi and Swift During the LIGO-Virgo Run O3b

We search for gravitational-wave signals associated with gamma-ray bursts detected by the Fermi and Swift satellites during the second half of the third observing run of Advanced LIGO and Advanced Virgo (1 November 2019 15:00 UTC-27 March 2020 17:00 UTC).We conduct two independent searches: a generic gravitational-wave transients search to analyze 86 gamma-ray bursts and an analysis to target binary mergers with at least one neutron star as short gamma-ray burst progenitors for 17 events. We find no significant evidence for gravitational-wave signals associated with any of these gamma-ray bursts. A weighted binomial test of the combined results finds no evidence for sub-threshold gravitational wave signals associated with this GRB ensemble either. We use several source types and signal morphologies during the searches, resulting in lower bounds on the estimated distance to each gamma-ray burst. Finally, we constrain the population of low luminosity short gamma-ray bursts using results from the first to the third observing runs of Advanced LIGO and Advanced Virgo. The resulting population is in accordance with the local binary neutron star merger rate.

preprint2021arXiv

Search for subsolar-mass binaries in the first half of Advanced LIGO and Virgo's third observing run

We report on a search for compact binary coalescences where at least one binary component has a mass between 0.2 $M_\odot$ and 1.0 $M_\odot$ in Advanced LIGO and Advanced Virgo data collected between 1 April 2019 1500 UTC and 1 October 2019 1500 UTC. We extend previous analyses in two main ways: we include data from the Virgo detector and we allow for more unequal mass systems, with mass ratio $q \geq 0.1$. We do not report any gravitational-wave candidates. The most significant trigger has a false alarm rate of 0.14 $\mathrm{yr}^{-1}$. This implies an upper limit on the merger rate of subsolar binaries in the range $[220-24200] \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$, depending on the chirp mass of the binary. We use this upper limit to derive astrophysical constraints on two phenomenological models that could produce subsolar-mass compact objects. One is an isotropic distribution of equal-mass primordial black holes. Using this model, we find that the fraction of dark matter in primordial black holes is $f_\mathrm{PBH} \equiv Ω_\mathrm{PBH} / Ω_\mathrm{DM} \lesssim 6\%$. The other is a dissipative dark matter model, in which fermionic dark matter can collapse and form black holes. The upper limit on the fraction of dark matter black holes depends on the minimum mass of the black holes that can be formed: the most constraining result is obtained at $M_\mathrm{min}=1 M_\odot$, where $f_\mathrm{DBH} \equiv Ω_\mathrm{PBH} / Ω_\mathrm{DM} \lesssim 0.003\%$. These are the tightest limits on spinning subsolar-mass binaries to date.

preprint2021arXiv

The Scientific Observation Campaign of the Hayabusa-2 Capsule Re-entry

On 5th December 2020 at 17:28 UTC, the Japan Aerospace Exploration Agency's Hayabusa-2 sample return capsule came back to the Earth. It re-entered the atmosphere over South Australia, visible for 53 seconds as a fireball from near the Northern Territory border toward Woomera where it landed in the the Woomera military test range. A scientific observation campaign was planned to observe the optical, seismo-acoustic, radio and high energy particle phenomena associated with the entry of an interplanetary object. A multi-institutional collaboration between Australian and Japanese universities resulted in the deployment of 49 instruments, with a further 13 permanent observation sites. The campaign successfully recorded optical, seismo-acoustic and spectral data for this event which will allow an in depth analysis of the effects produced by interplanetary objects impacting the Earth's atmosphere. This will allow future comparison and insights to be made with natural meteoroid objects.

preprint2020arXiv

Application of the independent component analysis to the iKAGRA data

We apply the independent component analysis (ICA) to the real data from a gravitational wave detector for the first time. Specifically we use the iKAGRA data taken in April 2016, and calculate the correlations between the gravitational wave strain channel and 35 physical environmental channels. Using a couple of seismic channels which are found to be strongly correlated with the strain, we perform ICA. Injecting a sinusoidal continuous signal in the strain channel, we find that ICA recovers correct parameters with enhanced signal-to-noise ratio, which demonstrates usefulness of this method. Among the two implementations of ICA used here, we find the correlation method yields the optimal result for the case environmental noises act on the strain channel linearly.

preprint2020arXiv

Overview of KAGRA: Detector design and construction history

KAGRA is a newly built gravitational-wave telescope, a laser interferometer comprising arms with a length of 3\,km, located in Kamioka, Gifu, Japan. KAGRA was constructed under the ground and it is operated using cryogenic mirrors that help in reducing the seismic and thermal noise. Both technologies are expected to provide directions for the future of gravitational-wave telescopes. In 2019, KAGRA finished all installations with the designed configuration, which we call the baseline KAGRA. In this occasion, we present an overview of the baseline KAGRA from various viewpoints in a series of of articles. In this article, we introduce the design configurations of KAGRA with its historical background.

preprint2020arXiv

Search for $η'$ bound nuclei in the $^{12}{\rm C}(γ,p)$ reaction with simultaneous detection of decay products

We measured missing mass spectrum of the $^{12}{\rm C}(γ,p)$ reaction for the first time in coincidence with potential decay products from $η'$ bound nuclei. We tagged an ($η+p$) pair associated with the $η'N\toηN$ process in a nucleus. After applying kinematical selections to reduce backgrounds, no signal events were observed in the bound-state region. An upper limit of the signal cross section in the opening angle $\cosθ^{ηp}_{lab}<-0.9$ was obtained to be 2.2 nb/sr at the 90$\%$ confidence level. It is compared with theoretical cross sections, whose normalization ambiguity is suppressed by measuring a quasifree $η'$ production rate. Our results indicate a small branching fraction of the $η'N\toηN$ process and/or a shallow $η'$-nucleus potential.

preprint2019arXiv

A Measurement of the Degree Scale CMB B-mode Angular Power Spectrum with POLARBEAR

We present a measurement of the $B$-mode polarization power spectrum of the cosmic microwave background (CMB) using taken from July 2014 to December 2016 with the POLARBEAR experiment. The CMB power spectra are measured using observations at 150 GHz with an instantaneous array sensitivity of $\mathrm{NET}_\mathrm{array}=23\, μ\mathrm{K} \sqrt{\mathrm{s}}$ on a 670 square degree patch of sky centered at (RA, Dec)=($+0^\mathrm{h}12^\mathrm{m}0^\mathrm{s},-59^\circ18^\prime$). A continuously rotating half-wave plate is used to modulate polarization and to suppress low-frequency noise. We achieve $32\,μ\mathrm{K}$-$\mathrm{arcmin}$ effective polarization map noise with a knee in sensitivity of $\ell = 90$, where the inflationary gravitational wave signal is expected to peak. The measured $B$-mode power spectrum is consistent with a $Λ$CDM lensing and single dust component foreground model over a range of multipoles $50 \leq \ell \leq 600$. The data disfavor zero $C_\ell^{BB}$ at $2.2σ$ using this $\ell$ range of POLARBEAR data alone. We cross-correlate our data with Planck high frequency maps and find the low-$\ell$ $B$-mode power in the combined dataset to be consistent with thermal dust emission. We place an upper limit on the tensor-to-scalar ratio $r < 0.90$ at 95% confidence level after marginalizing over foregrounds.

preprint2019arXiv

An arm length stabilization system for KAGRA and future gravitational-wave detectors

Modern ground-based gravitational wave (GW) detectors require a complex interferometer configuration with multiple coupled optical cavities. Since achieving the resonances of the arm cavities is the most challenging among the lock acquisition processes, the scheme called arm length stabilization (ALS) had been employed for lock acquisition of the arm cavities. We designed a new type of the ALS, which is compatible with the interferometers having long arms like the next generation GW detectors. The features of the new ALS are that the control configuration is simpler than those of previous ones and that it is not necessary to lay optical fibers for the ALS along the kilometer-long arms of the detector. Along with simulations of its noise performance, an experimental test of the new ALS was performed utilizing a single arm cavity of KAGRA. This paper presents the first results of the test where we demonstrated that lock acquisition of the arm cavity was achieved using the new ALS and residual noise was measured to be $8.2\,\mathrm{Hz}$ in units of frequency, which is smaller than the linewidth of the arm cavity and thus low enough to lock the full interferometer of KAGRA in a repeatable and reliable manner.

preprint2017arXiv

Construction of KAGRA: an Underground Gravitational Wave Observatory

Major construction and initial-phase operation of a second-generation gravitational-wave detector KAGRA has been completed. The entire 3-km detector is installed underground in a mine in order to be isolated from background seismic vibrations on the surface. This allows us to achieve a good sensitivity at low frequencies and high stability of the detector. Bare-bones equipment for the interferometer operation has been installed and the first test run was accomplished in March and April of 2016 with a rather simple configuration. The initial configuration of KAGRA is named {\it iKAGRA}. In this paper, we summarize the construction of KAGRA, including the study of the advantages and challenges of building an underground detector and the operation of the iKAGRA interferometer together with the geophysics interferometer that has been constructed in the same tunnel.

preprint2017arXiv

Prospects for CTA observations of the young SNR RX J1713.7-3946

We perform simulations for future Cherenkov Telescope Array (CTA) observations of RX~J1713.7$-$3946, a young supernova remnant (SNR) and one of the brightest sources ever discovered in very-high-energy (VHE) gamma rays. Special attention is paid to explore possible spatial (anti-)correlations of gamma rays with emission at other wavelengths, in particular X-rays and CO/H{\sc i} emission. We present a series of simulated images of RX J1713.7$-$3946 for CTA based on a set of observationally motivated models for the gamma-ray emission. In these models, VHE gamma rays produced by high-energy electrons are assumed to trace the non-thermal X-ray emission observed by {\it XMM-Newton}, whereas those originating from relativistic protons delineate the local gas distributions. The local atomic and molecular gas distributions are deduced by the NANTEN team from CO and H{\sc i} observations. Our primary goal is to show how one can distinguish the emission mechanism(s) of the gamma rays (i.e., hadronic vs leptonic, or a mixture of the two) through information provided by their spatial distribution, spectra, and time variation. This work is the first attempt to quantitatively evaluate the capabilities of CTA to achieve various proposed scientific goals by observing this important cosmic particle accelerator.

preprint2017arXiv

The status of KAGRA underground cryogenic gravitational wave telescope

KAGRA is a 3-km interferometric gravitational wave telescope located in the Kamioka mine in Japan. It is the first km-class gravitational wave telescope constructed underground to reduce seismic noise, and the first km-class telescope to use cryogenic cooling of test masses to reduce thermal noise. The construction of the infrastructure to house the interferometer in the tunnel, and the initial phase operation of the interferometer with a simple 3-km Michelson configuration have been completed. The first cryogenic operation is expected in 2018, and the observing runs with a full interferometer are expected in 2020s. The basic interferometer configuration and the current status of KAGRA are described.

preprint2016arXiv

A significant hardening and rising shape detected in the MeV/GeV nuFnu spectrum from the recently-discovered very-high-energy blazar S4 0954+65 during the bright optical flare in 2015 February

We report on Fermi Large Area Telescope (LAT) and multi-wavelength results on the recently-discovered very-high-energy (VHE, $E>$ 100 GeV) blazar S4 0954+65 ($z=0.368$) during an exceptionally bright optical flare in 2015 February. During the time period (2015 February, 13/14, or MJD 57067) when the MAGIC telescope detected VHE $γ$-ray emission from the source, the Fermi-LAT data indicated a significant spectral hardening at GeV energies, with a power-law photon index of $1.8 \pm 0.1$---compared with the 3FGL value (averaged over four years of observation) of $2.34 \pm 0.04$. In contrast, Swift/XRT data showed a softening of the X-ray spectrum, with a photon index of $1.72 \pm 0.08$ (compared with $1.38 \pm 0.03$ averaged during the flare from MJD 57066 to 57077), possibly indicating a modest contribution of synchrotron photons by the highest-energy electrons superposed on the inverse Compton component. Fitting of the quasi-simultaneous ($<1$ day) broadband spectrum with a one-zone synchrotron plus inverse-Compton model revealed that GeV/TeV emission could be produced by inverse-Compton scattering of external photons from the dust torus. We emphasize that a flaring blazar showing high flux of $\gtrsim 1.0 \times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$ ($E>$ 100 MeV) and a hard spectral index of $Γ_{\rm GeV} < 2.0$ detected by Fermi-LAT on daily time scales is a promising target for TeV follow-up by ground-based Cherenkov telescopes to discover high-redshift blazars, investigate their temporal variability and spectral features in the VHE band, and also constrain the intensity of the extragalactic background light.

preprint2016arXiv

Galactic-Center Hyper-Shell Model for the North Polar Spurs

The bipolar-hyper shell (BHS) model for the North Polar Spurs (NPS-E, -W, and Loop I) and counter southern spurs (SPS-E and -W) is revisited based on numerical hydrodynamical simulations. Propagations of shock waves produced by energetic explosive events in the Galactic Center are examined. Distributions of soft X-ray brightness on the sky at 0.25, 0.7, and 1.5 keV in a +/-50 deg x +/-50 deg region around the Galactic Center are modeled by thermal emission from high-temperature plasma in the shock-compressed shell considering shadowing by the interstellar HI and H2 gases. The result is compared with the ROSAT wide field X-ray images in R2, 4 and 6 bands. The NPS and southern spurs are well reproduced by the simulation as shadowed dumbbell-shaped shock waves. We discuss the origin and energetics of the event in relation to the starburst and/or AGN activities in the Galactic Center. [ High resolution pdf is available at http://www.ioa.s.u-tokyo.ac.jp/~sofue/htdocs/2016bhs/ ]

preprint2016arXiv

No Evidence of Intrinsic Optical/Near-Infrared Linear Polarization for V404 Cygni During its Bright Outburst in 2015: Broadband Modeling and Constraint on Jet Parameters

We present simultaneous optical and near-infrared (NIR) polarimetric results for the black hole binary V404 Cygni spanning the duration of its 7-day long optically-brightest phase of its 2015 June outburst. The simultaneous R and Ks-band light curves showed almost the same temporal variation except for the isolated (~30 min duration) orphan Ks-band flare observed at MJD 57193.54. We did not find any significant temporal variation of polarization degree (PD) and position angle (PA) in both R and Ks bands throughout our observations, including the duration of the orphan NIR flare. We show that the observed PD and PA are predominantly interstellar in origin by comparing the V404 Cyg polarimetric results with those of the surrounding sources within the 7'x7' field-of-view. The low intrinsic PD (less than a few percent) implies that the optical and NIR emissions are dominated by either disk or optically-thick synchrotron emission, or both. We also present the broadband spectra of V404 Cyg during the orphan NIR flare and a relatively faint and steady state by including quasi-simultaneous Swift/XRT and INTEGRAL fluxes. By adopting a single-zone synchrotron plus inverse-Compton model as widely used in modeling of blazars, we constrained the parameters of a putative jet. Because the jet synchrotron component cannot exceed the Swift/XRT disk/corona flux, the cutoff Lorentz factor in the electron energy distribution is constrained to be <10^2, suggesting particle acceleration is less efficient in this microquasar jet outburst compared to AGN jets. We also suggest that the loading of the baryon component inside the jet is inevitable based on energetic arguments.

preprint2016arXiv

Observation of gamma ray bursts at ground level under the thunderclouds

We observed three $γ$-ray bursts related to thunderclouds in winter using the prototype of anti-neutrino detector PANDA made of 360-kg plastic scintillator deployed at Ohi Power Station at the coastal area of the Japan Sea. The maximum rate of the events which deposited the energy higher than $3\,$MeV was $(5.5 \pm 0.1) \times 10^2 {\rm /s}$. Monte Carlo simulation showed that electrons with approximately monochromatic energy falling downwards from altitudes of order $100\,$m roughly produced the observed total energy spectra of the bursts. It is supposed that secondary cosmic-ray electrons, which act as seed, were accelerated in electric field of thunderclouds and multiplied by relativistic runaway electron avalanche. We actually found that the $γ$-rays of the bursts entered into the detector from the direction close to the zenith. The direction stayed constant during the burst within the detector resolution. In addition, taking advantage of the delayed coincidence detection of the detector, we found neutron events in one of the bursts at the maximum rate of $\sim 14\pm5\,{\rm /s}$.

preprint2016arXiv

Transverse energy production and charged-particle multiplicity at midrapidity in various systems from $\sqrt{s_{NN}}=7.7$ to 200 GeV

Measurements of midrapidity charged particle multiplicity distributions, $dN_{\rm ch}/dη$, and midrapidity transverse-energy distributions, $dE_T/dη$, are presented for a variety of collision systems and energies. Included are distributions for Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$, 130, 62.4, 39, 27, 19.6, 14.5, and 7.7 GeV, Cu$+$Cu collisions at $\sqrt{s_{_{NN}}}=200$ and 62.4 GeV, Cu$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV, U$+$U collisions at $\sqrt{s_{_{NN}}}=193$ GeV, $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV, $^{3}$He$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV, and $p$$+$$p$ collisions at $\sqrt{s_{_{NN}}}=200$ GeV. Centrality-dependent distributions at midrapidity are presented in terms of the number of nucleon participants, $N_{\rm part}$, and the number of constituent quark participants, $N_{q{\rm p}}$. For all $A$$+$$A$ collisions down to $\sqrt{s_{_{NN}}}=7.7$ GeV, it is observed that the midrapidity data are better described by scaling with $N_{q{\rm p}}$ than scaling with $N_{\rm part}$. Also presented are estimates of the Bjorken energy density, $\varepsilon_{\rm BJ}$, and the ratio of $dE_T/dη$ to $dN_{\rm ch}/dη$, the latter of which is seen to be constant as a function of centrality for all systems.

preprint2015arXiv

Gamma-ray flaring activity from the gravitationally lensed blazar PKS 1830-211 observed by Fermi LAT

The Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope routinely detects the highly dust-absorbed, reddened, and MeV-peaked flat spectrum radio quasar PKS 1830-211 (z=2.507). Its apparent isotropic gamma-ray luminosity (E>100 MeV) averaged over $\sim$ 3 years of observations and peaking on 2010 October 14/15 at 2.9 X 10^{50} erg s^{-1}, makes it among the brightest high-redshift Fermi blazars. No published model with a single lens can account for all of the observed characteristics of this complex system. Based on radio observations, one expects time delayed variability to follow about 25 days after a primary flare, with flux about a factor 1.5 less. Two large gamma-ray flares of PKS 1830-211 have been detected by the LAT in the considered period and no substantial evidence for such a delayed activity was found. This allows us to place a lower limit of about 6 on the gamma rays flux ratio between the two lensed images. Swift XRT observations from a dedicated Target of Opportunity program indicate a hard spectrum and with no significant correlation of X-ray flux with the gamma-ray variability. The spectral energy distribution can be modeled with inverse Compton scattering of thermal photons from the dusty torus. The implications of the LAT data in terms of variability, the lack of evident delayed flare events, and different radio and gamma-ray flux ratios are discussed. Microlensing effects, absorption, size and location of the emitting regions, the complex mass distribution of the system, an energy-dependent inner structure of the source, and flux suppression by the lens galaxy for one image path may be considered as hypotheses for understanding our results.

preprint2015arXiv

Global structure of isothermal X-ray emission along the Fermi bubbles

In our previous works (Kataoka et al. 2013, Tahara et al. 2015), we found absorbed thermal X-ray plasma with kT ~ 0.3 keV observed ubiquitously near the edges of the Fermi bubbles and interpreted this emission as weakly shock-heated Galactic halo (GH) gas. Here we present a systematic and uniform analysis of archival Suzaku (29 pointings; 6 newly presented) and Swift (68 pointings; 49 newly presented) data within Galactic longitudes |l| < 20 deg and latitude 5 deg < |b| < 60 deg, covering the whole extent of the Fermi bubbles. We show that the plasma temperature is constant at kT = 0.30+-0.07 keV, while the emission measure (EM) varies by an order of magnitude, increasing toward the Galactic center (i.e., low |b|) with enhancements at the north polar spur (NPS), SE-claw and NW-clump features. Moreover, the EM distribution of kT ~ 0.30 keV plasma is highly asymmetric in the northern and southern bubbles. Although the association of the X-ray emission with the bubbles is not conclusive, we compare the observed EM properties with simple models assuming (i) a filled halo without bubbles, whose gas density follows a hydrostatic isothermal model (King profile) and (ii) a bubble-in-halo in which two identical bubbles expand into the halo forming thick shells of swept halo gas. We argue that the EM profile in the north (b > 0 deg) favors (ii), whereas that of the south (b < 0 deg) is rather close to (i), but weak excess signature is clearly detected also in the south like NPS (South Polar Spur; SPS). Such an asymmetry, if due to the bubbles, cannot be fully understood only by the inclination of bubbles' axis against the Galactic disk normal, thus suggesting asymmetric outflow due to different environmental/initial condition.

preprint2015arXiv

Modeling atmospheric emission for CMB ground-based observations

Atmosphere is one of the most important noise sources for ground-based cosmic microwave background (CMB) experiments. By increasing optical loading on the detectors, it amplifies their effective noise, while its fluctuations introduce spatial and temporal correlations between detected signals. We present a physically motivated 3d-model of the atmosphere total intensity emission in the millimeter and sub-millimeter wavelengths. We derive a new analytical estimate for the correlation between detectors time-ordered data as a function of the instrument and survey design, as well as several atmospheric parameters such as wind, relative humidity, temperature and turbulence characteristics. Using an original numerical computation, we examine the effect of each physical parameter on the correlations in the time series of a given experiment. We then use a parametric-likelihood approach to validate the modeling and estimate atmosphere parameters from the POLARBEAR-I project first season data set. We derive a new 1.0% upper limit on the linear polarization fraction of atmospheric emission. We also compare our results to previous studies and weather station measurements. The proposed model can be used for realistic simulations of future ground-based CMB observations.

preprint2015arXiv

Six Years of Fermi-LAT and Multi-wavelength Monitoring of the Broad-Line Radio Galaxy 3C 120: Jet Dissipation at Sub-parsec Scales from the Central Engine

We present multi-wavelength monitoring results for the broad-line radio galaxy 3C 120 in the MeV/GeV, sub-millimeter, and 43 GHz bands over six years. Over the past two years, Fermi-LAT sporadically detected 3C 120 with high significance and the 230 GHz data also suggest an enhanced activity of the source. After the MeV/GeV detection from 3C 120 in MJD 56240-56300, 43 GHz VLBA monitoring revealed a brightening of the radio core, followed by the ejection of a superluminal knot. Since we observed the gamma-ray and VLBA phenomena in temporal proximity to each other, it is naturally assumed that they are physically connected. This assumption was further supported by the subsequent observation that the 43 GHz core brightened again after a gamma-ray flare occurred around MJD 56560. We can then infer that the MeV/GeV emission took place inside an unresolved 43 GHz core of 3C 120 and that the jet dissipation occurred at sub-parsec distances from the central black hole, if we take the distance of the 43 GHz core from the central black hole as ~ 0.5 pc, as previously estimated from the time lag between X-ray dips and knot ejections (Marscher et al. 2002; Chatterjee et al. 2009). Based on our constraints on the relative locations of the emission regions and energetic arguments, we conclude that the gamma rays are more favorably produced via the synchrotron self-Compton process, rather than inverse Compton scattering of external photons coming from the broad line region or hot dusty torus. We also derived the electron distribution and magnetic field by modeling the simultaneous broadband spectrum.

preprint2015arXiv

Suzaku X-ray Observations of the Fermi Bubbles: Northernmost Cap and Southeast Claw Discovered with MAXI-SSC

We report on Suzaku observations of large-scale X-ray structures possibly related with the Fermi Bubbles obtained in 2013 with a total duration of ~ 80 ks. The observed regions were the: (i) northern cap (N-cap; l ~ 0 deg, 45 deg < b < 55 deg) seen in the Mid-band (1.7-4.0 keV) map recently provided by MAXI-SSC and (ii) southeast claw (SE-claw; l ~ 10 deg, -20 deg < b < -10 deg) seen in the ROSAT all-sky map and MAXI-SSC Low-band (0.7-1.7 keV) map. In each region, we detected diffuse X-ray emissions which are represented by a three component plasma model consisting of an unabsorbed thermal component (kT ~ 0.1 keV) from the Local Bubble, absorbed kT = 0.30+/-0.05 keV emission representing the Galactic Halo, and a power-law component due to the isotropic cosmic X-ray background radiation. The emission measure of the GH component in the SE-claw shows an excess by a factor of ~ 2.5 over the surrounding emission at 2 deg away. We also found a broad excess in the 1.7-4.0 keV count rates across the N-cap after compiling other archival data from Suzaku and Swift. The spectral stacking analysis of the N-cap data indicates the presence of another thermal component with kT = 0.70 (+0.22,-0.11) keV. The temperature of kT ~ 0.3 keV of the Galactic Halo is higher than the ubiquitous value of kT ~ 0.2 keV near the Fermi Bubbles, and can be even higher (~ 0.7 keV). We discuss our findings in the context of bubble-halo interaction.

preprint2015arXiv

The Origin of the Extragalactic Gamma-Ray Background and Implications for Dark-Matter Annihilation

The origin of the extragalactic $γ$-ray background (EGB) has been debated for some time. { The EGB comprises the $γ$-ray emission from resolved and unresolved extragalactic sources, such as blazars, star-forming galaxies and radio galaxies, as well as radiation from truly diffuse processes.} This letter focuses on the blazar source class, the most numerous detected population, and presents an updated luminosity function and spectral energy distribution model consistent with the blazar observations performed by the {\it Fermi} Large Area Telescope (LAT). We show that blazars account for 50$^{+12}_{-11}$\,\% of the EGB photons ($>$0.1\,GeV), and that {\it Fermi}-LAT has already resolved $\sim$70\,\% of this contribution. Blazars, and in particular low-luminosity hard-spectrum nearby sources like BL Lacs, are responsible for most of the EGB emission above 100\,GeV. We find that the extragalactic background light, which attenuates blazars' high-energy emission, is responsible for the high-energy cut-off observed in the EGB spectrum. Finally, we show that blazars, star-forming galaxies and radio galaxies can naturally account for the amplitude and spectral shape of the background in the 0.1--820\,GeV range, leaving only modest room for other contributions. This allows us to set competitive constraints on the dark-matter annihilation cross section.

preprint2015arXiv

The Physics of the B Factories

This work is on the Physics of the B Factories. Part A of this book contains a brief description of the SLAC and KEK B Factories as well as their detectors, BaBar and Belle, and data taking related issues. Part B discusses tools and methods used by the experiments in order to obtain results. The results themselves can be found in Part C. Please note that version 3 on the archive is the auxiliary version of the Physics of the B Factories book. This uses the notation alpha, beta, gamma for the angles of the Unitarity Triangle. The nominal version uses the notation phi_1, phi_2 and phi_3. Please cite this work as Eur. Phys. J. C74 (2014) 3026.

preprint2015arXiv

The POLARBEAR-2 and the Simons Array Experiment

We present an overview of the design and status of the \Pb-2 and the Simons Array experiments. \Pb-2 is a Cosmic Microwave Background polarimetry experiment which aims to characterize the arc-minute angular scale B-mode signal from weak gravitational lensing and search for the degree angular scale B-mode signal from inflationary gravitational waves. The receiver has a 365~mm diameter focal plane cooled to 270~milli-Kelvin. The focal plane is filled with 7,588 dichroic lenslet-antenna coupled polarization sensitive Transition Edge Sensor (TES) bolometric pixels that are sensitive to 95~GHz and 150~GHz bands simultaneously. The TES bolometers are read-out by SQUIDs with 40 channel frequency domain multiplexing. Refractive optical elements are made with high purity alumina to achieve high optical throughput. The receiver is designed to achieve noise equivalent temperature of 5.8~$μ$K$_{CMB}\sqrt{s}$ in each frequency band. \Pb-2 will deploy in 2016 in the Atacama desert in Chile. The Simons Array is a project to further increase sensitivity by deploying three \Pb-2 type receivers. The Simons Array will cover 95~GHz, 150~GHz and 220~GHz frequency bands for foreground control. The Simons Array will be able to constrain tensor-to-scalar ratio and sum of neutrino masses to $σ(r) = 6\times 10^{-3}$ at $r = 0.1$ and $\sum m_ν(σ=1)$ to 40 meV.

preprint2014arXiv

Development and characterization of the readout system for POLARBEAR-2

POLARBEAR-2 is a next-generation receiver for precision measurements of the polarization of the cosmic microwave background (Cosmic Microwave Background (CMB)). Scheduled to deploy in early 2015, it will observe alongside the existing POLARBEAR-1 receiver, on a new telescope in the Simons Array on Cerro Toco in the Atacama desert of Chile. For increased sensitivity, it will feature a larger area focal plane, with a total of 7,588 polarization sensitive antenna-coupled Transition Edge Sensor (TES) bolometers, with a design sensitivity of 4.1 uKrt(s). The focal plane will be cooled to 250 milliKelvin, and the bolometers will be read-out with 40x frequency domain multiplexing, with 36 optical bolometers on a single SQUID amplifier, along with 2 dark bolometers and 2 calibration resistors. To increase the multiplexing factor from 8x for POLARBEAR-1 to 40x for POLARBEAR-2 requires additional bandwidth for SQUID readout and well-defined frequency channel spacing. Extending to these higher frequencies requires new components and design for the LC filters which define channel spacing. The LC filters are cold resonant circuits with an inductor and capacitor in series with each bolometer, and stray inductance in the wiring and equivalent series resistance from the capacitors can affect bolometer operation. We present results from characterizing these new readout components. Integration of the readout system is being done first on a small scale, to ensure that the readout system does not affect bolometer sensitivity or stability, and to validate the overall system before expansion into the full receiver. We present the status of readout integration, and the initial results and status of components for the full array.

preprint2014arXiv

Evidence for Gravitational Lensing of the Cosmic Microwave Background Polarization from Cross-correlation with the Cosmic Infrared Background

We reconstruct the gravitational lensing convergence signal from Cosmic Microwave Background (CMB) polarization data taken by the POLARBEAR experiment and cross-correlate it with Cosmic Infrared Background (CIB) maps from the Herschel satellite. From the cross-spectra, we obtain evidence for gravitational lensing of the CMB polarization at a statistical significance of 4.0$σ$ and evidence for the presence of a lensing $B$-mode signal at a significance of 2.3$σ$. We demonstrate that our results are not biased by instrumental and astrophysical systematic errors by performing null-tests, checks with simulated and real data, and analytical calculations. This measurement of polarization lensing, made via the robust cross-correlation channel, not only reinforces POLARBEAR auto-correlation measurements, but also represents one of the early steps towards establishing CMB polarization lensing as a powerful new probe of cosmology and astrophysics.

preprint2014arXiv

Inferred cosmic-ray spectrum from ${\it Fermi}$-LAT $γ$-ray observations of the Earth's limb

Recent accurate measurements of cosmic-ray (CR) species by ATIC-2, CREAM, and PAMELA reveal an unexpected hardening in the proton and He spectra above a few hundred GeV, a gradual softening of the spectra just below a few hundred GeV, and a harder spectrum of He compared to that of protons. These newly-discovered features may offer a clue to the origin of high-energy CRs. We use the ${\it Fermi}$ Large Area Telescope observations of the $γ$-ray emission from the Earth's limb for an indirect measurement of the local spectrum of CR protons in the energy range $\sim 90~$GeV-$6~$TeV (derived from a photon energy range $15~$GeV-$1~$TeV). Our analysis shows that single power law and broken power law spectra fit the data equally well and yield a proton spectrum with index $2.68 \pm 0.04$ and $2.61 \pm 0.08$ above $\sim 200~$GeV, respectively.

preprint2014arXiv

Measurement of the Cosmic Microwave Background Polarization Lensing Power Spectrum with the POLARBEAR experiment

Gravitational lensing due to the large-scale distribution of matter in the cosmos distorts the primordial Cosmic Microwave Background (CMB) and thereby induces new, small-scale $B$-mode polarization. This signal carries detailed information about the distribution of all the gravitating matter between the observer and CMB last scattering surface. We report the first direct evidence for polarization lensing based on purely CMB information, from using the four-point correlations of even- and odd-parity $E$- and $B$-mode polarization mapped over $\sim30$ square degrees of the sky measured by the POLARBEAR experiment. These data were analyzed using a blind analysis framework and checked for spurious systematic contamination using null tests and simulations. Evidence for the signal of polarization lensing and lensing $B$-modes is found at 4.2$σ$ (stat.+sys.) significance. The amplitude of matter fluctuations is measured with a precision of $27\%$, and is found to be consistent with the Lambda Cold Dark Matter ($Λ$CDM) cosmological model. This measurement demonstrates a new technique, capable of mapping all gravitating matter in the Universe, sensitive to the sum of neutrino masses, and essential for cleaning the lensing $B$-mode signal in searches for primordial gravitational waves.

preprint2014arXiv

Reactor antineutrino monitoring with a plastic scintillator array as a new safeguards method

We developed a segmented reactor-antineutrino detector made of plastic scintillators for application as a tool in nuclear safeguards inspection and performed mostly unmanned field operations at a commercial power plant reactor. At a position outside the reactor building, we measured the difference in reactor antineutrino flux above the ground when the reactor was active and inactive.

preprint2014arXiv

Suzaku observations of the diffuse x-ray emission across the fermi bubbles' edges

We present Suzaku X-ray observations along two edge regions of the Fermi Bubbles, with eight ~20 ksec pointings across the northern part of the North Polar Spur (NPS) surrounding the north bubble and six across the southernmost edge of the south bubble. After removing compact X-ray features, diffuse X-ray emission is clearly detected and is well reproduced by a three-component spectral model consisting of unabsorbed thermal emission (temperature kT ~0.1 keV from the Local Bubble (LB), absorbed kT ~0.3 keV thermal emission related to the NPS and/or Galactic Halo (GH), and a power-law component at a level consistent with the cosmic X-ray background. The emission measure (EM) of the 0.3 keV plasma decreases by ~50% toward the inner regions of the north-east bubble, with no accompanying temperature change. However, such a jump in the EM is not clearly seen in the south bubble data. While it is unclear if the NPS originates from a nearby supernova remnant or is related to previous activity within/around the Galactic Center, our Suzaku observations provide evidence suggestive of the latter scenario. In the latter framework, the presence of a large amount of neutral matter absorbing the X-ray emission as well as the existence of the kT ~ 0.3 keV gas can be naturally interpreted as a weak shock driven by the bubbles' expansion in the surrounding medium, with velocity v_exp ~300 km/s (corresponding to shock Mach number M ~1.5), compressing the GH gas to form the NPS feature. We also derived an upper limit for any non-thermal X-ray emission component associated with the bubbles and demonstrate, that in agreement with the findings above, the non-thermal pressure and energy estimated from a one-zone leptonic model of its broad-band spectrum, are in rough equilibrium with that of the surrounding thermal plasma.

preprint2013arXiv

Constraints on the Galactic Population of TEV Pulsar Wind Nebulae Using Fermi Large Area Telescope Observations

Pulsar wind nebulae (PWNe) have been established as the most populous class of TeV gamma-ray emitters. Since launch, the Fermi Large Area Telescope (LAT)identified five high-energy (100MeV <E< 100 GeV) gamma-ray sources as PWNe, and detected a large number of PWNe candidates, all powered by young and energetic pulsars. The wealth of multi-wavelength data available and the new results provided by Fermi-LAT give us an opportunity to find new PWNe and to explore the radiative processes taking place in known ones. The TeV gamma-ray unidentifiedsources (UNIDs) are the best candidates for finding new PWNe. Using 45 months of Fermi-LAT data for energies above 10 GeV, an analysis was performed near the position of 58TeV PWNe and UNIDs within 5deg of the Galactic Plane to establish new constraints on PWNe properties and find new clues on the nature of UNIDs. Of the 58 sources, 30 were detected, and this work provides their gamma-rayfluxes for energies above 10 GeV. The spectral energy distributions (SED) andupper limits, in the multi-wavelength context, also provide new information on the source nature and can help distinguish between emission scenarios, i.e. between classification as a pulsar candidate or as a PWN candidate. Six new GeV PWNe candidates are described in detail and compared with existing models. A population study of GeV PWNe candidates as a function of the pulsar/PWN system characteristics is presented.

preprint2013arXiv

Experimental search for solar hidden photons in the eV energy range using kinetic mixing with photons

We have searched for solar hidden photons in the eV energy range using a dedicated hidden photon detector. The detector consisted of a parabolic mirror with a diameter of 500mm and a focal length of 1007mm installed in a vacuum chamber, and a photomultiplier tube at its focal point. The detector was attached to the Tokyo axion helioscope, Sumico which has a mechanism to track the sun. From the result of the measurement, we found no evidence for the existence of hidden photons and set a limit on the photon-hidden photon mixing parameter χdepending on the hidden photon mass m_{γ'}.

preprint2013arXiv

Fermi Large Area Telescope Detection of Two Very-High-Energy (E>100 GeV) Gamma-ray Photons from the z = 1.1 Blazar PKS 0426-380

We report the Fermi Large Area Telescope (LAT) detection of two very-high-energy (VHE, E>100 GeV) gamma-ray photons from the directional vicinity of the distant (redshift, z = 1.1) blazar PKS 0426-380. The null hypothesis that both the 134 and 122 GeV photons originate from unrelated sources can be rejected at the 5.5 sigma confidence level. We therefore claim that at least one of the two VHE photons is securely associated with the blazar, making PKS 0426-380 the most distant VHE emitter known to date. The results are in agreement with the most recent Fermi-LAT constraints on the Extragalactic Background Light (EBL) intensity, which imply a $z \simeq 1$ horizon for $\simeq$ 100 GeV photons. The LAT detection of the two VHE gamma-rays coincided roughly with flaring states of the source, although we did not find an exact correspondence between the VHE photon arrival times and the flux maxima at lower gamma-ray energies. Modeling the gamma-ray continuum of PKS 0426-380 with daily bins revealed a significant spectral hardening around the time of detection of the first VHE event (LAT photon index Γ $\simeq$ 1.4) but on the other hand no pronounced spectral changes near the detection time of the second one. This combination implies a rather complex variability pattern of the source in gamma rays during the flaring epochs. An additional flat component is possibly present above several tens of GeV in the EBL-corrected Fermi-LAT spectrum accumulated over the ~8-month high state.

preprint2013arXiv

First Results from NuSTAR Observations of Mkn 421

Mkn 421 is a nearby active galactic nucleus dominated at all wavelengths by a very broad non-thermal continuum thought to arise from a relativistic jet seen at a small angle to the line of sight. Its spectral energy distribution peaks in the X-ray and TeV gamma-ray bands, where the energy output is dominated by cooling of high-energy electrons in the jet. In order to study the electron distribution and its evolution, we carried out a dedicated multi-wavelength campaign, including extensive observations by the recently launched highly sensitive hard X-ray telescope NuSTAR, between December 2012 and May 2013. Here we present some initial results based on NuSTAR data from January through March 2013, as well as calibration observations conducted in 2012. Although the observations cover some of the faintest hard X-ray flux states ever observed for Mkn 421, the sensitivity is high enough to resolve intra-day spectral variability. We find that in this low state the dominant flux variations are smooth on timescales of hours, with typical intra-hour variations of less than 5%. We do not find evidence for either a cutoff in the hard X-ray spectrum, or a rise towards a high-energy component, but rather that at low flux the spectrum assumes a power law shape with a photon index of approximately 3. The spectrum is found to harden with increasing brightness.

preprint2013arXiv

Mission design of LiteBIRD

LiteBIRD is a next-generation satellite mission to measure the polarization of the cosmic microwave background (CMB) radiation. On large angular scales the B-mode polarization of the CMB carries the imprint of primordial gravitational waves, and its precise measurement would provide a powerful probe of the epoch of inflation. The goal of LiteBIRD is to achieve a measurement of the characterizing tensor to scalar ratio $r$ to an uncertainty of $δr=0.001$. In order to achieve this goal we will employ a kilo-pixel superconducting detector array on a cryogenically cooled sub-Kelvin focal plane with an optical system at a temperature of 4~K. We are currently considering two detector array options; transition edge sensor (TES) bolometers and microwave kinetic inductance detectors (MKID). In this paper we give an overview of LiteBIRD and describe a TES-based polarimeter designed to achieve the target sensitivity of 2~$μ$K$\cdot$arcmin over the frequency range 50 to 320~GHz.

preprint2013arXiv

X-ray and Radio Follow-up Observations of High-Redshift Blazar Candidates in the Fermi-LAT Unassociated Source Population

We report on the results of X-ray and radio follow-up observations of two GeV gamma-ray sources 2FGL J0923.5+1508 and 2FGL J1502.1+5548, selected as candidates for high-redshift blazars from unassociated sources in the {\it Fermi} Large Area Telescope Second Source Catalog. We utilize the Suzaku satellite and the VLBI Exploration of Radio Astrometry (VERA) telescopes for X-ray and radio observations, respectively. For 2FGL J0923.5+1508, a possible radio counterpart NVSS J092357+150518 is found at 1.4 GHz from an existing catalog, but we do not detect any X-ray emission from it and derive a flux upper limit $F_{\rm 2-8 keV} <$ 1.37 $\times$ 10$^{-14}$ erg cm$^{-2}$ s$^{-1}$. Radio observations at 6.7 GHz also result in an upper limit of $S_{\rm 6.7 GHz}$ $<$ 19 mJy, implying a steep radio spectrum that is not expected for a blazar. On the other hand, we detect X-rays from NVSS J150229+555204, the potential 1.4 GHz radio counterpart of 2FGL J1502.1+5548. The X-ray spectrum can be fitted with an absorbed power-law model with a photon index $γ$=1.8$^{+0.3}_{-0.2}$ and the unabsorbed flux is $F_{\rm 2-8 keV}$=4.3$^{+1.1}_{-1.0}$ $\times$ 10$^{-14}$ erg cm$^{-2}$ s$^{-1}$. Moreover, we detect unresolved radio emission at 6.7 GHz with flux $S_{\rm 6.7 GHz}$=30.1 mJy, indicating a compact, flat-spectrum radio source. If NVSS J150229+555204 is indeed associated with 2FGL J1502.1+5548, we find that its multiwavelength spectrum is consistent with a blazar at redshift $z \sim 3-4$.

preprint2012arXiv

Medium modification of jet fragmentation in Au+Au collisions at sqrt(s_NN)=200 GeV measured in direct photon-hadron correlations

The jet fragmentation function is measured with direct photon-hadron correlations in p+p and Au+Au collisions at sqrt(s_NN)=200 GeV. The p_T of the photon is an excellent approximation to the initial p_T of the jet and the ratio z_T=p_T^h/p_T^γis used as a proxy for the jet fragmentation function. A statistical subtraction is used to extract the direct photon-hadron yields in Au+Au collisions while a photon isolation cut is applied in p+p. I_ AA, the ratio of jet fragment yield in Au+Au to that in p+p, indicates modification of the jet fragmentation function. Suppression, most likely due to energy loss in the medium, is seen at high z_T. The fragment yield at low z_T is enhanced at large angles. Such a trend is expected from redistribution of the lost energy into increased production of low-momentum particles.

preprint2012arXiv

Orbital fluctuation theory in iron-based superconductors: s-wave superconductivity, structure transition, and impurity-induced nematic order

The main features in iron-based superconductors would be (i) the orthorhombic transition accompanied by remarkable softening of shear modulus, (ii) high-Tc superconductivity close to the orthorhombic phase, and (iii) nematic transition in the tetragonal phase. In this paper, we present a unified explanation for them, based on the orbital fluctuation theory, considering both the e-ph and the Coulomb interaction. It is found that a small e-ph coupling constant ($λ~ 0.2$) is enough to produce large orbital (=charge quadrupole $O_{xz/yz}$) fluctuations, which causes the s-wave superconductivity without sign reversal ($s_{++}$-wave state). The derived orbital fluctuations also cause the instability toward the structure transition due to the bound state formation of two orbitons with opposite momenta, which is called the "two-orbiton process". Moreover, impurity-induced non-local orbital order with $C_2$-symmetry is obtained when the orbital fluctuations are strong. This "impurity-induced nematic state" explains the in-plane anisotropy of resistivity in detwinned samples. We stress that (i)-(iii) are reproducible only when orbital fluctuations with respect to $O_{xz}$ and $O_{yz}$ charge quadrupoles are the most divergent. This fact ensures the reliability of the present model Hamiltonian and calculation.

preprint2012arXiv

The Tokyo Axion Helioscope

The Tokyo Axion Helioscope experiment aims to detect axions which are produced in the solar core. The helioscope uses a strong magnetic field in order to convert axions into X-ray photons and has a mounting to follow the sun very accurately. The photons are detected by an X-ray detector which is made of 16 PIN-photodiodes. In addition, a gas container and a gas regulation system is adopted for recovering the coherence between axions and photons in the conversion region giving sensitivity to axions with masses up to 2 eV. In this paper, we report on the technical detail of the Tokyo Axion Helioscope.

preprint2011arXiv

Production of omega mesons in p+p, d+Au, Cu+Cu, and Au+Au collisions at sqrt(s_NN)=200 GeV

The PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) has measured omega meson production via leptonic and hadronic decay channels in p+p, d+Au, Cu+Cu, and Au+Au collisions at sqrt(s_NN) = 200 GeV. The invariant transverse momentum spectra measured in different decay modes give consistent results. Measurements in the hadronic decay channel in Cu+Cu and Au+Au collisions show that omega production has a suppression pattern at high transverse momentum, similar to that of pi^0 and eta in central collisions, but no suppression is observed in peripheral collisions. The nuclear modification factors, R_AA, are consistent in Cu+Cu and Au+Au collisions at similar numbers of participant nucleons.

preprint2010arXiv

Event Structure and Double Helicity Asymmetry in Jet Production from Polarized p+p Collisions at sqrt(s) = 200 GeV

We report on event structure and double helicity asymmetry ($A_LL$) of jet production in longitudinally polarized p+p collisions at $\sqrt{s}$=200 GeV. Photons and charged particles were measured at midrapidity $|η| < 0.35$ with the requirement of a high-momentum ($>2$ GeV/$c$) photon in each event. Measured event structure is compared with {\sc pythia} and {\sc geant} simulations. The shape of jets and the underlying event were well reproduced at this collision energy. For the measurement of jet $A_{LL}$, photons and charged particles were clustered with a seed-cone algorithm to obtain the cluster $p_T$ sum ($p_T^{\rm reco}$). The effect of detector response and the underlying events on $p_T^{\rm reco}$ was evaluated with the simulation. The production rate of reconstructed jets is satisfactorily reproduced with the NLO pQCD jet production cross section. For $4 < p_T^{\rm reco} < 12$ GeV/$c$ with an average beam polarization of $< P > = 49%$ we measured $A_{LL} = -0.0014 \pm 0.0037^{\rm stat}$ at the lowest $p_T^{\rm reco}$ bin (4-5 GeV/$c$) and $-0.0181 \pm 0.0282^{\rm stat}$ at the highest $p_T^{\rm reco}$ bin (10-12 GeV/$c$) with a beam polarization scale error of 9.4% and a $\pT$ scale error of 10%. Jets in the measured $p_T^{\rm reco}$ range arise primarily from hard-scattered gluons with momentum fraction $0.02 < x < 0.3$ according to {\sc pythia}. The measured $A_{LL}$ is compared with predictions that assume various $ΔG(x)$ distributions based on the GRSV parameterization. The present result imposes the limit $-1.1 < \int_{0.02}^{0.3}dx ΔG(x, μ^2 = 1 {\rm GeV}^2) < 0.4$ at 95% confidence level or $\int_{0.02}^{0.3}dx ΔG(x, μ^2 = 1 {\rm GeV}^2) < 0.5$ at 99% confidence level.

preprint2010arXiv

Heavy Quark Production in p+p and Energy Loss and Flow of Heavy Quarks in Au+Au Collisions at sqrt(s_NN)=200 GeV

Transverse momentum (p^e_T) spectra of electrons from semileptonic weak decays of heavy flavor mesons in the range of 0.3 < p^e_T < 9.0 GeV/c have been measured at mid-rapidity (|eta| < 0.35) by the PHENIX experiment at the Relativistic Heavy Ion Collider in p+p and Au+Au collisions at sqrt(s_NN)=200 GeV. The nuclear modification factor R_AA with respect to p+p collisions indicates substantial energy loss of heavy quarks in the produced medium. In addition, the azimuthal anisotropy parameter v_2 has been measured for 0.3 < p^e_T < 5.0 GeV/c in Au+Au collisions. Comparisons of R_AA and v_2 are made to various model calculations.

preprint2010arXiv

Measurement of neutral mesons in p+p collisions at sqrt(s) = 200 GeV and scaling properties of hadron production

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured the invariant differential cross section for production of K^0_S , ω, ηprime, and ϕmesons in p + p collisions at = 200 GeV. Measurements ωand ϕproduction in different decay channels give consistent results. New results for the ϕare in agreement with previously published data and extend the measured pT coverage. The spectral shapes of all hadron transverse momentum distributions measured by PHENIX are well described by a Tsallis distribution functional form with only two parameters, n and T, determining the high-pT and characterizing the low-pT regions of the spectra, respectively. The values of these parameters are very similar for all analyzed meson spectra, but with a lower parameter T extracted for protons. The integrated invariant cross sections calculated from the fitted distributions are found to be consistent with existing measurements and with statistical model predictions.

preprint2010arXiv

Nuclear modification factors of phi mesons in d+Au, Cu+Cu and Au+Au collisions at sqrt(S_NN)=200 GeV

The PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) has performed systematic measurements of phi meson production in the K+K- decay channel at midrapidity in p+p, d+Au, Cu+Cu and Au+Au collisions at sqrt(S_NN)=200 GeV. Results are presented on the phi invariant yield and the nuclear modification factor R_AA for Au+Au and Cu+Cu, and R_dA for d+Au collisions, studied as a function of transverse momentum (1<p_T<7 GeV/c) and centrality. In central and mid-central Au+Au collisions, the R_AA of phi exhibits a suppression relative to expectations from binary scaled p+p results. The amount of suppression is smaller than that of the neutral pion and the eta meson in the intermediate p_T range (2--5 GeV/c); whereas at higher p_T the phi, pi^0, and eta show similar suppression. The baryon (protons and anti-protons) excess observed in central Au+Au collisions at intermediate p_T is not observed for the phi meson despite the similar mass of the proton and the phi. This suggests that the excess is linked to the number of constituent quarks rather than the hadron mass. The difference gradually disappears with decreasing centrality and for peripheral collisions the R_AA values for both particles are consistent with binary scaling. Cu+Cu collisions show the same yield and suppression as Au+Au collisions for the same number of N_part. The R_dA of phi shows no evidence for cold nuclear effects within uncertainties.

preprint2010arXiv

Phonon dynamics in SrFe2As2 and SrFeAsF parent pnictide compounds: How different/similar when compared to Ca- and Ba-compounds?

We report detailed measurements of temperature dependence and ab initio lattice dynamical calculations of the phonon density-of-states of two Sr based iron pnictide parent compounds; SrFe2As2 and SrFeAsF. Although both compounds undergo a magnetic as well as a tetragonal-to-orthorhombic structural phase transitions, no drastic change of their phonon spectra is found in the temperature range from 2 to 300 K.. In both the systems the Fe-As stretching modes at 32 meV react to the temperature change and shift towards higher energies with decreasing temperature. Furthermore, we found that in the antiferromagnetic phase of SrFe2As2 at 180 K a peak at 25 meV, due to Fe vibrations, moves slightly to higher energies and become narrower on cooling from 300 K. However we did not find any evidence of softening of low energy phonon modes in both SrFe2As2 and SrFeAsF on decreasing the temperature from 300 K in their paramagnetic phase. This is contrary to our earlier observation of phonon softening in CaFe2As2, Ca0.6Na0.4Fe2As2 and Sr0.6K0.4Fe2As2. For the sake of comparison, the previously studied Ca- and Ba-cases are also presented and the role of the magnetic ordering on the calculated phonon spectra in these two families of pnictides is reported.

preprint2010arXiv

Tokyo axion helioscope experiment

A search for solar axions has been performed using an axion helioscope which is equipped with a 2.3m x 4T superconducting magnet, a gas container to hold dispersion-matching gas, PIN-photodiode X-ray detectors, and a telescope mount mechanism to track the sun. In the past measurements, axion mass up to 0.27eV have been scanned. It has been upgraded to handle dispersion-matching gas (He-4) of higher density to explore higher mass region. From December 2007 through April 2008, the axion mass region 0.84 < m < 1.00eV was newly explored, where the axions in the "photon-coupling vs. mass" parameter region of the preferred axion models were newly searched. From the absence of any evidence, a limit on axion-photon coupling constant was set to be g < 5.6-13.4e-10/GeV at 95% confidence level in the above mass region.

preprint2009arXiv

Neutron diffraction study on phase transition and thermal expansion of SrFeAsF

The magnetic ordering and crystal structure of iron pnictide SrFeAsF was investigated by using neutron powder diffraction method. With decreasing temperature, the tetragonal to orthorhombic phase transition is found at 180 K, while the paramagnetic to antiferromagnetic phase transition set in at 133 K. Similar to the parent compound of other iron pnictide system, the striped Fe magnetism is confirmed in antiferromagnetic phase and the Fe moment of 0.58(6) uB aligned along long a axis. The thermal expansion of orthorhombic phase of SrFeAsF is also investigated. Based on the Grueneisen approximation and Debye approximation for internal energy, the volume of SrFeAsF can be well fitted with Debye temperature of 347(5) K. The experimental atomic displacement parameters for different crystallographic sites in SrFeAsF are analyzed with Debye model. The results suggested that the expansion of FeAs layers plays an important role in determining the thermal expansion coefficient.

preprint2009arXiv

Prospects of Search for Solar Axions with Mass over 1 eV and Hidden Sector Photons

We present prospects of two experiments using the Tokyo Axion Helioscope. One is a search for solar axions. In the past measurements, axion mass from 0 to 0.27 eV and from 0.84 to 1.00 eV have been scanned and no positive evidence was seen. We are now actively preparing a new phase of the experiment aiming at axion mass over 1 eV. The other is a search for hidden sector photons from the Sun. We have been designing and testing some additional equipments, which have to be installed on the helioscope to search for hidden photons with mass of over $10^{-3}$ eV.

preprint2008arXiv

Search for solar axions with mass around 1 eV using coherent conversion of axions into photons

A search for solar axions has been performed using an axion helioscope which is equipped with a 2.3m-long 4T superconducting magnet, a gas container to hold dispersion-matching gas, PIN-photodiode X-ray detectors, and a telescope mount mechanism to track the sun. A mass region around m_a = 1eV was newly explored. From the absence of any evidence, analysis sets a limit on axion-photon coupling constant to be g < 5.6-13.4x10^{-10} GeV^{-1} for the axion mass of 0.84<m_a<1.00eV at 95% confidence level. It is the first result to search for the axion in the g-m_a parameter region of the preferred axion models with a magnetic helioscope.

preprint2008arXiv

Status report of the Tokyo axion helioscope experiment

We have searched for solar axions with a detector which consists of a 4T x 2.3m superconducting magnet, PIN-photodiode X-ray detectors, and an altazimuth mount to track the sun. The conversion region is filled with cold helium gas which modifies the axion mass at which coherent conversion occurs. In the past measurements, axion mass from 0 to 0.27eV have been scanned. Since no positive evidence was seen, an upper limit to the axion-photon coupling constant was set to be g < 6-10E-10/GeV (95%CL) depending on the axion masses. We are now actively preparing for a new stage of the experiment aiming at one to a few eV solar axions. In this mass region, our detector might be able to check parameter regions which are preferable to the axion models.

preprint2004arXiv

Dark Matter Search with Direction Sensitive Scintillators

We have carried out the dark matter search with a 116g direction-sensitive stilbene crystal in Kamioka Observatory. With the crystal fixed to the earth, we searched the modulation of thelight output. No modulation signal was found due to the small size of the detector crystal and the higher background rate yet to be eliminated. However, it demonstrated the effectiveness of the method of direction sensitive search for the dark matter with an implementation of the anisotropic organic scintillation crystal.

preprint2002arXiv

Directional scintillation detector for the detection of the wind of WIMPs

The quenching factor of the proton recoil in the stilbene scintillator was measured with a 252Cf neutron source and was found to be 0.1 - 0.17 in the recoil energy range between 300 keV and 3 MeV. It was found to depend on the direction of the recoil proton. The directional anisotropy of the quenching factor could be used to detect the wind of the WIMPs caused by the motion of the earth around the galactic center.

preprint2000arXiv

Recent results from the Tokyo axion helioscope experiment

We have searched for axions which could have been produced in the solar core using an axion helioscope with a 2.3m-long 4T superconducting magnet. Axion mass region up to m_a=0.26eV was newly explored by introducing dispersion-matching gas. Preliminary analysis sets a limit on axion-photon coupling constant to be less than (6.4--9.6)e-10 GeV^{-1} (95%CL) for this mass region from the absence of the axion signal. This is more stringent than the limit inferred from the solar age consideration and also more stringent than the recent helioseismological bound.