BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Shigehiro Nagataki

Shigehiro Nagataki contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEastro-ph.SRhep-phastro-ph.COastro-ph.GA

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 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