scispace - formally typeset
Search or ask a question
Journal ArticleDOI

3D collapse of rotating stellar iron cores in general relativity including deleptonization and a nuclear equation of state.

29 Jun 2007-Physical Review Letters (American Physical Society)-Vol. 98, Iss: 26, pp 261101
TL;DR: This work focuses on gravitational wave (GW) emission from rotating collapse, bounce, and early postbounce phases and indicates that the GW signature of these phases is much more generic than previously estimated.
Abstract: We present 2D and 3D simulations of the collapse of rotating stellar iron cores in general relativity employing a nuclear equation of state and an approximate treatment of deleptonization. We compare fully general relativistic and conformally flat evolutions and find that the latter treatment is sufficiently accurate for the core-collapse supernova problem. We focus on gravitational wave (GW) emission from rotating collapse, bounce, and early postbounce phases. Our results indicate that the GW signature of these phases is much more generic than previously estimated. We also track the growth of a nonaxisymmetric instability in one model, leading to strong narrow-band GW emission.

Content maybe subject to copyright    Report

3D Collapse of Rotating Stellar Iron Cores in General Relativity Including Deleptonization
and a Nuclear Equation of State
C. D. Ott,
1,5,
*
H. Dimmelmeier,
2
A. Marek,
2
H.-T. Janka,
2
I. Hawke,
3
B. Zink,
2,4
and E. Schnetter
4
1
Max-Planck-Institut fu
¨
r Gravitationsphysik, Albert-Einstein-Institut, Am Mu
¨
hlenberg 1, D-14476 Potsdam, Germany
2
Max-Planck-Institut fu
¨
r Astrophysik, Karl-Schwarzschild-Straße 1, D-85741 Garching, Germany
3
School of Mathematics, University of Southampton, Southampton SO17 1BJ, UK
4
Center for Computation and Technology, Louisiana State University, Baton Rouge, Louisiana 70803, USA
5
Steward Observatory and Department of Astronomy, University of Arizona, 933 N. Cherry Avenue, Tucson, AZ 85721, USA
(Received 29 September 2006; published 29 June 2007)
We present 2D and 3D simulations of the collapse of rotating stellar iron cores in general relativity
employing a nuclear equation of state and an approximate treatment of deleptonization. We compare fully
general relativistic and conformally flat evolutions and find that the latter treatment is sufficiently accurate
for the core-collapse supernova problem. We focus on gravitational wave (GW) emission from rotating
collapse, bounce, and early postbounce phases. Our results indicate that the GW signature of these phases
is much more generic than previously estimated. We also track the growth of a nonaxisymmetric
instability in one model, leading to strong narrow-band GW emission.
DOI: 10.1103/PhysRevLett.98.261101 PACS numbers: 04.25.Dm, 04.30.Db, 95.30.Sf, 97.60.Bw
Introduction.For more than two decades, astrophysi-
cists have struggled to compute the gravitational wave
(GW) signal produced by rotating stellar core collapse
and the subsequent supernova evolution. Besides the co-
alescence of black hole and neutron star binaries, core-
collapse events are considered to be among the most
promising sources of detectable GWs. Theoretical predic-
tions are still hampered by three major problems: (i) the
rotational configuration prior to gravitational collapse is
still uncertain since multi-D evolutionary calculations
of rotating massive stars have not yet been performed;
(ii) reliable waveform estimates require a general relativ-
istic (GR) treatment, since both high densities and veloc-
ities in combination with strong gravitational fields are
encountered in this problem; and (iii) an adequate treat-
ment of the nuclear equation of state (EOS) and the neu-
trino microphysics and radiative transfer is crucial for
obtaining realistic collapse, bounce, and postbounce dy-
namics and waveforms. GW emission from core-collapse
supernovae may arise from rotating collapse and bounce,
postbounce neutrino-driven convection, anisotropic neu-
trino emission, nonaxisymmetric rotational instabilities
of the protoneutron star (PNS), or from the recently pro-
posed PNS core g-mode oscillations. Previous estimates of
the GW signature of core-collapse supernovae have either
relied on Newtonian simulations [16] (to some extent
approximating GR effects [7]) or GR simulations with
simplified analytic (so-called hybrid) EOSs and no neu-
trino treatment [811]. Depending on rotation, softness of
the EOS at subnuclear densities, and inclusion of GR
effects, the collapse dynamics and accordingly the GW
signature can differ significantly.
Here, we present new results from GR simulations,
focusing on the rotating collapse, bounce, and early post-
bounce phases. These are the first-ever multi-D GR simu-
lations with presupernova models from stellar evolution
calculations, a microphysical nuclear EOS, and a simple
but effective treatment of electron capture and neutrino
radiation effects during collapse. In this way, we obtain the
most accurate estimates of the GW signature of rotating
stellar core collapse in full GR to date.
Method and initial models.We perform all 3D simu-
lations in 3 1 GR using
CACTUS [ 12], Cartesian coordi-
nates, and mesh refinement provided by
CARPET [13].
Spacetime is evolved using the BSSN formulation (see,
e.g., [14]) with 1 log slicing and a hyperbolic shift [15].
We use the finite-volume GR hydrodynamics code
WHISKY
[16]. Typical simulation grids extend to 3000 km and use 9
refinement levels. The central resolution is r 350 m.In
addition, we perform axisymmetric (2D) simulations for
all models using the
COCONUT code [8,17], approximating
GR by the conformal flatness condition (CFC).
COCONUT
utilizes spherical coordinates with 250 logarithmically-
spaced radial zones with r 250 m and 45 equidistant
angular zones. Resolution tests with both codes do not
yield significant qualitative or quantitative differences.
We extract GWs using a variant of the Newtonian quadru-
pole formula [9].
We employ the microphysical nuclear EOS of Shen et al.
[18] as implemented in [19]. Deleptonization is approxi-
mated as proposed by [20]: The electron fraction Y
e
is
parameterized as a function of density based on data from
spherically symmetric radiation-hydrodynamics calcula-
tions with standard electron capture rates [21]. After core
bounce, Y
e
is passively advected, and further lepton loss is
neglected, but neutrino pressure continues to be taken into
account above trapping density [20].
In this Letter, we focus on the collapse of massive
presupernova iron cores with at most moderate differential
rotation, and rotation rates that may be too fast to match
PRL 98, 261101 (2007)
PHYSICAL REVIEW LETTERS
week ending
29 JUNE 2007
0031-9007=07=98(26)=261101(4) 261101-1 © 2007 The American Physical Society

garden-variety pulsar birth spin estimates [22,23], but
could be relevant in the collapsar-type gamma-ray burst
context. As initial data, we use the nonrotating 20M
presupernova model s20 of Woosley et al. [24], which
we force to rotate according to the rotation law discussed
in [5,8]. We parameterize our models in terms of the
differential rotation parameter A and the initial ratio of
rotational kinetic to gravitational energy
i
T=jWj.In
addition, we perform a calculation with the 20M
model
E20A of [25], which includes an approximate 1D treatment
of rotation. In Table I, we list the model parameters.
Results.In Fig. 1, we compare GW signals computed
with
COCONUT in 2D-CFC and those computed with our
3D-full-GR approach. Model s20A2B2 (red lines) is a
moderate rotator with a
i
0:50%, rotating rigidly in
its central region. It stays axisymmetric throughout its
numerical evolution. The agreement of 2D-CFC with 3D-
full-GR is excellent for this model: Both waveforms match
almost perfectly at bounce and during the very early post-
bounce phase. A few ms after bounce, when convection in
the region behind the stalling shock sets in due to a
negative entropy gradient, the signals begin to differ quan-
titatively while remaining in phase. We attribute this small
mismatch to the choice of coordinate grids and to differ-
ences in the growth and scale of vortical postbounce mo-
tions between 2D and 3D. Model s20A1B5 rotates with
constant in the entire core. Despite its very large
i
4%, it remains essentially axisymmetric during the time
covered by our simulation, since most of its angular mo-
mentum is attached to material at large radii that falls
inward and spins up only slowly. The waveforms in CFC
and full GR agree very well. Again, both waveforms match
best for the strong burst related to core bounce during
which more than 90% of the total GW energy are emitted
in an axisymmetric model. The overall excellent agree-
ment of CFC with full GR confirms results of [9,11] and
proves that CFC is a very good approximation to full GR in
the core-collapse scenario.
In Fig. 2, we present waveforms of models with varying
initial degree of differential rotation A and rotation rate
i
.
We find that the inclusion of a microphysical EOS and of
electron capture yields results that differ considerably from
those obtained in previous, less sophisticated studies.
Figure 2 exemplifies that largely independent of the initial
rotational configuration, the GW signal of rotating collapse
has a generic shape: a slow rise in the prebounce phase, a
large negative amplitude at bounce when the motion of the
inner core is reversed, followed by a ring-down. This so-
called ‘Type I’ signature corresponds to a baryonic pres-
sure dominated bounce [1,2,5,8]. Thus, all our models
undergo core bounce dominated by the stiffening of the
EOS at nuclear density.
TABLE I. Summary.
b
is the density at bounce, the maximum characteristic GW strain
h
char;max
is at a distance of 10 kpc, and E
gw
is the energy emitted in GWs. Values for E20A
pb
include GW emission from late-time 3D dynamics.
Model A [10
8
cm]
i
[%]
b
[%]
b
10
14
g
cm
3
h
char;max
[10
21
] E
gw
[10
9
M
c
2
]
s20A1B1 50.0 0.25 0.90 3.29 1.46 0.6
s20A1B5 50.0 4.00 10.52 2.90 9.68 26.9
s20A2B2 1.0 0.50 6.72 3.07 8.77 22.0
s20A2B4 1.0 1.80 16.33 2.35 4.28 9.4
s20A3B3 0.5 0.90 16.57 2.33 4.58 12.4
E20A  0.37 11.31 2.79 12.18 36.9
E20A
pb
24.23 75.4
t-t
bounce
(ms)
h
+
(10
−21
at 10 kpc)
-5 0 5 10
-10
-8
-6
-4
-2
0
2
4
6
s20A2B2 3D full GR
s20A2B2 2D CFC
s20A1B5 3D full GR
s20A1B5 2D CFC
FIG. 1 (color). GW strain h
along the equator for models
s20A2B2 and s20A1B5. We compare 2D-CFC and 3D-full-GR
results.
-5 0 5 10 15
-12
-10
-8
-6
-4
-2
0
2
4
6
t-t
bounce
(ms)
h
+
(10
−21
at 10 kpc)
s20A2B1
s20A2B2
s20A2B4
s20A3B3
E20A
FIG. 2 (color). Equatorial GW strain h
for a representative
subset of the models listed in Table I. Note the generic shape of
the core bounce GW signal.
PRL 98, 261101 (2007)
PHYSICAL REVIEW LETTERS
week ending
29 JUNE 2007
261101-2

This is in stark contrast to the studies using the simple
hybrid EOS [2,810], where initial models with rotation
rates in the range investigated here develop sufficient
centrifugal support during contraction to stop collapse at
subnuclear densities, resulting in several consecutive cen-
trifugal bounces separated by phases of coherent reexpan-
sion of the inner core. While in GR, models exhibiting such
multiple centrifugal bounce and the corresponding ‘Type
II’ GW signals are only rarer compared to Newtonian
gravity [8], we do not observe any such model in this study.
An evident example is model s20A2B4: In previous studies
without detailed microphysics, the corresponding model
with identical initial rotation parameters (A2B4G1)
showed clear ‘Type II’ behavior in both Newtonian and
GR calculations [2,8].
The suppression of the multiple centrifugal bounce sce-
nario is due to two physical effects: (i) In contrast to the
simple hybrid EOS, in our case the mass and dynamics of
the inner core (which is most important for the GW emis-
sion) is not merely determined by the adiabatic index
d lnP=d ln (at constant entropy) of the EOS, but also by
deleptonization during collapse. This leads to considerably
smaller inner cores with less angular momentum and
weaker pressure support. (ii) Since multiple centrifugal
bounce was observed for a model with moderately fast
rotation in a previous Newtonian study [1] employing
both a microphysical EOS and deleptonization, the ab-
sence of this collapse type in our study is not only caused
by microphysical effects, but also by the effectively
stronger gravity in GR. This is in accordance with simu-
lations using the simple hybrid EOS [8]. For a more de-
tailed discussion of these effects, see [26].
Model E20A possesses the largest GW amplitude of all
our models. In addition, it reaches a high
b
at core bounce
(see Table I) and settles at a postbounce
f
of 9%.A
previous Newtonian study [27] has found a low-T=jWj
nonaxisymmetric instability for a PNS with similar
f
.
In order to verify their findings, we trace the evolution of
model E20A to 70 ms after bounce and perform an analysis
of azimuthal density modes / e
im’
by computing complex
Fourier amplitudes C
m
1
2
R
2
0
$; ’; z 0e
im’
d’
on rings of constant coordinate radius with respect to the
coordinate center of mass. The latter stays within the
innermost zone at all times. In the top panel of Fig. 3,we
display the normalized mode amplitudes extracted at
15 km radius. Without adding artificial seed perturbations
to model E20A, discretization errors trigger m f1; 2; 3g
modes, which rise to a level of 10
5
during the 220 ms
collapse phase. After bounce, the m 1 mode exhibits the
fastest growth. This growth on a dynamical time scale,
lasting over tens of ms until saturation, is closely followed
by a growth of m f2; 3g daughter modes [27,28]. Note
that after core bounce, model E20A remains dynamically
stable to the m 4 grid mode. In the lower panel of Fig. 3,
we plot the GW strains h
and h
as seen along the polar
axis. The rotational symmetry of E20A before and shortly
after bounce is reflected in the fact that h
and h
are
practically zero until E20A develops considerable nonax-
isymmetry when the m 1 mode becomes dominant and
its m 2, GW-emitting harmonic reaches a sizable ampli-
tude. In agreement with expectations for a spinning bar, h
and h
oscillate at the same frequency ( 930 Hz) and are
phase-shifted by a quarter cycle.
Discussion.Our results indicate that the GW signature
of the collapse, bounce, and early postbounce phases of the
core-collapse supernova evolution is much more generic
than previously thought. We find that the dynamics of core
bounce are dominated by mainly gravity and microphysics,
reducing the importance of centrifugal support for the
rotation rates considered here. Importantly, for our model
set, we do not observe rotationally induced multiple core
bounce at subnuclear density as proposed by previous
studies that did not include a microphysical EOS and
electron capture treatment in combination with GR.
t-t
bounce
(ms)
h
+/×,pole
(10
−21
at 10 kpc)
-10 0 10 20 30 40 50 60 70
-3
-2
-1
0
1
2
h
+ ,pole
h
×,pole
A
m
10
−5
10
−4
10
−3
10
−2
10
−1
10
0
m= 1
m= 2
m= 3
m= 4
FIG. 3 (color). Normalized mode amplitudes A
m
jC
m
j=C
0
at postbounce times (upper panel) and GW strains h
and h
along the poles (lower panel) for model E20A.
f (Hz)
h
char
at 10 kpc
LIGO I
advanced
LIGO
10
2
10
3
10
−22
10
−21
10
−20
E20A
s20A2B1
s20A2B2
s20A2B4
s20A1B5
s20A3B3
FIG. 4 (color). Spectra of the characteristic GW strain h
char
of
all models and the LIGO (optimal) rms noise curves [32].
PRL 98, 261101 (2007)
PHYSICAL REVIEW LETTERS
week ending
29 JUNE 2007
261101-3

Thus, we predict that the core-bounce waveform of models
in a large parameter space of initial rotation rate and degree
of differential rotation will likely both qualitatively and
quantitatively resemble those presented in Fig. 2.
Model E20A, which we evolve to later postbounce
times, exhibits the dynamical growth of a nonaxisymmetric
low-T=jWj corotation-type m 1 instability [2729]. We
also find m f2; 3g contributions and significant GW
emission from the quadrupole components. We emphasize
that we observe this instability not only in E20A, but also
in other models with comparable values of
f
. Our results,
which remove the limitations of previous studies
[10,27,30,31], demonstrate that the development of non-
axisymmetric structures is neither limited to Newtonian
gravity, simple matter models, equilibrium configurations,
nor high values of , but may rather be a phenomenon
occurring generically in differentially rotating compact
stars.
For the GW signals from the axisymmetric collapse and
core-bounce phase, we obtain peak amplitudes of up to h
10
20
at 10 kpc, while the nonaxisymmetric structures in
model E20A developing later emit GWs with h only a
factor 5 smaller. However, since the latter emission
process operates over several tens of ms, the total energy
E
gw
emitted in GWs is larger than that from the
core-bounce signal. This is evident in Fig. 4, where we
display the characteristic GW strain spectra h
char
R
1

2
2
dE
gw
=df
q
[5] for all models, evolving E20A
for 70 ms after bounce (see also Table I). Considering
only the core-bounce waveforms, h
char
has its maximum
between 300 and 800 Hz, while for model E20A it peaks at
930 Hz, the frequency of the GW-emitting component of
its nonaxisymmetric structures. We conclude that the core-
bounce GW signals of all models investigated here may be
detectable by current and future LIGO-class detectors from
anywhere in the Milky Way. Models that develop non-
axisymmetric instabilities may be detectable out to much
larger distances if the instability persists for a sufficiently
long time.
We point out that due to the nature of the approximation
used for the neutrino effects, we can only accurately model
the GW emission in the collapse, bounce, and early post-
bounce epoch of the core-collapse supernova scenario. In
future work, we plan to improve upon this and carry out
longer-term postbounce evolutions, where other GW emis-
sion mechanisms may be relevant [6,7].
We thank A. Burrows, E. Mu
¨
ller, S. Ou, L. Rezzolla,
E. Seidel, D. Shoemaker, N. Stergioulas, and J. Tohline for
help and stimulating discussions. This research was par-
tially supported by the DFG (Nos. SFB/TR7, SFB 375) and
by the NCSA (Grant No. AST050022N).
*cott@as.arizona.edu
[1] R. Mo
¨
nchmeyer, G. Scha
¨
fer, E. Mu
¨
ller, and R. Kates,
Astron. Astrophys. 246, 417 (1991).
[2] T. Zwerger and E. Mu
¨
ller, Astron. Astrophys. 320, 209
(1997).
[3] M. Rampp, E. Mu
¨
ller, and M. Ruffert, Astron. Astrophys.
332, 969 (1998).
[4] K. Kotake, S. Yamada, K. Sato, K. Sumiyoshi, H. Ono,
and H. Suzuki, Phys. Rev. D 69, 124004 (2004).
[5] C. D. Ott, A. Burrows, E. Livne, and R. Walder,
Astrophys. J. 600, 834 (2004).
[6] C. D. Ott, A. Burrows, L. Dessart, and E. Livne, Phys. Rev.
Lett. 96, 201102 (2006).
[7] E. Mu
¨
ller, M. Rampp, R. Buras, H.-T. Janka, and D. H.
Shoemaker, Astrophys. J. 603, 221 (2004).
[8] H. Dimmelmeier, J. A. Font, and E. Mu
¨
ller, Astron.
Astrophys. 388, 917 (2002); 393, 523 (2002).
[9] M. Shibata and Y.-I. Sekiguchi, Phys. Rev. D 69, 084024
(2004).
[10] M. Shibata and Y.-I. Sekiguchi, Phys. Rev. D 71, 024014
(2005).
[11] P. Cerda
´
-Dura
´
n, G. Faye, H. Dimmelmeier, J. A. Font,
J. M. Iba
´
n
˜
ez, E. Mu
¨
ller, and G. Scha
¨
fer, Astron.
Astrophys. 439, 1033 (2005).
[12] http://www.cactuscode.org.
[13] E. Schnetter, S. H. Hawley, and I. Hawke, Classical
Quantum Gravity 21, 1465 (2004).
[14] T. W. Baumgarte and S. L. Shapiro, Phys. Rep. 376,41
(2003).
[15] M. Shibata, Astrophys. J. 595, 992 (2003).
[16] L. Baiotti, I. Hawke, P. J. Montero, F. Lo
¨
ffler, L. Rezzolla,
N. Stergioulas, J. A. Font, and E. Seidel, Phys. Rev. D 71,
024035 (2005).
[17] H. Dimmelmeier, J. Novak, J. A. Font, J. M. Iba
´
n
˜
ez, and
E. Mu
¨
ller, Phys. Rev. D
71, 064023 (2005).
[18] H. Shen, H. Toki, K. Oyamatsu, and K. Sumiyoshi, Prog.
Theor. Phys. 100, 1013 (1998).
[19] A. Marek, H.-T. Janka, R. Buras, M. Liebendo
¨
rfer, and
M. Rampp, Astron. Astrophys. 443, 201 (2005).
[20] M. Liebendo
¨
rfer, Astrophys. J. 633, 1042 (2005).
[21] R. Buras, H.-T. Janka, M. Rampp, and K. Kifonidis,
Astron. Astrophys. 457, 281 (2006).
[22] A. Heger, S. E. Woosley, and H. C. Spruit, Astrophys. J.
626, 350 (2005).
[23] C. D. Ott, A. Burrows, T. A. Thompson, E. Livne, and
R. Walder, Astrophys. J. Suppl. Ser. 164, 130 (2006).
[24] S. E. Woosley, A. Heger, and T. A. Weaver, Rev. Mod.
Phys. 74, 1015 (2002).
[25] A. Heger, N. Langer, and S. E. Woosley, Astrophys. J. 528,
368 (2000).
[26] H. Dimmelmeier, C. D. Ott, H.-T. Janka, A. Marek, and
E. Mu
¨
ller, Phys. Rev. Lett. (to be published).
[27] C. D. Ott, S. Ou, J. E. Tohline, and A. Burrows, Astrophys.
J. Lett. 625, L119 (2005).
[28] S. Ou and J. E. Tohline, Astrophys. J. 651, 1068 (2006).
[29] M. Saijo and S. Yoshida, Mon. Not. R. Astron. Soc. 368,
1429 (2006).
[30] J. M. Centrella, K. C. B. New, L. L. Lowe, and J. D. Brown,
Astrophys. J. Lett. 550, L193 (2001).
[31] M. Shibata, S. Karino, and Y. Eriguchi, Mon. Not. R.
Astron. Soc. 334, L27 (2002).
[32] D. H. Shoemaker (private communication).
PRL 98, 261101 (2007)
PHYSICAL REVIEW LETTERS
week ending
29 JUNE 2007
261101-4
Citations
More filters
Journal ArticleDOI
TL;DR: In this paper, the effects of rotation, progenitor structure, and postbounce neutrino leakage on dynamics, GW, and neutrinos signals were investigated in rotating iron core collapse.
Abstract: We present results from a new set of 3D general-relativistic hydrodynamic simulations of rotating iron core collapse. We assume octant symmetry and focus on axisymmetric collapse, bounce, the early postbounce evolution, and the associated gravitational wave (GW) and neutrino signals. We employ a finite-temperature nuclear equation of state, parameterized electron capture in the collapse phase, and a multi-species neutrino leakage scheme after bounce. The latter captures the important effects of deleptonization, neutrino cooling and heating and enables approximate predictions for the neutrino luminosities in the early evolution after core bounce. We consider 12_⊙ and 40_⊙ presupernova models and systematically study the effects of (i) rotation, (ii) progenitor structure, and (iii) postbounce neutrino leakage on dynamics, GW, and, neutrino signals. We demonstrate, that the GW signal of rapidly rotating core collapse is practically independent of progenitor mass and precollapse structure. Moreover, we show that the effects of neutrino leakage on the GW signal are strong only in nonrotating or slowly rotating models in which GW emission is not dominated by inner core dynamics. In rapidly rotating cores, core bounce of the centrifugally-deformed inner core excites the fundamental quadrupole pulsation mode of the nascent protoneutron star. The ensuing global oscillations (f ~700-800 Hz) lead to pronounced oscillations in the GW signal and correlated strong variations in the rising luminosities of antineutrino and heavy-lepton neutrinos. We find these features in cores that collapse to protoneutron stars with spin periods ≾ 2.5 ms and rotational energies sufficient to drive hyper-energetic core-collapse supernova explosions. Hence, joint GW + neutrino observations of a core collapse event could deliver strong evidence for or against rapid core rotation. Our estimates suggest that the GW signal should be detectable throughout the Milky Way by advanced laser-interferometer GW observatories, but a water-Cherenkov neutrino detector would have to be of near-megaton size to observe the variations in the early neutrino luminosities from a core collapse event at 1 kpc.

113 citations

Journal ArticleDOI
TL;DR: In this article, a reformulation of the elliptic sector of CFC is presented, which has the fundamental property of overcoming local uniqueness problems, and an extension of these ideas to FCF.
Abstract: The fully constrained formulation (FCF) proposed by Bonazzola, Gourgoulhon, Grand‐clement, and Novak is one of the constrained formulations of Einstein’s equations. It contains as an approximation the conformal flatness condition (CFC). The elliptic part of the FCF basically shares the same differential operators as the elliptic equations in the CFC scheme. We present here a reformulation of the elliptic sector of CFC that has the fundamental property of overcoming local uniqueness problems, and an extension of these ideas to FCF.

111 citations

Journal ArticleDOI
TL;DR: In this paper, the authors present results from the first generation of multi-dimensional hydrodynamic core-collapse simulations in full general relativity (GR) that include an approximate treatment of neutrino transport.
Abstract: We present results from the first generation of multi-dimensional hydrodynamic core-collapse simulations in full general relativity (GR) that include an approximate treatment of neutrino transport. Using an M1 closure scheme with an analytic variable Eddington factor, we solve the energy-independent set of radiation energy and momentum based on the Thorne's momentum formalism. Our newly developed code is designed to evolve the Einstein field equation together with the GR radiation hydrodynamic equations. We follow the dynamics starting from the onset of gravitational core collapse of a 15 M ☉ star, through bounce, up to about 100 ms postbounce in this study. By computing four models that differ according to 1D to 3D and by switching from special relativistic (SR) to GR hydrodynamics, we study how the spacial multi-dimensionality and GR would affect the dynamics in the early postbounce phase. Our 3D results support the anticipation in previous 1D results that the neutrino luminosity and average neutrino energy of any neutrino flavor in the postbounce phase increase when switching from SR to GR hydrodynamics. This is because the deeper gravitational well of GR produces more compact core structures, and thus hotter neutrino spheres at smaller radii. By analyzing the residency timescale to the neutrino-heating timescale in the gain region, we show that the criterion to initiate neutrino-driven explosions can be most easily satisfied in 3D models, irrespective of SR or GR hydrodynamics. Our results suggest that the combination of GR and 3D hydrodynamics provides the most favorable condition to drive a robust neutrino-driven explosion.

109 citations

Journal ArticleDOI
TL;DR: In this paper, the potential sensitivity of prospective detection scenarios for GWs from CCSNe within 5 Mpc, using realistic noise at the predicted sensitivity of the Advanced LIGO and Advanced Virgo detectors for 2015, 2017, and 2019.
Abstract: The next galactic core-collapse supernova (CCSN) has already exploded, and its electromagnetic (EM) waves, neutrinos, and gravitational waves (GWs) may arrive at any moment. We present an extensive study on the potential sensitivity of prospective detection scenarios for GWs from CCSNe within 5 Mpc, using realistic noise at the predicted sensitivity of the Advanced LIGO and Advanced Virgo detectors for 2015, 2017, and 2019. We quantify the detectability of GWs from CCSNe within the Milky Way and Large Magellanic Cloud, for which there will be an observed neutrino burst. We also consider extreme GW emission scenarios for more distant CCSNe with an associated EM signature. We find that a three-detector network at design sensitivity will be able to detect neutrino-driven CCSN explosions out to ∼5.5 kpc, while rapidly rotating core collapse will be detectable out to the Large Magellanic Cloud at 50 kpc. Of the phenomenological models for extreme GW emission scenarios considered in this study, such as long-lived bar-mode instabilities and disk fragmentation instabilities, all models considered will be detectable out to M31 at 0.77 Mpc, while the most extreme models will be detectable out to M82 at 3.52 Mpc and beyond.

106 citations

Journal ArticleDOI
TL;DR: In this article, a series of two-dimensional, special relativistic magnetohydrodynamic (MHD) simulations were performed to study signatures of gravitational waves in the MHD-driven core-collapse supernovae.
Abstract: By performing a series of two-dimensional, special relativistic magnetohydrodynamic (MHD) simulations, we study signatures of gravitational waves (GWs) in the MHD-driven core-collapse supernovae. In order to extract the gravitational waveforms, we present a stress formula including contributions both from magnetic fields and special relativistic corrections. By changing the precollapse magnetic fields and initial angular momentum distributions parametrically, we compute 12 models. As for the microphysics, a realistic equation of state is employed and the neutrino cooling is taken into account via a multiflavor neutrino leakage scheme. With these computations, we find that the total GW amplitudes show a monotonic increase after bounce for models with a strong precollapse magnetic field (1012 G) and with a rapid rotation imposed. We show that this trend stems both from the kinetic contribution of MHD outflows with large radial velocities and also from the magnetic contribution dominated by the toroidal magnetic fields that predominantly trigger MHD explosions. For models with weaker initial magnetic fields, the total GW amplitudes after bounce stay close to zero because the contribution from the magnetic fields cancels with the contribution from the hydrodynamic counterpart. These features can be clearly understood with a careful analysis of the explosion dynamics. We point out that the GW signals with the increasing trend, possibly visible to the next-generation detectors for a Galactic supernova, would be associated with MHD explosions with explosion energies exceeding 1051 erg.

101 citations


Cites background from "3D collapse of rotating stellar iro..."

  • ...…likely to appear than type I, because a combination of general relativity (GR) and electron capture near core bounce suppresses multiple bounce in the type II waveforms 1http://www.ego-gw.it/ 2http://geo600.aei.mpg.de/ 3http://www.gravity.uwa.edu.au/ (Dimmelmeier et al. 2007; Ott et al. 2007a,b)....

    [...]