Source author record

T. Yamamoto

T. Yamamoto 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

118works
17topics
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

118 published item(s)

preprint2026arXiv

Measurement of muon neutrino induced charged current interactions without charged pions in the final state using a new T2K off-axis near detector WAGASCI-BabyMIND

We report a flux-integrated cross section measurement of muon neutrino interactions on water and hydrocarbon via charged current reactions without charged pions in the final state with the WAGASCI-BabyMIND detector which was installed in the T2K near detector hall in 2018. The detector is located 1.5$^\circ$ off-axis and is exposed to a more energetic neutrino flux than ND280, another T2K near detector, which is located at a different off-axis position. The total flux-integrated cross section is measured to be $1.26 \pm 0.18\,(stat.+syst.) \times 10^{-39} $ $\mathrm{cm^{2}/nucleon}$ on CH and $1.44 \pm 0.21\,(stat.+syst.) \times 10^{-39} $ $\mathrm{cm^{2}/nucleon}$ on H$_{2}$O. These results are compared to model predictions provided by the NEUT v5.3.2 and GENIE v2.8.0 MC generators and the measurements are compatible with these models. Differential cross sections in muon momentum and cosine of the muon scattering angle are also reported. This is the first such measurement reported with the WAGASCI-BabyMIND detector and utilizes the 2020 and 2021 datasets.

preprint2025arXiv

Results from the T2K experiment on neutrino mixing including a new far detector $μ$-like sample

We have made improved measurements of three-flavor neutrino mixing with 19.7(16.3)$\times 10^{20}$ protons on target in (anti-)neutrino-enhanced beam modes. A new sample of muon-neutrino events with tagged pions has been added at the far detector, as well as new proton and photon-tagged samples at the near detector. Significant improvements have been made to the flux and neutrino interaction modeling. T2K data continue to prefer the normal mass ordering and upper octant of $\sin^2θ_{23}$ with a near-maximal value of the charge-parity violating phase with best-fit values in the normal ordering of $δ_{\scriptscriptstyle\mathrm{CP}}=-2.18\substack{+1.22 \\ -0.47}$, $\sin^2θ_{23}=0.559\substack{+0.018 \\ -0.078}$ and $Δm^2_{32}=(+2.506\substack{+0.039 \\ -0.052})\times 10^{-3}$ eV$^{2}$.

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

Angular distribution of $γ$ rays from the p-wave resonance of $^{118}$Sn

The neutron energy-dependent angular distribution of $γ$ rays from $^{117}{\rm Sn}(n,γ)$ reaction was measured with germanium detectors and a pulsed neutron beam. The angular distribution was clearly observed in $γ$-ray emissions with an energy of 9327 keV which corresponds to the transition from a neutron resonance of $^{117}{\rm Sn}+n$ to the ground state of $^{118}{\rm Sn}$. The angular distribution causes an angular-dependent asymmetric resonance shape. An asymmetry $A_{\rm LH}$ was defined as $(N_{\rm L}-N_{\rm H})/(N_{\rm L}+N_{\rm H})$, where $N_{\rm L}$ and $N_{\rm H}$ are integrated values for lower- and higher-energy regions of a neutron resonance, respectively. We found that the $A_{\rm LH}$ has the angular dependence of $(A \cos θ_γ+B)$, where $θ_γ$ is the $γ$-ray emission angle with respect to the incident neutron momentum, with $A=0.394 \pm 0.073$ and $B = 0.118 \pm 0.029$ in the 1.33 eV p-wave resonance.

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

Investigating the blazar TXS 0506+056 through sharp multi-wavelength eyes during 2017-2019

The blazar TXS 0506+056 got into the spotlight of the astrophysical community in September 2017, when a high-energy neutrino detected by IceCube (IceCube-170922A) was associated at the 3 $σ$ level to a $γ$-ray flare from this source. This multi-messenger photon-neutrino association remains, as per today, the most significant one ever observed. TXS 0506+056 was a poorly studied object before the IceCube-170922A event. To better characterize its broad-band emission, we organized a multi-wavelength campaign lasting 16 months (November 2017 to February 2019), covering the radio-band (Metsähovi, OVRO), the optical/UV (ASAS-SN, KVA, REM, Swift/UVOT), the X-rays (Swift/XRT, NuSTAR), the high-energy $γ$ rays (Fermi/LAT) and the very-high-energy (VHE) $γ$ rays (MAGIC). In $γ$ rays, the behaviour of the source was significantly different from the 2017 one: MAGIC observations show the presence of flaring activity during December 2018, while the source only shows an excess at the 4$σ$ level during the rest of the campaign (74 hours of accumulated exposure); Fermi/LAT observations show several short (days-to-week timescale) flares, different from the long-term brightening of 2017. No significant flares are detected at lower energies. The radio light curve shows an increasing flux trend, not seen in other wavelengths. We model the multi-wavelength spectral energy distributions in a lepto-hadronic scenario, in which the hadronic emission emerges as Bethe-Heitler and pion-decay cascade in the X-rays and VHE $γ$ rays. According to the model presented here, the December 2018 $γ$-ray flare was connected to a neutrino emission that was too brief and not bright enough to be detected by current neutrino instruments.

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.

preprint2022arXiv

Spectroscopic observation of the crossover from a classical Duffing oscillator to a Kerr parametric oscillator

We study microwave response of a Josephson parametric oscillator consisting of a superconducting transmission-line resonator with an embedded dc-SQUID. The dc-SQUID allows to control the magnitude of a Kerr nonlinearity over the ranges where it is smaller or larger than the photon loss rate. Spectroscopy measurements reveal the change of the microwave response from a classical Duffing oscillator to a Kerr parametric oscillator in a single device. In the single-photon Kerr regime, we observe parametric oscillations with a well-defined phase of either $0$ or $π$, whose probability can be controlled by an externally injected signal.

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

Isotopic effect of proton conductivity in barium-zirconates for various hydrogen-containing atmospheres

Among various perovskite proton conducting oxides, Y-doped BaZrO3 perovskite is a promising material for electrochemical hydrogen devices due to the good chemical stability and higher proton conductivity at higher operating temperatures like 500-800 °C. For the practical application of the functional BaZrO3 proton conductor in the electrochemical hydrogen devices, its necessary to understand the isotopic effect of proton conductivity. To understand the isotopic effect of proton conductivity in the barium zirconates, in this study, the proton conductivity in the Ar, (Ar + 4% H2), (Ar + 4% D2), (Ar + H2O), (Ar + D2O), and O2 atmospheres were measured for two different compositions: BaZr0.9Y0.1O2.95 (BZY), and BaZr0.955Y0.03Co0.015O2.97 (BZYC) in the temperature range from 500 °C to 1000 °C. By comparing the obtained results, a significant difference in sinterability, conductivity, and the isotopic effect was observed due to the co-doping of the Co element in the BaZr1-xYxO3-a proton conductor.

preprint2021arXiv

Observation of the gamma-ray binary HESS J0632+057 with the H.E.S.S., MAGIC, and VERITAS telescopes

The results of gamma-ray observations of the binary system HESS J0632+057 collected during 450 hours over 15 years, between 2004 and 2019, are presented. Data taken with the atmospheric Cherenkov telescopes H.E.S.S., MAGIC, and VERITAS at energies above 350 GeV were used together with observations at X-ray energies obtained with Swift-XRT, Chandra, XMM-Newton, NuSTAR, and Suzaku. Some of these observations were accompanied by measurements of the Hα emission line. A significant detection of the modulation of the VHE gamma-ray fluxes with a period of 316.7+-4.4 days is reported, consistent with the period of 317.3+-0.7 days obtained with a refined analysis of X-ray data. The analysis of data of four orbital cycles with dense observational coverage reveals short timescale variability, with flux-decay timescales of less than 20 days at very high energies. Flux variations observed over the time scale of several years indicate orbit-to-orbit variability. The analysis confirms the previously reported correlation of X-ray and gamma-ray emission from the system at very high significance, but can not find any correlation of optical Hα parameters with X-ray or gamma-ray energy fluxes in simultaneous observations. The key finding is that the emission of HESS J0632+057 in the X-ray and gamma-ray energy bands is highly variable on different time scales. The ratio of gamma-ray to X-ray flux shows the equality or even dominance of the gamma-ray energy range. This wealth of new data is interpreted taking into account the insufficient knowledge of the ephemeris of the system, and discussed in the context of results reported on other gamma-ray binary systems.

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.

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

Development and application of a $^3$He Neutron Spin Filter at J-PARC

We are developing a neutron polarizer with polarized $^3$He gas, referred to as a $^3$He spin filter, based on the Spin Exchange Optical Pumping (SEOP) for polarized neutron scattering experiments at Materials and Life Science Experimental Facility (MLF) of Japan Proton Accelerator Research Complex (J-PARC). A $^3$He gas-filling station was constructed at J-PARC, and several $^3$He cells with long spin relaxation times have been fabricated using the gas-filling station. A laboratory has been prepared in the MLF beam hall for polarizing $^3$He cells, and compact pumping systems with laser powers of 30~W and 110~W, which can be installed onto a neutron beamline, have been developed. A $^3$He polarization of 85% was achieved at a neutron beamline by using the pumping system with the 110~W laser. Recently, the first user experiment utilizing the $^3$He spin filter was conducted, and there have been several more since then. The development and utilization of $^3$He spin filters at MLF of J-PARC are reported.

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.

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

Anomalous time delays and quantum weak measurements in optical micro-resonators

We study inelastic resonant scattering of a Gaussian wave packet with the parameters close to a zero of the complex scattering coefficient. We demonstrate, both theoretically and experimentally, that such near-zero scattering can result in anomalously-large time delays and frequency shifts of the scattered wave packet. Furthermore, we reveal a close analogy of these anomalous shifts with the spatial and angular Goos-Hänchen optical beam shifts, which are amplified via quantum weak measurements. However, in contrast to other beam-shift and weak-measurement systems, we deal with a one-dimensional scalar wave without any intrinsic degrees of freedom. It is the non-Hermitian nature of the system that produces its rich and non-trivial behaviour. Our results are generic for any scattering problem, either quantum or classical. As an example, we consider the transmission of an optical pulse through a nano-fiber with a side-coupled toroidal micro-resonator. The zero of the transmission coefficient corresponds to the critical coupling conditions. Experimental measurements of the time delays near the critical-coupling parameters verify our weak-measurement theory and demonstrate amplification of the time delay from the typical inverse resonator linewidth scale to the pulse duration scale.

preprint2016arXiv

Change of nuclear configurations in the neutrinoless double-$β$ decay of $^{130}$Te $\rightarrow$ $^{130}$Xe and $^{136}$Xe $\rightarrow$ $^{136}$Ba

The change in the configuration of valence protons between the initial and final states in the neutrinoless double-$β$ decay of $^{130}$Te $\rightarrow$ $^{130}$Xe and of $^{136}$Xe $\rightarrow$ $^{136}$Ba has been determined by measuring the cross sections of the ($d$,$^3$He) reaction with 101-MeV deuterons. Together with our recent determination of the relevant neutron configurations involved in the process, a quantitative comparison with the latest shell-model and interacting-boson-model calculations reveals significant discrepancies. These are the same calculations used to determine the nuclear matrix elements governing the rate of neutrinoless double-$β$ decay in these systems.

preprint2016arXiv

Charge-changing-cross-section measurements of $^{12-16}$C at around $45A$ MeV and development of a Glauber model for incident energies $10A-2100A$ MeV

We have measured for the first time the charge-changing cross sections ($σ_{\text{CC}}$) of $^{12-16}$C on a $^{12}$C target at energies below $100A$ MeV. To analyze these low-energy data, we have developed a finite-range Glauber model with a global parameter set within the optical-limit approximation which is applicable to reaction cross section ($σ_{\text{R}}$) and $σ_{\text{CC}}$ measurements at incident energies from 10$A$ to $2100A$ MeV. Adopting the proton-density distribution of $^{12}$C known from the electron-scattering data, as well as the bare total nucleon-nucleon cross sections, and the real-to-imaginary-part ratios of the forward proton-proton elastic scattering amplitude available in the literatures, we determine the energy-dependent slope parameter $β_{\rm pn}$ of the proton-neutron elastic differential cross section so as to reproduce the existing $σ_{\text{R}}$ and interaction-cross-section data for $^{12}$C+$^{12}$C over a wide range of incident energies. The Glauber model thus formulated is applied to calculate the $σ_{\text{\tiny R}}$'s of $^{12}$C on a $^9$Be and $^{27}$Al targets at various incident energies. Our calculations show excellent agreement with the experimental data. Applying our model to the $σ_{\text{\tiny R}}$ and $σ_{\text{\tiny CC}}$ for the "neutron-skin" $^{16}$C nucleus, we reconfirm the importance of measurements at incident energies below $100A$ MeV. The proton root-mean-square radii of $^{12-16}$C are extracted using the measured $σ_{\text{CC}}$'s and the existing $σ_{\text{R}}$ data. The results for $^{12-14}$C are consistent with the values from the electron scatterings, demonstrating the feasibility, usefulness of the $σ_{\text{CC}}$ measurement and the present Glauber model.

preprint2016arXiv

Displacement of propagating squeezed microwave states

Displacement of propagating quantum states of light is a fundamental operation for quantum communication. It enables fundamental studies on macroscopic quantum coherence and plays an important role in quantum teleportation protocols with continuous variables. In our experiments we have successfully implemented this operation for propagating squeezed microwave states. We demonstrate that, even for strong displacement amplitudes, there is no degradation of the squeezing level in the reconstructed quantum states. Furthermore, we confirm that path entanglement generated by using displaced squeezed states stays constant over a wide range of the displacement power.

preprint2016arXiv

High Resolution Spectroscopic Study of $^{10}_Λ$Be

Spectroscopy of a $^{10}_Λ$Be hypernucleus was carried out at JLab Hall C using the $(e,e^{\prime}K^{+})$ reaction. A new magnetic spectrometer system (SPL+HES+HKS), specifically designed for high resolution hypernuclear spectroscopy, was used to obtain an energy spectrum with a resolution of 0.78 MeV (FWHM). The well-calibrated spectrometer system of the present experiment using the $p(e,e^{\prime}K^{+})Λ,Σ^{0}$ reactions allowed us to determine the energy levels, and the binding energy of the ground state peak (mixture of 1$^{-}$ and 2$^{-}$ states) was obtained to be B$_Λ$=8.55$\pm$0.07(stat.)$\pm$0.11(sys.) MeV. The result indicates that the ground state energy is shallower than that of an emulsion study by about 0.5 MeV which provides valuable experimental information on charge symmetry breaking effect in the $ΛN$ interaction.

preprint2016arXiv

Spectroscopy of the neutron-rich hypernucleus $^{7}_Λ$He from electron scattering

The missing mass spectroscopy of the $^{7}_Λ$He hypernucleus was performed, using the $^{7}$Li$(e,e^{\prime}K^{+})^{7}_Λ$He reaction at the Thomas Jefferson National Accelerator Facility Hall C. The $Λ$ binding energy of the ground state (1/2$^{+}$) was determined with a smaller error than that of the previous measurement, being $B_Λ$ = 5.55 $\pm$ 0.10(stat.) $\pm$ 0.11(sys.) MeV. The experiment also provided new insight into charge symmetry breaking in p-shell hypernuclear systems. Finally, a peak at $B_Λ$ = 3.65 $\pm$ 0.20(stat.) $\pm$ 0.11(sys.) MeV was observed and assigned as a mixture of 3/2$^{+}$ and 5/2$^{+}$ states, confirming the "gluelike" behavior of $Λ$, which makes an unstable state in $^{6}$He stable against neutron emission.

preprint2015arXiv

Correlated emission lasing in harmonic oscillators coupled via a single three-level artificial atom

A single superconducting artificial atom provides a unique basis for coupling electromagnetic fields and photons hardly achieved with a natural atom. Bringing a pair of harmonic oscillators into resonance with transitions of the three-level atom converts atomic spontaneous processes into correlated emission dynamics. We demonstrate two-mode correlated emission lasing on harmonic oscillators coupled via the fully controllable three-level artificial atom. Correlation of two different color emissions reveals itself as equally narrowed linewiths and quench of their mutual phase-diffusion. The mutual linewidth is more than four orders of magnitude narrower than the Schawlow-Townes limit. The interference between the different color lasing fields demonstrates the two-mode fields are strongly correlated.

preprint2015arXiv

Measurement and comparison of individual external doses of high-school students living in Japan, France, Poland and Belarus -- the "D-shuttle" project --

Twelve high schools in Japan (of which six are in Fukushima Prefecture), four in France, eight in Poland and two in Belarus cooperated in the measurement and comparison of individual external doses in 2014. In total 216 high-school students and teachers participated in the study. Each participant wore an electronic personal dosimeter "D-shuttle" for two weeks, and kept a journal of his/her whereabouts and activities. The distributions of annual external doses estimated for each region overlap with each other, demonstrating that the personal external individual doses in locations where residence is currently allowed in Fukushima Prefecture and in Belarus are well within the range of estimated annual doses due to the background radiation level of other regions/countries.

preprint2015arXiv

Probing astrophysically important states in $^{26}$Mg nucleus to study neutron sources for the $s$-Process

The $^{22}$Ne($α$,n)$^{25}$Mg reaction is the dominant neutron source for the slow neutron capture process ($s$-process) in massive stars and contributes, together with the $^{13}$C($α$,n)$^{16}$O, to the production of neutrons for the $s$-process in Asymptotic Giant Branch (AGB) stars. However, the reaction is endothermic and competes directly with the $^{22}$Ne($α,γ)^{26}$Mg radiative capture. The uncertainties for both reactions are large owing to the uncertainty in the level structure of $^{26}$Mg near the alpha and neutron separation energies. These uncertainties are affecting the s-process nucleosynthesis calculations in theoretical stellar models. Indirect studies in the past have been successful in determining the energies, $γ$-ray and neutron widths of the $^{26}$Mg states in the energy region of interest. But, the high Coulomb barrier hinders a direct measurement of the resonance strengths, which are determined by the $α$-widths for these states. The goal of the present experiments is to identify the critical resonance states and to precisely measure the $α$-widths by $α$ transfer techniques . Hence, the $α$-inelastic scattering and $α$-transfer measurements were performed on a solid $^{26}$Mg target and a $^{22}$Ne gas target, respectively, using the Grand Raiden Spectrometer at RCNP, Osaka, Japan. Six levels (E$_x$ = 10717 keV , 10822 keV, 10951 keV, 11085 keV, 11167 keV and 11317 keV) have been observed above the $α$-threshold in the region of interest (10.61 - 11.32 MeV). The rates are dominated in both reaction channels by the resonance contributions of the states at E$_x$ = 10951, 11167 and 11317 keV. The E$_x$ =11167 keV has the most appreciable impact on the ($α,γ$) rate and therefore plays an important role for the prediction of the neutron production in s-process environments.

preprint2015arXiv

Quadratic Fermi Node in a 3D Strongly Correlated Semimetal

Strong spin-orbit coupling fosters exotic electronic states such as topological insulators and superconductors, but the combination of strong spin-orbit and strong electron-electron interactions is just beginning to be understood. Central to this emerging area are the 5d transition metal iridium oxides. Here, in the pyrochlore iridate Pr2Ir2O7, we identify a nontrivial state with a single point Fermi node protected by cubic and time-reversal symmetries, using a combination of angle-resolved photoemission spectroscopy and first principles calculations. Owing to its quadratic dispersion, the unique coincidence of four degenerate states at the Fermi energy, and strong Coulomb interactions, non-Fermi liquid behavior is predicted, for which we observe some evidence. Our discovery implies that Pr2Ir2O7 is a parent state that can be manipulated to produce other strongly correlated topological phases, such as topological Mott insulator, Weyl semi-metal, and quantum spin and anomalous Hall states.

preprint2015arXiv

Simultaneous Observation of Solar Neutrons from the ISS and High Mountain Observatories in association with a flare on July 8, 2014

An M6.5-class flare was observed at N12E56 of the solar surface at 16:06 UT on July 8, 2014. In association with this flare, solar neutron detectors located on two high mountains, Mt. Sierra Negra and Chacaltaya and at the space station observed enhancements in the neutral channel. The authors analysed these data and a possible scenario of enhancements produced by high-energy protons and neutrons is proposed, using the data from continuous observation of a solar surface by the ultraviolet telescope onboard the Solar Dynamical Observatory (SDO).

preprint2015arXiv

The Search for >35 MeV Neutrons from the June 3, 2012 Impulsive Flare

We analyzed a highly impulsive solar flare observed on June 3, 2012. In association with this flare, emissions of hard X-rays, high-energy gamma rays, and neutrons were detected by the detectors onboard the FERMI, RHESSI satellites and the International Space Station. We compared those results with the pictures taken by the UV telescope onboard the Solar Dynamics Observatory satellite and found the crossing structure of two magnetic ropes at two positions on the solar surface almost at the same time. High-energy gamma rays were detected by the Fermi Large Area Telescope satellite, implying that the impulsive flare was one of a major source of proton acceleration processes on the solar surface. At the beginning of research, impulsive solar flares were considered to be the main source of particle acceleration processes; our current observations have confirmed this hypothesis.

preprint2014arXiv

A search for point sources of EeV photons

Measurements of air showers made using the hybrid technique developed with the fluorescence and surface detectors of the Pierre Auger Observatory allow a sensitive search for point sources of EeV photons anywhere in the exposed sky. A multivariate analysis reduces the background of hadronic cosmic rays. The search is sensitive to a declination band from -85° to +20°, in an energy range from 10^17.3 eV to 10^18.5 eV. No photon point source has been detected. An upper limit on the photon flux has been derived for every direction. The mean value of the energy flux limit that results from this, assuming a photon spectral index of -2, is 0.06 eV cm^-2 s^-1, and no celestial direction exceeds 0.25 eV cm^-2 s^-1. These upper limits constrain scenarios in which EeV cosmic ray protons are emitted by non-transient sources in the Galaxy.

preprint2014arXiv

A Targeted Search for Point Sources of EeV Neutrons

A flux of neutrons from an astrophysical source in the Galaxy can be detected in the Pierre Auger Observatory as an excess of cosmic-ray air showers arriving from the direction of the source. To avoid the statistical penalty for making many trials, classes of objects are tested in combinations as nine "target sets", in addition to the search for a neutron flux from the Galactic Center or from the Galactic Plane. Within a target set, each candidate source is weighted in proportion to its electromagnetic flux, its exposure to the Auger Observatory, and its flux attenuation factor due to neutron decay. These searches do not find evidence for a neutron flux from any class of candidate sources. Tabulated results give the combined p-value for each class, with and without the weights, and also the flux upper limit for the most significant candidate source within each class. These limits on fluxes of neutrons significantly constrain models of EeV proton emission from non-transient discrete sources in the Galaxy.

preprint2014arXiv

Isotopic identification of engineered nitrogen-vacancy spin qubits in ultrapure diamond

Nitrogen impurities help to stabilize the negatively-charged-state of NV$^-$ in diamond, whereas magnetic fluctuations from nitrogen spins lead to decoherence of NV$^-$ qubits. It is not known what donor concentration optimizes these conflicting requirements. Here we used 10-MeV $^{15}$N$^{3+}$ ion implantation to create NV$^-$ in ultrapure diamond. Optically detected magnetic resonance of single centers revealed a high creation yield of $40\pm3$% from $^{15}$N$^{3+}$ ions and an additional yield of $56\pm3$% from $^{14}$N impurities. High-temperature anneal was used to reduce residual defects, and charge stable NV$^-$, even in a dilute $^{14}$N impurity concentration of 0.06 ppb were created with long coherence times.

preprint2014arXiv

Josephson parametric phase-locked oscillator and its application to dispersive readout of superconducting qubits

The parametric phase-locked oscillator (PPLO), also known as a parametron, is a resonant circuit in which one of the reactances is periodically modulated. It can detect, amplify, and store binary digital signals in the form of two distinct phases of self-oscillation. Indeed, digital computers using PPLOs based on a magnetic ferrite ring or a varactor diode as its fundamental logic element were successfully operated in 1950s and 1960s. More recently, basic bit operations have been demonstrated in an electromechanical resonator, and an Ising machine based on optical PPLOs has been proposed. Here, using a PPLO realized with Josephson-junction circuitry, we demonstrate the demodulation of a microwave signal digitally modulated by binary phase-shift keying. Moreover, we apply this demodulation capability to the dispersive readout of a superconducting qubit. This readout scheme enables a fast and latching-type readout, yet requires only a small number of readout photons in the resonator to which the qubit is coupled, thus featuring the combined advantages of several disparate schemes. We have achieved high-fidelity, single-shot, and non-destructive qubit readout with Rabi-oscillation contrast exceeding 90%, limited primarily by the qubit's energy relaxation.

preprint2014arXiv

Microwave Down-Conversion with an Impedance-Matched $Λ$ System in Driven Circuit QED

By driving a dispersively coupled qubit-resonator system, we realize an "impedance-matched" $Λ$ system that has two identical radiative decay rates from the top level and interacts with a semi-infinite waveguide. It has been predicted that a photon input from the waveguide deterministically induces a Raman transition in the system and switches its electronic state. We confirm this through microwave response to a continuous probe field, observing near-perfect ($99.7\%$) extinction of the reflection and highly efficient ($74\%$) frequency down-conversion. These proof-of-principle results lead to deterministic quantum gates between material qubits and microwave photons and open the possibility for scalable quantum networks interconnected with waveguide photons.

preprint2014arXiv

Reconstruction of inclined air showers detected with the Pierre Auger Observatory

We describe the method devised to reconstruct inclined cosmic-ray air showers with zenith angles greater than $60^\circ$ detected with the surface array of the Pierre Auger Observatory. The measured signals at the ground level are fitted to muon density distributions predicted with atmospheric cascade models to obtain the relative shower size as an overall normalization parameter. The method is evaluated using simulated showers to test its performance. The energy of the cosmic rays is calibrated using a sub-sample of events reconstructed with both the fluorescence and surface array techniques. The reconstruction method described here provides the basis of complementary analyses including an independent measurement of the energy spectrum of ultra-high energy cosmic rays using very inclined events collected by the Pierre Auger Observatory.

preprint2014arXiv

Resonant Spin Wave Excitation by Terahertz Magnetic Near-field Enhanced with Split Ring Resonator

Excitation of antiferromagnetic spin waves in HoFeO$_{3}$ crystal combined with a split ring resonator (SRR) is studied using terahertz (THz) electromagnetic pulses. The magnetic field in the vicinity of the SRR induced by the incident THz electric field component excites and the Faraday rotation of the polarization of a near-infrared probe pulse directly measures oscillations that correspond to the antiferromagnetic spin resonance mode. The good agreement of the temperature-dependent magnetization dynamics with the calculation using the two-lattice Landau-Lifshitz-Gilbert equation confirms that the spin wave is resonantly excited by the THz magnetic near-field enhanced at the LC resonance frequency of the SRR, which is 20 times stronger than the incident magnetic field.

preprint2014arXiv

Revisiting Jovian-Resonance Induced Chondrule Formation

It is proposed that planetesimals perturbed by Jovian mean-motion resonances are the source of shock waves that form chondrules. It is considered that this shock-induced chondrule formation requires the velocity of the planetesimal relative to the gas disk to be on the order of > 7 km/s at 1 AU. In previous studies on planetesimal excitation, the effects of Jovian mean-motion resonance together with the gas drag were investigated, but the velocities obtained were at most 8 km/s in the asteroid belt, which is insufficient to account for the ubiquitous existence of chondrules. In this paper, we reexamine the effect of Jovian resonances and take into account the secular resonance in the asteroid belt caused by the gravity of the gas disk. We find that the velocities relative to the gas disk of planetesimals a few hundred kilometers in size exceed 12 km/s, and that this is achieved around the 3:1 mean-motion resonance. The heating region is restricted to a relatively narrow band between 1.5 AU and 3.5 AU. Our results suggest that chondrules were produced effectively in the asteroid region after Jovian formation. We also find that many planetesimals are scattered far beyond Neptune. Our findings can explain the presence of crystalline silicate in comets if the scattered planetesimals include silicate dust processed by shock heating.

preprint2014arXiv

Searches for Large-Scale Anisotropy in the Arrival Directions of Cosmic Rays Detected above Energy of $10^{19}$ eV at the Pierre Auger Observatory and the Telescope Array

Spherical harmonic moments are well-suited for capturing anisotropy at any scale in the flux of cosmic rays. An unambiguous measurement of the full set of spherical harmonic coefficients requires full-sky coverage. This can be achieved by combining data from observatories located in both the northern and southern hemispheres. To this end, a joint analysis using data recorded at the Telescope Array and the Pierre Auger Observatory above $10^{19}$ eV is presented in this work. The resulting multipolar expansion of the flux of cosmic rays allows us to perform a series of anisotropy searches, and in particular to report on the angular power spectrum of cosmic rays above $10^{19}$ eV. No significant deviation from isotropic expectations is found throughout the analyses performed. Upper limits on the amplitudes of the dipole and quadrupole moments are derived as a function of the direction in the sky, varying between 7% and 13% for the dipole and between 7% and 10% for a symmetric quadrupole.

preprint2014arXiv

The experiments with the High Resolution Kaon Spectrometer at JLab Hall C and the new spectroscopy of ${}^{12}_Λ\text{B}$ hypernuclei

Since the pioneering experiment, E89-009 studying hypernuclear spectroscopy using the $(e,e^{\prime}K^+)$ reaction was completed, two additional experiments, E01-011 and E05-115, were performed at Jefferson Lab. These later experiments used a modified experimental design, the "tilt method", to dramatically suppress the large electromagnetic background, and allowed for a substantial increase in luminosity. Additionally, a new kaon spectrometer, HKS (E01-011), a new electron spectrometer, HES, and a new splitting magnet were added to produce precision, high-resolution hypernuclear spectroscopy. These two experiments, E01-011 and E05-115, resulted in two new data sets, producing sub-MeV energy resolution in the spectra of ${}^{7}_Λ\text{He}$, ${}^{12}_Λ\text{B}$ and ${}^{28}_Λ\text{Al}$ and ${}^{7}_Λ\text{He}$, ${}^{10}_Λ\text{Be}$, ${}^{12}_Λ\text{B}$ and ${}^{52}_Λ\text{V}$. All three experiments obtained a ${}^{12}_Λ\text{B}$, spectrum, which is the most characteristic $p$-shell hypernucleus and is commonly used for calibration. Independent analyses of these different experiments demonstrate excellent consistency and provide the clearest level structure to date of this hypernucleus as produced by the $(e,e^{\prime}K^+)$ reaction. This paper presents details of these experiments, and the extraction and analysis of the observed ${}^{12}_Λ\text{B}$ spectrum.

preprint2013arXiv

An evaluation of the exposure in nadir observation of the JEM-EUSO mission

We evaluate the exposure during nadir observations with JEM-EUSO, the Extreme Universe Space Observatory, on-board the Japanese Experiment Module of the International Space Station. Designed as a mission to explore the extreme energy Universe from space, JEM-EUSO will monitor the Earth's nighttime atmosphere to record the ultraviolet light from tracks generated by extensive air showers initiated by ultra-high energy cosmic rays. In the present work, we discuss the particularities of space-based observation and we compute the annual exposure in nadir observation. The results are based on studies of the expected trigger aperture and observational duty cycle, as well as, on the investigations of the effects of clouds and different types of background light. We show that the annual exposure is about one order of magnitude higher than those of the presently operating ground-based observatories.

preprint2013arXiv

Broadband Excitation by Chirped Pulses: Application to Single Electron Spins in Diamond

Pulsed excitation of broad spectra requires very high field strengths if monochromatic pulses are used. If the corresponding high power is not available or not desirable, the pulses can be replaced by suitable low-power pulses that distribute the power over a wider bandwidth. As a simple case, we use microwave pulses with a linear frequency chirp. We use these pulses to excite spectra of single NV-centers in a Ramsey experiment. Compared to the conventional Ramsey experiment, our approach increases the bandwidth by at least an order of magnitude. Compared to the conventional ODMR experiment, the chirped Ramsey experiment does not suffer from power broadening and increases the resolution by at least an order of magnitude. As an additional benefit, the chirped Ramsey spectrum contains not only `allowed' single quantum transitions, but also `forbidden' zero- and double quantum transitions, which can be distinguished from the single quantum transitions by phase-shifting the readout pulse with respect to the excitation pulse or by variation of the external magnetic field strength.

preprint2013arXiv

Collectivity evolution in the neutron-rich Pd isotopes towards the N=82 shell closure

The neutron-rich, even-even 122,124,126Pd isotopes has been studied via in-beam gamma-ray spectroscopy at the RIKEN Radioactive Isotope Beam Factory. Excited states at 499(9), 590(11), and 686(17) keV were found in the three isotopes, which we assign to the respective 2+ -> 0+ decays. In addition, a candidate for the 4+ state at 1164(20) keV was observed in 122Pd. The resulting Ex(2+) systematics are essentially similar to those of the Xe (Z=54) isotopic chain and theoretical prediction by IBM-2, suggesting no serious shell quenching in the Pd isotopes in the vicinity of N=82.

preprint2013arXiv

Direct Measurements of the Lifetime of Heavy Hypernuclei

The lifetime of a Lambda particle embedded in a nucleus (hypernucleus) decreases from that of free Lambda decay due to the opening of the Lambda N to NN weak decay channel. However, it is generally believed that the lifetime of a hypernucleus attains a constant value (saturation) for medium to heavy hypernuclear masses, yet this hypothesis has been difficult to verify. The present paper reports a direct measurement of the lifetime of medium-heavy hypernuclei produced with a photon-beam from Fe, Cu, Ag, and Bi targets. The recoiling hypernuclei were detected by a fission fragment detector using low-pressure multi-wire proportional chambers. The experiment agrees remarkably well with the only previously-measured single-species heavy-hypernucleus lifetime, that of Fe56_Lambda at KEK, and has significantly higher precision. The experiment disagrees with the measured lifetime of an unknown combination of heavy hypernuclei with 180<A<225 and, with a small statistical and systematic uncertainty, strongly favors the expected saturation of the lifetime decrease.

preprint2013arXiv

Experimental implementation of assisted quantum adiabatic passage in a single spin

Quantum adiabatic passages can be greatly accelerated by a suitable control field, called a counter-diabatic field, which varies during the scan through resonance. Here, we implement this technique on the electron spin of a single nitrogen-vacancy center in diamond. We demonstrate two versions of this scheme. The first follows closely the procedure originally proposed by Demirplak and Rice (J. Phys. Chem. A 107, 9937 (2003)). In the second scheme, we use a control field whose amplitude is constant, but its phase varies with time. This version, which we call the rapid-scan approach, allows an even faster passage through resonance and therefore makes it applicable also for systems with shorter decoherence times.

preprint2013arXiv

Extending spin coherence times of diamond qubits by high temperature annealing

Spins of negatively charged nitrogen-vacancy (NV$^-$) defects in diamond are among the most promising candidates for solid-state qubits. The fabrication of quantum devices containing these spin-carrying defects requires position-controlled introduction of NV$^-$ defects having excellent properties such as spectral stability, long spin coherence time, and stable negative charge state. Nitrogen ion implantation and annealing enable the positioning of NV$^-$ spin qubits with high precision, but to date, the coherence times of qubits produced this way are short, presumably because of the presence of residual radiation damage. In the present work, we demonstrate that a high temperature annealing at 1000$^\circ$C allows 2 millisecond coherence times to be achieved at room temperature. These results were obtained for implantation-produced NV$^-$ defects in a high-purity, 99.99% $^{12}$C enriched single crystal chemical vapor deposited diamond. We discuss these remarkably long coherence times in the context of the thermal behavior of residual defect spins. [Published in Physical Review B {\bf{88}}, 075206 (2013)]

preprint2013arXiv

Measurement of High-Energy Solar Neutrons with SEDA-FIB onboard the ISS

A new type of solar neutron detector (SEDA-FIB) was launched on board the Space Shuttle Endeavor on July 16 2009, and began collecting data at the International Space Station (ISS) on August 25 2009. This paper summarizes four years of observations with the solar neutron detector SEDA-FIB (Space Environment Data Acquisition using the FIBer detector). The solar neutron detector FIB can determine both the energy and arrival direction of solar neutrons. In this paper, we first present the angular distribution of neutron induced protons obtained in Monte Carlo simulations. The results are compared with the experimental results. Then we provide the angular distribution of background neutrons during one full orbit of the ISS (90 minutes). Next, the angular distribution of neutrons during the flare onset time from 20:02 to 20:10 UT on March 7 2011 is presented. It is compared with the distribution when a solar flare is not occurring. Observed solar neutrons possibly originated from the M-class solar flares that occurred on March 7 (M3.7), June 7 (M2.5), September 24 (M3.0) (weak signal) and November 3 (X1.9) of 2011 and January 23 of 2012 (M8.7). This marked the first time that neutrons have been observed from M-class solar flares. A possible interpretation of the neutron production process will be also provided.

preprint2013arXiv

Observation of three-state dressed states in circuit quantum electrodynamics

We have investigated the microwave response of a transmon qubit coupled directly to a transmission line. In a transmon qubit, owing to its weak anharmonicity, a single driving field may generate dressed states involving more than two bare states. We confirmed the formation of three-state dressed states by observing all of the six associated Rabi sidebands, which appear as either amplification or attenuation of the probe field. The experimental results are reproduced with good precision by a theoretical model incorporating the radiative coupling between the qubit and the microwave.

preprint2013arXiv

Powerful terahertz emission from Bi2Sr2CaCu2O8+delta mesa arrays

Stacks of intrinsic Josephson junctions (IJJs) in high-temperature superconductors enable the fabrication of compact sources of coherent terahertz radiation. Here we demonstrate that multiple stacks patterned on the same Bi2Sr2CaCu2O8+delta crystal can - under optimized conditions - be synchronized to emit high-power THz-radiation. For three synchronized stacks we achieved 610 microwatts of continuous-wave coherent radiation power at 0.51 THz. We suggest that synchronization is promoted by THz-waves in the base crystal. We note though that synchronization cannot be achieved in all samples. However, even in these cases powers on the 100-microwatt scale can be generated.

preprint2013arXiv

Single-shot readout of a superconducting flux qubit with a flux-driven Josephson parametric amplifier

We report single-shot readout of a superconducting flux qubit by using a flux-driven Josephson parametric amplifier (JPA). After optimizing the readout power, gain of the JPA and timing of the data acquisition, we observe the Rabi oscillations with a contrast of 74% which is mainly limited by the bandwidth of the JPA and the energy relaxation of the qubit. The observation of quantum jumps between the qubit eigenstates under continuous monitoring indicates the nondestructiveness of the readout scheme.

preprint2013arXiv

Squeezing with a flux-driven Josephson parametric amplifier

Josephson parametric amplifiers (JPA) are promising devices for applications in circuit quantum electrodynamics (QED) and for studies on propagating quantum microwaves because of their good noise performance. In this work, we present a systematic characterization of a flux-driven JPA at millikelvin temperatures. In particular, we study in detail its squeezing properties by two different detection techniques. With the homodyne setup, we observe squeezing of vacuum fluctuations by superposing signal and idler bands. For a quantitative analysis we apply dual-path cross-correlation techniques to reconstruct the Wigner functions of various squeezed vacuum and thermal states. At 10 dB signal gain, we find 4.9+-0.2 dB squeezing below vacuum. In addition, we discuss the physics behind squeezed coherent microwave fields. Finally, we analyze the JPA noise temperature in the degenerate mode and find a value smaller than the standard quantum limit for phase-insensitive amplifiers.

preprint2012arXiv

A search for ultra-high energy neutrinos in highly inclined events at the Pierre Auger Observatory

The Surface Detector of the Pierre Auger Observatory is sensitive to neutrinos of all flavours above 0.1 EeV. These interact through charged and neutral currents in the atmosphere giving rise to extensive air showers. When interacting deeply in the atmosphere at nearly horizontal incidence, neutrinos can be distinguished from regular hadronic cosmic rays by the broad time structure of their shower signals in the water-Cherenkov detectors. In this paper we present for the first time an analysis based on down-going neutrinos. We describe the search procedure, the possible sources of background, the method to compute the exposure and the associated systematic uncertainties. No candidate neutrinos have been found in data collected from 1 January 2004 to 31 May 2010. Assuming an E^-2 differential energy spectrum the limit on the single flavour neutrino is (E^2 * dN/dE) < 1.74x10^-7 GeV cm^-2 s^-1 sr^-1 at 90% C.L. in the energy range 1x10^17 eV < E < 1x10^20 eV.

preprint2012arXiv

Antennas for the Detection of Radio Emission Pulses from Cosmic-Ray induced Air Showers at the Pierre Auger Observatory

The Pierre Auger Observatory is exploring the potential of the radio detection technique to study extensive air showers induced by ultra-high energy cosmic rays. The Auger Engineering Radio Array (AERA) addresses both technological and scientific aspects of the radio technique. A first phase of AERA has been operating since September 2010 with detector stations observing radio signals at frequencies between 30 and 80 MHz. In this paper we present comparative studies to identify and optimize the antenna design for the final configuration of AERA consisting of 160 individual radio detector stations. The transient nature of the air shower signal requires a detailed description of the antenna sensor. As the ultra-wideband reception of pulses is not widely discussed in antenna literature, we review the relevant antenna characteristics and enhance theoretical considerations towards the impulse response of antennas including polarization effects and multiple signal reflections. On the basis of the vector effective length we study the transient response characteristics of three candidate antennas in the time domain. Observing the variation of the continuous galactic background intensity we rank the antennas with respect to the noise level added to the galactic signal.

preprint2012arXiv

Description of Atmospheric Conditions at the Pierre Auger Observatory using the Global Data Assimilation System (GDAS)

Atmospheric conditions at the site of a cosmic ray observatory must be known for reconstructing observed extensive air showers. The Global Data Assimilation System (GDAS) is a global atmospheric model predicated on meteorological measurements and numerical weather predictions. GDAS provides altitude-dependent profiles of the main state variables of the atmosphere like temperature, pressure, and humidity. The original data and their application to the air shower reconstruction of the Pierre Auger Observatory are described. By comparisons with radiosonde and weather station measurements obtained on-site in Malargüe and averaged monthly models, the utility of the GDAS data is shown.

preprint2012arXiv

Electron spin resonance detected by a superconducting qubit

A new method for detecting the magnetic resonance of electronic spins at low temperature is demonstrated. It consists in measuring the signal emitted by the spins with a superconducting qubit that acts as a single-microwave-photon detector, resulting in an enhanced sensitivity. We implement this new type of electron-spin resonance spectroscopy using a hybrid quantum circuit in which a transmon qubit is coupled to a spin ensemble consisting of NV centers in diamond. With this setup we measure the NV center absorption spectrum at 30mK at an excitation level of \thicksim15\,μ_{B} out of an ensemble of 10^{11} spins.

preprint2012arXiv

Large Dispersive Shift of Cavity Resonance Induced by a Superconducting Flux Qubit in the Straddling Regime

We demonstrate enhancement of the dispersive frequency shift in a coplanar waveguide resonator induced by a capacitively-coupled superconducting flux qubit in the straddling regime. The magnitude of the observed shift, 80 MHz for the qubit-resonator detuning of 5 GHz, is quantitatively explained by the generalized Jaynes-Cummings model which takes into account the contribution of the qubit higher energy levels. By applying the enhanced dispersive shift to the qubit readout, we achieved 90% contrast of the Rabi oscillations which is mainly limited by the energy relaxation of the qubit.

preprint2012arXiv

Measurement by FIB on the ISS: Two Emissions of Solar Neutrons Detected?

A new type of solar neutron detector (FIB) was launched onboard the Space Shuttle Endeavour on July 16, 2009, and it began collecting data at the International Space Station (ISS) on August 25, 2009. This paper summarizes the three years of observations obtained by the solar neutron detector FIB until the end of July 2012. The solar neutron detector FIB can determine both the energy and arrival direction of neutrons. We measured the energy spectra of background neutrons over the SAA region and elsewhere, and found the typical trigger rates to be 20 counts/sec and 0.22 counts/sec, respectively. It is possible to identify solar neutrons to within a level of 0.028 counts/sec, provided that directional information is applied. Solar neutrons were observed in association with the M-class solar flares that occurred on March 7 (M3.7) and June 7 (M2.5) of 2011. This marked the first time that neutrons were observed in M-class solar flares. A possible interpretaion of the prodcution process is provided.

preprint2012arXiv

Path Entanglement of Continuous-Variable Quantum Microwaves

Path entanglement constitutes an essential resource in quantum information and communication protocols. Here, we demonstrate frequency-degenerate entanglement between continuous-variable quantum microwaves propagating along two spatially separated paths. We combine a squeezed and a vacuum state using a microwave beam splitter. Via correlation measurements, we detect and quantify the path entanglement contained in the beam splitter output state. Our experiments open the avenue to quantum teleportation, quantum communication, or quantum radar with continuous variables at microwave frequencies.

preprint2012arXiv

The Rapid Atmospheric Monitoring System of the Pierre Auger Observatory

The Pierre Auger Observatory is a facility built to detect air showers produced by cosmic rays above 10^17 eV. During clear nights with a low illuminated moon fraction, the UV fluorescence light produced by air showers is recorded by optical telescopes at the Observatory. To correct the observations for variations in atmospheric conditions, atmospheric monitoring is performed at regular intervals ranging from several minutes (for cloud identification) to several hours (for aerosol conditions) to several days (for vertical profiles of temperature, pressure, and humidity). In 2009, the monitoring program was upgraded to allow for additional targeted measurements of atmospheric conditions shortly after the detection of air showers of special interest, e.g., showers produced by very high-energy cosmic rays or showers with atypical longitudinal profiles. The former events are of particular importance for the determination of the energy scale of the Observatory, and the latter are characteristic of unusual air shower physics or exotic primary particle types. The purpose of targeted (or "rapid") monitoring is to improve the resolution of the atmospheric measurements for such events. In this paper, we report on the implementation of the rapid monitoring program and its current status. The rapid monitoring data have been analyzed and applied to the reconstruction of air showers of high interest, and indicate that the air fluorescence measurements affected by clouds and aerosols are effectively corrected using measurements from the regular atmospheric monitoring program. We find that the rapid monitoring program has potential for supporting dedicated physics analyses beyond the standard event reconstruction.

preprint2011arXiv

A contribution of stellar flares to the GRXE -- based on MAXI observations --

Using unbiased observations of MAXI/GSC the potential contribution of stellar flares and CVs to GRXE luminosity is estimated in the energy range of 2 - 10 keV. Currently, a reasonable luminosity has been obtained extrapolating the number of stellar flares and that of CVs toward the Galactic ridge from those of the observed flares and CVs near the solar system. The ionized emission lines of Si to Fe are also simulated making the composite thermal spectrum which is based on MAXI observations of stellar flares and CVs. The present estimated result strongly supports a picture that the cumulative stellar flares contribute primarily to the GRXE in terms of the composite thermal spectrum with emission lines and secondary contribution is due to the thermal spectrum with high temperature from CVs.

preprint2011arXiv

Advanced functionality for radio analysis in the Offline software framework of the Pierre Auger Observatory

The advent of the Auger Engineering Radio Array (AERA) necessitates the development of a powerful framework for the analysis of radio measurements of cosmic ray air showers. As AERA performs "radio-hybrid" measurements of air shower radio emission in coincidence with the surface particle detectors and fluorescence telescopes of the Pierre Auger Observatory, the radio analysis functionality had to be incorporated in the existing hybrid analysis solutions for fluoresence and surface detector data. This goal has been achieved in a natural way by extending the existing Auger Offline software framework with radio functionality. In this article, we lay out the design, highlights and features of the radio extension implemented in the Auger Offline framework. Its functionality has achieved a high degree of sophistication and offers advanced features such as vectorial reconstruction of the electric field, advanced signal processing algorithms, a transparent and efficient handling of FFTs, a very detailed simulation of detector effects, and the read-in of multiple data formats including data from various radio simulation codes. The source code of this radio functionality can be made available to interested parties on request.

preprint2011arXiv

Anisotropy and chemical composition of ultra-high energy cosmic rays using arrival directions measured by the Pierre Auger Observatory

The Pierre Auger Collaboration has reported evidence for anisotropy in the distribution of arrival directions of the cosmic rays with energies $E>E_{th}=5.5\times 10^{19}$ eV. These show a correlation with the distribution of nearby extragalactic objects, including an apparent excess around the direction of Centaurus A. If the particles responsible for these excesses at $E>E_{th}$ are heavy nuclei with charge $Z$, the proton component of the sources should lead to excesses in the same regions at energies $E/Z$. We here report the lack of anisotropies in these directions at energies above $E_{th}/Z$ (for illustrative values of $Z=6,\ 13,\ 26$). If the anisotropies above $E_{th}$ are due to nuclei with charge $Z$, and under reasonable assumptions about the acceleration process, these observations imply stringent constraints on the allowed proton fraction at the lower energies.

preprint2011arXiv

Energy decay and frequency shift of a superconducting qubit from non-equilibrium quasiparticles

Quasiparticles are an important decoherence mechanism in superconducting qubits, and can be described with a complex admittance that is a generalization of the Mattis-Bardeen theory. By injecting non-equilibrium quasiparticles with a tunnel junction, we verify qualitatively the expected change of the decay rate and frequency in a phase qubit. With their relative change in agreement to within 4% of prediction, the theory can be reliably used to infer quasiparticle density. We describe how settling of the decay rate may allow determination of whether qubit energy relaxation is limited by non-equilibrium quasiparticles.

preprint2011arXiv

Finite-temperature phase diagram of two-component bosons in a cubic optical lattice: Three-dimensional t-J model of hard-core bosons

We study the three-dimensional bosonic t-J model, i.e., the t-J model of "bosonic electrons", at finite temperatures. This model describes the $s={1 \over 2}$ Heisenberg spin model with the anisotropic exchange coupling $J_{\bot}=-αJ_z$ and doped {\it bosonic} holes, which is an effective system of the Bose-Hubbard model with strong repulsions. The bosonic "electron" operator $B_{rσ}$ at the site $r$ with a two-component (pseudo-)spin $σ(=1,2)$ is treated as a hard-core boson operator, and represented by a composite of two slave particles; a "spinon" described by a Schwinger boson (CP$^1$ boson) $z_{rσ}$ and a "holon" described by a hard-core-boson field $ϕ_r$ as $B_{rσ}=ϕ^†_r z_{rσ}$. By means of Monte Carlo simulations, we study its finite-temperature phase structure including the $α$ dependence, the possible phenomena like appearance of checkerboard long-range order, super-counterflow, superfluid, and phase separation, etc. The obtained results may be taken as predictions about experiments of two-component cold bosonic atoms in the cubic optical lattice.

preprint2011arXiv

Fluctuating pancake vortices revealed by dissipation of Josephson vortex lattice

In strongly anisotropic layered superconductors in tilted magnetic fields the Josephson vortex lattice coexists with the lattice of pancake vortices. Due to the interaction between them, the dissipation of the Josephson-vortex lattice occurs to be very sensitive to the presence of the pancake vortices. If the c-axis magnetic field is smaller then the corresponding lower critical field, the pancake stacks are not formed but the individual pancakes may exist in the fluctuational regime either near surface in large-size samples or in the central region for small-size mesas. We calculate the contribution of such fluctuating pancake vortices to the c-axis conductivity of the Josephson vortex lattice and compare the theoretical results with measurements on small mesas fabricated out of Bi$_{2}$Sr$_{2}$CaCu$_{2}$O$_{8+δ}$ crystals. A fingerprint of fluctuating pancakes is characteristic exponential dependence of the c-axis conductivity observed experimentally. Our results provide strong evidence of the existence of the fluctuating pancakes and their influence on the Josephson-vortex-lattice dissipation.

preprint2011arXiv

Implementing the Quantum von Neumann Architecture with Superconducting Circuits

The von Neumann architecture for a classical computer comprises a central processing unit and a memory holding instructions and data. We demonstrate a quantum central processing unit that exchanges data with a quantum random-access memory integrated on a chip, with instructions stored on a classical computer. We test our quantum machine by executing codes that involve seven quantum elements: Two superconducting qubits coupled through a quantum bus, two quantum memories, and two zeroing registers. Two vital algorithms for quantum computing are demonstrated, the quantum Fourier transform, with 66% process fidelity, and the three-qubit Toffoli OR phase gate, with 98% phase fidelity. Our results, in combination especially with longer qubit coherence, illustrate a potentially viable approach to factoring numbers and implementing simple quantum error correction codes.

preprint2011arXiv

Phase qubits fabricated with trilayer junctions

We have developed a novel Josephson junction geometry with minimal volume of lossy isolation dielectric, being suitable for higher quality trilayer junctions implemented in qubits. The junctions are based on in-situ deposited trilayers with thermal tunnel oxide, have micron-sized areas and a low subgap current. In qubit spectroscopy only a few avoided level crossings are observed, and the measured relaxation time of $T_1\approx400\;\rm{nsec}$ is in good agreement with the usual phase qubit decay time, indicating low loss due to the additional isolation dielectric.

preprint2011arXiv

Photon shell game in three-resonator circuit quantum electrodynamics

The generation and control of quantum states of light constitute fundamental tasks in cavity quantum electrodynamics (QED). The superconducting realization of cavity QED, circuit QED, enables on-chip microwave photonics, where superconducting qubits control and measure individual photon states. A long-standing issue in cavity QED is the coherent transfer of photons between two or more resonators. Here, we use circuit QED to implement a three-resonator architecture on a single chip, where the resonators are interconnected by two superconducting phase qubits. We use this circuit to shuffle one- and two-photon Fock states between the three resonators, and demonstrate qubit-mediated vacuum Rabi swaps between two resonators. This illustrates the potential for using multi-resonator circuits as photon quantum registries and for creating multipartite entanglement between delocalized bosonic modes.

preprint2011arXiv

Search for First Harmonic Modulation in the Right Ascension Distribution of Cosmic Rays Detected at the Pierre Auger Observatory

We present the results of searches for dipolar-type anisotropies in different energy ranges above $2.5\times 10^{17}$ eV with the surface detector array of the Pierre Auger Observatory, reporting on both the phase and the amplitude measurements of the first harmonic modulation in the right-ascension distribution. Upper limits on the amplitudes are obtained, which provide the most stringent bounds at present, being below 2% at 99% $C.L.$ for EeV energies. We also compare our results to those of previous experiments as well as with some theoretical expectations.

preprint2011arXiv

Search for signatures of magnetically-induced alignment in the arrival directions measured by the Pierre Auger Observatory

We present the results of an analysis of data recorded at the Pierre Auger Observatory in which we search for groups of directionally-aligned events (or `multiplets') which exhibit a correlation between arrival direction and the inverse of the energy. These signatures are expected from sets of events coming from the same source after having been deflected by intervening coherent magnetic fields. The observation of several events from the same source would open the possibility to accurately reconstruct the position of the source and also measure the integral of the component of the magnetic field orthogonal to the trajectory of the cosmic rays. We describe the largest multiplets found and compute the probability that they appeared by chance from an isotropic distribution. We find no statistically significant evidence for the presence of multiplets arising from magnetic deflections in the present data.

preprint2011arXiv

Surface spin fluctuations probed with flux noise and coherence in Josephson phase qubits

We measure the dependence of qubit phase coherence and surface spin induced flux noise on inductor loop geometry. While wider inductor traces change neither the flux noise power spectrum nor the qubit dephasing time, increased inductance leads to a simultaneous increase in both. Using our protocol for measuring low frequency flux noise, we make a direct comparison between the flux noise spectrum and qubit phase decay, finding agreement within 10% of theory. The dependence of the measured flux noise on inductor geometry is consistent with a noise source correlation length between 6 and 400 um.

preprint2011arXiv

The effect of the geomagnetic field on cosmic ray energy estimates and large scale anisotropy searches on data from the Pierre Auger Observatory

We present a comprehensive study of the influence of the geomagnetic field on the energy estimation of extensive air showers with a zenith angle smaller than $60^\circ$, detected at the Pierre Auger Observatory. The geomagnetic field induces an azimuthal modulation of the estimated energy of cosmic rays up to the ~2% level at large zenith angles. We present a method to account for this modulation of the reconstructed energy. We analyse the effect of the modulation on large scale anisotropy searches in the arrival direction distributions of cosmic rays. At a given energy, the geomagnetic effect is shown to induce a pseudo-dipolar pattern at the percent level in the declination distribution that needs to be accounted for.

preprint2011arXiv

The Lateral Trigger Probability function for the Ultra-High Energy Cosmic Ray Showers detected by the Pierre Auger Observatory

In this paper we introduce the concept of Lateral Trigger Probability (LTP) function, i.e., the probability for an extensive air shower (EAS) to trigger an individual detector of a ground based array as a function of distance to the shower axis, taking into account energy, mass and direction of the primary cosmic ray. We apply this concept to the surface array of the Pierre Auger Observatory consisting of a 1.5 km spaced grid of about 1600 water Cherenkov stations. Using Monte Carlo simulations of ultra-high energy showers the LTP functions are derived for energies in the range between 10^{17} and 10^{19} eV and zenith angles up to 65 degs. A parametrization combining a step function with an exponential is found to reproduce them very well in the considered range of energies and zenith angles. The LTP functions can also be obtained from data using events simultaneously observed by the fluorescence and the surface detector of the Pierre Auger Observatory (hybrid events). We validate the Monte-Carlo results showing how LTP functions from data are in good agreement with simulations.

preprint2011arXiv

Trigger and Aperture of the Surface Detector Array of the Pierre Auger Observatory

The surface detector array of the Pierre Auger Observatory consists of 1600 water-Cherenkov detectors, for the study of extensive air showers (EAS) generated by ultra-high-energy cosmic rays. We describe the trigger hierarchy, from the identification of candidate showers at the level of a single detector, amongst a large background (mainly random single cosmic ray muons), up to the selection of real events and the rejection of random coincidences. Such trigger makes the surface detector array fully efficient for the detection of EAS with energy above $3\times 10^{18}$ eV, for all zenith angles between 0$^\circ$ and 60$^\circ$, independently of the position of the impact point and of the mass of the primary particle. In these range of energies and angles, the exposure of the surface array can be determined purely on the basis of the geometrical acceptance.

preprint2011arXiv

Tunable terahertz emission from Bi2Sr2CaCu2O8+δ mesa devices

We have measured coherent terahertz emission spectra from Bi2Sr2CaCu2O8+δ mesa devices, as a function of temperature and mesa bias voltage. The emission frequency is found to be tunable by up to 12% by varying the temperature and bias voltage. We attribute the appearance of tunability to asymmetric boundaries at the top and bottom, and the non-rectangular cross-section of the mesas. This interpretation is consistent with numerical simulations of the dynamics of intrinsic Josephson junctions in the mesa. Easily tunable emission frequency may have important implications for the design of terahertz devices based on stacked intrinsic Josephson junctions.

preprint2010arXiv

Demonstration of a Tuneable Coupler for Superconducting Qubits Using Coherent, Time Domain, Two-Qubit Operations

A major challenge in the field of quantum computing is the construction of scalable qubit coupling architectures. Here, we demonstrate a novel tuneable coupling circuit that allows superconducting qubits to be coupled over long distances. We show that the inter-qubit coupling strength can be arbitrarily tuned over nanosecond timescales within a sequence that mimics actual use in an algorithm. The coupler has a measured on/off ratio of 1000. The design is self-contained and physically separate from the qubits, allowing the coupler to be used as a module to connect a variety of elements such as qubits, resonators, amplifiers, and readout circuitry over long distances. Such design flexibility is likely to be essential for a scalable quantum computer.

preprint2010arXiv

Deterministic entanglement of photons in two superconducting microwave resonators

Quantum entanglement, one of the defining features of quantum mechanics, has been demonstrated in a variety of nonlinear spin-like systems. Quantum entanglement in linear systems has proven significantly more challenging, as the intrinsic energy level degeneracy associated with linearity makes quantum control more difficult. Here we demonstrate the quantum entanglement of photon states in two independent linear microwave resonators, creating N-photon NOON states as a benchmark demonstration. We use a superconducting quantum circuit that includes Josephson qubits to control and measure the two resonators, and we completely characterize the entangled states with bipartite Wigner tomography. These results demonstrate a significant advance in the quantum control of linear resonators in superconducting circuits.

preprint2010arXiv

Femtosecond core level photoemision spectroscopy on 1T-TaS2 using 60 eV laser source

Time-resolved photoelectron spectroscopy (trPES) can directly detect transient electronic structure, thus bringing out its promising potential to clarify nonequilibrium processes arising in condensed matters. Here we report the result of core-level (CL) trPES on 1T-TaS2, realized by developing a high-intensity 60 eV laser obtained by high-order harmonic (HH) generation. Ta4f CL-trPES offers the transient amplitude of the charge-density-wave (CDW), via the site-selective and real-time observation of Ta electrons. The present result indicates an ultrafast photoinduced melting and recovery of CDW amplitude, followed by a peculiar long-life oscillation (i.e. collective amplitudon excitation) accompanying the transfer of 0.01 electrons among adjacent Ta atoms. CL-trPES offers a broad range of opportunities for investigating the ultrafast atom-specific electron dynamics in photo-related phenomena of interest.

preprint2010arXiv

Formation of Cosmic Crystals in Highly-Supersaturated Silicate Vapor Produced by Planetesimal Bow Shocks

Several lines of evidence suggest that fine silicate crystals observed in primitive meteorite and interplanetary dust particles (IDPs) nucleated in a supersaturated silicate vapor followed by crystalline growth. We investigated evaporation of $μ$m-sized silicate particles heated by a bow shock produced by a planetesimal orbiting in the gas in the early solar nebula and condensation of crystalline silicate from the vapor thus produced. Our numerical simulation of shock-wave heating showed that these μm-sized particles evaporated almost completely when the bow shock is strong enough to cause melting of chondrule precursor dust particles. We found that the silicate vapor cools very rapidly with expansion into the ambient unshocked nebular region; the cooling rate is estimated, for instance, to be as high as 2000 K s$^{-1}$ for a vapor heated by a bow shock associated with a planetesimal of radius 1 km. The rapid cooling of the vapor leads to nonequilibrium gas-phase condensation of dust at temperatures much lower than those expected from the equilibrium condensation. It was found that the condensation temperatures are lower by a few hundred K or more than the equilibrium temperatures. This explains the results of the recent experimental studies of condensation from a silicate vapor that condensation in such large supercooling reproduces morphologies similar to those of silicate crystals found in meteorites. Our results suggest strongly that the planetesimal bow shock is one of the plausible sites for formation of not only chondrules but also other cosmic crystals in the early solar system.

preprint2010arXiv

Generation of Three-Qubit Entangled States using Superconducting Phase Qubits

Entanglement is one of the key resources required for quantum computation, so experimentally creating and measuring entangled states is of crucial importance in the various physical implementations of a quantum computer. In superconducting qubits, two-qubit entangled states have been demonstrated and used to show violations of Bell's Inequality and to implement simple quantum algorithms. Unlike the two-qubit case, however, where all maximally-entangled two-qubit states are equivalent up to local changes of basis, three qubits can be entangled in two fundamentally different ways, typified by the states $|\mathrm{GHZ}> = (|000> + |111>)/\sqrt{2}$ and $|\mathrm{W}> = (|001> + |010> + |100>)/\sqrt{3}$. Here we demonstrate the operation of three coupled superconducting phase qubits and use them to create and measure $|\mathrm{GHZ}>$ and $|\mathrm{W}>$ states. The states are fully characterized using quantum state tomography and are shown to satisfy entanglement witnesses, confirming that they are indeed examples of three-qubit entanglement and are not separable into mixtures of two-qubit entanglement.

preprint2010arXiv

Quantum process tomography of two-qubit controlled-Z and controlled-NOT gates using superconducting phase qubits

We experimentally demonstrate quantum process tomography of controlled-Z and controlled-NOT gates using capacitively-coupled superconducting phase qubits. These gates are realized by using the $|2\rangle$ state of the phase qubit. We obtain a process fidelity of 0.70 for the controlled-phase and 0.56 for the controlled-NOT gate, with the loss of fidelity mostly due to single-qubit decoherence. The controlled-Z gate is also used to demonstrate a two-qubit Deutsch-Jozsa algorithm with a single function query.

preprint2010arXiv

Resonance Fluorescence of a Single Artificial Atom

An atom in open space can be detected by means of resonant absorption and reemission of electromagnetic waves, known as resonance fluorescence, which is a fundamental phenomenon of quantum optics. We report on the observation of scattering of propagating waves by a single artificial atom. The behavior of the artificial atom, a superconducting macroscopic two-level system, is in a quantitative agreement with the predictions of quantum optics for a pointlike scatterer interacting with the electromagnetic field in one-dimensional open space. The strong atom-field interaction as revealed in a high degree of extinction of propagating waves will allow applications of controllable artificial atoms in quantum optics and photonics.

preprint2010arXiv

Update on the correlation of the highest energy cosmic rays with nearby extragalactic matter

Data collected by the Pierre Auger Observatory through 31 August 2007 showed evidence for anisotropy in the arrival directions of cosmic rays above the Greisen-Zatsepin-Kuz'min energy threshold, \nobreak{$6\times 10^{19}$eV}. The anisotropy was measured by the fraction of arrival directions that are less than $3.1^\circ$ from the position of an active galactic nucleus within 75 Mpc (using the Véron-Cetty and Véron $12^{\rm th}$ catalog). An updated measurement of this fraction is reported here using the arrival directions of cosmic rays recorded above the same energy threshold through 31 December 2009. The number of arrival directions has increased from 27 to 69, allowing a more precise measurement. The correlating fraction is $(38^{+7}_{-6})%$, compared with $21%$ expected for isotropic cosmic rays. This is down from the early estimate of $(69^{+11}_{-13})%$. The enlarged set of arrival directions is examined also in relation to other populations of nearby extragalactic objects: galaxies in the 2 Microns All Sky Survey and active galactic nuclei detected in hard X-rays by the Swift Burst Alert Telescope. A celestial region around the position of the radiogalaxy Cen A has the largest excess of arrival directions relative to isotropic expectations. The 2-point autocorrelation function is shown for the enlarged set of arrival directions and compared to the isotropic expectation.

preprint2009arXiv

$μ$SR studies of the frustrated quasi-2d square-lattice spin system Cu(Cl,Br)La(Nb,Ta)$_{2}$O$_{7}$: evolution from spin-gap to antiferromagnetic state

We report muon spin relaxation ($μ$SR) and magnetic susceptibility measurements on Cu(Cl,Br)La(Nb,Ta)$_{2}$O$_{7}$, which demonstrate: (a) the absence of static magnetism in (CuCl)LaNb$_{2}$O$_{7}$ down to 15 mK confirming a spin-gapped ground state; (b) phase separation between partial volumes with a spin-gap and static magnetism in (CuCl)La(Nb,Ta)$_{2}$O$_{7}$; (c) history-dependent magnetization in the (Nb,Ta) and (Cl,Br) substitution systems; (d) a uniform long-range collinear antiferromagnetic state in (CuBr)LaNb$_{2}$O$_{7}$; and (e) a decrease of Néel temperature with decreasing Br concentration $x$ in Cu(Cl$_{1-x}$Br$_{x}$)LaNb$_{2}$O$_{7}$ with no change in the ordered Cu moment size for $0.33 \leq x \leq 1$. Together with several other $μ$SR studies of quantum phase transitions in geometrically-frustrated spin systems, the present results reveal that the evolution from a spin-gap to a magnetically ordered state is often associated with phase separation and/or a first order phase transition.

preprint2009arXiv

Emission of Terahertz Waves from Stacks of Intrinsic Josephson Junctions

By patterning mesoscopic crystals of Bi2Sr2CaCu2O8 (BSCCO) into electromagnetic resonators the oscillations of a large number of intrinsic Josephson junctions can be synchronized into a macroscopic coherent state accompanied by the emission of strong continuous wave THz-radiation. The temperature dependence of the emission is governed by the interplay of self-heating in the resonator and by re-trapping of intrinsic Josephson junctions which can yield a strongly nonmonotonic temperature dependence of the emission power. Furthermore, proper shaping of the resonators yields THzsources with voltage-tunable emission frequencies.

preprint2009arXiv

Scaling behavior of the crossover to short-stack regimes of Josephson vortex lattices in Bi$_2$Sr$_2$CaCu$_2$O$_{8+δ}$ stacks

We report the systematic investigations of the oscillation of the Josephson vortex (JV) flow resistance in Bi$_2$Sr$_2$CaCu$_2$O$_{8+δ}$ micro-fabricated junctions with various geometries and superconducting anisotropy parameters. As the applied magnetic field parallel to the $ab$-plane is increased, oscillation with a period corresponding to a $ϕ_0/2$ par atomic Josephson junction changes to oscillation with a doubled period. This crossover is scaled by both the junction length and the anisotropy parameter, indicating that the bulk inductive coupling that favors the triangular JV lattice is replaced with the surface deformation energy as the dominant interaction for a JV lattice. These results suggest that the in-phase square JV lattice is pronounced at a higher magnetic field in a smaller and more anisotropic sample.

preprint2009arXiv

The Fluorescence Detector of the Pierre Auger Observatory

The Pierre Auger Observatory is a hybrid detector for ultra-high energy cosmic rays. It combines a surface array to measure secondary particles at ground level together with a fluorescence detector to measure the development of air showers in the atmosphere above the array. The fluorescence detector comprises 24 large telescopes specialized for measuring the nitrogen fluorescence caused by charged particles of cosmic ray air showers. In this paper we describe the components of the fluorescence detector including its optical system, the design of the camera, the electronics, and the systems for relative and absolute calibration. We also discuss the operation and the monitoring of the detector. Finally, we evaluate the detector performance and precision of shower reconstructions.

preprint2009arXiv

Thermal Management in Large Bi2212 Mesas used for Terahertz Sources

We present a thermal analysis of a patterned mesa on a Bi2Sr2CaCu2O8 (Bi2212) single crystal that is based on tunneling characteristics of the c-axis stack of ~800 intrinsic Josephson junctions in the mesa. Despite the large mesa volume (e.g., 40x300x1.2 micron3) and power dissipation that result in selfheating and backbending of the current-voltage curve (I-V), there are accessible bias conditions for which significant polarized THz-wave emission can be observed. We estimate the mesa temperature by equating the quasiparticle resistance, Rqp(T), to the ratio V/I over the entire I-V including the backbending region. These temperatures are used to predict the unpolarized black-body radiation reaching our bolometer and there is substantial agreement over the entire I-V. As such, backbending results from the particular Rqp(T) for Bi2212, as first discussed by Fenton, rather than a significant suppression of the energy gap. This model also correctly predicts the observed disappearance of backbending above ~60 K.

preprint2006arXiv

Cosmogenic Neutrinos from the propagation of Ultra High Energy Nuclei

We calculate the flux of neutrinos generated by the propagation of ultra-high energy nuclei over cosmological distances. The propagation takes into account the interactions with cosmic background radiations including the CMB and the most recent estimates of higher energy (infra-red, optical, and ultra violet) backgrounds. We assume that the composition of ultra-high energy cosmic rays (UHECRs) at the source is the same as the observed one at low energies. This assumption fits well the present data at the highest energies. We compare the cosmogenic neutrino flux from mixed composition sources to pure proton sources. We find that the neutrino flux in the mixed composition case has a high energy peak, mainly due to photopion production off CMB photons, of similar shape and amplitude to the proton case. At low energies both composition cases have significant neutrino flux with a peak around 1014.5 eV due to the higher energy backgrounds. The mixed composition case induces a higher flux of neutrinos at energies below 1013 eV due to the neutron decay component that extends down to low energies. Detection of diffuse neutrino fluxes at ultra high energies can strongly constrain the source distribution of UHECR whereas fluxes at lower energies could be used to constrain confinement of VHE and UHE cosmic rays if combined with composition analysis from cosmic ray experiments.

preprint2006arXiv

Near Infrared Polarization Images of the Orion Nebula

Wide-field (~8' x 8') and deep near-infrared (JHKs bands) polarization images of the Orion nebulae (IRNe) around young stellar objects (YSOs), both massive and low-mass. We found the IRNe around both IRc2 and BN to be very extensive, suggesting that there might be two extended (>0.7 pc) bipolar/monopolar IRNe in these sources. We discovered at least 13 smaller-scale (~0.01-0.1 pc) IRNe around less-massive YSOs including the famous source theta^2 Ori C. We also suggest the presence of many unresolved (<690 AU) systems around low-mass YSOs and young brown dwarfs showing possible intrinsic polarizations. Wide-field infrared polarimetry is thus demonstrated to be a powerful technique in revealing IRNe and hence potential disk/outflow systems among high-mass to substellar YSOs.

preprint2006arXiv

The First Scientific Results from the Pierre Auger Observatory

The southern site of the Pierre Auger Observatory is under the construction near Malargue in Argentina and now more than 60% of the detectors are completed. The observatory has been collecting data for over 1 year and the cumulative exposure is already similar to that of the largest forerunner experiments. The hybrid technique provides model-independent energy measurements from the Fluorescence Detector to calibrate the Surface Detector. Based on this technique, the first estimation of the energy spectrum above 3 EeV has been presented and is discussed in this paper.

preprint2002arXiv

Quantum oscillations in two coupled charge qubits

Despite an apparent progress in implementing individual solid-state qubits, there have been no experimental reports so far on multi-bit gates required for building a real quantum computer. Here we report a new circuit comprising two coupled charge qubits. Using a pulse technique, we coherently mix quantum states and observe quantum oscillations whose spectrum reflects interaction between the qubits. Our results demonstrate the feasibility of coupling of multiple solid-state qubits and indicate the existence of entangled two-qubit states.

preprint2000arXiv

Small-scale anisotropy of cosmic rays above 10^19eV observed with the Akeno Giant Air Shower Array

With the Akeno Giant Air Shower Array (AGASA), 581 cosmic rays above 10^19eV, 47 above 4 x 10^19eV, and 7 above 10^20eV are observed until August 1998. Arrival direction distribution of these extremely high energy cosmic rays has been studied. While no significant large-scale anisotropy is found on the celestial sphere, some interesting clusters of cosmic rays are observed. Above 4 x 10^19eV, there are one triplet and three doublets within separation angle of 2.5^o and the probability of observing these clusters by a chance coincidence under an isotropic distribution is smaller than 1 %. Especially the triplet is observed against expected 0.05 events. The cos(θ_GC) distribution expected from the Dark Matter Halo model fits the data as well as an isotropic distribution above 2 x 10^19eV and 4 x 10^19eV, but is a poorer fit than isotropy above 10^19eV. Arrival direction distribution of seven 10^20eV cosmic rays is consistent with that of lower energy cosmic rays and is uniform. Three of seven are members of doublets above about 4 x 10^19eV.

preprint1999arXiv

The Anisotropy of Cosmic Ray Arrival Directions around 10$^{18}$eV

Anisotropy in the arrival directions of cosmic rays with energies above 10$^{17}$eV is studied using data from the Akeno 20 km$^2$ array and the Akeno Giant Air Shower Array (AGASA), using a total of about 117,000 showers observed during 11 years. In the first harmonic analysis, we have found strong anisotropy of $\sim$ 4% around 10$^{18}$eV, corresponding to a chance probability of 0.2%. With two dimensional analysis in right ascension and declination, this anisotropy is interpreted as an excess of showers near the directions of the Galactic Center and the Cygnus region.