Showing 1–50 of 230 results for author: Nagataki, S

GRB Redshift Classifier to Follow-up High-Redshift GRBs Using Supervised Machine Learning

Authors: Maria Giovanna Dainotti, Shubham Bhardwaj, Christopher Cook, Joshua Ange, Nishan Lamichhane, Malgorzata Bogdan, Monnie McGee, Pavel Nadolsky, Milind Sarkar, Agnieszka Pollo, Shigehiro Nagataki

Abstract: Gamma-ray bursts (GRBs) are intense, short-lived bursts of gamma-ray radiation observed up to a high redshift ($z \sim 10$) due to their luminosities. Thus, they can serve as cosmological tools to probe the early Universe. However, we need a large sample of high$-z$ GRBs, currently limited due to the difficulty in securing time at the large aperture Telescopes. Thus, it is painstaking to determine… ▽ More Gamma-ray bursts (GRBs) are intense, short-lived bursts of gamma-ray radiation observed up to a high redshift ($z \sim 10$) due to their luminosities. Thus, they can serve as cosmological tools to probe the early Universe. However, we need a large sample of high$-z$ GRBs, currently limited due to the difficulty in securing time at the large aperture Telescopes. Thus, it is painstaking to determine quickly whether a GRB is high$z$ or low$-z$, which hampers the possibility of performing rapid follow-up observations. Previous efforts to distinguish between high$-$ and low$-z$ GRBs using GRB properties and machine learning (ML) have resulted in limited sensitivity. In this study, we aim to improve this classification by employing an ensemble ML method on 251 GRBs with measured redshifts and plateaus observed by the Neil Gehrels Swift Observatory. Incorporating the plateau phase with the prompt emission, we have employed an ensemble of classification methods to enhance the sensitivity unprecedentedly. Additionally, we investigate the effectiveness of various classification methods using different redshift thresholds, $z_{threshold}$=$z_t$ at $z_{t}=$ 2.0, 2.5, 3.0, and 3.5. We achieve a sensitivity of 87\% and 89\% with a balanced sampling for both $z_{t}=3.0$ and $z_{t}=3.5$, respectively, representing a 9\% and 11\% increase in the sensitivity over Random Forest used alone. Overall, the best results are at $z_{t} = 3.5$, where the difference between the sensitivity of the training set and the test set is the smallest. This enhancement of the proposed method paves the way for new and intriguing follow-up observations of high$-z$ GRBs. △ Less

Submitted 16 August, 2024; originally announced August 2024.

Comments: 36 pages, 14 Figures (6 Figures with single panel, 4 Figures with 8 panels, 2 Figures with 4 panels, 1 Figure with 12 panels, 1 Figure with 2 panels), 6 Tables, submitted to ApJS

Neutron-quark stars: Discerning viable alternatives for the higher-density part of the equation of state of compact stars

Authors: Sudipta Hensh, Yong-Jia Huang, Toru Kojo, Luca Baiotti, Kentaro Takami, Shigehiro Nagataki, Hajime Sotani

Abstract: By taking into account the latest observations and theoretical constraints, we investigate the merger and post-merger of binary neutron stars (NSs) with general-relativistic numerical simulations employing hadronic and hybrid equations of state (EOSs). We name our hybrid stars neutron-quark stars (NQS), because the transition from hadrons to quarks starts at a density lower than the central densit… ▽ More By taking into account the latest observations and theoretical constraints, we investigate the merger and post-merger of binary neutron stars (NSs) with general-relativistic numerical simulations employing hadronic and hybrid equations of state (EOSs). We name our hybrid stars neutron-quark stars (NQS), because the transition from hadrons to quarks starts at a density lower than the central density of $\sim 1 M_{\odot}$ stars. We address two viable scenarios for the transition to quark matter: a crossover or a strong first-order phase transition (1PT). We find that a crossover transition is in principle observable when both the inspiral and post-merger signals are detected because the post-merger gravitational-wave (GW) main frequency $f_2$ is generally lower than that of hadronic models with the same tidal deformability ($Λ$). Since it is viable according to current multi-messenger constraints, we also highlight the possibility of an EOS with a strong 1PT that takes place at 1.8 times the nuclear saturation density, with a stiff quark EOS after the transition. It is the first time that mergers of binary NQSs with a deconfined quark-matter core are studied numerically in full general relativity. In this case, although ($Λ$, compactness) lies significantly outside the hadronic relation, ($Λ$, $f_2$) is close to the relation valid for hadronic EOSs. We also point out a linear correlation, valid within the observational constraints and not sensitive to the presence of a hadron-quark transition, between the emitted energy in GWs and their frequency. △ Less

Submitted 2 August, 2024; v1 submitted 12 July, 2024; originally announced July 2024.

Comments: 6+7 pages, 3+10 figures

Report number: RIKEN-iTHEMS-Report-24

A detailed study of the very-high-energy Crab pulsar emission with the LST-1

Authors: CTA-LST Project, :, K. Abe, S. Abe, A. Abhishek, F. Acero, A. Aguasca-Cabot, I. Agudo, N. Alvarez Crespo, L. A. Antonelli, C. Aramo, A. Arbet-Engels, C. Arcaro, M. Artero, K. Asano, P. Aubert, A. Baktash, A. Bamba, A. Baquero Larriva, L. Baroncelli, U. Barres de Almeida, J. A. Barrio, I. Batkovic, J. Baxter, J. Becerra González , et al. (272 additional authors not shown)

Abstract: Context: There are currently three pulsars firmly detected by imaging atmospheric Cherenkov telescopes (IACTs), two of them reaching TeV energies, challenging models of very-high-energy (VHE) emission in pulsars. More precise observations are needed to better characterize pulsar emission at these energies. The LST-1 is the prototype of the Large-Sized Telescope, that will be part of the Cherenkov… ▽ More Context: There are currently three pulsars firmly detected by imaging atmospheric Cherenkov telescopes (IACTs), two of them reaching TeV energies, challenging models of very-high-energy (VHE) emission in pulsars. More precise observations are needed to better characterize pulsar emission at these energies. The LST-1 is the prototype of the Large-Sized Telescope, that will be part of the Cherenkov Telescope Array Observatory (CTAO). Its improved performance over previous IACTs makes it well suited for studying pulsars. Aims: To study the Crab pulsar emission with the LST-1, improving and complementing the results from other telescopes. These observations can also be used to characterize the potential of the LST-1 to study other pulsars and detect new ones. Methods: We analyzed a total of $\sim$103 hours of gamma-ray observations of the Crab pulsar conducted with the LST-1 in the period from September 2020 to January 2023. The observations were carried out at zenith angles less than 50 degrees. A new analysis of the Fermi-LAT data was also performed, including $\sim$14 years of observations. Results: The Crab pulsar phaseogram, long-term light-curve, and phase-resolved spectra are reconstructed with the LST-1 from 20 GeV to 450 GeV for P1 and up to 700 GeV for P2. The pulsed emission is detected with a significance of 15.2$σ$. The two characteristic emission peaks of the Crab pulsar are clearly detected (>10$σ$), as well as the so-called bridge emission (5.7$σ$). We find that both peaks are well described by power laws, with spectral indices of $\sim$3.44 and $\sim$3.03 respectively. The joint analysis of Fermi-LAT and LST-1 data shows a good agreement between both instruments in the overlapping energy range. The detailed results obtained in the first observations of the Crab pulsar with LST-1 show the potential that CTAO will have to study this type of sources. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: Accepted by A&A

Isotropy of cosmic rays beyond $10^{20}$ eV favors their heavy mass composition

Authors: Telescope Array Collaboration, R. U. Abbasi, Y. Abe, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, I. Buckland, B. G. Cheon, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, N. Globus, R. Gonzalez, W. Hanlon, N. Hayashida, H. He , et al. (118 additional authors not shown)

Abstract: We report an estimation of the injected mass composition of ultra-high energy cosmic rays (UHECRs) at energies higher than 10 EeV. The composition is inferred from an energy-dependent sky distribution of UHECR events observed by the Telescope Array surface detector by comparing it to the Large Scale Structure of the local Universe. In the case of negligible extra-galactic magnetic fields the resul… ▽ More We report an estimation of the injected mass composition of ultra-high energy cosmic rays (UHECRs) at energies higher than 10 EeV. The composition is inferred from an energy-dependent sky distribution of UHECR events observed by the Telescope Array surface detector by comparing it to the Large Scale Structure of the local Universe. In the case of negligible extra-galactic magnetic fields the results are consistent with a relatively heavy injected composition at E ~ 10 EeV that becomes lighter up to E ~ 100 EeV, while the composition at E > 100 EeV is very heavy. The latter is true even in the presence of highest experimentally allowed extra-galactic magnetic fields, while the composition at lower energies can be light if a strong EGMF is present. The effect of the uncertainty in the galactic magnetic field on these results is subdominant. △ Less

Submitted 3 July, 2024; v1 submitted 27 June, 2024; originally announced June 2024.

Comments: 8 pages, 3 figures, accepted for publication in PRL

Mass composition of ultra-high energy cosmic rays from distribution of their arrival directions with the Telescope Array

Authors: Telescope Array Collaboration, R. U. Abbasi, Y. Abe, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, I. Buckland, B. G. Cheon, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, N. Globus, R. Gonzalez, W. Hanlon, N. Hayashida, H. He , et al. (118 additional authors not shown)

Abstract: We use a new method to estimate the injected mass composition of ultrahigh cosmic rays (UHECRs) at energies higher than 10 EeV. The method is based on comparison of the energy-dependent distribution of cosmic ray arrival directions as measured by the Telescope Array experiment (TA) with that calculated in a given putative model of UHECR under the assumption that sources trace the large-scale struc… ▽ More We use a new method to estimate the injected mass composition of ultrahigh cosmic rays (UHECRs) at energies higher than 10 EeV. The method is based on comparison of the energy-dependent distribution of cosmic ray arrival directions as measured by the Telescope Array experiment (TA) with that calculated in a given putative model of UHECR under the assumption that sources trace the large-scale structure (LSS) of the Universe. As we report in the companion letter, the TA data show large deflections with respect to the LSS which can be explained, assuming small extra-galactic magnetic fields (EGMF), by an intermediate composition changing to a heavy one (iron) in the highest energy bin. Here we show that these results are robust to uncertainties in UHECR injection spectra, the energy scale of the experiment and galactic magnetic fields (GMF). The assumption of weak EGMF, however, strongly affects this interpretation at all but the highest energies E > 100 EeV, where the remarkable isotropy of the data implies a heavy injected composition even in the case of strong EGMF. This result also holds if UHECR sources are as rare as $2 \times 10^{-5}$ Mpc$^{-3}$, that is the conservative lower limit for the source number density. △ Less

Submitted 3 July, 2024; v1 submitted 27 June, 2024; originally announced June 2024.

Comments: 18 pages, 11 figures, accepted for publication in PRD

Observation of Declination Dependence in the Cosmic Ray Energy Spectrum

Authors: The Telescope Array Collaboration, R. U. Abbasi, T. Abu-Zayyad, M. Allen, J. W. Belz, D. R. Bergman, I. Buckland, W. Campbell, B. G. Cheon, K. Endo, A. Fedynitch, T. Fujii, K. Fujisue, K. Fujita, M. Fukushima, G. Furlich, Z. Gerber, N. Globus, W. Hanlon, N. Hayashida, H. He, K. Hibino, R. Higuchi, D. Ikeda, T. Ishii , et al. (101 additional authors not shown)

Abstract: We report on an observation of the difference between northern and southern skies of the ultrahigh energy cosmic ray energy spectrum with a significance of ${\sim}8σ$. We use measurements from the two largest experiments$\unicode{x2014}$the Telescope Array observing the northern hemisphere and the Pierre Auger Observatory viewing the southern hemisphere. Since the comparison of two measurements fr… ▽ More We report on an observation of the difference between northern and southern skies of the ultrahigh energy cosmic ray energy spectrum with a significance of ${\sim}8σ$. We use measurements from the two largest experiments$\unicode{x2014}$the Telescope Array observing the northern hemisphere and the Pierre Auger Observatory viewing the southern hemisphere. Since the comparison of two measurements from different observatories introduces the issue of possible systematic differences between detectors and analyses, we validate the methodology of the comparison by examining the region of the sky where the apertures of the two observatories overlap. Although the spectra differ in this region, we find that there is only a $1.8σ$ difference between the spectrum measurements when anisotropic regions are removed and a fiducial cut in the aperture is applied. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 8 pages, 6 figures

Probing Shocked Ejecta in SN 1987A with XRISM-Resolve: the effects of the gate valve closed

Authors: Vincenzo Sapienza, Marco Miceli, Aya Bamba, Salvatore Orlando, Shiu-Hang Lee, Shigehiro Nagataki, Masaomi Ono, Satoru Katsuda, Koji Mori, Makoto Sawada, Yukikatsu Terada, Roberta Giuffrida, Fabrizio Bocchino

Abstract: Supernova (SN) 1987A is widely regarded as an excellent candidate for leveraging the capabilities of the freshly launched XRISM satellite. Recent researches indicate that the X-ray emission from SN 1987A will increasingly originate from its ejecta in the years to come. In a previous study, we thoroughly examined the proficiency of XRISM-Resolve in identifying signatures of shocked ejecta in SN 198… ▽ More Supernova (SN) 1987A is widely regarded as an excellent candidate for leveraging the capabilities of the freshly launched XRISM satellite. Recent researches indicate that the X-ray emission from SN 1987A will increasingly originate from its ejecta in the years to come. In a previous study, we thoroughly examined the proficiency of XRISM-Resolve in identifying signatures of shocked ejecta in SN 1987A, synthesizing the XRISM-Resolve spectrum based on a state-of-the-art magneto-hydrodynamic simulation. However, following the satellite's launch, a technical issue arose with the XRISM instrument's gate valve, which failed to open, thereby affecting observations with the Resolve spectrometer. Here, we update our analysis, reevaluating our diagnostic approach under the assumption that the gate valve remains closed. We find that, even with the reduced instrumental capabilities, it will be possible to pinpoint the ejecta contribution through the study of the line profiles in the XRISM-Resolve spectrum of SN 1987A. △ Less

Submitted 31 May, 2024; originally announced June 2024.

Comments: 4 pages, 1 figure. Updates of the results described in Sapienza et al. 2024 (http://doi.org/10.3847/2041-8213/ad16e3). Submitted for publication in Research Notes of the AAS

EUSO-SPB1 Mission and Science

Authors: JEM-EUSO Collaboration, :, G. Abdellaoui, S. Abe, J. H. Adams. Jr., D. Allard, G. Alonso, L. Anchordoqui, A. Anzalone, E. Arnone, K. Asano, R. Attallah, H. Attoui, M. Ave Pernas, R. Bachmann, S. Bacholle, M. Bagheri, M. Bakiri, J. Baláz, D. Barghini, S. Bartocci, M. Battisti, J. Bayer, B. Beldjilali, T. Belenguer , et al. (271 additional authors not shown)

Abstract: The Extreme Universe Space Observatory on a Super Pressure Balloon 1 (EUSO-SPB1) was launched in 2017 April from Wanaka, New Zealand. The plan of this mission of opportunity on a NASA super pressure balloon test flight was to circle the southern hemisphere. The primary scientific goal was to make the first observations of ultra-high-energy cosmic-ray extensive air showers (EASs) by looking down on… ▽ More The Extreme Universe Space Observatory on a Super Pressure Balloon 1 (EUSO-SPB1) was launched in 2017 April from Wanaka, New Zealand. The plan of this mission of opportunity on a NASA super pressure balloon test flight was to circle the southern hemisphere. The primary scientific goal was to make the first observations of ultra-high-energy cosmic-ray extensive air showers (EASs) by looking down on the atmosphere with an ultraviolet (UV) fluorescence telescope from suborbital altitude (33~km). After 12~days and 4~hours aloft, the flight was terminated prematurely in the Pacific Ocean. Before the flight, the instrument was tested extensively in the West Desert of Utah, USA, with UV point sources and lasers. The test results indicated that the instrument had sensitivity to EASs of approximately 3 EeV. Simulations of the telescope system, telescope on time, and realized flight trajectory predicted an observation of about 1 event assuming clear sky conditions. The effects of high clouds were estimated to reduce this value by approximately a factor of 2. A manual search and a machine-learning-based search did not find any EAS signals in these data. Here we review the EUSO-SPB1 instrument and flight and the EAS search. △ Less

Submitted 12 January, 2024; originally announced January 2024.

Comments: 18 pages, 19 figures

Journal ref: Astropart Phys 154 (2024) 102891

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

Authors: Maria Manuela Saez, Ermal Rrapaj, Akira Harada, Shigehiro Nagataki, Yong-Zhong Qian

Abstract: 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… ▽ More 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. △ Less

Submitted 8 January, 2024; v1 submitted 4 January, 2024; originally announced January 2024.

Comments: 31 pages, 12 Figures

Report number: RIKEN-iTHEMS-Report-24

Non-Hermitian $p$-wave superfluid and effects of the inelastic three-body loss in a one-dimensional spin-polarized Fermi gas

Authors: Hiroyuki Tajima, Yuta Sekino, Daisuke Inotani, Akira Dohi, Shigehiro Nagataki, Tomoya Hayata

Abstract: We theoretically investigate non-Hermitian $p$-wave Fermi superfluidity in one-dimensional spin-polarized Fermi gases which is relevant to recent ultracold atomic experiments. Considering an imaginary atom-dimer coupling responsible for the three-body recombination process in the Lindblad formalism, we discuss the stability of the superfluid state against the atomic loss effect. Within the two-cha… ▽ More We theoretically investigate non-Hermitian $p$-wave Fermi superfluidity in one-dimensional spin-polarized Fermi gases which is relevant to recent ultracold atomic experiments. Considering an imaginary atom-dimer coupling responsible for the three-body recombination process in the Lindblad formalism, we discuss the stability of the superfluid state against the atomic loss effect. Within the two-channel non-Hermitian BCS-Leggett theory, the atomic loss is characterized by the product of the imaginary atom-dimer coupling and the $p$-wave effective range. Our results indicate that for a given imaginary atom-dimer coupling, a smaller magnitude of the effective ranges of $p$-wave interaction is crucial for reaching the non-Hermitian $p$-wave Fermi superfluid state. △ Less

Submitted 25 December, 2023; originally announced December 2023.

Comments: 9 pages, 3 figures

Report number: RIKEN-iTHEMS-Report-23

Probing Shocked Ejecta in SN 1987A: A novel diagnostic approach using XRISM-Resolve

Authors: Vincenzo Sapienza, Marco Miceli, Aya Bamba, Salvatore Orlando, Shiu-Hang Lee, Shigehiro Nagataki, Masaomi Ono, Satoru Katsuda, Koji Mori, Makoto Sawada, Yukikatsu Terada, Roberta Giuffrida, Fabrizio Bocchino

Abstract: Supernova (SN) 1987A is one of the best candidates to exploit the capabilities of the freshly launched XRISM satellite. This celestial object offers the unique opportunity to study the evolution of a SN into a young supernova remnant. To date, the X-ray emission has been dominated by the shocked circumstellar medium (CSM), with no shocked ejecta firmly detected. However, recent studies provide com… ▽ More Supernova (SN) 1987A is one of the best candidates to exploit the capabilities of the freshly launched XRISM satellite. This celestial object offers the unique opportunity to study the evolution of a SN into a young supernova remnant. To date, the X-ray emission has been dominated by the shocked circumstellar medium (CSM), with no shocked ejecta firmly detected. However, recent studies provide compelling evidence that in the forthcoming years the X-ray emission from SN 1987A will increasingly stem from the ejecta. Our aim is to assess the proficiency of XRISM-Resolve high resolution spectrometer in pinpointing signatures of the shocked ejecta in SN 1987A. Taking advantage of a self consistent state-of-art magneto-hydrodynamic simulation that describes the evolution from SN 1987A to its remnant, we synthesized the XRISM-Resolve spectrum of SN 1987A, as it would be collected in the allocated observation during the performance verification phase, which is foreseen for 2024. Our predictions clearly show the leading role of shocked ejecta in shaping the profile of the emission lines. The Doppler broadening associated with the bulk motion along the line of sight of the rapidly expanding ejecta is shown to increase the line widths well above the values observed so far. The quantitative comparison between our synthetic spectra and the XRISM spectra will enable us to establish a strong connection between the broadened line emission and the freshly shocked ejecta. This, in turn, will allow us to retrieve the ejecta dynamics and chemical composition from the X-ray emission. △ Less

Submitted 18 December, 2023; originally announced December 2023.

Comments: 9 pages, 4 figures, 1 table, Accepted for publication in The Astrophysical Journal Letters

JEM-EUSO Collaboration contributions to the 38th International Cosmic Ray Conference

Authors: S. Abe, J. H. Adams Jr., D. Allard, P. Alldredge, R. Aloisio, L. Anchordoqui, A. Anzalone, E. Arnone, M. Bagheri, B. Baret, D. Barghini, M. Battisti, R. Bellotti, A. A. Belov, M. Bertaina, P. F. Bertone, M. Bianciotto, F. Bisconti, C. Blaksley, S. Blin-Bondil, K. Bolmgren, S. Briz, J. Burton, F. Cafagna, G. Cambiè , et al. (133 additional authors not shown)

Abstract: This is a collection of papers presented by the JEM-EUSO Collaboration at the 38th International Cosmic Ray Conference (Nagoya, Japan, July 26-August 3, 2023) This is a collection of papers presented by the JEM-EUSO Collaboration at the 38th International Cosmic Ray Conference (Nagoya, Japan, July 26-August 3, 2023) △ Less

Submitted 13 December, 2023; originally announced December 2023.

Developments and results in the context of the JEM-EUSO program obtained with the ESAF Simulation and Analysis Framework

Authors: S. Abe, J. H. Adams Jr., D. Allard, P. Alldredge, L. Anchordoqui, A. Anzalone, E. Arnone, B. Baret, D. Barghini, M. Battisti, J. Bayer, R. Bellotti, A. A. Belov, M. Bertaina, P. F. Bertone, M. Bianciotto, P. L. Biermann, F. Bisconti, C. Blaksley, S. Blin-Bondil, P. Bobik, K. Bolmgren, S. Briz, J. Burton, F. Cafagna , et al. (150 additional authors not shown)

Abstract: JEM--EUSO is an international program for the development of space-based Ultra-High Energy Cosmic Ray observatories. The program consists of a series of missions which are either under development or in the data analysis phase. All instruments are based on a wide-field-of-view telescope, which operates in the near-UV range, designed to detect the fluorescence light emitted by extensive air showers… ▽ More JEM--EUSO is an international program for the development of space-based Ultra-High Energy Cosmic Ray observatories. The program consists of a series of missions which are either under development or in the data analysis phase. All instruments are based on a wide-field-of-view telescope, which operates in the near-UV range, designed to detect the fluorescence light emitted by extensive air showers in the atmosphere. We describe the simulation software ESAFin the framework of the JEM--EUSO program and explain the physical assumptions used. We present here the implementation of the JEM--EUSO, POEMMA, K--EUSO, TUS, Mini--EUSO, EUSO--SPB1 and EUSO--TA configurations in ESAF. For the first time ESAF simulation outputs are compared with experimental data. △ Less

Submitted 21 November, 2023; originally announced November 2023.

Journal ref: Eur. Phys. J. C 83, 1028 (2023)

Constraints on the Neutron-Star Structure from the Clocked X-Ray Burster 1RXS J180408.9$-$342058

Authors: Akira Dohi, Wataru Iwakiri, Nobuya Nishimura, Tsuneo Noda, Shigehiro Nagataki, Masa-aki Hashimoto

Abstract: Type-I X-ray bursts are rapid-brightening transient phenomena on the surfaces of accreting neutron stars (NSs). Some X-ray bursts, called {\it clocked bursters}, exhibit regular behavior with similar light curve profiles in their burst sequences. The periodic nature of clocked bursters has the advantage of constraining X-ray binary parameters and physics inside the NS. In the present study, we com… ▽ More Type-I X-ray bursts are rapid-brightening transient phenomena on the surfaces of accreting neutron stars (NSs). Some X-ray bursts, called {\it clocked bursters}, exhibit regular behavior with similar light curve profiles in their burst sequences. The periodic nature of clocked bursters has the advantage of constraining X-ray binary parameters and physics inside the NS. In the present study, we compute numerical models, based on different equations of state and NS masses, which are compared with the observation of a recently identified clocked burster, 1RXS J180408.9$-$342058. We find that the relation between accretion rate and recurrence time is highly sensitive to the NS mass and radius. We determine, in particular, that 1RXS J180408.9$-$342058 appears to possess a mass less than $1.7M_{\odot}$ and favors a stiffer nuclear equation of state (with an NS radius $\gtrsim12.7{\rm km}$). Consequently, the observations of this new clocked burster may provide additional constraints for probing the structure of NSs. △ Less

Submitted 5 November, 2023; originally announced November 2023.

Comments: 7 pages, 5 figures, accepted for publication in ApJ

Report number: RIKEN-iTHEMS-Report-23

Potential of the Cherenkov Telescope Array for studying the young nearby supernova remnants

Authors: Dmitry Prokhorov, Shiu-Hang Lee, Shigehiro Nagataki, Jacco Vink, Donald Ellison, Gilles Ferrand, Daniel Patnaude, Friedrich Roepke, Ivo Seitenzahlh, Patrick Slane

Abstract: Modern imaging atmospheric Cherenkov telescopes have extensively observed young nearby supernova remnants (SNRs), with ages of about 1000 years or less, in the very-high-energy (VHE) gamma-ray band. These efforts resulted in the detection of VHE emission from three young SNRs - Cassiopeia A, Tycho, and SN 1006 - and provided significant evidence for emission from the more distant Kepler's SNR. How… ▽ More Modern imaging atmospheric Cherenkov telescopes have extensively observed young nearby supernova remnants (SNRs), with ages of about 1000 years or less, in the very-high-energy (VHE) gamma-ray band. These efforts resulted in the detection of VHE emission from three young SNRs - Cassiopeia A, Tycho, and SN 1006 - and provided significant evidence for emission from the more distant Kepler's SNR. However, many questions on the production of VHE gamma rays in these remnants remain unanswered. Using detailed physical models for Tycho's SNR based on the CR-hydro-NEI code and physically motivated models for the other young nearby remnants, we simulated observations with the Cherenkov Telescope Array (CTA) of these gamma-ray sources. We highlight properties of these remnants accessible for investigation with future CTA observations and discuss which questions are expected to be answered. △ Less

Submitted 4 October, 2023; originally announced October 2023.

Comments: Presented at the 38th International Cosmic Ray Conference (ICRC 2023), 2023 (arXiv:2309.08219)

Report number: CTA-ICRC/2023/23

Performance of the joint LST-1 and MAGIC observations evaluated with Crab Nebula data

Authors: H. Abe, K. Abe, S. Abe, V. A. Acciari, A. Aguasca-Cabot, I. Agudo, N. Alvarez Crespo, T. Aniello, S. Ansoldi, L. A. Antonelli, C. Aramo, A. Arbet-Engels, C. Arcaro, M. Artero, K. Asano, P. Aubert, D. Baack, A. Babić, A. Baktash, A. Bamba, A. Baquero Larriva, L. Baroncelli, U. Barres de Almeida, J. A. Barrio, I. Batković , et al. (344 additional authors not shown)

Abstract: Aims. LST-1, the prototype of the Large-Sized Telescope for the upcoming Cherenkov Telescope Array Observatory, is concluding its commissioning in Observatorio del Roque de los Muchachos on the island of La Palma. The proximity of LST-1 (Large-Sized Telescope 1) to the two MAGIC (Major Atmospheric Gamma Imaging Cherenkov) telescopes permits observations of the same gamma-ray events with both syste… ▽ More Aims. LST-1, the prototype of the Large-Sized Telescope for the upcoming Cherenkov Telescope Array Observatory, is concluding its commissioning in Observatorio del Roque de los Muchachos on the island of La Palma. The proximity of LST-1 (Large-Sized Telescope 1) to the two MAGIC (Major Atmospheric Gamma Imaging Cherenkov) telescopes permits observations of the same gamma-ray events with both systems. Methods. We describe the joint LST-1+MAGIC analysis pipeline and use simultaneous Crab Nebula observations and Monte Carlo simulations to assess the performance of the three-telescope system. The addition of the LST-1 telescope allows the recovery of events in which one of the MAGIC images is too dim to survive analysis quality cuts. Results. Thanks to the resulting increase in the collection area and stronger background rejection, we find a significant improvement in sensitivity, allowing the detection of 30% weaker fluxes in the energy range between 200 GeV and 3 TeV. The spectrum of the Crab Nebula, reconstructed in the energy range ~60 GeV to ~10 TeV, is in agreement with previous measurements. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 680, A66 (2023)

Prospects for $γ$-ray observations of the Perseus galaxy cluster with the Cherenkov Telescope Array

Authors: The Cherenkov Telescope Array Consortium, :, K. Abe, S. Abe, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Alfaro, N. Alvarez-Crespo, R. Alves Batista, J. -P. Amans, E. Amato, E. O. Angüner, L. A. Antonelli, C. Aramo, M. Araya, C. Arcaro, L. Arrabito, K. Asano, Y. Ascasíbar, J. Aschersleben , et al. (542 additional authors not shown)

Abstract: Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster med… ▽ More Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster medium. We estimate the sensitivity of the Cherenkov Telescope Array (CTA) to detect diffuse gamma-ray emission from the Perseus galaxy cluster. We perform a detailed spatial and spectral modelling of the expected signal for the DM and the CRp components. For each, we compute the expected CTA sensitivity. The observing strategy of Perseus is also discussed. In the absence of a diffuse signal (non-detection), CTA should constrain the CRp to thermal energy ratio within the radius $R_{500}$ down to about $X_{500}<3\times 10^{-3}$, for a spatial CRp distribution that follows the thermal gas and a CRp spectral index $α_{\rm CRp}=2.3$. Under the optimistic assumption of a pure hadronic origin of the Perseus radio mini-halo and depending on the assumed magnetic field profile, CTA should measure $α_{\rm CRp}$ down to about $Δα_{\rm CRp}\simeq 0.1$ and the CRp spatial distribution with 10% precision. Regarding DM, CTA should improve the current ground-based gamma-ray DM limits from clusters observations on the velocity-averaged annihilation cross-section by a factor of up to $\sim 5$, depending on the modelling of DM halo substructure. In the case of decay of DM particles, CTA will explore a new region of the parameter space, reaching models with $τ_χ>10^{27}$s for DM masses above 1 TeV. These constraints will provide unprecedented sensitivity to the physics of both CRp acceleration and transport at cluster scale and to TeV DM particle models, especially in the decay scenario. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Comments: 93 pages (including author list, appendix and references), 143 figures. Submitted to JCAP

Numerical simulation of photospheric emission in long gamma-ray bursts: prompt correlations, spectral shapes, and polarizations

Authors: Hirotaka Ito, Jin Matsumoto, Shigehiro Nagataki, Donald C. Warren, Maxim V. Barkov, Daisuke Yonetoku

Abstract: We explore the properties of photospheric emission in the context of long gamma-ray bursts (LGRBs) using three numerical models that combine relativistic hydrodynamical simulations and Monte Carlo radiation transfer calculations in three dimensions. Our simulations confirm that the photospheric emission gives rise to correlations between the spectral peak energy and luminosity that agree with the… ▽ More We explore the properties of photospheric emission in the context of long gamma-ray bursts (LGRBs) using three numerical models that combine relativistic hydrodynamical simulations and Monte Carlo radiation transfer calculations in three dimensions. Our simulations confirm that the photospheric emission gives rise to correlations between the spectral peak energy and luminosity that agree with the observed Yonetoku, Amati, and Golenetskii correlations. It is also shown that the spectral peak energy and luminosity correlate with the bulk Lorentz factor, as indicated in the literature. On the other hand, synthetic spectral shapes tend to be narrower than those of the observations. The result indicates that an additional physical process that can provide non-thermal broadening is needed to reproduce the spectral features. Furthermore, the polarization analysis finds that, while the degree of polarization is low for the emission from the jet core ($Π< 4~\%$), it tends to increase with the viewing angle outside the core and can be as high as $Π\sim 20-40~\%$ in an extreme case. This suggests that the typical GRBs show systematically low polarization compared to softer, dimmer counterparts (X-ray-rich GRBs and X-ray flashes). Interestingly, our simulations indicate that photospheric emission exhibits large temporal variation in the polarization position angle ($Δψ\sim 90^{\circ}$), which may be compatible with those inferred in observations. A notable energy dependence of the polarization property is another characteristic feature found in the current study. Particularly, the difference in the position angle among different energy bands can be as large as $\sim 90^{\circ}$. △ Less

Submitted 19 July, 2023; originally announced July 2023.

Comments: 20 pages, 12 figures, accepted for publication in ApJ

Report number: RIKEN-iTHEMS-Report-23

Examining Neutrino-Matter Interactions in the Cassiopeia A Supernova

Authors: Toshiki Sato, Takashi Yoshida, Hideyuki Umeda, John P. Hughes, Keiichi Maeda, Shigehiro Nagataki, Brian J. Williams

Abstract: Neutrino interactions with stellar material are widely believed to be fundamental to the explosion of massive stars. However, this important process has remained difficult to confirm observationally. We propose a new method to verify it using X-ray observations of the supernova remnant Cassiopeia A. The elemental composition in its Fe-rich ejecta that could have been produced at the innermost regi… ▽ More Neutrino interactions with stellar material are widely believed to be fundamental to the explosion of massive stars. However, this important process has remained difficult to confirm observationally. We propose a new method to verify it using X-ray observations of the supernova remnant Cassiopeia A. The elemental composition in its Fe-rich ejecta that could have been produced at the innermost region of the supernova, where neutrinos are expected to interact, allows us to examine the presence of neutrino interactions. Here we demonstrate that the amount of Mn produced without neutrino nucleosynthesis processes (i.e., the $ν$- and $ν$p-process) is too small to explain the Mn/Fe mass ratio we measure (0.14--0.67\%). This result supports the operation of significant neutrino interactions in the Cassiopeia A supernova. If the observed Mn/Fe mass ratio purely reflects the production at the innermost region of the supernova, this would be the first robust confirmation of neutrino-matter interactions in an individual supernova. We further show that the Mn/Fe mass ratio has the potential to constrain supernova neutrino parameters (i.e., total neutrino luminosity, neutrino temperature). Future spatially-resolved, high-resolution X-ray spectroscopy will allow us to investigate the details of neutrino-supernova astrophysics through its signatures in elemental composition not only in Cassiopeia A but also in other remnants. △ Less

Submitted 7 July, 2023; originally announced July 2023.

Comments: 17 pages, 10 figures. Accepted for publication in ApJ

Observations of the Crab Nebula and Pulsar with the Large-Sized Telescope Prototype of the Cherenkov Telescope Array

Authors: CTA-LST Project, :, H. Abe, K. Abe, S. Abe, A. Aguasca-Cabot, I. Agudo, N. Alvarez Crespo, L. A. Antonelli, C. Aramo, A. Arbet-Engels, C. Arcaro, M. Artero, K. Asano, P. Aubert, A. Baktash, A. Bamba, A. Baquero Larriva, L. Baroncelli, U. Barres de Almeida, J. A. Barrio, I. Batkovic, J. Baxter, J. Becerra González, E. Bernardini , et al. (467 additional authors not shown)

Abstract: CTA (Cherenkov Telescope Array) is the next generation ground-based observatory for gamma-ray astronomy at very-high energies. The Large-Sized Telescope prototype (LST-1) is located at the Northern site of CTA, on the Canary Island of La Palma. LSTs are designed to provide optimal performance in the lowest part of the energy range covered by CTA, down to $\simeq 20$ GeV. LST-1 started performing a… ▽ More CTA (Cherenkov Telescope Array) is the next generation ground-based observatory for gamma-ray astronomy at very-high energies. The Large-Sized Telescope prototype (LST-1) is located at the Northern site of CTA, on the Canary Island of La Palma. LSTs are designed to provide optimal performance in the lowest part of the energy range covered by CTA, down to $\simeq 20$ GeV. LST-1 started performing astronomical observations in November 2019, during its commissioning phase, and it has been taking data since then. We present the first LST-1 observations of the Crab Nebula, the standard candle of very-high energy gamma-ray astronomy, and use them, together with simulations, to assess the basic performance parameters of the telescope. The data sample consists of around 36 hours of observations at low zenith angles collected between November 2020 and March 2022. LST-1 has reached the expected performance during its commissioning period - only a minor adjustment of the preexisting simulations was needed to match the telescope behavior. The energy threshold at trigger level is estimated to be around 20 GeV, rising to $\simeq 30$ GeV after data analysis. Performance parameters depend strongly on energy, and on the strength of the gamma-ray selection cuts in the analysis: angular resolution ranges from 0.12 to 0.40 degrees, and energy resolution from 15 to 50%. Flux sensitivity is around 1.1% of the Crab Nebula flux above 250 GeV for a 50-h observation (12% for 30 minutes). The spectral energy distribution (in the 0.03 - 30 TeV range) and the light curve obtained for the Crab Nebula agree with previous measurements, considering statistical and systematic uncertainties. A clear periodic signal is also detected from the pulsar at the center of the Nebula. △ Less

Submitted 19 July, 2023; v1 submitted 22 June, 2023; originally announced June 2023.

Comments: Accepted in ApJ. v3: updated author list and acknowledgements, fixed typos and other minor issues

Quasars: standard candles up to z=7.5 with the precision of Supernovae Ia

Authors: Maria Giovanna Dainotti, Giada Bargiacchi, Aleksander Łukasz Lenart, Shigehiro Nagataki, Salvatore Capozziello

Abstract: Currently, the $Λ$ Cold Dark Matter model, which relies on the existence of cold dark matter and a cosmological constant $Λ$, best describes the Universe. However, we lack information in the high-redshift ($z$) region between Type Ia Supernovae (SNe Ia) (up to $z=2.26$) and the Cosmic Microwave Background ($z=1100$), an interval crucial to test cosmological models and their possible evolution. We… ▽ More Currently, the $Λ$ Cold Dark Matter model, which relies on the existence of cold dark matter and a cosmological constant $Λ$, best describes the Universe. However, we lack information in the high-redshift ($z$) region between Type Ia Supernovae (SNe Ia) (up to $z=2.26$) and the Cosmic Microwave Background ($z=1100$), an interval crucial to test cosmological models and their possible evolution. We have defined a sample of 983 Quasars up to $z=7.54$ with reduced intrinsic dispersion $δ=0.007$ which determines the matter density parameter $Ω_M$ with the same precision of SNe Ia. Although previous analysis have been used Quasars as cosmological tools (e.g. Risaliti and Lusso 2019), this is the first time that high-redshift sources, in this case Quasars, as standalone cosmological probes yield such tight constraints on $Ω_M$. Our results show the importance of correcting cosmological relationships for selection biases and redshift evolution and how the choice of a golden sample reduces considerably the intrinsic scatter. This proves the reliability of Quasars as standard cosmological candles. △ Less

Submitted 31 May, 2023; originally announced May 2023.

Comments: Accepted for publication in ApJ

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

Authors: Masaomi Ono, Takaya Nozawa, Shigehiro Nagataki, Alexandra Kozyreva, Salvatore Orlando, Marco Miceli, Ke-Jung Chen

Abstract: 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 refl… ▽ More 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. △ Less

Submitted 2 January, 2024; v1 submitted 4 May, 2023; originally announced May 2023.

Comments: 85 pages, 32 figures, 8 tables; accepted for publication in ApJS

Report number: RIKEN-iTHEMS-Report-24

Investigating Time Evolution of Thermal Emission from the Putative Neutron Star in SN 1987A for 50+ Years

Authors: Akira Dohi, Emanuele Greco, Shigehiro Nagataki, Masaomi Ono, Marco Miceli, Salvatore Orlando, Barbara Olmi

Abstract: Observations collected with the Atacama Large Millimeter/submillimeter Array (ALMA) and analysis of broadband X-ray spectra have recently suggested the presence of a central compact object (CCO) in SN 1987A. However, no direct evidence of the CCO has been found yet. Here we analyze Chandra X-ray observations of SN 1987A collected in 2007 and 2018, and synthesize the 2027 Chandra and 2037 Lynx spec… ▽ More Observations collected with the Atacama Large Millimeter/submillimeter Array (ALMA) and analysis of broadband X-ray spectra have recently suggested the presence of a central compact object (CCO) in SN 1987A. However, no direct evidence of the CCO has been found yet. Here we analyze Chandra X-ray observations of SN 1987A collected in 2007 and 2018, and synthesize the 2027 Chandra and 2037 Lynx spectra of the faint inner region of SN 1987A. We estimate the temporal evolution of the upper limits of the intrinsic luminosity of the putative CCO in three epochs (2018, 2027 and 2037). We find that these upper limits are higher for higher neutron star (NS) kick velocities due to the increased absorption from the surrounding cold ejecta. We compare NS cooling models with both the intrinsic luminosity limits obtained from the X-ray spectra, and the ALMA constraints with the assumption that the observed blob of SN 1987A is primarily heated by thermal emission. We find that the synthetic Lynx spectra are crucial to constrain physical properties of the CCO, which will be confirmed by future observations in the 2040s. We draw our conclusions based on two scenarios, namely the non-detection and detection of NS by Lynx. If the NS is not detected, its kick velocity should be ~700 km/s. Furthermore, the non-detection of the NS would suggest rapid cooling processes around the age of 40 years, implying strong crust superfluidity. Conversely, in the case of NS detection, the mass of the NS envelope must be high. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 37 pages, 17 figures, 6 tables, accepted for publication in ApJ

Report number: RIKEN-iTHEMS-Report-23

Sensitivity of the Cherenkov Telescope Array to spectral signatures of hadronic PeVatrons with application to Galactic Supernova Remnants

Authors: The Cherenkov Telescope Array Consortium, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Aloisio, N. Álvarez Crespo, R. Alves Batista, L. Amati, E. Amato, G. Ambrosi, E. O. Angüner, C. Aramo, C. Arcaro, T. Armstrong, K. Asano, Y. Ascasibar, J. Aschersleben, M. Backes, A. Baktash, C. Balazs, M. Balbo , et al. (334 additional authors not shown)

Abstract: The local Cosmic Ray (CR) energy spectrum exhibits a spectral softening at energies around 3~PeV. Sources which are capable of accelerating hadrons to such energies are called hadronic PeVatrons. However, hadronic PeVatrons have not yet been firmly identified within the Galaxy. Several source classes, including Galactic Supernova Remnants (SNRs), have been proposed as PeVatron candidates. The pote… ▽ More The local Cosmic Ray (CR) energy spectrum exhibits a spectral softening at energies around 3~PeV. Sources which are capable of accelerating hadrons to such energies are called hadronic PeVatrons. However, hadronic PeVatrons have not yet been firmly identified within the Galaxy. Several source classes, including Galactic Supernova Remnants (SNRs), have been proposed as PeVatron candidates. The potential to search for hadronic PeVatrons with the Cherenkov Telescope Array (CTA) is assessed. The focus is on the usage of very high energy $γ$-ray spectral signatures for the identification of PeVatrons. Assuming that SNRs can accelerate CRs up to knee energies, the number of Galactic SNRs which can be identified as PeVatrons with CTA is estimated within a model for the evolution of SNRs. Additionally, the potential of a follow-up observation strategy under moonlight conditions for PeVatron searches is investigated. Statistical methods for the identification of PeVatrons are introduced, and realistic Monte--Carlo simulations of the response of the CTA observatory to the emission spectra from hadronic PeVatrons are performed. Based on simulations of a simplified model for the evolution for SNRs, the detection of a $γ$-ray signal from in average 9 Galactic PeVatron SNRs is expected to result from the scan of the Galactic plane with CTA after 10 hours of exposure. CTA is also shown to have excellent potential to confirm these sources as PeVatrons in deep observations with $\mathcal{O}(100)$ hours of exposure per source. △ Less

Submitted 27 March, 2023; originally announced March 2023.

Comments: 34 pages, 16 figures, Accepted for publication in Astroparticle Physics

Gamma-Ray Bursts, Quasars, Baryonic Acoustic Oscillations, and Supernovae Ia: new statistical insights and cosmological constraints

Authors: G. Bargiacchi, M. G. Dainotti, S. Nagataki, S. Capozziello

Abstract: The recent $\sim 4 \, σ$ Hubble constant, $H_{0}$, tension is observed between the value of $H_{0}$ from the Cosmic Microwave Background (CMB) and Type Ia Supernovae (SNe Ia). It is a decade since this tension is excruciating the modern astrophysical community. To shed light on this problem is key to consider probes at intermediate redshifts between SNe Ia and CMB and reduce the uncertainty on… ▽ More The recent $\sim 4 \, σ$ Hubble constant, $H_{0}$, tension is observed between the value of $H_{0}$ from the Cosmic Microwave Background (CMB) and Type Ia Supernovae (SNe Ia). It is a decade since this tension is excruciating the modern astrophysical community. To shed light on this problem is key to consider probes at intermediate redshifts between SNe Ia and CMB and reduce the uncertainty on $H_0$. Toward these goals, we fill the redshift gap by employing Gamma-Ray Bursts (GRBs) and Quasars (QSOs), reaching $z=9.4$ and $z=7.6$, respectively, combined with Baryonic Acoustic Oscillations (BAO) and SNe Ia. To this end, we employ the ``Dainotti GRB 3D relation" among the rest-frame end time of the X-ray plateau emission, its corresponding luminosity, and the peak prompt luminosity, and the ``Risaliti-Lusso" QSO relation between ultraviolet and X-ray luminosities. We inquire the commonly adopted Gaussianity assumption on GRBs, QSOs, and BAO. With the joint sample, we fit the flat $Λ$ Cold Dark Matter model with both the Gaussian and the newly discovered likelihoods. We also investigate the impact of the calibration assumed for \textit{Pantheon} and \textit{Pantheon +} SNe Ia on this analysis. Remarkably, we show that only GRBs fulfill the Gaussianity assumption. We achieve small uncertainties on the matter density parameter $Ω_M$ and $H_0$. We find $H_0$ values compatible within 2 $σ$ with the one from the Tip of the Red Giant Branch. Finally, we show that the cosmological results are heavily biased against the arbitrary calibration choice for SNe Ia. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: 17 pages, 6 figures, 1 table, accepted in MNRAS

Reduced uncertainties up to 43\% on the Hubble constant and the matter density with the SNe Ia with a new statistical analysis

Authors: Maria Giovanna Dainotti, Giada Bargiacchi, Malgorzata Bogdan, Salvatore Capozziello, Shigehiro Nagataki

Abstract: Type Ia Supernovae (SNe Ia) are considered the most reliable \textit{standard candles} and they have played an invaluable role in cosmology since the discovery of the Universe's accelerated expansion. During the last decades, the SNe Ia samples have been improved in number, redshift coverage, calibration methodology, and systematics treatment. These efforts led to the most recent \textit{``Pantheo… ▽ More Type Ia Supernovae (SNe Ia) are considered the most reliable \textit{standard candles} and they have played an invaluable role in cosmology since the discovery of the Universe's accelerated expansion. During the last decades, the SNe Ia samples have been improved in number, redshift coverage, calibration methodology, and systematics treatment. These efforts led to the most recent \textit{``Pantheon"} (2018) and \textit{``Pantheon +"} (2022) releases, which enable to constrain cosmological parameters more precisely than previous samples. In this era of precision cosmology, the community strives to find new ways to reduce uncertainties on cosmological parameters. To this end, we start our investigation even from the likelihood assumption of Gaussianity, implicitly used in this domain. Indeed, the usual practise involves constraining parameters through a Gaussian distance moduli likelihood. This method relies on the implicit assumption that the difference between the distance moduli measured and the ones expected from the cosmological model is Gaussianly distributed. In this work, we test this hypothesis for both the \textit{Pantheon} and \textit{Pantheon +} releases. We find that in both cases this requirement is not fulfilled and the actual underlying distributions are a logistic and a Student's t distribution for the \textit{Pantheon} and \textit{Pantheon +} data, respectively. When we apply these new likelihoods fitting a flat $Λ$CDM model, we significantly reduce the uncertainties on $Ω_M$ and $H_0$ of $\sim 40 \%$. This boosts the SNe Ia power in constraining cosmological parameters, thus representing a huge step forward to shed light on the current debated tensions in cosmology. △ Less

Submitted 10 January, 2024; v1 submitted 13 March, 2023; originally announced March 2023.

Comments: 15 pages with 7 figures (with multiple panels), 1 table. This work includes the suggestions of scientists from the community who reached out. The paper is now to be published in JHEAP

Non-Hermitian topological Fermi superfluid near the $p$-wave unitary limit

Authors: Hiroyuki Tajima, Yuta Sekino, Daisuke Inotani, Akira Dohi, Shigehiro Nagataki, Tomoya Hayata

Abstract: We theoretically discuss the non-Hermitian superfluid phase transition in one-dimensional two-component Fermi gases near the $p$-wave Feshbach resonance accompanied by the two-body loss associated with the dipolar relaxation. For the first time we point out that this system gives us an opportunity to explore the interplay among various non-trivial properties such as universal thermodynamics at div… ▽ More We theoretically discuss the non-Hermitian superfluid phase transition in one-dimensional two-component Fermi gases near the $p$-wave Feshbach resonance accompanied by the two-body loss associated with the dipolar relaxation. For the first time we point out that this system gives us an opportunity to explore the interplay among various non-trivial properties such as universal thermodynamics at divergent $p$-wave scattering length, topological phase transition at vanishing chemical potential, and non-Hermitian Bardeen-Cooper-Schrieffer(BCS) to Bose-Einstein condensate (BEC) transition, in a unified manner. In the BCS phase, the loss-induced superfluid-normal transition occurs when the exceptional point appears in the effective non-Hermitian Hamiltonian. In the BEC phase, the diffusive gapless mode can be regarded as a precursor of the instability of the superfluid state. Moreover, we show that the superfluid state is fragile against the two-body loss near the topological phase transition point. △ Less

Submitted 22 December, 2022; originally announced December 2022.

Comments: 7 pages, 5 figures

Report number: RIKEN-iTHEMS-Report-22

Journal ref: Phys. Rev. A 107, 033331 (2023)

Polarized radio emission unveils the structure of the pre-supernova circumstellar magnetic field and the radio emission in SN1987A

Authors: O. Petruk, V. Beshley, S. Orlando, F. Bocchino, M. Miceli, S. Nagataki, M. Ono, S. Loru, A. Pellizzoni, E. Egron

Abstract: The detected polarized radio emission from remnant of SN1987A opens the possibility to unveil the structure of the pre-supernova magnetic field in the circumstellar medium. Properties derived from direct measurements would be of importance for understanding the progenitor stars and their magnetic fields. As the first step to this goal, we adopted the hydrodynamic data from an elaborated three-dime… ▽ More The detected polarized radio emission from remnant of SN1987A opens the possibility to unveil the structure of the pre-supernova magnetic field in the circumstellar medium. Properties derived from direct measurements would be of importance for understanding the progenitor stars and their magnetic fields. As the first step to this goal, we adopted the hydrodynamic data from an elaborated three-dimensional (3-D) numerical model of SN1987A. We have developed an approximate method for `reconstruction' of 3-D magnetic field structure inside supernova remnant on the `hydrodynamic background'. This method uses the distribution of the magnetic field around the progenitor as the initial condition. With such a 3-D magneto-hydrodynamic model, we have synthesized the polarization maps for a number of SN1987A models and compared them to the observations. In this way, we have tested different initial configurations of the magnetic field as well as a structure of the synchrotron emission in SN987A. We have recovered the observed polarization pattern and we have found that the radial component of the ambient pre-supernova magnetic field should be dominant on the length-scale of the present-day radius of SN1987A. The physical reasons for such a field are discussed. △ Less

Submitted 1 December, 2022; originally announced December 2022.

Comments: 13 pages, 8 figures, accepted in MNRAS

A bias-free cosmological analysis with quasars alleviating $H_0$ tension

Authors: Aleksander Łukasz Lenart, Giada Bargiacchi, Maria Giovanna Dainotti, Shigehiro Nagataki, Salvatore Capozziello

Abstract: Cosmological models and their parameters are widely debated because of theoretical and observational mismatches of the standard cosmological model, especially the current discrepancy between the value of the Hubble constant, $H_{0}$, obtained by Type Ia supernovae (SNe Ia), and the Cosmic Microwave Background Radiation (CMB). Thus, considering high-redshift probes like quasars (QSOs), having inter… ▽ More Cosmological models and their parameters are widely debated because of theoretical and observational mismatches of the standard cosmological model, especially the current discrepancy between the value of the Hubble constant, $H_{0}$, obtained by Type Ia supernovae (SNe Ia), and the Cosmic Microwave Background Radiation (CMB). Thus, considering high-redshift probes like quasars (QSOs), having intermediate redshifts between SNe Ia and CMB, is a necessary step. In this work, we use SNe Ia and the most updated QSO sample, reaching redshifts up to $z\sim7.5$, applying the Risaliti-Lusso QSO relation based on a non-linear relation between ultraviolet and X-ray luminosities. We consider this relation both in its original form and corrected for selection biases and evolution in redshift through a reliable statistical method also accounting for the circularity problem. We also explore two approaches: with and without calibration on SNe Ia. We then investigate flat and non-flat standard cosmological models and a flat $w$CDM model, with a constant dark energy equation of state parameter $w$. Remarkably, when correcting for the evolution as a function of cosmology, we obtain closed constraints on $Ω_M$ using only non-calibrated QSOs. We find that considering non-calibrated QSOs combined with SNe Ia and accounting for the same correction, our results are compatible with a flat $Λ$CDM model with $Ω_M = 0.3$ and $H_0 = 70 \, \mathrm{km\,s^{-1}\,Mpc^{-1}}$. Intriguingly, the $H_0$ values obtained place halfway between the one from SNe Ia and CMB, paving the way for new insights into the $H_0$ tension. △ Less

Submitted 28 January, 2023; v1 submitted 19 November, 2022; originally announced November 2022.

Comments: Accepted to be published in ApJS; 24 pages, 10 figures

Report number: RIKEN-iTHEMS-Report-22

PANDORA project: photo-nuclear reactions below $A=60$

Authors: A. Tamii, L. Pellegri, P. -A. Söderström, D. Allard, S. Goriely, T. Inakura, E. Khan, E. Kido, M. Kimura, E. Litvinova, S. Nagataki, P. von Neumann-Cosel, N. Pietralla, N. Shimizu, N. Tsoneva, Y. Utsuno, S. Adachi, P. Adsley, A. Bahini, D. Balabanski, B. Baret, J. A. C. Bekker, S. D. Binda, E. Boicu, A. Bracco , et al. (56 additional authors not shown)

Abstract: Photo-nuclear reactions of light nuclei below a mass of $A=60$ are studied experimentally and theoretically by the PANDORA (Photo-Absorption of Nuclei and Decay Observation for Reactions in Astrophysics) project. Two experimental methods, virtual-photon excitation by proton scattering and real-photo absorption by a high-brilliance gamma-ray beam produced by laser Compton scattering, will be applie… ▽ More Photo-nuclear reactions of light nuclei below a mass of $A=60$ are studied experimentally and theoretically by the PANDORA (Photo-Absorption of Nuclei and Decay Observation for Reactions in Astrophysics) project. Two experimental methods, virtual-photon excitation by proton scattering and real-photo absorption by a high-brilliance gamma-ray beam produced by laser Compton scattering, will be applied to measure the photo-absorption cross sections and the decay branching ratio of each decay channel as a function of the photon energy. Several nuclear models, e.g. anti-symmetrized molecular dynamics, mean-field type models, a large-scale shell model, and ab initio models, will be employed to predict the photo-nuclear reactions. The uncertainty in the model predictions will be evaluated from the discrepancies between the model predictions and the experimental data. The data and the predictions will be implemented in a general reaction calculation code TALYS . The results will be applied to the simulation of the photo-disintegration process of ultra-high-energy cosmic rays in inter-galactic propagation. △ Less

Submitted 18 November, 2022; v1 submitted 7 November, 2022; originally announced November 2022.

The Gamma-ray Bursts fundamental plane correlation as a cosmological tool

Authors: M. G. Dainotti, A. Ł. Lenart, A. Chraya, G. Sarracino, S. Nagataki, N. Fraija, S. Capozziello, M. Bogdan

Abstract: Cosmological models and their corresponding parameters are widely debated because of the current discrepancy between the results of the Hubble constant, $H_{0}$, obtained by SNe Ia, and the Planck data from the Cosmic Microwave Background Radiation. Thus, considering high redshift probes like Gamma-Ray Bursts (GRBs) is a necessary step. However, using GRB correlations between their physical featur… ▽ More Cosmological models and their corresponding parameters are widely debated because of the current discrepancy between the results of the Hubble constant, $H_{0}$, obtained by SNe Ia, and the Planck data from the Cosmic Microwave Background Radiation. Thus, considering high redshift probes like Gamma-Ray Bursts (GRBs) is a necessary step. However, using GRB correlations between their physical features to infer cosmological parameters is difficult because GRB luminosities span several orders of magnitude. In our work, we use a 3-dimensional relation between the peak prompt luminosity, the rest-frame time at the end of the X-ray plateau, and its corresponding luminosity in X-rays: the so-called 3D Dainotti fundamental plane relation. We correct this relation by considering the selection and evolutionary effects with a reliable statistical method, obtaining a lower central value for the intrinsic scatter, $σ_{int}=0.18 \pm 0.07$ (47.1 \%) compared to previous results, when we adopt a particular set of GRBs with well-defined morphological features, called the platinum sample. We have used the GRB fundamental plane relation alone with both Gaussian and uniform priors on cosmological parameters and in combination with SNe Ia and BAO measurements to infer cosmological parameters like $H_{0}$, the matter density in the universe ($Ω_{M}$), and the dark energy parameter $w$ for a $w$CDM model. Our results are consistent with the parameters given by the $Λ$CDM model but with the advantage of using cosmological probes detected up to $z=5$, much larger than the one observed for the furthest SNe Ia. △ Less

Submitted 18 September, 2022; originally announced September 2022.

Comments: 43 pages, 29 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

Impacts of the direct URCA and Superfluidity inside a Neutron Star on Type-I X-Ray Bursts and X-Ray Superbursts

Authors: Akira Dohi, Nobuya Nishimura, Hajime Sotani, Tsuneo Noda, Helei Liu, Shigehiro Nagataki, Masa-aki Hashimoto

Abstract: We investigate the impacts of neutrino cooling mechanism inside the neutron star (NS) core on the light curves of type-I X-ray bursts and X-ray superbursts. From several observations of NS thermal evolution, physical processes of fast neutrino cooling, such as the direct Urca (DU) process, are indicated. They significantly decrease the surface temperature of NSs, though the cooling effect could be… ▽ More We investigate the impacts of neutrino cooling mechanism inside the neutron star (NS) core on the light curves of type-I X-ray bursts and X-ray superbursts. From several observations of NS thermal evolution, physical processes of fast neutrino cooling, such as the direct Urca (DU) process, are indicated. They significantly decrease the surface temperature of NSs, though the cooling effect could be suppressed by nucleon superfluidity. In the present study, focusing on the DU process and nucleon superfluidity, we investigate the effects of NS cooling on the X-ray bursts using a general-relativistic stellar-evolution code. We find that the DU process leads find the longer recurrence time and the higher peak luminosity, which could be obstructed by the neutrons superfluidity. We also apply our burst models to the comparison with {\it Clocked burster} GS 1826$-$24, and to the recurrence time of superburst triggered by carbon ignition. These effects are significant within a certain range of binary parameters and uncertainty of the NS equation of state. △ Less

Submitted 20 March, 2023; v1 submitted 31 August, 2022; originally announced August 2022.

Comments: 12 pages, 8 figures, 3 Tables

Report number: RIKEN-iTHEMS-Report-22

Journal ref: ApJ 937 124 (2022)

Evaluations of uncertainties in simulations of propagation of ultrahigh-energy cosmic-ray nuclei derived from microscopic nuclear models

Authors: E. Kido, T. Inakura, M. Kimura, N. Kobayashi, S. Nagataki, N. Shimizu, A. Tamii, Y. Utsuno

Abstract: Photodisintegration is a main energy loss process for ultrahigh-energy cosmic-ray (UHECR) nuclei in intergalactic space. Therefore, it is crucial to understand systematic uncertainty in photodisintegration when simulating the propagation of UHECR nuclei. In this work, we calculated the cross sections using the random phase approximation (RPA) of density functional theory (DFT), a microscopic nucle… ▽ More Photodisintegration is a main energy loss process for ultrahigh-energy cosmic-ray (UHECR) nuclei in intergalactic space. Therefore, it is crucial to understand systematic uncertainty in photodisintegration when simulating the propagation of UHECR nuclei. In this work, we calculated the cross sections using the random phase approximation (RPA) of density functional theory (DFT), a microscopic nuclear model. We calculated the $E1$ strength of 29 nuclei using three different density functionals. We obtained the cross sections of photonuclear reactions, including photodisintegration, with the $E1$ strength. Then, we implemented the cross sections in the cosmic-ray propagation code CRPropa. We found that assuming certain astrophysical parameter values, the difference between UHECR energy spectrum predictions using the RPA calculation and the default photodisintegration model in CRPropa can be more than the statistical uncertainty of the spectrum. We also found that the differences between the RPA calculations and CRPropa default in certain astrophysical parameters obtained by a combined fit of UHECR energy spectrum and composition data assuming a phenomenological model of UHECR sources can be more than the uncertainty of the data. △ Less

Submitted 24 May, 2023; v1 submitted 7 June, 2022; originally announced June 2022.

Comments: 37 pages, 16 figures

Report number: RIKEN-iTHEMS-Report-22

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

Authors: Oliver Just, Miguel A. Aloy, Martin Obergaulinger, Shigehiro Nagataki

Abstract: 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 des… ▽ More 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. △ Less

Submitted 4 August, 2022; v1 submitted 27 May, 2022; originally announced May 2022.

Comments: 7 pages, 4 figures, 1 table, slightly revised discussion, main results unchanged compared to v1, accepted to ApJL

Report number: RIKEN-iTHEMS-Report-22

First High-speed Video Camera Observations of a Lightning Flash Associated with a Downward Terrestrial Gamma-ray Flash

Authors: R. U. Abbasi, M. M. F. Saba, J. W. Belz, P. R. Krehbiel, W. Rison, N. Kieu, D. R. da Silva, Dan Rodeheffer, M. A. Stanley, J. Remington, J. Mazich, R. LeVon, K. Smout, A. Petrizze, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, D. R. Bergman, S. A. Blake, I. Buckland, B. G. Cheon, M. Chikawa, T. Fujii , et al. (127 additional authors not shown)

Abstract: In this paper, we present the first high-speed video observation of a cloud-to-ground lightning flash and its associated downward-directed Terrestrial Gamma-ray Flash (TGF). The optical emission of the event was observed by a high-speed video camera running at 40,000 frames per second in conjunction with the Telescope Array Surface Detector, Lightning Mapping Array, interferometer, electric-field… ▽ More In this paper, we present the first high-speed video observation of a cloud-to-ground lightning flash and its associated downward-directed Terrestrial Gamma-ray Flash (TGF). The optical emission of the event was observed by a high-speed video camera running at 40,000 frames per second in conjunction with the Telescope Array Surface Detector, Lightning Mapping Array, interferometer, electric-field fast antenna, and the National Lightning Detection Network. The cloud-to-ground flash associated with the observed TGF was formed by a fast downward leader followed by a very intense return stroke peak current of -154 kA. The TGF occurred while the downward leader was below cloud base, and even when it was halfway in its propagation to ground. The suite of gamma-ray and lightning instruments, timing resolution, and source proximity offer us detailed information and therefore a unique look at the TGF phenomena. △ Less

Submitted 9 August, 2023; v1 submitted 10 May, 2022; originally announced May 2022.

Journal ref: Geophysical Research Letters, 50, e2023GL102958 (2023)

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

Authors: Emanuele Greco, Marco Miceli, Salvatore Orlando, Barbara Olmi, Fabrizio Bocchino, Shigehiro Nagataki, Lei Sun, Jacco Vink, Vincenzo Sapienza, Masaomi Ono, Akira Dohi, Giovanni Peres

Abstract: 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 arisi… ▽ More 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. △ Less

Submitted 14 April, 2022; originally announced April 2022.

Comments: Accepted for publication in ApJ

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

Authors: Maria Giovanna Dainotti, Via Nielson, Giuseppe Sarracino, Enrico Rinaldi, Shigehiro Nagataki, Salvatore Capozziello, Oleg Y. Gnedin, Giada Bargiacchi

Abstract: 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… ▽ More 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. △ Less

Submitted 17 June, 2022; v1 submitted 29 March, 2022; originally announced March 2022.

Comments: 31 pages, 17 figures, 10 tables

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 514, Issue 2, August 2022, Pages 1828-1856

Merger and post-merger of binary neutron stars with a quark-hadron crossover equation of state

Authors: Yong-Jia Huang, Luca Baiotti, Toru Kojo, Kentaro Takami, Hajime Sotani, Hajime Togashi, Tetsuo Hatsuda, Shigehiro Nagataki, Yi-Zhong Fan

Abstract: Fully general-relativistic binary-neutron-star (BNS) merger simulations with quark-hadron crossover (QHC) equations of state (EOSs) are studied for the first time. In contrast to EOSs with purely hadronic matter or with a first-order quark-hadron phase transition (1PT), in the transition region QHC EOSs show a peak in sound speed, and thus a stiffening. We study the effects of such stiffening in t… ▽ More Fully general-relativistic binary-neutron-star (BNS) merger simulations with quark-hadron crossover (QHC) equations of state (EOSs) are studied for the first time. In contrast to EOSs with purely hadronic matter or with a first-order quark-hadron phase transition (1PT), in the transition region QHC EOSs show a peak in sound speed, and thus a stiffening. We study the effects of such stiffening in the merger and post-merger gravitational (GW) signals. Through simulations in the binary-mass range $2.5 < M/M_{\odot} < 2.75$, characteristic differences due to different EOSs appear in the frequency of the main peak of the post-merger GW spectrum ($f_2$), extracted through Bayesian inference. In particular, we found that (i) for lower-mass binaries, since the maximum baryon number density ($n_{\rm max}$) after the merger stays below $3\text{--}4$ times the nuclear-matter density ($n_0$), the characteristic stiffening of the QHC models in that density range results in a lower $f_2$ than that computed for the underlying hadronic EOS and thus also than that for EOSs with a 1PT, (ii) for higher-mass binaries, where $n_{\rm max}$ may exceed $4\text{--}5 n_0$ depending on the EOS model, whether $f_2$ in QHC models is higher or lower than that in the underlying hadronic model depends on the height of the sound-speed peak. Comparing the values of $f_2$ for different EOSs and BNS masses gives important clues on how to discriminate different types of quark dynamics in the high-density end of EOSs and is relevant to future kilohertz GW observations with third-generation GW detectors. △ Less

Submitted 17 November, 2022; v1 submitted 9 March, 2022; originally announced March 2022.

Comments: 6+8 pages, 5+6 figures

Report number: RIKEN-iTHEMS-Report-22

Journal ref: Phys. Rev. Lett. 129, 181101 (2022)

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

Authors: Gilles Ferrand, Ataru Tanikawa, Donald C. Warren, Shigehiro Nagataki, Samar Safi-Harb, Anne Decourchelle

Abstract: 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… ▽ More 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. △ Less

Submitted 8 February, 2022; originally announced February 2022.

Comments: Accepted by ApJ. Animated figures will be available on the online version

Report number: RIKEN-iTHEMS-Report-22

Evidence for past interaction with an asymmetric circumstellar shell in the young SNR Cassiopeia A

Authors: S. Orlando, A. Wongwathanarat, H. -T. Janka, M. Miceli, S. Nagataki, M. Ono, F. Bocchino, J. Vink, D. Milisavljevic, D. J. Patnaude, G. Peres

Abstract: Observations of the SNR Cassiopeia A (Cas A) show asymmetries in the reverse shock that cannot be explained by models describing a remnant expanding through a spherically symmetric wind of the progenitor star. We investigate whether a past interaction of Cas A with an asymmetric circumstellar shell can account for the observed asymmetries. We performed 3D MHD simulations that describe the remnant… ▽ More Observations of the SNR Cassiopeia A (Cas A) show asymmetries in the reverse shock that cannot be explained by models describing a remnant expanding through a spherically symmetric wind of the progenitor star. We investigate whether a past interaction of Cas A with an asymmetric circumstellar shell can account for the observed asymmetries. We performed 3D MHD simulations that describe the remnant evolution from the SN to its interaction with a circumstellar shell. The initial conditions are provided by a 3D neutrino-driven SN model whose morphology resembles Cas A. We explored the parameter space of the shell, searching for a set of parameters able to produce reverse shock asymmetries at the age of 350 years analogous to those observed in Cas A. The interaction of the remnant with the shell can produce asymmetries resembling those observed in the reverse shock if the shell was asymmetric with the densest portion in the nearside to the northwest (NW). The reverse shock shows the following asymmetries at the age of Cas A: i) it moves inward in the observer frame in the NW region, while it moves outward in other regions; ii) the geometric center of the reverse shock is offset to the NW from the geometric center of the forward shock; iii) the reverse shock in the NW region has enhanced nonthermal emission because, there, the ejecta enter the reverse shock with a higher velocity (between 4000 and 7000 km/s) than in other regions (below 2000 km/s). The asymmetries observed in the reverse shock of Cas A can be interpreted as signatures of the interaction of the remnant with an asymmetric circumstellar shell that occurred between 180 and 240 years after the SN event. We suggest that the shell was, most likely, the result of a massive eruption from the progenitor star that occurred between $10^4$ and $10^5$ years prior to core-collapse. We estimate a total mass of the shell of the order 2 Msun. △ Less

Submitted 12 August, 2022; v1 submitted 3 February, 2022; originally announced February 2022.

Comments: 20 pages, 11 figures, 2 Tables, 1 interactive graphics. Accepted for publication on A&A, high quality figures will be available on the published paper

Journal ref: A&A 666, A2 (2022)

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

Authors: Wei-Jian Li, Rui Xue, Guang-Bo Long, Ze-Rui Wang, Shigehiro Nagataki, Da-Hai Yan, Jian-Cheng Wang

Abstract: 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) mode… ▽ More 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. △ Less

Submitted 29 January, 2022; originally announced January 2022.

Comments: 6 pages, 1 figure. Accepted for Publication in A&A

Journal ref: A&A 659, A184 (2022)

JEM-EUSO Collaboration contributions to the 37th International Cosmic Ray Conference

Authors: G. Abdellaoui, S. Abe, J. H. Adams Jr., D. Allard, G. Alonso, L. Anchordoqui, A. Anzalone, E. Arnone, K. Asano, R. Attallah, H. Attoui, M. Ave Pernas, M. Bagheri, J. Baláz, M. Bakiri, D. Barghini, S. Bartocci, M. Battisti, J. Bayer, B. Beldjilali, T. Belenguer, N. Belkhalfa, R. Bellotti, A. A. Belov, K. Benmessai , et al. (267 additional authors not shown)

Abstract: Compilation of papers presented by the JEM-EUSO Collaboration at the 37th International Cosmic Ray Conference (ICRC), held on July 12-23, 2021 (online) in Berlin, Germany. Compilation of papers presented by the JEM-EUSO Collaboration at the 37th International Cosmic Ray Conference (ICRC), held on July 12-23, 2021 (online) in Berlin, Germany. △ Less

Submitted 28 January, 2022; originally announced January 2022.

Comments: html page with links to the JEM-EUSO Collaboration papers presented at ICRC-2021, Berlin, Germany

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

Authors: The ANTARES collaboration, A. Albert, S. Alves, M. André, M. Anghinolfi, M. Ardid, S. Ardid, J. -J. Aubert, J. Aublin, B. Baret, S. Basa, B. Belhorma, M. Bendahman, V. Bertin, S. Biagi, M. Bissinger, J. Boumaaza, M. Bouta, M. C. Bouwhuis, H. Brânzaş, R. Bruijn, J. Brunner, J. Busto, B. Caiffi, D. Calvo , et al. (1025 additional authors not shown)

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

Submitted 23 August, 2022; v1 submitted 18 January, 2022; originally announced January 2022.

Comments: 39 pages, 7 figures, 4 tables; updated source files including xml authorlist

Report number: FERMILAB-PUB-22-033-AD-PPD-SCD-TD

Journal ref: ApJ 934 164 (2022)

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

Authors: R. U. Abbasi, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, I. Buckland, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, R. Fukushima, G. Furlich, N. Globus, R. Gonzalez, W. Hanlon, M. Hayashi , et al. (140 additional authors not shown)

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

Submitted 18 November, 2021; originally announced November 2021.

Indications of a Cosmic Ray Source in the Perseus-Pisces Supercluster

Authors: Telescope Array Collaboration, R. U. Abbasi, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, I. Buckland, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, R. Fukushima, G. Furlich, N. Globus, R. Gonzalez, W. Hanlon , et al. (135 additional authors not shown)

Abstract: The Telescope Array Collaboration has observed an excess of events with $E \ge 10^{19.4} ~{\rm eV}$ in the data which is centered at (RA, dec) = ($19^\circ$, $35^\circ$). This is near the center of the Perseus-Pisces supercluster (PPSC). The PPSC is about $70 ~{\rm Mpc}$ distant and is the closest supercluster in the Northern Hemisphere (other than the Virgo supercluster of which we are a part). A… ▽ More The Telescope Array Collaboration has observed an excess of events with $E \ge 10^{19.4} ~{\rm eV}$ in the data which is centered at (RA, dec) = ($19^\circ$, $35^\circ$). This is near the center of the Perseus-Pisces supercluster (PPSC). The PPSC is about $70 ~{\rm Mpc}$ distant and is the closest supercluster in the Northern Hemisphere (other than the Virgo supercluster of which we are a part). A Li-Ma oversampling analysis with $20^\circ$-radius circles indicates an excess in the arrival direction of events with a local significance of about 4 standard deviations. The probability of having such excess close to the PPSC by chance is estimated to be 3.5 standard deviations. This result indicates that a cosmic ray source likely exists in that supercluster. △ Less

Submitted 27 October, 2021; originally announced October 2021.

Comments: 8 pages, 4 figures, 1 table

High-entropy ejecta plumes in Cassiopeia A from neutrino-driven convection

Authors: Toshiki Sato, Keiichi Maeda, Shigehiro Nagataki, Takashi Yoshida, Brian Grefenstette, Brian J. Williams, Hideyuki Umeda, Masaomi Ono, John P. Hughes

Abstract: Recent multi-dimensional simulations suggest that high-entropy buoyant plumes help massive stars to explode. Outwardly protruding iron-rich fingers in the galactic supernova remnant Cassiopeia A are uniquely suggestive of this picture. Detecting signatures of specific elements synthesized in the high-entropy nuclear burning regime (i.e., $α$-rich freeze out) would be among the strongest substantia… ▽ More Recent multi-dimensional simulations suggest that high-entropy buoyant plumes help massive stars to explode. Outwardly protruding iron-rich fingers in the galactic supernova remnant Cassiopeia A are uniquely suggestive of this picture. Detecting signatures of specific elements synthesized in the high-entropy nuclear burning regime (i.e., $α$-rich freeze out) would be among the strongest substantiating evidence. Here we report the discovery of such elements, stable Ti and Cr, at a confidence level greater than 5$σ$ in the shocked high-velocity iron-rich ejecta of Cassiopeia A. We found the observed Ti/Fe and Cr/Fe mass ratios require $α$-rich freeze out, providing the first observational demonstration for the existence of high-entropy ejecta plumes that boosted the shock wave at explosion. The metal composition of the plumes agrees well with predictions for strongly neutrino-processed proton-rich ejecta. These results support the operation of the convective supernova engine via neutrino heating in the supernova that produced Cassiopeia A. △ Less

Submitted 20 October, 2021; originally announced October 2021.

Comments: 34 pages, 13 figures

Journal ref: Nature, Volume 592, Issue 7855, p.537-540 (2021)

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

Authors: Donald C. Warren, Maria Dainotti, Maxim V. Barkov, Bjorn Ahlgren, Hirotaka Ito, Shigehiro Nagataki

Abstract: 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) e… ▽ More 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. △ Less

Submitted 15 September, 2021; originally announced September 2021.

Comments: 32 pages, 15 figures. Includes changes requested by referee

Report number: RIKEN-iTHEMS-Report-21

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

Authors: A. Dohi, N. Nishimura, M. Hashimoto, Y. Matsuo, T. Noda, S. Nagataki

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

Submitted 5 January, 2022; v1 submitted 28 May, 2021; originally announced May 2021.

Comments: 19 pages, 16 figures, 3 tables, updated to match published version (ApJ)

Journal ref: The Astrophysical Journal, 923:64, 2021

Closure relations during the plateau emission of Swift GRBs and the fundamental plane

Authors: M. G. Dainotti, A. Ł. Lenart, N. Fraija, S. Nagataki, D. C. Warren, B. De Simone, G. Srinivasaragavan, A. Mata

Abstract: The Neil Gehrels Swift observatory observe Gamma-Ray bursts (GRBs) plateaus in X-rays. We test the reliability of the closure relations through the fireball model when dealing with the GRB plateau emission. We analyze 455 X-ray lightcurves (LCs) collected by \emph{Swift} from 2005 (January) until 2019 (August) for which the redshift is both known and unknown using the phenomenological Willingale 2… ▽ More The Neil Gehrels Swift observatory observe Gamma-Ray bursts (GRBs) plateaus in X-rays. We test the reliability of the closure relations through the fireball model when dealing with the GRB plateau emission. We analyze 455 X-ray lightcurves (LCs) collected by \emph{Swift} from 2005 (January) until 2019 (August) for which the redshift is both known and unknown using the phenomenological Willingale 2007 model. Using these fits, we analyze the emission mechanisms and astrophysical environments of these GRBs through the closure relations within the time interval of the plateau emission. Finally, we test the 3D fundamental plane relation (Dainotti relation) which connects the prompt peak luminosity, the time at the end of the plateau (rest-frame), and the luminosity at that time, on the GRBs with redshift, concerning groups determined by the closure relations. This allows us to check if the intrinsic scatter σ_{int} of any of these groups is reduced compared to previous literature. The most fulfilled environments for the electron spectral distribution, p>2, are Wind Slow Cooling (SC) and ISM Slow Cooling for cases in which the parameter q, which indicates the flatness of the plateau emission and accounts for the energy injection, is =0 and =0.5, respectively, both in the cases with known and unknown redshifts. We also find that for the sGRBs All ISM Environments with $q=0$ have the smallest σ_{int}=0.04 \pm 0.15 in terms of the fundamental plane relation holding a probability of occurring by chance of p=0.005. We have shown that the majority of GRBs presenting the plateau emission fulfil the closure relations, including the energy injection, with a particular preference for the Wind SC environment. The subsample of GRBs that fulfil given relations can be used as possible standard candles and can suggest a way to reduce the intrinsic scatter of these studied relationships. △ Less

Submitted 22 May, 2021; originally announced May 2021.

Comments: 44 pages, 23 figures; Accepted to the PASJ to be published soon

A global numerical model of the prompt emission in short gamma-ray bursts

Authors: Hirotaka Ito, Oliver Just, Yuki Takei, Shigehiro Nagataki

Abstract: We present the first global model of prompt emission from a short gamma-ray burst that consistently describes the evolution of the central black-hole (BH) torus system, the propagation of the jet through multi-component merger ejecta, the transition into free expansion, and the photospheric emission from the relativistic jet. To this end, we perform a special relativistic neutrino-hydrodynamics si… ▽ More We present the first global model of prompt emission from a short gamma-ray burst that consistently describes the evolution of the central black-hole (BH) torus system, the propagation of the jet through multi-component merger ejecta, the transition into free expansion, and the photospheric emission from the relativistic jet. To this end, we perform a special relativistic neutrino-hydrodynamics simulation of a viscous BH-torus system, which is formed about 500ms after the merger and is surrounded by dynamical ejecta as well as neutron star winds, along with a jet that is injected in the vicinity of the central BH. In a post-processing step, we compute the photospheric emission using a relativistic Monte-Carlo radiative transfer code. It is found that the wind from the torus leaves a strong imprint on the jet as well as on the emission causing narrow collimation and rapid time variability. The viewing angle dependence of the emission gives rise to correlations among the spectral peak energy, E_p, isotropic energy, E_iso, and peak luminosity, L_p, which may provide natural explanations for the Amati- and Yonetoku-relations. We also find that the degree of polarization is small for the emission from the jet core (<2%), while it tends to increase with viewing angle outside of the core and can become as high as ~10-40% for energies larger than the peak energy. Finally, the comparison of our model with GRB170817A strongly disfavors the photospheric emission scenario and therefore supports alternative scenarios, such as the cocoon shock breakout. △ Less

Submitted 22 June, 2021; v1 submitted 19 May, 2021; originally announced May 2021.

Comments: 13 pages, 6 figures, accepted for publication in ApJ

Report number: RIKEN-iTHEMS-Report-21