Issue |
A&A
Volume 496, Number 2, March III 2009
|
|
---|---|---|
Page(s) | 475 - 494 | |
Section | Stellar structure and evolution | |
DOI | https://doi.org/10.1051/0004-6361/200810883 | |
Published online | 09 February 2009 |
Equation-of-state dependent features in shock-oscillation modulated neutrino and gravitational-wave signals from supernovae
A. Marek - H.-T. Janka - E. Müller
Max-Planck-Institut für Astrophysik, Karl-Schwarzschild-Str.1, 85748 Garching, Germany
Received 29 August 2008 / Accepted 28 January 2009
Abstract
We present two-dimensional (axisymmetric)
neutrino-hydrodynamic simulations of the long-time
accretion phase of a 15
progenitor star after core
bounce and before the launch of a supernova explosion, when
non-radial hydrodynamic instabilities like convection occur
in different regions of the collapsing stellar core and
the standing accretion shock instability (SASI) leads
to large-amplitude oscillations of the stalled shock with a
period of tens of milliseconds. Our simulations were
performed with the P ROMETHEUS-V ERTEX code, which includes a
multi-flavor, energy-dependent neutrino transport scheme and employs
an effective relativistic gravitational potential. Testing
the influence of a stiff and a soft equation of state for hot
neutron star matter, we find that the non-radial mass motions
in the supernova core impose a time variability on the neutrino
and gravitational-wave signals with larger amplitudes, as well as
higher frequencies in the case of a more compact nascent neutron
star.
After the prompt shock-breakout burst of electron neutrinos, a
more compact accreting remnant produces higher neutrino luminosities
and higher mean neutrino energies. The observable neutrino emission
in the SASI sloshing direction exhibits a modulation of several
ten percent in the luminosities and around 1 MeV in the mean
energies with most power at typical SASI frequencies
between roughly 20 and 100 Hz. The modulation
is caused by quasi-periodic variations in the mass accretion rate
of the neutron star in each hemisphere. At times later
than
50-100 ms after bounce, the gravitational-wave
amplitude is dominated by the growing low-frequency (
200 Hz)
signal associated with anisotropic neutrino emission.
A high-frequency wave signal results from nonradial
gas flows in the outer layers of the anisotropically accreting
neutron star. Right after bounce such nonradial mass motions occur
due to prompt post-shock convection in both considered cases
and contribute mostly to the early wave production around 100 Hz.
Later they are instigated by the SASI and by convective overturn
that vigorously stir the neutrino-heating and cooling layers,
and also by
convective activity developing below the neutrinosphere. The
gravitational-wave power then peaks at about 300-800 Hz, connected
to changes in the mass quadrupole moment on a timescale of milliseconds.
Distinctively higher spectral frequencies originate from the
more compact and more rapidly contracting neutron star. Both the
neutrino and gravitational-wave emission therefore carry information
that is characteristic of the properties of the nuclear equation of
state in the hot remnant. The detectability of the SASI effects
in the neutrino and gravitational-wave signals is briefly discussed.
Key words: stars: supernovae: general - hydrodynamics - neutrinos - gravitational waves - dense matter
1 Introduction
Neutrinos and gravitational waves are the most direct potential probes of the processes that occur deep inside of a dying star, accompanying or causing the initiation of the stellar explosion. Neutrinos were already detected in connection with supernova SN 1987A (Bionta et al. 1987; Hirata et al. 1987; Alexeyev et al. 1988), although with poor statistics so that the extraction of information for constraining the explosion mechanism was not possible. The simultaneous measurement of signals of both types remains a very realistic hope for the next Galactic supernova.
Capturing neutrinos and gravitational waves from the same source has the advantage of providing complementary insight into the conditions of the stellar core. While neutrino signals reflect the density structure and thermodynamic conditions in the high-density plasma of the collapsing core and forming neutron star, gravitational waves carry crucial information about the dynamics and nonradial motions of the stellar matter, e.g. of its rotational state or of hydrodynamic instabilities such as convection that deform and stir the condensing central compact remnant of the explosion.
The prompt burst of electron neutrinos, for example, signals the breakout of the supernova shock from the neutrinosphere and the shock heating of matter in this region. The neutrino emission after core bounce is a sensitive probe of the mass accretion rate on the forming neutron star, hence of the density structure of the infalling layers of the dying star (e.g., Liebendörfer et al. 2003; Buras et al. 2006b). The onset of the explosion is expected to show up as a more or less sudden drop in the mass accretion rate of the forming neutron star. A soft supernuclear equation of state will lead to a more compact and hotter remnant, radiating higher neutrino luminosities and more energetic neutrinos (e.g., Janka et al. 2005; Marek 2007). And a possible phase transition to non-nucleonic matter in the neutron star core may impose characteristic features like a second neutrino burst (Sagert et al. 2008) or, if triggering the collapse to a black hole, may cause an abrupt termination of the neutrino emission (e.g., Burrows 1988; Keil & Janka 1995; Sumiyoshi et al. 2006, 2007; Fischer et al. 2008).
Theoretical work on the gravitational-wave signals from stellar core collapse and explosion has a long history of successively refined numerical models. In particular the infall and bounce phases, which are theoretically relatively well understood parts of the evolution, have received a lot of interest, because a strong and characteristic signature could make them a promising source of a detectable gravitational-wave burst (for a review-like introduction to the subject and a nearly complete list of publications, see Dimmelmeier et al. 2008). For this to be the case, the core of the progenitor star must develop a sufficiently large deformation during its infall, and for that it must rotate enough rapidly. This, however, does not seem to be compatible with predictions from the latest generation of stellar evolution models for the vast majority of massive stars, which are expected to have lost most of their angular momentum before collapse (Heger et al. 2005). Only in rare, very special cases, possibly accounting for the few tenths of a percent of all stellar core collapses that produce gamma-ray bursts, such stars seem to be able to retain a high angular momentum in their core and to thus produce relativistic jets and highly asymmetric and extremely energetic, probably magnetohydrodynamically driven explosions (for a recent review, see e.g. Woosley & Bloom 2006). If this hypothetical connection was true, the core bounce phase would not really offer grand perspectives for the measurement of gravitational waves.
In contrast, only relatively little work has so far been done on determining the wave signals from the post-bounce evolution of a supernova, although these signals are likely to carry important information about the still incompletely understood explosion mechanism and the associated core dynamics (see the review by Ott 2008). Müller & Janka (1997) and later Müller et al. (2004) on the basis of significantly improved numerical models, showed that in the case of delayed, neutrino-driven explosions convective overturn behind the stalled shock as well as convection inside the nascent neutron star can account for sizable gravitational-wave signals, which should be detectable with a high probability from a Galactic supernova when the Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO II) is running (see also Fryer et al. 2002, 2004).
While the delayed neutrino-heating mechanism relies on the support by
strong nonradial hydrodynamic instabilities in the region of
neutrino-energy deposition, magnetohydrodynamic explosions, if
linked to rapid rotation, are expected to
develop relatively soon after core bounce and thus to occur faster
than the growth of the mentioned nonradial instabilities.
They tap the reservoir of differential
rotation in a collapsing stellar environment and therefore
require rapidly spinning cores (see Burrows et al. 2007b;
Thompson et al. 2005). Thus they are good candidates for sizable
gravitational-wave pulses from the moment of
core bounce (e.g., Kotake et al. 2006; Ott et al. 2004).
The newly proposed acoustic explosion mechanism, in which
large-amplitude core gravity-mode oscillations of the neutron star
convert accretion power to pressure and shock waves that feed the
supernova shock with acoustic energy (Burrows et al. 2006, 2007a),
is predicted to be associated with enormous and very characteristic
gravitational-wave activity due to the fast periodic movement of
roughly a solar mass of dense matter at late times (1 s)
after core bounce (Ott et al. 2006).
In the present paper our focus is on an analysis of the features
in the neutrino and gravitational-wave signals that might provide
evidence for the action of the so-called
standing accretion shock instability
(SASI), which has been shown to lead to low- mode (in terms of an
expansion in spherical harmonics with order
), large-amplitude
nonradial shock deformation and violent sloshing motions of the
stalled supernova shock (Blondin et al. 2003; Blondin &
Mezzacappa 2007; Bruenn et al. 2006;
Scheck et al. 2004, 2006, 2008; Ohnishi et al. 2006; Foglizzo et al. 2007; Yamasaki & Foglizzo 2008). This
SASI activity is found to reach the nonlinear regime at
roughly 100 ms
after core bounce and to grow in amplitude over possibly hundreds
of milliseconds (Scheck et al. 2008; Marek & Janka 2007;
Burrows et al. 2006, 2007a). It does not only act as seed of powerful
secondary convection but can also provide crucial aid for
the neutrino-heating mechanism by pushing the accretion
shock to larger radii and by thus stretching the time accreted
matter is exposed to neutrino heating in the gain layer
(Buras et al. 2006b;
Scheck et al. 2008;
Marek & Janka 2007;
Murphy & Burrows 2008). Moreover, the SASI is found to
play the driving force of the g-modes pulsations of the neutron
star core that are the essential ingredient of the acoustic
mechanism.
Observational signatures of the presence of the SASI would
therefore be extremely important for our understanding of how
massive stars begin their explosion. Supernova asymmetries and pulsar
kicks are one, yet not unambiguous observational hint.
Neutrinos and gravitational waves,
which originate directly from the region where the blast is
initiated, however, may remain the only way to obtain direct
information. It is therefore a highly relevant question to ask whether
any characteristic structures are imprinted on the neutrino and
gravitational-wave emission by the SASI activity in the supernova
core. In contrast to convection, which exhibits the fastest
growth for the higher- modes
(see Foglizzo et al. 2006), the SASI is
expected to possess the largest growth rates for the dipolar
and quadrupolar deformations (corresponding to
;
Blondin
& Mezzacappa 2006; Foglizzo et al. 2007;
Yamasaki & Foglizzo 2008;
Foglizzo 2001, 2002). Even in the fully nonlinear situation the
geometry and motion of the shock and post-shock layer are found
to be governed by these lowest modes.
We present here an analysis of the neutrino and gravitational-wave
signals that are calculated on the basis of the two-dimensional
post-bounce and
pre-explosion simulations of a 15
star recently published
by Marek & Janka (2007). We constrain ourselves to two
nonrotating models, which allow us to discuss the differences
that can be expected from a stiff and a soft nuclear equation
of state (results of the corresponding 1D simulations can also be found
in Marek & Janka 2007). We find that the SASI sloshing of
the shock and the associated quasi-periodic mass motions lead to a
time-modulation of the neutrino emission and to gravitational-wave
amplitudes whose size and characteristic frequency depend
on the compactness of the proto-neutron star during the first half of a
second after core bounce. Though the SASI contributions to the power
spectra are superimposed by a significant high-frequency ``noise''
due to convective fluctuations, at least the combination of
measurements should make it possible to identify the SASI activity
in the supernova core.
Our paper is organized as follows. In Sect. 2 we will briefly outline the main numerical and physics ingredients of our simulations, and the basic features of the two simulations we compare. In Sect. 3 we discuss our results with respect to the SASI effects on the shock motion, neutrino emission, and gravitational-wave signal, and in Sect. 4 we will summarize our findings and draw conclusions, including a discussion whether the SASI modulations are detectable by neutrino and gravitational-wave experiments.
2 Code, input, models
The 2D simulations discussed in this paper were performed with the P ROMETHEUS-V ERTEX code, whose numerical aspects and the implemented microphysics were described by Rampp & Janka (2002) and Buras et al. (2006a), and the publications quoted in those papers. The detailed list of ingredients was also provided by Marek & Janka (2007), where the simulations were already introduced as Models M15LS-2D and M15HW-2D and compared with other cases in a greater set of calculations. We therefore repeat only a few essential aspects of immediate relevance here and refer the reader to Sect. 2 of the Marek & Janka (2007) paper for more complete information.
The hydrodynamics module of the code is based on a conservative and explicit Eulerian implementation of a Godunov-type scheme with higher order spatial and temporal accuracy. It solves the nonrelativistic equations of motions for the stellar fluid, whose self-gravity is described by an ``effective relativistic potential'' for an approximative treatment of general relativistic gravity (see Marek et al. 2006, the discussed simulations were performed with the potential of Case A from this work). The neutrino transport, which is coupled to the hydrodynamics part via lepton number, energy, and momentum source terms, is computed with our ``ray-by-ray plus scheme'' (see Buras et al. 2006a). It accounts for the full neutrino-energy dependence of the transport but treats its dependence on the direction of the neutrino momentum in an approximative way, which is numerically less demanding and more efficient than a full multi-dimensional version of the transport. Recent multi-angle simulations in two-dimensional situations (Ott et al. 2008) - though done at the expense of a sophisticated description of the energy dependence - show that a detailed angular treatment produces considerably less lateral smearing of the outward directed radiation field than flux-limited diffusion. We suspect that our ray-by-ray description compares much better with multi-angle results than flux-limited diffusion does.
The progenitor star used for our simulations was
model s15s7b2 from Woosley & Weaver (1995), which is a standard
nonrotating
15
star widely used for supernova simulations. We
employed two different nuclear equations of state (EoS) for our
studies: (1) a soft version of the Lattimer & Swesty (1991)
EoS (``L&S EoS'') with
an incompressibility modulus of bulk nuclear matter of 180 MeV
and a symmetry energy parameter of 29.3 MeV, and (2)
the considerably stiffer EoS of Hillebrandt & Wolff
(1985, ``H&W EoS''; see also Hillebrandt et al. 1984), whose
parameter values are 263 MeV and 32.9 MeV,
respectively
,
.
The former leads to a radius of about 12 km for cold neutron stars
with a ``typical'' (gravitational) mass of 1.4
,
whereas
this radius is roughly 14 km in the second case.
The Hillebrandt & Wolff EoS is based on a largely different
modeling approach for inhomogeneous nucleon matter than the L&S EoS
and the more recent EoS of Shen et al. (1998). It employs a
Hartree-Fock calculation in contrast to the compressible liquid
drop model of the L&S EoS and the relativistic mean field
description of the Shen et al. EoS. The three EoSs yield significantly
different results in 1D core-collapse
simulations with respect to the shock formation
point, the luminosities and mean energies of the radiated neutrinos,
and the evolution of the shock radius and neutron star radius
after bounce (see Janka et al. 2005; Marek 2007; Figs. 6 and 7 in
Janka et al. 2007). Since the EoS of
Shen et al. yields intermediate values for many of these
quantities, we consider
the soft Lattimer & Swesty EoS on the one hand and the stiff
Hillebrandt & Wolff EoS on the other as two cases that roughly
span the range of extreme possibilities for baryonic matter around
and above nuclear saturation density near core bounce and shortly
afterwards.
![]() |
Figure 1:
Regions of mass motions due to hydrodynamic instabilities
in the supernova core. The plots show color coded
the lateral velocity ( top) and the radial velocity ( bottom)
in the equatorial plane of the polar grid
as functions of post-bounce time and radius for the 2D simulation with
the L&S EoS ( left) and for the simulation with the H&W EoS
( right). The displayed range of velocity values is limited
to |
Open with DEXTER |
![]() |
Figure 2:
Entropy (black) and electron fraction (red) profiles (the
plotted values correspond to
|
Open with DEXTER |
![]() |
Figure 3:
Four representative snapshots from the 2D simulation with
the L&S EoS at post-bounce times of 247 ms ( top left),
255 ms ( top right), 322 ms ( bottom left),
and 375 ms ( bottom right). The lefthand panel of each figure
shows color-coded the entropy distribution, the righthand panel the
radial velocity component with white and whitish hues denoting
matter at or near rest; black arrows in the righthand panel indicate
the direction of the velocity field in the post-shock region
(arrows were plotted only in regions where the absolute values of
the velocities were less than
|
Open with DEXTER |
3 Results
In the following we will compare our two simulations for these equations of state in detail, focussing on the observable signals whose properties depend on the asymmetries due to hydrodynamic instabilities in the supernova core after bounce.
3.1 Hydrodynamic instabilities and shock motion
Basically one can discriminate three regions and episodes of nonradial hydrodynamic instability in the stellar core after bounce: (i) prompt post-shock convection (Epstein 1979; Burrows & Fryxell 1993; Janka & Müller 1996; Müller & Janka 1997; Swesty & Myra 2005); (ii) convection inside the nascent neutron star, i.e. below the neutrinosphere of the electron neutrinos (Burrows 1987; Keil et al. 1996; Buras et al. 2006b; Dessart et al. 2006); and (iii) convective overturn in the neutrino-heating layer between the gain radius and the stalled supernova shock (Herant et al. 1994; Burrows et al. 1995; Janka & Müller 1996; Fryer & Warren 2002, 2004) as well as SASI activity (Blondin et al. 2003; Blondin & Mezzacappa 2006; Ohnishi et al. 2006; Scheck et al. 2008). These regions can be identified for both 2D simulations in Fig. 1.
Region (i) arises when the newly formed supernova shock
propagates outward through the infalling stellar iron core
and experiences massive energy losses by nuclear photo-disintegrations.
This weakens the shock so that it leaves behind a layer with
negative entropy gradient (before it encounters infalling
material with even lower densities and therefore produces
rising gas entropies again in spite of its further deceleration).
In addition,
a
minimum is established near an enclosed mass of roughly
1
when a large number of newly produced electron neutrinos
begins to stream away from the neutrinosphere in a luminous
shock-breakout burst. The electron fraction forms a trough
around the radial position
where the shock makes the transition from the neutrino-opaque to
the neutrino-transparent regime (Fig. 2).
The electron fraction decreases towards this minimum, because
electron-capture neutrinos escape more slowly from regions with
higher densities, and because the conversion of electrons to
neutrinos proceeds more
slowly at lower densities. In both simulations the negative
entropy gradient overlaps partly with the negative
gradient.
The layer where both the entropy and the
gradients are
negative is found to be Ledoux unstable (marked by the grey
vertical bars in Fig. 2). This region
lies between the first local maximum and the following local
minimum of the entropy profile. It
contains considerably more mass in the model computed with the
H&W EoS, where it is bounded by the mass coordinates of
0.63
and 0.79
.
At higher enclosed masses,
the entropy rises slowly out of its local minimum. In contrast,
in the L&S case the boundaries are at
0.60
and 0.71
,
and the increase of the
entropy in the overlying shells is much steeper. In the H&W model
the thicker layer in mass corresponds also to a much wider radial
domain at a larger distance from the stellar center.
Because of higher growth rates the convective activity sets in
more immediately after the shock passage and lateron becomes
stronger and encompasses a greater fraction of the stellar core.
This can be seen at
ms after bounce in
Fig. 1, where the righthand panels for the
H&W simulation in comparison to the lefthand panels show significantly
more violent convective overturn. While the L&S model develops
visible mass motions between
15 km and
40 km,
the H&W model shows stronger activity at radii of
km.
After some 10 ms the profiles of entropy and
electron fraction are flattened and with the disappearance of
the driving force the convective activity begins to calm down (this
can be recognized better in the two lower panels of
Fig. 1). Now, however, region (ii) is
being formed. Neutrino transport begins to reduce the electron
fraction deeper inside the nascent neutron star. As a consequence
of this the lower boundary of the
trough and thus the
negative
gradient moves gradually inward. At the
same time the energy losses by neutrinos reduce the level of
the entropy plateau at
/nucleon such that
the plateau becomes wider and its inner edge also moves inward.
These changes in the entropy and lepton
number profiles in combination lead to
the establishment (or/and maintenance) of a convectively
unstable layer inside the proto-neutron star.
About 50-60 ms after bounce, following
a short, intermediate period of relative quiescence, both of our
2D models indeed show the reappearance of vivid overturn activity
at densities above 1012 g cm-3
(Fig. 1). The convective
transport of lepton number and entropy supports and enhances
the mentioned structural changes so that a long-lasting zone with
convective mixing develops (compare our Fig. 1
with Figs. 9 and 10 of Dessart et al. 2006).
This zone encompasses a growing mass region within the broadening
trough.
Region (iii) with nonradial hydrodynamic mass motions develops
outside of the steep density gradient at the surface of the
nascent neutron star, where the neutrinospheres of all neutrinos are
located. First signs of low-
mode, low-amplitude SASI sloshing
can be seen here already some 10 ms after bounce. At about
100 ms after bounce this activity strengthens significantly
because at this time neutrino energy deposition in the gain layer
has begun to build up a negative entropy gradient between the gain
radius and the
stagnant supernova shock. The neutrino-heating region is unstable
to convection, and the presence of convective motions produces
vorticity and entropy perturbations that feed back into
the advective-acoustic cycle that is considered as an explanation of
the SASI phenomenon (see Foglizzo 2001, 2002;
Foglizzo et al. 2007).
Thus the onset of convection may lead to an amplification of the
SASI growth. On the other hand, the violent shock oscillations that
are characteristic of the fully developed SASI produce large entropy
variations in the downstream region. These act as seeds for secondary
convection (Scheck et al. 2008) and have been found to aid the
neutrino-heating mechanism for powering supernova explosions
(Scheck et al. 2008; Marek & Janka 2007;
see also Murphy & Burrows 2008).
In the nonlinear regime convective and SASI activity are
inseparably coupled, and a high-
convective mode pattern
occurs superimposed on the low-
mode SASI deformations of
the shock contour and of the post-shock region (for
snapshots, see Fig. 3 and Marek & Janka 2007;
Scheck et al. 2008).
![]() |
Figure 4: Left: maximum shock radii (solid lines) and proto-neutron star radii (dashed lines) as functions of post-bounce time for the 2D simulations with different nuclear equations of state. The neutron star radii are determined as the locations where the rest-mass density is equal to 1011 g cm-3. Right: shock contours at the different post-bounce times listed in the figure. The vertical axis of the plot is the symmetry axis of the simulation. |
Open with DEXTER |
Although both simulations show the same basic features
and evolutionary stages, the differences in details are interesting.
As mentioned above, the H&W model develops more vigorous prompt
post-shock convection in a region with higher mass and wider
radial extension. In contrast, the SASI activity in this model is
appreciably less strong than that in the calculation with the
L&S EoS until roughly 350 ms after bounce. This can be seen in
Fig. 4, where the maximum shock radius of
the L&S case exhibits significantly bigger SASI amplitudes and
correspondingly the shock contours show more extreme nonspherical
deformation (righthand panel of Fig. 4). These
differences grow during the nonlinear phase of the SASI between
150 and 350 ms until at
ms p.b. the
SASI sloshing in the H&W model also gains more power and the
maximum shock radii of both simulations become more similar again.
We suspect that the more compact proto-neutron star for the softer L&S EoS leads to conditions that favor strong SASI activity, possibly because of the more efficient neutrino heating (and more vigorous post-shock convection) as a consequence of the higher neutrino luminosities and harder neutrino spectra that are radiated from a more compact and hotter nascent neutron star (see Figs. 6 and 7 and Sect. 3.2). Another possible explanation in the context of the advective-acoustic cycle scenario may be different amplification factors of perturbations in the two models. Due to the lack of good theoretical insight into the behavior of the SASI in the fully nonlinear regime, we do not see a way how to facilitate deeper understanding by further analysis. For both suggested explanations, however, one might expect that in the H&W simulation at late post-bounce times, when the proto-neutron star radius has contracted (Fig. 4) and the neutrino-heating timescale has decreased (see Fig. 6 in Marek & Janka 2007), the conditions for violent SASI activity have the tendency to improve. This would be consistent with our observed growth of the SASI amplitudes near the end of this simulation.
![]() |
Figure 5:
Left: amplitudes for the decomposition of the
angle-dependent position of the
shock front,
|
Open with DEXTER |
In Fig. 5, lefthand panels,
we plot the time-evolution of the
amplitudes A1,0, A2,0, and A3,0, of the decomposition
of the angle-dependent shock radius
into
spherical harmonics:
where m = 0 because of the axial symmetry of our models and

![[*]](/icons/foot_motif.gif)

![]() |
Figure 6:
Isotropic equivalent luminosities of electron neutrinos
( top), electron antineutrinos ( middle), and one kind of heavy-lepton
neutrinos ( |
Open with DEXTER |
The righthand panels of Fig. 5 provide the
Fourier transforms of the time-dependent mode amplitudes. One
can see a broad peak of the Fourier spectra between about 10 and
roughly 100 Hz, followed by a steep decline towards higher
frequencies. The
amplitude exhibits the clearest,
jagged peak structure at a frequency of 20-60 Hz for
the L&S run, and a smoother peak at 30-40 Hz with a
secondary one around
65 Hz for the H&W model. A2,0 shows most power
at frequencies around 30-50 Hz and 70-90 Hz in both
simulations, and A3,0 possesses a broad Fourier maximum
below
100-130 Hz, with an indication of a peak at about
50 Hz in the H&W case. There is no really clear correlation
of the SASI oscillation period with the compactness of the
nascent neutron star or with the average radius of the stalled
shock except for, maybe, the slightly higher frequency of the
main
peak in the L&S run.
Because of the time-dependent structure of the
accretor and of the whole post-shock region the appearance
of well localized and very prominent peaks in the spectrum
might not be expected. In addition, shock wobbling
associated with convective mass motions
produces a background of short-wavelength ``noise'' at frequences
up to roughly 200 Hz (Fig. 3).
Accordingly, in particular the
amplitudes for higher
values show power nearly evenly
distributed over a wide range of frequencies. The Fourier
transforms for A1,0, A2,0, and A3,0 plotted
in Fig. 5 nevertheless confirm the
stronger SASI activity in the simulation with the L&S EoS:
only for a few, narrow frequency windows the Fourier amplitudes
of the H&W model become larger than those of the L&S run.
In the following two sections we will turn to an analysis of the consequences of the described nonradial hydrodynamic instabilities, in particular also of the SASI, for observable signals from the supernova core. One question of interest will be, whether the differences caused by the use of stiff or soft neutron star equations of state manifest themselves in some distinctive features of the signals.
![]() |
Figure 7: Mean energies of radiated neutrinos as functions of post-bounce time for our 1D simulations ( top) and 2D models ( middle and bottom) with both equations of state (the lefthand panels are for the L&S EoS, the right ones for the H&W EoS). The displayed data are defined as ratios of the energy flux to the number flux and correspond to the luminosities plotted with dashed and solid lines in Fig. 6. The panels in the middle show results for a lateral grid zone near the north polar axis, the bottom panels provide results that are averaged over the whole northern hemisphere of the computational grid. In all cases the evaluation has been performed in the laboratory frame at a distance of 400 km from the stellar center. |
Open with DEXTER |
3.2 Neutrino emission
The luminosities for neutrinos and antineutrinos of
all flavors and for both 2D simulations, compared to the results
of the corresponding 1D models, are given as functions
of post-bounce time in Fig. 6 (muon and tau
neutrinos and their antiparticles, which we sometimes denote
with ,
are treated in the same way,
and the plots show the luminosities of one individual type of
these heavy-lepton neutrinos). The luminosities
displayed there are actually not defined as the total rate of energy
loss of the supernova core through neutrinos radiated in all directions.
Instead, they represent the isotropic
equivalent luminosities inferred by an observer at rest relative to
the stellar center, viewing the source from the direction of the
polar axis of the grid and from a large distance above the north pole
(i.e., gravitational redshifting is taken into account). The results
for an observer above the south pole look qualitatively and
quantitatively very similar and lead to the same conclusions.
One of the goals of the present work is an analysis of luminosity
fluctuations that are associated with modulations of the mass
accretion rate of the nascent star due to SASI oscillations and
convective overturn in the post-shock layer,
which become very vigorous and create nonstationary conditions
at later post-bounce times (see Sect. 3.1 and
Figs. 4 and 5).
Strong shock retraction leads to a transient increase of the gas
flow towards the neutron star and to the compression and enhanced
cooling of the matter near the neutron star surface. In contrast,
shock expansion has the opposite
effect because it causes a deceleration of the infall or even outward
acceleration of material that is accreted through the shock. Thus
shock expansion stretches the time this matter remains in the gain
layer and does not cool by neutrino emission
(Marek & Janka 2007; Scheck et al. 2008;
Murphy & Burrows 2008). Corresponding quasi-periodic fluctuations of the neutrino
luminosities, strongest for electron neutrinos
and antineutrinos
(which are more abundantly produced in the layer of
freshly accreted, hot matter near the neutron star surface) can be
seen in all panels of Fig. 6.
In the lefthand panels of Fig. 6 we have constructed
the isotropic equivalent luminosities from the neutrino radiation
leaving the stellar core in one angular grid bin near the north pole
of the mesh, i.e.,
for


where the second line exploits axial symmetry and the sum is performed over the angular zones of one hemisphere. The former approach means that one assumes that an observer at a position near the polar axis receives the radiation emitted only from one pole-near lateral zone (in our treatment this means that the received radiation flux has the properties calculated in exactly one angular grid bin or angular ``ray''), whereas the second evaluation implies the assumption that each unit of area on the hemisphere oriented towards the observer contributes to the observable luminosity with a weight defined by the (radial) flux calculated for the corresponding angular bin (``ray'') of the polar grid in our ray-by-ray transport treatment. This assumption would be correct, of course, if a spherical source were radiating uniformly and isotropically in all directions, but in general the surface parts oriented with smaller angles to the observer direction contribute more strongly. For a sphere with a sharp radiating surface (in contrast to a radially extended neutrino-decoupling layer) that emits neutrinos with a (locally) isotropic intensity towards an observer at great distance on the axis of symmetry of the source, Eq. (3) would contain a factor

![[*]](/icons/foot_motif.gif)

Not solving the full multi-angle problem of neutrino transport in
the 2D geometry but using our ray-by-ray transport approximation
does not allow us to exactly compute the direction
dependence of the neutrino intensity radiated from the nascent
neutron star. The ray-by-ray treatment essentially implies
the construction of a 1D transport solution in every lateral bin
of the computational polar grid, assuming a spherically symmetric
transport problem for the stellar conditions that are present in
this angular bin. Therefore the ray-by-ray approach
tends to underestimate
the directional smearing of luminosity features that result from
local and time-dependent emission increase in the semi-transparent
accretion layer near the neutron star surface. Such features are,
for example, caused by the
SASI and convection modulated mass inflow and must be
expected to be more prominent in the plots showing the measurable
luminosity ,
because this quantity is the isotropic
equivalent luminosity
resulting from neutrino emission occurring in one angular grid bin
(see lefthand panels of Fig. 6). In contrast,
spatially localized effects will show up only with
significant damping in the luminosity
,
where neutrinos
coming from the whole hemisphere facing the observer are added up
with equal weighting (righthand panels of Fig. 6).
We consider these two evaluated cases
as the extrema that define the theoretical limits of the true
result. We suspect that
yields the better approximation
of the real situation, but we also provide
for a
pessimistic estimate of the magnitude of the luminosity
fluctuations. It is important to note that the total energy loss
of the supernova core in neutrinos (or one sort of neutrinos),
which is usually defined as neutrino luminosity and computed by
integrating the neutrino flux over all observer directions,
shows SASI induced variability only on a much reduced level (see
the corresponding plots in the Marek & Janka 2007 paper). This
can be understood from the fact that the SASI fluctuations in the
northern and southern hemispheres are phase shifted by roughly
half a SASI period so that the emission maxima on one side fill
the emission minima on the other. The combined signal is therefore
much smoother than the neutrino flux that can be received by an
observer in any chosen direction.
The modulations of the apparent luminosity due to SASI and convective
variations of the neutron star accretion rate in the 2D simulations
are superimposed on the general luminosity evolution that is
characteristic of the employed progenitor star. These modulations
add another feature on top of the effects that discriminate 1D from
2D results and models based on the use of the L&S EoS from those
with the H&W EoS. The 15
stellar model used in the
present simulations shows a distinct drop in the mass accretion rate
by roughly a factor of two when the steep density decline at the interface
between the silicon shell and the oxygen-enriched Si layer falls
through the stalled shock (see Appendices A and B of Buras et al. 2006b).
This happens between 160 ms and 180 ms
and leads to the pronounced decrease of the luminosities of all neutrinos
and antineutrinos at around this time.
In Figs. 6 and 7 we see
an effect already mentioned in Sect. 3.1, namely that the
proto-neutron star in the simulations with the softer L&S EoS is
more compact and hotter and therefore radiates higher luminosities
and significantly higher mean energies of all kinds of neutrinos
and antineutrinos (the mean neutrino energies
are defined here as
the ratio of the energy flux to the number flux).
This holds for 1D as well as 2D models. Besides the quasi-periodic
fluctuations, the 2D results exhibit the characteristic differences
compared to the 1D case that were discussed before in much detail by
Buras et al. (2006b): convective transport below the neutrinospheres
leads to an accelerated loss of lepton number and energy, which
manifests itself on the one hand in a slightly enhanced electron neutrino
flux relative to the electron antineutrino flux, and on the other hand in
significantly increased muon and tau neutrino luminosities compared
to the 1D results. These differences are visible in both the polar
luminosities,
,
as well as in the hemispherically
integrated data,
,
of Fig. 6.
On the other hand, the convective neutron star is less compact
than its 1D counterpart, for which reason the mean energies of the
radiated neutrinos are slightly reduced in the 2D simulations (compare
the upper two panels of Fig. 7 with the bottom
ones). This
phenomenon is particularly strong for the L&S EoS, in which case the
more compact neutron star structure seems to react more sensitively to
the redistribution of lepton number and entropy associated with the
convective activity in its interior.
An interesting effect is the crossing of the mean energies
of electron antineutrinos and heavy-lepton neutrinos (
becoming higher than
)
that occurs in five of the six panels of Fig. 7 at
about 200 ms after bounce, but which cannot be observed until
this post-bounce time for other energy moments; for example it is
absent in the case of the rms energies of the energy
spectra. These rms energies are ordered in the standard way at
all times, with heavy-lepton neutrinos
being more energetic
on average than
,
and the latter being more energetic
than
:
(this was pointed out already by Marek & Janka 2007).
The energy crossing therefore is a consequence of a subtle difference
in the spectral shape of
and
.
A close inspection
of our results reveals that this effect is connected with the
modification of the radiated neutrino fluxes in the accretion layer
around a successively more compact, contracting proto-neutron star.
In contrast to the mean energies at a radius of 400 km, i.e.
far outside of the neutrino source (see Fig. 7),
the corresponding data at a fixed density of
1011 g cm-3 always follow the relation
.
This holds for our 1D results and
for the time-dependent unidirectional as well as (hemi)spherically
averaged 2D results with both equations of state; only during
short phases of a strong increase of the accretion rate can
the mean energy of electron antineutrinos exceed the one of
heavy-lepton neutrinos also at
g cm-3.
Similarly, accretion around the poles of
the 2D models also causes the polar rms energies of
and
- measured outside of the accretion layer - to converge
(
),
whereas the laterally averaged values are
always well separated, in agreement with the hierarchy of the rms
energies found in the corresponding 1D simulations. Moreover,
the neutrino radiating ``naked'' neutron star in
spherically symmetric models with artifically initiated explosions
does not produce the inversion of the mean
and
energies, whereas the accreting neutron star in a corresponding
nonexploding model does (see Fig. 2 in Janka et al. 2007).
All these findings point to the importance of accretion
for an explanation of the phenomenon. Accretion modifies the
density profile outside of the neutrinosphere and thus is
responsible for
a subtle change in the
spectrum relative to the
spectrum. Two effects seem to contribute. On the one hand,
with the inclusion of
-pair production by
nucleon-nucleon bremsstrahlung and by
annihilation
in state-of-the-art supernova models,
the energysphere of heavy-lepton neutrinos is shifted to
lower temperatures than in older simulations where these processes
were ignored. The heavy-lepton neutrinos that leave
their average neutrino-energysphere also experience a
``spectral filter effect'' in being down-graded in energy space
by neutrino-nucleon and neutrino-electron collisions when diffusing
through the cooler, outer layers of the neutron star (modern
multi-energy transport calculations follow the corresponding
spectral changes with great accuracy).
For both reasons the spectra of
and
are now much more similar than in the older models (for a detailed
discussion, see Raffelt 2001; Keil et al. 2003). On the other
hand, electron neutrinos and antineutrinos, which are emitted in
rapid charged-current reactions of electrons and positrons with protons
and neutrons, are more copiously produced in the hot accretion layer
than muon and tau neutrinos. This leads to an enhancement of the
flux at medium and high energies (
10 MeV) compared
to the
flux at these energies. In combination, these effects
can raise the average energy of electron antineutrinos radiated from
an accreting neutron star above that of heavy-lepton neutrinos.
The described relative spectral changes become stronger at later
times after core bounce so that a crossing of
and
occurs. As the post-bounce time advances,
the average energies of muon and tau neutrinos grow slowly because the
neutrinospheric temperature increases in the shrinking neutron star.
This increase is stronger for the soft L&S EoS than for the
stiff H&W EoS (see Fig. 7). In contrast, the
and
created in the semitransparent accretion layer
escape with mean energies that rise relatively faster, because
accretion onto a more compact proto-neutron star releases more
gravitational binding energy. This leads to more extreme heating
of the infalling and compressed gas and more abundant production of
and
with higher energies. This reasoning also
explains why the polar-near data in Fig. 7
show the effect more clearly than the
hemispherically averaged results: the SASI sloshing motions of the
shock along the direction of the polar grid axis with alternating
periods of shock inflation and contraction lead to successive phases
of gas expansion and mass accumulation on the one hand, and gas
compression on the other, in particular around the north
and south poles of the neutron star (compare Figs. 3
and 6).
At first sight this argument seems in conflict with the observation that the accretion funnels in the sequence of snapshots of Fig. 3, in particular the most prominent and persistent ones that carry most mass, reach down to the neutron star typically at intermediate latitudinal angles, i.e., around the equator; still the fluctuations of the neutrino emission are highest close to the poles. We understand this as a consequence of the following two facts. Firstly, because of the extremely large infall velocities only little neutrino loss occurs within the downdrafts. Instead, most of the binding energy is radiated when the accreted gas spreads around the neutron star and begins to settle onto the neutron star surface (for a discussion of this aspect, see also Scheck et al. 2006). Secondly, due to the SASI sloshing of the whole accretion layer, much of the gas accreted through the shock is redirected alternatingly towards the north pole or the south pole. There the gas is decelerated near the grid (and sloshing) axis and deflected so that its duration of stay is longest in the polar regions. The gas that is sucked to a pole at an enhanced rate during the SASI half-cycle of shock expansion is compressed there during the next half-cycle of shock contraction. This is the reason why most of the losses of accretion energy happen around the poles. Future three-dimensional (3D) simulations with neutrino transport will have to show whether this behavior is a consequence of the assumption of axial symmetry or whether a similar effect can also be found in connection with the SASI in three dimensions.
![]() |
Figure 8: Fourier amplitudes of the time-dependent neutrino luminosities shown in Fig. 6 (the left/right column there corresponds to the left/right column in the present figure). In order to reduce the level of background noise, the evaluation was performed for post-bounce times t > 100 ms, for which the SASI modulation is clearly visible in the neutrino luminosities, and the signals were split in 17 overlapping segments of 128 ms, advanced in steps of 10 ms. The Fourier spectra of the 1D models show a weak contribution from low-amplitude numerical fluctuations of the transport results. |
Open with DEXTER |
Concerning the possibility of measuring the SASI (and convection)
induced luminosity variations by the detection of neutrinos from a future
Galactic supernova it is important to note that the strong compressional
heating of accreted gas during the SASI half-cycles with shock
contraction does not only boost the neutrino emission at the poles
but also leads to the production of more energetic
neutrinos. Therefore the variations of the neutrino luminosities
correlate and are in phase with modulations of the mean neutrino energies
(see in particular the polar data in Figs. 6
and 7). While (for the polar emission) the luminosity
differences between maxima
and minima can be several 10% up to about 50% of the values in
low-emission episodes and the mean energies show variations up to
roughly 1 MeV or 10%, the observationally relevant energy moments
(approximately scaling like
or
,
depending on the neutrino
detector) fluctuate with an amplitude of typically
50%
and extrema up to nearly 100% of the values in the minima.
The fluctuations are similarly strong for
and
,
which are both abundantly created in the accretion layer,
and somewhat reduced for heavy-lepton neutrinos, whose emission is
more indirectly affected by the perturbations of their decoupling
region near the neutron star surface, e.g. through the surface gravity
modes caused by the SASI activity and through changes in the density
structure and optical depth associated with the mass motions in the
SASI region.
In Fig. 8 we provide the Fourier spectra of
the neutrino luminosities displayed in Fig. 6;
the lefthand panels belong to the polar luminosities, the righthand
panels to the hemispherically averaged ones.
The 2D simulations exhibit a broad peak between about 20 Hz and
roughly 400 Hz with a clear power excess compared to 1D results.
Instead of the expected monotonic decline the Fourier transforms of
the 1D luminosities possess a flat shoulder at frequencies between
50 Hz and
100 Hz. This feature is a consequence
of low-amplitude numerical noise produced by the transport scheme,
which is so small that it is physically insignificant for the
simulations. As expected from
the time evolution of the shock radii and corresponding luminosity
variations, which both are larger in the case of the L&S run,
the Fourier spectra of the 2D simulation with the L&S
EoS dominate the result with the H&W EoS by a factor of 2-5 for
all kinds of neutrinos. The flat local maxima in the
25-80 Hz range, whose location roughly
agrees with the peak of the
Fourier amplitude of the
shock oscillation mode
(Fig. 5) and also shows the tendency of slightly
higher frequencies for the L&S case, are superimposed by a
high-frequency component that originates from the faster convective
modulations of the accretion flow and neutrino emission.
![]() |
Figure 9:
Gravitational-wave quadrupole amplitudes
AE220as functions of post-bounce time ( left) and
corresponding Fourier spectra ( right) associated with mass
motions ( top) and anisotropic neutrino (
|
Open with DEXTER |
3.3 Gravitational waves
Nonradial mass motions in the three regions of hydrodynamic instabilities as well as anisotropic neutrino emission are sources of gravitational-wave signals from stellar core-collapse events (e.g., Müller et al. 2004 and references therein). We will analyze our data for EoS-specific properties of these signals now.
The (quadrupole) gravitational-wave amplitudes, energy spectra, and
spectrograms resulting from anisotropic mass motion can be computed
for axisymmetric models as described in
Müller & Janka (1997; Eqs. (10)-(12)), using the Einstein
quadrupole formula in the numerically convenient form derived
by Blanchet et al. (1990), and by standard FFT techniques.
Assuming an observer that is located at an angle with respect to the symmetry axis of the source, the dimensionless
gravitational wave amplitude h(t) is related to the quadrupole wave
amplitude
AE220 (measured in units of cm), the lowest-order
non-vanishing term of a multipole expansion of the radiation field
into pure-spin tensor harmonics (see Eq. (9) of Müller 1997),
according to
where R is the distance to the source. In the following we will always assume



(the quadrupole wave amplitude AE220 for neutrinos then results from inversion of Eq. (5)), with
where













The quadrupole wave amplitudes AE220 of our two 2D models and the corresponding Fourier spectra are shown in Fig. 9, in the upper panels for the signals from aspherical matter flow, and in the lower panels for those from anisotropic neutrino emission.
The matter signal consists of a superposition of
quasi-periodic variations
on timescales ranging from a few to several ten milliseconds.
Two phases of enhanced activity can be distinguished. The first one
occurs within tens of milliseconds after core bounce and
is followed by a more quiescent episode before a long-lasting period
with an overall trend in growing amplitudes begins. The transient,
early phase of strong gravitational-wave emission is a consequence of
the prompt post-shock convection (cf. Sect. 3.1) and
exhibits larger amplitudes, slightly higher frequencies (spectral
peak at about 100 Hz instead of 70 Hz),
and more power in our simulation with the
stiffer H&W EoS (see Fig. 10). The reason
for these differences is the
fact that the region of prompt convection involves more mass and
extends to larger radii in the H&W case (Figs. 1
and 2), where also the growth rates
of the instability are higher and the mass motions become more violent
and faster. The prominent peaks of the Fourier spectra
in Fig. 10 correspond to a clearly dominant
frequency that can be identified in the wave train during the first
80 ms after bounce.
The later phase of intense gravitational-wave production is associated
with the development of convective overturn and large-amplitude SASI
oscillations in the post-shock layer as well as mass flows and
convective activity in the
proto-neutron star. Accordingly, the wave amplitudes at
ms
show variability with different frequencies. It is interesting to note
that the enveloping curves of the wave amplitudes at this time
follow roughly the evolution of the maximum shock radii in
Fig. 4. This means that more active phases and
stronger gravitational-wave emission correlate with larger shock
radii. The nonradial gas motions in the post-shock
region contribute mostly to the power radiated at
frequencies up to about 200 Hz in accordance with the SASI
Fourier spectra shown in Fig. 5. The main peak
of the matter gravitational-wave spectrum at 300-600 Hz (H&W EoS)
and 600-800 Hz (L&S EoS) as well as most of the lower ``background''
above
Hz, however, originates from changes
in the mass-quadrupole moment in deeper regions, namely at densities
between 1011 g cm-3 and 1013 g cm-3,
where more mass is involved than in the post-shock layer.
The main source of this signal are gravity modes and gas flows
instigated by the violent impact of accretion funnels
in the surface layers of the nascent neutron star, which
perturbs regions at densities up to
g cm-3. A substantial part of the
emission also comes from the convective zone inside the neutron
star at
g cm-3. The importance of the
interaction of the SASI and accretion flows with the proto-neutron
star surface can be concluded from a closer inspection of
the upper lefthand panel of Fig. 9 in comparison with
Fig. 1. The gravitational-wave amplitudes
of the L&S run at t > 100 ms p.b. are significantly larger than
those of the H&W model until near the end of our simulations. Only then the amplifying
SASI and accretion activity begins to affect shells at
g cm-3 also in the H&W case,
while earlier the involved regions in this model possess clearly lower
densities.
The maximum values of the dimensionless gravitational-wave strain h(t) for the matter signal are of the order of several 10-22for a source at a distance of 10 kpc, and the dominant peak may be marginally observable with LIGO I in the case of a Galactic supernova (for a detailed discussion of measurement aspects, see Müller et al. 2004). While the spectral peak is slightly higher if the nuclear EoS is soft and the proto-neutron star more compact, the stiffer H&W EoS leads to somewhat lower peak frequencies and thus to a shifting of the main emission closer to the region of highest sensitivity of the LIGO instrument. It is therefore not clear that a soft EoS is more favorable for a measurement of such gravitational-wave signals.
![]() |
Figure 10: Fourier spectra of the gravitational wave quadrupole amplitudes of the mass motions caused by prompt post-shock convection during the first 100 ms after core bounce (see upper lefthand panel of Fig. 9). |
Open with DEXTER |
Significantly larger amplitudes are generated by the anisotropic
neutrino emission, however most power of this component of the
gravitational-wave signal comes at frequencies below 100 Hz
(Fig. 9).
The lower lefthand panel of Fig. 9
shows negative amplitudes with the continuous trend in growing
absolute values and superimposed high-frequency modulations for both
2D simulations discussed in this paper.
The negative values indicate a time-dependent quadrupole moment that
is determined by an excess of neutrino emission from regions near the
equatorial plane, a fact that can be concluded easily from the
latitudinal variation of the function
in
the angular integral of the anisotropy parameter
of
Eq. (7). This function becomes negative in the equatorial
belt between
and
,
see Fig. 1 of Kotake et al. (2007).
A continuous decrease of the neutrino gravitational-wave amplitude to
negative values is no common feature of all 2D simulations. This
behavior depends strongly on the particular emission asymmetry that
develops within the nascent neutron star and in its accretion layer.
Müller et al. (2004) found long-period variations between positive
and negative values, for example,
for an explosion simulation of an 11.2star (cf. their Fig. 4) and for a 2D neutrino-cooling simulation of a
convective proto-neutron star (cf. their Fig. 5), and they obtained a
positive gravitational-wave amplitude with an overall trend in
growth and only shorter periods of decrease for the neutrino emission
of a centrifugally flattened, accreting neutron star at the center
of a collapsing, rotating 15
supernova progenitor (cf. their Fig. 3).
This diversity of possible behaviors is in conflict with results obtained by Kotake et al. (2007), who reported almost monotonically rising, positive wave amplitudes for all of their 2D simulations of SASI unstable accretion shocks in collapsing supernova cores. They explained their finding by the particular properties of the SASI and the corresponding post-shock accretion flow in models that are constrained by the assumption of axial symmetry. In this case the SASI oscillations must occur along the polar grid axis, and Kotake et al. (2007) argued that the influence of this on the accretion flow makes the neutrino emission biased towards the regions around the symmetry axis (see also Sect. 3.2), thus leading to positive and growing wave amplitudes.
Why do we come to different conclusions about the neutrino
emission asymmetry than Kotake et al. (2007) in spite of
our apparent agreement about the role of the SASI for the
and
production in the accretion layer?
In the first place one must note that the modeling approach taken
by Kotake et al. (2007) is fundamentally different from ours.
Besides many other differences, they, for example, did not simulate
the collapse of a stellar progenitor model but started from
a steady-state solution of an accretion shock with
a mass accretion rate of 1
s-1; they
excised the neutron star interior at a constant radius chosen as
the inner boundary of the computational grid, where they assumed a
fixed density of 1011 g cm-3; and they kept the
fairly high
mass accretion rate and the accretor mass constant during all
of the simulation (time-independent boundary conditions at the
outer grid boundary), using a spherically symmetric Newtonian
gravitational potential and ignoring self-gravity of the matter
in the accretion flow. Moreover and very important in the context
of the present discussion of gravitational-wave signals associated
with anisotropic neutrino emission, they did not solve the
energy-dependent, multi-flavor neutrino
transport problem but employed a light-bulb description with
prescribed, time-independent neutrino luminosities and spectra
from the nascent neutron
star and local source rates for lepton number and energy due to
neutrino reactions in the accretion flow.
![]() |
Figure 11:
Quantities that are relevant for the gravitational-wave
signal associated with the anisotropic neutrino emission in our 2D model
using the L&S EoS.
Top left: gravitational-wave amplitudes as functions of time
after core bounce. The black curve was obtained by considering the
neutrino radiation at a large distance from the source (identical with
the corresponding line in
the lower lefthand panel of Fig. 9, where, as in
Figs. 6 and 7, the evaluation
was performed in the laboratory frame at a radius of 400 km).
The red line shows the wave amplitude computed from the neutrino
luminosities that leave the neutrinospheric region
(rather, a sphere where the density is 1011 g cm-3).
The green line displays the difference between the black and red curves.
Top right: anisotropy parameter
|
Open with DEXTER |
Based on the possibilities available to us for comparison we cannot make a reasonable judgement whether any, and if so, which of these simplifications are causal for the discrepant gravitational-wave results obtained by us and by Kotake et al. (2007). However, we strongly suspect that the approach of the latter work is too simplistic for reliably estimating the emission asymmetry of neutrinos that is responsible for the gravitational-wave generation. This conclusion is supported by a close analysis of our neutrino transport data and the insight they yield into the formation of the neutrino emission asymmetry. Figure 11 provides the corresponding information for our 2D simulation with the L&S EoS.
In the upper lefthand panel of Fig. 11 we display
the time evolution of the neutrino gravitational-wave amplitude for
this model. Different from
Fig. 9 we have performed the evaluation now
not only at a radius of 400 km, where the neutrino
luminosity has reached its asymptotic level, but also at the
(time-dependent) radius where the density has a value of
1011 g cm-3, i.e. in the vicinity of the
electron neutrinosphere and thus below the accretion layer where
neutrino cooling still enhances the outgoing
and
fluxes. The upper righthand panel of Fig. 11
contains the corresponding parameters
of the
neutrino emission anisotropy, which exhibit rapid variations
between positive and negative values on timescales of milliseconds
to
10 ms and growing amplitudes at later times.
Both upper panels show an even clearer long-time
trend to negative gravitational-wave amplitudes for the neutrino
luminosity leaving the neutrinospheric layer:
measured at 1011 g cm-3 attains positive values only
during the very short time intervals when the high-frequency
modulations reach their largest amplitudes,
while this parameter at a radius of 400 km enters the positive
side for significantly more numerous and longer time intervals,
although
for all times
at both locations.
Correspondingly, the gravitational-wave amplitude at
g cm-3 is even more negative than the one
at r = 400 km, and their difference (green line in the upper
lefthand
panel of Fig. 11) suggests that the amplitude
associated with the neutrino emission in the accretion layer is
positive and reduces the neutrinospheric emission asymmetry.
This seems to mean that the final emission anisotropy of the neutrinos
is determined already by the flux that leaks out from the
neutrinospheric region and in which the equatorial region has
obtained a slight flux excess. Below
we will see, however, that this interpretation is not completely
correct and the true situation is more complex. A comparion of the
different lines in the upper lefthand panel of Fig. 11
also reveals that the fluctuations of the radiated neutrino
gravitational-wave signal with a period of 10-20 ms
(frequency
50-100 Hz)
originate from the accretion layer and are therefore connected
to the SASI and convective activity in the post-shock region and
its influence on the neutrino-cooling zone as described in
Sect. 3.2.
3.3.1 Neutrino emission anisotropy
In order to obtain deeper insight into the origin of the neutrino
emission anisotropy, we have calculated the cumulative loss of
energy in neutrinos (i.e., the time- and surface-integrated
neutrino flux) once for the polar cap regions and another time
for the equatorial belt, again at r = 400 km and at
g cm-3. The results are shown in
the lower lefthand panel of Fig. 11. The small
differences
between equatorial and
polar neutrino energy losses, scaled by a factor of 10, are
plotted with
dash-dotted lines. Again we see that at large distances as
well as near the neutrinospheres the equatorial emission
exceeds slightly the polar one during the whole evolution after
ms post bounce, and
this excess is reduced by the neutrino production that occurs
in the cooling layer outside of a density of
1011 g cm-3. Separating the analysis for the
different neutrino kinds, it becomes obvious that only the
and
luminosities grow in the cooling layer
above the neutrinospheres, and that these neutrinos are more
abundantly produced in the polar regions where the loss of accretion
energy mainly occurs, as discussed in Sect. 3.2.
For electron neutrinos the emission in the accretion
layer near the poles is so strong that their polar energy
loss at large distances can become bigger than the equatorial
one. In contrast, the flux of heavy-lepton neutrinos, which
are radiated from the neutrinospheric region with a persistent
surplus of the equatorial energy loss, gets a bit reduced in
the accretion and postshock layers between
g cm-3 and r = 400 km, whereby
their equatorial emission excess grows to roughly twice the
value near the neutrinospheres. Their asymmetry
determines the final emission asymmetry, because they contribute
the major part (
56%) of the total neutrino energy loss
(in fact, the pole-equator asymmetries of
and
at large distances have opposite signs and nearly cancel
each other). These findings are in clear contradiction to the
results obtained by Kotake et al. (2007), who reported that the
neutrino gravitational-wave signals in their models were dominated
by the
emission and the other neutrinos did only contribute
very small amplitudes with the same signs (see their Fig. 3).
What is the reason that slightly (2-3%) more energy is radiated from the neutrinospheric layer near the equator than near the poles? We find that the emergence of this asymmetry coincides with the development of the stable convective shell inside the neutron star. This might suggest that convective transport near the equatorial plane is slightly more effective than in the polar regions. We have no satisfactory explanation for this phenomenon, because one would expect convection to operate similarly in all directions and not to produce a global asymmetry on angular scales much larger than the size of the (nonstationary) convective cells. This, however, is a statement that is correct in three dimensions but might not hold for convection in the 2D environment of our simulations. Here convective plumes and downflows are toroidal features because of the axial symmetry, and the polar grid has a preferred direction along the polar axis. We suspect - without having performed a detailed analysis - that the operation of convective overturn with cells of such a toroidal structure near and along the polar axis might be slightly different from the in- and outflows near the equatorial plane with fluid motions perpendicular to the axis. If this were indeed the origin of the neutrino emission asymmetry, it would imply a serious warning that the computation of such anisotropies for estimating gravitational-wave signals (and also pulsar kicks based on the recoil associated with anisotropic neutrino emission) requires 3D modeling, because in 2D the discussed effects may be just a consequence of the symmetry constraints.
But is this fishy directional asymmetry really responsible for
our prediction of the neutrino gravitational-wave signal? Fortunately,
it is not! To see this, we have evaluated our results in yet another
way. In the lower righthand panel of Fig. 11 we
display the cumulative - i.e. time-integrated from the beginning
until the end of our 2D simulation - neutrino energy loss (again
summed up for neutrinos and antineutrinos of all flavors) as a
function of radius for the polar, equatorial, and total emission.
The fourth (blue) curve represents the difference (
)
between equatorial
and polar terms, scaled by a factor of 20. Consistent with what
we discussed above there is an excess of the equatorial losses until
about 45 km. This radius encompasses well the convection zone
inside of the nascent neutron star during the whole post-bounce
evolution of the 2D model with L&S EoS
(Fig. 1) and is larger than the location
where
g cm-3 for post-bounce times
ms (lefthand panel of Fig. 4).
At greater radii when one enters the accretion and SASI layers,
however, the polar energy loss first dominates before the equatorial
one wins again at
km. After a local maximum at a
radius slightly larger than 100 km the equator-to-pole
asymmetry settles to its final value. This radial variation is
found to be qualitatively very similar for all kinds of neutrinos,
although with different values of the local maxima and minima and
of the final, asymptotic level at large radii.
From the lower righthand panel of
Fig. 11 we learn that the emission anisotropy
carried by the outgoing neutrinos is not imposed near the
neutrinospheres and preserved outside, but it is completely determined
by the accretion region behind the supernova shock and the SASI and
convective inhomogeneities in this region. The final emission anisotropy
just by chance gets close to the value near the neutrinospheres.
Moreover, the separate analysis for the different kinds of neutrinos
reveals that it is the muon and tau neutrinos that govern
the final anisotropy of the radiated energy flux. In the two
models presented here the energy carried away by these neutrinos
in the equatorial belt is slightly
higher than the one escaping from the polar cap areas.
This as well as the total asymmetry, however, is likely to be a
time- and model-dependent result and to
depend strongly on the detailed structure of the
inhomogeneities that evolve behind the deformed accretion shock.
It could also depend
on other model aspects such as the value of the mass accretion
rate, the direction-dependent stagnation radius of the shock,
or the relative sizes of the fluxes of the different neutrinos.
Such conjectures are not only supported by the H&W model in
comparison to the L&S run: the neutrino gravitational-wave
amplitude of the latter model exhibits a significantly clearer trend
in a monotonic decrease to negative values, while in the H&W case
one can see long phases where the amplitude or its time derivative
become positive (Fig. 9).
These conjectures are also in agreement with the results of other
models described in the Müller et al. (2004) paper, and, last
but not least, they are suggested by the
fact that the quantity
in the lower
righthand panel of Fig. 11 is a quantity with
large radial variations in both positive and negative directions.
It is rather unlikely that such a complex radial behavior can
be generic for models in two dimensions so that in all cases and
at all times the same effects are recovered that ultimately determine
the emission asymmetry in our models, namely (i) a delicate
almost-cancellation of the combined asymmetry of
and
due to the production and absorption processes of these
neutrinos in the neutrino-cooling and heating regions; and (ii) a
final dominance of the anisotropy of the radiated muon and tau
neutrinos as a consequence of their asymmetric transport through
the semi-transparent accretion layer, in which the total heavy-lepton
neutrino luminosity hardly changes but the equator-to-pole
asymmetry grows.
Our analysis therefore reveals that the gravitational-wave results for neutrinos depend on physical effects that are either not included or treated in a much simplified manner in previous works. It further reveals that a reliable and quantitatively meaningful computation of the tiny neutrino emission asymmetry is extremely difficult and does not only require a full transport treatment and good numerical resolution in all relevant regions, but ultimately must be done on the basis of 3D models.
4 Summary and conclusions
We have presented results of two stellar core-collapse simulations
performed with the P ROMETHEUS-V ERTEX code for a nonrotating
15
progenitor star in axial symmetry. One of these models
was computed with a soft version of the Lattimer & Swesty (1991)
nuclear equation of state, and the other one with the stiffer
Hillebrandt & Wolff (1985) EoS, which produces a significantly
less compact proto-neutron star. For these simulations we have
compared the neutrino and gravitational-wave signals between core
bounce and 400 ms later. In particular we were interested in
signal features that may yield information about the action of the
standing accretion shock instability (SASI) in the supernova core
during the long phase of post-bounce accretion that precedes the
onset of the supernova explosion in all recent multi-dimensional
calculations of stellar core collapse.
Our main results can be summarized as follows:
- 1)
- In both simulations we found that the weakened bounce shock produces a negative entropy gradient that partly overlaps with the negative lepton number gradient left behind by the escaping shock breakout burst of electron neutrinos. This region is Ledoux unstable and develops prompt post-shock convection, which leads to an early burst of gravitational-wave emission lasting several ten milliseconds after core bounce. Its frequency spectrum shows a prominent peak around 70-100 Hz.
- 2)
- A later, long-lasting phase of strong gravitational-wave emission
with a long-time trend in growing wave amplitudes sets in at
ms after bounce. It is produced by nonradial mass motions and anisotropic neutrino emission that are caused by convective overturn inside the proto-neutron star and in the neutrino-heating layer behind the stalled shock, and in particular by vigorous sloshing of the accretion shock due to the development of low-
mode SASI activity. The wave train for the matter signal consists of a superposition of quasi-periodic variations on timescales from a few milliseconds to several ten milliseconds. The spectra of the dominant
2 SASI modes peak between about 20 and roughly 100 Hz, and the gravitational-wave emission also shows significant power at frequencies up to
200 Hz. The most prominent maxima of the gravitational-wave spectra, however, are located at higher frequencies of 300-800 Hz and originate mainly from the outer layers of the neutron star where rapid gas flows are triggered by the supersonic impact and deceleration of accretion funnels.
- 3)
- The high-frequency peak of the gravitational-wave spectrum and
less clearly also the secondary maximum near the SASI frequency
exhibit a dependence on the
compactness of the nascent neutron star and thus on the properties
of the high-density EoS: the more compact neutron star leads to more
powerful shock oscillations and significantly larger gravitational-wave
amplitudes earlier after core bounce. Also the characteristic frequencies
of the low-
SASI modes and of the gravitational-wave signal are higher. While the main peak of the wave spectrum is located at 300-600 Hz for the stiffer EoS, it can be found at 600-800 Hz for the softer EoS with the more compact remnant.
- 4)
- The more compact neutron star also radiates neutrinos with higher
luminosities and greater mean energies. Anisotropic neutrino emission
produces a low-frequency component of the gravitational-wave signal
at
-200 Hz with an amplitude that dominates the matter signal by up to a factor of two. The emission anisotropy is mostly established in the cooling and heating regions of the accretion layer outside of the neutrinospheres and is determined by the anisotropic transport of muon and tau neutrinos through this layer, which leads to a slightly higher flux of muon and tau neutrinos near the equatorial plane of the polar grid, i.e., perpendicular to the direction of the SASI oscillations of the supernova shock. Despite the stronger emission of
and
from the accretion layer in the polar regions (i.e. in the direction of the SASI shock expansion and contraction), the gravitational-wave amplitude associated with the anisotropic energy loss in neutrinos exhibits a nearly monotonic trend to negative values.
- 5)
- The neutrino luminosities and mean energies, in particular the
ones that can be received by an observer in the polar (and thus
SASI) direction, show a quasi-periodic time variability with an
amplitude of several ten percent (up to about 50%) of the minimum
values for the luminosities and of roughly 1 MeV (up to 10%)
for the mean energies. The luminosity fluctuations are correlated
and in phase with the energy variations and are somewhat bigger
for
and
than for heavy-lepton neutrinos. They originate from alternating phases of gas accumulation and compression in the accretion layer in particular around the poles, caused by the expansion and contraction of the shock in the course of the SASI oscillations. The frequency spectra of the neutrino luminosities show most power between about 20 and 200 Hz and significant power in a decaying tail up to about 400 Hz as a consequence of emission variability associated with faster, convective modulations of the accretion flow between supernova shock and nascent neutron star. Since SASI and convective activity occur inseparably in the gain layer and stir the cooling layer below, power is found to be distributed over a wide range of frequencies in the Fourier spectra of the luminosities.
Moreover, nonradial mass motions such as convection and the SASI in the supernova core after bounce produce gravitational-wave emission with significant intensity. We found that the wave spectrum exhibits a broad hump between about 20 and 200 Hz and a main maximum with clearly more power at considerably higher frequencies. The latter is associated with SASI and accretion induced gas flows in the neutron star surface layers, where much more mass is involved than in the fairly dilute post-shock region. Both the primary and the secondary maximum of the spectrum are sensitive to the structure of the neutron star: larger gravitational-wave amplitudes and correspondingly higher spectral maxima as well as higher frequencies are obtained with a softer equation of state and more compact remnant. The relative SASI variations of the neutrino luminosities do not appear to be very sensitive to the nuclear equation of state, and a clear EoS-dependent difference in the modulation frequencies is also not visible, but the absolute values of the luminosities and mean energies of the neutrinos radiated by the more compact and hotter, accreting proto-neutron star in our simulations are 10-20% higher.
Anisotropic neutrino emission is also an important source of
gravitational waves with amplitudes that can dominate
the matter signal. The corresponding spectrum
possesses most power at frequencies below 100-200 Hz.
The emission anisotropy is a function with rapid short-time
variability (with periods of milliseconds to 10 ms)
and with an overlying, nearly monotonic long-time trend toward
slightly stronger emission from equatorial regions, i.e.,
perpendicular to the direction of the SASI shock oscillations.
This trend decides about the sign of the neutrino
gravitational-wave amplitude.
A detailed analysis of our models reveals that it is
the result of complex transport physics
in which the neutrinospheric flux, which has a low
pole-to-equator anisotropy, is modified strongly in the
inhomogeneous environment of the accretion layer
around the forming neutron star. The
pole-to-equator asymmetry changes with radius several times
in both directions and the
contributions from
,
,
and heavy-lepton
neutrinos partly cancel each other. A reasonable
determination of these transport asymmetries, which is not
only important for gravitational-wave calculations but also
for pulsar kick estimates based on neutrino recoil, requires
detailed transport calculations. Even our sophisticated
ray-by-ray approach may not be sufficient for accurate
predictions of measurable effects, but a full multi-angle
treatment of the transport may be necessary.
From this discussion it is clear that any reliable assessment of the detectability of SASI effects in the neutrino and gravitational-wave emission will ultimately require results from three-dimensional models. The growth of SASI modes with nonvanishing azimuthal order m was shown to be of potential importance, and the presence of such nonaxisymmetric asymmetries can have a big influence on the structure and the dynamics of the accretion layer (Blondin & Mezzacappa 2007; Yamasaki & Foglizzo 2008; Iwakami et al. 2008). It is probable that also the emission asymmetries and time modulations of the emitted signals exhibit significant differences compared to 2D models, in particular because in 3D the SASI axis is not forced to coincide with an axis of symmetry and with the axis of the numerical grid.
The SASI and convective activity, its duration as well as
strength, must be expected to vary with the progenitor star
and to depend on the length of the postbounce accretion phase
until the onset of the explosion (for hints to that, see the
results for 11.2
and 15
explosions in
Marek & Janka 2007). In the nonlinear stage, the SASI and
convective motions are of chaotic nature and therefore
stochastic differences also in the large-scale pattern and
global asymmetry are
likely to develop from small initial differences even in
progenitor stars (models) of the same kind (see, e.g., the
different explosion asymmetries obtained in a large set of
2D simulations performed by Scheck et al. 2006). Nevertheless,
basic features of the instabilites, in particular their mere
presence and long-time growth behavior and the sizable
modulation they impose on the
emitted signals over a period of possibly hundreds of
milliseconds, should be common to different realizations.
Judged from our present understanding, however, a
template-like uniformity of the signal structure with a
well defined frequency (or frequency range) and amplitude
that exhibit a regular dependence on the explosion
and progenitor properties, appears not to be very likely. More
model calculations are therefore needed to find out which
signal properties such as the modulation frequency,
are sufficiently generic to yield information about the
neutron star equation of state.
4.1 Detectability of SASI modulations
Should the SASI signatures and their magnitude discussed in this work be confirmed by 3D results, the question arises whether they will be measurable with existing or emerging neutrino and gravitational-wave experiments in the case of the next Galactic supernova. A detailed prediction of supernova neutrino signals in terrestrial detectors is a highly complex problem, because the neutrinos radiated by the nascent neutron star undergo flavor conversions so that the flavor arriving at the Earth is not necessarily the same as the one leaving the source. On the one hand, the flavor evolution is affected by neutrino-matter oscillations during the neutrino propagation through the density gradients in the supernova mantle and envelope (where the H and L resonances are encountered at densities around 103 g cm-3 and 10 g cm-3, respectively) and through the Earth. On the other hand, neutrino conversion is also induced by the nonlinear self-interaction associated with neutrino-neutrino forward scattering in the dense neutrino background close to the neutron star (for an application of these phenomena in the context of the neutrino signals from oxygen-neon-magnesium core supernovae, see Lunardini et al. 2008, and from black hole formation in stellar collapse events, see Nakazato et al. 2008, for a discussion of the flavor conversion physics, see the references to original publications in the latter two papers). The flavor composition of the neutrino burst in the detector and thus the detector response therefore do not only depend on the type and structure of the progenitor star but also on a number of degrees of freedom associated with unknown neutrino properties and the location of the experiment, e.g. the neutrino-mass hierarchy (normal or inverted), the still undetermined mixing angle between the mass eigenstates 1 and 3, and the nadir angle of the neutrino path through the Earth.
Because of the involved complexity we defer an investigation
of the flavor conversion effects and their consequences for the
neutrino measurements to a future paper. However, since the SASI
modulations affect the luminosities and mean energies of neutrinos
and antineutrinos of all flavors with similar strength, we can
ignore the flavor oscillation effects for a simple
estimate to convince us that the SASI modulations of the radiated
signal should be detectable in principle. Super-Kamiokande, for
example, is expected to capture roughly 8000 electron antineutrinos
through their absorption by protons, if a Galactic supernova
occurs at a distance of 10 kpc and releases at total energy
of
erg in
with an average energy
of 15 MeV. A luminosity of
erg s-1,
which is a typical value for the emission maxima during the
SASI phase, should account for an event rate of approximately 7000
captures per second, corresponding to about 70 events in the
10 ms intervals of peak emission.
If
is about 50%
lower in luminosity minima, about one half of this event rate
and event number can be expected during the low-emission periods.
Therefore the SASI modulations will lead to a variability
of the measured signal that can be larger, though not much, than
the statistical
fluctuations,
being the number of
detected neutrinos in a time interval of 10 ms.
For Hyper-Kamiokande, whose fiducial
volume would be roughly 30 times bigger than that of
Super-Kamiokande, the Poisson noise would be reduced to a few
percent of the event number within 10 ms, so that the
SASI modulations should create a clear
signature. Even more powerful for this purpose is the
IceCube experiment. It will be able to identify the supernova
neutrino burst by the Cherenkov glow caused
by the neutrino energy deposition in the ice, which will show
up as time-correlated noise in all phototubes. For a supernova
at a distance of 10 kpc, the
luminosity maxima
during the SASI phase are expected to cause a photon count
rate of more than 106 per second in the 4800 optical
modules (see Dighe et al. 2003).
The
events within 10 ms intervals and
the variations of several 1000 events between periods of
luminosity peaks and minima (in the case of IceCube
these variations
scale with
,
because the crucial quantity
is the energy deposition in the detector)
dominate by far the Poisson fluctuations associated
with the background counting rate of 300 Hz per optical
module; the Poisson variability of the background
corresponds to
photon events in
all phototubes, which is
1% of the supernova signal in
a 10 ms interval (an impression of the detection capability
even for a supernova model with a lower fluctuation
amplitude of the neutrino emission than in the models
discussed here, can be obtained from
the upper lefthand panel of Fig. 2 in Dighe et al. 2003).
The detectability of gravitational-wave signals very similar to
the ones presented in this work was investigated in much
detail by Müller et al. (2004). Due to the restriction
to a 90-degree grid (one hemisphere) in the case of the
15
core-collapse model s15r considered in the
latter paper, the contributions of
dipolar SASI modes were not included in there.
Moreover, since the collapse evolution was followed to later
times in our most recent simulations (about 400 ms instead
of
270 ms after bounce) and the quadrupole amplitudes
exhibit a trend in rise during the later post-bounce phases,
the signal-to-noise (S/N) ratios and maximum detection distances
given for LIGO I and Advanced LIGO by Müller et al. (2004)
are likely to be only lower limits. Since the S/N ratio scales
roughly linearly with the Fourier transform of the wave
amplitude, the detection prospects grow with longer integration
periods of the signal and thus higher Fourier power. The
measurement of dominant gravitational-wave power at frequencies
of several hundred Hz with growing strength over a period of hundreds
of milliseconds would therefore be a strong indication of SASI and
convective activity around the nascent neutron star. Müller
et al. (2004) estimated a S/N ratio of roughly 70 in Advanced
LIGO for a 15
supernova at 10 kpc, corresponding
to the possibility of capturing gravitational waves from such
sources up to a distance of roughly 100 kpc. A third-generation
underground interferometer facility (Einstein Telescope) would
lead to a gain in sensitivity of approximately a factor of 20
in the most relevant frequency range between 100 Hz
and 1000 Hz (Spiering 2007). This would increase the
expected detection rate of supernovae by up to a factor
of 3 (Ando et al. 2005). A detailed analysis of the
gravitational-wave signals and their observational implications
from a greater sample of core-collapse models will be presented
in a forthcoming paper.
Acknowledgements
We thank C.D. Ott and G.G. Raffelt for useful discussions. The project was supported by the Deutsche Forschungsgemeinschaft through the Transregional Collaborative Research Centers SFB/TR 27 ``Neutrinos and Beyond'' and SFB/TR 7 ``Gravitational Wave Astronomy'', and the Cluster of Excellence EXC 153 ``Origin and Structure of the Universe'' (http://www.universe-cluster.de). Computer time grants at the John von Neumann Institute for Computing (NIC) in Jülich, the Höchstleistungsrechenzentrum of the Stuttgart University (HLRS) under grant number SuperN/12758, the Leibniz-Rechenzentrum München, and the RZG in Garching are acknowledged.
References
- Alexeyev, E. N., Alexeyeva, L. N., Krivosheina, I. V., & Volchenko, V. I. 1988, Phys. Lett., B205, 209 [NASA ADS] (In the text)
- Ando, S., Beacom, J. F., & Yüksel, H. 2005, , 95, 171101 (In the text)
- Bionta, R. M., Blewitt, G., Bratton, C. B., et al. 1987, , 58, 1494 (IMB collaboration) (In the text)
- Blanchet, L., Damour, T., & Schäfer, G. 1990, MNRAS, 242, 289 [NASA ADS] (In the text)
- Blondin, J. M., & Mezzacappa, A. 2006, ApJ, 642, 401 [NASA ADS] [CrossRef] (In the text)
- Blondin, J. M., & Mezzacappa, A. 2007, Nature, 445, 58 [NASA ADS] [CrossRef] (In the text)
- Blondin, J. M., Mezzacappa, A., & DeMarino, C. 2003, ApJ, 584, 971 [NASA ADS] [CrossRef] (In the text)
- Bruenn, S. W., Dirk, C. J., Mezzacappa, A., et al. 2006, in SciDAC, Scientific Discovery through Advanced Computing, Denver, Colorado, USA, 25-29 June 2006, ed. W. M. Tang, et al., Journ. Phys. Conf. Ser., 46, 393 [arXiv:0709.0537] (In the text)
- Buras, R., Rampp, M., Janka, H.-Th., & Kifonidis, K. 2006a, A&A, 447, 1049 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Buras, R., Janka, H.-Th., Rampp, M., & Kifonidis, K. 2006b, A&A, 457, 281 [NASA ADS] [CrossRef] [EDP Sciences]
- Burrows, A. 1987, ApJ, 318, L57 [NASA ADS] [CrossRef] (In the text)
- Burrows, A. 1988, ApJ, 334, 891 [NASA ADS] [CrossRef] (In the text)
- Burrows, A., & Fryxell, B. A. 1993, ApJ, 418, L33 [NASA ADS] [CrossRef] (In the text)
- Burrows, A., Hayes, J., & Fryxell, B. A. 1995, ApJ, 450, 830 [NASA ADS] [CrossRef] (In the text)
- Burrows, A., Livne, E., Dessart, L., Ott, C. D., & Murphy, J. 2006, ApJ, 640, 878 [NASA ADS] [CrossRef] (In the text)
- Burrows, A., Livne, E., Dessart, L., Ott, C. D., & Murphy, J. 2007a, ApJ, 655, 416 [NASA ADS] [CrossRef] (In the text)
- Burrows, A., Dessart, L., Livne, E., Ott, C. D., & Murphy, J. 2007b, ApJ, 664, 416 [NASA ADS] [CrossRef]
- Dessart, L., Burrows, A., Livne, E., & Ott, C. D. 2006, ApJ, 645, 534 [NASA ADS] [CrossRef] (In the text)
- Dighe, A. S., Keil, M. Th., & Raffelt, G. G. 2003, JCAP, 6, 005 [NASA ADS] (In the text)
- Dimmelmeier, H., Ott, C. D., Marek, A., & Janka, H.-Th. 2008, Phys. Rev. D, 78, 064056 [NASA ADS] [CrossRef] (In the text)
- Epstein, R. 1978, ApJ, 223, 1037 [NASA ADS] [CrossRef] (In the text)
- Epstein, R. 1979, MNRAS, 188, 305 [NASA ADS] (In the text)
- Fischer, T., Gögelein, P., Liebendörfer, M., Mezzacappa, A., & Thielemann, F.-K. 2008, in Origin of Matter and Evolution of Galaxies, ed. T. Suda, T. Nozawa, A. Ohnishi, K. Kato, M. Y. Fujimoto, T. Kajino, & S. Kubono, Amer. Inst. Phys. Conf. Proc., 1016, 277 (In the text)
- Foglizzo, T. 2001, A&A, 368, 311 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Foglizzo, T. 2002, A&A, 392, 353 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Foglizzo, T., Scheck, L., & Janka, H.-Th. 2006, ApJ, 652, 1436 [NASA ADS] [CrossRef] (In the text)
- Foglizzo, T., Galletti, P., Scheck, L., & Janka, H.-Th. 2007, ApJ, 654, 1006 [NASA ADS] [CrossRef] (In the text)
- Fryer, C. L., & Warren, M. S. 2002, ApJ, 574, L65 [NASA ADS] [CrossRef] (In the text)
- Fryer, C. L., & Warren, M. S. 2004, ApJ, 601, 391 [NASA ADS] [CrossRef]
- Fryer, C. L., Holz, D. E., & Hughes, S. A. 2002, ApJ, 565, 430 [NASA ADS] [CrossRef]
- Fryer, C. L., Holz, D. E., & Hughes, S. A. 2004, ApJ, 609, 288 [NASA ADS] [CrossRef] (In the text)
- Heger, A., Woosley, S. E., & Spruit, H. C. 2005, ApJ, 626, 350 [NASA ADS] [CrossRef] (In the text)
- Herant, M., Benz, W., Hix, W. R., Fryer, C. L., & Colgate, S. A. 1994, ApJ, 435, 339 [NASA ADS] [CrossRef] (In the text)
- Hillebrandt, W., & Wolff, R. G. 1985, in Nucleosynthesis: Challenges and New Developments, ed. W. D. Arnett, & J. W. Truran (Chicago: Univ. Chicago Press), 131 (In the text)
- Hillebrandt, W., Nomoto, K., & Wolff, R. G. 1984, A&A, 133, 175 [NASA ADS] (In the text)
- Hirata, K., Kajita, T., Koshiba, M., Nakahata, M., Oyama, Y., et al. 1987, , 58, 1490 (Kamiokande II collaboration) (In the text)
- Iwakami, W., Kotake, K., Ohnishi, N., Yamada, S., & Sawada, K. 2008, ApJ, 678, 1207 [NASA ADS] [CrossRef] (In the text)
- Janka, H.-Th., & Müller, E. 1996, A&A, 306, 167 [NASA ADS] (In the text)
- Janka, H.-Th., Buras, R., Kitaura, F. S., et al. 2005, Nucl. Phys., A758, 19c [NASA ADS] (In the text)
- Janka, H.-Th., Langanke, K., Marek, A., Martínez-Pinedo, G., & Müller, B. 2007, Phys. Rep., 442, 38 [NASA ADS] [CrossRef] (In the text)
- Keil, W., & Janka, H.-Th. 1995, A&A, 296, 145 [NASA ADS] (In the text)
- Keil, W., Janka, H.-Th., & Müller, E. 1996, ApJ, 473, L111 [NASA ADS] [CrossRef] (In the text)
- Keil, W., Raffelt, G., & Janka, H.-Th. 2003, ApJ, 590, 971 [NASA ADS] [CrossRef] (In the text)
- Kotake, K., Sato, K., & Takahashi, K. 2006, Rep. Prog. Phys., 69, 971 [NASA ADS] [CrossRef] (In the text)
- Kotake, K., Ohnishi, N., & Yamada, S. 2007, ApJ, 655, 406 [NASA ADS] [CrossRef] (In the text)
- Lattimer, J. M., & Swesty, F. D. 1991, Nucl. Phys., A535, 331 [NASA ADS] (In the text)
- Liebendörfer, M., Mezzacappa, A., Messer, O. E. B., et al. 2003, Nucl. Phys., A719, 144c [NASA ADS] (In the text)
- Lunardini, C., Müller, B., & Janka, H.-Th. 2008, Phys. Rev. D, 78, 023016 [NASA ADS] [CrossRef] (In the text)
- Marek, A. 2007, Ph.D. Thesis, Technische Universität München (In the text)
- Marek, A., & Janka, H.-Th. 2007, ApJ, in press [arXiv:0708.3372]
- Marek, A., Dimmelmeier, H., Janka, H.-Th., Müller, E., & Buras, R. 2006, A&A, 445, 273 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Müller, E., & Janka, H.-Th. 1997, A&A, 317, 140 [NASA ADS] (In the text)
- Müller, E., Rampp, M., Buras, R., Janka, H.-Th., & Shoemaker, D. H. 2004, ApJ, 603, 221 [NASA ADS] [CrossRef] (In the text)
- Murphy, J. W., & Burrows, A. 2008, ApJ, 688, 1159 [NASA ADS] [CrossRef] (In the text)
- Nakazato, K., Sumiyoshi, K., Suzuki, H., & Yamada, S. 2008, Phys. Rev. D, 78, 083014 [NASA ADS] [CrossRef] (In the text)
- Ohnishi, N., Kotake, K., & Yamada, S. 2006, ApJ, 641, 1018 [NASA ADS] [CrossRef] (In the text)
- Ott, C. D. 2008, Classical and Quantum Gravity, submitted [arXiv:0809.0695] (In the text)
- Ott, C. D., Burrows, A., Livne, E., & Walder, R. 2004, ApJ, 600, 834 [NASA ADS] [CrossRef] (In the text)
- Ott, C. D., Burrows, A., Dessart, L., & Livne, E. 2006, , 96, 201102 (In the text)
- Ott, C. D., Burrows, A., Dessart, L., & Livne, E. 2008, ApJ, 685, 1069 [NASA ADS] [CrossRef]
- Raffelt, G. G. 2001, ApJ, 561, 890 [NASA ADS] [CrossRef] (In the text)
- Rampp, M., & Janka, H.-Th. 2002, A&A, 396, 361 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Sagert, I., Hempel, M., Pagliara, G., et al. 2008, , submitted (In the text)
- Scheck, L., Plewa, T., Janka, H.-Th., Kifonidis, K., & Müller, E. 2004, , 92, 011103 (In the text)
- Scheck, L., Kifonidis, K., Janka, H.-Th., & Müller, E. 2006, A&A, 457, 963 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Scheck, L., Janka, H.-Th., Foglizzo, T., & Kifonidis, K. 2008, A&A, 477, 931 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Shen, H., Toki, H., Oyamatsu, K., & Sumiyoshi, K. 1998, Nucl. Phys., A637, 435 [NASA ADS] (In the text)
- Spiering, C. 2007, The Messenger, 129, 33 [NASA ADS] (In the text)
- Sumiyoshi, K., Yamada, S., Suzuki, H., & Chiba, S. 2006, , 97, 091101 (In the text)
- Sumiyoshi, K., Yamada, S., & Suzuki, H. 2007, ApJ, 667, 382 [NASA ADS] [CrossRef] (In the text)
- Swesty, F. D., & Myra, E. S. 2005, J. Phys. Conf. Ser., 16, 380 [NASA ADS] [CrossRef] (In the text)
- Thompson, T. A., Quataert, E., & Burrows, A. 2005, ApJ, 620, 861 [NASA ADS] [CrossRef] (In the text)
- Yamasaki, T., & Foglizzo, T. 2008, ApJ, 679, 607 [NASA ADS] [CrossRef] (In the text)
- Woosley, S. E., & Bloom, J. S. 2006, ARA&A, 44, 507 [NASA ADS] [CrossRef] (In the text)
- Woosley, S. E., & Weaver, T. A. 1995, ApJS, 101, 181 [NASA ADS] [CrossRef] (In the text)
Footnotes
- ...
respectively
- Note that below a certain density, which is typically chosen to be 1011 g cm-3 after core bounce, we replace the high-density EoSs by an ideal-gas equation of state with electrons, positrons, photons, and a mixture of classical, nonrelativistic Boltzmann gases for nucleons, alpha particles, and 14 kinds of heavier nuclei (see Marek & Janka 2007, for more information).
- ...
- In order to save computer time, in particular
during the simulation phase right after core bounce when the timesteps
are constrained to very low values, the model with the H&W EoS was
performed with the assumption of equatorial symmetry until 125.3 ms
post bounce. Only afterwards it was continued with a full 180
grid. Tests showed that this had no important influence on the results, neither for the evolution of the shock radius, nor for the growth and development of hydrodynamic instabilities. In particular, we did not find any qualitative difference in the SASI modes, also for odd values of
, at later times (see Marek 2007).
- ... polynomials
- Note
that we have renormalized our amplitudes
compared to the coefficients
used by Kotake et al. (2007) in order to directly display in Fig. 5 the radius variations associated with the different spherical harmonics components. The coefficient A0,0 is identical with the angular average of the shock radius:
, whose values and evolution are similar to those of the maximum shock radius plotted in Fig. 4.
- ... integrand
- The factor
accounts for the projection of the surface elements on the radiating sphere perpendicular to the line connecting source and observer.
All Figures
![]() |
Figure 1:
Regions of mass motions due to hydrodynamic instabilities
in the supernova core. The plots show color coded
the lateral velocity ( top) and the radial velocity ( bottom)
in the equatorial plane of the polar grid
as functions of post-bounce time and radius for the 2D simulation with
the L&S EoS ( left) and for the simulation with the H&W EoS
( right). The displayed range of velocity values is limited
to |
Open with DEXTER | |
In the text |
![]() |
Figure 2:
Entropy (black) and electron fraction (red) profiles (the
plotted values correspond to
|
Open with DEXTER | |
In the text |
![]() |
Figure 3:
Four representative snapshots from the 2D simulation with
the L&S EoS at post-bounce times of 247 ms ( top left),
255 ms ( top right), 322 ms ( bottom left),
and 375 ms ( bottom right). The lefthand panel of each figure
shows color-coded the entropy distribution, the righthand panel the
radial velocity component with white and whitish hues denoting
matter at or near rest; black arrows in the righthand panel indicate
the direction of the velocity field in the post-shock region
(arrows were plotted only in regions where the absolute values of
the velocities were less than
|
Open with DEXTER | |
In the text |
![]() |
Figure 4: Left: maximum shock radii (solid lines) and proto-neutron star radii (dashed lines) as functions of post-bounce time for the 2D simulations with different nuclear equations of state. The neutron star radii are determined as the locations where the rest-mass density is equal to 1011 g cm-3. Right: shock contours at the different post-bounce times listed in the figure. The vertical axis of the plot is the symmetry axis of the simulation. |
Open with DEXTER | |
In the text |
![]() |
Figure 5:
Left: amplitudes for the decomposition of the
angle-dependent position of the
shock front,
|
Open with DEXTER | |
In the text |
![]() |
Figure 6:
Isotropic equivalent luminosities of electron neutrinos
( top), electron antineutrinos ( middle), and one kind of heavy-lepton
neutrinos ( |
Open with DEXTER | |
In the text |
![]() |
Figure 7: Mean energies of radiated neutrinos as functions of post-bounce time for our 1D simulations ( top) and 2D models ( middle and bottom) with both equations of state (the lefthand panels are for the L&S EoS, the right ones for the H&W EoS). The displayed data are defined as ratios of the energy flux to the number flux and correspond to the luminosities plotted with dashed and solid lines in Fig. 6. The panels in the middle show results for a lateral grid zone near the north polar axis, the bottom panels provide results that are averaged over the whole northern hemisphere of the computational grid. In all cases the evaluation has been performed in the laboratory frame at a distance of 400 km from the stellar center. |
Open with DEXTER | |
In the text |
![]() |
Figure 8: Fourier amplitudes of the time-dependent neutrino luminosities shown in Fig. 6 (the left/right column there corresponds to the left/right column in the present figure). In order to reduce the level of background noise, the evaluation was performed for post-bounce times t > 100 ms, for which the SASI modulation is clearly visible in the neutrino luminosities, and the signals were split in 17 overlapping segments of 128 ms, advanced in steps of 10 ms. The Fourier spectra of the 1D models show a weak contribution from low-amplitude numerical fluctuations of the transport results. |
Open with DEXTER | |
In the text |
![]() |
Figure 9:
Gravitational-wave quadrupole amplitudes
AE220as functions of post-bounce time ( left) and
corresponding Fourier spectra ( right) associated with mass
motions ( top) and anisotropic neutrino (
|
Open with DEXTER | |
In the text |
![]() |
Figure 10: Fourier spectra of the gravitational wave quadrupole amplitudes of the mass motions caused by prompt post-shock convection during the first 100 ms after core bounce (see upper lefthand panel of Fig. 9). |
Open with DEXTER | |
In the text |
![]() |
Figure 11:
Quantities that are relevant for the gravitational-wave
signal associated with the anisotropic neutrino emission in our 2D model
using the L&S EoS.
Top left: gravitational-wave amplitudes as functions of time
after core bounce. The black curve was obtained by considering the
neutrino radiation at a large distance from the source (identical with
the corresponding line in
the lower lefthand panel of Fig. 9, where, as in
Figs. 6 and 7, the evaluation
was performed in the laboratory frame at a radius of 400 km).
The red line shows the wave amplitude computed from the neutrino
luminosities that leave the neutrinospheric region
(rather, a sphere where the density is 1011 g cm-3).
The green line displays the difference between the black and red curves.
Top right: anisotropy parameter
|
Open with DEXTER | |
In the text |
Copyright ESO 2009
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.