Protostellar collapse: A comparison between SPH and AMR calculations

The development of parallel supercomputers allows today the detailed study of the collapse and the fragmentation of prestellar cores with increasingly accurate numerical simulations. Thanks to the advances in sub-millimeter observations, a wide range of observed initial conditions enable us to study the different modes of low-mass star formation. The challenge for the simulations is to reproduce the observational results. Two main numerical methods, namely AMR and SPH, are widely used to simulate the collapse and the fragmentation of prestellar cores. We compare thoroughly these two methods within their standard framework. We use the AMR code RAMSES and the SPH code DRAGON. Our physical model is as simple as possible and consists of an isothermal sphere rotating around the z-axis. We first study the conservation of angular momentum as a function of the resolution. Then, we explore a wide range of simulation parameters to study the fragmentation of prestellar cores. There seems to be a convergence between the two methods, provided resolution in each case is sufficient. Resolution criteria adapted to our physical cases, in terms of resolution per Jeans mass, for an accurate description of the formation of protostellar cores are deduced from the present study. This convergence is encouraging for future work in simulations of low-mass star formation, providing the aforementioned criteria are fulfilled. Higher resolution figures can be downloaded at http://www-dapnia.cea.fr/Projets/COAST/paper_amrvssph.pdfComment: 16 pages, 16 figures, accepted for publication in A&

The supernova remnant CTB 37B and its associated magnetar CXOU J171405.7-381031: evidence for a magnetar-driven remnant

We discuss in this Letter the association of the candidate magnetar CXOU J171405.7-381031 with the supernova remnant CTB 37B. The recent detection of the period derivative of the object allowed an estimation of a young characteristic age of only $\sim 1000 yr$. This value is too small to be compatible even with the minimum radius of the remnant $\geq 10 pc$, the value corresponding to the {\it lower} limit of the estimated distance of $10.2 \pm 3.5 kpc$, unless the true distance happens to be even smaller than the lower limit. We argue that a consistent scenario for the remnant origin, in which the latter is powered by the energy injected by a young magnetar, is indeed more accurate to explain the young age, and points out to its non-standard (i.e. magnetar-driven) nature.Comment: 6 pp., 1 figure, to appear in RAA Letter

An adjustable law of motion for relativistic spherical shells

A classical and a relativistic law of motion for an advancing shell are deduced applying the thin layer approximation. A new parameter connected with the quantity of absorbed matter in the expansion is introduced; this allows of matching theory and observation.Comment: 15 pages, 10 figures and article in press; Central European Journal of Physics 201

The contribution of supernova remnants to the galactic cosmic ray spectrum

The supernova paradigm for the origin of galactic cosmic rays has been deeply affected by the development of the non-linear theory of particle acceleration at shock waves. Here we discuss the implications of applying such theory to the calculation of the spectrum of cosmic rays at Earth as accelerated in supernova remnants and propagating in the Galaxy. The spectrum is calculated taking into account the dynamical reaction of the accelerated particles on the shock, the generation of magnetic turbulence which enhances the scattering near the shock, and the dynamical reaction of the amplified field on the plasma. Most important, the spectrum of cosmic rays at Earth is calculated taking into account the flux of particles escaping from upstream during the Sedov-Taylor phase and the adiabatically decompressed particles confined in the expanding shell and escaping at later times. We show how the spectrum obtained in this way is well described by a power law in momentum with spectral index close to -4, despite the concave shape of the instantaneous spectra of accelerated particles. On the other hand we also show how the shape of the spectrum is sensible to details of the acceleration process and environment which are and will probably remain very poorly known.Comment: 19 pages, 8 figures, published version (references updated

Non-linear diffusive acceleration of heavy nuclei in supernova remnant shocks

We describe a semi-analytical approach to non-linear diffusive shock acceleration in the case in which nuclei other than protons are also accelerated. The structure of the shock is determined by the complex interplay of all nuclei, and in turn this shock structure determines the spectra of all components. The magnetic field amplification upstream is described as due to streaming instability of all nuclear species. The amplified magnetic field is then taken into account for its dynamical feedback on the shock structure as well as in terms of the induced modification of the velocity of the scattering centers that enters the particle transport equation. The spectra of accelerated particles are steep enough to be compared with observed cosmic ray spectra only if the magnetic field is sufficiently amplified and the scattering centers have high speed in the frame of the background plasma. We discuss the implications of this generalized approach on the structure of the knee in the all-particle cosmic ray spectrum, which we interpret as due to an increasingly heavier chemical composition above $10^{15}$eV. The effects of a non trivial chemical composition at the sources on the gamma ray emission from a supernova remnant when gamma rays are of hadronic origin are also discussed.Comment: 23 pages, 5 figures, minor changes to reflect the published versio

MHD Interaction of Pulsar Wind Nebulae with SNRs and the ISM

In the late 1960s the discovery of the Crab pulsar in its associated supernova remnant, launched a new field in supernova remnant research: the study of pulsar-driven or plerionic supernova remnants. In these type of remnants, the relativistic wind emitted by the pulsar, blows a pulsar wind nebula into the interior of its supernova remnant. Now, more then forty years after the discovery of the Crab pulsar, there are more then fifty plerionic supernova remnants known, due to the ever-increasing capacity of observational facilities. I will review our current understanding of the different evolutionary stages of a pulsar wind nebula as it is interacting with its associated supernova remnant.Therefore I will discuss both analytical and more recent numerical (M)HD models.The four main stages of a pulsar wind nebula are: the supersonic expansion stage, the reverse shock interaction stage, the subsonic expansion stage and ultimatelythe stage when the head of the bubble is bounded by a bow shock, due to the supersonic motion of the pulsar. Ultimately this pulsar wind nebula bow shock will break through its associated remnant, after which the pulsar-powered bow shock will interact directly with the interstellar medium. I will discuss recent numerical models from these type of pulsar wind nebulae and their morphology.Comment: 11 pages, 3 figures, Advances in Space Research, in pres

The Nucleosynthetic Imprint of 15-40 Solar Mass Primordial Supernovae on Metal-Poor Stars

The inclusion of rotationally-induced mixing in stellar evolution can alter the structure and composition of presupernova stars. We survey the effects of progenitor rotation on nucleosynthetic yields in Population III and II supernovae using the new adaptive mesh refinement (AMR) code CASTRO. We examine spherical explosions in 15, 25 and 40 solar mass stars at Z = 0 and 10^-4 solar metallicity with three explosion energies and two rotation rates. Rotation in the Z = 0 models resulted in primary nitrogen production and a stronger hydrogen burning shell which led all models to die as red supergiants. On the other hand, the Z=10^-4 solar metallicity models that included rotation ended their lives as compact blue stars. Because of their extended structure, the hydrodynamics favors more mixing and less fallback in the metal free stars than the Z = 10^-4 models. As expected, higher energy explosions produce more enrichment and less fallback than do lower energy explosions, and less massive stars produce more enrichment and leave behind smaller remnants than do more massive stars. We compare our nucleosynthetic yields to the chemical abundances in the three most iron-poor stars yet found and reproduce the abundance pattern of one, HE 0557-4840, with a zero metallicity 15 solar mass, 2.4 x 10^51 erg supernova. A Salpeter IMF averaged integration of our yields for Z=0 models with explosion energies of 2.4x10^51 ergs or less is in good agreement with the abundances observed in larger samples of extremely metal-poor stars, provided 15 solar mass stars are included. Since the abundance patterns of extremely metal-poor stars likely arise from a representative sample of progenitors, our yields suggest that low-mass supernovae contributed the bulk of the metals to the early universe.Comment: 16 pages, 11 figures; submitted to Ap

The Origin and Distribution of Cold Gas in the Halo of a Milky Way-Mass Galaxy

We analyze an adaptive mesh refinement hydrodynamic cosmological simulation of a Milky Way-sized galaxy to study the cold gas in the halo. HI observations of the Milky Way and other nearby spirals have revealed the presence of such gas in the form of clouds and other extended structures, which indicates on-going accretion. We use a high-resolution simulation (136-272 pc throughout) to study the distribution of cold gas in the halo, compare it with observations, and examine its origin. The amount (10^8 Msun in HI), covering fraction, and spatial distribution of the cold halo gas around the simulated galaxy at z=0 are consistent with existing observations. At z=0 the HI mass accretion rate onto the disk is 0.2 Msun/yr. We track the histories of the 20 satellites that are detected in HI in the redshift interval 0.5>z>0 and find that most of them are losing gas, with a median mass loss rate per satellite of 3.1 x 10^{-3} Msun/yr. This stripped gas is a significant component of the HI gas seen in the simulation. In addition, we see filamentary material coming into the halo from the IGM at all redshifts. Most of this gas does not make it directly to the disk, but part of the gas in these structures is able to cool and form clouds. The metallicity of the gas allows us to distinguish between filamentary flows and satellite gas. We find that the former accounts for at least 25-75% of the cold gas in the halo seen at any redshift analyzed here. Placing constraints on cloud formation mechanisms allows us to better understand how galaxies accrete gas and fuel star formation at z=0.Comment: 13 pages, 8 figures. Accepted for publication in Ap

The evolution of mass loaded supernova remnants: II. Temperature dependent mass injection rates

We investigate the evolution of spherically symmetric supernova remnants in which mass loading takes place due to conductively driven evaporation of embedded clouds. Numerical simulations reveal significant differences between the evolution of conductively mass loaded and the ablatively mass loaded remnants studied in Paper I. A main difference is the way in which conductive mass loading is extinguished at fairly early times, once the interior temperature of the remnant falls below ~ 107 K. Thus, at late times remnants that ablatively mass load are dominated by loaded mass and thermal energy, while those that conductively mass load are dominated by swept-up mass and kinetic energy. Simple approximations to the remnant evolution, complementary to those in Paper I, are given

The kinematics and chemical stratification of the Type Ia supernova remnant 0519-69.0

We present an analysis of the XMM-Newton and Chandra X-ray data of the young Type Ia supernova remnant 0519-69.0 in the Large Magellanic Cloud. We used data from both the Chandra ACIS and XMM-Newton EPIC-MOS instruments, and high resolution X-ray spectra obtained with the XMM-Newton Reflection Grating Spectrometer. The Chandra data show that there is a radial stratification of oxygen, intermediate mass elements and iron, with the emission from more massive elements more toward the center. Using a deprojection technique we measure a forward shock radius of 4.0(3) pc and a reverse shock radius of 2.7(4) pc. We took the observed stratification of the shocked ejecta into account in the modeling of the X-ray spectra with multi-component NEI models, with the components corresponding to layers dominated by one or two elements. An additional component was added in order to represent the ISM, which mostly contributed to the continuum emission. This model fits the data well, and was also employed to characterize the spectra of distinct regions extracted from the Chandra data. From our spectral analysis we find that the fractional masses of shocked ejecta for the most abundant elements are: M(O)=32%, M(Si/S)=7%/5%, M(Ar+Ca)=1%, and M(Fe) = 55%. From the continuum component we derive a circumstellar density of nH= 2.4(2)/cm^3. This density, together with the measurements of the forward and reverse shock radii suggest an age of 450+/-200 yr,somewhat lower than, but consistent with the estimate based on the optical light echo (600+/-200 yr). From the RGS spectra we measured a Doppler broadening of sigma=1873+/-50 km/s, from implying a forward shock velocity of vS = 2770+/-500 km/s. We discuss the results in the context of single degenerate explosion models, using semi-analytical and numerical modeling, and compare the characteristics of 0519-69.0 with those of other Type Ia supernova remnants.Comment: Astronomy and Astrophysics in press. This version is the A&A accepted version, which contains improved figures and an extended discussion sectio