A Case Study of Small Scale Structure Formation in 3D Supernova Simulations

It is suggested in observations of supernova remnants that a number of large- and small-scale structures form at various points in the explosion. Multidimensional modeling of core-collapse supernovae has been undertaken since SN1987A, and both simulations and observations suggest/show that Rayleigh-Taylor instabilities during the explosion is a main driver for the formation of structure in the remnants. We present a case study of structure formation in 3D in a \msol{15} supernova for different parameters. We investigate the effect of moderate asymmetries and different resolutions of the formation and morphology of the RT unstable region, and take first steps at determining typical physical quantities (size, composition) of arising clumps. We find that in this progenitor the major RT unstable region develops at the He/OC interface for all cases considered. The RT instabilities result in clumps that are overdense by 1-2 orders of magnitude with respect to the ambient gas, have size scales on the level of a few % of the remnant diameter, and are not diffused after the first $\sim30$ yrs of the remnant evolution, in the absence of a surrounding medium.Comment: 59 pages, 34 figure

Non-stationary Rayleigh-Taylor instability in supernovae ejecta

The Rayleigh-Taylor instability plays an important role in the dynamics of several astronomical objects, in particular, in supernovae (SN) evolution. In this paper we develop an analytical approach to study the stability analysis of spherical expansion of the SN ejecta by using a special transformation in the co-moving coordinate frame. We first study a non-stationary spherical expansion of a gas shell under the pressure of a central source. Then we analyze its stability with respect to a no radial, non spherically symmetric perturbation of the of the shell. We consider the case where the polytropic constant of the SN shell is $\gamma=5/3$ and we examine the evolution of a arbitrary shell perturbation. The dispersion relation is derived. The growth rate of the perturbation is found and its temporal and spatial evolution is discussed. The stability domain depends on the ejecta shell thickness, its acceleration, and the perturbation wavelength.Comment: 16 page

Mass-loaded spherical accretion flows

We have calculated the evolution of spherical accretion flows undergoing mass-loading from embedded clouds through either conduction or hydrodynamical ablation. We have observed the effect of varying the ratios of the mass-loading timescale and the cooling timescale to the ballistic crossing timescale through the mass-loading region. We have also varied the ratio of the potential energy of a particle injected into the flow near the outer region of mass-loading to the temperature at which a minimum occurs in the cooling curve. The two types of mass-loading produce qualitatively different types of behaviour in the accretion flow, since mass-loading through conduction requires the ambient gas to be hot, whereas mass ablation from clumps occurs throughout the flow. Higher ratios of injected to accreted mass typically occur with hydrodynamical ablation, in agreement with previous work on wind-blown bubbles and supernova remnants. We find that mass-loading damps the radiative overstability of such flows, in agreement with our earlier work. If the mass-loading is high enough it can stabilize the accretion shock at a constant radius, yielding an almost isothermal subsonic post-shock flow. Such solutions may be relevant to cooling flows onto massive galaxies. Mass-loading can also lead to the formation of isolated shells of high temperature material, separated by gas at cooler temperatures

INTEGRAL observations of TeV plerions

Amongst the sources seen in very high gamma-rays several are associated with Pulsar Wind Nebulae (``TeV plerions''). The study of hard X-ray/soft gamma-ray emission is providing an important insight into the energetic particle population present in these objects. The unpulsed emission from pulsar/pulsar wind nebula systems in the energy range accessible to the INTEGRAL satellite is mainly synchrotron emission from energetic and fast cooling electrons close to their acceleration site. Our analyses of public INTEGRAL data of known TeV plerions detected by ground based Cherenkov telescopes indicate a deeper link between these TeV plerions and INTEGRAL detected pulsar wind nebulae. The newly discovered TeV plerion in the northern wing of the Kookaburra region (G313.3+0.6 powered by the middle aged PSR J1420-6048) is found to have a previously unknown INTEGRAL counterpart which is besides the Vela pulsar the only middle aged pulsar detected with INTEGRAL. We do not find an INTEGRAL counterpart of the TeV plerion associated with the X-ray PWN ``Rabbit'' G313.3+0.1 which is possibly powered by a young pulsar.Comment: 4 pages, 6 figures, proceedings of conference "The Multi-Messenger Approach to High-Energy Gamma-ray Sources" Barcelona/Spain (2006

X-ray variation statistics and wind clumping in Vela X-1

We investigate the structure of the wind in the neutron star X-ray binary system Vela X-1 by analyzing its flaring behavior. Vela X-1 shows constant flaring, with some flares reaching fluxes of more than 3.0 Crab between 20-60 keV for several 100 seconds, while the average flux is around 250 mCrab. We analyzed all archival INTEGRAL data, calculating the brightness distribution in the 20-60 keV band, which, as we show, closely follows a log-normal distribution. Orbital resolved analysis shows that the structure is strongly variable, explainable by shocks and a fluctuating accretion wake. Analysis of RXTE ASM data suggests a strong orbital change of N_H. Accreted clump masses derived from the INTEGRAL data are on the order of 5 x 10^19 -10^21 g. We show that the lightcurve can be described with a model of multiplicative random numbers. In the course of the simulation we calculate the power spectral density of the system in the 20-100 keV energy band and show that it follows a red-noise power law. We suggest that a mixture of a clumpy wind, shocks, and turbulence can explain the measured mass distribution. As the recently discovered class of supergiant fast X-ray transients (SFXT) seems to show the same parameters for the wind, the link between persistent HMXB like Vela X-1 and SFXT is further strengthened.Comment: 8 pages, 6 figures, accepted for publication in A&

A cost-effective method to quantify biological surface sediment reworking

We propose a simple and inexpensive method to determine the rate and pattern of surface sediment reworking by benthic organisms. Unlike many existing methods commonly used in bioturbation studies, which usually require sediment sampling, our approach is fully non-destructive and is well suited for investigating non-cohesive fine sediments in streams and rivers. Optical tracer (e.g., luminophores or coloured sand) disappearance or appearance is assessed through time based on optical quantification of surfaces occupied by tracers. Data are used to calculate surface sediment reworking (SSR) coefficients depicting bioturbation intensities. Using this method, we evaluated reworking activity of stream organisms (three benthic invertebrates and a fish) in laboratory microcosms mimicking pool habitats or directly in the field within arenas set in depositional zones. Our method was sensitive enough to measure SSR as low as 0.2 cm2.d-1, such as triggered by intermediate density (774 m-2) of Gammarus fossarum (Amphipoda) in microcosms. In contrast, complex invertebrate community in the field and a fish (Barbatula barabatula) in laboratory microcosms were found to yield to excessively high SSR (>60 cm2.d-1). Lastly, we suggest that images acquired during experiments can be used for qualitative evaluation of species-specific effects on sediment distribution

8.4GHz VLBI observations of SN2004et in NGC6946

We report on 8.4GHz Very Long Baseline Interferometry (VLBI) observations of the type II-P supernova SN2004et in the spiral galaxy NGC 6946, made on 20 February 2005 (151 days after explosion). The Very Large Array (VLA) flux density was 1.23$\pm$0.07 mJy, corresponding to an isotropic luminosity at 8.4GHz of (4.45$\pm$0.3)$\times10^{25}$ erg s$^{-1}$ Hz$^{-1}$ and a brightness temperature of (1.3$\pm$0.3)$\times10^{8}$ K. We also provide an improved source position, accurate to about 0.5 mas in each coordinate. The VLBI image shows a clear asymmetry. From model fitting of the size of the radio emission, we estimate a minimum expansion velocity of 15,700$\pm$2,000 km s$^{-1}$. This velocity is more than twice the expected mean expansion velocity estimated from a synchrotron self-absorbed emission model, thus suggesting that synchrotron self-absorption is not relevant for this supernova. With the benefit of an optical spectrum obtained 12 days after explosion, we favor an emission model which consists of two hot spots on an underlying expanding shell of width comparable to that of SN 1993J.Comment: Accepted for publication in A&A (22/05/07

Modeling of the Vela complex including the Vela supernova remnant, the binary system gamma2 Velorum, and the Gum nebula

We study the geometry and dynamics of the Vela complex including the Vela supernova remnant (SNR), the binary system gamma2 Velorum and the Gum nebula. We show that the Vela SNR belongs to a subclass of non-Sedov adiabatic remnants in a cloudy interstellar medium (ISM), the dynamics of which is determined by the heating and evaporation of ISM clouds. We explain observable characteristics of the Vela SNR with a SN explosion with energy 1.4 x 10^50 ergs near the step-like boundary of the ISM with low intercloud densities (~ 10^{-3} cm^{-3}) and with a volume-averaged density of clouds evaporated by shock in the north-east (NE) part about four times higher than the one in the south-west (SW) part. The observed asymmetry between the NE and SW parts of the Vela SNR could be explained by the presence of a stellar wind bubble (SWB) blown by the nearest-to-the Earth Wolf-Rayet (WR) star in the gamma2 Velorum system. We show that the size and kinematics of gamma2 Velorum SWB agree with predictions of numerical calculations for the evolution of the SWB of M_ini = 35M* star. The low initial mass of the WR star in gamma2 Velorum implies that the luminosity of the nuclear line of 26Al, produced by gamma2 Velorum, is below the sensitivity of existing gamma-ray telescopes.Comment: 8 pages, 2 figures, accepted for publication in A&

WD + MS systems as the progenitor of SNe Ia

We show the initial and final parameter space for SNe Ia in a ($\log P^{\rm i}, M_{\rm 2}^{\rm i}$) plane and find that the positions of some famous recurrent novae, as well as a supersoft X-ray source (SSS), RX J0513.9-6951, are well explained by our model. The model can also explain the space velocity and mass of Tycho G, which is now suggested to be the companion star of Tycho's supernova. Our study indicates that the SSS, V Sge, might be the potential progenitor of supernovae like SN 2002ic if the delayed dynamical-instability model due to Han & Podsiadlowski (2006) is appropriate. Following the work of Meng, Chen & Han (2009), we found that the SD model (WD + MS) with an optically thick wind can explain the birth rate of supernovae like SN 2006X and reproduce the distribution of the color excess of SNe Ia. The model also predicts that at least 75% of all SNe Ia may show a polarization signal in their spectra.Comment: 6 pages, 2 figures, accepted for publication in Astrophysics & Space Science (Proceeding of the 4th Meeting on Hot Subdwarf Stars and Related Objects, edited by Zhanwen Han, Simon Jeffery & Philipp Podsiadlowski

Modeling of Photoionized Plasmas

In this paper I review the motivation and current status of modeling of plasmas exposed to strong radiation fields, as it applies to the study of cosmic X-ray sources. This includes some of the astrophysical issues which can be addressed, the ingredients for the models, the current computational tools, the limitations imposed by currently available atomic data, and the validity of some of the standard assumptions. I will also discuss ideas for the future: challenges associated with future missions, opportunities presented by improved computers, and goals for atomic data collection.Comment: 17 pages, 8 figures, to appear in the proceedings of Xray2010, Utrecht, the Netherlands, March 15-17 201