Source author record

S. Nagataki

S. Nagataki 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

53works
5topics
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

53 published item(s)

preprint2022arXiv

Effects of Nuclear Equation of State on Type-I X-ray Bursts: Interpretation of the X-ray Bursts from GS 1826-24

Type I X-ray bursts are thermonuclear explosions on the neutron star (NS) surface caused by mass accretion from a companion star. Observations of X-ray bursts provide valuable information on X-ray binary systems, e.g., binary parameters, the chemical composition of accreted matter, and the nuclear equation of state (EOS) of NSs. There have been several theoretical studies to constrain the physics of X-ray bursters. However, they have mainly focused on the burning layers above the solid crust of the NS, which brings up issues of the treatment of NS gravitation and internal energy. In this study, focusing on the microphysics inside NSs, we calculate a series of X-ray bursts using a general-relativistic stellar-evolution code with several NS EOSs. We compare the X-ray-burst models with the burst parameters of a clocked burster associated with GS 1826-24. We find a monotonic correlation between the NS radius and the light-curve profile. A larger radius shows a higher recurrence time and a large peak luminosity. In contrast, the dependence of light curves on the NS mass becomes more complicated, where the neutrino cooling suppresses the efficiency of nuclear ignition. We also constrain the EOS and mass of GS 1826-24, i.e., stiffer EOSs, corresponding to larger NS radii, are not preferred due to too-high peak luminosity. The EOS and the cooling and heating of NSs are important to discuss the theoretical and observational properties of X-ray bursts.

preprint2022arXiv

Search for Spatial Correlations of Neutrinos with Ultra-High-Energy Cosmic Rays

For several decades, the origin of ultra-high-energy cosmic rays (UHECRs) has been an unsolved question of high-energy astrophysics. One approach for solving this puzzle is to correlate UHECRs with high-energy neutrinos, since neutrinos are a direct probe of hadronic interactions of cosmic rays and are not deflected by magnetic fields. In this paper, we present three different approaches for correlating the arrival directions of neutrinos with the arrival directions of UHECRs. The neutrino data is provided by the IceCube Neutrino Observatory and ANTARES, while the UHECR data with energies above $\sim$50 EeV is provided by the Pierre Auger Observatory and the Telescope Array. All experiments provide increased statistics and improved reconstructions with respect to our previous results reported in 2015. The first analysis uses a high-statistics neutrino sample optimized for point-source searches to search for excesses of neutrinos clustering in the vicinity of UHECR directions. The second analysis searches for an excess of UHECRs in the direction of the highest-energy neutrinos. The third analysis searches for an excess of pairs of UHECRs and highest-energy neutrinos on different angular scales. None of the analyses has found a significant excess, and previously reported over-fluctuations are reduced in significance. Based on these results, we further constrain the neutrino flux spatially correlated with UHECRs.

preprint2021arXiv

Observation of Variations in Cosmic Ray Single Count Rates During Thunderstorms and Implications for Large-Scale Electric Field Changes

We present the first observation by the Telescope Array Surface Detector (TASD) of the effect of thunderstorms on the development of cosmic ray single count rate intensity over a 700 km$^{2}$ area. Observations of variations in the secondary low-energy cosmic ray counting rate, using the TASD, allow us to study the electric field inside thunderstorms, on a large scale, as it progresses on top of the 700 km$^{2}$ detector, without dealing with the limitation of narrow exposure in time and space using balloons and aircraft detectors. In this work, variations in the cosmic ray intensity (single count rate) using the TASD, were studied and found to be on average at the $\sim(0.5-1)\%$ and up to 2\% level. These observations were found to be both in excess and in deficit. They were also found to be correlated with lightning in addition to thunderstorms. These variations lasted for tens of minutes; their footprint on the ground ranged from 6 to 24 km in diameter and moved in the same direction as the thunderstorm. With the use of simple electric field models inside the cloud and between cloud to ground, the observed variations in the cosmic ray single count rate were recreated using CORSIKA simulations. Depending on the electric field model used and the direction of the electric field in that model, the electric field magnitude that reproduces the observed low-energy cosmic ray single count rate variations was found to be approximately between 0.2-0.4 GV. This in turn allows us to get a reasonable insight on the electric field and its effect on cosmic ray air showers inside thunderstorms.

preprint2020arXiv

Evidence for a Supergalactic Structure of Magnetic Deflection Multiplets of Ultra-High Energy Cosmic Rays

Evidence for a large-scale supergalactic cosmic ray multiplet (arrival directions correlated with energy) structure is reported for ultra-high energy cosmic ray (UHECR) energies above 10$^{19}$ eV using seven years of data from the Telescope Array (TA) surface detector and updated to 10 years. Previous energy-position correlation studies have made assumptions regarding magnetic field shapes and strength, and UHECR composition. Here the assumption tested is that, since the supergalactic plane is a fit to the average matter density of the local Large Scale Structure (LSS), UHECR sources and intervening extragalactic magnetic fields are correlated with this plane. This supergalactic deflection hypothesis is tested by the entire field-of-view (FOV) behavior of the strength of intermediate-scale energy-angle correlations. These multiplets are measured in spherical cap section bins (wedges) of the FOV to account for coherent and random magnetic fields. The structure found is consistent with supergalactic deflection, the previously published energy spectrum anisotropy results of TA (the hotspot and coldspot), and toy-model simulations of a supergalactic magnetic sheet. The seven year data post-trial significance of this supergalactic structure of multiplets appearing by chance, on an isotropic sky, is found by Monte Carlo simulation to be 4.2$σ$. The ten years of data post-trial significance is 4.1$σ$. Furthermore, the starburst galaxy M82 is shown to be a possible source of the TA Hotspot, and an estimate of the supergalactic magnetic field using UHECR measurements is presented.

preprint2020arXiv

Measurement of the Proton-Air Cross Section with Telescope Array's Black Rock Mesa and Long Ridge Fluorescence Detectors, and Surface Array in Hybrid Mode

Ultra high energy cosmic rays provide the highest known energy source in the universe to measure proton cross sections. Though conditions for collecting such data are less controlled than an accelerator environment, current generation cosmic ray observatories have large enough exposures to collect significant statistics for a reliable measurement for energies above what can be attained in the lab. Cosmic ray measurements of cross section use atmospheric calorimetry to measure depth of air shower maximum ($X_{\mathrm{max}}$), which is related to the primary particle's energy and mass. The tail of the $X_{\mathrm{max}}$ distribution is assumed to be dominated by showers generated by protons, allowing measurement of the inelastic proton-air cross section. In this work the proton-air inelastic cross section measurement, $σ^{\mathrm{inel}}_{\mathrm{p-air}}$, using data observed by Telescope Array's Black Rock Mesa and Long Ridge fluorescence detectors and surface detector array in hybrid mode is presented. $σ^{\mathrm{inel}}_{\mathrm{p-air}}$ is observed to be $520.1 \pm 35.8$[Stat.] $^{+25.0}_{-40}$[Sys.]~mb at $\sqrt{s} = 73$ TeV. The total proton-proton cross section is subsequently inferred from Glauber formalism and is found to be $σ^{\mathrm{tot}}_{\mathrm{pp}} = 139.4 ^{+23.4}_{-21.3}$ [Stat.]$ ^{+15.0}_{-24.0}$[Sys.]~mb.

preprint2020arXiv

Observations of the Origin of Downward Terrestrial Gamma-Ray Flashes

In this paper we report the first close, high-resolution observations of downward-directed terrestrial gamma-ray flashes (TGFs) detected by the large-area Telescope Array cosmic ray observatory, obtained in conjunction with broadband VHF interferometer and fast electric field change measurements of the parent discharge. The results show that the TGFs occur during strong initial breakdown pulses (IBPs) in the first few milliseconds of negative cloud-to-ground and low-altitude intracloud flashes, and that the IBPs are produced by a newly-identified streamer-based discharge process called fast negative breakdown. The observations indicate the relativistic runaway electron avalanches (RREAs) responsible for producing the TGFs are initiated by embedded spark-like transient conducting events (TCEs) within the fast streamer system, and potentially also by individual fast streamers themselves. The TCEs are inferred to be the cause of impulsive sub-pulses that are characteristic features of classic IBP sferics. Additional development of the avalanches would be facilitated by the enhanced electric field ahead of the advancing front of the fast negative breakdown. In addition to showing the nature of IBPs and their enigmatic sub-pulses, the observations also provide a possible explanation for the unsolved question of how the streamer to leader transition occurs during the initial negative breakdown, namely as a result of strong currents flowing in the final stage of successive IBPs, extending backward through both the IBP itself and the negative streamer breakdown preceding the IBP.

preprint2020arXiv

Search for Large-scale Anisotropy on Arrival Directions of Ultra-high-energy Cosmic Rays Observed with the Telescope Array Experiment

Motivated by the detection of a significant dipole structure in the arrival directions of ultrahigh-energy cosmic rays above 8 EeV reported by the Pierre Auger Observatory (Auger), we search for a large-scale anisotropy using data collected with the surface detector array of the Telescope Array Experiment (TA). With 11 years of TA data, a dipole structure in a projection of the right ascension is fitted with an amplitude of 3.3+- 1.9% and a phase of 131 +- 33 degrees. The corresponding 99% confidence-level upper limit on the amplitude is 7.3%. At the current level of statistics, the fitted result is compatible with both an isotropic distribution and the dipole structure reported by Auger.

preprint2020arXiv

Search for point sources of ultra-high energy photons with the Telescope Array surface detector

The surface detector (SD) of the Telescope Array (TA) experiment allows one to indirectly detect photons with energies of order $10^{18}$ eV and higher and to separate photons from the cosmic-ray background. In this paper we present the results of a blind search for point sources of ultra-high energy (UHE) photons in the Northern sky using the TA SD data. The photon-induced extensive air showers (EAS) are separated from the hadron-induced EAS background by means of a multivariate classifier based upon 16 parameters that characterize the air shower events. No significant evidence for the photon point sources is found. The upper limits are set on the flux of photons from each particular direction in the sky within the TA field of view, according to the experiment's angular resolution for photons. Average 95% C.L. upper limits for the point-source flux of photons with energies greater than $10^{18}$, $10^{18.5}$, $10^{19}$, $10^{19.5}$ and $10^{20}$ eV are $0.094$, $0.029$, $0.010$, $0.0073$ and $0.0058$ km$^{-2}$yr$^{-1}$ respectively. For the energies higher than $10^{18.5}$ eV, the photon point-source limits are set for the first time. Numerical results for each given direction in each energy range are provided as a supplement to this paper.

preprint2020arXiv

The fully developed remnant of a neutrino-driven supernova: Evolution of ejecta structure and asymmetries in SNR Cassiopeia A

Abridged. We aim at exploring to which extent the remnant keeps memory of the asymmetries that develop stochastically in the neutrino-heating layer due to hydrodynamic instabilities (e.g., convective overturn and the standing accretion shock instability) during the first second after core bounce. We coupled a 3D HD model of a neutrino-driven SN explosion with 3D MHD/HD simulations of the remnant formation. The simulations cover 2000 years of expansion and include all physical processes relevant to describe the complexities in the SN evolution and the subsequent interaction of the stellar debris with the wind of the progenitor star. The interaction of large-scale asymmetries left from the earliest phases of the explosion with the reverse shock produces, at the age of $\approx 350$~years, an ejecta structure and a remnant morphology which are remarkably similar to those observed in Cas A. Small-scale structures in the large-scale Fe-rich plumes created during the initial stages of the SN, combined with HD instabilities that develop after the passage of the reverse shock, naturally produce a pattern of ring- and crown-like structures of shocked ejecta. The consequence is a spatial inversion of the ejecta layers with Si-rich ejecta being physically interior to Fe-rich ejecta. The full-fledged remnant shows voids and cavities in the innermost unshocked ejecta resulting from the expansion of Fe-rich plumes and their inflation due to the decay of radioactive species. The asymmetric distributions of $^{44}$Ti and $^{56}$Fe and their abundance ratio are both compatible with those inferred from high-energy observations of Chandra and NuSTAR. The main asymmetries observed in the ejecta distribution of Cas A can be explained by the interaction of the reverse shock with the large-scale asymmetries that developed from stochastic processes that originate during the first seconds of the SN blast.

preprint2019arXiv

Hydrodynamic simulations unravel the progenitor-supernova-remnant connection in SN 1987A

(Abridged) We aim at linking the dynamical and radiative properties of the remnant of SN 1987A to the geometrical and physical characteristics of the parent aspherical SN explosion and to the internal structure of its progenitor star. We performed 3D hydrodynamic simulations which describe the long-term evolution of SN 1987A from the onset of the SN to the full-fledged remnant at the age of 50 years, accounting for the pre-SN structure of the progenitor star. The simulations include all physical processes relevant for the complex phases of SN evolution and for the interaction of the SNR with the highly inhomogeneous ambient environment around SN 1987A. From the simulations, we synthesize observables to be compared with observations. By comparing the model results with observations, we constrained the initial SN anisotropy causing Doppler shifts observed in emission lines of heavy elements from ejecta, and leading to the remnant evolution observed in the X-ray band in the last 30 years. In particular, we found that the high mixing of ejecta unveiled by high redshifts and broadenings of [FeII] and $^{44}$Ti lines require a highly asymmetric SN explosion channeling a significant fraction of energy along an axis almost lying in the plane of the central equatorial ring around SN 1987A, roughly along the line-of-sight but with an offset of 40 deg, with the lobe propagating away from the observer slightly more energetic than the other. We found unambiguously that the observed distribution of ejecta and the dynamical and radiative properties of the SNR can be best reproduced if the structure of the progenitor star was that of a blue supergiant resulted from the merging of two massive stars.

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.

preprint2016arXiv

First Upper Limits on the Radar Cross Section of Cosmic-Ray Induced Extensive Air Showers

TARA (Telescope Array Radar) is a cosmic ray radar detection experiment colocated with Telescope Array, the conventional surface scintillation detector (SD) and fluorescence telescope detector (FD) near Delta, Utah, U.S.A. The TARA detector combines a 40 kW, 54.1 MHz VHF transmitter and high-gain transmitting antenna which broadcasts the radar carrier over the SD array and within the FD field of view, towards a 250 MS/s DAQ receiver. TARA has been collecting data since 2013 with the primary goal of observing the radar signatures of extensive air showers (EAS). Simulations indicate that echoes are expected to be short in duration (~10 microseconds) and exhibit rapidly changing frequency, with rates on the order of 1 MHz/microsecond. The EAS radar cross-section (RCS) is currently unknown although it is the subject of over 70 years of speculation. A novel signal search technique is described in which the expected radar echo of a particular air shower is used as a matched filter template and compared to waveforms obtained by triggering the radar DAQ using the Telescope Array fluorescence detector. No evidence for the scattering of radio frequency radiation by EAS is obtained to date. We report the first quantitative RCS upper limits using EAS that triggered the Telescope Array Fluorescence Detector.

preprint2016arXiv

Search for correlations between the arrival directions of IceCube neutrino events and ultrahigh-energy cosmic rays detected by the Pierre Auger Observatory and the Telescope Array

This paper presents the results of different searches for correlations between very high-energy neutrino candidates detected by IceCube and the highest-energy cosmic rays measured by the Pierre Auger Observatory and the Telescope Array. We first consider samples of cascade neutrino events and of high-energy neutrino-induced muon tracks, which provided evidence for a neutrino flux of astrophysical origin, and study their cross-correlation with the ultrahigh-energy cosmic ray (UHECR) samples as a function of angular separation. We also study their possible directional correlations using a likelihood method stacking the neutrino arrival directions and adopting different assumptions on the size of the UHECR magnetic deflections. Finally, we perform another likelihood analysis stacking the UHECR directions and using a sample of through-going muon tracks optimized for neutrino point-source searches with sub-degree angular resolution. No indications of correlations at discovery level are obtained for any of the searches performed. The smallest of the p-values comes from the search for correlation between UHECRs with IceCube high-energy cascades, a result that should continue to be monitored.

preprint2015arXiv

A Northern Sky Survey for Point-Like Sources of EeV Neutral Particles with the Telescope Array Experiment

We report on the search for steady point-like sources of neutral particles around 10$^{18}$ eV between 2008 May and 2013 May with the scintillator surface detector of the Telescope Array experiment. We found overall no significant point-like excess above 0.5 EeV in the northern sky. Subsequently, we also searched for coincidence with the Fermi bright Galactic sources. No significant coincidence was found within the statistical uncertainty. Hence, we set an upper limit on the neutron flux that corresponds to an averaged flux of 0.07 km$^{-2}$ yr$^{-1}$ for $E>1$ EeV in the northern sky at the 95% confidence level. This is the most stringent flux upper limit in a northern sky survey assuming point-like sources. The upper limit at the 95% confidence level on the neutron flux from Cygnus X-3 is also set to 0.2 km$^{-2}$ yr$^{-1}$ for $E>0.5$ EeV. This is an order of magnitude lower than previous flux measurements.

preprint2015arXiv

Luminosity--time and luminosity--luminosity correlations for GRB prompt and afterglow plateau emissions

We present an analysis of 123 Gamma-ray bursts (GRBs) with known redshifts possessing an afterglow plateau phase. We reveal that $L_a-T^{*}_a$ correlation between the X-ray luminosity $L_a$ at the end of the plateau phase and the plateau duration, $T^*_a$, in the GRB rest frame has a power law slope different, within more than 2 $σ$, from the slope of the prompt $L_{f}-T^{*}_{f}$ correlation between the isotropic pulse peak luminosity, $L_{f}$, and the pulse duration, $T^{*}_{f}$, from the time since the GRB ejection. Analogously, we show differences between the prompt and plateau phases in the energy-duration distributions with the afterglow emitted energy being on average $10\%$ of the prompt emission. Moreover, the distribution of prompt pulse versus afterglow spectral indexes do not show any correlation. In the further analysis we demonstrate that the $L_{peak}-L_a$ distribution, where $L_{peak}$ is the peak luminosity from the start of the burst, is characterized with a considerably higher Spearman correlation coefficient, $ρ=0.79$, than the one involving the averaged prompt luminosity, $L_{prompt}-L_a$, for the same GRB sample, yielding $ρ=0.60$. Since some of this correlation could result from the redshift dependences of the luminosities, namely from their cosmological evolution we use the Efron-Petrosian method to reveal the intrinsic nature of this correlation. We find that a substantial part of the correlation is intrinsic. We apply a partial correlation coefficient to the new de-evolved luminosities showing that the intrinsic correlation exists.

preprint2015arXiv

Matter Mixing in Core-collapse Supernova Ejecta: Large Density Perturbations in the Progenitor Star?

Matter mixing is one important topic in the study of core-collapse supernova (CCSN) explosions. In this paper, we perform two-dimensional hydrodynamic simulations to reproduce the high velocity $^{56}$Ni clumps observed in SN 1987A. This is the first time that large density perturbation is proposed in the CCSN progenitor to generate Rayleigh-Taylor (RT) instability and make the effective matter mixing. In the case of a spherical explosion, RT instability is efficient at both C+O/He and He/H interfaces of the SN progenitor. Radial coherent structures shown in perturbation patterns are important for obtaining high velocity $^{56}$Ni clumps. We can also obtain matter mixing features and high velocity $^{56}$Ni clumps in some cases of aspherical explosion. We find that one of the most favorable models in our work has a combination of bipolar and equatorially asymmetric explosions in which at least 25\% of density perturbation is introduced at different composition interfaces of the CCSN progenitor. These simulation results are comparable to the observational findings of SN 1987A.

preprint2015arXiv

Measurement of the Proton-Air Cross Section with Telescope Array's Middle Drum Detector and Surface Array in Hybrid Mode

In this work we are reporting on the measurement of the proton-air inelastic cross section $σ^{\rm inel}_{\rm p-air}$ using the Telescope Array (TA) detector. Based on the measurement of the $σ^{\rm inel}_{\rm p-air}$ the proton-proton cross section $σ_{\rm p-p}$ value is also determined at $\sqrt{s} = 95_{-8}^{+5}$ TeV. Detecting cosmic ray events at ultra high energies with Telescope Array enables us to study this fundamental parameter that we are otherwise unable to access with particle accelerators. The data used in this report is the hybrid events observed by the Middle Drum fluorescence detector together with the surface array detector collected over five years. The value of the $σ^{\rm inel}_{\rm p-air}$ is found to be equal to $567.0 \pm 70.5 [{\rm Stat.}] ^{+29}_{-25} [{\rm Sys.}]$ mb. The total proton-proton cross section is subsequently inferred from Glauber Formalism and Block, Halzen and Stanev QCD inspired fit and is found to be equal to $170_{-44}^{+48} [{\rm Stat.}] _{-17}^{+19} [{\rm Sys.}] $mb.

preprint2015arXiv

The IceCube Neutrino Observatory, the Pierre Auger Observatory and the Telescope Array: Joint Contribution to the 34th International Cosmic Ray Conference (ICRC 2015)

We have conducted three searches for correlations between ultra-high energy cosmic rays detected by the Telescope Array and the Pierre Auger Observatory, and high-energy neutrino candidate events from IceCube. Two cross-correlation analyses with UHECRs are done: one with 39 cascades from the IceCube `high-energy starting events' sample and the other with 16 high-energy `track events'. The angular separation between the arrival directions of neutrinos and UHECRs is scanned over. The same events are also used in a separate search using a maximum likelihood approach, after the neutrino arrival directions are stacked. To estimate the significance we assume UHECR magnetic deflections to be inversely proportional to their energy, with values $3^\circ$, $6^\circ$ and $9^\circ$ at 100 EeV to allow for the uncertainties on the magnetic field strength and UHECR charge. A similar analysis is performed on stacked UHECR arrival directions and the IceCube sample of through-going muon track events which were optimized for neutrino point-source searches.

preprint2014arXiv

A CR-hydro-NEI Model of the Structure and Broadband Emission from Tycho's SNR

Tycho's supernova remnant (SNR) is well-established as a source of particle acceleration to very high energies. Constraints from numerous studies indicate that the observed gamma-ray emission results primarily from hadronic processes, providing direct evidence of highly relativistic ions that have been accelerated by the SNR. Here we present an investigation of the dynamical and spectral evolution of Tycho's SNR by carrying out hydrodynamical simulations that include diffusive shock acceleration of particles in the amplified magnetic field at the forward shock of the SNR. Our simulations provide a consistent view of the shock positions, the nonthermal emission, the thermal X-ray emission from the forward shock, and the brightness profiles of the radio and X-ray emission. We compare these with the observed properties of Tycho to determine the density of the ambient material, the particle acceleration efficiency and maximum energy, the accelerated electron to-proton ratio, and the properties of the shocked gas downstream of the expanding SNR shell. We find that evolution of a typical Type Ia supernova in a low ambient density (n_0 ~ 0.3 cm^{-3}), with an upstream magnetic field of ~5\ microGauss, and with ~16% of the SNR kinetic energy being converted into relativistic electrons and ions through diffusive shock acceleration, reproduces the observed properties of Tycho. Under such a scenario, the bulk of observed gamma-ray emission at high energies is produced by pi^0-decay resulting from the collisions of energetic hadrons, while inverse-Compton emission is significant at lower energies, comprising roughly half of the flux between 1 and 10 GeV.

preprint2014arXiv

Indications of Intermediate-Scale Anisotropy of Cosmic Rays with Energy Greater Than 57 EeV in the Northern Sky Measured with the Surface Detector of the Telescope Array Experiment

We have searched for intermediate-scale anisotropy in the arrival directions of ultrahigh-energy cosmic rays with energies above 57~EeV in the northern sky using data collected over a 5 year period by the surface detector of the Telescope Array experiment. We report on a cluster of events that we call the hotspot, found by oversampling using 20$^\circ$-radius circles. The hotspot has a Li-Ma statistical significance of 5.1$σ$, and is centered at R.A.=146.7$^{\circ}$, Dec.=43.2$^{\circ}$. The position of the hotspot is about 19$^{\circ}$ off of the supergalactic plane. The probability of a cluster of events of 5.1$σ$ significance, appearing by chance in an isotropic cosmic-ray sky, is estimated to be 3.7$\times$10$^{-4}$ (3.4$σ$).

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

Study of Ultra-High Energy Cosmic Ray Composition Using Telescope Array's Middle Drum Detector and Surface Array in Hybrid Mode

Previous measurements of the composition of Ultra-High Energy Cosmic Rays(UHECRs) made by the High Resolution Fly's Eye(HiRes) and Pierre Auger Observatory(PAO) are seemingly contradictory, but utilize different detection methods, as HiRes was a stereo detector and PAO is a hybrid detector. The five year Telescope Array(TA) Middle Drum hybrid composition measurement is similar in some, but not all, respects in methodology to PAO, and good agreement is evident between data and a light, largely protonic, composition when comparing the measurements to predictions obtained with the QGSJetII-03 and QGSJet-01c models. These models are also in agreement with previous HiRes stereo measurements, confirming the equivalence of the stereo and hybrid methods. The data is incompatible with a pure iron composition, for all models examined, over the available range of energies. The elongation rate and mean values of Xmax are in good agreement with Pierre Auger Observatory data. This analysis is presented using two methods: data cuts using simple geometrical variables and a new pattern recognition technique.

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

Energy Spectrum of Ultra-High Energy Cosmic Rays Observed with the Telescope Array Using a Hybrid Technique

We measure the spectrum of cosmic rays with energies greater than $10^{18.2}$ eV with the Fluorescence Detectors (FDs) and the Surface Detectors (SDs) of the Telescope Array Experiment using the data taken in our first 2.3-year observation from May 27 2008 to September 7 2010. A hybrid air shower reconstruction technique is employed to improve accuracies in determination of arrival directions and primary energies of cosmic rays using both FD and SD data. The energy spectrum presented here is in agreement with our previously published spectra and the HiRes results.

preprint2013arXiv

Search for correlations of the arrival directions of ultra-high energy cosmic ray with extragalactic objects as observed by the telescope array experiment

We search for correlations between positions of extragalactic objects and arrival directions of Ultra-High Energy Cosmic Rays (UHECRs) with primary energy $E \ge 40$ EeV as observed by the surface detector array of the Telescope Array (TA) experiment during the first 40 months of operation. We examined several public astronomical object catalogs, including the Veron-Cetty and Veron catalog of active galactic nuclei. We counted the number of TA events correlated with objects in each catalog as a function of three parameters: the maximum angular separation between a TA event and an object, the minimum energy of the events, and the maximum redshift of the objects. We determine combinations of these parameters which maximize the correlations, and calculate the chance probabilities of having the same levels of correlations from an isotropic distribution of UHECR arrival directions. No statistically significant correlations are found when penalties for scanning over the above parameters and for searching in several catalogs are taken into account.

preprint2013arXiv

The Cosmic Ray Energy Spectrum Observed with the Surface Detector of the Telescope Array Experiment

The Telescope Array (TA) collaboration has measured the energy spectrum of ultra-high energy cosmic rays with primary energies above 1.6 x 10^(18) eV. This measurement is based upon four years of observation by the surface detector component of TA. The spectrum shows a dip at an energy of 4.6 x 10^(18) eV and a steepening at 5.4 x 10^(19) eV which is consistent with the expectation from the GZK cutoff. We present the results of a technique, new to the analysis of ultra-high energy cosmic ray surface detector data, that involves generating a complete simulation of ultra-high energy cosmic rays striking the TA surface detector. The procedure starts with shower simulations using the CORSIKA Monte Carlo program where we have solved the problems caused by use of the "thinning" approximation. This simulation method allows us to make an accurate calculation of the acceptance of the detector for the energies concerned.

preprint2013arXiv

The Energy Spectrum of Ultra-High-Energy Cosmic Rays Measured by the Telescope Array FADC Fluorescence Detectors in Monocular Mode

We present a measurement of the energy spectrum of ultra-high-energy cosmic rays performed by the Telescope Array experiment using monocular observations from its two new FADC-based fluorescence detectors. After a short description of the experiment, we describe the data analysis and event reconstruction procedures. Since the aperture of the experiment must be calculated by Monte Carlo simulation, we describe this calculation and the comparisons of simulated and real data used to verify the validity of the aperture calculation. Finally, we present the energy spectrum calculated from the merged monocular data sets of the two FADC-based detectors, and also the combination of this merged spectrum with an independent, previously published monocular spectrum measurement performed by Telescope Array's third fluorescence detector (Abu-Zayyad {\it et al.}, {Astropart. Phys.} 39 (2012), 109). This combined spectrum corroborates the recently published Telescope Array surface detector spectrum (Abu-Zayyad {\it et al.}, {Astrophys. Journ.} 768 (2013), L1) with independent systematic uncertainties.

preprint2013arXiv

Upper limit on the flux of photons with energies above 10^19 eV using the Telescope Array surface detector

We search for ultra-high energy photons by analyzing geometrical properties of shower fronts of events registered by the Telescope Array surface detector. By making use of an event-by-event statistical method, we derive upper limits on the absolute flux of primary photons with energies above 10^19, 10^19.5 and 10^20 eV based on the first three years of data taken.

preprint2012arXiv

Modeling high-energy light curves of the PSR B1259-63/LS 2883 binary based on 3-D SPH simulations

Temporal changes of X-ray to very-high-energy gamma-ray emissions from the pulsar-Be star binary PSR B1259-63/LS 2883 are studied based on 3-D SPH simulations of pulsar wind interaction with Be-disk and wind. We focus on the periastron passage of the binary and calculate the variation of the synchrotron and inverse-Compton emissions using the simulated shock geometry and pressure distribution of the pulsar wind. The characteristic double-peaked X-ray light curve from observations is reproduced by our simulation under a dense Be disk condition (base density ~10^{-9} g cm^{-3}). We interpret the pre- and post-periastron peaks as being due to a significant increase in the conversion efficiency from pulsar spin down power to the shock-accelerated particle energy at orbital phases when the pulsar crosses the disk before periastron passage, and when the pulsar wind creates a cavity in the disk gas after periastron passage, respectively. On the contrary, in the model TeV light curve, which also shows a double peak feature, the first peak appears around the periastron phase. The possible effects of cooling processes on the TeV light curve are briefly discussed.

preprint2012arXiv

Search for Anisotropy of Ultra-High Energy Cosmic Rays with the Telescope Array Experiment

We study the anisotropy of Ultra-High Energy Cosmic Ray (UHECR) events collected by the Telescope Array (TA) detector in the first 40 months of operation. Following earlier studies, we examine event sets with energy thresholds of 10 EeV, 40 EeV, and 57 EeV. We find that the distributions of the events in right ascension and declination are compatible with an isotropic distribution in all three sets. We then compare with previously reported clustering of the UHECR events at small angular scales. No significant clustering is found in the TA data. We then check the events with E>57 EeV for correlations with nearby active galactic nuclei. No significant correlation is found. Finally, we examine all three sets for correlations with the large-scale structure of the Universe. We find that the two higher-energy sets are compatible with both an isotropic distribution and the hypothesis that UHECR sources follow the matter distribution of the Universe (the LSS hypothesis), while the event set with E>10 EeV is compatible with isotropy and is not compatible with the LSS hypothesis at 95% CL unless large deflection angles are also assumed. We show that accounting for UHECR deflections in a realistic model of the Galactic magnetic field can make this set compatible with the LSS hypothesis.

preprint2012arXiv

The Energy Spectrum of Telescope Array's Middle Drum Detector and the Direct Comparison to the High Resolution Fly's Eye Experiment

The Telescope Array's Middle Drum fluorescence detector was instrumented with telescopes refurbished from the High Resolution Fly's Eye's HiRes-1 site. The data observed by Middle Drum in monocular mode was analyzed via the HiRes-1 profile-constrained geometry reconstruction technique and utilized the same calibration techniques enabling a direct comparison of the energy spectra and energy scales between the two experiments. The spectrum measured using the Middle Drum telescopes is based on a three-year exposure collected between December 16, 2007 and December 16, 2010. The calculated difference between the spectrum of the Middle Drum observations and the published spectrum obtained by the data collected by the HiRes-1 site allows the HiRes-1 energy scale to be transferred to Middle Drum. The HiRes energy scale is applied to the entire Telescope Array by making a comparison between Middle Drum monocular events and hybrid events that triggered both Middle Drum and the Telescope Array's scintillator Ground Array.

preprint2011arXiv

A high-frequency study of the Sunyaev-Zel'dovich effect morphology in galaxy clusters

High-frequency, high-resolution imaging of the Sunyaev-Zel'dovich (SZ) effect is an important technique to study the complex structures of the atmospheres of merging galaxy clusters. Such observations are sensitive to the details of the electron spectrum. We show that the morphology of the SZ intensity maps in simulated galaxy clusters observed at 345 GHz, 600 GHz, and 857 GHz are significantly different because of SZ relativistic corrections. These differences can be revealed by high-resolution imaging instruments. We calculate relativistically corrected SZ intensity maps of a simulated, massive, merging galaxy cluster and of the massive, merging clusters 1E0657-558 (the Bullet Cluster) and Abell 2219. The morphologies of the SZ intensity maps are remarkably different between 345 GHz and 857 GHz for each merging cluster. We show that high-resolution imaging observations of the SZ intensity maps at these frequencies, obtainable with the LABOCA and HERSCHEL-SPIRE instruments, allow to fully exploit the astrophysical relevance of the predicted SZ morphological effect.

preprint2011arXiv

An Off-Axis Relativistic Jet Model for the Type Ic supernova SN 2007gr

We propose an off-axis relativistic jet model for the Type Ic supernova SN 2007gr. Most of the energy ($\sim2\times10^{51}$ erg) in the explosion is contained in non-relativistic ejecta which produces the supernova. The optical emission is coming from the decay process of $\rm ^{56}Ni$ synthesized in the bulk SN ejecta. Only very little energy ($\sim10^{48}$ erg) is contained in the relativistic jet with initial velocity about 0.94 times the speed of light. The radio and X-ray emission comes from this relativistic jet. With some typical parameters of a Wolf-Rayet star (progenitor of Type Ic SN), i.e., the mass loss rate $\dot{M}=1.0 \times10^{-5} M_{\odot} \rm yr^{-1}$ and the wind velocity $v_{\rm w}=1.5\times10^{3} \rm km s^{-1}$ together with an observing angle of $θ_{\rm obs} = 63.3^{\circ}$, we can obtain the multiband light curves that fit the observations well. All the observed data are consistent with our model. Thus we conclude that SN 2007gr contains a weak relativistic jet and we are observing the jet from off-axis.

preprint2011arXiv

Can electron distribution functions be derived through the Sunyaev-Zel'dovich effect?

Measurements of the Sunyaev-Zel'dovich (hereafter SZ) effect distortion of the cosmic microwave background provide methods to derive the gas pressure and temperature of galaxy clusters. Here we study the ability of SZ effect observations to derive the electron distribution function (DF) in massive galaxy clusters. Our calculations of the SZ effect include relativistic corrections considered within the framework of the Wright formalism and use a decomposition technique of electron DFs into Fourier series. Using multi-frequency measurements of the SZ effect, we find the solution of a linear system of equations that is used to derive the Fourier coefficients; we further analyze different frequency samples to decrease uncertainties in Fourier coefficient estimations. We propose a method to derive DFs of electrons using SZ multi-frequency observations of massive galaxy clusters. We found that the best frequency sample to derive an electron DF includes high frequencies $ν$=375, 600, 700, 857 GHz. We show that it is possible to distinguish a Juttner DF from a Maxwell-Bolzman DF as well as from a Juttner DF with the second electron population by means of SZ observations for the best frequency sample if the precision of SZ intensity measurements is less than 0.1%. We demonstrate by means of 3D hydrodynamic numerical simulations of a hot merging galaxy cluster that the morphologies of SZ intensity maps are different for frequencies $ν$=375, 600, 700, 857 GHz. We stress that measurements of SZ intensities at these frequencies are a promising tool for studying electron distribution functions in galaxy clusters.

preprint2011arXiv

Failed Gamma-Ray Bursts: Thermal UV/Soft X-ray Emission Accompanied by Peculiar Afterglows

We show that the photospheres of "failed" Gamma-Ray Bursts (GRBs), whose bulk Lorentz factors are much lower than 100, can be outside of internal shocks. The resulting radiation from the photospheres is thermal and bright in UV/Soft X-ray band. The photospheric emission lasts for about one thousand seconds with luminosity about several times 10^46 erg/s. These events can be observed by current and future satellites. It is also shown that the afterglows of failed GRBs are peculiar at the early stage, which makes it possible to distinguish failed GRBs from ordinary GRBs and beaming-induced orphan afterglows.

preprint2011arXiv

Unveiling the 3D temperature structure of galaxy clusters by means of the thermal SZ effect

The Sunyaev-Zel'dovich (hereafter SZ) effect is a promising tool to derive the gas temperature of galaxy clusters. The approximation of a spherically symmetric gas distribution is usually used to determine the temperature structure of galaxy clusters, but this approximation cannot properly describe merging galaxy clusters. The methods used so far, which do not assume the spherically symmetric distribution, permit us to derive 2D temperature maps of merging galaxy clusters. In this paper, we propose a method to derive the standard temperature deviation and temperature variance along the line-of-sight, which permits us to analyze the 3D temperature structure of galaxy clusters by means of the thermal SZ effect. We also propose a method to reveal merger shock waves in galaxy clusters by analyzing the presence of temperature inhomogeneities along the line-of-sight.

preprint2010arXiv

An analysis of the temperature structure of galaxy clusters by means of the thermal Sunyaev-Zel'dovich effect

Measurements of the Sunyaev-Zel'dovich (hereafter SZ) effect distortion of the cosmic microwave background provide us with an independent method to derive the gas temperature of galaxy clusters. In merging galaxy clusters the gas distribution is inhomogeneous and, therefore, the method of temperature measuring based on the SZ effect should be more relevant than that based on an X-ray emission analysis. Here we study a method for measuring the gas temperature in merging clusters by means of the SZ effect. Our calculations of intensity maps of the SZ effect include relativistic corrections considered within the framework of the Wright formalism and utilize a cosmological numerical simulation of a merging galaxy cluster evolved with its baryon physics. We found that the gas temperature in merging clusters can be measured by means of the ratio of the SZ intensity at a low frequency (128 GHz) to that at a high frequency (369 GHz). This SZ intensity ratio permits us to reveal prominent features of the temperature structure caused by violent merger shock waves. Therefore, measurements of the ratio of the SZ intensities are a promising tool for measuring gas temperature in merging galaxy clusters.

preprint2009arXiv

Propagation of ultrahigh energy nuclei in clusters of galaxies: resulting composition and secondary emissions

We study the survival of ultrahigh energy nuclei injected in clusters of galaxies, as well as their secondary neutrino and photon emissions, using a complete numerical propagation method and a realistic modeling of the magnetic, baryonic and photonic backgrounds. It is found that the survival of heavy nuclei highly depends on the injection position and on the profile of the magnetic field. Taking into account the limited lifetime of the central source could also lead in some cases to the detection of a cosmic ray afterglow, temporally decorrelated from neutrino and gamma ray emissions. We calculate that the diffusive neutrino flux around 1 PeV coming from clusters of galaxies may have a chance to be detected by current instruments. The observation of single sources in neutrinos and in gamma rays produced by ultrahigh energy cosmic rays will be more difficult. Signals coming from lower energy cosmic rays (E < 1 PeV), if they exist, might however be detected by Fermi, for reasonable sets of parameters.

preprint2008arXiv

High Energy neutrino signals from the Epoch of Reionization

We perform a new estimate of the high energy neutrinos expected from GRBs associated with the first generation of stars in light of new models and constraints on the epoch of reionization and a more detailed evaluation of the neutrino emission yields. We also compare the diffuse high energy neutrino background from Population III stars with the one from "ordinary stars" (Population II), as estimated consistently within the same cosmological and astrophysical assumptions. In disagreement with previous literature, we find that high energy neutrinos from Population III stars will not be observable with current or near future neutrino telescopes, falling below both IceCube sensitivity and atmospheric neutrino background under the most extreme assumptions for the GRB rate. This rules them out as a viable diagnostic tool for these still elusive metal-free stars.

preprint2002arXiv

Revision of the Selection Function of the Optical Redshift Survey using the Sloan Digital Sky Survey Early Data Release: Toward an Accurate Estimate of Source Number Density of Ultra-High Energy Cosmic Rays

If Ultra-High Energy Cosmic Rays (UHECRs) are originated from nearby galaxies, modeling of the distribution of nearby galaxies is important to an accurate estimate the source number density of UHECRs. We investigate uncertainty of the selection function of the Optical Redshift Survey (ORS), which we used to construct a source model of UHECRs. The investigation is based on a comparison of numbe counts of ORS galaxies with those of the spectroscopic sample of the Sloan Digital Sky Survey (SDSS) Early Data Release (EDR). We carefully count galaxies in the same absolute magnitude bin from the two samples. We find a slight systematic overestimate of the ORS counts outside 5000 km s$^{-1}$ by about a factor of 2. We revise the selection function of the ORS assuming that the SDSS counts are correct. Our revision is based on the absorption given in the ORS catalog as well as that computed from Schlegel et al. (1998), which is systematically larger than the former by $A_B \sim 0.1$ mag in the region of low absorption. It is found that introduction of Schlegel et al.'s absorption changes one of the parameters of the ORS selection function by more than 10%. The revision should be taken into account in the future analysis of the source number density of UHECRs based on the ORS. Using the revised selection function, we determine the global structure of the Local Supercluster (LSC) with a source model of UHECRs, that is, a number-density model consisting of a uniform spherical halo and an exponential disk. We find that the revision is insignificant in terms of the structure of the LSC. However, the revised selection function will be useful to other studies such as peculiar velocity and correlation function.

preprint2001arXiv

Propagation of UHECRs from the Sources in the Super-Galactic Plane

We have performed the detailed numerical simulations on the propagation of the UHE protons in the energy range $E=(10^{19.5} - 10^{22.0}$) eV in the relatively strong extra-galactic magnetic field with strength $B= (10, 100)$ nG within about 40 Mpc. In this case, the deflection angles of UHECRs become so large that the no counterparts problem is simply solved. As for the source distribution, we assumed that it is proportional to the number distribution of galaxies within the GZK sphere. We have found many clusters, which mean the small-scale anisotropy, in our simulations. It has been also shown that the observed energy spectrum is well reproduced in our models without any fine-tuned parameter. We have used the correlation value in order to investigate statistically the similarity between the distribution of arrival directions of UHECRs and that of galaxies. We have found that each correlation value for each parameter set begins to converge when the number of the detected events becomes ${\large O}(10^3)$. Since the expected number counts by the experiment of the next generation such as TA, HiRes, Auger, and EUSO are thought to be the order of $10^3$, we will be able to determine the source distribution and values of the parameters in this study in the very near future. Compared with the AGASA data, the significant anisotropy on the arrival directions of UHECRs are found in the analysis of first and second harmonics. This may originate from the incompleteness of the ORS database. This problem may be also solved if the source distribution is slightly changed. For example, this problem may be solved if we assume that UHECRs come from some of the galaxies such as AGNs and radio galaxies.

preprint2001arXiv

Rapid-Process Nucleosynthesis in Neutrino-Magneto-Centrifugally Driven Winds

We have studied whether the rotation and magnetic fields in neutrino-driven winds can be key processes for the rapid-process (r-process) nucleosynthesis. We have examined the features of a steady and subsonic wind solutions which extend the model of Weber and Davis (1967), which is a representative solar wind model. As a result, we found that the entropy per baryon becomes lower and the dynamical timescale becomes longer as the angular velocity becomes higher. These results are inappropriate for the production of the r-process nuclei. As for the effects of magnetic fields, we found that a solution as a steady wind from the surface of the proto-neutron star can not be obtained when the strength of the magnetic field becomes $\ge$ $10^{11}$ G. Since the magnetic field in normal pulsars is of order $10^{12}$ G, a steady wind solution might not be realized there, which means that the models in this study may not be adopted for normal proto-neutron stars. In this situation, we have little choice but to conclude that it is difficult to realize a successful r-process nucleosynthesis in the wind models in this framework.

preprint2001arXiv

Specifying the Environments around GRB, Explaining the Fe line in the X-Ray Afterglow of GRB000214

We present a model explaining the Fe K alpha line and the continuum in the afterglow of GRB000214. We pose the importance to seek the physically natural environment around GRB000214. For the reproduction of the observation, we need the ring-like remnant around the progenitor like that of SN 1987A produced by the mass-loss of the progenitor and the fireball spread over in every directions. The observation of GRB000214, in which the continuum power law spectrum decreased faster than the line, motivated us to consider the two independent systems for the line emission and the continuum spectrum. At first, the continuum spectrum can be fitted by the afterglow emission of the fireball pointing toward the observer which does not collide with the ring because the emission of GRB and the afterglow are highly collimated to the observer by the relativistic beaming effect. Secondly, the line can be fitted by the fluorescence of the Fe atoms in the ring illuminated by the X-ray afterglow. The significance of this study is that our model may constrain strongly the GRB model. Although the Supranova model assumes the extreme-ring-like remnant produced by the usual supernova explosion, this may not be probable. It is because the supernova remnants are known to be shell-like. The model also assumes two steps of explosions, on the other hand, we need only one explosion of the progenitor. In this sense, our scenario is more natural. Moreover, in the numerical simulations of Hypernova, the jet of the opening angle of only 1 degree is generated. In our model, the fireball which spreads over in every directions reconciles with the observation of 1 percent of the polarization in the observation of SN1998bw which showed the explosion might not be so collimated.

preprint2000arXiv

Chemical Evolution in the Large Magellanic Cloud

We present a new input parameter set of the Pagel model (Pagel & Tautvai$\rm \breve{s}$ien$\rm \dot{e}$ 1998) for the Large Magellanic Cloud (LMC) in order to reproduce the observations, including the star formation rate (SFR) history. It is concluded that the probability for (3-8)$M_{\odot}$ stars to explode as SNe Ia has to be quite high ($\sim 0.17$) in the LMC. As a result, a steep initial mass function (IMF) slope and existence of the outflow are not needed in order to attain the low [O/Fe] ratio in the LMC. As for the current supernova ratio, a high ratio ($\sim 1.3$) is concluded by the new parameter set, which is consistent with the recent X-ray observations.

preprint2000arXiv

Effects of Jet-like Explosion in SN 1987A

We study the effects of jet-like explosion in SN 1987A. Calculations of the explosive nucleosynthesis and the matter mixing in a jet-like explosion are performed and their results are compared with the observations of SN 1987A. It is shown that the jet-like explosion model is favored because the radioactive nuclei 44Ti is produced in a sufficient amount to explain the observed bolometric luminosity at 3600 days after the explosion. It is also shown that the observed line profiles of Fe[II] are well reproduced by the jet-like explosion model. In particular, the fast moving component travelling at (3000-4000) km/s is well reproduced, which has not been reproduced by the spherical explosion models. Moreover, we conclude that the favored degree of a jet-like explosion to explain the tail of the light curve is consistent with the one favored in the calculation of the matter mixing. The concluded ratio of the velocity along to the polar axis relative to that in the equatorial plane at the Si/Fe interface is about 2:1. This conclusion will give good constraints on the calculations of the dynamics of the collapse-driven supernova. The asymmetry of the observed line profiles of Fe[II] can be explained when the assumption of the equatorial symmetry of the system is removed, which can be caused by the asymmetry of the jet-like explosion with respect to the equatorial plane. In the case of SN 1987A, the jet on the north pole has to be stronger than that on the south pole in order to reproduce the observed asymmetric line profiles. Such an asymmetry may also be the origin of the pulsar kick. When we believe some theories that cause such an asymmetric explosion, the proto -neutron star born in SN 1987A is moving in the southern part of the remnant.

preprint1998arXiv

Dispersion of Inferred SNe Ia/SNe II Ratios for Current Models of Supernova Nucleosynthesis

We have estimated the dispersion of the inferred relative frequencies of Type-Ia and Type- II supernovae ($N_{Ia}/ N_{II}$) in our Galaxy by fitting the numerical results of supernova nucleosynthesis to the solar-system abundances. The ratio $N_{Ia}/ N_{II}$ is estimated to be 0.056--0.14, which is consistent with the observation, if the model of Woosley and Weaver (1995, WW95) is adopted for Type-II supernovae (SNe II). On the other hand, the upper limit of $N_{Ia}/ N_{II}$ becomes too large for the model of Hashimoto (1995, Ha95). However, Ha95 can fit the solar values better than WW95 as far as the abundant nuclei are concerned. These results mean that Ha95 can reproduce the main nuclei of the solar-system well and that WW95 can reproduce the solar values over a wide mass-number range. We also note that $N_{Ia}/ N_{II}$ tends to become smaller if the delayed detonation model is adopted for Type-Ia supernovae (SNe Ia). The dispersion of $N_{Ia}/ N_{II}$ obtained in the present investigation is larger than that concluded by Tsujimoto et al. (1995), which will have an influence on the estimate of the average life time of SNe Ia's progenitors and/or the star-formation rate history in our galaxy.

preprint1998arXiv

Initial Shock Waves for Explosive Nucleosynthesis in Type II Supernova

We have performed 1-dimensional calculations for explosive nucleosynthesis in collapse-driven supernova and investigated its sensitivity to the initial form of the shock wave. We have found the tendency that the peak temperature becomes higher around the mass cut if the input energy is injected more in the form of kinetic energy rather than internal energy. Then, the mass cut becomes larger, and, as a result, neutron-rich matter is less included in the ejecta; this is favorable for producing the observational data compared with a previous model. Our results imply that the standard method to treat various processes for stellar evolution, such as convection and electron capture during the silicon burning stage, are still compatible with the calculation of explosive nucleosynthesis.