Normalized to: Kotake, K.
[1]
oai:arXiv.org:2005.05574 [pdf] - 2127475
Overview of KAGRA: Detector design and construction history
Akutsu, T.;
Ando, M.;
Arai, K.;
Arai, Y.;
Araki, S.;
Araya, A.;
Aritomi, N.;
Aso, Y.;
Bae, S. -W.;
Bae, Y. -B.;
Baiotti, L.;
Bajpai, R.;
Barton, M. A.;
Cannon, K.;
Capocasa, E.;
Chan, M. -L.;
Chen, C. -S.;
Chen, K. -H.;
Chen, Y. -R.;
Chu, H. -Y.;
Chu, Y-K.;
Eguchi, S.;
Enomoto, Y.;
Flaminio, R.;
Fujii, Y.;
Fukunaga, M.;
Fukushima, M.;
Ge, G. -G.;
Hagiwara, A.;
Haino, S.;
Hasegawa, K.;
Hayakawa, H.;
Hayama, K.;
Himemoto, Y.;
Hiranuma, Y.;
Hirata, N.;
Hirose, E.;
Hong, Z.;
Hsieh, B. -H.;
Huang, G. -Z.;
Huang, P. -W.;
Huang, Y. -J.;
Ikenoue, B.;
Imam, S.;
Inayoshi, K.;
Inoue, Y.;
Ioka, K.;
Itoh, Y.;
Izumi, K.;
Jung, K.;
Jung, P.;
Kajita, T.;
Kamiizumi, M.;
Kanda, N.;
Kang, G. -W.;
Kawaguchi, K.;
Kawai, N.;
Kawasaki, T.;
Kim, C.;
Kim, J.;
Kim, W.;
Kim, Y. -M.;
Kimura, N.;
Kita, N.;
Kitazawa, H.;
Kojima, Y.;
Kokeyama, K.;
Komori, K.;
Kong, A. K. H.;
Kotake, K.;
Kozakai, C.;
Kozu, R.;
Kumar, R.;
Kume, J.;
Kuo, C. -M.;
Kuo, H. -S.;
Kuroyanagi, S.;
Kusayanagi, K.;
Kwak, K.;
Lee, H. -K.;
Lee, H. -W.;
Lee, R. -K.;
Leonardi, M.;
Lin, C. -Y.;
Lin, F. -L.;
Lin, L. C. -C.;
Liu, G. -C.;
Luo, L. -W.;
Marchio, M.;
Michimura, Y.;
Mio, N.;
Miyakawa, O.;
Miyamoto, A.;
Miyazaki, Y.;
Miyo, K.;
Miyoki, S.;
Morisaki, S.;
Moriwaki, Y.;
Nagano, K.;
Nagano, S.;
Nakamura, K.;
Nakano, H.;
Nakano, M.;
Nakashima, R.;
Narikawa, T.;
Negishi, R.;
Ni, W. -T.;
Nishizawa, A.;
Obuchi, Y.;
Ogaki, W.;
Oh, J. J.;
Oh, S. -H.;
Ohashi, M.;
Ohishi, N.;
Ohkawa, M.;
Okutomi, K.;
Oohara, K.;
Ooi, C. -P.;
Oshino, S.;
Pan, K. -C.;
Pang, H. -F.;
Park, J.;
Arellano, F. E. Peña;
Pinto, I.;
Sago, N.;
Saito, S.;
Saito, Y.;
Sakai, K.;
Sakai, Y.;
Sakuno, Y.;
Sato, S.;
Sato, T.;
Sawada, T.;
Sekiguchi, T.;
Sekiguchi, Y.;
Shibagaki, S.;
Shimizu, R.;
Shimoda, T.;
Shimode, K.;
Shinkai, H.;
Shishido, T.;
Shoda, A.;
Somiya, K.;
Son, E. J.;
Sotani, H.;
Sugimoto, R.;
Suzuki, T.;
Suzuki, T.;
Tagoshi, H.;
Takahashi, H.;
Takahashi, R.;
Takamori, A.;
Takano, S.;
Takeda, H.;
Takeda, M.;
Tanaka, H.;
Tanaka, K.;
Tanaka, K.;
Tanaka, T.;
Tanaka, T.;
Tanioka, S.;
Martin, E. N. Tapia San;
Telada, S.;
Tomaru, T.;
Tomigami, Y.;
Tomura, T.;
Travasso, F.;
Trozzo, L.;
Tsang, T. T. L.;
Tsubono, K.;
Tsuchida, S.;
Tsuzuki, T.;
Tuyenbayev, D.;
Uchikata, N.;
Uchiyama, T.;
Ueda, A.;
Uehara, T.;
Ueno, K.;
Ueshima, G.;
Uraguchi, F.;
Ushiba, T.;
van Putten, M. H. P. M.;
Vocca, H.;
Wang, J.;
Wu, C. -M.;
Wu, H. -C.;
Wu, S. -R.;
Xu, W. -R.;
Yamada, T.;
Yamamoto, Ka.;
Yamamoto, Ko.;
Yamamoto, T.;
Yokogawa, K.;
Yokoyama, J.;
Yokozawa, T.;
Yoshioka, T.;
Yuzurihara, H.;
Zeidler, S.;
Zhao, Y.;
Zhu, Z. -H.
Submitted: 2020-05-12, last modified: 2020-07-02
KAGRA is a newly built gravitational-wave telescope, a laser interferometer
comprising arms with a length of 3\,km, located in Kamioka, Gifu, Japan. KAGRA
was constructed under the ground and it is operated using cryogenic mirrors
that help in reducing the seismic and thermal noise. Both technologies are
expected to provide directions for the future of gravitational-wave telescopes.
In 2019, KAGRA finished all installations with the designed configuration,
which we call the baseline KAGRA. In this occasion, we present an overview of
the baseline KAGRA from various viewpoints in a series of of articles. In this
article, we introduce the design configurations of KAGRA with its historical
background.
[2]
oai:arXiv.org:1911.10656 [pdf] - 2124569
Detectability of SASI activity in supernova neutrino signals
Submitted: 2019-11-24, last modified: 2020-06-02
We introduce a novel methodology for establishing the presence of Standing
Accretion Shock Instabilities (SASI) in the dynamics of a core collapse
supernova from the observed neutrino event rate at water- or ice-based neutrino
detectors. The methodology uses a likelihood ratio in the frequency domain as a
test-statistics; it is also employed to assess the potential to estimate the
frequency and the amplitude of the SASI modulations of the neutrino signal. The
parameter estimation errors are consistent with the minimum possible errors as
evaluated from the inverse of the Fisher information matrix, and close to the
theoretical minimum for the SASI amplitude. Using results from a core-collapse
simulation of a 15 solar-mass star by Kuroda $\it {et\, al.}$ (2017) as a test
bed for the method, we find that SASI can be identified with high confidence
for a distance to the supernova of up to $\sim 6$ kpc for IceCube and and up to
$\sim 3$ kpc for a 0.4 Mt mass water Cherenkov detector. This methodology will
aid the investigation of a future galactic supernova.
[3]
oai:arXiv.org:1908.03013 [pdf] - 2105542
Application of the independent component analysis to the iKAGRA data
KAGRA Collaboration;
Akutsu, T.;
Ando, M.;
Arai, K.;
Arai, Y.;
Araki, S.;
Araya, A.;
Aritomi, N.;
Asada, H.;
Aso, Y.;
Atsuta, S.;
Awai, K.;
Bae, S.;
Bae, Y.;
Baiotti, L.;
Bajpai, R.;
Barton, M. A.;
Cannon, K.;
Capocasa, E.;
Chan, M.;
Chen, C.;
Chen, K.;
Chen, Y.;
Chu, H.;
Chu, Y-K.;
Craig, K.;
Creus, W.;
Doi, K.;
Eda, K.;
Eguchi, S.;
Enomoto, Y.;
Flaminio, R.;
Fujii, Y.;
Fujimoto, M. -K.;
Fukunaga, M.;
Fukushima, M.;
Furuhata, T.;
Ge, G.;
Hagiwara, A.;
Haino, S.;
Hasegawa, K.;
Hashino, K.;
Hayakawa, H.;
Hayama, K.;
Himemoto, Y.;
Hiranuma, Y.;
Hirata, N.;
Hirobayashi, S.;
Hirose, E.;
Hong, Z.;
Hsieh, B. H.;
Huang, G-Z.;
Huang, P.;
Huang, Y.;
Ikenoue, B.;
Imam, S.;
Inayoshi, K.;
Inoue, Y.;
Ioka, K.;
Itoh, Y.;
Izumi, K.;
Jung, K.;
Jung, P.;
Kaji, T.;
Kajita, T.;
Kakizaki, M.;
Kamiizumi, M.;
Kanbara, S.;
Kanda, N.;
Kanemura, S.;
Kaneyama, M.;
Kang, G.;
Kasuya, J.;
Kataoka, Y.;
Kawaguchi, K.;
Kawai, N.;
Kawamura, S.;
Kawasaki, T.;
Kim, C.;
Kim, J. C.;
Kim, W. S.;
Kim, Y. -M.;
Kimura, N.;
Kinugawa, T.;
Kirii, S.;
Kita, N.;
Kitaoka, Y.;
Kitazawa, H.;
Kojima, Y.;
Kokeyama, K.;
Komori, K.;
Kong, A. K. H.;
Kotake, K.;
Kozakai, C.;
Kozu, R.;
Kumar, R.;
Kume, J.;
Kuo, C.;
Kuo, H-S.;
Kuroyanagi, S.;
Kusayanagi, K.;
Kwak, K.;
Lee, H. K.;
Lee, H. M.;
Lee, H. W.;
Lee, R.;
Leonardi, M.;
Lin, C.;
Lin, C-Y.;
Lin, F-L.;
Liu, G. C.;
Liu, Y.;
Luo, L.;
Majorana, E.;
Mano, S.;
Marchio, M.;
Matsui, T.;
Matsushima, F.;
Michimura, Y.;
Mio, N.;
Miyakawa, O.;
Miyamoto, A.;
Miyamoto, T.;
Miyazaki, Y.;
Miyo, K.;
Miyoki, S.;
Morii, W.;
Morisaki, S.;
Moriwaki, Y.;
Morozumi, T.;
Musha, M.;
Nagano, K.;
Nagano, S.;
Nakamura, K.;
Nakamura, T.;
Nakano, H.;
Nakano, M.;
Nakao, K.;
Nakashima, R.;
Narikawa, T.;
Naticchioni, L.;
Negishi, R.;
Quynh, L. Nguyen;
Ni, W. -T.;
Nishizawa, A.;
Obuchi, Y.;
Ochi, T.;
Ogaki, W.;
Oh, J. J.;
Oh, S. H.;
Ohashi, M.;
Ohishi, N.;
Ohkawa, M.;
Okutomi, K.;
Oohara, K.;
Ooi, C. P.;
Oshino, S.;
Pan, K.;
Pang, H.;
Park, J.;
Arellano, F. E. Pena;
Pinto, I.;
Sago, N.;
Saijo, M.;
Saito, S.;
Saito, Y.;
Sakai, K.;
Sakai, Y.;
Sakai, Y.;
Sakuno, Y.;
Sasaki, M.;
Sasaki, Y.;
Sato, S.;
Sato, T.;
Sawada, T.;
Sekiguchi, T.;
Sekiguchi, Y.;
Seto, N.;
Shibagaki, S.;
Shibata, M.;
Shimizu, R.;
Shimoda, T.;
Shimode, K.;
Shinkai, H.;
Shishido, T.;
Shoda, A.;
Somiya, K.;
Son, E. J.;
Sotani, H.;
Suemasa, A.;
Sugimoto, R.;
Suzuki, T.;
Suzuki, T.;
Tagoshi, H.;
Takahashi, H.;
Takahashi, R.;
Takamori, A.;
Takano, S.;
Takeda, H.;
Takeda, M.;
Tanaka, H.;
Tanaka, K.;
Tanaka, K.;
Tanaka, T.;
Tanaka, T.;
Tanioka, S.;
Martin, E. N. Tapia San;
Tatsumi, D.;
Telada, S.;
Tomaru, T.;
Tomigami, Y.;
Tomura, T.;
Travasso, F.;
Trozzo, L.;
Tsang, T.;
Tsubono, K.;
Tsuchida, S.;
Tsuzuki, T.;
Tuyenbayev, D.;
Uchikata, N.;
Uchiyama, T.;
Ueda, A.;
Uehara, T.;
Ueki, S.;
Ueno, K.;
Ueshima, G.;
Uraguchi, F.;
Ushiba, T.;
van Putten, M. H. P. M.;
Vocca, H.;
Wada, S.;
Wakamatsu, T.;
Wang, J.;
Wu, C.;
Wu, H.;
Wu, S.;
Xu, W-R.;
Yamada, T.;
Yamamoto, A.;
Yamamoto, K.;
Yamamoto, K.;
Yamamoto, S.;
Yamamoto, T.;
Yokogawa, K.;
Yokoyama, J.;
Yokozawa, T.;
Yoon, T. H.;
Yoshioka, T.;
Yuzurihara, H.;
Zeidler, S.;
Zhao, Y.;
Zhu, Z. -H.
Submitted: 2019-08-08, last modified: 2020-06-01
We apply the independent component analysis (ICA) to the real data from a
gravitational wave detector for the first time. Specifically we use the iKAGRA
data taken in April 2016, and calculate the correlations between the
gravitational wave strain channel and 35 physical environmental channels. Using
a couple of seismic channels which are found to be strongly correlated with the
strain, we perform ICA. Injecting a sinusoidal continuous signal in the strain
channel, we find that ICA recovers correct parameters with enhanced
signal-to-noise ratio, which demonstrates usefulness of this method. Among the
two implementations of ICA used here, we find the correlation method yields the
optimal result for the case environmental noises act on the strain channel
linearly.
[4]
oai:arXiv.org:1912.11489 [pdf] - 2101341
Time of Flight and Supernova Progenitor Effects on the Neutrino Halo
Submitted: 2019-12-24, last modified: 2020-05-25
We argue that the neutrino halo, a population of neutrinos that have
undergone direction-changing scattering in the stellar envelope of a
core-collapse supernova (CCSNe), is sensitive to neutrino emission history
through time of flight. We show that the constant time approximation, commonly
used in calculating the neutrino halo, does not capture the spatiotemporal
evolution of the halo neutrino population and that correcting for time of
flight can produce conditions which may trigger fast neutrino flavor
conversion. We also find that there exists a window of time early in all CCSNe
where the neutrino halo population is sufficiently small that it may be
negligible. This suggests that collective neutrino oscillation calculations
which neglect the Halo may be well founded at sufficiently early times.
[5]
oai:arXiv.org:1908.10594 [pdf] - 2115093
Neutrino halo effect on collective neutrino oscillation in iron
core-collapse supernova model of a 9.6 $M_{\odot}$ star
Submitted: 2019-08-28, last modified: 2020-05-17
We extend the multi-angle computational framework and investigate the time
evolution of the neutrino halo on collective neutrino oscillation in the core
collapse of an iron core progenitor. We find that in the case of the $9.6\, \rm
M_\odot$ progenitor adopted in this work, there are windows of time when the
effects of neutrino halo and collective neutrino oscillation are not
simultaneously large. Inside the shock, the impact of the inward-scattered halo
neutrino cannot in general be neglected compared to the outward-propagating
neutrino flux. However, during early epochs, collective neutrino oscillation is
effectively shut down by multi-angle matter suppression. During the
intermediate epoch, collective neutrino oscillation is not suppressed, but its
onset radius is beyond the still relatively small explosion shock front where
the halo is prominent. We also find in the case of the $9.6\, \rm M_\odot$
progenitor the halo neutrinos induce a delay in the onset of collective
neutrino oscillations. This causes novel flavor conversions which sharpen
collective neutrino oscillation spectral features. We predict that the
inclusion of neutrino halo effects makes neutrino signals that are more clearly
distinct from thermal emission that when halo neutrinos are omitted.
[6]
oai:arXiv.org:2003.02004 [pdf] - 2119901
Magnetorotational Explosion of A Massive Star Supported by Neutrino
Heating in General Relativistic Three Dimensional Simulations
Submitted: 2020-03-04, last modified: 2020-05-12
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.
[7]
oai:arXiv.org:1909.09730 [pdf] - 2049166
A New Gravitational Wave Signature of Low-$T/|W|$ Instability in Rapidly
Rotating Stellar Core Collapse
Submitted: 2019-09-20, last modified: 2020-02-06
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.
[8]
oai:arXiv.org:2001.08452 [pdf] - 2036387
Deep Core Convection and Dynamo in Newly-born Neutron Star
Submitted: 2020-01-23
We report results from a convection dynamo simulation of proto-neutron star
(PNS), with a nuclear equation of state (EOS) and the initial hydrodynamic
profile taken from a neutrino radiation-hydrodynamics simulation of a massive
stellar core-collapse. A moderately-rotating PNS with the spin period of $170$
ms in the lepton-driven convection stage is focused. We find that large-scale
flow and thermodynamic fields with north-south asymmetry develop in the
turbulent flow, as a consequence of the convection in the central part of the
PNS, which we call as a "deep core convection". Intriguingly, even with such a
moderate rotation, large-scale, $10^{15}$ G, magnetic field with dipole
symmetry is spontaneously built up in the PNS. The turbulent electro-motive
force arising from rotationally-constrained core convection is shown to play a
key role in the large-scale dynamo. The large-scale structures organized in the
PNS may impact the explosion dynamics of supernovae and subsequent evolution to
the neutron stars.
[10]
oai:arXiv.org:1912.08571 [pdf] - 2017329
The Recent Developmental Status of SNEGRAF: a Web-Based Gravitational
Wave Signal Analyzer
Submitted: 2019-12-18
Unveiling physical processes in a supernova is one of challenging topics of
modern physics and astrophysics since that event is due to particle physics on
a stellar scale and tightly related to nucleosynthesis in Universe.
Multi-messenger astronomy, a combination, such as of electromagnetic-wave,
gravitational-wave, and neutrino observations, will be a breakthrough to the
puzzle. To boost the research, we released a web-based gravitational wave
signal analyzer "SuperNova Event Gravitational-wave-display in Fukuoka
(SNEGRAF)" last year (Eguchi et al. 2019). We are now working on an integration
of the application with the RIDGE pipeline, which is for a coherent network
analysis between the LIGO, VIRGO, and KAGRA observations (Hayama et al. 2007),
and implemented in MATLAB. In the basic design phase, we decided to wrap RIDGE
with a simple Python script and make it listen for connections from SNEGRAF.
This design enables these two programs to be hosted on different servers
independently, and minimizes the cyber risks of RIDGE. In this paper, we report
the current developmental status of our system.
[11]
oai:arXiv.org:1910.00955 [pdf] - 2030547
An arm length stabilization system for KAGRA and future
gravitational-wave detectors
Akutsu, T.;
Ando, M.;
Arai, K.;
Arai, K.;
Arai, Y.;
Araki, S.;
Araya, A.;
Aritomi, N.;
Aso, Y.;
Bae, S.;
Bae, Y.;
Baiotti, L.;
Bajpai, R.;
Barton, M. A.;
Cannon, K.;
Capocasa, E.;
Chan, M.;
Chen, C.;
Chen, K.;
Chen, Y.;
Chu, H.;
Chu, Y-K.;
Doi, K.;
Eguchi, S.;
Enomoto, Y.;
Flaminio, R.;
Fujii, Y.;
Fukunaga, M.;
Fukushima, M.;
Ge, G.;
Hagiwara, A.;
Haino, S.;
Hasegawa, K.;
Hayakawa, H.;
Hayama, K.;
Himemoto, Y.;
Hiranuma, Y.;
Hirata, N.;
Hirose, E.;
Hong, Z.;
Hsieh, B. H.;
Huang, G-Z.;
Huang, P.;
Huang, Y.;
Ikenoue, B.;
Imam, S.;
Inayoshi, K.;
Inoue, Y.;
Ioka, K.;
Itoh, Y.;
Izumi, K.;
Jung, K.;
Jung, P.;
Kajita, T.;
Kamiizumi, M.;
Kanbara, S.;
Kanda, N.;
Kang, G.;
Kawaguchi, K.;
Kawai, N.;
Kawasaki, T.;
Kim, C.;
Kim, J. C.;
Kim, W. S.;
Kim, Y. -M.;
Kimura, N.;
Kita, N.;
Kitazawa, H.;
Kojima, Y.;
Kokeyama, K.;
Komori, K.;
Kong, A. K. H.;
Kotake, K.;
Kozakai, C.;
Kozu, R.;
Kumar, R.;
Kume, J.;
Kuo, C.;
Kuo, H-S.;
Kuroyanagi, S.;
Kusayanagi, K.;
Kwak, K.;
Lee, H. K.;
Lee, H. W.;
Lee, R.;
Leonardi, M.;
Lin, L. C. -C.;
Lin, C-Y.;
Lin, F-L.;
Liu, G. C.;
Luo, L. -W.;
Marchio, M.;
Michimura, Y.;
Mio, N.;
Miyakawa, O.;
Miyamoto, A.;
Miyazaki, Y.;
Miyo, K.;
Miyoki, S.;
Morisaki, S.;
Moriwaki, Y.;
Musha, M.;
Nagano, K.;
Nagano, S.;
Nakamura, K.;
Nakano, H.;
Nakano, M.;
Nakashima, R.;
Narikawa, T.;
Negishi, R.;
Ni, W. -T.;
Nishizawa, A.;
Obuchi, Y.;
Ogaki, W.;
Oh, J. J.;
Oh, S. H.;
Ohashi, M.;
Ohishi, N.;
Ohkawa, M.;
Ohmae, N.;
Okutomi, K.;
Oohara, K.;
Ooi, C. P.;
Oshino, S.;
Pan, K.;
Pang, H.;
Park, J.;
Arellano, F. E. Peña;
Pinto, I.;
Sago, N.;
Saito, S.;
Saito, Y.;
Sakai, K.;
Sakai, Y.;
Sakuno, Y.;
Sato, S.;
Sato, T.;
Sawada, T.;
Sekiguchi, T.;
Sekiguchi, Y.;
Shibagaki, S.;
Shimizu, R.;
Shimoda, T.;
Shimode, K.;
Shinkai, H.;
Shishido, T.;
Shoda, A.;
Somiya, K.;
Son, E. J.;
Sotani, H.;
Sugimoto, R.;
Suzuki, T.;
Suzuki, T.;
Tagoshi, H.;
Takahashi, H.;
Takahashi, R.;
Takamori, A.;
Takano, S.;
Takeda, H.;
Takeda, M.;
Tanaka, H.;
Tanaka, K.;
Tanaka, K.;
Tanaka, T.;
Tanaka, T.;
Tanioka, S.;
Martin, E. N. Tapia San;
Tatsumi, D.;
Telada, S.;
Tomaru, T.;
Tomigami, Y.;
Tomura, T.;
Travasso, F.;
Trozzo, L.;
Tsang, T.;
Tsubono, K.;
Tsuchida, S.;
Tsuzuki, T.;
Tuyenbayev, D.;
Uchikata, N.;
Uchiyama, T.;
Ueda, A.;
Uehara, T.;
Ueno, K.;
Ueshima, G.;
Uraguchi, F.;
Ushiba, T.;
van Putten, M. H. P. M.;
Vocca, H.;
Wang, J.;
Wu, C.;
Wu, H.;
Wu, S.;
Xu, W-R.;
Yamada, T.;
Yamamoto, K.;
Yamamoto, K.;
Yamamoto, T.;
Yokogawa, K.;
Yokoyama, J.;
Yokozawa, T.;
Yoshioka, T.;
Yuzurihara, H.;
Zeidler, S.;
Zhao, Y.;
Zhu, Z. -H.
Submitted: 2019-10-02, last modified: 2019-11-28
Modern ground-based gravitational wave (GW) detectors require a complex
interferometer configuration with multiple coupled optical cavities. Since
achieving the resonances of the arm cavities is the most challenging among the
lock acquisition processes, the scheme called arm length stabilization (ALS)
had been employed for lock acquisition of the arm cavities. We designed a new
type of the ALS, which is compatible with the interferometers having long arms
like the next generation GW detectors. The features of the new ALS are that the
control configuration is simpler than those of previous ones and that it is not
necessary to lay optical fibers for the ALS along the kilometer-long arms of
the detector. Along with simulations of its noise performance, an experimental
test of the new ALS was performed utilizing a single arm cavity of KAGRA. This
paper presents the first results of the test where we demonstrated that lock
acquisition of the arm cavity was achieved using the new ALS and residual noise
was measured to be $8.2\,\mathrm{Hz}$ in units of frequency, which is smaller
than the linewidth of the arm cavity and thus low enough to lock the full
interferometer of KAGRA in a repeatable and reliable manner.
[12]
oai:arXiv.org:1904.08088 [pdf] - 1983844
Long-term Simulations of Multi-Dimensional Core-collapse Supernovae:
Implications for Neutron Star Kicks
Submitted: 2019-04-17, last modified: 2019-08-05
Core-collapse supernovae (CCSNe) are the final stage of massive stars,
marking the birth of neutron stars (NSs). The aspherical mass ejection drives a
natal kick of the forming NS. In this work we study the properties of the NS
kick based on our long-term hydrodynamics CCSN simulations. We perform
two-dimensional (2D) simulations for ten progenitors from a 10.8 to 20 $\,
M_{\odot}$ star covering a wide range of the progenitor's compactness
parameter, and two three-dimensional (3D) simulations for an 11.2 $\,
M_{\odot}$ star. Our 2D models present a variety of explosion energies between
$\sim 1.3 \times 10^{50}$ erg and $\sim 1.2 \times 10^{51}$ erg, and NS kick
velocities between $\sim 100$ km s$^{-1}$ and $\sim 1500$ km s$^{-1}$. For the
2D exploding models, we find that the kick velocities tend to become higher
with the progenitor's compactness. This is because the high progenitor
compactness results in high neutrino luminosity from the proto-neutron star
(PNS), leading to more energetic explosions. Since high-compactness progenitors
produce massive PNSs, we point out that the NS masses and the kick velocities
can be correlated, which is moderately supported by observation. Comparing 2D
and 3D models of the 11.2 $\, M_{\odot}$ star, the diagnostic explosion energy
in 3D is, as previously identified, higher than that in 2D, whereas the 3D
model results in a smaller asymmetry in the ejecta distribution and a smaller
kick velocity than in 2D. Our results confirm the importance of self-consistent
CCSN modeling covering a long-term postbounce evolution in 3D for a
quantitative prediction of the NS kicks.
[13]
oai:arXiv.org:1903.07811 [pdf] - 1938370
One-, Two-, and Three-dimensional Simulations of Oxygen Shell Burning
Just Before the Core-Collapse of Massive Stars
Submitted: 2019-03-18, last modified: 2019-06-19
We perform two- (2D) and three-dimensional (3D) hydrodynamics simulations of
convective oxygen shell-burning that takes place deep inside a massive
progenitor star of a core-collapse supernova. Using one dimensional (1D)
stellar evolution code, we first calculate the evolution of massive stars with
an initial mass of 9-40 $M_\odot$. Four different overshoot parameters are
applied, and CO core mass trend similar to previous works is obtained in the 1D
models. Selecting eleven 1D models that have a silicon and oxygen coexisting
layer, we perform 2D hydrodynamics simulations of the evolution $\sim$100 s
until the onset of core-collapse. We find that convection with large-scale
eddies and the turbulent Mach number $\sim$0.1 is obtained in the models having
a Si/O layer with a scale of 10$^8$ cm, whereas most models that have an
extended O/Si layer up to a few $\times 10^9$ cm exhibit lower turbulent
velocity. Our results indicate that the supernova progenitors that possess a
thick Si/O layer could provide a preferable condition for perturbation-aided
explosions. We perform 3D simulation of a 25 $M_\odot$ model, which exhibits
large-scale convection in the 2D models. The 3D model develops large ($\ell =
2$) convection similar to the 2D model, however, the turbulent velocity is
lower. By estimating the neutrino emission properties of the 3D model, we point
out that a time modulation of the event rates, if observed in KamLAND and
Hyper-Kamiokande, would provide an important information about structural
changes in the presupernova convective layer.
[14]
oai:arXiv.org:1906.04354 [pdf] - 1908729
Dependence of outer boundary condition on protoneutron star
asteroseismology with gravitational-wave signatures
Submitted: 2019-06-10
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.
[15]
oai:arXiv.org:1903.09224 [pdf] - 1854203
The Yet-Unobserved Multi-Messenger Gravitational-Wave Universe
Kalogera, Vassiliki;
Bizouard, Marrie-Anne;
Burrows, Adam;
Janka, Thomas;
Kotake, Kei;
Messer, Bronson;
Mezzacappa, Tony;
Mueller, Bernhard;
Mueller, Ewald;
Papa, Maria Alessandra;
Reddy, Sanjay;
Rosswog, Stephan
Submitted: 2019-03-21
Observations with next-generation ground-based detectors further enhanced
with multi-messenger (electromagnetic and neutrino) detections will allow us to
probe new extreme astrophysics. Target sources included: core-collapse
supernovae, continuous emission from isolated or accreting neutron stars, and
bursts from magnetars and other pulsars.
[16]
oai:arXiv.org:1901.10692 [pdf] - 1824490
Prototype Implementation of a Web-Based Gravitational Wave Signal
Analyzer: SNEGRAF
Submitted: 2019-01-30
A direct detection of gravitational waves is one of the most exciting
frontiers for modern astronomy and astrophysics. Gravitational wave signals
combined with classical electro-magnetic observations, known as multi-messenger
astronomy, promise newer and deeper insights about the cosmic evolution of
astrophysical objects such as neutron starts and black holes. To this end, we
have been developing an original data processing pipeline for KAGRA, a Japanese
gravitational wave telescope, for optimal detections of supernova events. As a
part of our project, we released a web application named SuperNova Event
Gravitational-wave-display in Fukuoka (SNEGRAF) in autumn 2018. SNEGRAF accepts
the users' theoretical waveforms as a plain text file consisting of a time
series of $h_{+}$ and $h_{\times}$ (the plus and cross mode of gravitational
waves, respectively), then displays the input, a corresponding spectrogram, and
power spectrum together with KAGRA sensitivity curve and the signal-to-noise
ratio; we adopt Google Visualization API for the interactive visualization of
the input waveforms. However, it is a time-consuming task to draw more than
$\sim 10^{5}$ data points directly with JavaScript, although the number can be
typical for a supernova hunt by assuming a typical duration of the event and
sampling rate of the detectors; a combination of recursive decimations of the
original in the server-side program and an appropriate selection of them
depending on the time duration requested by the user in a web browser achieves
an acceptable latency. In this paper, we present the current design,
implementation and optimization algorithms of SNEGRAF, and its future
perspectives.
[17]
oai:arXiv.org:1901.03569 [pdf] - 1979528
First cryogenic test operation of underground km-scale
gravitational-wave observatory KAGRA
KAGRA Collaboration;
Akutsu, T.;
Ando, M.;
Arai, K.;
Arai, Y.;
Araki, S.;
Araya, A.;
Aritomi, N.;
Asada, H.;
Aso, Y.;
Atsuta, S.;
Awai, K.;
Bae, S.;
Baiotti, L.;
Barton, M. A.;
Cannon, K.;
Capocasa, E.;
Chen, C-S.;
Chiu, T-W.;
Cho, K.;
Chu, Y-K.;
Craig, K.;
Creus, W.;
Doi, K.;
Eda, K.;
Enomoto, Y.;
Flaminio, R.;
Fujii, Y.;
Fujimoto, M. -K.;
Fukunaga, M.;
Fukushima, M.;
Furuhata, T.;
Hagiwara, A.;
Haino, S.;
Hasegawa, K.;
Hashino, K.;
Hayama, K.;
Hirobayashi, S.;
Hirose, E.;
Hsieh, B. H.;
Huang, C-Z.;
Ikenoue, B.;
Inoue, Y.;
Ioka, K.;
Itoh, Y.;
Izumi, K.;
Kaji, T.;
Kajita, T.;
Kakizaki, M.;
Kamiizumi, M.;
Kanbara, S.;
Kanda, N.;
Kanemura, S.;
Kaneyama, M.;
Kang, G.;
Kasuya, J.;
Kataoka, Y.;
Kawai, N.;
Kawamura, S.;
Kawasaki, T.;
Kim, C.;
Kim, J.;
Kim, J. C.;
Kim, W. S.;
Kim, Y. -M.;
Kimura, N.;
Kinugawa, T.;
Kirii, S.;
Kitaoka, Y.;
Kitazawa, H.;
Kojima, Y.;
Kokeyama, K.;
Komori, K.;
Kong, A. K. H.;
Kotake, K.;
Kozu, R.;
Kumar, R.;
Kuo, H-S.;
Kuroyanagi, S.;
Lee, H. K.;
Lee, H. M.;
Lee, H. W.;
Leonardi, M.;
Lin, C-Y.;
Lin, F-L.;
Liu, G. C.;
Liu, Y.;
Majorana, E.;
Mano, S.;
Marchio, M.;
Matsui, T.;
Matsushima, F.;
Michimura, Y.;
Mio, N.;
Miyakawa, O.;
Miyamoto, A.;
Miyamoto, T.;
Miyo, K.;
Miyoki, S.;
Morii, W.;
Morisaki, S.;
Moriwaki, Y.;
Morozumi, T.;
Murakami, I.;
Musha, M.;
Nagano, K.;
Nagano, S.;
Nakamura, K.;
Nakamura, T.;
Nakano, H.;
Nakano, M.;
Nakao, K.;
Namai, Y.;
Narikawa, T.;
Naticchioni, L.;
Quynh, L. Nguyen;
Ni, W. -T.;
Nishizawa, A.;
Obuchi, Y.;
Ochi, T.;
Oh, J. J.;
Oh, S. H.;
Ohashi, M.;
Ohishi, N.;
Ohkawa, M.;
Okutomi, K.;
Ono, K.;
Oohara, K.;
Ooi, C. P.;
Pan, S-S.;
Park, J.;
Arellano, F. E. Peña;
Pinto, I.;
Sago, N.;
Saijo, M.;
Saito, Y.;
Saitou, S.;
Sakai, K.;
Sakai, Y.;
Sakai, Y.;
Sasai, M.;
Sasaki, M.;
Sasaki, Y.;
Sato, N.;
Sato, S.;
Sato, T.;
Sekiguchi, Y.;
Seto, N.;
Shibata, M.;
Shimoda, T.;
Shinkai, H.;
Shishido, T.;
Shoda, A.;
Somiya, K.;
Son, E. J.;
Suemasa, A.;
Suzuki, T.;
Suzuki, T.;
Tagoshi, H.;
Tahara, H.;
Takahashi, H.;
Takahashi, R.;
Takamori, A.;
Takeda, H.;
Tanaka, H.;
Tanaka, K.;
Tanaka, T.;
Tanioka, S.;
Martin, E. N. Tapia San;
Tatsumi, D.;
Terashima, S.;
Tomaru, T.;
Tomura, T.;
Travasso, F.;
Tsubono, K.;
Tsuchida, S.;
Uchikata, N.;
Uchiyama, T.;
Ueda, A.;
Uehara, T.;
Ueki, S.;
Ueno, K.;
Uraguchi, F.;
Ushiba, T.;
van Putten, M. H. P. M.;
Vocca, H.;
Wada, S.;
Wakamatsu, T.;
Watanabe, Y.;
Xu, W-R.;
Yamada, T.;
Yamamoto, A.;
Yamamoto, K.;
Yamamoto, K.;
Yamamoto, S.;
Yamamoto, T.;
Yokogawa, K.;
Yokoyama, J.;
Yokozawa, T.;
Yoon, T. H.;
Yoshioka, T.;
Yuzurihara, H.;
Zeidler, S.;
Zhu, Z. -H.
Submitted: 2019-01-11
KAGRA is a second-generation interferometric gravitational-wave detector with
3-km arms constructed at Kamioka, Gifu in Japan. It is now in its final
installation phase, which we call bKAGRA (baseline KAGRA), with scientific
observations expected to begin in late 2019. One of the advantages of KAGRA is
its underground location of at least 200 m below the ground surface, which
brings small seismic motion at low frequencies and high stability of the
detector. Another advantage is that it cools down the sapphire test mass
mirrors to cryogenic temperatures to reduce thermal noise. In April-May 2018,
we have operated a 3-km Michelson interferometer with a cryogenic test mass for
10 days, which was the first time that km-scale interferometer was operated at
cryogenic temperatures. In this article, we report the results of this "bKAGRA
Phase 1" operation. We have demonstrated the feasibility of 3-km interferometer
alignment and control with cryogenic mirrors.
[18]
oai:arXiv.org:1811.08079 [pdf] - 1979518
KAGRA: 2.5 Generation Interferometric Gravitational Wave Detector
Akutsu, T.;
Ando, M.;
Arai, K.;
Arai, Y.;
Araki, S.;
Araya, A.;
Aritomi, N.;
Asada, H.;
Aso, Y.;
Atsuta, S.;
Awai, K.;
Bae, S.;
Baiotti, L.;
Barton, M. A.;
Cannon, K.;
Capocasa, E.;
Chen, C-S.;
Chiu, T-W.;
Cho, K.;
Chu, Y-K.;
Craig, K.;
Creus, W.;
Doi, K.;
Eda, K.;
Enomoto, Y.;
Flaminio, R.;
Fujii, Y.;
Fujimoto, M. -K.;
Fukunaga, M.;
Fukushima, M.;
Furuhata, T.;
Haino, S.;
Hasegawa, K.;
Hashino, K.;
Hayama, K.;
Hirobayashi, S.;
Hirose, E.;
Hsieh, B. H.;
Huang, C-Z.;
Ikenoue, B.;
Inoue, Y.;
Ioka, K.;
Itoh, Y.;
Izumi, K.;
Kaji, T.;
Kajita, T.;
Kakizaki, M.;
Kamiizumi, M.;
Kanbara, S.;
Kanda, N.;
Kanemura, S.;
Kaneyama, M.;
Kang, G.;
Kasuya, J.;
Kataoka, Y.;
Kawai, N.;
Kawamura, S.;
Kawasaki, T.;
Kim, C.;
Kim, J.;
Kim, J. C.;
Kim, W. S.;
Kim, Y. -M.;
Kimura, N.;
Kinugawa, T.;
Kirii, S.;
Kitaoka, Y.;
Kitazawa, H.;
Kojima, Y.;
Kokeyama, K.;
Komori, K.;
Kong, A. K. H.;
Kotake, K.;
Kozu, R.;
Kumar, R.;
Kuo, H-S.;
Kuroyanagi, S.;
Lee, H. K.;
Lee, H. M.;
Lee, H. W.;
Leonardi, M.;
Lin, C-Y.;
Lin, F-L.;
Liu, G. C.;
Liu, Y.;
Majorana, E.;
Mano, S.;
Marchio, M.;
Matsui, T.;
Matsushima, F.;
Michimura, Y.;
Mio, N.;
Miyakawa, O.;
Miyamoto, A.;
Miyamoto, T.;
Miyo, K.;
Miyoki, S.;
Morii, W.;
Morisaki, S.;
Moriwaki, Y.;
Morozumi, T.;
Musha, M.;
Nagano, K.;
Nagano, S.;
Nakamura, K.;
Nakamura, T.;
Nakano, H.;
Nakano, M.;
Nakao, K.;
Narikawa, T.;
Naticchioni, L.;
Quynh, L. Nguyen;
Ni, W. -T.;
Nishizawa, A.;
Ochi, T.;
Oh, J. J.;
Oh, S. H.;
Ohashi, M.;
Ohishi, N.;
Ohkawa, M.;
Okutomi, K.;
Ono, K.;
Oohara, K.;
Ooi, C. P.;
Pan, S-S.;
Park, J.;
Arellano, F. E. Peña;
Pinto, I.;
Sago, N.;
Saijo, M.;
Saito, Y.;
Sakai, K.;
Sakai, Y.;
Sakai, Y.;
Sasai, M.;
Sasaki, M.;
Sasaki, Y.;
Sato, S.;
Sato, T.;
Sekiguchi, Y.;
Seto, N.;
Shibata, M.;
Shimoda, T.;
Shinkai, H.;
Shishido, T.;
Shoda, A.;
Somiya, K.;
Son, E. J.;
Suemasa, A.;
Suzuki, T.;
Suzuki, T.;
Tagoshi, H.;
Tahara, H.;
Takahashi, H.;
Takahashi, R.;
Takamori, A.;
Takeda, H.;
Tanaka, H.;
Tanaka, K.;
Tanaka, T.;
Tanioka, S.;
Martin, E. N. Tapia San;
Tatsumi, D.;
Tomaru, T.;
Tomura, T.;
Travasso, F.;
Tsubono, K.;
Tsuchida, S.;
Uchikata, N.;
Uchiyama, T.;
Uehara, T.;
Ueki, S.;
Ueno, K.;
Ushiba, T.;
van Putten, M. H. P. M.;
Vocca, H.;
Wada, S.;
Wakamatsu, T.;
Watanabe, Y.;
Xu, W-R.;
Yamada, T.;
Yamamoto, A.;
Yamamoto, K.;
Yamamoto, K.;
Yamamoto, S.;
Yamamoto, T.;
Yokogawa, K.;
Yokoyama, J.;
Yokozawa, T.;
Yoon, T. H.;
Yoshioka, T.;
Yuzurihara, H.;
Zeidler, S.;
Zhu, Z. -H.
Submitted: 2018-11-20
The recent detections of gravitational waves (GWs) reported by LIGO/Virgo
collaborations have made significant impact on physics and astronomy. A global
network of GW detectors will play a key role to solve the unknown nature of the
sources in coordinated observations with astronomical telescopes and detectors.
Here we introduce KAGRA (former name LCGT; Large-scale Cryogenic Gravitational
wave Telescope), a new GW detector with two 3-km baseline arms arranged in the
shape of an "L", located inside the Mt. Ikenoyama, Kamioka, Gifu, Japan.
KAGRA's design is similar to those of the second generations such as Advanced
LIGO/Virgo, but it will be operating at the cryogenic temperature with sapphire
mirrors. This low temperature feature is advantageous for improving the
sensitivity around 100 Hz and is considered as an important feature for the
third generation GW detector concept (e.g. Einstein Telescope of Europe or
Cosmic Explorer of USA). Hence, KAGRA is often called as a 2.5 generation GW
detector based on laser interferometry. The installation and commissioning of
KAGRA is underway and its cryogenic systems have been successfully tested in
May, 2018. KAGRA's first observation run is scheduled in late 2019, aiming to
join the third observation run (O3) of the advanced LIGO/Virgo network. In this
work, we describe a brief history of KAGRA and highlights of main feature. We
also discuss the prospects of GW observation with KAGRA in the era of O3. When
operating along with the existing GW detectors, KAGRA will be helpful to locate
a GW source more accurately and to determine the source parameters with higher
precision, providing information for follow-up observations of a GW trigger
candidate.
[19]
oai:arXiv.org:1805.10419 [pdf] - 1771554
Chiral magnetohydrodynamic turbulence in core-collapse supernovae
Submitted: 2018-05-25, last modified: 2018-10-15
Macroscopic evolution of relativistic charged matter with chirality imbalance
is described by the chiral magnetohydrodynamics (chiral MHD). One such
astrophysical system is high-density lepton matter in core-collapse supernovae
where the chirality imbalance of leptons is generated by the parity-violating
weak processes. After developing the chiral MHD equations for this system, we
perform numerical simulations for the real-time evolutions of magnetic and flow
fields, and study the properties of the chiral MHD turbulence. In particular,
we observe the inverse cascade of the magnetic energy and the fluid kinetic
energy. Our results suggest that the chiral effects that have been neglected so
far can reverse the turbulent cascade direction from direct to inverse cascade,
which would impact the magnetohydrodynamics evolution in the supernova core
toward explosion.
[20]
oai:arXiv.org:1810.00334 [pdf] - 1783938
A Linear and Quadratic Time-Frequency Analysis of Gravitational Waves
from Core-Collapse Supernovae
Submitted: 2018-09-30
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.
[21]
oai:arXiv.org:1806.04175 [pdf] - 1755942
Global Comparison of Core-Collapse Supernova Simulations in Spherical
Symmetry
O'Connor, Evan;
Bollig, Robert;
Burrows, Adam;
Couch, Sean;
Fischer, Tobias;
Janka, Hans-Thomas;
Kotake, Kei;
Lentz, Eric J.;
Liebendörfer, Matthias;
Messer, O. E. Bronson;
Mezzacappa, Anthony;
Takiwaki, Tomoya;
Vartanyan, David
Submitted: 2018-06-11
We present a comparison between several simulation codes designed to study
the core-collapse supernova mechanism. We pay close attention to controlling
the initial conditions and input physics in order to ensure a meaningful and
informative comparison. Our goal is three-fold. First, we aim to demonstrate
the current level of agreement between various groups studying the
core-collapse supernova central engine. Second, we desire to form a strong
basis for future simulation codes and methods to compare to. Lastly, we want
this work to be a stepping stone for future work exploring more complex
simulations of core-collapse supernovae, i.e., simulations in multiple
dimensions and simulations with modern neutrino and nuclear physics. We compare
the early (first ~500ms after core bounce) spherically-symmetric evolution of a
20 solar mass progenitor star from six different core-collapse supernovae
codes: 3DnSNe-IDSA, AGILE-BOLTZTRAN, FLASH, F{\sc{ornax}}, GR1D, and
PROMETHEUS-VERTEX. Given the diversity of neutrino transport and hydrodynamic
methods employed, we find excellent agreement in many critical quantities,
including the shock radius evolution and the amount of neutrino heating. Our
results provide an excellent starting point from which to extend this
comparison to higher dimensions and compare the development of hydrodynamic
instabilities that are crucial to the supernova explosion mechanism, such as
turbulence and convection.
[22]
oai:arXiv.org:1805.00050 [pdf] - 1675131
Diagnosing the Structure of Massive Stars with Galactic Supernova
Neutrinos
Submitted: 2018-04-30
It has been suggested that whether a star explodes or not, and what kind of
explosion properties it shows, is strongly dependent on the progenitor's core
structure. We present the results from 101 axisymmetric core-collapse supernova
simulations performed with progenitors spanning initial masses in the range
from 10.8 to 75 solar masses, and focus on their connections to the compactness
of the progenitor's core. Our simulations confirm a correlation between the
neutrinos emitted during the accretion phase and the progenitor's compactness.
We suggest that the ratio of observed neutrino events during the first hundreds
of milliseconds can be used to infer the progenitor's inner mass density
structure.
[23]
oai:arXiv.org:1801.01293 [pdf] - 1682539
A full general relativistic neutrino radiation-hydrodynamics simulation
of a collapsing very massive star and the formation of a black hole
Submitted: 2018-01-04, last modified: 2018-04-04
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.
[24]
oai:arXiv.org:1710.10372 [pdf] - 1659567
Intermediate-Mass-Elements in Young Supernova Remnants Reveal Neutron
Star Kicks by Asymmetric Explosions
Katsuda, Satoru;
Morii, Mikio;
Janka, Hans-Thomas;
Wongwathanarat, Annop;
Nakamura, Ko;
Kotake, Kei;
Mori, Koji;
Müller, Ewald;
Takiwaki, Tomoya;
Tanaka, Masaomi;
Tominaga, Nozomu;
Tsunemi, Hiroshi
Submitted: 2017-10-27, last modified: 2018-02-18
The birth properties of neutron stars yield important information on the
still debated physical processes that trigger the explosion and on intrinsic
neutron-star physics. These properties include the high space velocities of
young neutron stars with average values of several 100 km/s, whose underlying
"kick" mechanism is not finally clarified. There are two competing
possibilities that could accelerate NSs during their birth: anisotropic
ejection of either stellar debris or neutrinos. We here present new evidence
from X-ray measurements that chemical elements between silicon and calcium in
six young gaseous supernova remnants are preferentially expelled opposite to
the direction of neutron star motion. There is no correlation between the kick
velocities and magnetic field strengths of these neutron stars. Our results
support a hydrodynamic origin of neutron-star kicks connected to asymmetric
explosive mass ejection, and they conflict with neutron-star acceleration
scenarios that invoke anisotropic neutrino emission caused by particle and
nuclear physics in combination with very strong neutron-star magnetic fields.
[25]
oai:arXiv.org:1802.03842 [pdf] - 1682576
Circular polarization of gravitational waves from non-rotating supernova
cores: a new probe into the pre-explosion hydrodynamics
Submitted: 2018-02-11
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.
[26]
oai:arXiv.org:1801.02703 [pdf] - 1631472
Impact of Neutrino Opacities on Core-Collapse Supernova Simulations
Submitted: 2018-01-08, last modified: 2018-02-07
Accurate description of neutrino opacities is central both to the
core-collapse supernova (CCSN) phenomenon and to the validity of the explosion
mechanism itself. In this work, we study in a systematic fashion the role of a
variety of well-selected neutrino opacities in CCSN simulations where
multi-energy, three-flavor neutrino transport is solved by the isotropic
diffusion source approximation (IDSA) scheme. To verify our code, we first
present results from one-dimensional (1D) simulations following core-collapse,
bounce, and up to ~ 250 ms postbounce of a 15 $M_{\odot}$ star using a standard
set of neutrino opacities by Bruenn (1985). Detailed comparison with published
results supports the reliability of our three-flavor IDSA scheme using the
standard opacity set. We then investigate in 1D simulations how the individual
opacity update leads to the difference from the base-line run with the standard
opacity set. By making a detailed comparison with previous work, we check the
validity of our implementation of each update in a step-by-step manner.
Individual neutrino opacities with the largest impact on the overall evolution
in 1D simulations are selected for a systematic comparison in our
two-dimensional (2D) simulations. Special emphasis is devoted to the criterion
of explodability in the 2D models. We discuss the implications of these results
as well as the limitations and requirements for future towards more elaborate
CCSN modeling.
[27]
oai:arXiv.org:1709.06567 [pdf] - 1630189
Diffuse Supernova Neutrino Background from extensive core-collapse
simulations of $8$-$100 {\rm M}_\odot$ progenitors
Submitted: 2017-09-19, last modified: 2018-01-15
We revisit the diffuse supernova neutrino background in light of recent
systematic studies of stellar core collapse that reveal the quantitative
impacts of the progenitor conditions on the collapse process. In general, the
dependence of the progenitor on the core-collapse neutrino emission is not
monotonic in progenitor initial mass, but we show that it can, at first order,
be characterized by the core compactness. For the first time, we incorporate
the detailed variations in the neutrino emission over the entire mass range
$8$-$100 {\rm M}_\odot$, based on (i) a long-term simulation of the core
collapse of a $8.8 {\rm M}_\odot$ O-Ne-Mg core progenitor, (ii) over 100
simulations of iron core collapse to neutron stars, and (iii) half a dozen
simulations of core collapse to black holes (the "failed channel"). The
fraction of massive stars that undergo the failed channel remains uncertain,
but in view of recent simulations which reveal high compactness to be conducive
to collapse to black holes, we characterize the failed fraction by considering
a threshold compactness above which massive stars collapse to black holes and
below which the final remnant is a neutron star. We predict that future
detections of the diffuse supernova neutrino background may have the power to
reveal this threshold compactness, if its value is relatively small as
suggested by interpretations of several recent astronomical observations.
[28]
oai:arXiv.org:1711.01905 [pdf] - 1622320
Anisotropic emission of neutrino and gravitational-wave signals from
rapidly rotating core-collapse supernovae
Submitted: 2017-11-06, last modified: 2018-01-06
We present analysis of neutrino and GW signals based on three-dimensional
(3D) core-collapse supernova simulations of a rapidly rotating 27 $M_{\odot}$
star. We find a new neutrino signature that is produced by a lighthouse effect
where the spinning of strong neutrino emission regions around the rotational
axis leads to quasi-periodic modulation in the neutrino signal. Depending on
the observer's viewing angle, the time modulation will be clearly detectable in
IceCube and the future Hyper-Kamiokande. The GW emission is also anisotropic
where the GW signal is emitted, as previously identified, most strongly toward
the equator at rotating core-collapse and bounce, and the non-axisymmetric
instabilities in the postbounce phase lead to stronger GW emission toward the
spin axis. We show that these GW signals can be a target of LIGO-class
detectors for a Galactic event. The origin of the postbounce GW emission
naturally explains why the peak GW frequency is about twice of the neutrino
modulation frequency. We point out that the simultaneous detection of the
rotation-induced neutrino and GW signatures could provide a smoking-gun
signature of a rapidly rotating proto-neutron star at the birth.
[29]
oai:arXiv.org:1712.00148 [pdf] - 2127370
Construction of KAGRA: an Underground Gravitational Wave Observatory
Akutsu, T.;
Ando, M.;
Araki, S.;
Araya, A.;
Arima, T.;
Aritomi, N.;
Asada, H.;
Aso, Y.;
Atsuta, S.;
Awai, K.;
Baiotti, L.;
Barton, M. A.;
Chen, D.;
Cho, K.;
Craig, K.;
DeSalvo, R.;
Doi, K.;
Eda, K.;
Enomoto, Y.;
Flaminio, R.;
Fujibayashi, S.;
Fujii, Y.;
Fujimoto, M. -K.;
Fukushima, M.;
Furuhata, T.;
Hagiwara, A.;
Haino, S.;
Harita, S.;
Hasegawa, K.;
Hasegawa, M.;
Hashino, K.;
Hayama, K.;
Hirata, N.;
Hirose, E.;
Ikenoue, B.;
Inoue, Y.;
Ioka, K.;
Ishizaki, H.;
Itoh, Y.;
Jia, D.;
Kagawa, T.;
Kaji, T.;
Kajita, T.;
Kakizaki, M.;
Kakuhata, H.;
Kamiizumi, M.;
Kanbara, S.;
Kanda, N.;
Kanemura, S.;
Kaneyama, M.;
Kasuya, J.;
Kataoka, Y.;
Kawaguchi, K.;
Kawai, N.;
Kawamura, S.;
Kawazoe, F.;
Kim, C.;
Kim, J.;
Kim, J. C.;
Kim, W.;
Kimura, N.;
Kitaoka, Y.;
Kobayashi, K.;
Kojima, Y.;
Kokeyama, K.;
Komori, K.;
Kotake, K.;
Kubo, K.;
Kumar, R.;
Kume, T.;
Kuroda, K.;
Kuwahara, Y.;
Lee, H. -K.;
Lee, H. -W.;
Lin, C. -Y.;
Liu, Y.;
Majorana, E.;
Mano, S.;
Marchio, M.;
Matsui, T.;
Matsumoto, N.;
Matsushima, F.;
Michimura, Y.;
Mio, N.;
Miyakawa, O.;
Miyake, K.;
Miyamoto, A.;
Miyamoto, T.;
Miyo, K.;
Miyoki, S.;
Morii, W.;
Morisaki, S.;
Moriwaki, Y.;
Muraki, Y.;
Murakoshi, M.;
Musha, M.;
Nagano, K.;
Nagano, S.;
Nakamura, K.;
Nakamura, T.;
Nakano, H.;
Nakano, M.;
Nakano, M.;
Nakao, H.;
Nakao, K.;
Narikawa, T.;
Ni, W. -T.;
Nonomura, T.;
Obuchi, Y.;
Oh, J. J.;
Oh, S. -H.;
Ohashi, M.;
Ohishi, N.;
Ohkawa, M.;
Ohmae, N.;
Okino, K.;
Okutomi, K.;
Ono, K.;
Ono, Y.;
Oohara, K.;
Ota, S.;
Park, J.;
Arellano, F. E. Peña;
Pinto, I. M.;
Principe, M.;
Sago, N.;
Saijo, M.;
Saito, T.;
Saito, Y.;
Saitou, S.;
Sakai, K.;
Sakakibara, Y.;
Sasaki, Y.;
Sato, S.;
Sato, T.;
Sato, Y.;
Sekiguchi, T.;
Sekiguchi, Y.;
Shibata, M.;
Shiga, K.;
Shikano, Y.;
Shimoda, T.;
Shinkai, H.;
Shoda, A.;
Someya, N.;
Somiya, K.;
Son, E. J.;
Starecki, T.;
Suemasa, A.;
Sugimoto, Y.;
Susa, Y.;
Suwabe, H.;
Suzuki, T.;
Tachibana, Y.;
Tagoshi, H.;
Takada, S.;
Takahashi, H.;
Takahashi, R.;
Takamori, A.;
Takeda, H.;
Tanaka, H.;
Tanaka, K.;
Tanaka, T.;
Tatsumi, D.;
Telada, S.;
Tomaru, T.;
Tsubono, K.;
Tsuchida, S.;
Tsukada, L.;
Tsuzuki, T.;
Uchikata, N.;
Uchiyama, T.;
Uehara, T.;
Ueki, S.;
Ueno, K.;
Uraguchi, F.;
Ushiba, T.;
van Putten, M. H. P. M.;
Wada, S.;
Wakamatsu, T.;
Yaginuma, T.;
Yamamoto, K.;
Yamamoto, S.;
Yamamoto, T.;
Yano, K.;
Yokoyama, J.;
Yokozawa, T.;
Yoon, T. H.;
Yuzurihara, H.;
Zeidler, S.;
Zhao, Y.;
Zheng, L.;
Agatsuma, K.;
Akiyama, Y.;
Arai, N.;
Asano, M.;
Bertolini, A.;
Fujisawa, M.;
Goetz, R.;
Guscott, J.;
Hashimoto, Y.;
Hayashida, Y.;
Hennes, E.;
Hirai, K.;
Hirayama, T.;
Ishitsuka, H.;
Kato, J.;
Khalaidovski, A.;
Koike, S.;
Kumeta, A.;
Miener, T.;
Morioka, M.;
Mueller, C. L.;
Narita, T.;
Oda, Y.;
Ogawa, T.;
Sekiguchi, T.;
Tamura, H.;
Tanner, D. B.;
Tokoku, C.;
Toritani, M.;
Utsuki, T.;
Uyeshima, M.;
Brand, J. van den;
van Heijningen, J.;
Yamaguchi, S.;
Yanagida, A.
Submitted: 2017-11-30, last modified: 2017-12-11
Major construction and initial-phase operation of a second-generation
gravitational-wave detector KAGRA has been completed. The entire 3-km detector
is installed underground in a mine in order to be isolated from background
seismic vibrations on the surface. This allows us to achieve a good sensitivity
at low frequencies and high stability of the detector. Bare-bones equipment for
the interferometer operation has been installed and the first test run was
accomplished in March and April of 2016 with a rather simple configuration. The
initial configuration of KAGRA is named {\it iKAGRA}. In this paper, we
summarize the construction of KAGRA, including the study of the advantages and
challenges of building an underground detector and the operation of the iKAGRA
interferometer together with the geophysics interferometer that has been
constructed in the same tunnel.
[30]
oai:arXiv.org:1708.05252 [pdf] - 1608427
Correlated Signatures of Gravitational-Wave and Neutrino Emission in
Three-Dimensional General-Relativistic Core-Collapse Supernova Simulations
Submitted: 2017-08-17, last modified: 2017-11-04
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.
[31]
oai:arXiv.org:1710.04823 [pdf] - 2104587
The status of KAGRA underground cryogenic gravitational wave telescope
KAGRA Collaboration;
Akutsu, T.;
Ando, M.;
Araya, A.;
Aritomi, N.;
Asada, H.;
Aso, Y.;
Atsuta, S.;
Awai, K.;
Barton, M. A.;
Cannon, K.;
Craig, K.;
Creus, W.;
Doi, K.;
Eda, K.;
Enomoto, Y.;
Flaminio, R.;
Fujii, Y.;
Fujimoto, M. -K.;
Furuhata, T.;
Haino, S.;
Hasegawa, K.;
Hashino, K.;
Hayama, K.;
Hirobayashi, S.;
Hirose, E.;
Hsieh, B. H.;
Inoue, Y.;
Ioka, K.;
Itoh, Y.;
Kaji, T.;
Kajita, T.;
Kakizaki, M.;
Kamiizumi, M.;
Kambara, S.;
Kanda, N.;
Kanemura, S.;
Kaneyama, M.;
Kang, G.;
Kasuya, J.;
Kataoka, Y.;
Kawai, N.;
Kawamura, S.;
Kim, C.;
Kim, H.;
Kim, J.;
Kim, Y.;
Kimura, N.;
Kinugawa, T.;
Kirii, S.;
Kitaoka, Y.;
Kojima, Y.;
Kokeyama, K.;
Komori, K.;
Kotake, K.;
Kumar, R.;
Lee, H.;
Lee, H.;
Liu, Y.;
Luca, N.;
Majorana, E.;
Mano, S.;
Marchio, M.;
Matsui, T.;
Matsushima, F.;
Michimura, Y.;
Miyakawa, O.;
Miyamoto, T.;
Miyamoto, A.;
Miyo, K.;
Miyoki, S.;
Morii, W.;
Morisaki, S.;
Moriwaki, Y.;
Morozumi, T.;
Musha, M.;
Nagano, S.;
Nagano, K.;
Nakamura, K.;
Nakamura, T.;
Nakano, H.;
Nakano, M.;
Nakao, K.;
Narikawa, T.;
Quynh, L. Nguyen;
Ni, W. -T.;
Ochi, T.;
Oh, J.;
Oh, S.;
Ohashi, M.;
Ohishi, N.;
Ohkawa, M.;
Okutomi, K.;
Oohara, K.;
Alleano, F. E. Peña;
Pinto, I.;
Sago, N.;
Saijo, M.;
Saito, Y.;
Sakai, K.;
Sakai, Y.;
Sasaki, Y.;
Sasaki, M.;
Sato, S.;
Sato, T.;
Sekiguchi, Y.;
Seto, N.;
Shibata, M.;
Shimoda, T.;
Shinkai, H.;
Shoda, A.;
Somiya, K.;
Son, E.;
Suemasa, A.;
Suzuki, T.;
Suzuki, T.;
Tagoshi, H.;
Takahashi, H.;
Takahashi, R.;
Takamori, A.;
Takeda, H.;
Tanaka, H.;
Tanaka, K.;
Tanaka, T.;
Tatsumi, D.;
Tomaru, T.;
Tomura, T.;
Travasso, F.;
Tsubono, K.;
Tsuchida, S.;
Uchikata, N.;
Uchiyama, T.;
Uehara, T.;
Ueki, S.;
Ueno, K.;
Ushiba, T.;
van Putten, M. H. P. M.;
Vocca, H.;
Wada, S.;
Wakamatsu, T.;
Yamada, T.;
Yamamoto, S.;
Yamamoto, T.;
Yamamoto, K.;
Yamamoto, A.;
Yokoyama, J.;
Yokozawa, T.;
Yoon, T. H.;
Yuzurihara, H.;
Zeidler, S.;
Zhu, Z. -H.
Submitted: 2017-10-13
KAGRA is a 3-km interferometric gravitational wave telescope located in the
Kamioka mine in Japan. It is the first km-class gravitational wave telescope
constructed underground to reduce seismic noise, and the first km-class
telescope to use cryogenic cooling of test masses to reduce thermal noise. The
construction of the infrastructure to house the interferometer in the tunnel,
and the initial phase operation of the interferometer with a simple 3-km
Michelson configuration have been completed. The first cryogenic operation is
expected in 2018, and the observing runs with a full interferometer are
expected in 2020s. The basic interferometer configuration and the current
status of KAGRA are described.
[32]
oai:arXiv.org:1708.08513 [pdf] - 1625375
Estimating the core compactness of massive stars with Galactic supernova
neutrinos
Submitted: 2017-08-28, last modified: 2017-09-05
We suggest the future detection of neutrinos from a Galactic core-collapse
supernova can be used to infer the progenitor's inner mass density structure.
We present the results from 20 axisymmetric core-collapse supernova simulations
performed with progenitors spanning initial masses in the range 11--30Msun, and
focus on their connections to the progenitor compactness. The compactness is a
measure of the mass density profile of the progenitor core and recent
investigations have suggested its salient connections to the outcomes of core
collapse. Our simulations confirm a correlation between the neutrinos emitted
during the accretion phase and the progenitor's compactness, and that the ratio
of observed neutrino events during the first hundreds of milliseconds provides
a promising handle on the progenitor's inner structure. Neutrino flavor mixing
during the accretion phase remains a large source of uncertainty.
[33]
oai:arXiv.org:1708.08393 [pdf] - 1587568
Nucleosynthesis in 2D Core-Collapse Supernovae of 11.2 and 17.0
M$_{\odot}$ Progenitors: Implications for Mo and Ru Production
Submitted: 2017-08-28
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
$\nu$p-process. Furthermore, the production of Ru in the $\nu$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.
[34]
oai:arXiv.org:1708.03738 [pdf] - 1587042
Probing mass-radius relation of protoneutron stars from
gravitational-wave asteroseismology
Submitted: 2017-08-11
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.
[35]
oai:arXiv.org:1602.03028 [pdf] - 1544073
Multi-messenger signals of long-term core-collapse supernova simulations
: synergetic observation strategies
Submitted: 2016-02-09, last modified: 2017-03-08
The next Galactic supernova is expected to bring great opportunities for the
direct detection of gravitational waves (GW), full flavor neutrinos, and
multi-wavelength photons. To maximize the science return from such a rare
event, it is essential to have established classes of possible situations and
preparations for appropriate observations. To this end, we use a long-term
numerical simulation of the core-collapse supernova (CCSN) of a 17 solar-mass
red supergiant progenitor to self-consistently model the multi-messenger
signals expected in GW, neutrino, and electromagnetic messengers. This
supernova model takes into account the formation and evolution of a
protoneutron star, neutrino-matter interaction, and neutrino transport, all
within a two-dimensional shock hydrodynamics simulation. With this, we
separately discuss three situations: (i) a CCSN at the Galactic Center, (ii) an
extremely nearby CCSN within hundreds of parsecs, and (iii) a CCSN in nearby
galaxies within several Mpc. These distance regimes necessitate different
strategies for synergistic observations. In a Galactic CCSN, neutrinos provide
strategic timing and pointing information. We explore how these in turn deliver
an improvement in the sensitivity of GW analyses and help to guarantee
observations of early electromagnetic signals. To facilitate the detection of
multi-messenger signals of CCSNe in extremely nearby and extragalactic
distances, we compile a list of nearby red supergiant candidates and a list of
nearby galaxies with their expected CCSN rates. By exploring the sequential
multi-messenger signals of a nearby CCSN, we discuss preparations for
maximizing successful studies of such an unprecedented stirring event.
[36]
oai:arXiv.org:1605.09215 [pdf] - 1486864
A New Gravitational-Wave Signature from Standing Accretion Shock
Instabilities in Supernovae
Submitted: 2016-05-30, last modified: 2016-08-24
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.
[37]
oai:arXiv.org:1602.06759 [pdf] - 1436307
Three-dimensional simulations of rapidly rotating core-collapse
supernovae: finding a neutrino-powered explosion aided by non-axisymmetric
flows
Submitted: 2016-02-22, last modified: 2016-05-19
We report results from a series of three-dimensional (3D) rotational
core-collapse simulations for $11.2$ and 27 $M_{/odot}$ stars employing
neutrino transport scheme by the isotropic diffusion source approximation. By
changing the initial strength of rotation systematically, we find a
rotation-assisted explosion for the 27$M_{/odot}$ progenitor, which fails in
the absence of rotation. The unique feature was not captured in previous
two-dimensional (2D) self-consistent rotating models because the growing
non-axisymmetric instabilities play a key role. In the rapidly rotating case,
strong spiral flows generated by the so-called low $T/|W|$ instability enhance
the energy transport from the proto-neutron star (PNS) to the gain region,
which makes the shock expansion more energetic. The explosion occurs more
strongly in the direction perpendicular to the rotational axis, which is
different from previous 2D predictions.
[38]
oai:arXiv.org:1501.06330 [pdf] - 1370163
A New Multi-Energy Neutrino Radiation-Hydrodynamics Code in Full General
Relativity and Its Application to Gravitational Collapse of Massive Stars
Submitted: 2015-01-26, last modified: 2016-01-08
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.
[39]
oai:arXiv.org:1406.6414 [pdf] - 1347422
The criterion of supernova explosion revisited: the mass accretion
history
Submitted: 2014-06-24, last modified: 2015-11-19
By performing neutrino-radiation hydrodynamic simulations in spherical
symmetry (1D) and axial symmetry (2D) with different progenitor models by
Woosley & Heger (2007) from 12 $M_{\odot}$ to 100 $M_{\odot}$, we find that all
1D runs fail to produce an explosion and several 2D runs succeed. The
difference in the shock evolutions for different progenitors can be interpreted
by the difference in their mass accretion histories, which are in turn
determined by the density structures of progenitors. The mass accretion history
has two phases in the majority of the models: the earlier phase in which the
mass accretion rate is high and rapidly decreasing and the later phase with a
low and almost constant accretion rate. They are separated by the so-called
turning point, the origin of which is a change of the accreting layer. We argue
that shock revival will most likely occur around the turning point and hence
that its location in the $\dot M$-$L_\nu$ plane will be a good measure for the
possibility of shock revival: if the turning point lies above the critical
curve and the system stays there for a long time, shock revival will obtain. In
addition, we develop a phenomenological model to approximately evaluate the
trajectories in the $\dot M$-$L_\nu$ plane, which, after calibrating free
parameters by a small number of 1D simulations, reproduces the location of the
turning point reasonably well by using the initial density structure of
progenitor alone. We suggest the application of the phenomenological model to a
large collection of progenitors in order to infer without simulations which
ones are more likely to explode.
[40]
oai:arXiv.org:1501.00966 [pdf] - 1334075
Coherent Network Analysis of Gravitational Waves from Three-Dimensional
Core-Collapse Supernova Models
Submitted: 2015-01-02, last modified: 2015-11-15
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.
[41]
oai:arXiv.org:1406.2415 [pdf] - 1300121
Systematic Features of Axisymmetric Neutrino-Driven Core-Collapse
Supernova Models in Multiple Progenitors
Submitted: 2014-06-10, last modified: 2015-07-28
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 $\xi$ that characterizes the structure of the
progenitors is also a key in 2D to diagnose the properties of neutrino-driven
explosions. Models with high $\xi$ 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 $\xi$,
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
$\xi$.
[42]
oai:arXiv.org:1411.6705 [pdf] - 1223326
Magnetohydrodynamic Turbulence Powered by Magnetorotational Instability
in Nascent Proto-Neutron Stars
Submitted: 2014-11-24, last modified: 2014-12-11
Magnetorotational instability (MRI) in a convectively-stable layer around the
neutrinosphere is simulated by a three-dimensional model of supernova core. To
resolve MRI-unstable modes, a thin layer approximation considering only the
radial global stratification is adopted. Our intriguing finding is that the
convectively-stable layer around the neutrinosphere becomes fully-turbulent due
to the MRI and its nonlinear penetration into the strongly-stratified
MRI-stable region. The intensity of the MRI-driven turbulence increases with
magnetic flux threading the core, but is limited by a free energy stored in the
differential rotation. The turbulent neutrinosphere is a natural consequence of
rotating core-collapse and could exert a positive impact on the supernova
mechanism.
[43]
oai:arXiv.org:1409.0006 [pdf] - 873916
The red supergiant and supernova rate problems: implications for
core-collapse supernova physics
Submitted: 2014-08-29, last modified: 2014-09-30
Mapping supernovae to their progenitors is fundamental to understanding the
collapse of massive stars. We investigate the red supergiant problem, which
concerns why red supergiants with masses $\sim16$-$30 M_\odot$ have not been
identified as progenitors of Type IIP supernovae, and the supernova rate
problem, which concerns why the observed cosmic supernova rate is smaller than
the observed cosmic star formation rate. We find key physics to solving these
in the compactness parameter, which characterizes the density structure of the
progenitor. If massive stars with compactness above $\xi_{2.5} \sim 0.2$ fail
to produce canonical supernovae, (i) stars in the mass range $16$-$30 M_\odot$
populate an island of stars that have high $\xi_{2.5}$ and do not produce
canonical supernovae, and (ii) the fraction of such stars is consistent with
the missing fraction of supernovae relative to star formation. We support this
scenario with a series of two- and three-dimensional radiation hydrodynamics
core-collapse simulations. Using more than 300 progenitors covering initial
masses $10.8$-$75 M_\odot$ and three initial metallicities, we show that high
compactness is conducive to failed explosions. We then argue that a critical
compactness of $\sim 0.2$ as the divide between successful and failed
explosions is consistent with state-of-the-art three-dimensional core-collapse
simulations. Our study implies that numerical simulations of core collapse need
not produce robust explosions in a significant fraction of compact massive star
initial conditions.
[44]
oai:arXiv.org:1403.7290 [pdf] - 1208673
Impacts of Rotation on Three-dimensional Hydrodynamics of Core-collapse
Supernovae
Submitted: 2014-03-28, last modified: 2014-07-08
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.
[45]
oai:arXiv.org:1304.4372 [pdf] - 785175
Gravitational Wave Signatures from Low-mode Spiral Instabilities in
Rapidly Rotating Supernova Cores
Submitted: 2013-04-16, last modified: 2014-02-14
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.
[46]
oai:arXiv.org:1401.2723 [pdf] - 1202696
An Investigation into Surface Temperature Distributions of High-B
Pulsars
Submitted: 2014-01-13
Bearing in mind the application to high-magnetic-field (high-B) radio
pulsars, we investigate two-dimensional (2D) thermal evolutions of neutron
stars (NSs). We pay particular attention to the influence of different
equilibrium configurations on the surface temperature distributions. The
equilibrium configurations are constructed in a systematic manner, in which
both toroidal and poloidal magnetic fields are determined self-consistently
with the inclusion of general relativistic effects. To solve the 2D heat
transfer inside the NS interior out to the crust, we have developed an implicit
code based on a finite-difference scheme that deals with anisotropic thermal
conductivity and relevant cooling processes in the context of a standard
cooling scenario. In agreement with previous studies, the surface temperatures
near the pole become higher than those in the vicinity of the equator as a
result of anisotropic heat transfer. Our results show that the ratio of the
highest to the lowest surface temperatures changes maximally by one order of
magnitude, depending on the equilibrium configurations. Despite such
difference, we find that the area of such hot and cold spots is so small that
the simulated X-ray spectrum could be well reproduced by a single temperature
blackbody fitting.
[47]
oai:arXiv.org:1207.5955 [pdf] - 791652
Revisiting Impacts of Nuclear Burning for Reviving Weak Shocks in
Neutrino-Driven Supernovae
Submitted: 2012-07-25, last modified: 2014-01-08
We revisit potential impacts of nuclear burning on the onset of the
neutrino-driven explosions of core-collapse supernovae. By changing the
neutrino luminosity and its decay time to obtain parametric explosions in
one-(1D) and two-dimensional (2D) models with or without a 13-isotope alpha
network, we study how the inclusion of nuclear burning could affect the
postbounce dynamics for four progenitor models; three for 15.0 Msun stars, one
for an 11.2 Msun star. We find that the energy supply due to nuclear burning of
infalling material behind the shock can energize the shock expansion especially
for models that produce only marginal explosions in the absence of nuclear
burning. These models are energized by nuclear energy deposition when the shock
front passes through the silicon-rich layer and/or later it touches the
oxygen-rich layer. Depending on the neutrino luminosity and its decay time, a
diagnostic energy of explosion increases up to a few times 10^50 erg for models
with nuclear burning compared to the corresponding models without. We point out
that these features are most remarkable for the Limongi-Chieffi progenitor in
both 1D and 2D, because the progenitor model possesses a massive oxygen layer
with its inner-edge radius being smallest among the employed progenitors, so
that the shock can touch the rich fuel on a shorter timescale after bounce. The
energy difference is generally smaller (~0.1-0.2 times 10^51 erg) in 2D than in
1D (at most ~0.6 times 10^51 erg). This is because neutrino-driven convection
and the shock instability in 2D models enhance the neutrino heating efficiency,
which makes the contribution of nuclear burning relatively smaller compared to
1D models. Considering uncertainties in progenitor models, our results indicate
that nuclear burning should remain as one of the important ingredients to
foster the onset of neutrino-driven explosions.
[48]
oai:arXiv.org:1308.5755 [pdf] - 1178799
A Comparison of Two- and Three-dimensional Neutrino-hydrodynamics
simulations of Core-collapse Supernovae
Submitted: 2013-08-27, last modified: 2014-01-08
We present numerical results on two- (2D) and three-dimensional (3D)
hydrodynamic core-collapse simulations of an 11.2$M_\odot$ star. By changing
numerical resolutions and seed perturbations systematically, we study how the
postbounce dynamics is different in 2D and 3D. The calculations were performed
with an energy-dependent treatment of the neutrino transport based on the
isotropic diffusion source approximation scheme, which we have updated to
achieve a very high computational efficiency. All the computed models in this
work including nine 3D models and fifteen 2D models exhibit the revival of the
stalled bounce shock, leading to the possibility of explosion. All of them are
driven by the neutrino-heating mechanism, which is fostered by neutrino-driven
convection and the standing-accretion-shock instability (SASI). Reflecting the
stochastic nature of multi-dimensional (multi-D) neutrino-driven explosions,
the blast morphology changes from models to models. However, we find that the
final fate of the multi-D models whether an explosion is obtained or not, is
little affected by the explosion stochasticity. In agreement with some previous
studies, higher numerical resolutions lead to slower onset of the shock revival
in both 3D and 2D. Based on the self-consistent supernova models leading to the
possibility of explosions, our results systematically show that the revived
shock expands more energetically in 2D than in 3D.
[49]
oai:arXiv.org:1203.5192 [pdf] - 732488
Multimessenger astronomy with gravitational waves and high-energy
neutrinos
Ando, S.;
Baret, B.;
Bouhou, B.;
Chassande-Mottin, E.;
Kouchner, A.;
Moscoso, L.;
Van Elewyck, V.;
Bartos, I.;
Márka, S.;
Márka, Z.;
Corsi, A.;
Di Palma, I.;
Papa, M. A.;
Dietz, A.;
Donzaud, C.;
Eichler, D.;
Finley, C.;
Guetta, D.;
Halzen, F.;
Jones, G.;
Sutton, P. J.;
Kandhasamy, S.;
Mandic, V.;
Thrane, E.;
Kotake, K.;
Piran, T.;
Pradier, T.;
Romero, G. E.;
Waxman, E.
Submitted: 2012-03-23, last modified: 2013-07-25
Many of the astrophysical sources and violent phenomena observed in our
Universe are potential emitters of gravitational waves (GW) and high-energy
neutrinos (HEN). Both GWs and HENs may escape very dense media and travel
unaffected over cosmological distances, carrying information from the innermost
regions of the astrophysical engines. Such messengers could also reveal new,
hidden sources that have not been observed by conventional photon-based
astronomy. Coincident observation of GWs and HENs may thus play a critical role
in multimessenger astronomy. This is particularly true at the present time
owing to the advent of a new generation of dedicated detectors: IceCube,
ANTARES, VIRGO and LIGO. Given the complexity of the instruments, a successful
joint analysis of this data set will be possible only if the expertise and
knowledge of the data is shared between the two communities. This review aims
at providing an overview of both theoretical and experimental state-of-the-art
and perspectives for such a GW+HEN multimessenger astronomy.
[50]
oai:arXiv.org:1211.6817 [pdf] - 1158146
Effects of Resistivity on Magnetized Core-Collapse Supernovae
Submitted: 2012-11-29, last modified: 2013-01-23
We studied roles of a turbulent resistivity in the core-collapse of a
strongly magnetized massive star, carrying out 2D-resistive-MHD simulations.
The three cases with different initial strengths of magnetic field and rotation
are investigated; 1. strongly magnetized rotating core; 2.moderately magnetized
rotating core; 3. very strongly magnetized non-rotating core. In each case,
both an ideal-MHD model and resistive-MHD models are computed. As a result of
computations, each model shows a matter eruption helped by a magnetic
acceleration (and also by a centrifugal acceleration in the rotating cases). We
found that a resistivity attenuates the explosion in case~1 and 2, while it
enhances the explosion in case~3. We also found that in the rotating cases,
main mechanisms for the amplification of a magnetic field in the post-bounce
phase are an outward advection of magnetic field and a winding of poloidal
magnetic field-lines by differential rotation, which are somewhat dampened down
with the presence of a resistivity. Although the magnetorotational instability
seems to occur in the rotating models, it will play only a minor role in a
magnetic field amplification. Another impact of resistivity is that on the
aspect ratio. In the rotating cases, a large aspect ratio of the ejected
matters, $> 2.5$, attained in a ideal-MHD model is reduced to some extent in a
resistive model. These results indicate that a resistivity possibly plays an
important role in the dynamics of strongly magnetized supernovae.
[51]
oai:arXiv.org:1206.6101 [pdf] - 968391
On the Importance of the Equation of State for the Neutrino-Driven
Supernova Explosion Mechanism
Submitted: 2012-06-26, last modified: 2012-12-11
By implementing widely-used equations of state (EOS) from Lattimer & Swesty
(LS) and H. Shen et al. (SHEN) in core-collapse supernova simulations, we
explore possible impacts of these EOS on the post-bounce dynamics prior to the
onset of neutrino-driven explosions. Our spherically symmetric (1D) and axially
symmetric (2D) models are based on neutrino radiation hydrodynamics including
spectral transport, which is solved by the isotropic diffusion source
approximation. We confirm that in 1D simulations neutrino-driven explosions
cannot be obtained for any of the employed EOS. Impacts of the EOS on the
post-bounce hydrodynamics are more clearly visible in 2D simulations. In 2D
models of a 15 M_sun progenitor using the LS EOS, the stalled bounce shock
expands to increasingly larger radii, which is not the case using the SHEN EOS.
Keeping in mind that the omission of the energy drain by heavy-lepton neutrinos
in the present scheme could facilitate explosions, we find that 2D models of an
11.2 M_sun progenitor produce neutrino-driven explosions for all the EOS under
investigation. Models using the LS EOS are slightly more energetic compared to
those with the SHEN EOS. The more efficient neutrino heating in the LS models
coincides with a higher electron antineutrino luminosity and a larger mass that
is enclosed within the gain region. The models based on the LS EOS also show a
more vigorous and aspherical downflow of accreting matter to the surface of the
protoneutron star (PNS). The accretion pattern is essential for the production
and strength of outgoing pressure waves, that can push in turn the shock to
larger radii and provide more favorable conditions for the explosion.
[abbreviated]
[52]
oai:arXiv.org:1203.6488 [pdf] - 584414
Explosive Nucleosynthesis in Magnetohydrodynamical Jets from Collapsars
II. Heavy-Element Nucleosynthesis of s, r, p-Processes
Submitted: 2012-03-29, last modified: 2012-10-16
We investigate the nucleosynthesis in a massive star of 70 M_solar with solar
metallicity in the main sequence stage. The helium core mass after hydrogen
burning corresponds to 32 M_solar. Nucleosynthesis calculations have been
performed during the stellar evolution and the jetlike supernova explosion of a
collapsar model, where the weak s-, p-, and r-processes are taken into account.
We confirm that s-elements of 60 < A < 90 are highly overproduced relative to
the solar abundances in the hydrostatic nucleosynthesis. During oxygen burning,
p-elements of A > 90 are produced via photodisintegrations of seed s-elements.
However, the produced p-elements are disintegrated in later stages except for
^{180}Ta. In the explosive nucleosynthesis, elements of 90 < A < 160 are
significantly overproduced relative to the solar values owing to the r-process.
Only heavy p-elements (N > 50) are overproduced via the p-process. Compared
with the previous study of r-process nucleosynthesis calculations in the
collapsar model of 40 M_solar by Fujimoto et al. 2007, 2008, our jet model
cannot contribute to the third peak of the solar r-elements and intermediate
p-elements. Averaging the overproduction factors over the progenitor masses
with the use of Salpeter's IMF, we suggest that the 70 M_solar star could
contribute to the solar weak s-elements of 60 < A < 90 and neutron-rich
elements of 90 < A < 160. We confirm the primary synthesis of light p-elements
in the ejected matter of high peak temperature. The ejected matter has [Sr/Eu]
\sim -0.4, which is different from that of a typical r-process-enriched star
CS22892-052 ([Sr/Eu] \sim -1). We find that Sr-Y-Zr isotopes are primarily
synthesized in the explosive nucleosynthesis in a similar process of the
primary production of light p-elements, which has been considered as one of the
sites of a lighter element primary process (LEPP).
[53]
oai:arXiv.org:1209.2360 [pdf] - 1151317
Local Simulations of the Magneto-rotational Instability in Core-Collapse
Supernovae
Submitted: 2012-09-11
Bearing in mind the application to core-collapse supernovae, we study
nonlinear properties of the magneto-rotational instability (MRI) by means of
three- dimensional simulations in the framework of a local shearing box
approximation. By changing systematically the shear rates that symbolize the
degree of differential rotation in nascent proto-neutron stars (PNSs), we
derive a scaling relation between the turbulent stress sustained by the MRI and
the shear- vorticity ratio. Our parametric survey shows a power-law scaling
between the turbulent stress ($<< w_{\rm tot}>>$) and the shear- vorticity
ratio ($g_q$) as $<<w_{\rm tot}>> \propto g_q^{\delta}$ with its index $\delta
\sim 0.5$. The MRI-amplified magnetic energy has a similar scaling relative to
the turbulent stress, while the Maxwell stress has slightly smaller power-law
index ($\sim 0.36$). By modeling the effect of viscous heating rates due to the
MRI turbulence, we show that the stronger magnetic fields or the larger shear
rates initially imposed lead to the higher dissipation rates. For a rapidly
rotating PNS with the spin period in milliseconds and with strong magnetic
fields of $10^{15}$ G, the energy dissipation rate is estimated to exceed
$10^{51} {\rm erg\ sec^{-1}}$. Our results suggest that the conventional
magnetohydrodynamic (MHD) mechanism of core-collapse supernovae is likely to be
affected by the MRI-driven turbulence, which we speculate, on one hand, could
harm the MHD-driven explosions due to the dissipation of the shear rotational
energy at the PNS surface, on the other hand the energy deposition there might
be potentially favorable for the working of the neutrino-heating mechanism.
[54]
oai:arXiv.org:1205.6284 [pdf] - 656859
Core-Collapse Supernovae as Supercomputing Science: a status report
toward 6D simulations with exact Boltzmann neutrino transport in full general
relativity
Submitted: 2012-05-29, last modified: 2012-06-26
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.
[55]
oai:arXiv.org:1202.2487 [pdf] - 1116573
Fully General Relativistic Simulations of Core-Collapse Supernovae with
An Approximate Neutrino Transport
Submitted: 2012-02-11, last modified: 2012-05-29
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.
[56]
oai:arXiv.org:1205.6061 [pdf] - 1123705
Gravitational Wave Signatures of Hyperaccreting Collapsar Disks
Submitted: 2012-05-28
By performing two-dimensional special relativistic (SR) magnetohydrodynamic
simulations, we study possible signatures of gravitational waves (GWs) in the
context of the collapsar model for long-duration gamma-ray bursts. In our SR
simulations, the central black hole is treated as an absorbing boundary. By
doing so, we focus on the GWs generated by asphericities in neutrino emission
and matter motions in the vicinity of the hyperaccreting disks. We compute nine
models by adding initial angular momenta and magnetic fields parametrically to
a precollapse core of a $35 M_{\odot}$ progenitor star. As for the
microphysics, a realistic equation of state is employed and the neutrino
cooling is taken into account via a multiflavor neutrino leakage scheme. To
accurately estimate GWs produced by anisotropic neutrino emission, we perform a
ray-tracing analysis in general relativity by a post-processing procedure. By
employing a stress formula that includes contributions both from magnetic
fields and special relativistic corrections, we study also the effects of
magnetic fields on the gravitational waveforms. We find that the GW amplitudes
from anisotropic neutrino emission show a monotonic increase with time, whose
amplitudes are much larger than those from matter motions of the accreting
material. We show that the increasing trend of the neutrino GWs stems from the
excess of neutrino emission in the direction near parallel to the spin axis
illuminated from the hyperaccreting disks. We point out that a recently
proposed future space-based interferometer like Fabry-Perot type DECIGO would
permit the detection of these GW signals within $\approx$ 100 Mpc.
[57]
oai:arXiv.org:1204.2330 [pdf] - 1117889
Multimessengers from core-collapse supernovae: multidimensionality as a
key to bridge theory and observation
Submitted: 2012-04-10, last modified: 2012-04-24
Core-collapse supernovae are dramatic explosions marking the catastrophic end
of massive stars. The only means to get direct information about the supernova
engine is from observations of neutrinos emitted by the forming neutron star,
and through gravitational waves which are produced when the hydrodynamic flow
or the neutrino flux is not perfectly spherically symmetric. The
multidimensionality of the supernova engine, which breaks the sphericity of the
central core such as convection, rotation, magnetic fields, and hydrodynamic
instabilities of the supernova shock, is attracting great attention as the most
important ingredient to understand the long-veiled explosion mechanism. Based
on our recent work, we summarize properties of gravitational waves, neutrinos,
and explosive nucleosynthesis obtained in a series of our multidimensional
hydrodynamic simulations and discuss how the mystery of the central engines can
be unraveled by deciphering these multimessengers produced under the thick
veils of massive stars.
[58]
oai:arXiv.org:1110.5107 [pdf] - 1085083
Multiple physical elements to determine the gravitational-wave
signatures of core-collapse supernovae
Submitted: 2011-10-23, last modified: 2012-04-10
We review recent progress in the theoretical predictions of gravitational
waves (GWs) of core-collapse supernovae. Following a brief summary of the
methods in the numerical modeling, we summarize multiple physical elements that
determine the GW signatures which have been considered to be important in
extracting the information of the long-veiled explosion mechanism from the
observation of the GWs. We conclude with a summary of the most urgent tasks to
make the dream come true.
[59]
oai:arXiv.org:1204.1655 [pdf] - 497572
R-process Nucleosynthesis during the Magnetohydrodynamics Explosions of
a Massive Star
Submitted: 2012-04-07
We investigate the possibility of the r-process during the
magnetohydrohynamical explosion of supernova in a massive star of 13 solar mass
with the effects of neutrinos induced. We adopt five kinds of initial models
which include properties of rotation and the toroidal component of the magnetic
field . The simulations which succeed the explosions are limitted to a
concentrated magnetic field and strong differential rotation. Low $Y_{e}$
ejecta produce heavy elements and the third peak can be reprocuced. However,
the second peak is low because $Y_{e}$ distribution as a function of radius is
steep and ejecta corresponding to middle $Y_{e}$ is very few.
[60]
oai:arXiv.org:1108.3989 [pdf] - 1083495
Three-dimensional Hydrodynamic Core-Collapse Supernova Simulations for
an $11.2 M_{\odot}$ Star with Spectral Neutrino Transport
Submitted: 2011-08-19, last modified: 2012-01-25
We present numerical results on three-dimensional (3D) hydrodynamic
core-collapse simulations of an $11.2 M_{\odot}$ star. By comparing one-(1D)
and two-dimensional(2D) results with those of 3D, we study how the increasing
spacial multi-dimensionality affects the postbounce supernova dynamics. The
calculations were performed with an energy-dependent treatment of the neutrino
transport that is solved by the isotropic diffusion source approximation
scheme. By performing a tracer-particle analysis, we show that the maximum
residency time of material in the gain region is shown to be longer for 3D due
to non-axisymmetric flow motions than 2D, which is one of advantageous aspects
of 3D models to obtain neutrino-driven explosions. Our results show that
convective matter motions below the gain radius become much more violent in 3D
than 2D, making the neutrino luminosity larger for 3D. Nevertheless the emitted
neutrino energies are made smaller due to the enhanced cooling. Our results
indicate whether these advantages for driving 3D explosions could or could not
overwhelm the disadvantages is sensitive to the employed numerical resolutions.
An encouraging finding is that the shock expansion tends to become more
energetic for models with finer resolutions. To draw a robust conclusion, 3D
simulations with much more higher numerical resolutions and also with more
advanced treatment of neutrino transport as well as of gravity is needed, which
could be hopefully practicable by utilizing forthcoming Petaflops-class
supercomputers.
[61]
oai:arXiv.org:1004.2896 [pdf] - 1026368
Gravitational Wave Signatures of Magnetohydrodynamically-Driven
Core-Collapse Supernova Explosions
Submitted: 2010-04-16, last modified: 2011-08-19
By performing a series of two-dimensional, special relativistic
magnetohydrodynamic (MHD) simulations, we study signatures of gravitational
waves (GWs) in the magnetohydrodynamically-driven core-collapse supernovae. In
order to extract the gravitational waveforms, we present a stress formula
including contributions both from magnetic fields and special relativistic
corrections. By changing the precollapse magnetic fields and initial angular
momentum distributions parametrically, we compute twelve models. As for the
microphysics, a realistic equation of state is employed and the neutrino
cooling is taken into account via a multiflavor neutrino leakage scheme.
With these computations, we find that the total GW amplitudes show a
monotonic increase after bounce for models with a strong precollapse magnetic
field ($10^{12}$G) also with a rapid rotation imposed.
We show that this trend stems both from the kinetic contribution of MHD
outflows with large radial velocities and also from the magnetic contribution
dominated by the toroidal magnetic fields that predominantly trigger MHD
explosions. For models with weaker initial magnetic fields, the total GW
amplitudes after bounce stay almost zero, because the contribution from the
magnetic fields cancels with the one from the hydrodynamic counterpart.
These features can be clearly understood with a careful analysis on the
explosion dynamics. We point out that the GW signals with the increasing trend,
possibly visible to the next-generation detectors for a Galactic supernova,
would be associated with MHD explosions with the explosion energies exceeding
$10^{51}$ erg.
[62]
oai:arXiv.org:1106.5487 [pdf] - 1077577
Impacts of Collective Neutrino Oscillations on Supernova Explosions
Submitted: 2011-06-27
By performing a series of one- and two-dimensional (1-, 2D) hydrodynamic
simulations with spectral neutrino transport, we study possible impacts of
collective neutrino oscillations on the dynamics of core-collapse supernovae.
To model the spectral swapping which is one of the possible outcome of the
collective neutrino oscillations, we parametrize the onset time when the
spectral swap begins, the radius where the spectral swap occurs, and the
threshold energy above which the spectral interchange between heavy-lepton
neutrinos and electron/anti-electron neutrinos takes place, respectively. By
doing so, we systematically study how the neutrino heating enhanced by the
spectral swapping could affect the shock evolution as well as the matter
ejection. We also investigate the progenitor dependence using a suite of
progenitor models (13, 15, 20, and 25 $M_\odot$). We find that there is a
critical heating rate induced by the spectral swapping to trigger explosions,
which significantly differs between the progenitors. The critical heating rate
is generally smaller for 2D than 1D due to the multidimensionality that
enhances the neutrino heating efficiency. For the progenitors employed in this
paper, the final remnant masses are estimated to range in 1.1-1.5$M_\odot$. For
our 2D model of the $15M_\odot$ progenitor, we find a set of the oscillation
parameters that could account for strong supernova explosions ($\sim 10^{51}$
erg), simultaneously leaving behind the remnant mass close to $\sim 1.4
M_\odot$.
[63]
oai:arXiv.org:1106.2606 [pdf] - 1077256
Explosive nucleosynthesis in the neutrino-driven aspherical supernova
explosion of a non-rotating 15$M_{\odot}$ star with solar metallicity
Submitted: 2011-06-14
We investigate explosive nucleosynthesis in a non-rotating 15$M_\odot$ star
with solar metallicity that explodes by a neutrino-heating supernova (SN)
mechanism aided by both standing accretion shock instability (SASI) and
convection. To trigger explosions in our two-dimensional hydrodynamic
simulations, we approximate the neutrino transport with a simple light-bulb
scheme and systematically change the neutrino fluxes emitted from the
protoneutron star. By a post-processing calculation, we evaluate abundances and
masses of the SN ejecta for nuclei with the mass number $\le 70$ employing a
large nuclear reaction network. Aspherical abundance distributions, which are
observed in nearby core-collapse SN remnants, are obtained for the non-rotating
spherically-symmetric progenitor, due to the growth of low-mode SASI. Abundance
pattern of the supernova ejecta is similar to that of the solar system for
models whose masses ranges $(0.4-0.5) \Ms$ of the ejecta from the inner region
($\le 10,000\km$) of the precollapse core. For the models, the explosion
energies and the \nuc{Ni}{56} masses are $ \simeq 10^{51} \rm erg$ and
$(0.05-0.06) \Ms$, respectively; their estimated baryonic masses of the neutron
star are comparable to the ones observed in neutron-star binaries. These
findings may have little uncertainty because most of the ejecta is composed by
matter that is heated via the shock wave and has relatively definite
abundances. The abundance ratios for Ne, Mg, Si and Fe observed in Cygnus loop
are well reproduced with the SN ejecta from an inner region of the $15\Ms$
progenitor.
[64]
oai:arXiv.org:1106.0544 [pdf] - 1077049
Effects of Rotation on Stochasticity of Gravitational Waves in Nonlinear
Phase of Core-Collapse Supernovae
Submitted: 2011-06-02
By performing three-dimensional (3D) simulations that demonstrate the
neutrino-driven core-collapse supernovae aided by the standing accretion shock
instability (SASI), we study how the spiral modes of the SASI can have impacts
on the properties of the gravitational-wave (GW) emission. To see the effects
of rotation in the non-linear postbounce phase, we give a uniform rotation on
the flow advecting from the outer boundary of the iron core, whose specific
angular momentum is assumed to agree with recent stellar evolution models. We
compute fifteen 3D models in which the initial angular momentum as well as the
input neutrino luminosities from the protoneutron star are changed in a
systematic manner. By performing a ray-tracing analysis, we accurately estimate
the GW amplitudes generated by anisotropic neutrino emission. Our results show
that the gravitational waveforms from neutrinos in models that include rotation
exhibit a common feature otherwise they vary much more stochastically in the
absence of rotation. The breaking of the stochasticity stems from the excess of
the neutrino emission parallel to the spin axis. This is because the
compression of matter is more enhanced in the vicinity of the equatorial plane
due to the growth of the spiral SASI modes, leading to the formation of spiral
flows circulating around the spin axis with higher temperatures. We point out
that a recently proposed future space interferometers like Fabry-Perot type
DECIGO would permit detection of these signals for a Galactic supernova.
[65]
oai:arXiv.org:1012.4080 [pdf] - 1042655
Magnetic Field Decay Due to the Wave-Particle Resonances in the Outer
Crust of the Neutron Star
Submitted: 2010-12-18
Bearing in mind the application to the outer crust of the neutron stars
(NSs), we investigate the magnetic field decay by means of the fully
relativistic Particle-In-Cell simulations. Numerical computations are carried
out in 2-dimensions, in which the initial magnetic fields are set to be
composed both of the uniform magnetic fields that model the global fields
penetrating the NS and of the turbulent magnetic fields that would be
originated from the Hall cascade of the large-scale turbulence. Our results
show that the whistler cascade of the turbulence transports the magnetic energy
preferentially in the direction perpendicular to the uniform magnetic fields.
It is also found that the distribution function of electrons becomes
anisotropic because electrons with lower energies are predominantly heated in
the direction parallel to the uniform magnetic fields due to the Landau
resonance, while electrons with higher energies are heated mainly by the
cyclotron resonance that makes the distribution function isotropic for the high
energy tails. Furthermore we point out that the degree of anisotropy takes
maximum as a function of the initial turbulent magnetic energy. As an
alternative to the conventional ohmic dissipation, we propose that the magnetic
fields in the outer crust of NSs, cascading down to the electron inertial scale
via the whistler turbulence, would decay predominantly by the dissipation
processes through the Landau damping and the cyclotron resonance.
[66]
oai:arXiv.org:0912.1157 [pdf] - 1018643
Explosion geometry of a rotating 13 $M_{\odot}$ star driven by the
SASI-aided neutrino-heating supernova mechanism
Submitted: 2009-12-06, last modified: 2010-09-28
By performing axisymmetric hydrodynamic simulations of core-collapse
supernovae with spectral neutrino transport based on the isotropic diffusion
source approximation scheme, we support the assumption that the
neutrino-heating mechanism aided by the standing accretion shock instability
and convection can initiate an explosion of a 13 $M_{\odot}$ star. Our results
show that bipolar explosions are more likely to be associated with models which
include rotation. We point out that models, which form a north-south symmetric
bipolar explosion, can lead to larger explosion energies than for the
corresponding unipolar explosions.
[67]
oai:arXiv.org:1007.3165 [pdf] - 1033757
General Relativistic Ray-Tracing Method for Estimating the Energy and
Momentum Deposition by Neutrino Pair Annihilation in Collapsars
Submitted: 2010-07-19, last modified: 2010-07-23
Bearing in mind the application to the collapsar models of gamma-ray bursts
(GRBs), we develop a numerical scheme and code for estimating the deposition of
energy and momentum due to the neutrino pair annihilation ($\nu + {\bar \nu}
\rightarrow e^{-} + e^{+}$) in the vicinity of accretion tori around a Kerr
black hole. Our code is designed to solve the general relativistic neutrino
transfer by a ray-tracing method. To solve the collisional Boltzmann equation
in curved spacetime, we numerically integrate the so-called rendering equation
along the null geodesics. For the neutrino opacity, the charged-current
$\beta$-processes are taken into account, which are dominant in the vicinity of
the accretion tori. The numerical accuracy of the developed code is
certificated by several tests, in which we show comparisons with the
corresponding analytic solutions. Based on the hydrodynamical data in our
collapsar simulation, we estimate the annihilation rates in a post-processing
manner. Increasing the Kerr parameter from 0 to 1, it is found that the general
relativistic effect can increase the local energy deposition rate by about one
order of magnitude, and the net energy deposition rate by several tens of
percents. After the accretion disk settles into a stationary state (typically
later than $\sim 9$ s from the onset of gravitational collapse), we point out
that the neutrino-heating timescale in the vicinity of the polar funnel region
can be shorter than the dynamical timescale. Our results suggest the neutrino
pair annihilation has a potential importance equal to the conventional
magnetohydrodynamic mechanism for igniting the GRB fireballs.
[68]
oai:arXiv.org:0912.2590 [pdf] - 902735
Neutrino Pair Annihilation in Collapsars; Ray-Tracing Method in Special
Relativity
Submitted: 2009-12-14, last modified: 2010-02-26
We develop a numerical scheme and code for estimating the energy and momentum
transfer via neutrino pair annihilation ($\nu + {\bar \nu} \to e^{-}+ e^{+}$),
bearing in mind the application to the collapsar models of gamma-ray bursts
(GRBs). To calculate the neutrino flux illuminated from the accretion disk, we
perform a ray-tracing calculation in the framework of special relativity. The
numerical accuracy of the developed code is certificated by several tests, in
which we show comparisons with the corresponding analytical solutions. Using
hydrodynamical data in our collapsar simulation, we estimate the annihilation
rates in a post-processing manner. We show that the neutrino energy deposition
and momentum transfers are strongest near the inner edge of the accretion disk.
The beaming effects of special relativity are found to change the annihilation
rates by several factors in the polar funnel region. After the accretion disk
settles into a stationary state (typically later than
$\sim 9$ s from the onset of gravitational collapse), we find that the
neutrino-heating timescale in the vicinity of the polar funnel ($\lesssim 80$
km) can become shorter than the hydrodynamical timescale, indicating that the
neutrino-heated outflows can be launched there. We point out that the momentum
transfer can play as important role as the energy deposition for the efficient
acceleration of neutrino-driven outflows.
Our results suggest that the neutrino pair annihilation has a potential
importance equal to the conventional magnetohydrodynamic mechanism for igniting
the GRB fireballs.
[69]
oai:arXiv.org:0912.3577 [pdf] - 31943
Long-Term Evolution of Collapsar: Mechanism of Outflow Production
Submitted: 2009-12-17, last modified: 2009-12-20
We present our numerical results of two-dimensional hydrodynamic (HD)
simulations and magnetohydrodynamic (MHD) simulations of the collapse of
rotating massive stars in light of the collapsar model of gamma-ray bursts
(GRBs). Pushed by recent evolution calculations of GRB progenitors, we focus on
lower angular momentum of the central core than the ones taken mostly in
previous studies. By performing special relativistic simulations including both
realistic equation of state and neutrino cooling, we follow a long-term
evolution of the slowly rotating collapsars up to 10 s, accompanied by the
formation of jets and accretion disks. We find such outfows can be launched
both by MHD process and neutrino process. We investigeate the properties of
these jets whether it can become GRBs or remains primary weak outflow.
[70]
oai:arXiv.org:0911.3968 [pdf] - 30878
The pasta structure in the hadron-quark phase transition and the effects
on magnetized compact stars
Submitted: 2009-11-20, last modified: 2009-11-23
We report the effects of quark-hadron phase transition on the structures of
general relativistic stars with purely toroidal magnetic field. For the mixed
phase, we take into account of the finite-size effects, which lead to
non-uniform "Pasta" structures. Our study is based on axisymetric and
stationary formalism including purely toroidal magnetic field. For hybrid
stars, we find the characteristic distribution of magnetic field, which has a
discontinuity originated in the quark-hadron mixed phase. These distributions
of magnetic field will change astrophysical phenomena, such as cooling
processes.
[71]
oai:arXiv.org:0910.0327 [pdf] - 1017805
Relativistic hybrid stars with super-strong toroidal magnetic fields: An
evolutionary track with QCD phase transition
Submitted: 2009-10-02
We investigate structures of hybrid stars, which feature quark core
surrounded by a hadronic matter mantle, with super-strong toroidal magnetic
fields in full general relativity. Modeling the equation of state (EOS) with a
first order transition by bridging the MIT bag model for the description of
quark matter and the nuclear EOS by Shen et al., we numerically construct
thousands of the equilibrium configurations to study the effects of the phase
transition. It is found that the appearance of the quark phase can affect
distributions of the magnetic fields inside the hybrid stars, making the
maximum field strength about up to 30 % larger than for the normal neutron
stars. Using the equilibrium configurations, we explore the possible
evolutionary paths to the formation of hybrid stars due to the spin-down of
magnetized rotating neutron stars. We find that the energy release by the phase
transition to the hybrid stars is quite large ($\la 10^{52} \rm erg$) even for
super strongly magnetized compact stars. Our results suggest that the strong
gravitational-wave emission and the sudden spin-up signature could be
observable signals of the QCD phase transition, possibly for a source out to
Megaparsec distances.
[72]
oai:arXiv.org:0909.3622 [pdf] - 901786
Ray-Tracing Analysis of Anisotropic Neutrino Radiation for Estimating
Gravitational Waves in Core-Collapse Supernovae
Submitted: 2009-09-20
We propose a ray-tracing method to estimate gravitational waves (GWs)
generated by anisotropic neutrino emission in supernova cores. To calculate the
gravitational waveforms, we derive analytic formulae in a useful form, which
are applicable also for three-dimensional computations. Pushed by evidence of
slow rotation prior to core-collapse, we focus on asphericities in neutrino
emission and matter motions outside the protoneutron star. Based on the
two-dimensional (2D) models, which mimic SASI-aided neutrino heating
explosions, we compute the neutrino anisotropies via the ray-tracing method in
a post-processing manner and calculate the resulting waveforms. With these
computations, it is found that the waveforms exhibit more variety in contrast
to the ones previously estimated by the ray-by-ray analysis (e.g., Kotake et
al. (2007)). In addition to a positively growing feature, which was predicted
to determine the total wave amplitudes predominantly, the waveforms are shown
to exhibit large negative growth for some epochs during the growth of SASI.
These features are found to stem from the excess of neutrino emission in
lateral directions, which can be precisely captured by the ray-tracing
calculation. Due to the negative contributions and the neutrino absorptions
appropriately taken into account by the ray-tracing method, the wave amplitudes
become more than one-order-of magnitude smaller than the previous estimation,
thus making their detections very hard for a galactic source.On the other hand,
it is pointed out that the GW spectrum from matter motions have its peak near
$\sim 100$ Hz, which could be characteristic for the SASI-induced supernova
explosions.(abridged)
[73]
oai:arXiv.org:0906.3180 [pdf] - 25363
Neutrino oscillations in magnetically driven supernova explosions
Submitted: 2009-06-17, last modified: 2009-09-19
We investigate neutrino oscillations from core-collapse supernovae that
produce magnetohydrodynamic (MHD) explosions. By calculating numerically the
flavor conversion of neutrinos in the highly non-spherical envelope, we study
how the explosion anisotropy has impacts on the emergent neutrino spectra
through the Mikheyev-Smirnov-Wolfenstein effect. In the case of the inverted
mass hierarchy with a relatively large theta_(13), we show that survival
probabilities of electron type neutrinos and antineutrinos seen from the
rotational axis of the MHD supernovae (i.e., polar direction), can be
significantly different from those along the equatorial direction. The event
numbers of electron type antineutrinos observed from the polar direction are
predicted to show steepest decrease, reflecting the passage of the
magneto-driven shock to the so-called high-resonance regions. Furthermore we
point out that such a shock effect, depending on the original neutrino spectra,
appears also for the low-resonance regions, which leads to a noticeable
decrease in the electron type neutrino signals. This reflects a unique nature
of the magnetic explosion featuring a very early shock-arrival to the resonance
regions, which is in sharp contrast to the neutrino-driven delayed supernova
models. Our results suggest that the two features in the electron type
antineutrinos and neutrinos signals, if visible to the Super-Kamiokande for a
Galactic supernova, could mark an observational signature of the magnetically
driven explosions, presumably linked to the formation of magnetars and/or
long-duration gamma-ray bursts.
[74]
oai:arXiv.org:0905.2006 [pdf] - 24220
Long-Term Evolution of Slowly Rotating Collapsar in Special Relativistic
Magnetohydrodynamics
Submitted: 2009-05-12, last modified: 2009-08-27
We present our numerical results of two-dimensional magnetohydrodynamic (MHD)
simulations of the collapse of rotating massive stars in light of the collapsar
model of gamma-ray bursts (GRBs). Pushed by recent evolution calculations of
GRB progenitors, we focus on lower angular momentum of the central core than
the ones taken mostly in previous studies. By performing special relativistic
simulations including both realistic equation of state and neutrino coolings,
we follow a unprecedentedly long-term evolution of the slowly rotating
collapsars up to $\sim$ 10 s, accompanied by the formation of jets and
accretion disks. Our results show that for the GRB progenitors to function as
collapsars, there is a critical initial angular momentum, below which matter is
quickly swallowed to the central objects, no accretion disks and no MHD
outflows are formed. When larger than the criteria, we find the launch of the
MHD jets in the following two ways. For models with stronger initial magnetic
fields, the magnetic pressure amplified inside the accretion disk can drive the
MHD outflows, which makes the strong magnetic explosions like a 'magnetic
tower' (type II). For the models with weaker initial magnetic fields, the
magnetic tower stalls first and the subsequent MHD outflows are produced by
some eventual inflows of the accreting material from the equator to the polar
regions (type I). Regardless of type I or II, the jets can attain only mildly
relativistic speeds with the explosion energy less than $10^{49} \erg$, which
could possibly be related to the X-ray flashes. Due to high opacity for
neutrinos inside the disk, we find that the luminosities of $\nu_e$ and
$\bar{\nu}_e$ become almost comparable, which is advantageous for making the
energy deposition rate larger.
[75]
oai:arXiv.org:0906.4957 [pdf] - 25712
Joint searches between gravitational-wave interferometers and
high-energy neutrino telescopes: science reach and analysis strategies
Van Elewyck, V.;
Ando, S.;
Aso, Y.;
Baret, B.;
Barsuglia, M.;
Bartos, I.;
Chassande-Mottin, E.;
Di Palma, I.;
Dwyer, J.;
Finley, C.;
Kotake, K.;
Kouchner, A.;
Marka, S.;
Marka, Z.;
Rollins, J.;
Ott, C. D.;
Pradier, T.;
Searle, A.
Submitted: 2009-06-26
Many of the astrophysical sources and violent phenomena observed in our
Universe are potential emitters of gravitational waves (GWs) and high-energy
neutrinos (HENs). A network of GW detectors such as LIGO and Virgo can
determine the direction/time of GW bursts while the IceCube and ANTARES
neutrino telescopes can also provide accurate directional information for HEN
events. Requiring the consistency between both, totally independent, detection
channels shall enable new searches for cosmic events arriving from potential
common sources, of which many extra-galactic objects.
[76]
oai:arXiv.org:0904.4300 [pdf] - 23725
Stochastic Nature of Gravitational Waves from Supernova Explosions with
Standing Accretion Shock Instability
Submitted: 2009-04-28
We study properties of gravitational waves based on the three-dimensional
simulations, which demonstrate the neutrino-driven explosions aided by the
standing accretion shock instability (SASI). Pushed by evidence supporting slow
rotation prior to core-collapse, we focus on the asphericities in neutrino
emissions and matter motions outside the protoneutron star. By performing a
ray-tracing calculation in 3D, we estimate accurately the gravitational
waveforms from anisotropic neutrino emissions. In contrast to the previous work
assuming axisymmetry, we find that the gravitational waveforms vary much more
stochastically because the explosion anisotropies depend sensitively on the
growth of the SASI which develops chaotically in all directions. Our results
show that the gravitational-wave spectrum has its peak near $\sim 100$ Hz,
reflecting the SASI-induced matter overturns of $\sim O(10)$ ms. We point out
that the detection of such signals, possibly visible to the LIGO-class
detectors for a Galactic supernova, could be an important probe into the
long-veiled explosion mechanism.
[77]
oai:arXiv.org:0904.3203 [pdf] - 900708
General relativistic compact stars with exotic matter
Submitted: 2009-04-21
We study structures of general relativistic compact stars with exotic matter.
Our study is based on axisymetric and stationary formalism including purely
toroidal magnetic field. We also study the finite size effects of quark-hadron
mixed phase on structures of magnetars. For hybrid stars, we find a
characteristic distribution of magnetic field, which has a discontinuity
originated in the quark-hadron mixed phase. These distributions of magnetic
field will change astrophysical phenomena, such as cooling processes.
[78]
oai:arXiv.org:0904.2044 [pdf] - 550873
Relativistic stars with purely toroidal magnetic fields with realistic
equations of state
Submitted: 2009-04-14
We investigate equilibrium sequences of relativistic stars containing purely
toroidal magnetic fields with four kinds of realistic equations of state (EOSs)
of SLy (Douchin et al.), FPS (Pandharipande et al.), Shen (Shen et al.), and LS
(Lattimer & Swesty). We numerically construct thousands of equilibrium
configurations. Particularly we pay attention to the equilibrium sequences of
constant baryon mass and/or constant magnetic flux, which model evolutions of
an isolated neutron star. Important properties obtained in this study are
summarized as follows ; (1) The dependence of the mass-shedding angular
velocity on the EOSs is determined from that of the non-magnetized case. The
stars with Shen(FPS) EOS reach the mass-shedding limit at the smallest(largest)
angular velocity, while the stars with SLy or Lattimer-Swesty EOSs take the
moderate values. (2) For the supramassive sequences, the equilibrium
configurations are found to be generally oblate for the realistic EOSs in sharp
contrast to the polytropic stars. For FPS(LS) EOS, the parameter region which
permits the prolately deformed stars is widest(narrowest). For SLy and Shen
EOS, it is in medium. Furthermore, the angular velocities $\Omega_{\rm up}$,
above which the stars start to spin up as they lose angular momentum, are found
to depend sharply on the realistic EOSs. Our analysis indicates that the
hierarchy of this spin up angular velocity is $\Omega_{\rm up,SLy} >
\Omega_{\rm up,FPS} > \Omega_{\rm up,LS}>\Omega_{\rm up,Shen}$ and this
relation holds even if the sequences have strong magnetic fields. Our results
suggest the EOSs within the relativistic stars containing purely toroidal
magnetic fields can be constrained by observing the angular velocity, the
gravitational wave, and the signature of the spin up.
[79]
oai:arXiv.org:0712.1949 [pdf] - 8057
Special Relativistic Simulations of Magnetically-dominated Jets in
Collapsing Massive Stars
Submitted: 2007-12-12, last modified: 2008-12-10
We perform a series of two-dimensional magnetohydrodynamic core-collapse
simulations of rapidly rotating and strongly magnetized massive stars. To study
the properties of magnetic explosions for a longer time stretch of postbounce
evolution, we develop a new code under the framework of special relativity
including a realistic equation of state with a multiflavor neutrino leakage
scheme. Our results show the generation of the magnetically-dominated jets in
the two ways. One is launched just after the core-bounce in a prompt way and
another is launched at $ \sim 100 $ ms after the stall of the prompt shock. We
find that the shock-revival occurs when the magnetic pressure becomes strong,
due to the field wrapping, enough to overwhelm the ram pressure of the
accreting matter. The critical toroidal magnetic fields for the magnetic
shock-revival are found to be universal of $\sim 10^{15}\mathrm{G}$ behind the
jets. We point out that the time difference before the shock-revival has a
strong correlation with the explosions energies. Our results suggest that the
magnetically dominated jets are accompanied by the formation of the magnetars.
Since the jets are mildly relativistic, we speculate that they might be the
origin of some observed X-ray flashes.
[80]
oai:arXiv.org:0811.0651 [pdf] - 18185
Effects of Rotation on Standing Accretion Shock Instability in Nonlinear
Phase for Core-Collapse Supernovae
Submitted: 2008-11-05
We studied the effects of rotation on standing accretion shock instability
(SASI) by performing three-dimensional hydrodynamics simulations. Taking into
account a realistic equation of state and neutrino heating/cooling, we prepared
a spherically symmmetric and steady accretion flow through a standing shock
wave onto a proto-neutron star (PNS). When the SASI entered the nonlinear
phase, we imposed uniform rotation on the flow advecting from the outer
boundary of the iron core, whose specific angular momentum was assumed to agree
with recent stellar evolution models. Using spherical harmonics in space and
Fourier decompositions in time, we performed mode analysis of the nonspherical
deformed shock wave to observe rotational effects on the SASI in the nonlinear
phase. We found that rotation imposed on the axisymmetric SASI did not make any
spiral modes and hardly affected sloshing modes, except for steady l=2, m=0
modes. In contrast, rotation imposed on the non-axisymmetric flow increased the
amplitude of spiral modes so that some spiral flows accreting on the PNS were
more clearly formed inside the shock wave than without rotation. The amplitudes
of spiral modes increased significantly with rotation in the progressive
direction.
[81]
oai:arXiv.org:0806.1072 [pdf] - 13308
Impact of Rotation on Neutrino Emission and Relic Neutrino Background
from Population III Stars
Submitted: 2008-06-05
We study the effects of rotation on the neutrino emission from Population III
(Pop III) stars by performing a series of two-dimensional rotational collapse
simulations of Pop III stellar cores. Our results show that rotation enhances
the neutrino luminosities and the average energies of emitted neutrinos. This
is because the thermalized inner core, which is the dominant neutrino source
from Pop III stars, can be enlarged, due to rotational flattening, enough to
extend the inner core outside the neutrinospheres. This is in sharp contrast to
the case of spherical collapse, in which the case of inner core shrinks deeper
inside the neutrinospheres before black hole formation, which hinders the
efficient neutrino emission. In the case of rotational core-collapse, the
emitted neutrino energies are found to become larger in the vicinity near the
pole than the ones near the equatorial plane. These factors make the emergent
neutrino spectrum broader and harder than the spherical collapse case. By
computing the overall neutrino signals produced by the ensemble of individual
rotating Pop III stars, we find that the amplitudes of the relic neutrinos,
depending on their star formation rates, can dominate over the contributions
from ordinary core-collapse supernovae below a few MeV. A detection of this
signal could be an important tool to probe star formation history in the early
universe.
[82]
oai:arXiv.org:0711.0107 [pdf] - 6613
Search method for coincident events from LIGO and IceCube detectors
Submitted: 2007-11-01, last modified: 2008-01-31
We present a coincidence search method for astronomical events using
gravitational wave detectors in conjunction with other astronomical
observations. We illustrate our method for the specific case of the LIGO
gravitational wave detector and the IceCube neutrino detector. LIGO
trigger-events and IceCube events which occur within a given time window are
selected as time-coincident events. Then the spatial overlap of the
reconstructed event directions is evaluated using an unbinned maximum
likelihood method. Our method was tested with Monte Carlo simulations based on
realistic LIGO and IceCube event distributions. We estimated a typical false
alarm rate for the analysis to be 1 event per 435 years. This is significantly
smaller than the false alarm rates of the ndividual detectors.
[83]
oai:arXiv.org:0710.2191 [pdf] - 712792
Three-Dimensional Simulations of Standing Accretion Shock Instability in
Core-Collapse Supernovae
Submitted: 2007-10-11
We have studied non-axisymmetric standing accretion shock instability, or
SASI, by 3D hydrodynamical simulations. This is an extention of our previous
study on axisymmetric SASI. We have prepared a spherically symmetric and steady
accretion flow through a standing shock wave onto a proto-neutron star, taking
into account a realistic equation of state and neutrino heating and cooling.
This unperturbed model is supposed to represent approximately the typical
post-bounce phase of core-collapse supernovae. We then have added a small
perturbation (~1%) to the radial velocity and computed the ensuing evolutions.
Not only axisymmetric but non-axisymmetric perturbations have been also
imposed. We have applied mode analysis to the non-spherical deformation of the
shock surface, using the spherical harmonics. We have found that (1) the growth
rates of SASI are degenerate with respect to the azimuthal index m of the
spherical harmonics Y_l^m, just as expected for a spherically symmetric
background, (2) nonlinear mode couplings produce only m=0 modes for the
axisymmetric perturbations, whereas m=!0 modes are also generated in the
non-axisymmetric cases according to the selection rule for the quadratic
couplings, (3) the nonlinear saturation level of each mode is lower in general
for 3D than for 2D because a larger number of modes are contributing to
turbulence in 3D, (4) low l modes are dominant in the nonlinear phase, (5) the
equi-partition is nearly established among different m modes in the nonlinear
phase, (6) the spectra with respect to l obey power laws with a slope slightly
steeper for 3D, and (7) although these features are common to the models with
and without a shock revival at the end of simulation, the dominance of low l
modes is more remarkable in the models with a shock revival.
[84]
oai:arXiv.org:0709.1795 [pdf] - 4853
Numerical Simulations of Equatorially-Asymmetric Magnetized Supernovae:
Formation of Magnetars and Their Kicks
Submitted: 2007-09-12
A series of numerical simulations on magnetorotational core-collapse
supernovae are carried out. Dipole-like configurations which are offset
northward are assumed for the initially strong magnetic fields together with
rapid differential rotations. Aims of our study are to investigate effects of
the offset magnetic field on magnetar kicks and on supernova dynamics. Note
that we study a regime where the proto-neutron star formed after collapse has a
large magnetic field strength approaching that of a ``magnetar'', a highly
magnetized slowly rotating neutron star. As a result, equatorially-asymmetric
explosions occur with a formation of the bipolar jets. Resultant magnetar's
kick velocities are $\sim 300-1000$ km s$^{-1}$. We find that the acceleration
is mainly due to the magnetic pressure while the somewhat weaker magnetic
tension works toward the opposite direction, which is due to stronger magnetic
field in the northern hemisphere. Noted that observations of magnetar's proper
motions are very scarce, our results supply a prediction for future
observations. Namely, magnetars possibly have large kick velocities, several
hundred km s$^{-1}$, as ordinary neutron stars do, and in an extreme case they
could have those up to 1000 km s$^{-1}$.
[85]
oai:arXiv.org:0708.3597 [pdf] - 354390
Equilibrium Configurations of Strongly Magnetized Neutron Stars with
Realistic Equations of State
Submitted: 2007-08-27
We investigate equilibrium sequences of magnetized rotating stars with four
kinds of realistic equations of state (EOSs) of SLy (Douchin et al.), FPS
(Pandharipande et al.), Shen (Shen et al.), and LS (Lattimer & Swesty).
Employing the Tomimura-Eriguchi scheme to construct the equilibrium
configurations. we study the basic physical properties of the sequences in the
framework of Newton gravity. In addition we newly take into account a general
relativistic effect to the magnetized rotating configurations. With these
computations, we find that the properties of the Newtonian magnetized stars,
e.g., structure of magnetic field, highly depends on the EOSs.
The toroidal magnetic fields concentrate rather near the surface for Shen and
LS EOSs than those for SLy and FPS EOSs. The poloidal fields are also affected
by the toroidal configurations. Paying attention to the stiffness of the EOSs,
we analyze this tendency in detail. In the general relativistic stars, we find
that the difference due to the EOSs becomes small because all the employed EOSs
become sufficiently stiff for the large maximum density, typically greater than
$10^{15}\rm{g} \rm{cm}^{-3}$. The maximum baryon mass of the magnetized stars
with axis ratio $q\sim 0.7$ increases about up to twenty percents for that of
spherical stars. We furthermore compute equilibrium sequences at finite
temperature, which should serve as an initial condition for the hydrodynamic
study of newly-born magnetars. Our results suggest that we may obtain
information about the EOSs from the observation of the masses of magnetars.
[86]
oai:arXiv.org:0706.3495 [pdf] - 260066
Gravitational Wave Background from Population III Stars
Submitted: 2007-06-24
We estimate the gravitational wave (GW) background from Population III (Pop
III) stars using the results from our hydrodynamic simulations (Suwa et al.
2007). We calculate gravitational waveforms from matter motions and anisotropic
neutrino emissions for single Pop III stars. We find that the GW amplitudes
from matter motions are dominant until just after bounce, but those from
neutrinos dominate later on at frequencies below $\sim 10$ Hz in the GW
spectrum. Computing the overall signal produced by the ensemble of such Pop III
stars, we find that the resultant density parameter of the GW background peaks
at the amplitude of $\Omega_\mathrm{GW}\sim 10^{-10}$ in the frequency interval
$\sim 1-10$ Hz. We show that such signals, depending on the formation rate of
Pop III stars, can be within the detection limits of future planned
interferometers such as DECIGO and BBO in the frequency interval of $\sim
0.1-1$ Hz. Our results suggest that the detection of the GW background from Pop
III stars can be an important tool to supply the information about the star
formation history in the early universe.
[87]
oai:arXiv.org:0704.1945 [pdf] - 1000367
Magnetorotational Collapse of Population III Stars
Submitted: 2007-04-16
We perform a series of two-dimensional magnetorotational core-collapse
simulations of Pop III stars. Changing the initial distributions of rotation
and magnetic fields prior to collapse in a parametric manner, we compute 19
models. By so doing, we systematically investigate how rotation and magnetic
fields affect the collapse dynamics and explore how the properties of the
black-hole formations and neutrino emissions could be affected. As for the
microphysics, we employ a realistic equation of state and approximate the
neutrino transfer by a multiflavour leakage scheme. With these computations, we
find that the jet-like explosions are obtained by the magnetodriven shock waves
if the initial magnetic field is as large as $10^{12}$ G. We point out that the
black-hole masses at the formation decrease with the initial field strength, on
the other hand, increase with the initial rotation rates. As for the neutrino
properties, we point out that the degree of the differential rotation plays an
important role to determine which species of the neutrino luminosity is more
dominant than the others. Furthermore, we find that the stronger magnetic
fields make the peak neutrino luminosities smaller, because the magnetic
pressure acts to halt the collapse in the central regions, leading to the
suppression of the releasable gravitational binding energies.
[88]
oai:arXiv.org:astro-ph/0702476 [pdf] - 89530
Effects of QCD phase transition on gravitational radiation from
two-dimensional collapse and bounce of massive stars
Submitted: 2007-02-19
We perform two-dimensional, magnetohydrodynamical core-collapse simulations
of massive stars accompanying the QCD phase transition. We study how the phase
transition affects the gravitational waveforms near the epoch of core-bounce.
As for initial models, we change the strength of rotation and magnetic fields.
Particularly, the degree of differential rotation in the iron core (Fe-core) is
changed parametrically. As for the microphysics, we adopt a phenomenological
equation of state above the nuclear density, including two parameters to change
the hardness before the transition. We assume the first order phase transition,
where the conversion of bulk nuclear matter to a chirally symmetric quark-gluon
phase is described by the MIT bag model. Based on these computations, we find
that the phase transition can make the maximum amplitudes larger up to $\sim$
10 percents than the ones without the phase transition. On the other hand, the
maximum amplitudes become smaller up to $\sim$ 10 percents owing to the phase
transition, when the degree of the differential rotation becomes larger. We
find that even extremely strong magnetic fields $\sim 10^{17}$ G in the
protoneutron star do not affect these results.
[89]
oai:arXiv.org:astro-ph/0607224 [pdf] - 1469051
Gravitational Radiation from Standing Accretion Shock Instability in
Core-Collapse Supernovae
Submitted: 2006-07-11, last modified: 2006-09-22
We present the results of numerical experiments, in which we study how the
asphericities induced by the growth of the standing accretion shock instability
(SASI) produce the gravitational waveforms in the postbounce phase of
core-collapse supernovae. To obtain the neutrino-driven explosions, we
parameterize the neutrino fluxes emitted from the central protoneutron star and
approximate the neutrino transfer by a light-bulb scheme. We find that the
waveforms due to the anisotropic neutrino emissions show the monotonic increase
with time, whose amplitudes are up to two order-of-magnitudes larger than the
ones from the convective matter motions outside the protoneutron stars.
We point out that the amplitudes begin to become larger when the growth of
the SASI enters the nonlinear phase, in which the deformation of the shocks and
the neutrino anisotropy become large. From the spectrum analysis of the
waveforms, we find that the amplitudes from the neutrinos are dominant over the
ones from the matter motions at the frequency below $\sim 100$ Hz, which are
suggested to be within the detection limits of the detectors in the next
generation such as LCGT and the advanced LIGO for a supernova at 10 kpc.
As a contribution to the gravitational wave background, we show that the
amplitudes from this source could be larger at the frequency above $\sim$ 1 Hz
than the primordial gravitational wave backgrounds, but unfortunately,
invisible to the proposed space-based detectors.
[90]
oai:arXiv.org:astro-ph/0606187 [pdf] - 316380
Inelastic Neutrino-Helium Scatterings and Standing Accretion Shock
Instability in Core-Collapse Supernovae
Submitted: 2006-06-08
We present the results of numerical experiments, in which we have
investigated the influence of the inelastic neutrino-helium interactions on the
standing accretion shock instability supposed to occur in the post-bounce
supernova core. The axisymmetric hydrodynamical simulations of accretion flows
through the standing accretion shock wave onto the protoneutron star show that
the interactions are relatively minor and the linear growth of the shock
instability is hardly affected. The extra heating given by the inelastic
reactions becomes important for the shock revival after the instability enters
the non-linear regime, but only when the neutrino luminosity is very close to
the critical value, at which the shock would be revived without the
interactions. We have also studied the dependence of the results on the initial
amplitudes of perturbation and the temperatures of mu and tau neutrinos.
[91]
oai:arXiv.org:astro-ph/0602457 [pdf] - 80033
Magnetohydrodynamic Simulations of A Rotating Massive Star Collapsing to
A Black Hole
Submitted: 2006-02-21
We perform two-dimensional, axisymmetric, magnetohydrodynamic simulations of
the collapse of a rotating star of 40 Msun and in the light of the collapsar
model of gamma-ray burst. Considering two distributions of angular momentum, up
to \sim 10^{17} cm^2/s, and the uniform vertical magnetic field, we investigate
the formation of an accretion disk around a black hole and the jet production
near the hole. After material reaches to the black hole with the high angular
momentum, the disk is formed inside a surface of weak shock. The disk becomes
in a quasi-steady state for stars whose magnetic field is less than 10^{10} G
before the collapse. We find that the jet can be driven by the magnetic fields
even if the central core does not rotate as rapidly as previously assumed and
outer layers of the star has sufficiently high angular momentum. The magnetic
fields are chiefly amplified inside the disk due to the compression and the
wrapping of the field. The fields inside the disk propagate to the polar region
along the inner boundary near the black hole through the Alfv{\'e}n wave, and
eventually drive the jet. The quasi-steady disk is not an advection-dominated
disk but a neutrino cooling-dominated one. Mass accretion rates in the disks
are greater than 0.01 Msun/sec with large fluctuations. The disk is transparent
for neutrinos. The dense part of the disk, which locates near the hole, emits
neutrino efficiently at a constant rate of < 8 \times 10^{51} erg/s. The
neutrino luminosity is much smaller than those from supernovae after the
neutrino burst.
[92]
oai:arXiv.org:astro-ph/0602460 [pdf] - 80036
Heavy element nucleosynthesis in a collapsar
Submitted: 2006-02-21
We investigate synthesis of heavy elements in a collapsar. We have calculated
detailed composition of magnetically driven jets ejected from a collapsar,
which is based on long-term, magneto-hydrodynamic simulations of a rapidly
rotating massive star of 40Msun during core collapse. We follow evolution of
abundances of about 4000 nuclides from the collapse phase to the ejection phase
through the jet generation phase with use of two large nuclear reaction
networks. We find that the r-process successfully operates in the jets, so that
U and Th are synthesized abundantly when the progenitor has large magnetic
field of 10^{12} G and rapidly rotating core. Abundance pattern inside the jets
is similar compared to that of r-elements in the solar system. Heavy
neutron-rich nuclei \sim 0.01Msun can be ejected from the collapsar. The
detailed abundances depend on nuclear properties of mass model, beta-decay
rate, and fission, for nuclei near the neutron drip line. Furthermore, we find
that p-nuclei are produced without seed nuclei: not only light p-nuclei, such
as Se74, Kr78, Sr84, and Mo92, but also heavy p-nuclei, In113, Sn115, and
La138, can be abundantly synthesized in the jets. The amounts of p-nuclei in
the ejecta are much greater than those in core-collapse supernovae (SNe). In
particular, Mo92, In113, Sn115, and La138 deficient in the SNe, are
significantly produced in the ejecta.
[93]
oai:arXiv.org:astro-ph/0510855 [pdf] - 77351
Biermann Mechanism in Primordial Supernova Remnant and Seed Fields
Submitted: 2005-10-31, last modified: 2006-01-24
We have studied the generation of magnetic fields by the Biermann mechanism
in the pair-instability supernovae explosions of the first stars. The Biermann
mechanism produces magnetic fields in the shocked region between the bubble and
interstellar medium (ISM), even if magnetic fields are absent initially. We
have performed a series of two-dimensional magnetohydrodynamic simulations with
the Biermann term and estimate the amplitude and total energy of the produced
magnetic fields. We find that magnetic fields with amplitude
$10^{-14}-10^{-17}$ G are generated inside the bubble, though the amount of
magnetic fields generated depend on specific values of initial conditions. This
corresponds to magnetic fields of $10^{28}-10^{31}$ ergs per each supernova
remnant, which is strong enough to be the seed magnetic field for a galactic
and/or interstellar dynamo.
[94]
oai:arXiv.org:astro-ph/0511826 [pdf] - 78184
Toward Radiation-Magnetohydrodynamic Simulations in Core-Collapse
Supernovae
Submitted: 2005-11-30
We report a current status of our radiation-magnetohydrodynamic code for the
study of core-collapse supernovae. In this contribution, we discuss the
accuracy of our newly developed numerical code by presenting the test problem
in a static background model. We also present the application to the
spherically symmetric core-collapse simulations. Since close comparison with
the previously published models is made, we are now applying it for the study
of magnetorotational core-collapse supernovae.
[95]
oai:arXiv.org:astro-ph/0509765 [pdf] - 386221
Numerical Analysis on Standing Accretion Shock Instability with Neutrino
Heating in the Supernova Cores
Submitted: 2005-09-26
We have numerically studied the instability of the spherically symmetric
standing accretion shock wave against non-spherical perturbations. We have in
mind the application to the collapse-driven supernovae in the post bounce
phase, where the prompt shock wave generated by core bounce is commonly
stalled. We take an experimental stand point in this paper. Using spherically
symmetric, completely steady, shocked accretion flows as unperturbed states, we
have clearly observed both the linear growth and the subsequent nonlinear
saturation of the instability. In so doing, we have employed a realistic
equation of state together with heating and cooling via neutrino reactions with
nucleons. We have done a mode analysis based on the spherical harmonics
decomposition and found that the modes with l=1, 2 are dominant not only in the
linear regime, but also after the nonlinear couplings generate various modes
and the saturation occurs. Varying the neutrino luminosity, we have constructed
the unperturbed states both with and without a negative entropy-gradient. We
have found that in both cases the growth of the instability is similar,
suggesting the convection does not play a dominant role, which also appears to
be supported by the recent linear analysis of the convection in accretion flows
by Foglizzo et al. The real part of the eigen frequency seems to be mainly
determined by the advection time rather than by the sound-crossing time.
Whatever the cause may be, the instability is favorable for the shock revival.
[96]
oai:arXiv.org:astro-ph/0509787 [pdf] - 76368
Gravitational Wave Background from Neutrino-Driven Gamma-Ray Bursts
Submitted: 2005-09-26
We discuss the gravitational wave background (GWB) from a cosmological
population of gamma-ray bursts (GRBs). Among various emission mechanisms for
the gravitational waves (GWs), we pay a particular attention to the vast
anisotropic neutrino emissions from the accretion disk around the black hole
formed after the so-called failed supernova explosions. The produced GWs by
such mechanism are known as burst with memory, which could dominate over the
low-frequency regime below \sim 10Hz. To estimate their amplitudes, we derive
general analytic formulae for gravitational waveform from the axisymmetric
jets. Based on the formulae, we first quantify the spectrum of GWs from a
single GRB. Then, summing up its cosmological population, we find that the
resultant value of the density parameter becomes roughly \Omega_{GW} \approx
10^{-20} over the wide-band of the low-frequency region, f\sim 10^{-4}-10^1Hz.
The amplitude of GWB is sufficiently smaller than the primordial GWBs
originated from an inflationary epoch and far below the detection limit.
[97]
oai:arXiv.org:astro-ph/0509456 [pdf] - 76038
Explosion Mechanism, Neutrino Burst, and Gravitational Wave in
Core-Collapse Supernovae
Submitted: 2005-09-15, last modified: 2005-09-16
Core-collapse supernovae are among the most energetic explosions in the
universe marking the catastrophic end of massive stars. In spite of rigorous
studies for several decades, we still don't understand the explosion mechanism
completely. Since they are related to many astrophysical phenomena such as
nucleosynthesis, gamma-ray bursts and acceleration of cosmic rays,
understanding of their physics has been of wide interest to the astrophysical
community.
In this article, we review recent progress in the study of core-collapse
supernovae focusing on the explosion mechanism, supernova neutrinos, and the
gravitational waves. As for the explosion mechanism, we present a review paying
particular attention to the roles of multidimensional aspects, such as
convection, rotation, and magnetic fields, on the neutrino heating mechanism.
Next, we discuss supernova neutrinos, which is a powerful tool to probe not
only deep inside of the supernovae but also intrinsic properties of neutrinos.
For this purpose, it is necessary to understand neutrino oscillation which has
been established recently by a lot of experiments. Gravitational astronomy is
now also becoming reality.
We present an extensive review on the physical foundations and the emission
mechanism of gravitational waves in detail, and discuss the possibility of
their detections.
[98]
oai:arXiv.org:astro-ph/0501538 [pdf] - 70650
Biermann Mechanism in Primordial Supernova Remnant and Seed Magnetic
Fields
Submitted: 2005-01-25, last modified: 2005-07-11
We study generation of magnetic fields by the Biermann mechanism in the
pair-instability supernovae explosions of first stars. The Biermann mechanism
produces magnetic fields in the shocked region between the bubble and
interstellar medium (ISM), even if magnetic fields are absent initially. We
perform a series of two-dimensional magnetohydrodynamic simulations with the
Biermann term and estimate the amplitude and total energy of the produced
magnetic fields. We find that magnetic fields with amplitude
$10^{-14}-10^{-17}$ G are generated inside the bubble, though the amount of
magnetic fields generated depend on specific values of initial conditions. This
corresponds to magnetic fields of $10^{28}-10^{31}$ erg per each supernova
remnant, which is strong enough to be the seed magnetic field for galactic
and/or interstellar dynamo.
[99]
oai:arXiv.org:astro-ph/0505611 [pdf] - 73381
The Core-Collapse Supernova with "Non-Uniform" Magnetic Fields
Submitted: 2005-05-31, last modified: 2005-06-13
We perform two-dimensional numerical simulations on the core-collapse of a
massive star with strong magnetic fields and differential rotations using a
numerical code ZEUS-2D. Changing field configurations and laws of differential
rotation parametrically, we compute 14 models and investigate effects of these
parameters on the dynamics. In our models, we do not solve the neutrino
transport and instead employ a phenomenological parametric EOS that takes into
account the neutrino emissions. As a result of the calculations, we find that
the field configuration plays a significant role in the dynamics of the core if
the initial magnetic field is large enough. Models with initially concentrated
fields produce more energetic explosions and more prolate shock waves than the
uniform field. Quadrapole-like fields produce remarkably collimated and fast
jet, which might be important for gamma-ray bursts(GRB). The Lorentz forces
exerted in the region where the plasma-beta is less than unity are responsible
for these dynamics. The pure toroidal field, on the other hand, does not lead
to any explosion or matter ejection. This suggests the presupernova models of
Heger et al.(2003), in which toroidal fields are predominant, is
disadvantageous for the magnetorotation-induced supernova considered here.
Models with initially weak magnetic fields do not lead to explosion or matter
ejection, either. In these models magnetic fields play no role as they do not
grow on the timescale considered in this paper so that the magnetic pressure
could be comparable to the matter pressure. This is because the exponential
field growth as expected in MRI is not seen in our models. The magnetic field
is amplified mainly by field-compression and field-wrapping in our simulations.
[100]
oai:arXiv.org:astro-ph/0504100 [pdf] - 72197
R-Process Nucleosynthesis in MHD Jet Explosions of Core-Collapse
Supernovae
Submitted: 2005-04-05
We investigate $r$-process nucleosynthesis during the magnetohydrodynamical
(MHD) explosion of supernova in a massive star of 13 $M_{\odot}$. Contrary to
the case of the spherical explosion, jet-like explosion due to the combined
effects of the rotation and magnetic field lowers the electron fraction
significantly inside the layers above the iron core. We find that the ejected
material of low electron fraction responsible for the $r$-process comes out
from the silicon rich layer of the presupernova model. This leads to the
production up to the third peak in the solar $r$-process elements. We examine
whether the fission affects the $r$-process paths by using the full nuclear
reaction network with both the spontaneous and $\beta$-delayed fission
included. Moreover, we pay particular attention how the mass formula affects
the $r$-process peaks with use of two mass formulae. It is found that both
formulae can reproduce the global abundance pattern up to the third peak though
detailed distributions are rather different. We point out that there are
variations in the $r$-process nucleosynthesis if the MHD effects play an
important role in the supernova explosion.
[101]
oai:arXiv.org:astro-ph/0409244 [pdf] - 265340
North-South Neutrino Heating Asymmetry in Strongly Magnetized and
Rotating Stellar Cores
Submitted: 2004-09-10
We perform a series of two-dimensional magnetohydrodynamic simulations of
supernova cores. Since the distributions of the angular momentum and the
magnetic fields of strongly magnetized stars are quite uncertain, we
systematically change the combinations of the strength of the angular momentum,
the rotations law, the degree of differential rotation, and the profiles of the
magnetic fields to construct the initial conditions. By so doing, we estimate
how the rotation-induced anisotropic neutrino heating are affected by the
strong magnetic fields through parity-violating effects and first investigate
how the north-south asymmetry of the neutrino heating in a strongly magnetized
supernova core could be. As for the microphysics, we employ a realistic
equation of state based on the relativistic mean field theory and take into
account electron captures and the neutrino transport via the neutrino leakage
scheme. With these computations, we find that the parity-violating corrections
reduce $ \lesssim 0.5 %$ of the neutrino heating rate than that without the
magnetic fields in the vicinity of the north pole of a star, on the other hand,
enhance about $ \lesssim 0.5 %$ in the vicinity of the south pole. If the
global asymmetry of the neutrino heating in the both of the poles develops in
the later phases, the newly born neutron star might be kicked toward the north
pole in the subsequent time.
[102]
oai:arXiv.org:astro-ph/0408388 [pdf] - 66879
Magneto-driven Shock Waves in Core-Collapse Supernova
Submitted: 2004-08-21, last modified: 2004-08-25
We perform a series of two-dimensional magnetohydrodynamic simulations of the
rotational core-collapse of a magnetized massive star. We employ a realistic
equation of state and take into account the neutrino cooling by the so-called
leakage scheme. In this study we systematically investigate how the strong
magnetic field and the rapid rotation affect the propagation of the shock
waves. Our results show that in the case of the strong initial poloidal
magnetic field, the toroidal magnetic field amplified by the differential
rotation, becomes strong enough to generate a tightly collimated shock wave
along the rotational axis. On the other hand, in the case of the weak initial
magnetic field, although the differential rotation amplifies toroidal magnetic
field over the long rotational period, the launched shock wave is weak and the
shape of it becomes wider. The former case is expected to be accompanied by the
formation of the so-called magnetar. Our models with rapid rotation and strong
magnetic field can create a nozzle formed by the collimated shock wave. This
might be the analogous situation of the collapsar that is plausible for the
central engine of the Gamma-Ray Bursts.
[103]
oai:arXiv.org:astro-ph/0401563 [pdf] - 265329
Gravitational Radiation from Rotational Core Collapse: Effects of
Magnetic Fields and Realistic Equation of States
Submitted: 2004-01-27
We perform a series of two-dimensional, axisymmetric, magnetohydrodynamic
simulations of the rotational collapse of a supernova core. In order to
calculate the waveforms of the gravitational wave, we derive the quadrupole
formula including the contributions from the electromagnetic fields. Recent
stellar evolution calculations imply that the magnetic fields of the toroidal
components are much stronger than those of the poloidal ones at the
presupernova stage. Thus, we systematically investigate the effects of the
toroidal magnetic fields on the amplitudes and waveforms. Furthermore, we
employ the two kinds of the realistic equation of states, which are often used
in the supernova simulations. Then, we investigate the effects of the equation
of states on the gravitational wave signals. With these computations, we find
that the peak amplitudes are lowered by an order of 10% for the models with the
strongest toroidal magnetic fields. However, the peak amplitudes are mostly
within sensitivity range of laser interferometers such as TAMA and the first
LIGO for a source at a distance of 10 kpc. Furthermore, we point out that the
amplitudes of second peaks are still within the detection limit of the first
LIGO for the source, although the characteristics of second peaks are reduced
by the magnetic fields. We stress the importance of the detection, since it
will give us information about the angular momentum distribution of massive
evolved stars. When we compare the gravitational waves from the two realistic
equation of states, significant differences are not found, except that the
typical frequencies of the gravitational wave become slightly higher for the
softer equation of state.
[104]
oai:arXiv.org:astro-ph/0306020 [pdf] - 57101
Decaying Cold Dark Matter and the Evolution of the Cluster Abundance
Submitted: 2003-06-01, last modified: 2003-07-15
The cluster abundance and its redshift evolution are known to be a powerful
tool to constrain the amplitude of mass fluctuations \sigma_8 and the mass
density parameter \Omega_m0. We study the impact of the finite decay rate of
cold dark matter particles on the cluster abundances. On the basis of spherical
model in decaying cold dark matter universe, we calculate the mass function of
clusters and compare it with observed cluster abundance. We find the decay of
cold dark matter particles significantly changes the evolution of the cluster
abundance. In particular, we point out that the lifetime of dark matter
particles comparable to the age of the universe lowers the ratio of the local
cluster abundance to the high-redshift cluster abundance and can account for
the observed evolution of the cluster abundance quite well. The strong
dependence of the cluster abundance on the decay rate of dark matter suggests
that distant cluster surveys may offer clues to the nature of dark matter.
[105]
oai:arXiv.org:astro-ph/0306430 [pdf] - 1468541
Gravitational Radiation from Axisymmetric Rotational Core Collapse
Submitted: 2003-06-21
We have done a series of two-dimensional hydrodynamic simulations of the
rotational collapse of a supernova core in axisymmetry. We have employed a
realistic equation of state (EOS) and taken into account electron captures and
neutrino transport by the so-called leakage scheme. It is an important progress
to apply the realistic EOS coupled with the microphysics to 2-D simulations for
computing gravitational radiation in rotational core collapse.
We have used the quadrupole formula to calculate the amplitudes and the
waveforms of gravitational wave assuming the Newtonian gravity. From these
computations, we have extended the conventional category of the gravitational
waveforms. Our results have shown that the peak amplitudes of gravitational
wave are mostly within the sensitivity range of the laser interferometers such
as TAMA and first
LIGO for a source at a distance of 10 kpc. Furthermore we have found that the
amplitudes of the second peaks are within the detection limit of first LIGO for
the source and first pointed out the importance of the detections, since they
will give us the information as to the angular momentum distribution of evolved
massive stars.
[106]
oai:arXiv.org:astro-ph/0305260 [pdf] - 56733
Probing dark energy with gamma-ray bursts
Submitted: 2003-05-15, last modified: 2003-06-01
We propose a new method to use gamma-ray bursts (GRBs) as an alternative
probe of the dark energy. By calibrating luminosity-variability and
luminosity-lag time relations at low redshift where distance-redshift relations
have been already determined from type Ia supernova (SNIa), GRBs at high
redshift can be used as a distance indicator. We investigate the potential
impact of future GRB data on determining the current matter density
$\Omega_{m}$ and the dark energy equation of state $w$ which is assumed to be
constant. We show that a combined analysis of a data set expected from the {\it
Swift} and the current SNIa data results in excellent determination of both
$\Omega_{m}$ and $w$.
[107]
oai:arXiv.org:astro-ph/0104485 [pdf] - 1232673
Specifying the Environments around GRB, Explaining the Fe line in the
X-Ray Afterglow of GRB000214
Submitted: 2001-04-30, last modified: 2001-06-23
We present a model explaining the Fe K alpha line and the continuum in the
afterglow of GRB000214. We pose the importance to seek the physically natural
environment around GRB000214. For the reproduction of the observation, we need
the ring-like remnant around the progenitor like that of SN 1987A produced by
the mass-loss of the progenitor and the fireball spread over in every
directions. The observation of GRB000214, in which the continuum power law
spectrum decreased faster than the line, motivated us to consider the two
independent systems for the line emission and the continuum spectrum. At first,
the continuum spectrum can be fitted by the afterglow emission of the fireball
pointing toward the observer which does not collide with the ring because the
emission of GRB and the afterglow are highly collimated to the observer by the
relativistic beaming effect. Secondly, the line can be fitted by the
fluorescence of the Fe atoms in the ring illuminated by the X-ray afterglow.
The significance of this study is that our model may constrain strongly the GRB
model. Although the Supranova model assumes the extreme-ring-like remnant
produced by the usual supernova explosion, this may not be probable. It is
because the supernova remnants are known to be shell-like. The model also
assumes two steps of explosions, on the other hand, we need only one explosion
of the progenitor. In this sense, our scenario is more natural. Moreover, in
the numerical simulations of Hypernova, the jet of the opening angle of only 1
degree is generated. In our model, the fireball which spreads over in every
directions reconciles with the observation of 1 percent of the polarization in
the observation of SN1998bw which showed the explosion might not be so
collimated.
