Evolution of Magnetic Fields in Supernova Remnants

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.Comment: 5 pages, 3 figures. To appear in conference proceedings of "Magnetic Fields in the Universe II" (2008), RevMexA

Simulation and theory of fluid demixing and interfacial tension of mixtures of colloids and non-ideal polymers

An extension of the Asakura-Oosawa-Vrij model of hard sphere colloids and non-adsorbing polymers, that takes polymer non-ideality into account through a repulsive stepfunction pair potential between polymers, is studied with grand canonical Monte Carlo simulations and density functional theory. Simulation results validate previous theoretical findings for the shift of the bulk fluid demixing binodal upon increasing strength of polymer-polymer repulsion, promoting the tendency to mix. For increasing strength of the polymer-polymer repulsion, simulation and theory consistently predict the interfacial tension of the free colloidal liquid-gas interface to decrease significantly for fixed colloid density difference in the coexisting phases, and to increase for fixed polymer reservoir packing fraction.Comment: 10 pages, 4 figure

The relation between post-shock temperature, cosmic-ray pressure and cosmic-ray escape for non-relativistic shocks

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.Comment: 9 pages, 6 color figures. This is updated with a corrected figure 5a and 5b, reflecting an ApJ erratu

Jets as diagnostics of the circumstellar medium and the explosion energetics of supernovae: the case of Cas A

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.Comment: 9 pages, 5 figures, ApJ accepted. Version with high resolution figures available at http://www.phys.uu.nl/~schure/CasA_jet.pd

The loss-limited electron energy in SN 1006: effects of the shock velocity and of the diffusion process

The spectral shape of the synchrotron X-ray emission from SN 1006 reveals the fundamental role played by radiative losses in shaping the high-energy tail of the electron spectrum. We analyze data from the XMM-Newton SN 1006 Large Program and confirm that in both nonthermal limbs the loss-limited model correctly describes the observed spectra. We study the physical origin of the observed variations of the synchrotron cutoff energy across the shell. We investigate the role played by the shock velocity and by the electron gyrofactor. We found that the cutoff energy of the syncrotron X-ray emission reaches its maximum value in regions where the shock has experienced its highest average speed. This result is consistent with the loss-limited framework. We also find that the electron acceleration in both nonthermal limbs of SN 1006 proceeds close to the Bohm diffusion limit, the gyrofactor being in the range 1.5-4. We finally investigate possible explanations for the low values of cutoff energy measured in thermal limbs.Comment: Accepted for publication in Astronomische Nachrichten. Proceedings of the XMM-Newton Science Workshop 201

Progenitor's signatures in Type Ia supernova remnants

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.Comment: 4 pages, 9 figures, proceedings for IAU Symposium 281, Padova, July 201

Modeling the interaction of thermonuclear supernova remnants with circumstellar structures: The case of Tycho's supernova remnant

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.Comment: 12 pages, 7 figures, accepted for publication in MNRA

Capillary Waves in a Colloid-Polymer Interface

The structure and the statistical fluctuations of interfaces between coexisting phases in the Asakura-Oosawa (AO) model for a colloid--polymer mixture are analyzed by extensive Monte Carlo simulations. We make use of a recently developed grand canonical cluster move with an additional constraint stabilizing the existence of two interfaces in the (rectangular) box that is simulated. Choosing very large systems, of size LxLxD with L=60 and D=120, measured in units of the colloid radius, the spectrum of capillary wave-type interfacial excitations is analyzed in detail. The local position of the interface is defined in terms of a (local) Gibbs surface concept. For small wavevectors capillary wave theory is verified quantitatively, while for larger wavevectors pronounced deviations show up. For wavevectors that correspond to the typical distance between colloids in the colloid-rich phase, the interfacial fluctuations exhibit the same structure as observed in the bulk structure factor. When one analyzes the data in terms of the concept of a wavevector-dependent interfacial tension, a monotonous decrease of this quantity with increasing wavevector is found. Limitations of our analysis are critically discussed.Comment: 12 pages, 15 figure

Revealing the obscured supernova remnant Kes 32 with Chandra

I report here on the analysis and interpretation of a Chandra observation of the supernova remnant Kes 32. Kes 32 is rather weak in X-rays due to a large interstellar absorption, which is found to be ~4E22 cm^-2, larger than previously reported. Spectral analysis indicates that the ionization age of this object is very young, with n_e t ~ 4E9 cm^-3s, and a temperature of kT_e ~ 1 keV. The X-ray emission peaks at a smaller radius than in the radio. The low ionization age suggests that Kes 32 is a young remnant. However, a young age is in contradiction with the relatively large apparent size, which indicates an age of several thousand years, instead of a few hundred years. This problem is discussed in connection with Kes 32's unknown distance and its possible association with the Norma galactic arm.Comment: Accepted for publication in the Astrophysical Journal. 7 pages, 7 figure