BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Shigehiro Nagataki

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

ResearcherUnclaimed source record
astro-ph.HEastro-phastro-ph.COhep-phastro-ph.SRastro-ph.GAgr-qc

Catalog footprint

What is connected

56works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

56 published item(s)

preprint2024arXiv

Correlations and Distinguishability Challenges in Supernova Models: Insights from Future Neutrino Detectors

This paper explores core-collapse supernovae as crucial targets for neutrino telescopes, addressing uncertainties in their simulation results. We comprehensively analyze eighteen modern simulations and discriminate among supernova models using realistic detectors and interactions. A significant correlation between the total neutrino energy and cumulative counts, driven by massive lepton neutrinos and oscillations, is identified, particularly noticeable with the DUNE detector. Bayesian techniques indicate strong potential for model differentiation during a Galactic supernova event, with HK excelling in distinguishing models based on equation of state, progenitor mass, and mixing scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

The impact of effective matter mixing based on three-dimensional hydrodynamical models on the molecule formation in the ejecta of SN 1987A

To investigate the impact of matter mixing on the formation of molecules in the ejecta of SN 1987A, time-dependent rate equations for chemical reactions are solved for one-zone and one-dimensional ejecta models of SN 1987A. The latter models are based on the one-dimensional profiles obtained by angle-averaging of the three-dimensional hydrodynamical models (Ono et al. 2020), which effectively reflect the 3D matter mixing; the impact is demonstrated, for the first time, based on three-dimensional hydrodynamical models. The distributions of initial seed atoms and radioactive $^{56}$Ni influenced by the mixing could affect the formation of molecules. By comparing the calculations for spherical cases and for several specified directions in the bipolar-like explosions in the three-dimensional hydrodynamical models, the impact is discussed. The decay of $^{56}$Ni, practically $^{56}$Co at later phases, could heat the gas and delay the molecule formation. Additionally, Compton electrons produced by the decay could ionize atoms and molecules and could destruct molecules. Several chemical reactions involved with ions such as H$^+$ and He$^+$ could also destruct molecules. The mixing of $^{56}$Ni plays a non-negligible role in both the formation and destruction of molecules through the processes above. The destructive processes of carbon monoxide and silicon monoxide due to the decay of $^{56}$Ni generally reduce the amounts. However, if the molecule formation is sufficiently delayed under a certain condition, the decay of $^{56}$Ni could locally increase the amounts through a sequence of reactions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Additional evidence for a pulsar wind nebula in the heart of SN 1987A from multi-epoch X-ray data and MHD modeling

Since the day of its explosion, supernova (SN) 1987A has been closely monitored to study its evolution and to detect its central compact relic. In fact, the formation of a neutron star is strongly supported by the detection of neutrinos from the SN. However, besides the detection in the Atacama Large Millimeter/submillimeter Array (ALMA) data of a feature that is compatible with the emission arising from a proto-pulsar wind nebula (PWN), the only hint for the existence of such elusive compact object is provided by the detection of hard emission in NuSTAR data up to ~ 20 keV. We report on the simultaneous analysis of multi-epoch observations of SN 1987A performed with Chandra, XMM-Newton and NuSTAR. We also compare the observations with a state-of-the-art 3D magnetohydrodynamic (MHD) simulation of SN 1987A. A heavily absorbed power-law, consistent with the emission from a PWN embedded in the heart of SN 1987A, is needed to properly describe the high-energy part of the observed spectra. The spectral parameters of the best-fit power-law are in agreement with the previous estimate, and exclude diffusive shock acceleration as a possible mechanism responsible for the observed non-thermal emission. The information extracted from our analysis are used to infer the physical characteristics of the pulsar and the broad-band emission of its nebula, in agreement with the ALMA data. Analysis of the synthetic spectra also show that, in the near future, the main contribution to Fe K emission line will originate in the outermost shocked ejecta of SN 1987A.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Can one-zone hadronuclear model explain the hard-TeV spectrum of BL Lac objects?

Context. The intrinsic TeV emission of some BL Lacs are characterized by a hard spectrum (the hard-TeV spectrum) after correcting for the extragalactic background light. The hard-TeV spectra pose a challenge to conventional one-zone models, including the leptonic model, the photohadronic model, the proton synchrotron model, etc. Aims. In this work, we study if the one-zone hadronuclear (pp) model can be used to interpret the hard-TeV spectra of BL Lacs without introducing extreme parameters. Methods. We give analytical calculations to study if there is a parameter space and the charge neutrality condition of jet can be satisfied when interpreting the hard-TeV spectra of BL Lacs without introducing a super-Eddington jet power. Results. We find that in a sample of hard-TeV BL Lacs collected by Xue et al. (2019a), only the hard-TeV spectrum of 1ES 0229+200 could be explained by gamma-ray from pi-0 decay produced in the pp interactions, but at the cost of setting a small radius of the radiation region that comparable to the Schwarzschild radius of the central black hole. Combining with previous studies of other one-zone models, we suggest that the hard-TeV spectra of BL Lacs cannot be explained by any one-zone models without introducing extreme parameters, and should originate from the multiple radiation regions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Optical and X-ray GRB Fundamental Planes as Cosmological Distance Indicators

Gamma-Ray Bursts (GRBs), can be employed as standardized candles, extending the distance ladder beyond Supernovae Type Ia (SNe Ia, $z=2.26$). We standardize GRBs using the 3D fundamental plane relation (the Dainotti relation) among the rest-frame end time of the X-ray plateau emission, its corresponding luminosity, and the peak prompt luminosity. Combining SNe Ia and GRBs, we constrain $Ω_{\text{M}}= 0.299 \pm 0.009$ assuming a flat $Λ$CDM cosmology with and without correcting GRBs for selection biases and redshift evolution. Using a 3D optical Dainotti correlation, we find this sample is as efficacious in the determination of $Ω_{\text{M}}$ as the X-ray sample. We trimmed our GRB samples to achieve tighter planes to simulate additional GRBs. We determined how many GRBs are needed as standalone probes to achieve a comparable precision on $Ω_{\text{M}}$ to the one obtained by SNe Ia only. We reach the same error measurements derived using SNe Ia in 2011 and 2014 with 142 and 284 simulated optical GRBs, respectively, considering the errorbars on the variables halved. These error limits will be reached in 2038 and in 2047, respectively. Using a doubled sample (obtained by future machine learning approaches allowing a lightcurve reconstruction and the estimates of GRB redhifts when z is unknown) compared to the current sample, with errorbars halved we will reach the same precision as SNe Ia in 2011 and 2014, now and in 2026, respectively. If we consider the current SNe precision, this will be reached with 390 optical GRBs by 2054.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

R-process viable outflows are suppressed in global alpha-viscosity models of collapsar disks

Collapsar disks have been proposed to be rich factories of heavy elements, but the major question of whether their outflows are neutron-rich, and could therefore represent significant sites of the rapid neutron-capture (r-) process, or dominated by iron-group elements remains unresolved. We present the first global models of collapsars that start from a stellar progenitor and self-consistently describe the evolution of the disk, its composition, and its outflows in response to the imploding stellar mantle, using energy-dependent M1 neutrino transport and an alpha-viscosity to approximate turbulent angular-momentum transport. We find that a neutron-rich, neutrino-dominated accretion flow (NDAF) is established only marginally--either for short times or relatively low viscosities--because the disk tends to disintegrate into an advective disk (ADAF) already at relatively high mass-accretion rates, launching powerful outflows but preventing it from developing a hot, dense, and therefore neutron-rich core. Viscous outflows disrupt the star within ~100s with explosion energies close to that of hypernovae. If viscosity is neglected, a stable NDAF with disk mass of about 1Msun is formed but is unable to release neutron-rich ejecta, while it produces a relatively mild explosion powered by a neutrino-driven wind blown off its surface. With ejecta electron fractions close to 0.5, all models presumably produce large amounts of Ni56. Our results suggest that collapsar models based on the alpha-viscosity are inefficient r-process sites and that genuinely magnetohydrodynamic effects may be required to generate neutron-rich outflows. A relatively weak effective viscosity generated by magnetohydrodynamic turbulence would improve the prospects for obtaining neutron-rich ejecta.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The double detonation of a double degenerate system, from Type Ia supernova explosion to its supernova remnant

Type Ia supernovae (SNe) are believed to be caused by the thermonuclear explosion of a white dwarf (WD), but the nature of the progenitor system(s) is still unclear. Recent theoretical and observational developments have led to renewed interest in double degenerate models, in particular the "helium-ignited violent merger" or "dynamically-driven double-degenerate double-detonation" (D$^6$). In this paper we take the output of an existing D$^6$ SN model and carry it into the supernova remnant (SNR) phase up to 4000 years after the explosion, past the time when all the ejecta have been shocked. Assuming a uniform ambient medium, we reveal specific signatures of the explosion mechanism and spatial variations intrinsic to the ejecta. The first detonation produces an ejecta tail visible at early times, while the second detonation leaves a central density peak in the ejecta that is visible at late times. The SNR shell is off-centre at all times, because of an initial velocity shift due to binary motion. The companion WD produces a large conical shadow in the ejecta, visible in projection as a dark patch surrounded by a bright ring. This is a clear and long-lasting feature that is localized, and its impact on the observed morphology is dependent on the viewing angle of the SNR. These results offer a new way to diagnose the explosion mechanism and progenitor system using observations of a Type Ia SNR.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A semi-analytic afterglow with thermal electrons and synchrotron self-Compton emission

We extend previous work on gamma-ray burst (GRB) afterglows involving hot thermal electrons at the base of a shock-accelerated tail. Using a physically-motivated electron distribution based on first-principles simulations, we compute broadband emission from radio to TeV gamma-rays. For the first time, we present the effects of a thermal distribution of electrons on synchrotron self-Compton (SSC) emission. The presence of thermal electrons causes temporal and spectral structure across the entire observable afterglow, which is substantively different from models that assume a pure power-law distribution for the electrons. We show that early-time TeV emission is enhanced by more than an order of magnitude for our fiducial parameters, with a time-varying spectral index that does not occur for a pure power law of electrons. We further show that the X-ray "closure relations" take a very different, also time-dependent, form when thermal electrons are present; the shape traced out by the X-ray afterglows is a qualitative match to observations of the traditional decay phase.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Indication of a Pulsar Wind Nebula in the hard X-ray emission from SN 1987A

Since the day of its explosion, SN 1987A (SN87A) was closely monitored with the aim to study its evolution and to detect its central compact relic. The detection of neutrinos from the supernova strongly supports the formation of a neutron star (NS). However, the constant and fruitless search for this object has led to different hypotheses on its nature. Up to date, the detection in the ALMA data of a feature somehow compatible with the emission arising from a proto Pulsar Wind Nebula (PWN) is the only hint of the existence of such elusive compact object. Here we tackle this 33-years old issue by analyzing archived observations of SN87A performed Chandra and NuSTAR in different years. We firmly detect nonthermal emission in the $10-20$ kev energy band, due to synchrotron radiation. The possible physical mechanism powering such emission is twofold: diffusive shock acceleration (DSA) or emission arising from an absorbed PWN. By relating a state-of-the-art magneto-hydrodynamic simulation of SN87A to the actual data, we reconstruct the absorption pattern of the PWN embedded in the remnant and surrounded by cold ejecta. We found that, even though the DSA scenario cannot be firmly excluded, the most likely scenario that well explains the data is the PWN emission.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Sub-Solar Metallicity Progenitor for Cassiopeia A, the remnant of a Type IIb Supernova

We report, for the first time, the detection of the Mn-K$α$ line in the Type IIb supernova (SN IIb) remnant, Cassiopeia A. Manganese ($^{55}$Mn after decay of $^{55}$Co), a neutron-rich element, together with chromium ($^{52}$Cr after decay of $^{52}$Fe), is mainly synthesized at the explosive incomplete Si burning regime. Therefore, the Mn/Cr mass ratio with its neutron excess reflects the neutronization at the relevant burning layer during the explosion. Chandra's archival X-ray data of Cassiopeia A indicate a low Mn/Cr mass ratio with values in the range 0.10--0.66, which, when compared to one-dimensional SN explosion models, requires that the electron fraction be 0.4990 $\lesssim Y_{\rm e} \lesssim$ 0.5 at the incomplete Si burning layer. An explosion model assuming a solar-metallicity progenitor with a typical explosion energy ($1 \times 10^{51}$ erg) fails to reproduce such a high electron fraction. In such models, the explosive Si-burning regime extends only to the Si/O layer established during the progenitor's hydrostatic evolution; the $Y_e$ in the Si/O layer is lower than the value required by our observational constraints. We can satisfy the observed Mn/Cr mass ratio if the explosive Si-burning regime were to extend into the O/Ne hydrostatic layer, which has a higher $Y_{\rm e}$. This would require an energetic ($> 2 \times 10^{51}$ erg) and/or asymmetric explosion of a sub-solar metallicity progenitor ($Z \lesssim 0.5Z_{\odot}$) for Cassiopeia A. The low initial metallicity can be used to rule out a single-star progenitor, leaving the possibility of a binary progenitor with a compact companion (white dwarf, neutron star or black hole). We discuss the detectability of X-rays from Bondi accretion onto such a compact companion around the explosion site. We also discuss other possible mass-loss scenarios for the progenitor system of Cassiopeia A.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Follow-up Observations for IceCube-170922A: Detection of Rapid Near-Infrared Variability and Intensive Monitoring of TXS 0506+056

We present our follow-up observations to search for an electromagnetic counterpart of the IceCube high-energy neutrino, IceCube-170922A. Monitoring observations of a likely counterpart, TXS 0506+056, are also described. First, we quickly took optical and near-infrared images of 7 flat-spectrum radio sources within the IceCube error region right after the neutrino detection and found a rapid flux decline of TXS 0506+056 in Kanata/HONIR J-band data. Motivated by this discovery, intensive follow-up observations of TXS 0506+056 are continuously done, including our monitoring imaging observations, spectroscopic observations, and polarimetric observations in optical and near-infrared wavelengths. TXS 0506+056 shows a large amplitude (~1.0 mag) variability in a time scale of several days or longer, although no significant variability is detected in a time scale of a day or shorter. TXS 0506+056 also shows a bluer-when-brighter trend in optical and near-infrared wavelengths. Structure functions of variabilities are examined and indicate that TXS 0506+056 is not a special blazar in terms of optical variability. Polarization measurement results of TXS 0506+056 are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

From supernova to supernova remnant: comparison of thermonuclear explosion models

Progress in the three-dimensional modeling of supernovae (SN) prompts us to revisit the supernova remnant (SNR) phase. We continue our study of the imprint of a thermonuclear explosion on the SNR it produces, that we started with a delayed-detonation model of a Chandrasekhar-mass white dwarf. Here we compare two different types of explosion models, each with two variants: two delayed detonation models (N100ddt, N5ddt) and two pure deflagration models (N100def, N5def), where the N number parametrizes the ignition. The output of each SN simulation is used as input of a SNR simulation carried on until 500 yr after the explosion. While all SNR models become more spherical over time and overall display the theoretical structure expected for a young SNR, clear differences are visible amongst the models, depending on the geometry of the ignition and on the presence or not of detonation fronts. Compared to N100 models, N5 models have a strong dipole component, and produce asymmetric remnants. N5def produces a regular-looking, but offset remnant, while N5ddt produces a two-sided remnant. Pure deflagration models exhibit specific traits: a central over-density, because of the incomplete explosion, and a network of seam lines across the surface, boundaries between burning cells. Signatures from the SN dominate the morphology of the SNR up to 100 yr to 300 yr after the explosion, depending on the model, and are still measurable at 500 yr, which may provide a way of testing explosion models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Photohadronic model for the neutrino and $γ$-ray emission from TXS 0506+056

The detection of a high energy muon neutrino on 22 September 2017 by IceCube neutrino detector coincides with the multiwavelength flaring from the BL Lac object TXS 0506+056, most likely confirming AGN as a source of high energy cosmic rays and neutrinos. Using the photohadronic model, we have explained the very high energy $γ$-rays observed by MAGIC telescopes few days after the neutrino event and extend the model to calculate the neutrino flux at different windows consistent with the flaring period of TXS 0506+056 and compared with the IceCube and MAGIC estimates. We also use this model to estimate the neutrino flux from the flaring of FSRQ PKS B1424-418 which is believed to be associated with the 2 PeV neutrino event observed by IceCube.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Matter Mixing in Aspherical Core-collapse Supernovae: Three-dimensional Simulations with Single Star and Binary Merger Progenitor Models for SN 1987A

We perform three-dimensional hydrodynamic simulations of aspherical core-collapse supernovae focusing on the matter mixing in SN 1987A. The impacts of four progenitor (pre-supernova) models and parameterized aspherical explosions are investigated. The four pre-supernova models include a blue supergiant (BSG) model based on a slow merger scenario developed recently for the progenitor of SN 1987A (Urushibata et al. 2018). The others are a BSG model based on a single star evolution and two red supergiant (RSG) models. Among the investigated explosion (simulation) models, a model with the binary merger progenitor model and with an asymmetric bipolar-like explosion, which invokes a jetlike explosion, best reproduces constraints on the mass of high velocity $^{56}$Ni, as inferred from the observed [Fe II] line profiles. The advantage of the binary merger progenitor model for the matter mixing is the flat and less extended $ρ\,r^3$ profile of the C+O core and the helium layer, which may be characterized by the small helium core mass. From the best explosion model, the direction of the bipolar explosion axis (the strongest explosion direction), the neutron star (NS) kick velocity, and its direction are predicted. Other related implications and future prospects are also given.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Monte-Carlo simulations of relativistic radiation mediated shocks: II. photon starved regime

Radiation mediated shocks (RMS) play a key role in shaping the early emission observed in many transients. In most cases, e.g., shock breakout in supernovae, llGRBs and neutron star mergers, the upstream plasma is devoid of radiation, and the photons that ultimately reach the observer are generated predominantly inside and downstream of the shock. Predicting the observed spectrum requires detailed calculations of the shock structure and thermodynamic state that account properly for the shock microphysics. We present results of self-consistent Monte-Carlo simulations of photon-starved RMS, that yield the shock structure and emission for a broad range of shock velocities, from sub-relativistic ($β_{sh} = 0.1$) to highly relativistic ($Γ_{sh} = 20$). Our simulations confirm that in relativistic RMS the immediate downstream temperature is regulated by exponential pair creation, ranging from $50$ keV at $β_{sh}=0.5$ to $200$ keV at $Γ_{sh}=20$. At lower velocities the temperature becomes sensitive to the shock velocity, with $kT\sim 0.5$ keV at $β_{sh}=0.1$. We also confirm that in relativistic shocks the opacity is completely dominated by newly created pairs, which has important implications for the breakout physics. We find the transition to pair dominance to occur at $β_{sh}=0.5$ roughly. In all cases examined, the spectrum below the $νF_ν$ peak has been found to be substantially softer than the Planck distribution. This has important implications for the optical emission in fast and relativistic breakouts, and their detection. The applications to GRB 060218 and GRB 170817A are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The VHE SED modelling of Markarian 501 in 2009

The high energy blazar, Markarian 501 was observed as a part of multi-instrument and multiwavelength campaign spanning the whole electromagnetic spectrum for 4.5 months during March 15 to August 1, 2009. On May 1, Whipple 10m telescope observed a very strong $γ$-ray flare in a time interval of about 0.5 h. Apart from this flare, high state and low state emissions were also observed by Whipple, VERITAS and MAGIC telescopes. Using the photohadronic model and accounting for the absorption of the extragalactic background light (EBL) to these very high energy $γ$-rays, excellent fits are obtained for the observed spectra. We have shown that the intrinsic spectrum for low state emission is flat, however, for high and very high states this is a power-law with slowly increasing exponent.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Monte Carlo Bayesian Search for the Plausible Source of the Telescope Array Hotspot

The Telescope Array (TA) collaboration has reported a hotspot of 19 ultrahigh-energy cosmic rays (UHECRs). Using a universal model with one source and energy-dependent magnetic deflections, we show that the distribution of the TA hotspot events is consistent with a single source hypothesis, although multiple sources cannot be ruled out. The chance probability of this distribution arising from a homogeneous distribution is $0.2\%$. We describe a Monte Carlo Bayesian (MCB) inference approach, which can be used to derive parameters of the magnetic fields as well as the source coordinates, and we apply this method to the TA hotspot data, inferring the location of the likely source. We discuss possible applications of the same approach to future data.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Study of the Gamma-Ray Burst Fundamental Plane

A class of long gamma-ray bursts (GRBs) with a plateau phase in their X-ray afterglows obeys a three-dimensional (3D) relation (Dainotti et al. 2016), between the rest-frame time at the end of the plateau, Ta, its corresponding X-ray luminosity, La, and the peak luminosity in the prompt emission, Lpeak. We extended the original analysis with X-ray data from July 2014 to July 2016 achieving a total sample of 183 Swift GRBs with afterglow plateaus and known redshifts. We added the most recent GRBs to the previous 'gold sample' (now including 45 GRBs) and obtained a relation plane with intrinsic scatter compatible within one sigma with the previous result. We compared several GRB categories, such as short with extended emission, X-ray Flashes, GRBs associated with SNe, long-duration GRBs, and the gold sample, composed only by GRBs with light curves with good data coverage and relatively flat plateaus and evaluated their relation planes. We found that they are not statistically different from the fundamental plane derived from the gold sample and that the fundamental plane still has the smallest scatter. We compared the jet opening angles tabulated in literature with the angles derived using the Eiso-Egamma relation of the method in Pescalli et al. (2015) and calculated the relation plane for a sample of long GRBs accounting for the different jet opening angles. We observed that this correction does not significantly reduce the scatter. In an extended analysis, we found that the fundamental plane is independent from several prompt and afterglow parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Shock wave structure in astrophysical flows with an account of photon transfer

For an accurate treatment of the shock wave propagation in high-energy astrophysical phenomena, such as supernova shock breakouts, gamma-ray bursts and accretion disks, knowledge of radiative transfer plays a crucial role. In this paper we consider one-dimensional (1D) special relativistic radiation hydrodynamics by solving the Boltzmann equation for radiative transfer. The structure of a radiative shock is calculated for a number of shock tube problems, including strong shock waves, and relativistic- and radiation-dominated cases. Calculations are performed using an iterative technique that consistently solves the equations of relativistic hydrodynamics and relativistic comoving radiative transfer. A comparison of radiative transfer solutions with the Eddington approximation and the M1 closure is made. A qualitative analysis of moment equations for radiation is performed and the conditions for the existence of jump discontinuity for non-relativistic cases are investigated numerically.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Are Models for Core-Collapse Supernova Progenitors Consistent with the Properties of Supernova Remnants?

The recent discovery that the Fe-K line luminosities and energy centroids observed in nearby SNRs are a strong discriminant of both progenitor type and circumstellar environment has implications for our understanding of supernova progenitor evolution. Using models for the chemical composition of core-collapse supernova ejecta, we model the dynamics and thermal X-ray emission from shocked ejecta and circumstellar material, modeled as an $r^{-2}$ wind, to ages of 3000 years. We compare the X-ray spectra expected from these models to observations made with the Suzaku satellite. We also model the dynamics and X-ray emission from Type Ia progenitor models. We find a clear distinction in Fe-K line energy centroid between core-collapse and Type Ia models. The core-collapse supernova models predict higher Fe-K line centroid energies than the Type Ia models, in agreement with observations. We argue that the higher line centroids are a consequence of the increased densities found in the circumstellar environment created by the expansion of the slow-moving wind from the massive progenitors.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Evidence for Thermal X-Ray Line Emission from the Synchrotron-Dominated Supernova Remnant RX J1713.7-3946

We report the first detection of thermal X-ray line emission from the supernova remnant (SNR) RX J1713.7-3946, the prototype of the small class of synchrotron dominated SNRs. A softness-ratio map generated using XMM-Newton data shows that faint interior regions are softer than bright shell regions. Using Suzaku and deep XMM-Newton observations, we have extracted X-ray spectra from the softest area, finding clear line features at 1 keV and 1.35 keV. These lines can be best explained as Ne Ly-alpha and Mg He-alpha from a thermal emission component. Since the abundance ratios of metals to Fe are much higher than solar values in the thermal component, we attribute the thermal emission to reverse-shocked SN ejecta. The measured Mg/Ne, Si/Ne, and Fe/Ne ratios of 2.0-2.6, 1.5-2.0, and <0.05 solar suggest that the progenitor star of RX J1713.7-3946 was a relatively low-mass star (<~20 M_sun), consistent with a previous inference based on the effect of stellar winds of the progenitor star on the surrounding medium. Since the mean blastwave speed of ~6000 km/s (the radius of 9.6 pc divided by the age of 1600 yr) is relatively fast compared with other core-collapse SNRs, we propose that RX J1713.7-3946 is a result of a Type Ib/c supernova whose progenitor was a member of an interacting binary. While our analysis provides strong evidence for X-ray line emission, our interpretation of its nature as thermal emission from SN ejecta requires further confirmation especially through future precision spectroscopic measurements using ASTRO-H.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Modeling Bright Gamma-ray and Radio Emission at Fast Cloud Shocks

Recent observations by the Large Area Telescope (LAT) onboard the Fermi satellite have revealed bright gamma-ray emission from middle-aged supernova remnants (SNRs) inside our Galaxy. These remnants, which also possess bright non-thermal radio shells, are often found to be interacting directly with surrounding gas clouds. We explore the non-thermal emission mechanism at these dynamically evolved SNRs by constructing a hydrodynamical model. Two scenarios of particle acceleration, either a re-acceleration of Galactic cosmic rays (CRs) or an efficient nonlinear diffusive shock acceleration (NLDSA) of particles injected from downstream, are considered. Using parameters inferred from observations, our models are contrasted with the observed spectra of SNR W44. For the re-acceleration case, we predict a significant enhancement of radio and GeV emission as the SNR undergoes a transition into the radiative phase. If sufficiently strong magnetic turbulence is present in the molecular cloud, the re-acceleration scenario can explain the observed broadband spectral properties. The NLDSA scenario also succeeds in explaining the $γ$-ray spectrum but fails to reproduce the radio spectral index. Efficient NLDSA also results in a significant post-shock non-thermal pressure that limits the compression during cooling and prevents the formation of a prominent dense shell. Some other interesting differences between the two models in hydrodynamical behavior and resulting spectral features are illustrated.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Particle acceleration in superluminal strong waves

We calculate the electron acceleration in random superluminal strong waves (SLSWs) and radiation from them by using numerical methods in the context of the termination shock of the pulsar wind nebulae. We pursue the electrons by solving the equation of motion in the analytically expressed electromagnetic turbulences. These consist of primary SLSW and isotropically distributed secondary electromagnetic waves. Under the dominance of the secondary waves, all electrons gain nearly equal energy. On the other hand, when the primary wave is dominant, selective acceleration occurs. The phase of the primary wave felt by the electrons moving nearly along the wavevector changes very slowly compared to the oscillation of the wave, which is called "phase locked", and such electrons are continuously accelerated. This acceleration by SLSWs may play a crucial role in the pre-acceleration for the shock acceleration. In general, the radiation from the phase-locked population is different from the synchro-Compton radiation. However, when the amplitude of the secondary waves is not extremely weaker than that of the primary wave, the typical frequency can be estimated from the synchro-Compton theory by using the secondary waves. The primary wave does not contribute to the radiation, because the SLSW accelerates electrons almost linearly. This radiation can be observed as a radio knot at the upstream of the termination shock of the pulsar wind nebulae without counter parts in higher frequency range.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Photospheric Emission from Collapsar Jets in 3D Relativistic Hydrodynamics

We explore the photospheric emission from a relativistic jet breaking out from a massive stellar envelope based on relativistic hydrodynamical simulations and post-process radiation transfer calculations in three dimensions. To investigate the impact of three-dimensional (3D) dynamics on the emission, two models of injection conditions are considered for the jet at the center of the progenitor star: one with periodic precession and another without precession. We show that structures developed within the jet due to the interaction with the stellar envelope, as well as due to the precession, have a significant imprint on the resulting emission. Particularly, we find that the signature of precession activity by the central engine is not smeared out and can be directly observed in the light curve as a periodic signal. We also show non-thermal features, which can account for observations of gamma-ray bursts, are produced in the resulting spectra, even though only thermal photons are injected initially and the effect of non-thermal particles is not considered.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Selection effects in Gamma Ray Bursts correlations: consequences on the ratio between GRB and star formation rates

Gamma Ray Bursts (GRBs) visible up to very high redshift have become attractive targets as potential new distance indicators. It is still not clear whether the relations proposed so far originate from an unknown GRB physics or result from selection effects. We investigate this issue in the case of the $L_X-T^*_a$ correlation (hereafter LT) between the X-ray luminosity $L_X (T_a)$ at the end of the plateau phase, $T_a$, and the rest frame time $T^{*}_a$. We devise a general method to build mock data sets starting from a GRB world model and taking into account selection effects on both time and luminosity. This method shows how not knowing the efficiency function could influence the evaluation of the intrinsic slope of any correlation and the GRB density rate. We investigate biases (small offsets in slope or normalization) that would occur in the LT relation as a result of truncations, possibly present in the intrinsic distributions of $L_X$ and $T^*_a$. We compare these results with the ones in Dainotti et al. (2013) showing that in both cases the intrinsic slope of the LT correlation is $\approx -1.0$. This method is general, therefore relevant to investigate if any other GRB correlation is generated by the biases themselves. Moreover, because the farthest GRBs and star-forming galaxies probe the reionization epoch, we evaluate the redshift-dependent ratio $Ψ(z)=(1+z)^α$ of the GRB rate to star formation rate (SFR). We found a modest evolution $-0.2\leq α\leq 0.5$ consistent with Swift GRB afterglow plateau in the redshift range $0.99<z<9.4$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Anisotropy expectations for ultra-high-energy cosmic rays with future high statistics experiments

UHECRs have attracted a lot of attention due to their challengingly high energies and their potential value to constrain physical processes and astrophysical parameters in the most energetic sources of the universe. Current detectors have failed to detect significant anisotropies which had been expected to allow source identification. Some indications about the UHECR composition, which may become heavier at the highest energies, has even put into question the possibility that such a goal could be achieved soon. We investigate the potential value of a new-generation detector, with 10 times larger exposure, to overcome the current situation and make significant progress in the detection of anisotropies and thus in the study of UHECRs. We take as an example the expected performances of the JEM-EUSO, assuming a uniform full-sky coverage with a total exposure of 300,000 km2 sr yr. We simulate realistic UHECR sky maps for a wide range of possible astrophysical scenarios allowed by the current constraints, taking into account the energy losses and photo-dissociation of the UHECRs, as well as their deflections by magnetic fields. These sky maps, built for the expected statistics of JEM-EUSO as well as for the current Auger statistics, as a reference, are analyzed from the point of view of their intrinsic anisotropies, using the two-point correlation function. A statistical study of the resulting anisotropies is performed for each astrophysical scenario, varying the UHECR source composition and spectrum as well as the source density. We find that significant anisotropies are expected to be detected by a next-generation UHECR detector, for essentially all the astrophysical scenarios studied, and give precise, quantitative meaning to this statement. Our results show that a gain of one order of magnitude in exposure would make a significant difference compared to the existing detectors.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Reverse and Forward Shock X-ray Emission in an Evolutionary Model of Supernova Remnants undergoing Efficient Diffusive Shock Acceleration

We present new models for the forward and reverse shock thermal X-ray emission from core-collapse and Type Ia supernova remnants (SNRs) which include the efficient production of cosmic rays via non-linear diffusive shock acceleration (DSA). Our CR-hydro-NEI code takes into account non-equilibrium ionization (NEI), hydrodynamic effects of efficient CR production on the SNR evolution, and collisional temperature equilibration among heavy ions and electrons in both the shocked supernova (SN) ejecta and the shocked circumstellar material. While X-ray emission is emphasized here, our code self-consistently determines both thermal and non-thermal broadband emission from radio to TeV energies. We include Doppler broadening of the spectral lines by thermal motions of the ions and by the remnant expansion. We study, in general terms, the roles which the ambient environment, progenitor models, temperature equilibration, and processes related to DSA have on the thermal and non-thermal spectra. The study of X-ray line emission from young SNRs is a powerful tool for determining specific SN elemental contributions, and for providing critical information that helps to understand the type and energetics of the explosion, the composition of the ambient medium in which the SN exploded, and the ionization and dynamics of the hot plasma in the shocked SN ejecta and interstellar medium. With the approaching launch of the next-generation X-ray satellite Astro-H, observations of spectral lines with unprecedented high resolution will become a reality. Our self-consistent calculations of the X-ray spectra from various progenitors will help interpret future observations of SNRs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spectral and Polarization Properties of Photospheric Emission From Stratified Jets

We explore the spectral and polarization properties of photospheric emissions from stratified jets in which multiple components, separated by a sharp velocity shear regions, are distributed in lateral direction. Propagation of thermal photons injected at high optical depth region are calculated until they escape from the photosphere. It is found that presence of the lateral structure within the jet leads to non-thermal feature of the spectra and significant polarization signal in the resulting emission. The deviation from thermal spectra as well as the polarization degree tends to be enhanced as the velocity gradient in the shear region increases. In particular, we show that emissions from multi-component jet can reproduce the typical observed spectra of gamma-ray bursts (GRBs) irrespective to the position of the observer when a velocity shear region is closely spaced in various lateral ($θ$) positions. The degree of polarization associated in the emission is significant (> few%) at wide range of observer angles and can be higher than 30%.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A CR-hydro-NEI Model of Multi-wavelength Emission from the Vela Jr. Supernova Remnant (SNR RX J0852.0-4622)

Based largely on energy budget considerations and the observed cosmic-ray (CR) ionic composition, supernova remnant (SNR) blast waves are the most likely sources of CR ions with energies at least up to the "knee" near 3 PeV. Shocks in young shell-type TeV-bright SNRs are surely producing TeV particles, but the emission could be dominated by ions producing neutral pion-decay emission or electrons producing inverse-Compton gamma-rays. Unambiguously identifying the GeV-TeV emission process in a particular SNR will not only help pin down the origin of CRs, it will add significantly to our understanding of the diffusive shock acceleration (DSA) mechanism and improve our understanding of supernovae and the impact SNRs have on the circumstellar medium. In this study, we investigate the Vela Jr. SNR, an example of TeV-bright non-thermal SNRs. We perform hydrodynamic simulations coupled with non-linear DSA and non-equilibrium ionization near the forward shock (FS) to confront currently available multi-wavelength data. We find, with an analysis similar to that used earlier for SNR RX J1713.7-3946, that self-consistently modeling the thermal X-ray line emission with the non-thermal continuum in our one-dimensional model strongly constrains the fitting parameters, and this leads convincingly to a leptonic origin for the GeV-TeV emission for Vela Jr. This conclusion is further supported by applying additional constraints from observation, including the radial brightness profiles of the SNR shell in TeV gamma-rays, and the spatial variation of the X-ray synchrotron spectral index. We will discuss implications of our models on future observations by the next-generation telescopes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Matter mixing in aspherical core-collapse supernovae: a search for possible conditions for conveying $^{56}$Ni into high velocity regions

We perform two-dimensional axisymmetric hydrodynamic simulations of matter mixing in aspherical core-collapse supernova explosions of a 16.3 $M_{\odot}$ star with a compact hydrogen envelope. Observations of SN 1987A have provided evidence that $^{56}$Ni synthesized by explosive nucleosynthesis is mixed into fast moving matter ($\gtrsim$ 3,500 km s$^{-1}$) in the exploding star. In order to clarify the key conditions for reproducing such high velocity of $^{56}$Ni, we revisit matter mixing in aspherical core-collapse supernova explosions. Explosions are initiated artificially by injecting thermal and kinetic energies around the interface between the iron core and the silicon-rich layer. Perturbations of 5% or 30% amplitude in the radial velocities are introduced at several points in time. We found that no high velocity $^{56}$Ni can be obtained if we consider bipolar explosions with perturbations (5% amplitude) of pre-supernova origins. If large perturbations (30% amplitude) are introduced or exist due to some unknown mechanism in a later phase just before the shock wave reaches the hydrogen envelope, $^{56}$Ni with a velocity of 3,000 km s$^{-1}$ can be obtained. Aspherical explosions that are asymmetric across the equatorial plane with clumpy structures in the initial shock waves are investigated. We found that the clump sizes affect the penetration of $^{56}$Ni. Finally, we report that an aspherical explosion model that is asymmetric across the equatorial plane with multiple perturbations of pre-supernova origins can cause the penetration of $^{56}$Ni clumps into fast moving matter of 3,000 km s$^{-1}$. We show that both aspherical explosion with clumpy structures and perturbations of pre-supernova origins may be necessary to reproduce the observed high velocity of $^{56}$Ni. To confirm this, more robust three-dimensional simulations are required.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Photospheric emission from stratified jets

We explore photospheric emissions from stratified two-component jets, wherein a highly relativistic spine outflow is surrounded by a wider and less relativistic sheath outflow. Thermal photons are injected in regions of high optical depth and propagated until they escape at the photosphere. Due to the presence of shear in velocity (Lorentz factor) at the boundary of the spine and sheath region, a fraction of the injected photons are accelerated via a Fermi-like acceleration mechanism such that a high energy power-law tail is formed in the resultant spectrum. We show, in particular, that if a velocity shear with a considerable variance in the bulk Lorentz factor is present, the high energy part of observed Gamma-ray Bursts (GRBs) photon spectrum can be explained by this photon acceleration mechanism. We also show that the accelerated photons may also account for the origin of the extra hard power-law component above the bump of the thermal-like peak seen in some peculiar bursts (e.g., GRB 090510, 090902B, 090926A). It is demonstrated that time-integrated spectra can also reproduce the low energy spectrum of GRBs consistently due to a multi-temperature effect when time evolution of the outflow is considered. Finally, we show that the empirical Ep-Lp relation can be explained by differences in the outflow properties of individual sources.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Generalized Model of Nonlinear Diffusive Shock Acceleration Coupled to an Evolving Supernova Remnant

To better model the efficient production of cosmic rays (CRs) in supernova remnants (SNRs) with the associated coupling between CR production and SNR dynamics, we have generalized an existing cr-hydro-NEI code (i.e., Ellison et al. 2012) to include the following processes: (1) an explicit calculation of the upstream precursor structure including the position dependent flow speed, density, temperature, and magnetic field strength; (2) a momentum and space dependent CR diffusion coefficient; (3) an explicit calculation of magnetic field amplification (MFA); (4) calculation of the maximum CR momentum using the amplified magnetic field; (5) a finite Alfven speed for the particle scattering centers; and (6) the ability to accelerate a superthermal seed population of CRs as well as the ambient thermal plasma. While a great deal of work has been done modeling SNRs, most work has concentrated on either the continuum emission from relativistic electrons or ions, or the thermal emission from the shock heated plasma. Our generalized code combines these elements and describes the interplay between CR production and SNR evolution, including the nonlinear coupling of efficient diffusive shock acceleration (DSA), based mainly on the work of P. Blasi and co-workers, and a non-equilibrium ionization (NEI) calculation of thermal X-ray line emission. We believe our generalized model will provide a consistent modeling platform for SNRs, including those interacting with molecular clouds, and improve the interpretation of current and future observations, including the high-quality spectra expected from Astro-H. SNR RX J1713.7-3946 is modeled as an example.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

IceCube Non-detection of GRBs: Constraints on the Fireball Properties

The increasingly deep limit on the neutrino emission from gamma-ray bursts (GRBs) with IceCube observations has reached the level that could put useful constraints on the fireball properties. We first present a revised analytic calculation of the neutrino flux, which predicts a flux an order of magnitude lower than that obtained by the IceCube collaboration. For benchmark model parameters (e.g. the bulk Lorentz factor is Γ=10^{2.5}, the observed variability time for long GRBs is t_v=0.01 s and the ratio between the energy in accelerated protons and in radiation is η_p=10 for every burst) in the standard internal shock scenario, the predicted neutrino flux from 215 bursts during the period of the 40-string and 59-string configurations is found to be a factor of ~3 below the IceCube sensitivity. However, if we accept the recently found inherent relation between the bulk Lorentz factor and burst energy, the expected neutrino flux increases significantly and the spectral peak shifts to lower energy. In this case, the non-detection then implies that the baryon loading ratio should be η_p<10 if the variability time of long GRBs is fixed to t_v=0.01 s. Instead, if we relax the standard internal shock scenario but keep to assume η_p=10, the non-detection constrains the dissipation radius to be R>4x10^{12} cm assuming the same dissipation radius for every burst and benchmark parameters for fireballs. We also calculate the diffuse neutrino flux from GRBs for different luminosity functions existing in the literature. The expected flux exceeds the current IceCube limit for some luminosity functions, and thus the non-detection constrains η_p<10 in such cases when the variability time of long GRBs is fixed to t_v=0.01 s.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

GeV Gamma-ray Emission from the Binary PSR B1259-63/SS2883 During the 2010 Periastron Passage

PSR B1259-63/SS2883 is a binary system which consists of a 48-ms radio pulsar and a massive star in a highly eccentric orbit with a period of about 3.4 years. Non-pulsed and non-thermal emissions from this binary have been reported in the radio, X-ray and TeV gamma-ray energy ranges. Light curves in the radio and X-ray bands showed characteristic double-peaked features which can be attibuted to the interactions of the pulsar wind and the Be disk during the crossings by the pulsar. The TeV light curves around periastron differ between 2004 and 2007 observations, and the feature is not conclusive. We report a detection of GeV gamma-ray emission around the periastron passage in December 2010 with Fermi-LAT.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hydrodynamic Interaction between the Be Star and the Pulsar in the TeV Binary PSR B1259-63/LS 2883

We study the interaction between the Be star and the pulsar in the TeV binary PSR B1259-63/LS 2883, using 3-D SPH simulations of the tidal and wind interactions in this Be-pulsar system. We first run a simulation without pulsar wind nor Be wind, taking into account only the gravitational effect of the pulsar on the Be disk. In this simulation, the gas particles are ejected at a constant rate from the equatorial surface of the Be star, which is tilted in a direction consistent with multi-waveband observations. We run the simulation until the Be disk is fully developed and starts to repeat a regular tidal interaction with the pulsar. Then, we turn on the pulsar wind and the Be wind. We run two simulations with different wind mass-loss rates for the Be star, one for a B2V type and the other for a significantly earlier spectral type. Although the global shape of the interaction surface between the pulsar wind and the Be wind agrees with the analytical solution, the effect of the pulsar wind on the Be disk is profound. The pulsar wind strips off an outer part of the Be disk, truncating the disk at a radius significantly smaller than the pulsar orbit. Our results, therefore, rule out the idea that the pulsar passes through the Be disk around periastron, which has been assumed in the previous studies. It also turns out that the location of the contact discontinuity can be significantly different between phases when the pulsar wind directly hits the Be disk and those when the pulsar wind collides with the Be wind. It is thus important to adequately take into account the circumstellar environment of the Be star, in order to construct a satisfactory model for this prototypical TeV binary.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Rotating BHs as Central Engines of Long GRBs: Faster is Better

We performed simulations of collapsars with different Kerr parameters a=0, 0.5, 0.9, 0.95. It is shown that a more rapidly rotating black hole is driving a more energetic jet. No jet is seen for the case of Schwartzschild black hole case, while the total energy of the jet is as large as 10^50 erg for a rapidly rotating Kerr black hole case (a=0.95). In order to explain the high luminosity of a GRB, it is concluded that a rapidly rotating black hole is favored ('faster is better'). We also find in the case of a=0.95 that (i) the stagnation region is clearly found in the jet region, (ii) the ordered poloidal field lines are seen in the jet, (iii) the jet region is surrounded by a 'Wall-like' structure that has a higher pressure than the jet region and contains strong vertical magnetic fields, and (iv) the jet is initiated by outgoing Poynting flux from the outer horizon of the black hole (Blandford-Znajek effect). The bulk Lorentz factor of the jet is still of the order of unity. However, energy density of electro-magnetic fields dominates the one of rest-mass in the jet. It can be expected that a relativistic jet will be seen if we perform a simulation for a longer time scale (of the order of 10-100 sec).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Thermal Radiation from GRB Jets

In this study, the light curves and spectrum of the photospheric thermal radiation from ultrarelativistic gamma-ray burst (GRB) jets are calculated using 2D relativistic hydrodynamic simulations of jets from a collapsar. As the jet advances, the density around the head of the jet decreases, and its Lorentz factor reaches as high as 200 at the photosphere and 400 inside the photosphere. For an on-axis observer, the photosphere appears concave shaped due to the low density and high beaming factor of the jet. The luminosity varies because of the abrupt change in the position of the photosphere due to the internal structure of the jet. Comparing our results with GRB090902B, the flux level of the thermal-like component is similar to our model, although the peak energy looks a little bit higher (but still within a factor of 2). From the comparison, we estimate that the bulk Lorentz factor of GRB090902B is $Γ\sim 2.4 \times 10^2 (r/10^{12} \rm cm$) where $r$ is the radius of the photosphere. The spectrum for an on-axis observer is harder than that for an off-axis observer. There is a time lag of a few seconds for high energy bands in the light curve. This may be the reason for the delayed onset of GeV emission seen in GRB080916C. The spectrum below the peak energy is a power law and the index is $2.3 \sim 2.6$ which is softer than that of single temperature plank distribution but still harder than that of typical value of observed one.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Origins of Short Gamma-Ray Bursts Deduced from Offsets to Their Host Galaxies Revisited

The spatial distribution of short Gamma-ray bursts (GRBs) in their host galaxies provide us an opportunity to investigate their origins. Based on the currently observed distribution of short GRBs relative to their host galaxies, we obtain the fraction of the component that traces the mergers of binary compact objects and the one that traces star formation rate (such as massive stars) in early- and late-type host galaxies. We find that the fraction of massive star component is 0.37\pm 0.13 from the analysis of projected offset distribution. This suggests that a good fraction of short GRBs still originate from merger events. From our analysis, we also conclude that the fraction of late-type hosts among the elliptical, starburst and spiral galaxy is 0.82\pm 0.05, which is consistent with the observed early- to late-type number ratio of host galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Role of Galactic sources and magnetic fields in forming the observed energy-dependent composition of ultrahigh-energy cosmic rays

Recent results from Pierre Auger Observatory, showing energy dependent chemical composition of ultrahigh-energy cosmic rays (UHECRs) with a growing fraction of heavy elements at high energies, suggest a possible non-negligible contribution of the Galactic sources. We show that, in the case of UHECRs produced by gamma-ray bursts (GRBs) or rare types of supernova explosions that took place in Milky Way in the past, the change in UHECR composition can result from the difference in diffusion times for different species. The anisotropy in the direction of the Galactic Center is expected to be a few per cent on average, but the locations of the most recent/closest bursts can be associated with observed clusters of UHECRs.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Solar-type Magnetic Reconnection Model for Magnetar Giant Flare

We present a theoretical model describing magnetar giant flares on the basis of solar flare/coronal mass ejection theory. In our model, a preflare activity plays a crucial role in driving evaporating flows and supplying baryonic matters into the magnetosphere. The loaded baryonic matter, that is called "prominence", is then gradually uplifted via crustal cracking with maintaining a quasi-force-free equilibrium of the magnetosphere. Finally the prominence is erupted by the magnetic pressure force due to the loss of equilibrium triggered by the explosive magnetic reconnection. The giant flare should be induced as a final outcome of the prominence eruption accompanied by large-scale field reconfigurations. An essential difference between solar and magnetar flares is the control process of their evolutionary dynamics. The flaring activity on magnetars is mainly controlled by the radiative process unlike the solar flare governed by the electron conduction. It is highly suggestive that our model is accountable for the physical properties of the extraordinary giant flare observed on 2004 December 27 from SGR1806-20, including the source of baryonic matters loaded in the expanding ejecta observed after the giant burst.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A peculiar jet and arc of molecular gas toward the rich and young stellar cluster Westerlund 2 and a TeV gamma ray source

We have discovered remarkable jet- and arc-like molecular features toward the rich and young stellar cluster Westerlund2. The jet has a length of ~100 pc and a width of ~10 pc, while the arc shows a crescent shape with a radius of ~30 pc. These molecular features each have masses of ~10000 solar mass and show spatial correlations with the surrounding lower density HI gas. The jet also shows an intriguing positional alignment with the core of the TeV gamma ray source HESS J1023-575 and with the MeV/GeV gamma-ray source recently reported by the Fermi collaboration. We argue that the jet and arc are caused by an energetic event in Westerlund 2, presumably due to an anisotropic supernova explosion of one of the most massive member stars. While the origin of the TeV and GeV gamma-ray sources is uncertain, one may speculate that they are related to the same event via relativistic particle acceleration by strong shock waves produced at the explosion or by remnant objects such as a pulsar wind nebula or microquasar.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

High-energy emission as a test of the prior emission model for gamma-ray burst afterglows

We study high-energy gamma-ray afterglow emission from gamma-ray bursts (GRBs) in the prior emission model, which is proposed to explain the plateau phase of the X-ray afterglow. This model predicts the high-energy gamma-ray emission when the prompt GRB photons from the main flow are up-scattered by relativistic electrons accelerated at the external shock due to the prior flow. The expected spectrum has the peak of 10-100 GeV at around the end time of the plateau phase for typical GRBs, and high-energy gamma rays from nearby and/or energetic GRBs can be detected by current and future Cherenkov telescopes such as MAGIC, VERITAS, CTA, and possibly Fermi. Multi-wavelength observations by ground-based optical telescopes as well as Fermi and/or Swift sattelites are important to constrain the model. Such external inverse-Compton emission may even lead to GeV-TeV gamma-ray signals with the delay time of 10-100 s, only if the plateau phase is short-lived.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Cosmic Rays Above the Second Knee from Clusters of Galaxies and Associated High-Energy Neutrino Emission

Accretion and merger shocks in clusters of galaxies are potential accelerators of high-energy protons, which can give rise to high-energy neutrinos through pp interactions with the intracluster gas. We discuss the possibility that protons from cluster shocks make a significant contribution to the observed cosmic rays in the energy range between the second knee around 10^17.5 eV and the ankle around 10^18.5 eV. The accompanying cumulative neutrino background above PeV may be detectable by upcoming neutrino telescopes such as IceCube or KM3NeT, providing a test of this scenario as well as a probe of cosmic-ray confinement properties in clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

High-energy cosmic-ray nuclei from high- and low-luminosity gamma-ray bursts and implications for multi-messenger astronomy

Gamma-ray bursts (GRBs) are one of the candidates of ultra-high-energy (around 10^18.5 eV) cosmic-ray (UHECR) sources. We investigate high-energy cosmic-ray acceleration including heavy nuclei in GRBs by using Geant 4, and discuss its various implications, taking both of high-luminosity (HL) and low-luminosity (LL) GRBs into account. This is because LL GRBs may also make a significant contribution to the observed UHECR flux if they form a distinct population. We show that not only protons but also heavier nuclei can be accelerated up to ultra-high energies in the internal, (external) reverse and forward shock models. We also show that the condition for ultra-high-energy heavy nuclei such as iron to survive is almost the same as that for about TeV gamma-rays to escape from the source and for high-energy neutrinos not to be much produced. The multi-messenger astronomy by neutrino and GeV-TeV gamma-ray telescopes such as IceCube and KM3Net, GLAST and MAGIC will be important to see whether GRBs can be accelerators of ultra-high-energy heavy nuclei. We also demonstrate expected spectra of high-energy neutrinos and gamma rays, and discuss their detectabilities. In addition, we discuss implictaions of the GRB-UHECR hypothesis. We point out, since the number densities of HL-GRBs and LL-GRBs are quite different, its detemination by UHECR observations is also important.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Effects of Cosmic Infrared Background on High Energy Delayed Gamma-Rays from Gamma-Ray Bursts

Regenerated high-energy emissions from gamma-ray bursts (GRBs) are studied in detail. If the intrinsic primary spectrum extends to the TeV range, these very high-energy photons are absorbed by the cosmic infrared background (CIB). Created high-energy electron-positron pairs up-scatter mainly cosmic microwave background (CMB) photons, and secondary photons are generated in the GeV-TeV range. These secondary delayed photons may be observed in the near future, which are useful for a consistency check of the intrinsic primary spectrum. In this paper, we focus on effects of the CIB on delayed secondary emissions. In particular, we show that not only up-scattered CMB photons but also up-scattered CIB ones are important, especially for low redshift bursts. They also give us additional information on the CIB, whose photon density is not definitely determined so far.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

The effect of energy amplification variance on the shock-acceleration

The shock-acceleration theory predicts a power-law energy spectrum in the test particle approximation, and there are two ways to calculate a power-law index, Peacock's approximation and Vietri's formulation. In Peacock's approximation, it is assumed that particles cross a shock front many times and energy-gains for each step are fully uncorrelated. On the other hand, correlation of the distribution of an energy-gain factor for a particle is considered in Vietri's formulation. We examine how Peacock's approximation differs from Vietri's formulation. It is useful to know when we can use Peacock's approximation because Peacock's approximation is simple to derive the power-law index. In addition, we focus on how the variance of the energy-gain factor has an influence on the difference between Vietri's formulation and Peacock's approximation. The effect of the variance has not been examined well until now. For demonstration, we consider two cases for the scattering in the upstream: the large-angle scattering (model A) and the regular deflection by large-scale magnetic fields (model B). Especially there is no correlation among the distribution of an energy-gain factor for every step in model A. In this model, we see the power-law index derived from Peacock's approximation differs from the one derived from Vietri's formulation when we consider the mildly-relativistic shock, and the variance of the energy-gain factor affects this difference. We can use Peacock's approximation for a non-relativistic shock and a highly-relativistic shock because the effect of the variance is hidden. In model B, we see the difference of the power-law converging along the shock velocity.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Unstable GRB photospheres and electron-positron annihilation lines

We propose an emission mechanism of prompt gamma-ray bursts (GRBs) that can reproduce the observed non-thermal spectra with high radiative efficiencies, >50%. Internal dissipation below a photosphere can create a radiation-dominated thermal fireball. If electron-positron pairs outnumber protons, radiative acceleration of pairs drives the two-stream instabilities between pairs and protons, leading to the ``proton sedimentation'' in the accelerating pair frame. Pairs are continuously shock heated by proton clumps, scattering the thermal photons into the broken power-law shape, with the non-thermal energy that is comparable to the proton kinetic energy, consistent with observations. Pair photospheres become unstable around the radius of the progenitor star where strong thermalization occurs, if parameters satisfy the observed spectral (Yonetoku) relation. Pair annihilation lines are predicted above continua, which could be verified by GLAST.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Explosive Nucleosynthesis in GRB Jets Accompanied by Hypernovae

Two-dimensional hydrodynamic simulations are performed to investigate explosive nucleosynthesis in a collapsar using the model of MacFadyen and Woosley (1999). It is shown that 56Ni is not produced in the jet of the collapsar sufficiently to explain the observed amount of a hypernova when the duration of the explosion is \sim 10 sec, which is considered to be the typical timescale of explosion in the collapsar model. Even though a considerable amount of 56Ni is synthesized if all explosion energy is deposited initially, the opening angles of the jets become too wide to realize highly relativistic outflows and gamma-ray bursts in such a case. From these results, it is concluded that the origin of 56Ni in hypernovae associated with GRBs is not the explosive nucleosynthesis in the jet. We consider that the idea that the origin is the explosive nucleosynthesis in the accretion disk is more promising. We also show that the explosion becomes bi-polar naturally due to the effect of the deformed progenitor. This fact suggests that the 56Ni synthesized in the accretion disk and conveyed as outflows are blown along to the rotation axis, which will explain the line features of SN 1998bw and double peaked line features of SN 2003jd. Some fraction of the gamma-ray lines from 56Ni decays in the jet will appear without losing their energies because the jet becomes optically thin before a considerable amount of 56Ni decays as long as the jet is a relativistic flow. We show that abundance of nuclei whose mass number \sim 40 in the ejecta depends sensitively on the energy deposition rate. So it may be determined by observations of chemical composition in metal poor stars which model is the proper one as a model of a gamma-ray burst accompanied by a hypernova.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

High Energy Neutrinos and Cosmic-Rays from Low-Luminosity Gamma-Ray Bursts?

The recently discovered gamma-ray burst (GRB) GRB 060218/SN 2006aj is classified as an X-ray Flash with very long duration driven possibly by a neutron star. Since GRB 060218 is very near about 140 Mpc and very dim, one-year observation by Swift suggests that the rate of GRB 060218-like events might be very high so that such low luminosity GRBs (LL-GRBs) might form a different population from the cosmological high luminosity GRBs (HL-GRBs). We found that the high energy neutrino background from LL-GRBs could be comparable with that from HL-GRBs. If each neutrino event is detected by IceCube, later optical-infrared follow-up observations such as by Subaru and HST have possibilities to identify a Type Ibc supernova associated with LL-GRBs, even if gamma- and X-rays are not observed by Swift. This is in a sense a new window from neutrino astronomy, which might enable us to confirm the existence of LL-GRBs and to obtain information about their rate and origin. We also argue LL-GRBs as high energy gamma-ray and cosmic-ray sources.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Alfven Wave-Driven Proto-Neutron Star Winds And R-Process Nucleosynthesis

We propose magnetic proto-neutron star (PNS) winds driven by \Alfven waves as well as the neutrino heating as an appropriate site for the r-process nucleosynthesis. \Alfven waves excited by surface motions of a PNS propagate outwardly, and they heat and accelerate the wind by dissipation. Compared with the wind purely driven by the neutrino heating, larger entropy per baryon and shorter dynamical time scale are achieved, which favors the r-process. We study reasonable cases that a wave amplitude is 10% of \Alfven speed at the surface to find that a PNS with surface field strength, $\gtrsim 5\times 10^{14}$G, gives suitable wind properties for the r-process, provided that dissipation length of the wave is at most $\sim$ 10 times of the PNS radius. We also compare properties of transcritical and subcritical winds in light of the r-process. We finally discuss possibilities of detections of $γ$-rays from radioactive nuclei and absorption lines due to Ba in supernova remnants which possess magnetars.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Effects of Magnetic Fields on Proto-Neutron Star Winds

We discuss effects of magnetic fields on proto-neutron star winds by performing numerical simulation. We assume that the atmosphere of proto-neutron star has a homogenous magnetic field (ranging from ~10^{12} G to ~10^{15} G) perpendicular to the radial direction and examine the dependence of the three key quantities (dynamical time scale, electron fraction, and entropy per baryon) for the successful r-process on the magnetic field strength. Our results show that even with a magneter-class field strength, ~10^{15} G, the feature of the wind dynamics varies only little from that of non-magnetic winds, and that the condition for successful r-process is not realized.

0 citations0 reviewsNo code linkedOpen access