Showing 1–40 of 40 results for author: Kuroda, T

Numerical relativity simulations of black hole and relativistic jet formation

Authors: Takami Kuroda, Masaru Shibata

Abstract: We investigate impacts of stellar rotation and magnetic fields on black hole (BH) formation and its subsequent explosive activities, by conducting axisymmetric radiation-magnetohydrodynamics simulations of gravitational collapse of a 70 $M_\odot$ star with two-moment multi energy neutrino transport in numerical relativity. Due to its dense stellar structure, all models cannot avoid the eventual BH… ▽ More We investigate impacts of stellar rotation and magnetic fields on black hole (BH) formation and its subsequent explosive activities, by conducting axisymmetric radiation-magnetohydrodynamics simulations of gravitational collapse of a 70 $M_\odot$ star with two-moment multi energy neutrino transport in numerical relativity. Due to its dense stellar structure, all models cannot avoid the eventual BH formation even though a strongly magnetized model experiences the so-called magnetorotational explosion prior to the BH formation. One intriguing phenomenon observed in the strongly magnetized model is the formation of a relativistic jet in the post-BH formation. The relativistic jet is the outcome of a combination of strong magnetic fields and low-density materials above the BH. The jet further enhances the explosion energy beyond $\sim10^{52}$ erg, which is well exceeding the gravitational overburden ahead of the shock. Our self-consistent supernova models demonstrate that rotating magnetized massive stars at the high-mass end of supernova progenitors could be a potential candidate of hypernova and long gamma-ray burst progenitors. △ Less

Submitted 3 April, 2024; originally announced April 2024.

Comments: 6 pages, 4 figures, submitted to MNRAS Letter

Three-dimensional GRMHD Simulations of Rapidly Rotating Stellar Core-Collapse

Authors: Shota Shibagaki, Takami Kuroda, Kei Kotake, Tomoya Takiwaki, Tobias Fischer

Abstract: We present results from fully general relativistic (GR), three-dimensional (3D), neutrino-radiation magneto-hydrodynamic (MHD) simulations of stellar core collapse of a 20 M$_\odot$ star with spectral neutrino transport. Our focus is to study the gravitational-wave (GW) signatures from the magnetorotationally (MR)-driven models. By parametrically changing the initial angular velocity and the stren… ▽ More We present results from fully general relativistic (GR), three-dimensional (3D), neutrino-radiation magneto-hydrodynamic (MHD) simulations of stellar core collapse of a 20 M$_\odot$ star with spectral neutrino transport. Our focus is to study the gravitational-wave (GW) signatures from the magnetorotationally (MR)-driven models. By parametrically changing the initial angular velocity and the strength of the magnetic fields in the core, we compute four models. Our results show that the MHD outflows are produced only for models (two out of four), to which magnetic field strengths of 10$^{12}$ G and rotation rates of 1 or 2 rad s$^{-1}$ are initially imposed in the core. Seen from the direction perpendicular to the rotational axis, a characteristic waveform is obtained exhibiting a monotonic time increase in the wave amplitude. As previously identified, this stems from the propagating MHD outflows along the axis. We show that the GW amplitude from anisotropic neutrino emission becomes more than one order-of-magnitude bigger than that from the matter contribution, whereas seen from the rotational axis, both of the two components are in the same order-of-magnitudes. Due to the memory effect, the frequency of the neutrino GW from our full-fledged 3D-MHD models is in the range less than $\sim$10 Hz. Toward the future GW detection for a Galactic core-collapse supernova, if driven by the MR mechanism, the planned next-generation detector as DECIGO is urgently needed to catch the low-frequency signals. △ Less

Submitted 17 November, 2023; v1 submitted 10 September, 2023; originally announced September 2023.

Comments: 12 pages, 10 figures

Failed supernova simulations beyond black hole formation

Authors: Takami Kuroda, Masaru Shibata

Abstract: We present an axisymmetric failed supernova simulation beyond black hole formation, for the first time with numerical relativity and two-moment multi energy neutrino transport. To ensure stable numerical evolution, we use an excision method for neutrino radiation-hydrodynamics within the inner part of black hole domain. We demonstrate that our excision method is capable to stably evolve the radiat… ▽ More We present an axisymmetric failed supernova simulation beyond black hole formation, for the first time with numerical relativity and two-moment multi energy neutrino transport. To ensure stable numerical evolution, we use an excision method for neutrino radiation-hydrodynamics within the inner part of black hole domain. We demonstrate that our excision method is capable to stably evolve the radiation-hydrodynamics in dynamical black hole spacetime. As a remarkable signature of the final moment of PNS, we find the emergence of high energy neutrinos. Those high energy neutrinos are associated with the proto-neutron star shock surface being swallowed by the central black hole and could be a possible observable of failed supernovae. △ Less

Submitted 6 September, 2023; v1 submitted 12 July, 2023; originally announced July 2023.

Comments: 8 pages, 5 figures, accepted by MNRAS

Spontaneous Scalarization as a New Core-Collapse Supernova Mechanism and its Multi-Messenger Signals

Authors: Takami Kuroda, Masaru Shibata

Abstract: We perform multi-dimensional core-collapse supernova (CCSN) simulations in a massive scalar-tensor theory for the first time with a realistic equation of state and multi-energy neutrino radiation. Among the set of our models varying the scalar mass and the coupling strength between the scalar and gravitational fields, a particular model allows for recurrent spontaneous scalarizations (SSs) in the… ▽ More We perform multi-dimensional core-collapse supernova (CCSN) simulations in a massive scalar-tensor theory for the first time with a realistic equation of state and multi-energy neutrino radiation. Among the set of our models varying the scalar mass and the coupling strength between the scalar and gravitational fields, a particular model allows for recurrent spontaneous scalarizations (SSs) in the proto-neutron star (PNS). Each SS induces the PNS collapse and subsequent bounce, from which devastating shock waves emanate and eject the PNS envelope. The explosion energy can easily exceed $\mathcal O(10^{51})$ erg. This study reveals new aspects of SS as the explosion mechanism of CCSNe. We also discuss its characteristic multi-messenger signals: neutrinos and gravitational waves. △ Less

Submitted 15 May, 2023; v1 submitted 20 February, 2023; originally announced February 2023.

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

Photochemical and RadiatiOn Transport model for Extensive USe (PROTEUS)

Authors: Yuki Nakamura, Naoki Terada, Shungo Koyama, Tatsuya Yoshida, Hiroki Karyu, Kaori Terada, Takeshi Kuroda, Arihiro Kamada, Isao Murata, Shotaro Sakai, Yuhei Suzuki, Mirai Kobayashi, François Leblanc

Abstract: We introduce a new flexible one-dimensional photochemical model named Photochemical and RadiatiOn Transport model for Extensive USe (PROTEUS), which consists of a Python graphical user interface (GUI) program and Fortran 90 modules. PROTEUS is designed for adaptability to many planetary atmospheres, for flexibility to deal with thousands of or more chemical reactions with high efficiency, and for… ▽ More We introduce a new flexible one-dimensional photochemical model named Photochemical and RadiatiOn Transport model for Extensive USe (PROTEUS), which consists of a Python graphical user interface (GUI) program and Fortran 90 modules. PROTEUS is designed for adaptability to many planetary atmospheres, for flexibility to deal with thousands of or more chemical reactions with high efficiency, and for intuitive operation with GUI. Chemical reactions can be easily implemented into the Python GUI program in a simple string format, and users can intuitively select a planet and chemical reactions on GUI. Chemical reactions selected on GUI are automatically analyzed by string parsing functions in the Python GUI program, then applied to the Fortran 90 modules to simulate with the selected chemical reactions on a selected planet. PROTEUS can significantly save the time for those who need to develop a new photochemical model; users just need to write chemical reactions in the Python GUI program and just select them on GUI to run a new photochemical model. △ Less

Submitted 6 January, 2023; originally announced January 2023.

Climate of high obliquity exo-terrestrial planets with a three-dimensional cloud system resolving climate model

Authors: Takanori Kodama, Daisuke Takasuka, Sam Sherriff-Tadano, Takeshi Kuroda, Tomoki Miyakawa, Ayako Abe-Ouchi, Masaki Satoh

Abstract: Planetary climates are strongly affected by planetary orbital parameters such as obliquity, eccentricity, and precession. In exoplanetary systems, exo-terrestrial planets should have various obliquities. High-obliquity planets would have extreme seasonal cycles due to the seasonal change of the distribution of the insolation. Here, we introduce the Non-hydrostatic ICosahedral Atmospheric Model(NIC… ▽ More Planetary climates are strongly affected by planetary orbital parameters such as obliquity, eccentricity, and precession. In exoplanetary systems, exo-terrestrial planets should have various obliquities. High-obliquity planets would have extreme seasonal cycles due to the seasonal change of the distribution of the insolation. Here, we introduce the Non-hydrostatic ICosahedral Atmospheric Model(NICAM), a global cloud-resolving model, to investigate the climate of high-obliquity planets. This model can explicitly simulate a three-dimensional cloud distribution and vertical transports of water vapor. We simulated exo-terrestrial climates with high resolution using the supercomputer FUGAKU. We assumed aqua-planet configurations with 1 bar of air as a background atmosphere, with four different obliquities ($0^{\circ}$, $23.5^{\circ}$, $45^{\circ}$, and $60^{\circ}$). We ran two sets of simulations: 1) low-resolution (~ 220 km-mesh as the standard resolution of a general circulation model for exoplanetary science) with parametrization for cloud formation, and 2) high-resolution (~ 14 km-mesh) with an explicit cloud microphysics scheme. Results suggest that high-resolution simulations with an explicit treatment of cloud microphysics reveal warmer climates due to less low cloud fraction and a large amount of water vapor in the atmosphere. It implies that treatments of cloud-related processes lead to a difference between different resolutions in climatic regimes in cases with high obliquities. △ Less

Submitted 10 October, 2022; originally announced October 2022.

Comments: 18 pages, 9 figure, 2 tables. Accepted for publication in The Astrophysical Journal

3D template-based $Fermi$-LAT constraints on the diffuse supernova axion-like particle background

Authors: Francesca Calore, Pierluca Carenza, Christopher Eckner, Tobias Fischer, Maurizio Giannotti, Joerg Jaeckel, Kei Kotake, Takami Kuroda, Alessandro Mirizzi, Francesco Sivo

Abstract: Axion-like particles (ALPs) may be abundantly produced in core-collapse (CC) supernovae (SNe), hence the cumulative signal from all past SN events can create a diffuse flux peaked at energies of about 25~MeV. We improve upon the modeling of the ALPs flux by including a set of CC SN models with different progenitor masses, as well as the effects of failed CC SNe -- which yield the formation of blac… ▽ More Axion-like particles (ALPs) may be abundantly produced in core-collapse (CC) supernovae (SNe), hence the cumulative signal from all past SN events can create a diffuse flux peaked at energies of about 25~MeV. We improve upon the modeling of the ALPs flux by including a set of CC SN models with different progenitor masses, as well as the effects of failed CC SNe -- which yield the formation of black holes instead of explosions. Relying on the coupling strength of ALPs to photons and the related Primakoff process, the diffuse SN ALP flux is converted into gamma rays while traversing the magnetic field of the Milky Way. The spatial morphology of this signal is expected to follow the shape of the Galactic magnetic field lines. We make use of this via a template-based analysis that utilizes 12 years of $Fermi$-LAT data in the energy range from 50 MeV to 500 GeV. In our benchmark case of the realization of astrophysical and cosmological parameters, we find an upper limit of $g_{aγ} \lesssim 3.76\times10^{-11}\;\mathrm{GeV}^{-1}$ at 95$\%$ confidence level for $m_a \ll 10^{-11}$ eV, while we find that systematic deviations from this benchmark scenario induce an uncertainty as large as about a factor of two. Our result slightly improves the CAST bound, while still being a factor of six (baseline scenario) weaker than the SN1987A gamma-ray burst limit. △ Less

Submitted 7 March, 2022; v1 submitted 7 October, 2021; originally announced October 2021.

Comments: 26 pages, 17 figures; v2 changes: inclusion of gravitational redshift on CC SN ALP emission spectra, impact of the creation of alpha particles during a CC SN on ALP emission spectra

Report number: LAPTH-038/21

Core-collapse supernova simulations and the formation of neutron stars, hybrid stars, and black holes

Authors: Takami Kuroda, Tobias Fischer, Tomoya Takiwaki, Kei Kotake

Abstract: We investigate observable signatures of a first-order quantum chromodynamics (QCD) phase transition in the context of core collapse supernovae. To this end, we conduct axially symmetric numerical relativity simulations with multi-energy neutrino transport, using a hadron-quark hybrid equation of state (EOS). We consider four non-rotating progenitor models, whose masses range from $9.6$ to $70$ M… ▽ More We investigate observable signatures of a first-order quantum chromodynamics (QCD) phase transition in the context of core collapse supernovae. To this end, we conduct axially symmetric numerical relativity simulations with multi-energy neutrino transport, using a hadron-quark hybrid equation of state (EOS). We consider four non-rotating progenitor models, whose masses range from $9.6$ to $70$ M$_\odot$. We find that the two less massive progenitor stars (9.6 and 11.2 M$_\odot$) show a successful explosion, which is driven by the neutrino heating. They do not undergo the QCD phase transition and leave behind a neutron star (NS). As for the more massive progenitor stars (50 and 70 M$_\odot$), the proto-neutron star (PNS) core enters the phase transition region and experiences the second collapse. Because of a sudden stiffening of the EOS entering to the pure quark matter regime, a strong shock wave is formed and blows off the PNS envelope in the 50 M$_\odot$ model. Consequently, the remnant becomes a quark core surrounded by hadronic matters, leading to the formation of the hybrid star. However for the 70 M$_\odot$ model, the shock wave cannot overcome the continuous mass accretion and it readily becomes a black hole. We find that the neutrino and gravitational wave (GW) signals from supernova explosions driven by the hadron-quark phase transition are detectable for the present generation of neutrino and GW detectors. Furthermore, the analysis of the GW detector response reveals unique kHz signatures, which will allow us to distinguish this class of supernova explosions from failed and neutrino-driven explosions. △ Less

Submitted 20 October, 2021; v1 submitted 3 September, 2021; originally announced September 2021.

Comments: 24 pages, 16 figures, Accepted for publication in ApJ

Can we constrain the origin of Mars' recurring slope lineae using atmospheric observations?

Authors: H. Kurokawa, T. Kuroda, S. Aoki, H. Nakagawa

Abstract: Flowing water and brine have been proposed to cause seasonally reappearing dark streaks called recurring slope lineae (RSL) on steep warm slopes on Mars, along with other formation mechanisms that do not involve water. This study aims to examine whether the evaporation of water vapor from the RSL, whether from fresh water or brine, is detectable by observing water vapor and/or clouds. In this stud… ▽ More Flowing water and brine have been proposed to cause seasonally reappearing dark streaks called recurring slope lineae (RSL) on steep warm slopes on Mars, along with other formation mechanisms that do not involve water. This study aims to examine whether the evaporation of water vapor from the RSL, whether from fresh water or brine, is detectable by observing water vapor and/or clouds. In this study, we summarize the possible rate and duration of water-vapor emission from RSL in different scenarios, simulate how the emitted water vapor behaves in a global climate model, and discuss the detectability of water vapor in nadir observations during existing and future explorations. We found that, in typical cases, rapid horizontal dissipation within the planetary boundary layer (PBL) following the release of water vapor prohibits cloud formation and the excess water vapor from being distinguished from the background with existing observations. Thus, we conclude that the lack of correlation between the RSL activities and the overlying water-vapor column density does not necessarily rule out the wet origin of RSL. Nevertheless, we also found that water vapor tends to accumulate in basins and valleys in some cases due to the combined effects of topography and low PBL; we suggest the locations of such configuration as targets for future atmospheric studies of Mars dedicated to quantifying water-vapor release (associated with RSL) to elucidate the formation mechanism(s) of the RSL on the planet. △ Less

Submitted 27 August, 2021; originally announced August 2021.

Comments: 52 pages, 15 figures, accepted for publication in Icarus. Supplementary movies will become available from the journal

Application of the Hilbert-Huang transform for analyzing standing-accretion-shock-instability induced gravitational waves in a core-collapse supernova

Authors: M. Takeda, Y. Hiranuma, N. Kanda, K. Kotake, T. Kuroda, R. Negishi, K. Oohara, K. Sakai, Y. Sakai, T. Sawada, H. Takahashi, S. Tsuchida, Y. Watanabe, T. Yokozawa

Abstract: Through numerical simulations, it is predicted that the gravitational waves (GWs) reflect the characteristics of the core-collapse supernova (CCSN) explosion mechanism. There are multiple GW excitation processes that occur inside a star before its explosion, and it is suggested that the GWs originating from the CCSN have a mode for each excitation process in terms of time-frequency representation.… ▽ More Through numerical simulations, it is predicted that the gravitational waves (GWs) reflect the characteristics of the core-collapse supernova (CCSN) explosion mechanism. There are multiple GW excitation processes that occur inside a star before its explosion, and it is suggested that the GWs originating from the CCSN have a mode for each excitation process in terms of time-frequency representation. Therefore, we propose an application of the Hilbert-Huang Transform (HHT), which is a high-resolution time-frequency analysis method, to analyze these GW modes for theoretically probing and increasing our understanding of the explosion mechanism. The HHT defines frequency as a function of time, and is not bound by the trade-off between time and frequency resolutions. In this study, we analyze a gravitational waveform obtained from a three-dimensional general-relativistic CCSN model that showed a vigorous activity of the standing-accretion-shock-instability (SASI). We succeed in extracting the SASI induced GWs with high resolution on a time-frequency representation using the HHT and we examine their instantaneous frequencies. △ Less

Submitted 19 October, 2021; v1 submitted 12 July, 2021; originally announced July 2021.

Comments: 12 pages, 9 figures

Report number: JGW-P2112746

Journal ref: Phys. Rev. D 104, 084063 (2021)

Post-explosion evolution of core-collapse supernovae

Authors: M. Witt, A. Psaltis, H. Yasin, C. Horn, M. Reichert, T. Kuroda, M. Obergaulinger, S. M. Couch, A. Arcones

Abstract: We investigate the post-explosion phase in core-collapse supernovae with 2D hydrodynamical simulations and a simple neutrino treatment. The latter allows us to perform 46 simulations and follow the evolution of the 32 successful explosions during several seconds. We present a broad study based on three progenitors (11.2 $M_\odot$, 15 $M_\odot$, and 27 $M_\odot$), different neutrino-heating efficie… ▽ More We investigate the post-explosion phase in core-collapse supernovae with 2D hydrodynamical simulations and a simple neutrino treatment. The latter allows us to perform 46 simulations and follow the evolution of the 32 successful explosions during several seconds. We present a broad study based on three progenitors (11.2 $M_\odot$, 15 $M_\odot$, and 27 $M_\odot$), different neutrino-heating efficiencies, and various rotation rates. We show that the first seconds after shock revival determine the final explosion energy, remnant mass, and properties of ejected matter. Our results suggest that a continued mass accretion increases the explosion energy even at late times. We link the late-time mass accretion to initial conditions such as rotation strength and shock deformation at explosion time. Only some of our simulations develop a neutrino-driven wind that survives for several seconds. This indicates that neutrino-driven winds are not a standard feature expected after every successful explosion. Even if our neutrino treatment is simple, we estimate the nucleosynthesis of the exploding models for the 15 $M_\odot$ progenitor after correcting the neutrino energies and luminosities to get a more realistic electron fraction. △ Less

Submitted 1 July, 2021; originally announced July 2021.

Comments: 16 pages, 15 figures, submitted to ApJ

Axion-like Particles from Hypernovae

Authors: Andrea Caputo, Pierluca Carenza, Giuseppe Lucente, Edoardo Vitagliano, Maurizio Giannotti, Kei Kotake, Takami Kuroda, Alessandro Mirizzi

Abstract: It was recently pointed out that very energetic subclasses of supernovae (SNe), like hypernovae and superluminous SNe, might host ultra-strong magnetic fields in their core. Such fields may catalyze the production of feebly interacting particles, changing the predicted emission rates. Here we consider the case of axion-like particles (ALPs) and show that the predicted large scale magnetic fields i… ▽ More It was recently pointed out that very energetic subclasses of supernovae (SNe), like hypernovae and superluminous SNe, might host ultra-strong magnetic fields in their core. Such fields may catalyze the production of feebly interacting particles, changing the predicted emission rates. Here we consider the case of axion-like particles (ALPs) and show that the predicted large scale magnetic fields in the core contribute significantly to the ALP production, via a coherent conversion of thermal photons. Using recent state-of-the-art SN simulations including magnetohydrodynamics, we find that if ALPs have masses $m_a \sim {\mathcal O}(10)\, \rm MeV$, their emissivity via magnetic conversions is over two orders of magnitude larger than previously estimated. Moreover, the radiative decay of these massive ALPs would lead to a peculiar delay in the arrival times of the daughter photons. Therefore, high-statistics gamma-ray satellites can potentially discover MeV ALPs in an unprobed region of the parameter space and shed light on the magnetohydrodinamical nature of the SN explosion. △ Less

Submitted 12 April, 2021; originally announced April 2021.

Comments: 6 pages, 3 Figures

Characteristic Time Variability of Gravitational-Wave and Neutrino Signals from Three-dimensional Simulations of Non-Rotating and Rapidly Rotating Stellar Core-Collapse

Authors: Shota Shibagaki, Takami Kuroda, Kei Kotake, Tomoya Takiwaki

Abstract: We present results from full general relativistic three-dimensional hydrodynamics simulations of stellar core collapse of a 70 M$_\odot$ star with spectral neutrino transport. To investigate the impact of rotation on non-axisymmetric instabilities, we compute three models by parametrically changing the initial strength of rotation. The most rapidly rotating model exhibits a transient development o… ▽ More We present results from full general relativistic three-dimensional hydrodynamics simulations of stellar core collapse of a 70 M$_\odot$ star with spectral neutrino transport. To investigate the impact of rotation on non-axisymmetric instabilities, we compute three models by parametrically changing the initial strength of rotation. The most rapidly rotating model exhibits a transient development of the low-$T/|W|$ instability with one-armed spiral flow at the early postbounce phase. Subsequently, the two-armed spiral flow appears, which persists during the simulation time. The moderately rotating model also shows the growth of the low-$T/|W|$ instability, but only with the two-armed spiral flow. In the nonrotating model, a vigorous activity of the standing accretion-shock instability (SASI) is only observed. The SASI is first dominated by the sloshing mode, which is followed by the spiral SASI until the black hole formation. We present a spectrogram analysis of the gravitational waves (GWs) and neutrinos, focusing on the time correlation. Our results show that characteristic time modulations in the GW and neutrino signals can be linked to the growth of the non-axisymmetric instabilities. We find that the degree of the protoneutron star (PNS) deformation, depending upon which modes of the non-axisymmetric instabilities develop, predominantly affects the characteristic frequencies of the correlated GW and neutrino signals. We point out that these signals would be simultaneously detectable by the current-generation detectors up to $\sim10$ kpc. Our findings suggest that the joint observation of GWs and neutrinos is indispensable for extracting information on the PNS evolution preceding the black hole formation. △ Less

Submitted 13 January, 2021; v1 submitted 8 October, 2020; originally announced October 2020.

Comments: 20 pages, 19 figures

Impact of magnetic field on neutrino-matter interactions in core-collapse supernova

Authors: Takami Kuroda

Abstract: We explore the impact of magnetic field on neutrino-matter interactions in core-collapse supernova. We first derive the modified source terms for neutrino-nucleon scattering and neutrino absorption and emission processes in the moment formalism. Then we perform full relativistic three-dimensional, magnetorotational core-collapse supernova simulations of a 20 $M_\odot$ star with spectral neutrino t… ▽ More We explore the impact of magnetic field on neutrino-matter interactions in core-collapse supernova. We first derive the modified source terms for neutrino-nucleon scattering and neutrino absorption and emission processes in the moment formalism. Then we perform full relativistic three-dimensional, magnetorotational core-collapse supernova simulations of a 20 $M_\odot$ star with spectral neutrino transport. Our simulations treat self-consistently the parity violation effects of weak interaction in the presence of external magnetic field. The result shows a significant global asymmetry, mostly confined in the proto-neutron star, with clearly reflecting the magnetic field structure. The asymmetric property arises from two factors: the angle between the neutrino flux and magnetic field, and the term, which is parallel to the magnetic field and is also proportional to the deviation of distribution function of neutrinos from thermal equilibrium. The typical correction value amounts to $\sim1$ % relative to the total neutrino-matter interaction rate for the magnetic field strength of $\sim 10^{15-16}$~G. Although these asymmetric properties do not immediately affect the explosion dynamics, our results imply that they would be significant once the neutrinos diffuse out the proto-neutron star core carrying those asymmetries away. We also show that, during our simulation time of $\sim370$ ms after bounce, our results indicate that the correction value due to the modified inelastic scattering process dominates over that of the modified neutrino absorption and emission process. △ Less

Submitted 27 November, 2020; v1 submitted 16 September, 2020; originally announced September 2020.

Comments: 21 pages, 8 figures; accepted for publication in ApJ

Magnetorotational Explosion of A Massive Star Supported by Neutrino Heating in General Relativistic Three Dimensional Simulations

Authors: Takami Kuroda, Almudena Arcones, Tomoya Takiwaki, Kei Kotake

Abstract: We present results of three-dimensional (3D), radiation-magnetohydrodynamics (MHD) simulations of core-collapse supernovae in full general relativity (GR) with spectral neutrino transport. In order to study the effects of progenitor's rotation and magnetic fields, we compute three models, where the precollapse rotation rate and magnetic fields are included parametrically to a 20 M$_{\odot}$ star.… ▽ More We present results of three-dimensional (3D), radiation-magnetohydrodynamics (MHD) simulations of core-collapse supernovae in full general relativity (GR) with spectral neutrino transport. In order to study the effects of progenitor's rotation and magnetic fields, we compute three models, where the precollapse rotation rate and magnetic fields are included parametrically to a 20 M$_{\odot}$ star. While we find no shock revival in our two non-magnetized models during our simulation times ($\sim500$ ms after bounce), the magnetorotationally (MR) driven shock expansion immediately initiates after bounce in our rapidly rotating and strongly magnetized model. We show that the expansion of the MR-driven flows toward the polar directions is predominantly driven by the magnetic pressure, whereas the shock expansion toward the equatorial direction is supported by neutrino heating. Our detailed analysis indicates that the growth of the so-called kink instability may hinder the collimation of jets, resulting in the formation of broader outflows. Furthermore we find a dipole emission of lepton number, only in the MR explosion model, whose asymmetry is consistent with the explosion morphology. Although it is similar to the lepton-number emission self-sustained asymmetry (LESA), our analysis shows that the dipole emission occurs not from the protoneutron star convection zone but from above the neutrino sphere indicating that it is not associated with the LESA. We also report several unique neutrino signatures, which are significantly dependent on both the time and the viewing angle, if observed, possibly providing a rich information regarding the onset of the MR-driven explosion. △ Less

Submitted 12 May, 2020; v1 submitted 4 March, 2020; originally announced March 2020.

Comments: 19 pages, 14 figures, accepted for publication in the ApJ

A New Gravitational Wave Signature of Low-$T/|W|$ Instability in Rapidly Rotating Stellar Core Collapse

Authors: Shota Shibagaki, Takami Kuroda, Kei Kotake, Tomoya Takiwaki

Abstract: We present results from a full general relativistic three-dimensional hydrodynamics simulation of rapidly rotating core-collapse of a 70 M$_{\odot}$ star with three-flavor spectral neutrino transport. We find a strong gravitational wave (GW) emission that originates from the growth of the one- and two-armed spiral waves extending from the nascent proto-neutron star (PNS). The GW spectrogram shows… ▽ More We present results from a full general relativistic three-dimensional hydrodynamics simulation of rapidly rotating core-collapse of a 70 M$_{\odot}$ star with three-flavor spectral neutrino transport. We find a strong gravitational wave (GW) emission that originates from the growth of the one- and two-armed spiral waves extending from the nascent proto-neutron star (PNS). The GW spectrogram shows several unique features that are produced by the non-axisymmetric instabilities. After bounce, the spectrogram first shows a transient quasi-periodic time modulation at $\sim$ 450 Hz. In the second active phase, it again shows the quasi-periodic modulation but with the peak frequency increasing with time, which continues until the final simulation time. From our detailed analysis, such features can be well explained by a combination of the so-called low-$T/|W|$ instability and the PNS core contraction. △ Less

Submitted 6 February, 2020; v1 submitted 20 September, 2019; originally announced September 2019.

Comments: 6 pages, 5 figures, accepted for publication in MNRAS Letters

Dependence of outer boundary condition on protoneutron star asteroseismology with gravitational-wave signatures

Authors: Hajime Sotani, Takami Kuroda, Tomoya Takiwaki, Kei Kotake

Abstract: To obtain the eigenfrequencies of a protoneutron star (PNS) in the postbounce phase of core-collapse supernovae (CCSNe), we perform a linear perturbation analysis of the angle-averaged PNS profiles using results from a general relativistic CCSN simulation of a $15 M_{\odot}$ star. In this work, we investigate how the choice of the outer boundary condition could affect the PNS oscillation modes in… ▽ More To obtain the eigenfrequencies of a protoneutron star (PNS) in the postbounce phase of core-collapse supernovae (CCSNe), we perform a linear perturbation analysis of the angle-averaged PNS profiles using results from a general relativistic CCSN simulation of a $15 M_{\odot}$ star. In this work, we investigate how the choice of the outer boundary condition could affect the PNS oscillation modes in the linear analysis. By changing the density at the outer boundary of the PNS surface in a parametric manner, we show that the eigenfrequencies strongly depend on the surface density. By comparing with the gravitational wave (GW) signatures obtained in the hydrodynamics simulation, the so-called surface $g$-mode of the PNS can be well ascribed to the fundamental oscillations of the PNS. The frequency of the fundamental oscillations can be fitted by a function of the mass and radius of the PNS similar to the case of cold neutron stars. In the case that the position of the outer boundary is chosen to cover not only the PNS but also the surrounding postshock region, we obtain the eigenfrequencies close to the modulation frequencies of the standing accretion-shock instability (SASI). However, we point out that these oscillation modes are unlikely to have the same physical origin of the SASI modes seen in the hydrodynamics simulation. We discuss possible limitations of applying the angle-averaged, linear perturbation analysis to extract the full ingredients of the CCSN GW signatures. △ Less

Submitted 10 June, 2019; originally announced June 2019.

Comments: accepted for publication in PRD

Journal ref: Phys. Rev. D 99, 123024 (2019)

A Linear and Quadratic Time-Frequency Analysis of Gravitational Waves from Core-Collapse Supernovae

Authors: Hajime Kawahara, Takami Kuroda, Tomoya Takiwaki, Kazuhiro Hayama, Kei Kotake

Abstract: Recent core-collapse supernova (CCSN) simulations have predicted several distinct features in gravitational-wave (GW) spectrograms, including a ramp-up signature due to the g-mode oscillation of the proto-neutron star (PNS) and an excess in the low-frequency domain (100-300 Hz) potentially induced by the standing accretion shock instability (SASI). These predictions motivated us to perform a sophi… ▽ More Recent core-collapse supernova (CCSN) simulations have predicted several distinct features in gravitational-wave (GW) spectrograms, including a ramp-up signature due to the g-mode oscillation of the proto-neutron star (PNS) and an excess in the low-frequency domain (100-300 Hz) potentially induced by the standing accretion shock instability (SASI). These predictions motivated us to perform a sophisticated time-frequency analysis (TFA) of the GW signals, aimed at preparation for future observations. By reanalyzing a gravitational waveform obtained in a three-dimensional general-relativistic CCSN simulation, we show that both the spectrogram with an adequate window and the quadratic TFA separate the multimodal GW signatures much more clearly compared with the previous analysis. We find that the observed low-frequency excess during the SASI active phase is divided into two components, a stronger one at 130 Hz and an overtone at 260 Hz, both of which evolve quasi-statically during the simulation time. We also identify a new mode whose frequency varies from 700 to 600 Hz. Furthermore, we develop the quadratic TFA for the Stokes I, Q, U, and V parameters as a new tool to investigate the GW circular polarization. We demonstrate that the polarization states that randomly change with time after bounce are associated with the PNS g-mode oscillation, whereas a slowly changing polarization state in the low-frequency domain is connected to the PNS core oscillation. This study demonstrates the capability of the sophisticated TFA for diagnosing the polarized CCSN GWs in order to explore their complex nature. △ Less

Submitted 30 September, 2018; originally announced October 2018.

Comments: 14 pages,18 figures. Accepted for publication in ApJ. The code for the time-frequency analysis used in the paper is available here: https://github.com/HajimeKawahara/juwvid

Core-collapse supernovae in the hall of mirrors. A three-dimensional code-comparison project

Authors: Rubén M. Cabezón, Kuo-Chuan Pan, Matthias Liebendörfer, Takami Kuroda, Kevin Ebinger, Oliver Heinimann, Friedrich-Karl Thielemann, Albino Perego

Abstract: Modeling core-collapse supernovae (CCSNe) with neutrino transport in three dimensions (3D) requires tremendous computing resources and some level of approximation. We present a first comparison study of CCSNe in 3D with different physics approximations and hydrodynamics codes. We aim to assess the impact of the hydrodynamics code, approximations for the neutrino and gravity treatments, and rotatio… ▽ More Modeling core-collapse supernovae (CCSNe) with neutrino transport in three dimensions (3D) requires tremendous computing resources and some level of approximation. We present a first comparison study of CCSNe in 3D with different physics approximations and hydrodynamics codes. We aim to assess the impact of the hydrodynamics code, approximations for the neutrino and gravity treatments, and rotation on the simulation of CCSNe in 3D. We use four different hydrodynamics codes in this work (ELEPHANT, FLASH, fGR1, and SPHYNX) in combination with two different neutrino treatments, the isotropic diffusion source approximation (IDSA) and two-moment M1, and three different gravity treatments: Newtonian, 1D General Relativity (GR) correction, and full GR). Additional parameters discussed in this study are the inclusion of neutrino-electron scattering via a parametrized deleptonization (PD) and the influence of rotation. The four codes compared in this work include Eulerian and fully Lagrangian (smoothed particle hydrodynamics) codes for the first time. They show agreement in the overall evolution of the collapse phase and early post-bounce within the range of 10% (20% in some cases). The comparison of the different neutrino treatments highlights the need to further investigate the antineutrino luminosities in IDSA, which tend to be relatively high. We also demonstrate the requirement for a more detailed heavy-lepton neutrino leakage. When comparing with a full GR code, including an M1 transport method, we confirm the influence of neutrino-electron scattering during the collapse phase, which is adequately captured by the PD scheme. Also, the effective GR potential reproduces the overall dynamic evolution correctly in all Newtonian codes. Additionally, we verify that rotation aids the shock expansion and estimate the overall angular momentum losses for each code in rotating scenarios. △ Less

Submitted 18 December, 2018; v1 submitted 24 June, 2018; originally announced June 2018.

Comments: 20 pages, 11 figures, 4 Tables. Published in A&A

Journal ref: A&A 619, A118 (2018)

Circular polarization of gravitational waves from non-rotating supernova cores: a new probe into the pre-explosion hydrodynamics

Authors: Kazuhiro Hayama, Takami Kuroda, Kei Kotake, Tomoya Takiwaki

Abstract: We present an analysis of the circular polarization of gravitational-waves (GWs) using results from three-dimensional (3D), general relativistic (GR) core-collapse simulations of a non-rotating $15 M_{\odot}$ star. For the signal detection, we perform a coherent network analysis taking into account the four interferometers of LIGO Hanford, LIGO Livingston, VIRGO, and KAGRA. We focus on the Stokes… ▽ More We present an analysis of the circular polarization of gravitational-waves (GWs) using results from three-dimensional (3D), general relativistic (GR) core-collapse simulations of a non-rotating $15 M_{\odot}$ star. For the signal detection, we perform a coherent network analysis taking into account the four interferometers of LIGO Hanford, LIGO Livingston, VIRGO, and KAGRA. We focus on the Stokes $V$ parameter, which directly characterizes the asymmetry of the GW circular polarization. We find that the amplitude of the GW polarization becomes bigger for our 3D-GR model that exhibits strong activity of the standing accretion shock instability (SASI). Our results suggest that the SASI-induced accretion flows to the proto-neutron star (PNS) lead to a characteristic, low-frequency modulation (100 $\sim$ 200 Hz) in both the waveform and the GW circular polarization. During the vigorous SASI phase, we observe that the GW polarization switches from the right- to left-handed mode, which is clearly visible in the spectrogram. By estimating the signal-to-noise ratio of the GW polarization, we demonstrate that the detection horizon of the circular polarization extends by more than a factor of several times farther comparing to that of the GW amplitude. Our results suggest that the GW circular polarization, if detected, could provide a new probe into the pre-explosion hydrodynamics such as the SASI activity and the $g$-mode oscillation of the PNS. △ Less

Submitted 11 February, 2018; originally announced February 2018.

Comments: 5 pages, 4 figures

A full general relativistic neutrino radiation-hydrodynamics simulation of a collapsing very massive star and the formation of a black hole

Authors: Takami Kuroda, Kei Kotake, Tomoya Takiwaki, Friedrich-Karl Thielemann

Abstract: We study the final fate of a very massive star by performing full general relativistic (GR), three-dimensional (3D) simulation with three-flavor multi-energy neutrino transport. Utilizing a 70 solar mass zero metallicity progenitor, we self-consistently follow the radiation-hydrodynamics from the onset of gravitational core-collapse until the second collapse of the proto-neutron star (PNS), leadin… ▽ More We study the final fate of a very massive star by performing full general relativistic (GR), three-dimensional (3D) simulation with three-flavor multi-energy neutrino transport. Utilizing a 70 solar mass zero metallicity progenitor, we self-consistently follow the radiation-hydrodynamics from the onset of gravitational core-collapse until the second collapse of the proto-neutron star (PNS), leading to black hole (BH) formation. Our results show that the BH formation occurs at a post-bounce time of ~300 ms for the 70 Msun star. This is significantly earlier than those in the literature where lower mass progenitors were employed. At a few ~10 ms before BH formation, we find that the stalled bounce shock is revived by intense neutrino heating from the very hot PNS, which is aided by violent convection behind the shock. In the context of 3D-GR core-collapse modeling with multi-energy neutrino transport, our numerical results present the first evidence to validate a fallback BH formation scenario of the 70 Msun star. △ Less

Submitted 4 April, 2018; v1 submitted 4 January, 2018; originally announced January 2018.

Comments: 5 pages, 4 figures, accepted for publication in MNRAS letter

Nucleosynthesis in 2D Core-Collapse Supernovae of 11.2 and 17.0 M$_{\odot}$ Progenitors: Implications for Mo and Ru Production

Authors: M. Eichler, K. Nakamura, T. Takiwaki, T. Kuroda, K. Kotake, M. Hempel, R. Cabezón, M. Liebendörfer, F-K. Thielemann

Abstract: Core-collapse supernovae are the first polluters of heavy elements in the galactic history. As such, it is important to study the nuclear compositions of their ejecta, and understand their dependence on the progenitor structure (e.g., mass, compactness, metallicity). Here, we present a detailed nucleosynthesis study based on two long-term, two-dimensional core-collapse supernova simulations of a 1… ▽ More Core-collapse supernovae are the first polluters of heavy elements in the galactic history. As such, it is important to study the nuclear compositions of their ejecta, and understand their dependence on the progenitor structure (e.g., mass, compactness, metallicity). Here, we present a detailed nucleosynthesis study based on two long-term, two-dimensional core-collapse supernova simulations of a 11.2 M$_{\odot}$ and a 17.0 M$_{\odot}$ star. We find that in both models nuclei well beyond the iron group (up to $Z \approx 44$) can be produced, and discuss in detail also the nucleosynthesis of the p-nuclei $^{92,94}$Mo and $^{96,98}$Ru. While we observe the production of $^{92}$Mo and $^{94}$Mo in slightly neutron-rich conditions in both simulations, $^{96,98}$Ru can only be produced efficiently via the $ν$p-process. Furthermore, the production of Ru in the $ν$p-process heavily depends on the presence of very proton-rich material in the ejecta. This disentanglement of production mechanisms has interesting consequences when comparing to the abundance ratios between these isotopes in the solar system and in presolar grains. △ Less

Submitted 28 August, 2017; originally announced August 2017.

Comments: 48 pages, 19 figures, accepted for publication in: J. Phys. G: Nucl. Part. Phys

Correlated Signatures of Gravitational-Wave and Neutrino Emission in Three-Dimensional General-Relativistic Core-Collapse Supernova Simulations

Authors: Takami Kuroda, Kei Kotake, Kazuhiro Hayama, Tomoya Takiwaki

Abstract: We present results from general-relativistic (GR) three-dimensional (3D) core-collapse simulations with approximate neutrino transport for three non-rotating progenitors (11.2, 15, and 40 Msun) using different nuclear equations of state (EOSs). We find that the combination of progenitor's higher compactness at bounce and the use of softer EOS leads to stronger activity of the standing accretion sh… ▽ More We present results from general-relativistic (GR) three-dimensional (3D) core-collapse simulations with approximate neutrino transport for three non-rotating progenitors (11.2, 15, and 40 Msun) using different nuclear equations of state (EOSs). We find that the combination of progenitor's higher compactness at bounce and the use of softer EOS leads to stronger activity of the standing accretion shock instability (SASI). We confirm previous predications that the SASI produces characteristic time modulations both in neutrino and gravitational-wave (GW) signals. By performing a correlation analysis of the SASI-modulated neutrino and GW signals, we find that the correlation becomes highest when we take into account the time-delay effect due to the advection of material from the neutrino sphere to the proto-neutron star core surface. Our results suggest that the correlation of the neutrino and GW signals, if detected, would provide a new signature of the vigorous SASI activity in the supernova core, which can be hardly seen if neutrino-convection dominates over the SASI. △ Less

Submitted 4 November, 2017; v1 submitted 17 August, 2017; originally announced August 2017.

Comments: 24 pages, 10 figures, Accepted for publication in ApJ

Probing mass-radius relation of protoneutron stars from gravitational-wave asteroseismology

Authors: Hajime Sotani, Takami Kuroda, Tomoya Takiwaki, Kei Kotake

Abstract: The gravitational-wave (GW) asteroseismology is a powerful technique for extracting interior information of compact objects. In this work, we focus on spacetime modes, the so-called $w$-modes, of GWs emitted from a proto-neutron star (PNS) in the postbounce phase of core-collapse supernovae. Using results from recent three-dimensional supernova models, we study how to infer the properties of the P… ▽ More The gravitational-wave (GW) asteroseismology is a powerful technique for extracting interior information of compact objects. In this work, we focus on spacetime modes, the so-called $w$-modes, of GWs emitted from a proto-neutron star (PNS) in the postbounce phase of core-collapse supernovae. Using results from recent three-dimensional supernova models, we study how to infer the properties of the PNS based on a quasi-normal mode analysis in the context of the GW asteroseismology. We find that the $w_1$-mode frequency multiplied by the PNS radius is expressed as a linear function with respect to the ratio of the PNS mass to the PNS radius. This relation is insensitive to the nuclear equation of state (EOS) employed in this work. Combining with another universal relation of the $f$-mode oscillations, we point out that the time dependent mass-radius relation of the PNS can be obtained by observing both the $f$- and $w_1$-mode GWs simultaneously. Our results suggest that the simultaneous detection of the two modes could provide a new probe into finite-temperature nuclear EOS that predominantly determines the PNS evolution. △ Less

Submitted 11 August, 2017; originally announced August 2017.

Comments: accepted for publication in PRD

Journal ref: Phys. Rev. D 96, 063005 (2017)

Circular polarizations of gravitational waves from core-collapse supernovae: a clear indication of rapid rotation

Authors: Kazuhiro Hayama, Takami Kuroda, Ko Nakamura, Shoichi Yamada

Abstract: We propose to employ the circular polarization of gravitational waves emitted by core-collapse supernovae as an unequivocal indication of rapid rotation deep in their cores. It has been demonstrated by three dimensional simulations that non-axisymmetric accretion flows may develop spontaneously via hydrodynamical instabilities in the post-bounce cores. It is not surprising then that the gravitatio… ▽ More We propose to employ the circular polarization of gravitational waves emitted by core-collapse supernovae as an unequivocal indication of rapid rotation deep in their cores. It has been demonstrated by three dimensional simulations that non-axisymmetric accretion flows may develop spontaneously via hydrodynamical instabilities in the post-bounce cores. It is not surprising then that the gravitational waves emitted by such fluid motions are circularly polarized. We show in this letter that a network of the second generation detectors of gravitational waves worldwide may be able to detect such polarizations up to the opposite side of Galaxy as long as the rotation period is shorter than a few seconds prior to collapse. △ Less

Submitted 7 June, 2016; v1 submitted 5 June, 2016; originally announced June 2016.

Comments: 5 pages, 3 figures

Journal ref: Physical Review Letters, Volume 116, Issue 15, id.151102, 2016

A New Gravitational-Wave Signature from Standing Accretion Shock Instabilities in Supernovae

Authors: Takami Kuroda, Kei Kotake, Tomoya Takiwaki

Abstract: We present results from fully relativistic three-dimensional core-collapse supernova (CCSN) simulations of a non-rotating 15 Msun star using three different nuclear equations of state (EoSs). From our simulations covering up to ~350 ms after bounce, we show that the development of the standing accretion shock instability (SASI) differs significantly depending on the stiffness of nuclear EoS. Gener… ▽ More We present results from fully relativistic three-dimensional core-collapse supernova (CCSN) simulations of a non-rotating 15 Msun star using three different nuclear equations of state (EoSs). From our simulations covering up to ~350 ms after bounce, we show that the development of the standing accretion shock instability (SASI) differs significantly depending on the stiffness of nuclear EoS. Generally, the SASI activity occurs more vigorously in models with softer EoS. By evaluating the gravitational-wave (GW) emission, we find a new GW signature on top of the previously identified one, in which the typical GW frequency increases with time due to an accumulating accretion to the proto-neutron star (PNS). The newly observed quasi-periodic signal appears in the frequency range from ~100 to 200 Hz and persists for ~150 ms before neutrino-driven convection dominates over the SASI. By analyzing the cycle frequency of the SASI sloshing and spiral modes as well as the mass accretion rate to the emission region, we show that the SASI frequency is correlated with the GW frequency. This is because the SASI-induced temporary perturbed mass accretion strike the PNS surface, leading to the quasi-periodic GW emission. Our results show that the GW signal, which could be a smoking-gun signature of the SASI, is within the detection limits of LIGO, advanced Virgo, and KAGRA for Galactic events. △ Less

Submitted 24 August, 2016; v1 submitted 30 May, 2016; originally announced May 2016.

Comments: 7 pages, 5 figures, Accepted for publication in ApJL

A global view of gravity waves in the Martian atmosphere inferred from a high-resolution general circulation model

Authors: Takeshi Kuroda, Alexander S. Medvedev, Erdal Yiğit, Paul Hartogh

Abstract: Global characteristics of the small-scale gravity wave (GW) field in the Martian atmosphere obtained from a high-resolution general circulation model (GCM) are presented for the first time. The simulated GW-induced temperature variances are in a good agreement with available radio occultation data in the lower atmosphere between 10 and 30 km. The model reveals a latitudinal asymmetry with stronger… ▽ More Global characteristics of the small-scale gravity wave (GW) field in the Martian atmosphere obtained from a high-resolution general circulation model (GCM) are presented for the first time. The simulated GW-induced temperature variances are in a good agreement with available radio occultation data in the lower atmosphere between 10 and 30 km. The model reveals a latitudinal asymmetry with stronger wave generation in the winter hemisphere, and two distinctive sources of GWs: mountainous regions and the meandering winter polar jet. Orographic GWs are filtered while propagating upward, and the mesosphere is primarily dominated by harmonics with faster horizontal phase velocities. Wave fluxes are directed mainly against the local wind. GW dissipation in the upper mesosphere generates body forces of tens of m~s$^{-1}$~sol$^{-1}$, which tend to close the simulated jets. The results represent a realistic surrogate for missing observations, which can be used for constraining GW parameterizations and validating GCM simulations. △ Less

Submitted 1 November, 2015; originally announced November 2015.

Comments: Accepted for publication in Geophysical Research Letters

High-altitude gravity waves in the Martian thermosphere observed by MAVEN/NGIMS and modeled by a gravity wave scheme

Authors: Erdal Yiğit, Scott L. England, Guiping Liu, Alexander S. Medvedev, Paul R. Mahaffy, Takeshi Kuroda, Bruce M. Jakosky

Abstract: First high-altitude observations of gravity wave (GW)-induced CO$_2$ density perturbations in the Martian thermosphere retrieved from NASA's NGIMS instrument on board the MAVEN satellite are presented and interpreted using the extended GW parameterization of Yiğit et al. [2008] and the Mars Climate Database as an input. Observed relative density perturbations between 180-220 km of 20-40 % demonstr… ▽ More First high-altitude observations of gravity wave (GW)-induced CO$_2$ density perturbations in the Martian thermosphere retrieved from NASA's NGIMS instrument on board the MAVEN satellite are presented and interpreted using the extended GW parameterization of Yiğit et al. [2008] and the Mars Climate Database as an input. Observed relative density perturbations between 180-220 km of 20-40 % demonstrate appreciable local time, latitude, and altitude variations. Modeling for the spatiotemporal conditions of the MAVEN observations suggests that GWs can directly propagate from the lower atmosphere to the thermosphere, produce appreciable dynamical effects, and likely contribute to the observed fluctuations. Modeled effects are somewhat smaller than the observed but their highly variable nature is in qualitative agreement with observations. Possible reasons for discrepancies between modeling and measurements are discussed. △ Less

Submitted 12 August, 2015; originally announced August 2015.

Comments: Accepted for publication in Geophysical Research Letters (GRL). Special section: First Results from the MAVEN Mission to Mars

A New Multi-Energy Neutrino Radiation-Hydrodynamics Code in Full General Relativity and Its Application to Gravitational Collapse of Massive Stars

Authors: Takami Kuroda, Tomoya Takiwaki, Kei Kotake

Abstract: We present a new multi-dimensional radiation-hydrodynamics code for massive stellar core-collapse in full general relativity (GR). Employing an M1 analytical closure scheme, we solve spectral neutrino transport of the radiation energy and momentum based on a truncated moment formalism. Regarding neutrino opacities, we take into account a baseline set in state-of-the-art simulations, in which inela… ▽ More We present a new multi-dimensional radiation-hydrodynamics code for massive stellar core-collapse in full general relativity (GR). Employing an M1 analytical closure scheme, we solve spectral neutrino transport of the radiation energy and momentum based on a truncated moment formalism. Regarding neutrino opacities, we take into account a baseline set in state-of-the-art simulations, in which inelastic neutrinoelectron scattering, thermal neutrino production via pair annihilation and nucleonnucleon bremsstrahlung are included. While the Einstein field equations and the spatial advection terms in the radiation-hydrodynamics equations are evolved explicitly, the source terms due to neutrino-matter interactions and energy shift in the radiation moment equations are integrated implicitly by an iteration method. To verify our code, we first perform a series of standard radiation tests with analytical solutions that include the check of gravitational redshift and Doppler shift. A good agreement in these tests supports the reliability of the GR multi-energy neutrino transport scheme. We then conduct several test simulations of core-collapse, bounce, and shock-stall of a 15Msun star in the Cartesian coordinates and make a detailed comparison with published results. Our code performs quite well to reproduce the results of full-Boltzmann neutrino transport especially before bounce. In the postbounce phase, our code basically performs well, however, there are several differences that are most likely to come from the insufficient spatial resolution in our current 3D-GR models. For clarifying the resolution dependence and extending the code comparison in the late postbounce phase, we discuss that next-generation Exaflops-class supercomputers are at least needed. △ Less

Submitted 8 January, 2016; v1 submitted 26 January, 2015; originally announced January 2015.

Comments: 61 pages, 20 figures, accepted for publication in ApJS

Coherent Network Analysis of Gravitational Waves from Three-Dimensional Core-Collapse Supernova Models

Authors: Kazuhiro Hayama, Takami Kuroda, Kei Kotake, Tomoya Takiwaki

Abstract: Using predictions from three-dimensional (3D) hydrodynamics simulations of core-collapse supernovae (CCSNe), we present a coherent network analysis to detection, reconstruction, and the source localization of the gravitational-wave (GW) signals. We use the {\tt RIDGE} pipeline for the analysis, in which the network of LIGO Hanford, LIGO Livingston, VIRGO, and KAGRA is considered. By combining with… ▽ More Using predictions from three-dimensional (3D) hydrodynamics simulations of core-collapse supernovae (CCSNe), we present a coherent network analysis to detection, reconstruction, and the source localization of the gravitational-wave (GW) signals. We use the {\tt RIDGE} pipeline for the analysis, in which the network of LIGO Hanford, LIGO Livingston, VIRGO, and KAGRA is considered. By combining with a GW spectrogram analysis, we show that several important hydrodynamics features in the original waveforms persist in the waveforms of the reconstructed signals. The characteristic excess in the spectrograms originates not only from rotating core-collapse, bounce and the subsequent ring down of the proto-neutron star (PNS) as previously identified, but also from the formation of magnetohydrodynamics jets and non-axisymmetric instabilities in the vicinity of the PNS. Regarding the GW signals emitted near at the rotating core bounce, the horizon distance extends up to $\sim$ 18 kpc for the most rapidly rotating 3D model in this work. Following the rotating core bounce, the dominant source of the GW emission shifts to the non-axisymmetric instabilities. The horizon distances extend maximally up to $\sim$ 40 kpc seen from the spin axis. With an increasing number of 3D models trending towards explosion recently, our results suggest that in addition to the best studied GW signals due to rotating core-collapse and bounce, the time is ripe to consider how we can do science from GWs of CCSNe much more seriously than before. Particularly the quasi-periodic signals due to the non-axisymmetric instabilities and the detectability should deserve further investigation to elucidate the inner-working of the rapidly rotating CCSNe. △ Less

Submitted 15 November, 2015; v1 submitted 2 January, 2015; originally announced January 2015.

Comments: PRD in press

Journal ref: Phys. Rev. D 92, 122001, 2015

Systematic Features of Axisymmetric Neutrino-Driven Core-Collapse Supernova Models in Multiple Progenitors

Authors: Ko Nakamura, Tomoya Takiwaki, Takami Kuroda, Kei Kotake

Abstract: We present an overview of two-dimensional (2D) core-collapse supernova simulations employing neutrino transport scheme by the isotropic diffusion source approximation. We study 101 solar-metallicity, 247 ultra metal-poor, and 30 zero-metal progenitors covering zero-age main sequence mass from $10.8 M_{\odot} $ to $75.0 M_{\odot} $. Using the 378 progenitors in total, we systematically investigate… ▽ More We present an overview of two-dimensional (2D) core-collapse supernova simulations employing neutrino transport scheme by the isotropic diffusion source approximation. We study 101 solar-metallicity, 247 ultra metal-poor, and 30 zero-metal progenitors covering zero-age main sequence mass from $10.8 M_{\odot} $ to $75.0 M_{\odot} $. Using the 378 progenitors in total, we systematically investigate how the differences in the structures of these multiple progenitors impact the hydrodynamics evolution. By following a long-term evolution over 1.0 s after bounce, most of the computed models exhibit neutrino-driven revival of the stalled bounce shock at $\sim$ 200 - 800 ms postbounce, leading to the possibility of explosion. Pushing the boundaries of expectations in previous one-dimensional (1D) studies, our results confirm that the compactness parameter $ξ$ that characterizes the structure of the progenitors is also a key in 2D to diagnose the properties of neutrino-driven explosions. Models with high $ξ$ undergo high ram pressure from the accreting matter onto the stalled shock, which affects the subsequent evolution of the shock expansion and the mass of the protoneutron star under the influence of neutrino-driven convection and the standing accretion-shock instability. We show that the accretion luminosity becomes higher for models with high $ξ$, which makes the growth rate of the diagnostic explosion energy higher and the synthesized nickel mass bigger. We find that these explosion characteristics tend to show a monotonic increase as a function of the compactness parameter $ξ$. △ Less

Submitted 28 July, 2015; v1 submitted 10 June, 2014; originally announced June 2014.

Comments: 15 pages, 22 figures, 277 progenitors added, accepted to PASJ

Impacts of Rotation on Three-dimensional Hydrodynamics of Core-collapse Supernovae

Authors: Ko Nakamura, Takami Kuroda, Tomoya Takiwaki, Kei Kotake

Abstract: We perform a series of simplified numerical experiments to explore how rotation impacts on the three-dimensional (3D) hydrodynamics of core-collapse supernovae. For the sake of our systematic study, we employ a light-bulb scheme to trigger explosions and a three-flavor neutrino leakage scheme to treat deleptonization effects and neutrino losses from proto-neutron star interior. Using a 15 solar ma… ▽ More We perform a series of simplified numerical experiments to explore how rotation impacts on the three-dimensional (3D) hydrodynamics of core-collapse supernovae. For the sake of our systematic study, we employ a light-bulb scheme to trigger explosions and a three-flavor neutrino leakage scheme to treat deleptonization effects and neutrino losses from proto-neutron star interior. Using a 15 solar mass progenitor, we compute thirty models in 3D with a wide variety of initial angular momentum and light-bulb neutrino luminosity. We find that the rotation can help onset of neutrino-driven explosions for the models in which the initial angular momentum is matched to that obtained in recent stellar evolutionary calculations (0.3-3 rad/s at the center). For the models with larger initial angular momentum, the shock surface deforms to be more oblate due to larger centrifugal force. This makes not only a gain region more concentrated around the equatorial plane, but also the mass in the gain region bigger. As a result, buoyant bubbles tend to be coherently formed and rise in the equatorial region, which pushes the revived shock ever larger radii until a global explosion is triggered. We find that these are the main reasons that the preferred direction of explosion in 3D rotating models is often perpendicular to the spin axis, which is in sharp contrast to the polar explosions around the axis that was obtained in previous 2D simulations. △ Less

Submitted 8 July, 2014; v1 submitted 28 March, 2014; originally announced March 2014.

Comments: 15pages, 13 figures, 1 table, revised and submitted to Astrophysical Journal after referee report

Gravitational Wave Signatures from Low-mode Spiral Instabilities in Rapidly Rotating Supernova Cores

Authors: Takami Kuroda, Tomoya Takiwaki, Kei Kotake

Abstract: We study properties of gravitational waves (GWs) from rotating core-collapse of a 15M_odot star by performing three-dimensional general-relativistic hydrodynamic simulations with an approximate neutrino transport. By parametrically changing the precollapse angular momentum, we focus on the effects of rotation on the GW signatures in the early postbounce evolution. Regarding three-flavor neutrino t… ▽ More We study properties of gravitational waves (GWs) from rotating core-collapse of a 15M_odot star by performing three-dimensional general-relativistic hydrodynamic simulations with an approximate neutrino transport. By parametrically changing the precollapse angular momentum, we focus on the effects of rotation on the GW signatures in the early postbounce evolution. Regarding three-flavor neutrino transport, we solve the energy-averaged set of radiation energy and momentum based on the Thorne's momentum formalism. In addition to the gravitational quadrupole radiation from matter motions, we take into account GWs from anisotropic neutrino emission. With these computations, our results present several supporting evidences for the previous anticipation that non-axisymmetric instabilities play an essential role in determining the postbounce GW signatures. During prompt convection, we find that the waveforms show narrow-band and highly quasi-periodic signals which persist until the end of simulations. We point out that such feature reflects the growth of the one-armed spiral modes. The typical frequency of the quasi-periodic waveforms can be well explained by the propagating acoustic waves between the stalled shock and the rotating proto-neutron star surface, which suggests the appearance of the standing-accretion-shock instability. Although the GW signals exhibit strong variability between the two polarizations and different observer directions, they are within the detection limits of next generation detectors such as by KAGRA and Advanced LIGO, if the source with sufficient angular momentum is located in our Galaxy. △ Less

Submitted 14 February, 2014; v1 submitted 16 April, 2013; originally announced April 2013.

Comments: 47 pages, 20 figures, Accepted by PRD

Journal ref: Phys. Rev. D 89, 044011 (2014)

Rapid Spin Deceleration of Magnetized Proto-Neutron Stars via Asymmetric Neutrino Emission

Authors: T. Maruyama, J. Hidaka, T. Kajino, N. Yasutake, T. Kuroda, M. K. Cheoun, C. Y. Ryu, Grant. J. Mathews

Abstract: We calculate the spin deceleration of proto-neutron stars by asymmetric neutrino absorption in the context of a fully relativistic mean field theory. We calculate for the first time the spin deceleration of neutron stars due to asymmetric neutrino absorption in a toroidal magnetic field configurations. We find a surprising effect that the deceleration can much larger for asymmetric neutrino absorp… ▽ More We calculate the spin deceleration of proto-neutron stars by asymmetric neutrino absorption in the context of a fully relativistic mean field theory. We calculate for the first time the spin deceleration of neutron stars due to asymmetric neutrino absorption in a toroidal magnetic field configurations. We find a surprising effect that the deceleration can much larger for asymmetric neutrino absorption in a toroidal magnetic field than the usually presumed braking due to magnetic dipole radiation. This may explain the observation that magnetars appear to have had a more rapid deceleration in the past. △ Less

Submitted 5 March, 2014; v1 submitted 30 January, 2013; originally announced January 2013.

Comments: 13 pages, 1 figure, 1 table

Journal ref: Physical Review C89, 035801 (2014)

Core-Collapse Supernovae as Supercomputing Science: a status report toward 6D simulations with exact Boltzmann neutrino transport in full general relativity

Authors: Kei Kotake, Kohsuke Sumiyoshi, Shoichi Yamada, Tomoya Takiwaki, Takami Kuroda, Yudai Suwa, Hiroki Nagakura

Abstract: This is a status report on our endeavor to reveal the mechanism of core-collapse supernovae (CCSNe) by large-scale numerical simulations. Multi-dimensionality of the supernova engine, general relativistic magnetohydrodynamics, energy and lepton number transport by neutrinos emitted from the forming neutron star as well as nuclear interactions there, are all believed to play crucial roles in repell… ▽ More This is a status report on our endeavor to reveal the mechanism of core-collapse supernovae (CCSNe) by large-scale numerical simulations. Multi-dimensionality of the supernova engine, general relativistic magnetohydrodynamics, energy and lepton number transport by neutrinos emitted from the forming neutron star as well as nuclear interactions there, are all believed to play crucial roles in repelling infalling matter and producing energetic explosions. These ingredients are nonlinearly coupled with one another in the dynamics of core-collapse, bounce, and shock expansion. Serious quantitative studies of CCSNe hence make extensive numerical computations mandatory. Since neutrinos are neither in thermal nor in chemical equilibrium in general, their distributions in the phase space should be computed. This is a six dimensional (6D) neutrino transport problem and quite a challenge even for those with an access to the most advanced numerical resources such as the "K computer". To tackle this problem, we have embarked on multi-front efforts. In particular we report in this paper our recent progresses in the treatments of multi-dimensional (multi-D) radiation-hydrodynamics. We are currently proceeding on two different paths to the ultimate goal; in one approach we employ an approximate but highly efficient scheme for neutrino transport and treat 3D hydrodynamics and/or general relativity rigorously; some neutrino-driven explosions will be presented and comparisons will be made between 2D and 3D models quantitatively; in the second approach, on the other hand, exact but so far Newtonian Boltzmann equations are solved in two and three spatial dimensions; we will show some demonstrative test simulations. We will also address the perspectives of exa-scale computations on the next generation supercomputers. △ Less

Submitted 26 June, 2012; v1 submitted 29 May, 2012; originally announced May 2012.

Comments: 41 pages, 21 figures, accepted by Progress of Theoretical and Experimental Physics

Journal ref: Prog. Theor. Exp. Phys. (2012) 01A301

Fully General Relativistic Simulations of Core-Collapse Supernovae with An Approximate Neutrino Transport

Authors: Takami Kuroda, Kei Kotake, Tomoya Takiwaki

Abstract: We present results from the first generation of multi-dimensional hydrodynamic core-collapse simulations in full general relativity (GR) that include an approximate treatment of neutrino transport. Using a M1 closure scheme with an analytic variable Eddington factor, we solve the energy-independent set of radiation energy and momentum based on the Thorne's momentum formalism. To simplify the sourc… ▽ More We present results from the first generation of multi-dimensional hydrodynamic core-collapse simulations in full general relativity (GR) that include an approximate treatment of neutrino transport. Using a M1 closure scheme with an analytic variable Eddington factor, we solve the energy-independent set of radiation energy and momentum based on the Thorne's momentum formalism. To simplify the source terms of the transport equations, a methodology of multiflavour neutrino leakage scheme is partly employed. Our newly developed code is designed to evolve the Einstein field equation together with the GR radiation hydrodynamic equations. We follow the dynamics starting from the onset of gravitational core-collapse of a 15 $M_{\odot}$ star, through bounce, up to about 100 ms postbounce in this study to study how the spacial multi-dimensionality and GR would affect the dynamics in the early postbounce phase. Our 3D results support the anticipation in previous 1D results that the neutrino luminosity and average neutrino energy of any neutrino flavor in the postbounce phase increase when switching from SR to GR hydrodynamics. This is because the deeper gravitational well of GR produces more compact core structures, and thus hotter neutrino spheres at smaller radii. By analyzing the residency timescale to the neutrino-heating timescale in the gain region, we show that the criterion to initiate neutrino-driven explosions can be most easily satisfied in 3D models, irrespective of SR or GR hydrodynamics. Our results suggest that the combination of GR and 3D hydrodynamics provides the most favorable condition to drive a robust neutrino-driven explosion. △ Less

Submitted 29 May, 2012; v1 submitted 11 February, 2012; originally announced February 2012.

Comments: 50pages, 20 figures, Accepted by ApJ. Latest version with following the referee's suggestions and comments

Three Dimensional Magneto Hydrodynamical Simulations of Gravitational Collapse of a 15Msun Star

Authors: Takami Kuroda, Hideyuki Umeda

Abstract: We introduce our newly developed two different, three dimensional magneto hydrodynamical codes in detail. One of our codes is written in the Newtonian limit (NMHD) and the other is in the fully general relativistic code (GRMHD). Both codes employ adaptive mesh refinement and, in GRMHD, the metric is evolved with the "Baumgarte-Shapiro-Shibata-Nakamura" formalism known as the most stable method at… ▽ More We introduce our newly developed two different, three dimensional magneto hydrodynamical codes in detail. One of our codes is written in the Newtonian limit (NMHD) and the other is in the fully general relativistic code (GRMHD). Both codes employ adaptive mesh refinement and, in GRMHD, the metric is evolved with the "Baumgarte-Shapiro-Shibata-Nakamura" formalism known as the most stable method at present. We did several test problems and as for the first practical test, we calculated gravitational collapse of a $15M_\odot$ star. Main features found from our calculations are; (1) High velocity bipolar outflow is driven from the proto-neutronstar and moves through along the rotational axis in strongly magnetized models; (2) A one-armed spiral structure appears which is originated from the low-$|T/W|$ instability; (3) By comparing GRMHD and NMHD models, the maximum density increases about $\sim30%$ in GRMHD models due to the stronger gravitational effect. These features agree very well with previous studies and our codes are thus reliable to numerical simulation of gravitational collapse of massive stars. △ Less

Submitted 5 November, 2010; v1 submitted 7 August, 2010; originally announced August 2010.

Comments: Accepted by ApJS, 55 pages, 34 figures

The r-Process in Supersonic Neutrino-Driven Winds: The Roll of Wind Termination Shock

Authors: Takami Kuroda, Shinya Wanajo, Ken'ichi Nomoto

Abstract: Recent hydrodynamic studies of core-collapse supernovae imply that the neutrino-heated ejecta from a nascent neutron star develops to supersonic outflows. These supersonic winds are influenced by the reverse shock from the preceding supernova ejecta, forming the wind termination shock. We investigate the effects of the termination shock in neutrino-driven winds and its roll on the r-process. Sup… ▽ More Recent hydrodynamic studies of core-collapse supernovae imply that the neutrino-heated ejecta from a nascent neutron star develops to supersonic outflows. These supersonic winds are influenced by the reverse shock from the preceding supernova ejecta, forming the wind termination shock. We investigate the effects of the termination shock in neutrino-driven winds and its roll on the r-process. Supersonic outflows are calculated with a semi-analytic neutrino-driven wind model. Subsequent termination-shocked, subsonic outflows are obtained by applying the Rankine-Hugoniot relations. We find a couple of effects that can be relevant for the r-process. First is the sudden slowdown of the temperature decrease by the wind termination. Second is the entropy jump by termination-shock heating, up to several 100NAk. Nucleosynthesis calculations in the obtained winds are performed to examine these effects on the r-process. We find that 1) the slowdown of the temperature decrease plays a decisive roll to determine the r-process abundance curves. This is due to the strong dependences of the nucleosynthetic path on the temperature during the r-process freezeout phase. Our results suggest that only the termination-shocked winds with relatively small shock radii (~500km) are relevant for the bulk of the solar r-process abundances (A~100-180). The heaviest part in the solar r-process curve (A~180-200), however, can be reproduced both in shocked and unshocked winds. These results may help to constrain the mass range of supernova progenitors relevant for the r-process. We find, on the other hand, 2) negligible roles of the entropy jump on the r-process. This is a consequence that the sizable entropy increase takes place only at a large shock radius (~10,000km) where the r-process has already ceased. △ Less

Submitted 17 September, 2007; v1 submitted 15 July, 2007; originally announced July 2007.

Comments: 11 pages, 7 figures, submitted to ApJ, revised following referee's comments,Accepted for publication in ApJ

Journal ref: AIP Conf.Proc.1016:436-438,2008

SN 2006aj Associated with XRF 060218 At Late Phases: Nucleosynthesis-Signature of A Neutron Star-Driven Explosion

Authors: Keiichi Maeda, Koji Kawabata, Masaomi Tanaka, Ken'ichi Nomoto, Nozomu Tominaga, Takashi Hattori, Takeo Minezaki, Takami Kuroda, Tomoharu Suzuki, Jinsong Deng, Paolo A. Mazzali, Elena Pian

Abstract: Optical spectroscopy and photometry of SN 2006aj have been performed with the Subaru telescope at t > 200 days after GRB060218, the X-ray Flash with which it was associated. Strong nebular emission-lines with an expansion velocity of v ~ 7,300 km/s were detected. The peaked but relatively broad [OI]6300,6363 suggests the existence of ~ 2 Msun of materials in which ~1.3 Msun is oxygen. The core m… ▽ More Optical spectroscopy and photometry of SN 2006aj have been performed with the Subaru telescope at t > 200 days after GRB060218, the X-ray Flash with which it was associated. Strong nebular emission-lines with an expansion velocity of v ~ 7,300 km/s were detected. The peaked but relatively broad [OI]6300,6363 suggests the existence of ~ 2 Msun of materials in which ~1.3 Msun is oxygen. The core might be produced by a mildly asymmetric explosion. The spectra are unique among SNe Ic in (1) the absence of [CaII]7291,7324 emission, and (2) a strong emission feature at ~ 7400A, which requires ~ 0.05 Msun of newly-synthesized 58Ni. Such a large amount of stable neutron-rich Ni strongly indicates the formation of a neutron star. The progenitor and the explosion energy are constrained to 18 Msun < Mms < 22 Msun and E ~ (1 - 3) 10^{51} erg, respectively. △ Less

Submitted 8 March, 2007; v1 submitted 6 March, 2007; originally announced March 2007.

Comments: Accepted for Publication in the Astrophysical Journal Letters (2007, ApJ, 658, L5). 8 pages, including 1 table and 3 figures. Typos corrected

Journal ref: Astrophys.J.658:L5-L8,2007

Asteroseismological Observations of the Central Star of the Planetary Nebula NGC 1501

Authors: H. E. Bond, S. Kawaler, R. Ciardullo, R. Stover, T. Kuroda, T. Ishida, T. Ono, S. Tamura, H. Malasan, A. Yamasaki, O. Hashimoto, E. Kambe, M. Takeuti, T. Kato, M. Kato, J. -S. Chen, E. M. Leibowitz, M. M. Roth, T. Soffner, W. Mitsch

Abstract: We report on a global CCD time-series photometric campaign to decode the pulsations of the nucleus of the planetary nebula NGC1501. The star is hot and hydrogen-deficient, similar to the pre-white-dwarf PG 1159 stars. NGC1501 shows pulsational brightness variations of a few percent with periods ranging from 19 to 87 minutes. The variations are very complex, suggesting a pulsation spectrum that r… ▽ More We report on a global CCD time-series photometric campaign to decode the pulsations of the nucleus of the planetary nebula NGC1501. The star is hot and hydrogen-deficient, similar to the pre-white-dwarf PG 1159 stars. NGC1501 shows pulsational brightness variations of a few percent with periods ranging from 19 to 87 minutes. The variations are very complex, suggesting a pulsation spectrum that requires a long unbroken time series to resolve. Our CCD photometry of the star covers a two-week period in 1991 November, and used a global network of observatories. We obtained nearly continuous coverage over an interval of one week in the middle of the run. We have identified 10 pulsation periods, ranging from 5235 s down to 1154 s. We find strong evidence that the modes are indeed nonradial g-modes. The ratios of the frequencies of the largest-amplitude modes agree with those expected for modes that are trapped by a density discontinuity in the outer layers. We offer a model for the pulsation spectrum that includes a common period spacing of 22.3 s and a rotation period of 1.17 days; the period spacing allows us to assign a seismological mass of 0.55+/-0.03 Msun. △ Less

Submitted 6 September, 1996; v1 submitted 3 September, 1996; originally announced September 1996.

Comments: 12 pages, AASTEX, 7 tables, 6 EPS figures, to appear in AJ, 12/96 Corrected version repairs table formatting and adds missing Table 7