Normalized to: Schure, K.
[3]
oai:arXiv.org:1310.7027 [pdf] - 1180250
From cosmic ray source to the Galactic pool
Submitted: 2013-10-25
The Galactic cosmic ray spectrum is a remarkably straight power law. Our
current understanding is that the dominant sources that accelerate cosmic rays
up to the knee ($3 \times 10^{15}$ eV) or perhaps even the ankle ($3 \times
10^{18}$ eV), are young Galactic supernova remnants. In theory, however, there
are various reasons why the spectrum may be different for different sources,
and may not even be a power law if nonlinear shock acceleration applies during
the most efficient stages of acceleration. We show how the spectrum at the
accelerator translates to the spectrum that make up the escaping cosmic rays
that replenish the Galactic pool of cosmic rays. We assume that cosmic ray
confinement, and thus escape, is linked to the level of magnetic field
amplification, and that the magnetic field is amplified by streaming cosmic
rays according to the non-resonant hybrid or resonant instability. When a fixed
fraction of the energy is transferred to cosmic rays, it turns out that a
source spectrum that is flatter than $E^{-2}$ will result in a $E^{-2}$ escape
spectrum, whereas a steeper source spectrum will result in an escape spectrum
with equal steepening. This alleviates some of the concern that may arise from
expected flat or concave cosmic ray spectra associated with nonlinear shock
modification.
[4]
oai:arXiv.org:1008.4367 [pdf] - 721088
The relation between post-shock temperature, cosmic-ray pressure and
cosmic-ray escape for non-relativistic shocks
Submitted: 2010-08-25, last modified: 2013-09-18
Supernova remnants are thought to be the dominant source of Galactic cosmic
rays. This requires that at least 5% of the available energy is transferred to
cosmic rays, implying a high cosmic-ray pressure downstream of supernova
remnant shocks. Recently, it has been shown that the downstream temperature in
some remnants is low compared to the measured shock velocities, implying that
additional pressure support by accelerated particles is present.
Here we use a two-fluid thermodynamic approach to derive the relation between
post-shock fractional cosmic-ray pressure and post-shock temperature, assuming
no additional heating beyond adiabatic heating in the shock precursor and with
all non-adiabatic heating occurring at the subshock. The derived relations show
that a high fractional cosmic-ray pressure is only possible, if a substantial
fraction of the incoming energy flux escapes from the system. Recently a shock
velocity and a downstream proton temperature were measured for a shock in the
supernova remnant RCW 86. We apply the two-fluid solutions to these
measurements and find that the the downstream fractional cosmic-ray pressure is
at least 50% with a cosmic-ray energy flux escape of at least 20%. In general,
in order to have 5% of the supernova energy go into accelerating cosmic rays,
on average the post-shock cosmic-ray pressure needs to be 30% for an effective
cosmic-ray adiabatic index of 4/3.
[5]
oai:arXiv.org:1307.8123 [pdf] - 1173081
Modeling the interaction of thermonuclear supernova remnants with
circumstellar structures: The case of Tycho's supernova remnant
Submitted: 2013-07-30
The well-established Type Ia remnant of Tycho's supernova (SN 1572) reveals
discrepant ambient medium density estimates based on either the measured
dynamics or on the X-ray emission properties. This discrepancy can potentially
be solved by assuming that the supernova remnant (SNR) shock initially moved
through a stellar wind bubble, but is currently evolving in the uniform
interstellar medium with a relatively low density.
We investigate this scenario by combining hydrodynamical simulations of the
wind-loss phase and the supernova remnant evolution with a coupled X-ray
emission model, which includes non-equilibrium ionization. For the explosion
models we use the well-known W7 deflagration model and the delayed detonation
model that was previously shown to provide good fits to the X-ray emission of
Tycho's SNR.
Our simulations confirm that a uniform ambient density cannot simultaneously
reproduce the dynamical and X-ray emission properties of Tycho. In contrast,
models that considered that the remnant was evolving in a dense, but small,
wind bubble reproduce reasonably well both the measured X-ray emission spectrum
and the expansion parameter of Tycho's SNR. Finally, we discuss possible mass
loss scenarios in the context of single- and double-degenerate models which
possible could form such a small dense wind bubble.
[7]
oai:arXiv.org:1307.6575 [pdf] - 1172962
Cosmic ray acceleration in young supernova remnants
Submitted: 2013-07-24
We investigate the appearance of magnetic field amplification resulting from
a cosmic ray escape current in the context of supernova remnant shock waves.
The current is inversely proportional to the maximum energy of cosmic rays, and
is a strong function of the shock velocity. Depending on the evolution of the
shock wave, which is drastically different for different circumstellar
environments, the maximum energy of cosmic rays as required to generate enough
current to trigger the non-resonant hybrid instability that confines the cosmic
rays follows a different evolution and reaches different values. We find that
the best candidates to accelerate cosmic rays to ~few PeV energies are young
remnants in a dense environment, such as a red supergiant wind, as may be
applicable to Cassiopeia A. We also find that for a typical background magnetic
field strength of 5 microG the instability is quenched in about 1000 years,
making SN1006 just at the border of candidates for cosmic ray acceleration to
high energies.
[8]
oai:arXiv.org:1301.7264 [pdf] - 1159378
Cosmic ray acceleration and escape from supernova remnants
Submitted: 2013-01-30
Galactic cosmic ray (CR) acceleration to the knee in the spectrum at a few
PeV is only possible if the magnetic field ahead of a supernova remnant (SNR)
shock is strongly amplified by CR escaping the SNR. A model formulated in terms
of the electric charge carried by escaping CR predicts the maximum CR energy
and the energy spectrum of CR released into the surrounding medium. We find
that historical SNR such as Cas A, Tycho and Kepler may be expanding too slowly
to accelerate CR to the knee at the present time.
[9]
oai:arXiv.org:1209.3043 [pdf] - 563422
Magnetic field amplification by cosmic rays in supernova remnants
Submitted: 2012-09-13
Magnetic field amplification is needed to accelerate cosmic cays to PeV
energies in supernova remants. Escaping cosmic rays trigger a return current in
the plasma that drives a non-resonant hybrid instability. We run simulations in
which we represent the escaping cosmic rays with the plasma return current,
keeping the maximum cosmic ray energy fixed, and evaluate its effects on the
upstream medium. In addition to magnetic field amplification, density
perturbations arise that, when passing through the shock, further increase
amplification levels downstream. As the growth rate of the instability is most
rapid for the smaller scales, the resolution is a limiting factor in the
amplification that can be reached with these simulations.
[10]
oai:arXiv.org:1203.1637 [pdf] - 1117143
Diffusive shock acceleration and magnetic field amplification
Submitted: 2012-03-07
Diffusive shock acceleration is the theory of particle acceleration through
multiple shock crossings. In order for this process to proceed at a rate that
can be reconciled with observations of high-energy electrons in the vicinity of
the shock, and for cosmic rays protons to be accelerated to energies up to
observed galactic values, significant magnetic field amplification is required.
In this review we will discuss various theories on how magnetic field
amplification can proceed in the presence of a cosmic ray population. On both
short and long scales, cosmic ray streaming can induce instabilities that act
to amplify the magnetic field. Developments in this area that have occurred
over the past decade are the main focus of this paper.
[11]
oai:arXiv.org:1111.0904 [pdf] - 1579385
Progenitor's signatures in Type Ia supernova remnants
Submitted: 2011-11-03
The remnants of Type Ia supernovae can provide important clues about their
progenitor-histories. We discuss two well-observed supernova remnants (SNRs)
that are believed to result from a Type Ia SN and use various tools to shed
light on the possible progenitor history. We find that Kepler's SNR is
consistent with a symbiotic binary progenitor consisted of a white dwarf and an
AGB star. Our hydrosimulations can reproduce the observed kinematic and
morphological properties. For Tycho's remnant we use the characteristics of the
X-ray spectrum and the kinematics to show that the ejecta has likely interacted
with dense circumstellar gas.
[12]
oai:arXiv.org:1110.1096 [pdf] - 422173
Confining the high-energy cosmic rays
Submitted: 2011-10-05
Diffusive shock acceleration is the prime candidate for efficient
acceleration of cosmic rays. Galactic cosmic rays are believed to originate
predominantly from this process in supernova remnant shock waves. Confinement
of the cosmic rays in the shock region is key in making the mechanism
effective. It has been known that on small scales (smaller than the typical
gyroradius) high-amplitude non-resonant instabilities arise due to cosmic ray
streaming ahead of the shock. For the efficiency of scattering of the highest
energy cosmic rays it is of interest to determine the type of instabilities
that act on longer length scales, i.e. larger than the cosmic ray gyroradius.
We will present the results of our analysis of an instability that acts in this
regime and will discuss its driving mechanism and typical growth times.
[13]
oai:arXiv.org:1108.0582 [pdf] - 1083147
Cosmic ray acceleration at oblique shocks
Submitted: 2011-08-02
We show that the diffusion approximation breaks down for particle
acceleration at oblique shocks with velocities typical of young supernova
remnants. Higher order anisotropies flatten the spectral index at
quasi-parallel shocks and steepen the spectral index at quasi-perpendicular
shocks. We compare the theory with observed spectral indices.
[14]
oai:arXiv.org:1107.5817 [pdf] - 1083064
A long-wavelength instability involving the stress tensor
Submitted: 2011-07-28
Cosmic ray acceleration through first-order Fermi acceleration in a
collisionless plasma relies on efficient scattering off magnetic field
fluctuations. Scattering is most efficient for magnetic field fluctuations with
wavelengths on the order of the gyroradius of the particles. In order to
determine the highest energy to which cosmic rays can be accelerated, it is
important to understand the growth of the magnetic field on these large scales.
We derive the growth rate of the long-wavelength fluctuations in the linear
regime, using the kinetic equation coupled to Maxwell's equations for the
background plasma. The instability, driven by the cosmic ray current, acts on
large scales due to the stress tensor and efficient scattering on small scales,
and operates for both left- and right circular polarisations. This
long-wavelength instability is potentially important in determining the
acceleration efficiency and maximum energy of cosmic rays around shock waves
such as in supernova remnants.
[15]
oai:arXiv.org:1103.5487 [pdf] - 1053003
The imprint of a symbiotic binary progenitor on the properties of
Kepler's supernova remnant
Submitted: 2011-03-28
We present a model for the Type Ia supernova remnant (SNR) of SN 1604, also
known as Kepler's SNR. We find that its main features can be explained by a
progenitor model of a symbiotic binary consisting of a white dwarf and an AGB
donor star with an initial mass of 4-5 M_sun. The slow, nitrogen rich wind
emanating from the donor star has partially been accreted by the white dwarf,
but has also created a circumstellar bubble. Based on observational evidence,
we assume that the system moves with a velocity of 250 km/s. Due to the
systemic motion the interaction between the wind and the interstellar medium
has resulted in the formation of a bow shock, which can explain the presence of
a one-sided, nitrogen rich shell. We present two-dimensional hydrodynamical
simulations of both the shell formation and the SNR evolution. The SNR
simulations show good agreement with the observed kinematic and morphological
properties of Kepler's SNR. Specifically, the model reproduces the observed
expansion parameters (m=V/(R/t)) of m=0.35 in the north and m=0.6 in the south
of Kepler's SNR. We discuss the variations among our hydrodynamical simulations
in light of the observations, and show that part of the blast wave may have
traversed through the one-sided shell completely. The simulations suggest a
distance to Kepler's SNR of 6 kpc, or otherwise require that SN 1604 was a
sub-energetic Type Ia explosion. Finally, we discuss the possible implications
of our model for Type Ia supernovae and their remnants in general.
[16]
oai:arXiv.org:1103.3049 [pdf] - 332113
A more accurate numerical scheme for diffusive shock acceleration
Submitted: 2011-03-15
We present a more accurate numerical scheme for the calculation of diffusive
shock acceleration of cosmic rays using Stochastic Differential Equations. The
accuracy of this scheme is demonstrated using a simple analytical flow profile
that contains a shock of finite width and a varying diffusivity of the cosmic
rays, where the diffusivity decreases across the shock. We compare the results
for the slope of the momentum distribution with those obtained from a
perturbation analysis valid for finite but small shock width. These
calculations show that this scheme, although computationally more expensive,
provides a significantly better performance than the Cauchy-Euler type schemes
that were proposed earlier in the case where steep gradients in the cosmic ray
diffusivity occur. For constant diffusivity the proposed scheme gives similar
results as the Cauchy-Euler scheme.
[17]
oai:arXiv.org:1004.2766 [pdf] - 1026349
Time-dependent particle acceleration in supernova remnants in different
environments
Submitted: 2010-04-16
We simulate time-dependent particle acceleration in the blast wave of a young
supernova remnant (SNR), using a Monte Carlo approach for the diffusion and
acceleration of the particles, coupled to an MHD code. We calculate the
distribution function of the cosmic rays concurrently with the hydrodynamic
evolution of the SNR, and compare the results with those obtained using simple
steady-state models. The surrounding medium into which the supernova remnant
evolves turns out to be of great influence on the maximum energy to which
particles are accelerated. In particular, a shock going through a $\rho \propto
r^{-2}$ density profile causes acceleration to typically much higher energies
than a shock going through a medium with a homogeneous density profile. We find
systematic differences between steady-state analytical models and our
time-dependent calculation in terms of spectral slope, maximum energy, and the
shape of the cut-off of the particle spectrum at the highest energies. We also
find that, provided that the magnetic field at the reverse shock is
sufficiently strong to confine particles, cosmic rays can be easily
re-accelerated at the reverse shock.
[18]
oai:arXiv.org:0909.5204 [pdf] - 28896
A new radiative cooling curve based on an up to date plasma emission
code
Submitted: 2009-09-28, last modified: 2009-11-26
This work presents a new plasma cooling curve that is calculated using the
SPEX package. We compare our cooling rates to those in previous works, and
implement the new cooling function in the grid-adaptive framework `AMRVAC'.
Contributions to the cooling rate by the individual elements are given, to
allow for the creation of cooling curves tailored to specific abundance
requirements. In some situations, it is important to be able to include
radiative losses in the hydrodynamics. The enhanced compression ratio can
trigger instabilities (such as the Vishniac thin-shell instability) that would
otherwise be absent. For gas with temperatures below 10,000 K, the cooling time
becomes very long and does not affect the gas on the timescales that are
generally of interest for hydrodynamical simulations of circumstellar plasmas.
However, above this temperature, a significant fraction of the elements is
ionised, and the cooling rate increases by a factor 1000 relative to lower
temperature plasmas.
[19]
oai:arXiv.org:0905.1134 [pdf] - 24054
Evolution of Magnetic Fields and Cosmic Ray Acceleration in Supernova
Remnants
Submitted: 2009-05-07
Observations show that the magnetic field in young supernova remnants (SNRs)
is significantly stronger than can be expected from the compression of the
circumstellar medium (CSM) by a factor of four expected for strong blast waves.
Additionally, the polarization is mainly radial, which is also contrary to
expectation from compression of the CSM magnetic field. Cosmic rays (CRs) may
help to explain these two observed features. They can increase the compression
ratio to factors well over those of regular strong shocks by adding a
relativistic plasma component to the pressure, and by draining the shock of
energy when CRs escape from the region. The higher compression ratio will also
allow for the contact discontinuity, which is subject to the Rayleigh-Taylor
(R-T) instability, to reach much further out to the forward shock. This could
create a preferred radial polarization of the magnetic field. With an adaptive
mesh refinement MHD code (AMRVAC), we simulate the evolution of SNRs with three
different configurations of the initial CSM magnetic field, and look at two
different equations of state in order to look at the possible influence of a CR
plasma component. The spectrum of CRs can be simulated using test particles, of
which we also show some preliminary results that agree well with available
analytical solutions.
[20]
oai:arXiv.org:0810.5150 [pdf] - 17938
Evolution of Magnetic Fields in Supernova Remnants
Submitted: 2008-10-28
Supernova remnants (SNR) are now widely believed to be a source of cosmic
rays (CRs) up to an energy of 1 PeV. The magnetic fields required to accelerate
CRs to sufficiently high energies need to be much higher than can result from
compression of the circumstellar medium (CSM) by a factor 4, as is the case in
strong shocks. Non-thermal synchrotron maps of these regions indicate that
indeed the magnetic field is much stronger, and for young SNRs has a dominant
radial component while for old SNRs it is mainly toroidal. How these magnetic
fields get enhanced, or why the field orientation is mainly radial for young
remnants, is not yet fully understood. We use an adaptive mesh refinement MHD
code, AMRVAC, to simulate the evolution of supernova remnants and to see if we
can reproduce a mainly radial magnetic field in early stages of evolution. We
follow the evolution of the SNR with three different configurations of the
initial magnetic field in the CSM: an initially mainly toroidal field, a
turbulent magnetic field, and a field parallel to the symmetry axis. Although
for the latter two topologies a significant radial field component arises at
the contact discontinuity due to the Rayleigh-Taylor instability, no radial
component can be seen out to the forward shock. Ideal MHD appears not
sufficient to explain observations. Possibly a higher compression ratio and
additional turbulence due to dominant presence of CRs can help us to better
reproduce the observations in future studies.
[21]
oai:arXiv.org:0806.4617 [pdf] - 13997
Jets as diagnostics of the circumstellar medium and the explosion
energetics of supernovae: the case of Cas A
Submitted: 2008-06-27
We present hydrodynamical models for the Cassiopeia A (Cas A) supernova
remnant and its observed jet / counter-jet system. We include the evolution of
the progenitor's circumstellar medium, which is shaped by a slow red supergiant
wind that is followed by a fast Wolf-Rayet (WR) wind.
The main parameters of the simulations are the duration of the WR phase and
the jet energy. We find that the jet is destroyed if the WR phase is
sufficiently long and a massive circumstellar shell has formed. We therefore
conclude that the WR phase must have been short (a few thousand yr), if present
at all. Since the actual jet length of Cas A is not known we derive a lower
limit for the jet energy, which is ~10^{48} erg. We discuss the implications
for the progenitor of Cas A and the nature of its explosion.
