Normalized to: Sekiguchi, Y.
[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:2007.00474 [pdf] - 2126356
Post-merger Mass Ejection of Low-mass Binary Neutron Stars
Submitted: 2020-07-01
We study the post-merger mass ejection of low-mass binary neutron stars (NSs)
with the system mass of $2.5\, M_\odot$, and subsequent nucleosynthesis by
performing general-relativistic, neutrino-radiation viscous-hydrodynamics
simulations in axial symmetry. We find that the merger remnants are long-lived
massive NSs surviving more than several seconds, irrespective of the nuclear
equations of state (EOSs) adopted. The ejecta masses of our fiducial models are
$\sim 0.06$-$0.1\, M_\odot$ (depending on the EOS), being $\sim 30\%$ of the
initial disk masses ($\sim 0.15$-$0.3\, M_\odot$). Post-processing
nucleosynthesis calculations indicate that the ejecta is composed mainly of
light $r$-process nuclei with small amounts of lanthanides (mass fraction $\sim
0.002$-$0.004$) and heavier species due to the modest average electron fraction
($\sim 0.32$-$0.34$) for a reasonable value of the viscous coefficient. Such
abundance distributions are incompatible with the solar $r$-process-like
abundance patterns found in all measured $r$-process-enhanced metal-poor stars.
Therefore, low-mass binary NS mergers should be rare. If such low-mass NS
mergers occur, their electromagnetic counterparts, kilonovae, will be
characterized by an early bright blue emission because of the large ejecta mass
as well as the small lanthanide fraction. We also show, however, that if the
effective turbulent viscosity is very high, or there is an efficient mass
ejection working in the early post-merger phase, the electron fraction of the
ejecta could be low enough that the solar $r$-process-like abundance pattern is
reproduced and the lanthanide fraction becomes so high that the kilonova would
be characterized by early bright blue and late bright red emissions.
[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:1907.03790 [pdf] - 2076596
Sub-radian-accuracy gravitational waves from coalescing binary neutron
stars II: Systematic study on the equation of state, binary mass, and mass
ratio
Submitted: 2019-07-08, last modified: 2020-03-12
We report results of numerical relativity simulations for {\it new} 26
non-spinning binary neutron star systems with 6 grid resolutions using an
adaptive mesh refinement numerical re\ lativity code {\tt SACRA-MPI}. The
finest grid spacing is $\approx 64$--$85$ m, depending on the systems. First,
we derive long-term high-precision inspiral gravitational waveforms and show
that the accumulated gravitational-wave phase error due to the finite grid
resolution is less than $0.5$ rad during more than $200$ rad phase evolution
irrespective of the systems. We also find that the gravitational-wave phase
error for a binary system with a tabulated equation of state (EOS) is
comparable to that for a piecewise polytropic EOS. Then we validate the SACRA
inspiral gravitational waveform template, which will be used to extract tidal
deformability from gravitational wave observation, and find that accuracy of \
our waveform modeling is $\lesssim 0.1$ rad in the gravitational-wave phase and
$\lesssim 20 \%$ in the gravitational-wave amplitude up to the
gravitational-wave frequency $1000$ Hz.\
Finally, we calibrate the proposed universal relations between a post-merger
gravitational wave signal and tidal deformability/neutron star radius in the
literature and show that th\ ey suffer from systematics and many relations
proposed as universal are not very universal. Improved fitting formulae are
also proposed.
[6]
oai:arXiv.org:2001.04467 [pdf] - 2085105
Mass ejection from disks surrounding a low-mass black hole: Viscous
neutrino-radiation hydrodynamics simulation in full general relativity
Submitted: 2020-01-13
New viscous neutrino-radiation hydrodynamics simulations are performed for
accretion disks surrounding a spinning black hole with low mass $3M_\odot$ and
dimensionless spin 0.8 or 0.6 in full general relativity, aiming at modeling
the evolution of a merger remnant of massive binary neutron stars or low-mass
black hole-neutron star binaries. We reconfirm the following results found by
previous studies of other groups: 15-30% of the disk mass is ejected from the
system with the average velocity of $\sim $5-10% of the speed of light for the
plausible profile of the disk as merger remnants. In addition, we find that for
the not extremely high viscous coefficient case, the neutron richness of the
ejecta does not become very high, because weak interaction processes enhance
the electron fraction during the viscous expansion of the disk before the onset
of the mass ejection, resulting in the suppression of the lanthanide synthesis.
For high-mass disks, the viscous expansion timescale is increased by a
longer-term neutrino emission, and hence, the electron fraction of the ejecta
becomes even higher. We also confirm that the mass distribution of the electron
fraction depends strongly on the magnitude of the given viscous coefficient.
This demonstrates that a first-principle magnetohydrodynamics simulation is
necessary for black hole-disk systems with sufficient grid resolution and with
sufficiently long timescale (longer than seconds) to clarify the
nucleosynthesis and electromagnetic signals from them.
[7]
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.
[8]
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.
[9]
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.
[10]
oai:arXiv.org:1711.02093 [pdf] - 1705157
Mass Ejection from the Remnant of a Binary Neutron Star Merger:
Viscous-Radiation Hydrodynamics Study
Submitted: 2017-11-06, last modified: 2018-06-05
We perform long-term general relativistic neutrino radiation hydrodynamics
simulations (in axisymmetry) for a massive neutron star (MNS) surrounded by a
torus, which is a canonical remnant formed after the binary neutron star
merger. We take into account the effects of viscosity, which is likely to arise
in the merger remnant due to magnetohydrodynamical turbulence. As the initial
condition, we employ the azimuthally averaged data of the MNS-torus system
derived in a three-dimensional, numerical-relativity simulation for the binary
neutron star merger. The viscous effect plays key roles for the remnant
evolution and mass ejection from it in two phases of the evolution. In the
first $t\lesssim10$ ms, a differential rotation state of the MNS is changed to
a rigidly rotating state, and as a result, a sound wave, which subsequently
becomes a shock wave, is formed in the vicinity of the MNS due to the variation
of the quasi-equilibrium state of the MNS. The shock wave induces significant
mass ejection of mass $\sim(0.5-2.0)\times 10^{-2}M_\odot$ for the alpha
viscosity parameter of $0.01-0.04$. For the longer-term evolution with $\sim
0.1-10$ s, a significant fraction of the torus material is ejected. The ejecta
mass is likely to be of order $10^{-2}M_\odot$, so that the total mass of the
viscosity-driven ejecta could dominate that of the dynamical ejecta of mass
$\lesssim 10^{-2}M_\odot$. The electron fraction, $Y_e$, of the ejecta is
always high enough ($Y_e\gtrsim0.25$) that this post-merger ejecta is
lanthanide-poor; hence, the opacity of the ejecta is likely to be $\sim 10-100$
times lower than that of the dynamical ejecta. This indicates that the
electromagnetic signal from the ejecta would be rapidly evolving, bright, and
blue if it is observed from a small viewing angle ($\lesssim 45^\circ$) for
which the effect of the dynamical ejecta is minor.
[11]
oai:arXiv.org:1710.01311 [pdf] - 1701512
Global simulations of strongly magnetized remnant massive neutron stars
formed in binary neutron star mergers
Submitted: 2017-10-03, last modified: 2018-05-09
We perform a general-relativistic magnetohydrodynamics simulation for
$\approx 30$ ms after merger of a binary neutron star to a remnant massive
neutron star (RMNS) with a high spatial resolution of the finest grid
resolution $12.5$ m. First, we estimate that the Kelvin-Helmholtz instability
at merger could amplify the magnetic-field energy up to $\sim 1\%$ of the
thermal energy. Second, we find that the magnetorotational instability in the
RMNS envelope and torus with $\rho < 10^{13}~{\rm g~cm^{-3}}$ sustains
magneto-turbulent state and the effective viscous parameter in these regions is
likely to converge to $\approx 0.01$--$0.02$ with respect to the grid
resolution. Third, the current grid resolution is not still fine enough to
sustain magneto-turbulent state in the RMNS with $\rho \ge 10^{13}~{\rm
g~cm^{-3}}$.
[12]
oai:arXiv.org:1802.06518 [pdf] - 1659716
Frequency-domain gravitational waveform models for inspiraling binary
neutron stars
Submitted: 2018-02-19
We develop a model for frequency-domain gravitational waveforms from
inspiraling binary neutron stars. Our waveform model is calibrated by
comparison with hybrid waveforms constructed from our latest high-precision
numerical-relativity waveforms and the SEOBNRv2T waveforms in the frequency
range of $10$--$1000\,{\rm Hz}$. We show that the phase difference between our
waveform model and the hybrid waveforms is always smaller than $0.1\, {\rm
rad}$ for the binary tidal deformability, ${\tilde \Lambda}$, in the range
$300\lesssim{\tilde \Lambda}\lesssim1900$ and for the mass ratio between 0.73
and 1. We show that, for $10$--$1000\,{\rm Hz}$, the distinguishability for the
signal-to-noise ratio $\lesssim50$ and the mismatch between our waveform model
and the hybrid waveforms are always smaller than 0.25 and $1.1\times10^{-5}$,
respectively. The systematic error of our waveform model in the measurement of
${\tilde \Lambda}$ is always smaller than $20$ with respect to the hybrid
waveforms for $300\lesssim{\tilde \Lambda}\lesssim1900$. The statistical error
in the measurement of binary parameters is computed employing our waveform
model, and we obtain results consistent with the previous studies. We show that
the systematic error of our waveform model is always smaller than $20\%$
(typically smaller than $10\%$) of the statistical error for events with the
signal-to-noise ratio of $50$.
[13]
oai:arXiv.org:1710.00827 [pdf] - 1621341
Neutrino transport in black hole-neutron star binaries: neutrino
emission and dynamical mass ejection
Submitted: 2017-10-02, last modified: 2018-01-21
We study the merger of black hole-neutron star binaries by fully
general-relativistic neutrino-radiation-hydrodynamics simulations throughout
the coalescence, particularly focusing on the role of neutrino irradiation in
dynamical mass ejection. Neutrino transport is incorporated by an approximate
transfer scheme based on the truncated moment formalism. While we fix the mass
ratio of the black hole to the neutron star to be 4 and the dimensionless spin
parameter of the black hole to be 0.75, the equations of state for
finite-temperature neutron-star matter are varied. The hot accretion disk
formed after tidal disruption of the neutron star emits a copious amount of
neutrinos with the peak total luminosity ~1--3x10^53 erg s^(-1) via thermal
pair production and subsequent electron/positron captures on free nucleons.
Nevertheless, the neutrino irradiation does not modify significantly the
electron fraction of the dynamical ejecta from the neutrinoless
beta-equilibrium value at zero temperature of initial neutron stars. The mass
of the wind component driven by neutrinos from the remnant disk is negligible
compared to the very neutron-rich dynamical component, throughout our
simulations performed until a few tens milliseconds after the onset of merger,
for the models considered in this study. These facts suggest that the ejecta
from black hole-neutron star binaries are very neutron rich and are expected to
accommodate strong r-process nucleosynthesis, unless magnetic or viscous
processes contribute substantially to the mass ejection from the disk. We also
find that the peak neutrino luminosity does not necessarily increase as the
disk mass increases, because tidal disruption of a compact neutron star can
result in a remnant disk with a small mass but high temperature.
[14]
oai:arXiv.org:1710.05865 [pdf] - 1648658
Subaru Hyper Suprime-Cam Survey for An Optical Counterpart of GW170817
Tominaga, Nozomu;
Tanaka, Masaomi;
Morokuma, Tomoki;
Utsumi, Yousuke;
Yamaguchi, Masaki S.;
Yasuda, Naoki;
Tanaka, Masayuki;
Yoshida, Michitoshi;
Fujiyoshi, Takuya;
Furusawa, Hisanori;
Kawabata, Koji S.;
Lee, Chien-Hsiu;
Motohara, Kentaro;
Ohsawa, Ryou;
Ohta, Kouji;
Terai, Tsuyoshi;
Abe, Fumio;
Aoki, Wako;
Asakura, Yuichiro;
Barway, Sudhanshu;
Bond, Ian A.;
Fujisawa, Kenta;
Honda, Satoshi;
Ioka, Kunihito;
Itoh, Youichi;
Kawai, Nobuyuki;
Kim, Ji Hoon;
Koshimoto, Naoki;
Matsubayashi, Kazuya;
Miyazaki, Shota;
Saito, Tomoki;
Sekiguchi, Yuichiro;
Sumi, Takahiro;
Tristram, Paul J.
Submitted: 2017-10-16, last modified: 2018-01-07
We perform a $z$-band survey for an optical counterpart of a binary neutron
star coalescence GW170817 with Subaru/Hyper Suprime-Cam. Our untargeted
transient search covers $23.6$ deg$^2$ corresponding to the $56.6\%$ credible
region of GW170817 and reaches the $50\%$ completeness magnitude of $20.6$ mag
on average. As a result, we find 60 candidates of extragalactic transients,
including J-GEM17btc (a.k.a. SSS17a/DLT17ck). While J-GEM17btc is associated
with NGC 4993 that is firmly located inside the 3D skymap of GW170817, the
other 59 candidates do not have distance information in the GLADE v2 catalog or
NASA/IPAC Extragalactic Database (NED). Among 59 candidates, 58 are located at
the center of extended objects in the Pan-STARRS1 catalog, while one candidate
has an offset. We present location, $z$-band apparent magnitude, and time
variability of the candidates and evaluate the probabilities that they are
located inside of the 3D skymap of GW170817. The probability for J-GEM17btc is
$64\%$ being much higher than those for the other 59 candidates
($9.3\times10^{-3}-2.1\times10^{-1}\%$). Furthermore, the possibility, that at
least one of the other 59 candidates is located within the 3D skymap, is only
$3.2\%$. Therefore, we conclude that J-GEM17btc is the most-likely and
distinguished candidate as the optical counterpart of GW170817.
[15]
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.
[16]
oai:arXiv.org:1710.07579 [pdf] - 1608526
Modeling GW170817 based on numerical relativity and its implications
Submitted: 2017-10-20, last modified: 2017-12-04
Gravitational-wave observation together with a large number of
electromagnetic observations shows that the source of the latest
gravitational-wave event, GW170817, detected primarily by advanced LIGO, is the
merger of a binary neutron star. We attempt to interpret this observational
event based on our results of numerical-relativity simulations performed so far
paying particular attention to the optical and infra-red observations. We
finally reach a conclusion that this event is described consistently by the
presence of a long-lived hypermassive or supramassive neutron star as the
merger remnant, because (i) significant contamination by lanthanide elements
along our line of sight to this source can be avoided by the strong neutrino
irradiation from it and (ii) it could play a crucial role to produce an ejecta
component of appreciable mass with fast motion in the post-merger phase. We
also point out that (I) the neutron-star equation of state has to be
sufficiently stiff (i.e., the maximum mass of cold spherical neutron stars,
M_max, has to be appreciably higher than 2M_sun in order that a long-lived
massive neutron star can be formed as the merger remnant for the binary systems
of GW170817, for which the initial total mass is >~ 2.73M_sun and (II) no
detection of relativistic optical counterpart suggests a not-extremely high
value of M_max approximately as 2.15-2.25M_sun.
[17]
oai:arXiv.org:1710.05848 [pdf] - 1589836
J-GEM observations of an electromagnetic counterpart to the neutron star
merger GW170817
Utsumi, Yousuke;
Tanaka, Masaomi;
Tominaga, Nozomu;
Yoshida, Michitoshi;
Barway, Sudhanshu;
Nagayama, Takahiro;
Zenko, Tetsuya;
Aoki, Kentaro;
Fujiyoshi, Takuya;
Furusawa, Hisanori;
Kawabata, Koji S.;
Koshida, Shintaro;
Lee, Chien-Hsiu;
Morokuma, Tomoki;
Motohara, Kentaro;
Nakata, Fumiaki;
Ohsawa, Ryou;
Ohta, Kouji;
Okita, Hirofumi;
Tajitsu, Akito;
Tanaka, Ichi;
Terai, Tsuyoshi;
Yasuda, Naoki;
Abe, Fumio;
Asakura, Yuichiro;
Bond, Ian A.;
Miyazaki, Shota;
Sumi, Takahiro;
Tristram, Paul J.;
Honda, Satoshi;
Itoh, Ryosuke;
Itoh, Yoichi;
Kawabata, Miho;
Morihana, Kumiko;
Nagashima, Hiroki;
Nakaoka, Tatsuya;
Ohshima, Tomohito;
Takahashi, Jun;
Takayama, Masaki;
Aoki, Wako;
Baar, Stefan;
Doi, Mamoru;
Finet, François;
Kanda, Nobuyuki;
Kawai, Nobuyuki;
Kim, Ji Hoon;
Kuroda, Daisuke;
Liu, Wei;
Matsubayashi, Kazuya;
Murata, Katsuhiro L.;
Nagai, Hiroshi;
Saito, Tomoki;
Saito, Yoshihiko;
Sako, Shigeyuki;
Sekiguchi, Yuichiro;
Tamura, Yoichi;
Tanaka, Masayuki;
Uemura, Makoto;
Yamaguchi, Masaki S.;
collaboration, the J-GEM
Submitted: 2017-10-16
The first detected gravitational wave from a neutron star merger was
GW170817. In this study, we present J-GEM follow-up observations of SSS17a, an
electromagnetic counterpart of GW170817. SSS17a shows a 2.5-mag decline in the
$z$-band from 1.7 days to 7.7 days after the merger. Such a rapid decline is
not comparable with supernovae light curves at any epoch. The color of SSS17a
also evolves rapidly and becomes redder for later epochs; the $z-H$ color
changed by approximately 2.5 mag in the period of 0.7 days to 7.7 days. The
rapid evolution of both the optical brightness and the color are consistent
with the expected properties of a kilonova that is powered by the radioactive
decay of newly synthesized $r$-process nuclei. Kilonova models with Lanthanide
elements can reproduce the aforementioned observed properties well, which
suggests that $r$-process nucleosynthesis beyond the second peak takes place in
SSS17a. However, the absolute magnitude of SSS17a is brighter than the expected
brightness of the kilonova models with the ejecta mass of 0.01 $\Msun$, which
suggests a more intense mass ejection ($\sim 0.03 \Msun$) or possibly an
additional energy source.
[18]
oai:arXiv.org:1710.05850 [pdf] - 1589838
Kilonova from post-merger ejecta as an optical and near-infrared
counterpart of GW170817
Tanaka, Masaomi;
Utsumi, Yousuke;
Mazzali, Paolo A.;
Tominaga, Nozomu;
Yoshida, Michitoshi;
Sekiguchi, Yuichiro;
Morokuma, Tomoki;
Motohara, Kentaro;
Ohta, Kouji;
Kawabata, Koji S.;
Abe, Fumio;
Aoki, Kentaro;
Asakura, Yuichiro;
Baar, Stefan;
Barway, Sudhanshu;
Bond, Ian A.;
Doi, Mamoru;
Fujiyoshi, Takuya;
Furusawa, Hisanori;
Honda, Satoshi;
Itoh, Yoichi;
Kawabata, Miho;
Kawai, Nobuyuki;
Kim, Ji Hoon;
Lee, Chien-Hsiu;
Miyazaki, Shota;
Morihana, Kumiko;
Nagashima, Hiroki;
Nagayama, Takahiro;
Nakaoka, Tatsuya;
Nakata, Fumiaki;
Ohsawa, Ryou;
Ohshima, Tomohito;
Okita, Hirofumi;
Saito, Tomoki;
Sumi, Takahiro;
Tajitsu, Akito;
Takahashi, Jun;
Takayama, Masaki;
Tamura, Yoichi;
Tanaka, Ichi;
Terai, Tsuyoshi;
Tristram, Paul J.;
Yasuda, Naoki;
Zenko, Tetsuya
Submitted: 2017-10-16
Recent detection of gravitational waves from a neutron star (NS) merger event
GW170817 and identification of an electromagnetic counterpart provide a unique
opportunity to study the physical processes in NS mergers. To derive properties
of ejected material from the NS merger, we perform radiative transfer
simulations of kilonova, optical and near-infrared emissions powered by
radioactive decays of r-process nuclei synthesized in the merger. We find that
the observed near-infrared emission lasting for > 10 days is explained by 0.03
Msun of ejecta containing lanthanide elements. However, the blue optical
component observed at the initial phases requires an ejecta component with a
relatively high electron fraction (Ye). We show that both optical and
near-infrared emissions are simultaneously reproduced by the ejecta with a
medium Ye of ~ 0.25. We suggest that a dominant component powering the emission
is post-merger ejecta, which exhibits that mass ejection after the first
dynamical ejection is quite efficient. Our results indicate that NS mergers
synthesize a wide range of r-process elements and strengthen the hypothesis
that NS mergers are the origin of r-process elements in the Universe.
[19]
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.
[20]
oai:arXiv.org:1708.08926 [pdf] - 1587642
Sub-radian-accuracy gravitational waveforms of coalescing binary neutron
stars in numerical relativity
Submitted: 2017-08-29
Extending our previous studies, we perform high-resolution simulations of
inspiraling binary neutron stars in numerical relativity. We thoroughly carry
through a convergence study in our currently available computational resources
with the smallest grid spacing of $\approx 63$--86~meter for the neutron-star
radius 10.9--13.7\,km. The estimated total error in the gravitational-wave
phase is of order 0.1~rad for the total phase of $\gtrsim 210$\,rad in the last
$\sim 15$--16 inspiral orbits. We then compare the waveforms (without
resolution extrapolation) with those calculated by the latest
effective-one-body formalism (tidal SEOBv2 model referred to as TEOB model). We
find that for any of our models of binary neutron stars, the waveforms
calculated by the TEOB formalism agree with the numerical-relativity waveforms
up to $\approx 3$\,ms before the peak of the gravitational-wave amplitude is
reached: For this late inspiral stage, the total phase error is $\lesssim
0.1$\,rad. Although the gravitational waveforms have an inspiral-type feature
for the last $\sim 3$\,ms, this stage cannot be well reproduced by the current
TEOB formalism, in particular, for neutron stars with large tidal deformability
(i.e., lager radius). The reason for this is described.
[21]
oai:arXiv.org:1708.09101 [pdf] - 1626297
Properties of Kilonovae from Dynamical and Post-Merger Ejecta of Neutron
Star Mergers
Tanaka, Masaomi;
Kato, Daiji;
Gaigalas, Gediminas;
Rynkun, Pavel;
Radziute, Laima;
Wanajo, Shinya;
Sekiguchi, Yuichiro;
Nakamura, Nobuyuki;
Tanuma, Hajime;
Murakami, Izumi;
Sakaue, Hiroyuki A.
Submitted: 2017-08-29
Ejected material from neutron star mergers give rise to electromagnetic
emission powered by radioactive decays of r-process nuclei, which is so called
kilonova or macronova. While properties of the emission are largely affected by
opacities in the ejected material, available atomic data for r-process elements
are still limited. We perform atomic structure calculations for r-process
elements: Se (Z=34), Ru (Z=44), Te (Z=52), Ba (Z=56), Nd (Z=60), and Er (Z=68).
We confirm that the opacities from bound-bound transitions of open f-shell,
Lanthanide elements (Nd and Er) are higher than those of the other elements
over a wide wavelength range. The opacities of open s-shell (Ba), p-shell (Se
and Te), and d-shell (Ru) elements are lower than those of open f-shell
elements and their transitions are concentrated in the ultraviolet wavelengths.
We show that the optical brightness can be different by >2 mag depending on the
element abundances in the ejecta such that post-merger, Lanthanide-free ejecta
produce brighter and bluer optical emission. Such blue emission from
post-merger ejecta can be observed from the polar directions if the mass of the
preceding dynamical ejecta in these regions is small. For the ejecta mass of
0.01 Msun, observed magnitudes of the blue emission will reach 21.0 mag (100
Mpc) and 22.5 mag (200 Mpc) in g and r bands within a few days after the
merger, which are detectable with 1m or 2m-class telescopes.
[22]
oai:arXiv.org:1703.10191 [pdf] - 1720947
Properties of Neutrino-driven Ejecta from the Remnant of Binary Neutron
Star Merger : Purely Radiation Hydrodynamics Case
Submitted: 2017-03-29, last modified: 2017-07-17
We performed general relativistic, long-term, axisymmetric neutrino radiation
hydrodynamics simulations for the remnant formed after the binary neutron star
merger, which consist of a massive neutron star and a torus surrounding it. As
an initial condition, we employ the result derived in a three-dimensional,
numerical relativity simulation for the binary neutron star merger. We
investigate the properties of neutrino-driven ejecta. Due to the
pair-annihilation heating, the dynamics of the neutrino-driven ejecta is
significantly modified. The kinetic energy of the ejecta is about two times
larger than that in the absence of the pair-annihilation heating. This suggests
that the pair-annihilation heating plays an important role in the evolution of
the merger remnants. The relativistic outflow, which is required for driving
gamma-ray bursts, is not observed because the specific heating rate around the
rotational axis is not sufficiently high due to the baryon loading caused by
the neutrino-driven ejecta from the massive neutron star. We discuss the
condition for launching the relativistic outflow and the nucleosynthesis in the
ejecta.
[23]
oai:arXiv.org:1704.00433 [pdf] - 1582303
Gravitational Collapse of Rotating Supermassive Stars including Nuclear
Burning Effects
Submitted: 2017-04-03
Supermassive stars (SMSs) of mass $\gtrsim 10^5 M_\odot$ are candidates for
seeds of supermassive black holes found in the center of many massive galaxies.
We simulate the gravitational collapse of a rigidly rotating SMS core including
nuclear burning effects in axisymmetric numerical-relativity simulation. We
find that for realistic initial conditions, the nuclear burning does not play
an important role. After the collapse, a torus surrounding a rotating black
hole is formed and a fraction of the torus material is ejected. We
quantitatively study the relation between the properties of these objects and
rotation. We find that if a SMS core is sufficiently rapidly rotating, the
torus and outflow mass have approximately $6\%$ and $1\%$ of the initial mass,
respectively. The typical average velocity and the total kinetic energy of the
outflow are $0.2~c$ and $10^{54-56}$ erg where $c$ is the speed of light.
Finally, we briefly discuss the possibility for observing the torus and
outflow.
[24]
oai:arXiv.org:1703.10303 [pdf] - 1567753
General relativistic viscous hydrodynamics of differentially rotating
neutron stars
Submitted: 2017-03-29
Employing a simplified version of the Israel-Stewart formalism for
general-relativistic shear-viscous hydrodynamics, we perform axisymmetric
general-relativistic simulations for a rotating neutron star surrounded by a
massive torus, which can be formed from differentially rotating stars. We show
that with our choice of a shear-viscous hydrodynamics formalism, the
simulations can be stably performed for a long time scale. We also demonstrate
that with a possibly high shear-viscous coefficient, not only viscous angular
momentum transport works but also an outflow could be driven from a hot
envelope around the neutron star for a time scale $\gtrsim 100\,$ms with the
ejecta mass $\gtrsim 10^{-2}M_\odot$ which is comparable to the typical mass
for dynamical ejecta of binary neutron star mergers. This suggests that massive
neutron stars surrounded by a massive torus, which are typical outcomes formed
after the merger of binary neutron stars, could be the dominant source for
providing neutron-rich ejecta, if the effective shear viscosity is sufficiently
high, i.e., if the viscous $\alpha$ parameter is $\gtrsim 10^{-2}$. The present
numerical result indicates the importance of a future high-resolution
magnetohydrodynamics simulation that is the unique approach to clarify the
viscous effect in the merger remnants of binary neutron stars by the
first-principle manner.
[25]
oai:arXiv.org:1606.07147 [pdf] - 1447824
Gravitational waves from supermassive stars collapsing to a supermassive
black hole
Submitted: 2016-06-22
We derive the gravitational waveform from the collapse of a rapidly rotating
supermassive star (SMS) core leading directly to a seed of a supermassive black
hole (SMBH) in axisymmetric numerical-relativity simulations. We find that the
peak strain amplitude of gravitational waves emitted during the black-hole
formation is $\approx 5 \times 10^{-21}$ at the frequency $f \approx 5$\,mHz
for an event at the cosmological redshift $z=3$, if the collapsing SMS core is
in the hydrogen-burning phase. Such gravitational waves will be detectable by
space laser interferometric detectors like eLISA with signal-to-noise ratio
$\approx 10$, if the sensitivity is as high as LISA for $f=1$--10\,mHz. The
detection of the gravitational-wave signal will provide a potential opportunity
for testing the direct-collapse scenario for the formation of a seed of SMBHs.
[26]
oai:arXiv.org:1605.04992 [pdf] - 1443976
Formation of Overheated Regions and Truncated Disks around Black Holes;
Three-dimensional General Relativistic Radiation-magnetohydrodynamics
Simulations
Submitted: 2016-05-16
Using three-dimensional general relativistic radiation magnetohydrodynamics
simulations of accretion flows around stellar mass black holes, we report that
the relatively cold disk ($\gtrsim 10^{7}$K) is truncated near the black hole.
Hot and less-dense regions, of which the gas temperature is $ \gtrsim 10^9$K
and more than ten times higher than the radiation temperature (overheated
regions), appear within the truncation radius. The overheated regions also
appear above as well as below the disk, and sandwich the cold disk, leading to
the effective Compton upscattering. The truncation radius is $\sim 30 r_{\rm
g}$ for $\dot{M} \sim L_{\rm Edd}/c^2$, where $r_{\rm g}, \dot M,
L_\mathrm{Edd}, c$ are the gravitational radius, mass accretion rate, Eddington
luminosity, and light speed. Our results are consistent with observations of
very high state, whereby the truncated disk is thought to be embedded in the
hot rarefied regions. The truncation radius shifts inward to $\sim 10 r_{\rm
g}$ with increasing mass accretion rate $\dot{M} \sim 100 L_{\rm Edd}/c^2$,
which is very close to an innermost stable circular orbit. This model
corresponds to the slim disk state observed in ultra luminous X-ray sources.
Although the overheated regions shrink if the Compton cooling effectively
reduces the gas temperature, the sandwich-structure does not disappear at the
range of $\dot{M} \lesssim 100L_{\rm Edd}/c^2$. Our simulations also reveal
that the gas temperature in the overheated regions depends on black hole spin,
which would be due to efficient energy transport from black hole to disks
through the Poynting flux, resulting gas heating.
[27]
oai:arXiv.org:1604.00643 [pdf] - 1389011
Stability of rigidly rotating supermassive stars against gravitational
collapse
Submitted: 2016-04-03
We revisit secular stability against quasi-radial collapse for rigidly
rotating supermassive stars (SMSs) in general relativity. We suppose that the
SMSs are in a nuclear-burning phase and can be modeled by polytropic equations
of state with the polytropic index $n_p$ slightly smaller than $3$. The
stability is determined in terms of the turning-point method. We find a fitting
formula of the stability condition for the plausible range of $n_p$ ($2.95 \alt
n_p \alt 3$) for SMSs. This condition reconfirms that, while non-rotating SMSs
with mass $\sim 10^5M_\odot$--$10^6M_\odot$ may undergo a
general-relativistically induced quasi-radial collapse, rigidly rotating SMSs
with a ratio of rotational to gravitational potential energy ($\beta$) of $\sim
10^{-2}$ are likely to be stable against collapse unless they are able to
accrete $\sim 5$ times more mass during the (relatively brief) hydrogen-burning
phase of their evolution. We discuss implications of our results.
[28]
oai:arXiv.org:1603.01918 [pdf] - 1426275
Dynamical mass ejection from the merger of asymmetric binary neutron
stars: Radiation-hydrodynamics study in general relativity
Submitted: 2016-03-06
We perform neutrino radiation-hydrodynamics simulations for the merger of
asymmetric binary neutron stars in numerical relativity. Neutron stars are
modeled by soft and moderately stiff finite-temperature equations of state
(EOS). We find that the properties of the dynamical ejecta such as the total
mass, neutron richness profile, and specific entropy profile depend on the mass
ratio of the binary systems for a given EOS in a unique manner. For the soft
EOS (SFHo), the total ejecta mass depends weakly on the mass ratio, but the
average of electron number per baryon ($Y_e$) and specific entropy ($s$) of the
ejecta decreases significantly with the increase of the degree of mass
asymmetry. For the stiff EOS (DD2), with the increase of the mass asymmetry
degree, the total ejecta mass significantly increases while the average of
$Y_e$ and $s$ moderately decreases. We find again that only for the soft EOS
(SFHo), the total ejecta mass exceeds $0.01M_\odot$ irrespective of the mass
ratio chosen in this paper. The ejecta have a variety of electron number per
baryon with its average approximately between $Y_e \sim 0.2$ and $\sim 0.3$
irrespective of the EOS employed, which is well-suited for the production of
the r-process heavy elements (second and third peaks), although its averaged
value decreases with the increase of the degree of mass asymmetry.
[29]
oai:arXiv.org:1601.03398 [pdf] - 1359218
Alpha-constrained QSE Nucleosynthesis in High-entropy and Fast-expanding
Material
Submitted: 2016-01-13
We investigate the nucleosynthesis process in high-entropy ($s/k_{\rm
B}\gtrsim100$) and very fast-expanding ($\tau_{\rm exp}\sim10^{-3}\ {\rm s}$)
materials. In such a material with the electron fraction near 0.5, an
interesting nucleosynthesis process occurs. In this process, the abundance
distribution of heavy-nuclei of $A>100$ achieve quasi-statistical equilibrium
(QSE) at high temperature and the abundances are frozen at the end of the
nucleosynthesis. We explain this abundance distribution using the
"alpha-constrained QSE" abundances formulated in this paper. We demonstrate
that this nucleosynthesis would occur in neutrino-driven winds from massive
proto-neutron stars in hypernovae, where $A\sim140$ $p$-nuclei are synthesized.
[30]
oai:arXiv.org:1509.09205 [pdf] - 1331032
Efficient magnetic-field amplification due to the Kelvin-Helmholtz
instability in binary neutron star mergers
Submitted: 2015-09-30, last modified: 2015-11-13
We explore magnetic-field amplification due to the Kelvin-Helmholtz
instability during binary neutron star mergers. By performing high-resolution
general relativistic magnetohydrodynamics simulations with a resolution of
$17.5$ m for $4$--$5$ ms after the onset of the merger on the Japanese
supercomputer "K", we find that an initial magnetic field of moderate maximum
strength $10^{13}$ G is amplified at least by a factor of $\approx 10^3$. We
also explore the saturation of the magnetic-field energy and our result shows
that it is likely to be $\gtrsim 4 \times 10^{50}$ erg, which is $\gtrsim
0.1\%$ of the bulk kinetic energy of the merging binary neutron stars.
[31]
oai:arXiv.org:1507.05945 [pdf] - 1280930
Nucleosynthesis in neutrino-driven winds in hypernovae
Submitted: 2015-07-21, last modified: 2015-09-13
We investigate the nucleosynthesis in the neutrino-driven winds blown off
from a $3M_\odot$ massive proto-neutron star (mPNS) temporarily formed during
the collapse of a $100M_\odot$ star. Such mPNSs would be formed in hypernovae.
We construct steady and spherically symmetric wind solutions. We set large
neutrino luminosities of $\sim 10^{53}\ {\rm erg\ s^{-1}}$ and average energies
of electron neutrinos and antineutrinos in the ranges of
$\epsilon_{\nu_e}=9-16\ {\rm MeV}$ and $\epsilon_{\bar{\nu}_e}=11-18\ {\rm
MeV}$ based on a recent numerical relativity simulation. The wind solutions
indicate much shorter temperature-decrease timescale than that of the winds
from ordinary PNSs and, depending on $\epsilon_\nu$, the winds can be both
neutron-rich and proton-rich. In the neutron-rich wind, the $r$-process occurs
and the abundance distribution of a fiducial wind model of the mPNS gives an
approximate agreement with the abundance pattern of metal-poor weak $r$ star
HD~122563, although the third-peak elements are produced only when the
$\bar{\nu}_e$ energy is much larger than the $\nu_e$ energy. In the proton-rich
wind, the strong $\nu p$-process occurs and $A>100$ nuclides are synthesized.
The synthesized nuclei can be neutron-rich in some cases because the large
neutrino luminosity of the mPNS supplies a sufficient amount of neutrons.
[32]
oai:arXiv.org:1506.06811 [pdf] - 1284855
High-resolution magnetohydrodynamics simulation of black hole-neutron
star merger: Mass ejection and short gamma-ray burst
Submitted: 2015-06-22, last modified: 2015-09-03
We report results of a high-resolution numerical-relativity simulation for
the merger of black hole-magnetized neutron star binaries on Japanese
supercomputer "K". We focus on a binary that is subject to tidal disruption and
subsequent formation of a massive accretion torus. We find the launch of
thermally driven torus wind, subsequent formation of a funnel wall above the
torus and a magnetosphere with collimated poloidal magnetic field, and high
Blandford-Znajek luminosity. We show for the first time this picture in a
self-consistent simulation. The turbulence-like motion induced by the
non-axisymmetric magnetorotational instability as well as the Kelvin-Helmholtz
instability inside the accretion torus works as an agent to drive the mass
accretion and converts the accretion energy to thermal energy, which results in
the generation of a strong wind. By an in-depth resolution study, we reveal
that high resolution is essential to draw such a picture. We also discuss the
implication for the r-process nucleosynthesis, the radioactively-powered
transient emission, and short gamma-ray bursts.
[33]
oai:arXiv.org:1502.06660 [pdf] - 1231849
The dynamical mass ejection from binary neutron star mergers:
Radiation-hydrodynamics study in general relativity
Submitted: 2015-02-23, last modified: 2015-03-03
We perform radiation-hydrodynamics simulations of binary neutron star mergers
in numerical relativity on the Japanese "K" supercomputer, taking into account
neutrino cooling and heating by an updated leakage-plus-transfer scheme for the
first time. Neutron stars are modeled by three modern finite-temperature
equations of state (EOS) developed by Hempel and his collaborators. We find
that the properties of the dynamical ejecta of the merger such as total mass,
average electron fraction, and thermal energy depend strongly on the EOS. Only
for a soft EOS (the so-called SFHo), the ejecta mass exceeds $0.01M_{\odot}$.
In this case, the distribution of the electron fraction of the ejecta becomes
broad due to the shock heating during the merger. These properties are
well-suited for the production of the solar-like $r$-process abundance. For the
other stiff EOS (DD2 and TM1), for which a long-lived massive neutron star is
formed after the merger, the ejecta mass is smaller than $0.01M_{\odot}$,
although broad electron-fraction distributions are achieved by the positron
capture and the neutrino heating.
[34]
oai:arXiv.org:1407.2660 [pdf] - 1215573
High resolution numerical-relativity simulations for the merger of
binary magnetized neutron stars
Submitted: 2014-07-09, last modified: 2014-07-14
We perform high-resolution magnetohydrodynamics simulations of binary neutron
star mergers in numerical relativity on the Japanese supercomputer K. The
neutron stars and merger remnants are covered by a grid spacing of 70\,m, which
yields the highest-resolution results among those derived so far. By an
in-depth resolution study, we clarify several amplification mechanisms of
magnetic fields during the binary neutron star merger for the first time.
First, the Kelvin-Helmholtz instability developed in the shear layer at the
onset of the merger significantly amplifies the magnetic fields. A hypermassive
neutron star (HMNS) formed after the merger is then subject to the
nonaxisymmetric magnetorotational instability, which amplifies the magnetic
field in the HMNS. These two amplification mechanisms cannot be found with
insufficient-resolution runs. We also show that the HMNS eventually collapses
to a black hole surrounded by an accretion torus which is strongly magnetized
at birth.
[35]
oai:arXiv.org:1402.7317 [pdf] - 1203641
Production of all the r-process nuclides in the dynamical ejecta of
neutron star mergers
Submitted: 2014-02-28, last modified: 2014-06-10
Recent studies suggest that binary neutron star (NS-NS) mergers robustly
produce the heavy r-process nuclei above the atomic mass number A ~ 130 because
of their ejecta consisting of almost pure neutrons (electron fraction of Y_e <
0.1). However, little production of the lighter r-process nuclei (A = 90-120)
conflicts with the spectroscopic results of r-process-enhanced Galactic halo
stars. We present, for the first time, the result of nucleosynthesis
calculations based on the fully general-relativistic simulation of a NS-NS
merger with approximate neutrino transport. It is found that the bulk of the
dynamical ejecta are appreciably shock-heated and neutrino-processed, resulting
in a wide range of Y_e (= 0.09-0.45). The mass-averaged abundance distribution
of calculated nucleosynthesis yields is in reasonable agreement with the
full-mass range (A = 90-240) of the solar r-process curve. This implies, if our
model is representative of such events, that the dynamical ejecta of NS-NS
mergers can be the origin of the Galactic r-process nuclei. Our result also
shows that the radioactive heating after ~ 1 day from the merging, giving rise
to r-process-powered transient emission, is dominated by the beta-decays of
several species close to stability with precisely measured half-lives. This
implies that the total radioactive heating rate for such an event can be well
constrained within about a factor of two if the ejected material has a
solar-like r-process pattern.
[36]
oai:arXiv.org:1403.0956 [pdf] - 1208173
Jet collimation in the ejecta of double neutron star merger: new
canonical picture of short gamma-ray bursts
Submitted: 2014-03-04
The observations of jet breaks in the afterglows of short gamma-ray bursts
(SGRBs) indicate that the jet has a small opening angle of < 10{\deg}. The
collimation mechanism of the jet is a longstanding theoretical problem. We
numerically analyze the jet propagation in the material ejected by double
neutron star merger, and demonstrate that if the ejecta mass is > 10^{-2}
M_{sun}, the jet is well confined by the cocoon and emerges from the ejecta
with the required collimation angle. Our results also suggest that there are
some populations of choked (failed) SGRBs or low-luminous new types of event.
By constructing a model for SGRB 130603B, which is associated with the first
kilonova/macronova can- didate, we infer that the equation-of-state of neutron
stars would be soft enough to provide sufficient ejecta to collimate the jet,
if this event was associated with a double neutron star merger.
[37]
oai:arXiv.org:1310.2774 [pdf] - 1179857
Radioactively Powered Emission from Black Hole-Neutron Star Mergers
Submitted: 2013-10-10
Detection of electromagnetic counterparts of gravitational wave (GW) sources
is important to unveil the nature of compact binary coalescences. We perform
three-dimensional, time-dependent, multi-frequency radiative transfer
simulations for radioactively powered emission from the ejecta of black hole
(BH) - neutron star (NS) mergers. Depending on the BH to NS mass ratio, spin of
the BH, and equations of state of dense matter, BH-NS mergers can eject more
material than NS-NS mergers. In such cases, radioactively powered emission from
the BH-NS merger ejecta can be more luminous than that from NS-NS mergers. We
show that, in spite of the expected larger distances to BH-NS merger events,
observed brightness of BH-NS mergers can be comparable to or even higher than
that of NS-NS mergers. We find that, when the tidally disrupted BH-NS merger
ejecta are confined to a small solid angle, the emission from BH-NS merger
ejecta tends to be bluer than that from NS-NS merger ejecta for a given total
luminosity. Thanks to this property, we might be able to distinguish BH-NS
merger events from NS-NS merger events by multi-band observations of the
radioactively powered emission. In addition to the GW observations, such
electromagnetic observations can potentially provide independent information on
the nature of compact binary coalescences.
[38]
oai:arXiv.org:1310.1623 [pdf] - 1179743
Progenitor models of the electromagnetic transient associated with the
short GRB 130603B
Submitted: 2013-10-06
An electromagnetic transient powered by the radioactive decay of r-process
elements, a so-called kilonova/macronova, is one of the possible observable
consequences of compact binary mergers including at east one neutron star.
Recent observations strongly suggest the first discovery of the electromagnetic
transient, which is associated with the short GRB 130603B. We explore a
possible progenitor of this event combining the numerical-relativity
simulations and radiative transfer simulations of the dynamical ejecta of
binary neutron star and black hole - neutron star mergers. We show that the
ejecta models within a realistic parameter range consistently reproduce the
observed near-infrared excess. We also show that the soft equation of state
models for binary neutron star mergers and the stiff equation of state models
for black hole - neutron star mergers are favored to reproduce the observed
luminosity.
[39]
oai:arXiv.org:1307.5888 [pdf] - 1172900
Remnant massive neutron stars of binary neutron star mergers: Evolution
process and gravitational waveform
Submitted: 2013-07-22
Massive (hypermassive and supramassive) neutron stars are likely to be often
formed after the merger of binary neutron stars. We explore the evolution
process of the remnant massive neutron stars and gravitational waves emitted by
them, based on numerical-relativity simulations for binary neutron star mergers
employing a variety of equations of state and choosing a plausible range of the
neutron-star mass of binaries. We show that the lifetime of remnant
hypermassive neutron stars depends strongly on the total binary mass and also
on the equations of state. Gravitational waves emitted by the remnant massive
neutron stars universally have a quasiperiodic nature of an approximately
constant frequency although the frequency varies with time. We also show that
the frequency and time-variation feature of gravitational waves depend strongly
on the equations of state. We derive a fitting formula for the quasiperiodic
gravitational waveforms, which may be used for the data analysis of a
gravitational-wave signal.
[40]
oai:arXiv.org:1212.0905 [pdf] - 1158277
The mass ejection from the merger of binary neutron stars
Submitted: 2012-12-04, last modified: 2013-01-30
Numerical-relativity simulations for the merger of binary neutron stars are
performed for a variety of equations of state (EOSs) and for a plausible range
of the neutron-star mass, focusing primarily on the properties of the material
ejected from the system. We find that a fraction of the material is ejected as
a mildly relativistic and mildly anisotropic outflow with the typical and
maximum velocities $\sim 0.15$ -- $0.25c$ and $\sim 0.5$ -- $0.8c$ (where $c$
is the speed of light), respectively, and that the total ejected rest mass is
in a wide range $10^{-4}$ -- $10^{-2}M_{\odot}$, which depends strongly on the
EOS, the total mass, and the mass ratio. The total kinetic energy ejected is
also in a wide range between $10^{49}$ and $10^{51} {\rm ergs}$. The numerical
results suggest that for a binary of canonical total mass $2.7M_{\odot}$, the
outflow could generate an electromagnetic signal observable by the planned
telescopes through the production of heavy-element unstable nuclei via the
$r$-process or through the formation of blast waves during the interaction with
the interstellar matter, if the EOS and mass of the binary are favorable ones.
[41]
oai:arXiv.org:1212.4910 [pdf] - 1158634
Explicit-Implicit Scheme for Relativistic Radiation Hydrodynamics
Submitted: 2012-12-19
We propose an explicit-implicit scheme for numerically solving Special
Relativistic Radiation Hydrodynamic (RRHD) equations, which ensures a
conservation of total energy and momentum (matter and radiation). In our
scheme, 0th and 1st moment equations of the radiation transfer equation are
numerically solved without employing a flux-limited diffusion (FLD)
approximation. For an hyperbolic term, of which the time scale is the light
crossing time when the flow velocity is comparable to the speed of light, is
explicitly solved using an approximate Riemann solver. Source terms describing
an exchange of energy and momentum between the matter and the radiation via the
gas-radiation interaction are implicitly integrated using an iteration method.
The implicit scheme allows us to relax the Courant-Friedrichs-Lewy condition in
optically thick media, where heating/cooling and scattering timescales could be
much shorter than the dynamical timescale. We show that our numerical code can
pass test problems of one- and two-dimensional radiation energy transport, and
one-dimensional radiation hydrodynamics. Our newly developed scheme could be
useful for a number of relativistic astrophysical problems. We also discuss how
to extend our explicit-implicit scheme to the relativistic radiation
magnetohydrodynamics.
[42]
oai:arXiv.org:1206.5911 [pdf] - 1530093
Radiation Magnetohydrodynamics for Black Hole-Torus System in Full
General Relativity: A Step toward Physical Simulation
Submitted: 2012-06-26
A radiation-magnetohydrodynamic simulation for the black hole-torus system is
performed in the framework of full general relativity for the first time. A
truncated moment formalism is employed for a general relativistic neutrino
radiation transport. Several systems in which the black hole mass is $M_{\rm
BH}=3$ or $6M_{\odot}$, the black hole spin is zero, and the torus mass is
$\approx 0.14$--$0.38M_{\odot}$ are evolved as models of the remnant formed
after the merger of binary neutron stars or black hole-neutron star binaries.
The equation of state and microphysics for the high-density and
high-temperature matter are phenomenologically taken into account in a
semi-quantitative manner. It is found that the temperature in the inner region
of the torus reaches $\agt 10$ MeV which enhances a high luminosity of
neutrinos $\sim 10^{51}$ ergs/s for $M_{\rm BH}=6M_{\odot}$ and $\sim 10^{52}$
ergs/s for $M_{\rm BH}=3M_{\odot}$. It is shown that neutrinos are likely to be
emitted primarily toward the outward direction in the vicinity of the
rotational axis and their energy density may be high enough to launch a
low-energy short gamma-ray burst via the neutrino-antineutrino
pair-annihilation process with the total energy deposition $\sim
10^{47}$--$10^{49}$ ergs. It is also shown in our model that for $M_{\rm
BH}=3M_{\odot}$, the neutrino luminosity is larger than the electromagnetic
luminosity while for $M_{\rm BH}=6M_{\odot}$, the neutrino luminosity is
comparable to or slightly smaller than the electromagnetic luminosity.
[43]
oai:arXiv.org:1206.5927 [pdf] - 529443
Current Status of Numerical-Relativity Simulations in Kyoto
Submitted: 2012-06-26
We describe the current status of our numerical simulations for the collapse
of a massive stellar core to a BH and the BNS mergers, performed in the
framework of full general relativity incorporating finite-temperature EOS and
neutrino cooling. For the stellar core collapse simulation, we present the
latest numerical results. We employed a purely nucleonic EOS (Shen-EOS). As an
initial condition, we adopted a 100 $M_{\odot}$ presupernova model calculated
by Umeda and Nomoto. Changing the degree of rotation for the initial condition,
we clarify the strong dependence of the outcome of the collapse on this. When
the rotation is rapid enough, the shock wave formed at the core bounce is
deformed to be a torus-like shape. Then, the infalling matter is accumulated in
the central region due to the oblique shock at the torus surface, hitting the
PNS and dissipating the kinetic energy there. As a result, outflows can be
launched. The PNS eventually collapses to a BH and an accretion torus is formed
around it. We also found that the evolution of the BH and torus depends
strongly on the rotation initially given. In the BNS merger simulations, we in
addition employ an EOS incorporating a degree of freedom for hyperons. The
numerical simulations show that for the purely nucleonic EOS, a HMNS with a
long lifetime ($\gg 10$ ms) is the outcome for the total mass $M \lesssim
3.0M_{\odot}$. By contrast, the formed HMNS collapses to a BH in a shorter time
scale with the hyperonic EOS for $M \gtrsim 2.7M_{\odot}$. It is shown that the
typical total neutrino luminosity of the HMNS is $\sim 3$--$10\times 10^{53}$
ergs/s and the effective amplitude of gravitational waves from the HMNS is
2--$6 \times 10^{-22}$ at $f\approx 2$--2.5 kHz for a source distance of 100
Mpc.
[44]
oai:arXiv.org:1206.0509 [pdf] - 1123858
Gravitational waves, neutrino emissions, and effects of hyperons in
binary neutron star mergers
Submitted: 2012-06-03
Numerical simulations for the merger of binary neutron stars are performed in
full general relativity incorporating both nucleonic and hyperonic
finite-temperature equations of state (EOS) and neutrino cooling. It is found
that for the nucleonic and hyperonic EOS, a hyper massive neutron star (HMNS)
with a long lifetime $(t_{\rm life}\gtrsim 10 {\rm ms})$ is the outcome for the
total mass $\approx 2.7 M_\odot$. For the total mass $\approx 3 M_\odot$, a
long-lived (short-lived with $t_{\rm life}\approx 3 {\rm ms}$) HMNS is the
outcome for the nucleonic (hyperonic) EOS. It is shown that the typical total
neutrino luminosity of the HMNS is $\sim 3$ -- $6 \times 10^{53} {\rm erg /s}$
and the effective amplitude of gravitational waves from the HMNS is 1 --
$4\times 10^{-22}$ at $f\approx 2$ -- $3.2 {\rm kHz}$ for a source of distance
of 100 Mpc. During the HMNS phase, characteristic frequencies of gravitational
waves shift to a higher frequency for the hyperonic EOS in contrast to the
nucleonic EOS in which they remain constant approximately. Our finding suggests
that the effects of hyperons are well imprinted in gravitational wave and its
detection will give us a potential opportunity to explore the composition of
the neutron star matter. We present the neutrino luminosity curve when a black
hole is formed as well.
[45]
oai:arXiv.org:1110.4442 [pdf] - 1085007
Effects of hyperons in binary neutron star mergers
Submitted: 2011-10-20
Numerical simulations for the merger of binary neutron stars are performed in
full general relativity incorporating both nucleonic and hyperonic
finite-temperature equations of state (EOS) and neutrino cooling for the first
time. It is found that even for the hyperonic EOS, a hypermassive neutron star
is first formed after the merger for the typical total mass $\approx$
2.7M\bigodot, and subsequently collapses to a black hole (BH). It is shown that
hyperons play a substantial role in the post-merger dynamics, torus formation
around the BH, and emission of gravitational waves (GWs). In particular, the
existence of hyperons is imprinted in GWs. Therefore, GW observations will
provide a potential opportunity to explore the composition of the neutron star
matter.
[46]
oai:arXiv.org:1105.2125 [pdf] - 400352
Gravitational waves and neutrino emission from the merger of binary
neutron stars
Submitted: 2011-05-11, last modified: 2011-06-23
Numerical simulations for the merger of binary neutron stars are performed in
full general relativity incorporating a finite-temperature (Shen's) equation of
state (EOS) and neutrino cooling for the first time. It is found that for this
stiff EOS, a hypermassive neutron star (HMNS) with a long lifetime ($\gg 10$
ms) is the outcome for the total mass $\alt 3.0M_{\odot}$. It is shown that the
typical total neutrino luminosity of the HMNS is $\sim 3$--$8\times 10^{53}$
ergs/s and the effective amplitude of gravitational waves from the HMNS is
4--$6 \times 10^{-22}$ at $f=2.1$--2.5 kHz for a source distance of 100 Mpc. We
also present the neutrino luminosity curve when a black hole is formed for the
first time.
[47]
oai:arXiv.org:1009.5303 [pdf] - 1040863
Formation of black hole and accretion disk in a massive high-entropy
stellar core collapse
Submitted: 2010-09-27, last modified: 2011-05-24
We present the first numerical result of fully general relativistic
axisymmetric simulations for the collapse of a rotating high-entropy stellar
core to a black hole and an accretion disk. The simulations are performed
taking into account the relevant microphysics. We adopt as initial condition a
spherical core with constant electron fraction ($Y_e = 0.5$) and entropy per
baryon $s$ = 8 $k_B$, and angular velocity is superimposed. In the early phase,
the core collapses in a homologous manner. Then, it experiences a weak bounce
due to the gas pressure of free nucleons. Because the bounce is weak, the core
collapses eventually to a black hole. Subsequent evolution depends on initial
angular velocity. When the rotation is not fast, a geometrically thin (but
optically thick) accretion disk is formed, and shock waves are formed in the
inner part of the disk. For the moderately rotating case, the thin accretion
disk expands eventually to be a geometrically thick torus after sufficient
accumulation of the thermal energy generated at the shocks. Furthermore,
convection occurs inside the torus. Neutrino luminosities vary violently with
time because of the convective motion. For the rapidly rotating case, by
contrast, a geometrically thick torus is formed soon after the black hole
formation, and convective activity is weak due to the presence of epicyclic
mode.
[48]
oai:arXiv.org:1104.3937 [pdf] - 379763
Truncated Moment Formalism for Radiation Hydrodynamics in Numerical
Relativity
Submitted: 2011-04-20
A truncated moment formalism for general relativistic radiation
hydrodynamics, based on the Thorne's moment formalism, is derived. The fluid
rest frame is chosen to be the fiducial frame for defining the radiation
moments. Then, zeroth-, first-, and second-rank radiation moments are defined
from the distribution function with a physically reasonable assumption for it
in the optically thin and thick limits. The source terms are written, focusing
specifically on the neutrino transfer and neglecting higher harmonic angular
dependence of the reaction angle. Finally, basic equations for a truncated
moment formalism for general relativistic radiation hydrodynamics in a closed
covariant form are derived assuming a closure relation among the radiation
stress tensor, energy density, and energy flux, and a variable Eddington
factor, which works well.
[49]
oai:arXiv.org:1009.3358 [pdf] - 1040653
An implementation of the microphysics in full general relativity :
General relativistic neutrino leakage scheme
Submitted: 2010-09-17
Performing fully general relativistic simulations taking account of
microphysical processes is one of long standing problems in numerical
relativity. One of main difficulties in implementation of weak interactions in
the general relativistic framework lies on the fact that the characteristic
timescale of weak interaction processes (the WP timescale) in hot dense matters
is much shorter than the dynamical timescale. Numerically this means that stiff
source terms appears in the equations so that an implicit scheme is in general
necessary to stably solve the relevant equations. Otherwise a very short
timestep will be required to solve them explicitly. Furthermore, in the
relativistic framework, the Lorentz factor is coupled with the rest mass
density and the energy density. The specific enthalpy is also coupled with the
momentum. Due to these couplings, it is very complicated to recover the
primitive variables and the Lorentz factor from conserved quantities. At the
current status, no implicit procedure have been proposed except for the case of
the spherical symmetry. Therefore, an approximate, explicit procedure is
developed in the fully general relativistic framework in this paper as an first
implementation of the microphysics toward a more realistic sophisticated model.
The procedure is based on the so-called neutrino leakage schemes which is based
on the property that the characteristic timescale in which neutrinos leak out
of the system (the leakage timescale) is much longer than the WP timescale. In
this paper, I present a detailed neutrino leakage scheme and a simple and
stable method for solving the equations explicitly in the fully general
relativistic framework. I also perform a test simulation to check the validity
of the present method, showing that it works fairly well.
[50]
oai:arXiv.org:1009.3320 [pdf] - 275616
Stellar core collapse in full general relativity with microphysics -
Formulation and Spherical collapse test -
Submitted: 2010-09-16
One of the longstanding issues in numerical relativity is to enable a
simulation taking account of microphysical processes (e.g., weak interactions
and neutrino cooling). We develop an approximate and explicit scheme in the
fully general relativistic framework as a first implementation of the
microphysics toward a more realistic and sophisticated modeling. In this paper,
we describe in detail a method for implementation of a realistic equation of
state, the electron capture and the neutrino cooling in a multidimensional,
fully general relativistic code. The procedure is based on the so-called
neutrino leakage scheme. To check the validity of the code, we perform a two
dimensional (2D) simulation of spherical stellar core collapse. Until the
convective activities set in, our results approximately agree, or at least are
consistent, with those in the previous so-called state-of-the-art simulations.
In particular, the radial profiles of thermodynamical quantities and the time
evolution of the neutrino luminosities agree quantitatively. The convection is
driven by negative gradients of the entropy per baryon and the electron
fraction as in the previous 2D Newtonian simulations. We clarify which gradient
is more responsible for the convection. Gravitational waves from the convection
are also calculated. We find that the characteristic frequencies of the
gravitational-wave spectra are distributed for higher frequencies than those in
Newtonian simulations due to the general relativistic effects.
[51]
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.
[52]
oai:arXiv.org:1002.2689 [pdf] - 146883
Exploring binary-neutron-star-merger scenario of short-gamma-ray bursts
by gravitational-wave observation
Submitted: 2010-02-13
We elucidate the feature of gravitational waves (GWs) from binary neutron
star merger collapsing to a black hole by general relativistic simulation. We
show that GW spectrum imprints the coalescence dynamics, formation process of
disk, equation of state for neutron stars, total masses, and mass ratio. A
formation mechanism of the central engine of short $\gamma$-ray bursts, which
are likely to be composed of a black hole and surrounding disk, therefore could
be constrained by GW observation.
[53]
oai:arXiv.org:0904.4551 [pdf] - 23770
Longterm general relativistic simulation of binary neutron stars
collapsing to a black hole
Submitted: 2009-04-29, last modified: 2009-08-31
General relativistic simulations for the merger of binary neutron stars are
performed as an extension of a previous work\cite{Shibata:2006nm}. We prepare
binary neutron stars with a large initial orbital separation and employ the
moving-puncture formulation, which enables to follow merger and ringdown phases
for a long time, even after black hole formation. For modeling inspiraling
neutron stars, which should be composed of cold neutron stars, the
Akmal-Pandhalipande-Ravenhall (APR) equation of state (EOS) is adopted. After
the onset of the merger, the hybrid-type EOS is used; i.e., the cold and
thermal parts are given by the APR and $\Gamma$-law EOSs, respectively. Three
equal-mass binaries each with mass $1.4M_\odot,1.45M_\odot,1.5M_\odot$ and two
unequal-mass binaries with mass 1.3--$1.6M_\odot$, 1.35--$1.65M_\odot$ are
prepared. We focus primarily on the black hole formation case, and explore mass
and spin of the black hole, mass of disks which surround the black hole, and
gravitational waves emitted during the black hole formation. We find that (i)
for the systems of $m_0=2.9$--$3.0M_\odot$ and of mass ratio $\approx 0.8$, the
mass of disks which surround the formed black hole is 0.006--$0.02M_{\odot}$;
(ii) the spin of the formed black hole is $0.78 \pm 0.02$ when a black hole is
formed after the merger in the dynamical time scale. This value depends weakly
on the total mass and mass ratio, and is about 0.1 larger than that of a black
hole formed from nonspinning binary black holes; (iii) for the black-hole
formation case, Fourier spectrum shape of gravitational waves emitted in the
merger and ringdown phases has a universal qualitative feature irrespective of
the total mass and mass ratio, but quantitatively, the spectrum reflects the
parameters of the binary neutron stars.
[54]
oai:arXiv.org:0709.1766 [pdf] - 4847
Magnetohydrodynamics of Neutrino-Cooled Accretion Tori around a Rotating
Black Hole in General Relativity
Submitted: 2007-09-12
We present our first numerical results of axisymmetric magnetohydrodynamic
simulations for neutrino-cooled accretion tori around rotating black holes in
general relativity. We consider tori of mass $\sim 0.1$--0.4$M_{\odot}$ around
a black hole of mass $M=4M_{\odot}$ and spin $a=0$--$0.9M$; such systems are
candidates for the central engines of gamma-ray bursts (GRBs) formed after the
collapse of massive rotating stellar cores and the merger of a black hole and a
neutron star. In this paper, we consider the short-term evolution of a torus
for a duration of $\approx 60$ ms, focusing on short-hard GRBs. Simulations
were performed with a plausible microphysical equation of state that takes into
account neutronization, the nuclear statistical equilibrium of a gas of free
nucleons and $\alpha$-particles, black body radiation, and a relativistic Fermi
gas (neutrinos, electrons, and positrons). Neutrino-emission processes, such as
$e^{\pm}$ capture onto free nucleons, $e^{\pm}$ pair annihilation, plasmon
decay, and nucleon-nucleon bremsstrahlung are taken into account as cooling
processes. Magnetic braking and the magnetorotational instability in the
accretion tori play a role in angular momentum redistribution, which causes
turbulent motion, resultant shock heating, and mass accretion onto the black
hole. The mass accretion rate is found to be $\dot M_* \sim 1$--$10
M_{\odot}$/s, and the shock heating increases the temperature to $\sim 10^{11}$
K. This results in a maximum neutrino emission rate of $L_{\nu}=$ several
$\times 10^{53}$ ergs/s and a conversion efficiency $L_{\nu}/\dot M_* c^2$ on
the order of a few percent for tori with mass $M_{\rm t} \approx
0.1$--0.4$M_{\odot}$ and for moderately high black hole spins.
[55]
oai:arXiv.org:0706.4154 [pdf] - 2647
A Formation Mechanism of Collapsar Black Hole -- early evolution phase
Submitted: 2007-06-28
The latest studies of massive star evolution indicate that an initially
rapidly rotating star with sufficiently low metallicity can produce a rapidly
rotating, massive stellar core that could be a progenitor of long-soft
gamma-ray bursts (LGRBs). Motivated by these studies, we follow the collapse of
a rapidly rotating massive stellar core to a 'collapsar' black hole (BH)
surrounded by a massive, hot accretion disk performing fully general
relativistic simulations. We focus on the general relativistic dynamics of the
collapse, and the relevant microphysics is treated in a qualitative manner. The
simulations are performed until the system consisting of the BH and the disk
has relaxed to a quasi-stationary state. A novel mechanism found in this study
is that strong shock waves are formed at the inner part of the disk after the
formation of the BH. These shock waves propagate mainly along the rotation
axis, heating the disk and sweeping materials around the rotational axis, and
thereby forming a low density region. The temperature of the disk is high
enough for copious neutrino emission. All these features indicate that the
direct formation of a rapidly rotating BH is a promising source of LGRBs even
in the absence of strong magnetic fields.
[56]
oai:arXiv.org:astro-ph/0507383 [pdf] - 74552
Magnetohydrodynamics in full general relativity: Formulation and tests
Submitted: 2005-07-16
A new implementation for magnetohydrodynamics (MHD) simulations in full
general relativity (involving dynamical spacetimes) is presented. In our
implementation, Einstein's evolution equations are evolved by a BSSN formalism,
MHD equations by a high-resolution central scheme, and induction equation by a
constraint transport method. We perform numerical simulations for standard test
problems in relativistic MHD, including special relativistic magnetized shocks,
general relativistic magnetized Bondi flow in stationary spacetime, and a
longterm evolution for self-gravitating system composed of a neutron star and a
magnetized disk in full general relativity. In the final test, we illustrate
that our implementation can follow winding-up of the magnetic field lines of
magnetized and differentially rotating accretion disks around a compact object
until saturation, after which magnetically driven wind and angular momentum
transport inside the disk turn on.
[57]
oai:arXiv.org:astro-ph/0504567 [pdf] - 72664
Axisymmetric collapse simulations of rotating massive stellar cores in
full general relativity: Numerical study for prompt black hole formation
Submitted: 2005-04-26
We perform axisymmetric simulations for gravitational collapse of a massive
iron core to a black hole in full general relativity. The iron cores are
modeled by $\Gamma = 4/3$ equilibrium polytrope for simplicity and a parametric
equation of state is adopted during the collapse. Performing simulations for a
wide variety of initial conditions changing the mass, the angular momentum, the
rotational velocity profile of the core, and the parameters of the equations of
state (which are chosen so that the maximum mass of the cold spherical
polytrope is $\approx 1.6M_{\odot}$), a criterion for the prompt black hole
formation is clarified in terms of them. It is found that (i) with the increase
of the thermal energy generated by shocks, the threshold mass for the prompt
black hole formation is increased by 20--40%, (ii) the rotational centrifugal
force increases the threshold mass by $\alt 25%$, (iii) with the increase of
the degree of differential rotation, the threshold mass is also increased, and
(iv) the amplification factors shown in the results (i)--(iii) depend
sensitively on the equation of state. We also find that the collapse dynamics
and the structure of the shock formed at the bounce depend strongly on the
stiffness of the adopted equation of state. In particular, as a new feature, a
strong bipolar explosion is observed for the collapse of rapidly rotating iron
cores with an equation of state which is stiff in subnuclear density and soft
in supranuclear density. Gravitational waves are computed in terms of a
quadrupole formula. It is also found that the waveform depends sensitively on
the equations of state.
[58]
oai:arXiv.org:astro-ph/0412243 [pdf] - 69647
Three-dimensional simulations of stellar core collapse in full general
relativity: Nonaxisymmetric dynamical instabilities
Submitted: 2004-12-10, last modified: 2005-01-11
We perform fully general relativistic simulations of rotating stellar core
collapse in three spatial dimension. A parametric equation of state is adopted
following Dimmelmeier et al. The early stage of the collapse is followed by an
axisymmetric code. When the stellar core becomes compact enough, we start a
3-dimensional simulation adding a bar-mode nonaxisymmetric density
perturbation. In the axisymmetric simulations, it is clarified that the maximum
value of $\beta \equiv T/W$ achieved during the stellar collapse and depends
sensitively on the velocity profile and total mass of the initial core, and
equations of state. It is also found that for all the models with high degree
of differential rotation, a funnel structure is formed around the rotational
axis after the formation of neutron stars. For selected models in which the
maximum value of $\beta$ is larger than $\sim 0.27$, 3-dimensional simulations
are performed. It is found that the bar-mode dynamical instability sets in for
the case that the following conditions are satisfied: (i) the progenitor should
be rapidly rotating with the initial value of $0.01 \alt \beta \alt 0.02$, (ii)
the degree of differential rotation of the initial condition should be
sufficiently high, and (iii) a depletion factor of pressure in an early stage
of collapse should be large enough to induce a significant contraction for
which an efficient spin-up can be achieved. As a result of the onset of the
bar-mode dynamical instabilities, the amplitude of gravitational waves can be
by a factor of $\sim 10$ larger than that in the axisymmetric collapse. It is
found that a dynamical instability with the $m=1$ mode is also induced, but the
perturbation does not grow significantly and, hence, it does not contribute to
an outstanding amplification of gravitational waves.
[59]
oai:arXiv.org:gr-qc/0403036 [pdf] - 111471
New criterion for direct black hole formation in rapidly rotating
stellar collapse
Submitted: 2004-03-08
We study gravitational collapse of rapidly rotating relativistic polytropes
of the adiabatic index $\Gamma = 1.5$ and 2, in which the spin parameter $q
\equiv J/M^{2} > 1$ where $J$ and $M$ are total angular momentum and
gravitational mass, in full general relativity.
First, analyzing initial distributions of the mass and the spin parameter
inside stars, we predict the final outcome after the collapse. Then, we perform
fully general relativistic simulations on assumption of axial and equatorial
symmetries and confirm our predictions. As a result of simulations, we find
that in contrast with the previous belief, even for stars with $q > 1$, the
collapse proceeds to form a seed black hole at central region, and the seed
black hole subsequently grows as the ambient fluids accrete onto it. We also
find that growth of angular momentum and mass of the seed black hole can be
approximately determined from the initial profiles of the density and the
specific angular momentum. We define an effective spin parameter at the central
region of the stars, $q_{c}$, and propose a new criterion for black hole
formation as $q_{c} \alt 1$. Plausible reasons for the discrepancy between our
and previous results are clarified.
[60]
oai:arXiv.org:astro-ph/0402184 [pdf] - 62698
Gravitational waves from axisymmetrically oscillating neutron stars in
general relativistic simulations
Submitted: 2004-02-08
Gravitational waves from oscillating neutron stars in axial symmetry are
studied performing numerical simulations in full general relativity. Neutron
stars are modeled by a polytropic equation of state for simplicity. A
gauge-invariant wave extraction method as well as a quadrupole formula are
adopted for computation of gravitational waves. It is found that the
gauge-invariant variables systematically contain numerical errors generated
near the outer boundaries in the present axisymmetric computation. We clarify
their origin, and illustrate it possible to eliminate the dominant part of the
systematic errors. The best corrected waveforms for oscillating and rotating
stars currently contain errors of magnitude $\sim 10^{-3}$ in the local wave
zone. Comparing the waveforms obtained by the gauge-invariant technique with
those by the quadrupole formula, it is shown that the quadrupole formula yields
approximate gravitational waveforms besides a systematic underestimation of the
amplitude of $O(M/R)$ where $M$ and $R$ denote the mass and the radius of
neutron stars. However, the wave phase and modulation of the amplitude can be
computed accurately. This indicates that the quadrupole formula is a useful
tool for studying gravitational waves from rotating stellar core collapse to a
neutron star in fully general relativistic simulations. Properties of the
gravitational waveforms from the oscillating and rigidly rotating neutron stars
are also addressed paying attention to the oscillation associated with
fundamental modes.
[61]
oai:arXiv.org:gr-qc/0402040 [pdf] - 386270
Gravitational waves from axisymmetric rotating stellar core collapse to
a neutron star in full general relativity
Submitted: 2004-02-08
Axisymmetric numerical simulations of rotating stellar core collapse to a
neutron star are performed in the framework of full general relativity. The
so-called Cartoon method, in which the Einstein field equations are solved in
the Cartesian coordinates and the axisymmetric condition is imposed around the
$y=0$ plane, is adopted. The hydrodynamic equations are solved in the
cylindrical coordinates (on the $y=0$ plane in the Cartesian coordinates) using
a high-resolution shock-capturing scheme with the maximum grid size
$(2500,2500)$. A parametric equation of state is adopted to model collapsing
stellar cores and neutron stars following Dimmelmeier et al. It is found that
the evolution of central density during the collapse, bounce, and formation of
protoneutron stars agree well with those in the work of Dimmelmeier et al. in
which an approximate general relativistic formulation is adopted. This
indicates that such approximation is appropriate for following axisymmetric
stellar core collapses and subsequent formation of protoneutron stars.
Gravitational waves are computed using a quadrupole formula. It is found that
the waveforms are qualitatively in good agreement with those by Dimmelmeier et
al. However, quantitatively, two waveforms do not agree well. Possible reasons
for the disagreement are discussed.
