Simulation of the growth of the 3D Rayleigh-Taylor instability in Supernova Remnants using an expanding reference frame

Context: The Rayleigh-Taylor instabilities generated by the deceleration of a supernova remnant during the ejecta-dominated phase are known to produce finger-like structures in the matter distribution which modify the geometry of the remnant. The morphology of supernova remnants is also expected to be modified when efficient particle acceleration occurs at their shocks. Aims: The impact of the Rayleigh-Taylor instabilities from the ejecta-dominated to the Sedov-Taylor phase is investigated over one octant of the supernova remnant. We also study the effect of efficient particle acceleration at the forward shock on the growth of the Rayleigh-Taylor instabilities. Methods: We modified the Adaptive Mesh Refinement code RAMSES to study with hydrodynamic numerical simulations the evolution of supernova remnants in the framework of an expanding reference frame. The adiabatic index of a relativistic gas between the forward shock and the contact discontinuity mimics the presence of accelerated particles. Results: The great advantage of the super-comoving coordinate system adopted here is that it minimizes numerical diffusion at the contact discontinuity, since it is stationary with respect to the grid. We propose an accurate expression for the growth of the Rayleigh-Taylor structures that connects smoothly the early growth to the asymptotic self-similar behaviour. Conclusions: The development of the Rayleigh-Taylor structures is affected, although not drastically, if the blast wave is dominated by cosmic rays. The amount of ejecta that makes it into the shocked interstellar medium is smaller in the latter case. If acceleration occurs at both shocks the extent of the Rayleigh-Taylor structures is similar but the reverse shock is strongly perturbed.Comment: 15 pages, 12 figures, accepted for publication in Astronomy and Astrophysics with minor editorial changes. Version with full resolution images can be found at http://www.lpl.arizona.edu/~ffrasche/~12692.pd

XMM-Newton observation of Kepler's supernova remnant

We present the first results coming from the observation of Kepler's supernova remnant obtained with the EPIC instruments on board the XMM-Newton satellite. We focus on the images and radial profiles of the emission lines (Si K, Fe L, Fe K) and of the high energy continuum. Chiefly, the Fe L and Si K emission-line images are generally consistent with each other and the radial profiles show that the Si K emission extends to a larger radius than the Fe L emission (distinctly in the southern part of the remnant). Therefore, in contrast to Cas A, no inversion of the Si- and Fe-rich ejecta layers is observed in Kepler. Moreover, the Fe K emission peaks at a smaller radius than the Fe L emission, which implies that the temperature increases inwards in the ejecta. The 4-6 keV high energy continuum map shows the same distribution as the asymmetric emission-line images except in the southeast where there is a strong additional emission. A two color image of the 4-6 keV and 8-10 keV high energy continuum illustrates that the hardness variations of the continuum are weak all along the remnant except in a few knots. The asymmetry in the Fe K emission-line is not associated with any asymmetry in the Fe K equivalent width map. The Si K maps lead to the same conclusions. Hence, abundance variations do not cause the north-south brightness asymmetry. The strong emission in the north may be due to overdensities in the circumstellar medium. In the southeastern region of the remnant, the lines have a very low equivalent width and the X-ray emission is largely nonthermal.Comment: 15 pages, 15 figures, accepted for publication in A&

A Chandra Survey of the Nearest ULIRGs: Obscured AGN or Super-Starbursts?

We present initial results from a Chandra survey of a complete sample of the 8 nearest (z <= 0.04) ultraluminous IR galaxies (ULIRGs), and also include the IR-luminous galaxy NGC 6240 for comparison. In this paper we use the hard X-rays (2-8 kev) to search for the possible presence of an obscured AGN. In every case, a hard X-ray source is detected in the nuclear region. If we divide the sample according to the optical/IR spectroscopic classification (starburst vs. AGN), we find that the 5 ``starburst'' ULIRGs have hard X-ray luminosities about an order-of-magnitude smaller than the 3 ``AGN'' ULIRGs. NGC 6240 has an anomalously high hard X-ray luminosity compared to the ``starburst'' ULIRGs. The Fe Kalpha line is convincingly detected in only two ULIRGs. The weakness of the Fe-K emission in these ULIRGs generally suggests that the hard X-ray spectrum is not dominated by reflection from high N_H neutral material. The hard X-ray continuum flux ranges from a few X 10^3 to a few X 10^-5 of the far-IR flux, similar to values in pure starbursts, and several orders-of-magnitude smaller than in Compton-thin AGN. The upper limits on the ratio of the Fe Kalpha to far-IR flux are below the values measured in Compton-thick type 2 Seyfert galaxies. While very large column densities of molecular gas are observed in the nuclei of these galaxies, we find no evidence that the observed X-ray sources are obscured by Compton-thick material. Thus, our new hard X-ray data do not provide direct evidence that powerful ``buried quasars'' dominate the overall energetics of most ultraluminous infrared galaxies.Comment: 22 pages, 8 figures, formatted with emulateapj.sty, accepted for publication in August 2003 Ap

On the spherical-axial transition in supernova remnants

A new law of motion for supernova remnant (SNR) which introduces the quantity of swept matter in the thin layer approximation is introduced. This new law of motion is tested on 10 years observations of SN1993J. The introduction of an exponential gradient in the surrounding medium allows to model an aspherical expansion. A weakly asymmetric SNR, SN1006, and a strongly asymmetric SNR, SN1987a, are modeled. In the case of SN1987a the three observed rings are simulated.Comment: 19 figures and 14 pages Accepted for publication in Astrophysics & Space Science in the year 201

Gas and Dark Matter Spherical Dynamics

We investigate the formation of spherical cosmological structures following both dark matter and gas components. We focus on the dynamical aspect of the collapse assuming an adiabatic, $\gamma = 5/3$, fully ionized primordial plasma. We use for that purpose a fully Lagrangian hydrodynamical code designed to describe highly compressible flows in spherical geometry. We investigate also a "fluid approach" to describe the mean physical quantities of the dark matter flow. We test its validity for a wide range of initial density contrast. We show that an homogeneous isentropic core forms in the gas distribution, surrounded by a self-similar hydrostatic halo, with much higher entropy generated by shock dissipation. We derive analytical expressions for the size, density and temperature of the core, as well as for the surrounding halo. We show that, unless very efficient heating processes occur in the intergalactic medium, we are unable to reproduce within adiabatic models the typical core sizes in X-ray clusters. We also show that, for dynamical reasons only, the gas distribution is naturally antibiased relative to the total mass distribution, without invoking any reheating processes. This could explain why the gas fraction increases with radius in very large X-ray clusters. As a preparation for the next study devoted to the thermodynamical aspect of the collapse, we investigate the initial entropy level required to solve the core problem in X-ray clusters.Comment: 26 pages, 5 figures, accepted for publication in The Astrophysical Journa

XMM-Newton observation of the Tycho Supernova Remnant

We present the observation of the Tycho supernova remnant obtained with the EPIC and RGS instruments onboard the XMM-Newton satellite. We compare images and azimuthally averaged radial profiles in emission lines from different elements (silicon and iron) and different transition lines of iron (Fe L and Fe K). While the Fe XVII L line and Si XIII K line images are globally spatially coincident, the Fe K emission clearly peaks at a smaller radius, indicating a higher temperature toward the reverse shock. This is qualitatively the profile expected when the reverse shock, after travelling through the outer power-law density profile, has entered the central plateau of the ejecta. The high energy continuum map has an overall smooth distribution, with a similar extent to the radio emission. Its radial profile peaks further out than the lines emission. Brighter and harder continuum regions are observed with a rough bipolar symmetry in the eastern and western edges. The spectral analysis of the southeastern knots supports spatial variations of the relative abundance of silicon and iron, which implies an incomplete mixing of the silicon and iron layers.Comment: 6 pages, 6 figures, accepted for publication in A&A Special Letters Issue on XMM-Newto

Cosmic ray diffusion near the Bohm limit in the Cassiopeia A supernova remnant

Supernova remnants (SNRs) are believed to be the primary location of the acceleration of Galactic cosmic rays, via diffusive shock (Fermi) acceleration. Despite considerable theoretical work the precise details are still unknown, in part because of the difficulty in directly observing nucleons that are accelerated to TeV energies in, and affect the structure of, the SNR shocks. However, for the last ten years, X-ray observatories ASCA, and more recently Chandra, XMM-Newton, and Suzaku have made it possible to image the synchrotron emission at keV energies produced by cosmic-ray electrons accelerated in the SNR shocks. In this article, we describe a spatially-resolved spectroscopic analysis of Chandra observations of the Galactic SNR Cassiopeia A to map the cutoff frequencies of electrons accelerated in the forward shock. We set upper limits on the electron diffusion coefficient and find locations where particles appear to be accelerated nearly as fast as theoretically possible (the Bohm limit).Comment: 18 pages, 5 figures. Accepted for publication in Nature Physics (DOI below), final version available week of August 28, 2006 at http://www.nature.com/nphy

The European Photon Imaging Camera on XMM-Newton: The MOS Cameras

The EPIC focal plane imaging spectrometers on XMM-Newton use CCDs to record the images and spectra of celestial X-ray sources focused by the three X-ray mirrors. There is one camera at the focus of each mirror; two of the cameras contain seven MOS CCDs, while the third uses twelve PN CCDs, defining a circular field of view of 30 arcmin diameter in each case. The CCDs were specially developed for EPIC, and combine high quality imaging with spectral resolution close to the Fano limit. A filter wheel carrying three kinds of X-ray transparent light blocking filter, a fully closed, and a fully open position, is fitted to each EPIC instrument. The CCDs are cooled passively and are under full closed loop thermal control. A radio-active source is fitted for internal calibration. Data are processed on-board to save telemetry by removing cosmic ray tracks, and generating X-ray event files; a variety of different instrument modes are available to increase the dynamic range of the instrument and to enable fast timing. The instruments were calibrated using laboratory X-ray beams, and synchrotron generated monochromatic X-ray beams before launch; in-orbit calibration makes use of a variety of celestial X-ray targets. The current calibration is better than 10% over the entire energy range of 0.2 to 10 keV. All three instruments survived launch and are performing nominally in orbit. In particular full field-of-view coverage is available, all electronic modes work, and the energy resolution is close to pre-launch values. Radiation damage is well within pre-launch predictions and does not yet impact on the energy resolution. The scientific results from EPIC amply fulfil pre-launch expectations.Comment: 9 pages, 11 figures, accepted for publication in the A&A Special Issue on XMM-Newto

Magnetic fields in supernova remnants and pulsar-wind nebulae

We review the observations of supernova remnants (SNRs) and pulsar-wind nebulae (PWNe) that give information on the strength and orientation of magnetic fields. Radio polarimetry gives the degree of order of magnetic fields, and the orientation of the ordered component. Many young shell supernova remnants show evidence for synchrotron X-ray emission. The spatial analysis of this emission suggests that magnetic fields are amplified by one to two orders of magnitude in strong shocks. Detection of several remnants in TeV gamma rays implies a lower limit on the magnetic-field strength (or a measurement, if the emission process is inverse-Compton upscattering of cosmic microwave background photons). Upper limits to GeV emission similarly provide lower limits on magnetic-field strengths. In the historical shell remnants, lower limits on B range from 25 to 1000 microGauss. Two remnants show variability of synchrotron X-ray emission with a timescale of years. If this timescale is the electron-acceleration or radiative loss timescale, magnetic fields of order 1 mG are also implied. In pulsar-wind nebulae, equipartition arguments and dynamical modeling can be used to infer magnetic-field strengths anywhere from about 5 microGauss to 1 mG. Polarized fractions are considerably higher than in SNRs, ranging to 50 or 60% in some cases; magnetic-field geometries often suggest a toroidal structure around the pulsar, but this is not universal. Viewing-angle effects undoubtedly play a role. MHD models of radio emission in shell SNRs show that different orientations of upstream magnetic field, and different assumptions about electron acceleration, predict different radio morphology. In the remnant of SN 1006, such comparisons imply a magnetic-field orientation connecting the bright limbs, with a non-negligible gradient of its strength across the remnant.Comment: 20 pages, 24 figures; to be published in SpSciRev. Minor wording change in Abstrac

No cold dust within the supernova remnant Cassiopeia A

A large amount (about three solar masses) of cold (18 K) dust in the prototypical type II supernova remnant Cassiopeia A was recently reported. It was concluded that dust production in type II supernovae can explain how the large quantities (10^8 solar masses) of dust observed in the most distant quasars could have been produced within only 700 million years after the Big Bang. Foreground clouds of interstellar material, however, complicate the interpretation of the earlier submillimetre observations of Cas A. Here we report far-infrared and molecular line observations that demonstrate that most of the detected submillimetre emission originates from interstellar dust in a molecular cloud complex located in the line of sight between the Earth and Cas A, and is therefore not associated with the remnant. The argument that type II supernovae produce copious amounts of dust is not supported by the case of Cas A, which previously appeared to provide the best evidence for this possibility.Comment: 13 pages, 3 figure