Search

Observations of several Fanaroff-Riley (FR) type I sources reveal outflowing bipolar bubbles of hot gas surrounded by a weak forward shock. We consider the possibility that these bubbles were driven by choked relativistic jets which failed to penetrate the ambient intracluster medium (ICM). Using new results on choked jets linking the geometry of the forward shock to the jet properties, we infer robust limits on the radius $R_{ch}$ at which the jet was quenched in 5 well-studied FRI sources, finding typically $R_{ch}\sim 10$ kpc. We further show that, in order to reach this radius in less than the current age of the system, the jet must have been tightly collimated, with the jet head subtending an angle of $\theta_h < 2^\circ$. The ambient pressure is not high enough to explain this collimation, suggesting that the jet was collimated by interaction with its own cocoon. Although the choking radius is well-constrained, we find a degeneracy between the initial jet opening angle before collimation, $\theta_0$, and the duration of jet activity, $t_b$, with $(t_b/1\rm{Myr})(\theta_0/5^\circ)^{-2} \sim 0.1$. We speculate that the working time and/or opening angle of the jet may be important factors contributing to the FR type I/type II morphology in galaxy clusters, with short-lived or wide jets being choked to form bipolar bubbles filled with diffuse radio emission, and longer-lived or narrow jets successfully escaping the cluster core to produce cocoons with radio hotspots., Comment: 12 pages, 6 figures, submitted to MNRAS

G21.5-0.9 is a plerionic supernova remnant (SNR) used as a calibration target for the Chandra X-ray telescope. The first observations found an extended halo surrounding the bright central pulsar wind nebula (PWN). A 2005 study discovered that this halo is limb-brightened and suggested the halo to be the missing SNR shell. In 2010 the spectrum of the limb-brightened shell was found to be dominated by non-thermal X-rays. In this study, we combine 15 years of Chandra observations comprising over 1~Msec of exposure time (796.1~ks with the Advanced CCD Imaging Spectrometer (ACIS) and 306.1~ks with the High Resolution Camera (HRC)) to provide the deepest-to-date imaging and spectroscopic study. The emission from the limb is primarily non-thermal and is described by a power-law model with a photon index $\Gamma = 2.22 \, (2.04-2.34)$, plus a weak thermal component characterized by a temperature $kT = 0.37\, (0.20-0.64)$ keV and a low ionization timescale of $n_{e}t < 2.95 \times 10^{10}$ cm$^{-3}$s. The northern knot located in the halo is best fitted with a two-component power-law + non-equilibrium ionization thermal model characterized by a temperature of 0.14 keV and an enhanced abundance of silicon, confirming its nature as ejecta. We revisit the spatially resolved spectral study of the PWN and find that its radial spectral profile can be explained by diffusion models. The best fit diffusion coefficient is $D \sim 2.1\times 10^{27}\rm cm^2/s$ assuming a magnetic field $B =130 \mu G$, which is consistent with recent 3D MHD simulation results., Comment: 13 pages, 14 figures, accepted to be published in MNRAS, updated to match published version

Nearby pulsars have been suggested as sources of $\sim$TeV $e^+/e^-$ Cosmic Ray (CR) excess on Earth. The High-Altitude Water Cherenkov Observatory (HAWC) detected extended TeV emission regions in the direction of two nearby middle-aged pulsars, Geminga and PSR B0656+14. By modeling the TeV emission as inverse Compton emission from electron-positron pairs diffusing in the interstellar medium (ISM), the HAWC collaboration derives a diffusion coefficient much smaller than the standard value in the vicinity of the two pulsars, which make them unlikely the origin of the positron excess. We propose that the observed $\gamma$-ray emission originate from the relic pulsar wind nebula. A two zone diffusion model with a slow diffusion in the nebula and a fast diffusion in the ISM can explain the HAWC surface brightness profile and the positron excess simultaneously. Inefficient diffusion in the $\gamma$-ray emission region surrounding a middle-aged pulsar maybe a common phenomenon that can be tested by future observation. The implied diffusion coefficient in the ISM is smaller than the one suggested by the standard CR propagation model, but it is fully consistent with the predictions of the spiral arm model., Comment: 12 pages, 9 figures and 4 tables, submitted to MNRAS

We report on a broadband study of a complex X-ray source (1SAX J0618.0+2227) associated with the interaction site of the supernova remnant (SNR) IC 443 and ambient molecular cloud (MC) using NuSTAR, XMM_Newton, and Chandra observations. Its X-ray spectrum is composed of both thermal and non-thermal components. The thermal component can be equally well represented by either a thin plasma model with kT=0.19 keV or a blackbody model with kT=0.11 keV. The non-thermal component can be fit with either a power-law with Gamma~1.7 or a cutoff power-law with Gamma~1.5 and a cutoff energy at E_cut~18 keV. Using the newly obtained NuSTAR dataset, we test three possible scenarios for isolated X-ray sources in the SNR-MC interaction site: 1) pulsar wind nebula (PWN); 2) SNR ejecta fragment; 3) shocked molecular clump. We conclude that this source is most likely composed of a SNR ejecta (or a PWN) and surrounding shocked molecular clumps. The nature of this hard X-ray source in the SNR-MC interaction site of IC 443 may shed light on unidentified X-ray sources with hard X-ray spectra in rich environments for star forming regions, such as the Galactic center., Comment: 10 pages, 6 figures. Accepted by ApJ

In this work, we compare the $\gamma$-ray spectra available in literature from 11 middle aged supernova remnants (SNRs) interacting with molecular clouds (MCs). It is found that 5 remnants prefer a smoothly broken power law proton spectrum with similar power law index but different break energy. The rest of the SNRs need updated data to test whether a spectral break is preferred in the proton spectrum. Then we compare the $\gamma$-ray spectra from all 11 SNRs with the prediction from widely accepted escaping scenario and direct interaction scenario. We show that current $\gamma$-ray data is inconsistent with the escaping model statistically, as it predicts a diversity of $\gamma$-ray spectra which is not detected in the observation. We also find that ambient CRs can be very important for the $\gamma$-ray emission in the MCs external to W28 and W44, which requires further investigation. In the direct interaction scenario, we focus on re-acceleration of pre-existing ambient CRs. The model can produce the overall profile of $\gamma$-ray data with different acceleration time, but it suggests a transition of seed particles in the evolution of SNR. Whether such transition indeed exists has to be tested by future observation. In the end, we propose that radiative SNR without MC interaction can also produce a significant amount of $\gamma$-ray emission. One good candidate is S147. With accumulated Fermi data and CTA in future we expect to detect more remnants like S147., Comment: Proceedings of the 7th International Fermi Symposium, Garmisch, Oct 15-20 2017

We analyze the impact of thermal conduction on the appearance of a shock-heated gas shell which is produced when a spherically symmetric outburst of a supermassive black hole inflates bubbles of relativistic plasma at the center of a galaxy cluster. The presence of the hot and low-density shell can be used as an ancillary indicator for a high rate of energy release during the outburst, which is required to drive strong shocks into the gas. Here we show that conduction can effectively erase such shell, unless the diffusion of electrons is heavily suppressed. We conclude that a more robust proxy to the energy release rate is the ratio between the shock radius and bubble radius. We also revisited the issue of sound waves dissipation induced by thermal conduction in a scenario, where characteristic wavelength of the sound wave is set by the total energy of the outburst. For a fiducial short outburst model, the dissipation length does not exceed the cooling radius in a typical cluster, provided that the conduction is suppressed by a factor not larger than $\sim$100. For quasi-continuous energy injection neither the shock-heated shell nor the outgoing sound wave are important and the role of conduction is subdominant., Comment: 12 pages, 10 figures, MNRAS in press

We compare the $\gamma$-ray spectra from 10 middle-aged supernova remnants (SNRs), which are interacting with molecular clouds (MCs), with the model prediction from widely used escaping scenario and direct interaction scenario. It is found that the $\gamma$-ray data is inconsistent with the escaping scenario statistically, as it predicts a diversity of spectral shape which is not observed. The inconsistency suggests that the free escape boundary adopted in the escaping model is not a good approximation, which challenges our understanding of cosmic ray (CR) escaping in SNRs. In addition, we show that ambient CRs is potentially important for the $\gamma$-ray emission of illuminated MCs external to W28 and W44. In direct interaction scenario, the model involving re-acceleration of pre-existing CRs and adiabatic compression is able to explain the emission from most SNRs. The dispersion shown in the TeV data is naturally explained by different acceleration time of CR particles in SNRs. Re-acceleration of pre-existing CRs suggests a transition of seed particles, which is from thermal injected seed particle in young SNRs to ambient CRs in old SNRs. The transition needs to be tested by future multi-wavelength observation. In the end, we propose that radiative SNR without MC interaction is also able to produce a significant amount of $\gamma$-ray emission. A good candidate is S147. With accumulated Fermi data and CTA in future we expect to detect more remnants like S147., Comment: 15 pages and 7 figures. Accepted by MNRAS

We consider the evolution of an outburst in a uniform medium under spherical symmetry, having in mind AGN feedback in the intra cluster medium (ICM). For a given density and pressure of the medium, the spatial structure and energy partition at a given time $t_{age}$ (since the onset of the outburst) are fully determined by the total injected energy $E_{inj}$ and the duration $t_b$ of the outburst. We are particularly interested in the late phase evolution when the strong shock transforms into a sound wave. We studied the energy partition during such transition with different combinations of $E_{inj}$ and $t_b$. For an instantaneous outburst with $t_b\rightarrow 0$, which corresponds to the extension of classic Sedov-Taylor solution with counter-pressure, the fraction of energy that can be carried away by sound waves is $\lesssim$12% of $E_{inj}$. As $t_b$ increases, the solution approaches the "slow piston" limit, with the fraction of energy in sound waves approaching zero. We then repeat the simulations using radial density and temperature profiles measured in Perseus and M87/Virgo clusters. We find that the results with a uniform medium broadly reproduce an outburst in more realistic conditions once proper scaling is applied. We also develop techniques to map intrinsic properties of an outburst $(E_{inj}, t_b$ and $t_{age})$ to the observables like the Mach number of the shock and radii of the shock and ejecta. For the Perseus cluster and M87, the estimated $(E_{inj}, t_b$ and $t_{age})$ agree with numerical simulations tailored for these objects with $20-30\%$ accuracy., Comment: Accepted by MNRAS, add one figure in appendix and minor changes in text based on referee's comment

Recent gamma-ray observations show that middle aged supernova remnants (SNRs) interacting with molecular clouds (MCs) can be sources of both GeV and TeV emission. Based on the MC association, two scenarios have been proposed to explain the observed gamma-ray emission. In one, energetic cosmic ray (CR) particles escape from the SNR and then illuminate nearby MCs, producing gamma-ray emission, while the other involves direct interaction between the SNR and MC. In the direct interaction scenario, re-acceleration of pre-existing CRs in the ambient medium is investigated while particles injected from the thermal pool are neglected in view of the slow shock speeds in middle aged SNRs. However, standard diffusive shock acceleration (DSA) theory produces a steady state particle spectrum that is too flat compared to observations, which suggests that the high energy part of the observed spectrum has not yet reached a steady state. We derive a time dependent DSA solution in the test particle limit for re-acceleration of pre-existing CRs case and show that it is capable of reproducing the observed gamma-ray emission in SNRs like IC 443 and W44, in the context of a MC interaction model. We also provide a simple physical picture to understand the time dependent DSA spectrum. A spatially averaged diffusion coefficient around the SNR can be estimated through fitting the gamma-ray spectrum. The spatially averaged diffusion coefficient in middle aged SNRs like IC 443 and W44 is estimated to be ~ 10^25 cm^2/s at ~ 1GeV, which is between the Bohm limit and interstellar value., Comment: Proceedings of the conference "Supernova Remnants: an Odyssey in Space after Stellar Death"

We present a new analytical approach to derive approximate solutions describing the shock evolution in non-radiative supernova remnants (SNRs). We focus on the study of the forward shock and contact discontinuity while application to the reverse shock is only discussed briefly. The spherical shock evolution of a SNR in both the interstellar medium with a constant density profile and a circumstellar medium with a wind density profile is investigated. We compared our new analytical solution with numerical simulations and found that a few percent accuracy is achieved. For the evolution of the forward shock, we also compared our new solution to previous analytical models. In a uniform ambient medium, the accuracy of our analytical approximation is comparable to that in Truelove&McKee (1999). In a wind density profile medium, our solution performs better than that in Micelotta et al. (2016), especially when the ejecta envelope has a steep density profile. The new solution is significantly simplified compared to previous analytical models, as it only depends on the asymptotic behaviors of the remnant during its evolution., Comment: Add discussion for contact discontinuity

Recent gamma ray observations show that middle aged supernova remnants interacting with molecular clouds can be sources of both GeV and TeV emission. Models involving re-acceleration of pre-existing cosmic rays in the ambient medium and direct interaction between supernova remnant and molecular clouds have been proposed to explain the observed gamma ray emission. For the re-acceleration process, standard DSA theory in the test particle limit produces a steady state particle spectrum that is too flat compared to observations, which suggests that the high energy part of the observed spectrum has not yet reached a steady state. We derive a time dependent DSA solution in the test particle limit for situations involving re-acceleration of pre-existing cosmic rays in the preshock medium. Simple estimates with our time dependent DSA solution plus a molecular cloud interaction model can reproduce the overall shape of the spectra of IC 443 and W44 from GeV to TeV energies through pure $\pi^0$-decay emission. We allow for a power law momentum dependence of the diffusion coefficient, finding that a power law index of 0.5 is favored., Comment: 26 pages, 5 figures, accepted by ApJ

Observations of the middle-aged supernova remnants IC 443, W28 and W51C indicate that the brightnesses at GeV and TeV energies are correlated with each other and with regions of molecular clump interaction, but not with the radio synchrotron brightness. We suggest that the radio emission is primarily associated with a radiative shell in the interclump medium of a molecular cloud, while the gamma-ray emission is primarily associated with the interaction of the radiative shell with molecular clumps. The shell interaction produces a high pressure region, so that the gamma-ray luminosity can be approximately reproduced even if shock acceleration of particles is not efficient, provided that energetic particles are trapped in the cooling region. In this model, the spectral shape \ga 2 GeV is determined by the spectrum of cosmic ray protons. Models in which diffusive shock acceleration determines the spectrum tend to underproduce TeV emission because of the limiting particle energy that is attained., Comment: 15 pages, ApJ Letters, accepted

The model for pulsar wind nebulae (PWNe) as the result of the magnetohydrodynamic (MHD) downstream flow from a shocked, relativistic pulsar wind has been successful in reproducing many features of the nebulae observed close to the central pulsars. However, observations of well-studied young nebulae like the Crab Nebula, 3C 58, and G21.5--0.9 do not show the toroidal magnetic field on a larger scale that might be expected in the MHD flow model; in addition, the radial variation of spectral index due to synchrotron losses is smoother than expected in the MHD flow model. We find that pure diffusion models can reproduce the basic data on nebular size and spectral index variation for the Crab, 3C 58, and G21.5--0.9. Most of our models use an energy independent diffusion coefficient; power law variations of the coefficient with energy are degenerate with variation in the input particle energy distribution index in the steady state, transmitting boundary case. Energy dependent diffusion is a possible reason for the smaller diffusion coefficient inferred for the Crab. Monte Carlo simulations of the particle transport allowing for advection and diffusion of particles suggest that diffusion dominates over much of the total nebular volume of the Crab. Advection dominates close to the pulsar and is likely to play a role in the X-ray half-light radius. The source of diffusion and mixing of particles is uncertain, but may be related to the Rayleigh-Taylor instability at the outer boundary of a young PWN or to instabilities in the toroidal magnetic field structure., Comment: 13 pages, ApJ, in press, corrected typo