Non-Stationary Rayleigh-Taylor Instabilities in Pulsar Wind Interaction with a Supernova Shell

The Rayleigh-Taylor instability (RTI) plays an important role in the dynamics of several astronomical objects, in particular, in the supernovae (SN) evolution. In the present paper we examine the dynamics of a shell (representing a type II SN remnant) blown-up by a wind emitted by a central pulsar. The shell is accelerated by the pulsar wind and its inner surface experiences the RTI. We develop an analytical approach by using a specific transformation into the coordinate frame co-moving with the SN ejecta. We first derive a non-stationary spherically symmetric solution describing an expansion of a gas shell under the pressure of a central source (pulsar). Then, we analyze its 3D stability with respect to a small perturbation on the inner shell surface. The dispersion relation is derived in the co-moving reference frame. The growth rate of the perturbation is found and its temporal evolution is discussed. We compare our result with the previous published studies and apply it to the Crab nebula evolution.     

Systematics in the interpretation of aggregated neutrino flux limits and flavor ratios from gamma-ray bursts

Gamma-ray burst analyses at neutrino telescopes are typically based on diffuse or stacked (i.e., aggregated) neutrino fluxes, because the number of events expected from a single burst is small. The interpretation of aggregated flux limits implies new systematics not present for a single burst, such as by the integration over parameter distributions (diffuse fluxes), or by the low statistics in small burst samples (stacked fluxes). We simulate parameter distributions with a Monte Carlo method computing the spectra burst by burst, as compared to a conventional Monte Carlo integration. With this approach, we can predict the behavior of the flux in the diffuse limit as well as in low statistics stacking samples, such as used in recent IceCube data analyses. We also include the flavor composition at the detector (ratio between muon tracks and cascades) into our considerations. We demonstrate that the spectral features, such as a characteristic multi-peak structure coming from photohadronic interactions, flavor mixing, and magnetic field effects, are typically present even in diffuse neutrino fluxes if only the redshift distribution of the sources is considered, with z ? 1 dominating the neutrino flux. On the other hand, we show that variations of the Lorentz boost can only be interpreted in a model-dependent way, and can be used as a model discriminator. For example, we illustrate that the observation of spectral features in aggregated fluxes will disfavor the commonly used assumption that bursts with small Lorentz factors dominate the neutrino flux, whereas it will be consistent with the hypothesis that the bursts have similar properties in the comoving frame.     

Collective effects of stellar winds and unidentified gamma-ray sources

We study collective wind configurations produced by a number of massive stars, and obtain densities and expansion velocities of the stellar wind gas that is to be target, in this model, of hadronic interactions. We study the expected γ-ray emission from these regions, considering in an approximate way the effect of cosmic ray modulation. We compute secondary particle production (electrons from knock-on interactions and electrons and positrons from charged pion decay), and solve the loss equation with ionization, synchrotron, bremsstrahlung, inverse Compton, and expansion losses. We provide examples where configurations can produce sources for GLAST satellite, and the MAGIC, HESS, or VERITAS telescopes in non-uniform ways, i.e., with or without the corresponding counterparts. We show that in all cases we studied no EGRET source is expected.     

Population statistics study of radio and gamma-ray pulsars in the Galactic plane

We present results of our pulsar population synthesis of normal and millisecond pulsars in the Galactic plane. Over the past several years, a program has been developed to simulate pulsar birth, evolution and emission using Monte Carlo techniques. We have added to the program the capability to simulate millisecond pulsars, which are old, recycled pulsars with extremely short periods. We model the spatial distribution of the simulated pulsars by assuming that they start with a random kick velocity and then evolve through the Galactic potential. We use a polar cap/slot gap model for γ-ray emission from both millisecond and normal pulsars. From our studies of radio pulsars that have clearly identifiable core and cone components, in which we fit the polarization sweep as well as the pulse profiles in order to constrain the viewing geometry, we develop a model describing the ratio of radio core-to-cone peak fluxes. In this model, short period pulsars are more cone-dominated than in our previous studies. We present the preliminary results of our recent study and the implications for observing these pulsars with GLAST and AGILE.      

Very high energy γ-rays from galactic sources

The field of Very High Energy (VHE) gamma ray astronomy using the Atmospheric Cerenkov Technique has entered an interesting phase with detection of various galactic and extragalactic sources. Among galactic sources, only the Crab nebula has been established as a standard candle. Most observations on pulsars are in agreement as to the necessity for the GeV spectra to steepen at < 200 GeV. While the Imaging method for increase of sensitivity has been successful with many results, an alternate technique —Wavefront Sampling Technique- is also being used by an increasing number of experiments. The recently commissioned experiment at Pachmarhi (PACT) in India is presented as an example of this technique. Preliminary results from this experiment show detection of VHE γ-ray emission from (a) the Crab nebula at a high significance and (b) Crab and Geminga pulsars at > 1.5 TeV which could be the second component of the Outer Gap pulsar models.