Showing papers in "arXiv: High Energy Astrophysical Phenomena in 2010"

Predictions for the Rates of Compact Binary Coalescences Observable by Ground-based Gravitational-wave Detectors

Abstract: We present an up-to-date, comprehensive summary of the rates for all types of compact binary coalescence sources detectable by the Initial and Advanced versions of the ground-based gravitational-wave detectors LIGO and Virgo Astrophysical estimates for compact-binary coalescence rates depend on a number of assumptions and unknown model parameters, and are still uncertain The most confident among these estimates are the rate predictions for coalescing binary neutron stars which are based on extrapolations from observed binary pulsars in our Galaxy These yield a likely coalescence rate of 100 per Myr per Milky Way Equivalent Galaxy (MWEG), although the rate could plausibly range from 1 per Myr per MWEG to 1000 per Myr per MWEG We convert coalescence rates into detection rates based on data from the LIGO S5 and Virgo VSR2 science runs and projected sensitivities for our Advanced detectors Using the detector sensitivities derived from these data, we find a likely detection rate of 002 per year for Initial LIGO-Virgo interferometers, with a plausible range between 00002 and 02 per year The likely binary neutron-star detection rate for the Advanced LIGO-Virgo network increases to 40 events per year, with a range between 04 and 400 per year

Giant Gamma-ray Bubbles from Fermi-LAT: AGN Activity or Bipolar Galactic Wind?

Abstract: Data from the Fermi-LAT reveal two large gamma-ray bubbles, extending 50 degrees above and below the Galactic center, with a width of about 40 degrees in longitude. The gamma-ray emission associated with these bubbles has a significantly harder spectrum (dN/dE ~ E^-2) than the IC emission from electrons in the Galactic disk, or the gamma-rays produced by decay of pions from proton-ISM collisions. There is no significant spatial variation in the spectrum or gamma-ray intensity within the bubbles, or between the north and south bubbles. The bubbles are spatially correlated with the hard-spectrum microwave excess known as the WMAP haze; the edges of the bubbles also line up with features in the ROSAT X-ray maps at 1.5-2 keV. We argue that these Galactic gamma-ray bubbles were most likely created by some large episode of energy injection in the Galactic center, such as past accretion events onto the central massive black hole, or a nuclear starburst in the last ~10 Myr. Dark matter annihilation/decay seems unlikely to generate all the features of the bubbles and the associated signals in WMAP and ROSAT; the bubbles must be understood in order to use measurements of the diffuse gamma-ray emission in the inner Galaxy as a probe of dark matter physics. Study of the origin and evolution of the bubbles also has the potential to improve our understanding of recent energetic events in the inner Galaxy and the high-latitude cosmic ray population.

The Proto-Magnetar Model for Gamma-Ray Bursts

Abstract: Long duration Gamma-Ray Bursts (GRBs) originate from the core collapse of massive stars, but the identity of the central engine remains elusive. Previous work has shown that rapidly spinning, strongly magnetized proto-neutron stars (`millisecond proto-magnetars') produce outflows with energies, timescales, and magnetizations sigma_0 (maximum Lorentz factor) that are consistent with those required to produce long GRBs. Here we extend this work in order to construct a self-consistent model that directly connects the properties of the central engine to the observed prompt emission. Just after the launch of the supernova shock, a wind heated by neutrinos is driven from the proto-magnetar. The outflow is collimated into a bipolar jet by its interaction with the star. As the magnetar cools, the wind becomes ultra-relativistic and Poynting-flux dominated (sigma_0 >> 1) on a timescale comparable to that required for the jet to clear a cavity through the star. Although the site and mechanism of the prompt emission are debated, we calculate the emission predicted by two models: magnetic dissipation and internal shocks. Our results favor the magnetic dissipation model in part because it predicts a relatively constant `Band' spectral peak energy E_peak with time during the GRB. The jet baryon loading decreases abruptly when the neutron star becomes transparent to neutrinos at t ~ 10-100 seconds. Jets with ultra-high magnetization cannot effectively accelerate and dissipate their energy, suggesting this transition ends the prompt emission and may explain the steep decay phase that follows. We assess several phenomena potentially related to magnetar birth, including low luminosity GRBs, thermal-rich GRBs/X-ray Flashes, very luminous supernovae, and short duration GRBs with extended emission.

The First Catalog of Active Galactic Nuclei Detected by the Fermi Large Area Telescope

Abstract: We present the first catalog of active galactic nuclei (AGN) detected by the LAT, corresponding to 11 months of data collected in scientific operation mode. The First LAT AGN Catalog (1LAC) includes 671 gamma-ray sources located at high Galactic latitudes (|b| > 10 deg) that are detected with a test statistic (TS) greater than 25 and associated statistically with AGNs. Some LAT sources are associated with multiple AGNs, and consequently, the catalog includes 709 AGNs, comprising 300 BL Lacertae objects (BL Lacs), 296 flat-spectrum radio quasars (FSRQs), 41 AGNs of other types, and 72 AGNs of unknown type. We also classify the blazars based on their spectral energy distributions (SEDs) as archival radio, optical, and X-ray data permit. In addition to the format 1LAC sample, we provide AGN associations for 51 low-latitude LAT sources and AGN "affiliations" (unquantified counterpart candidates) for 104 high-latitude LAT sources without AGN associations. The overlap of the 1LAC with existing gamma-ray AGN catalogs (LBAS, EGRET, AGILE, Swift, INTEGRAL, TeVCat) is briefly discussed. Various properties--such as gamma-ray fluxes and photon power law spectral indices, redshifts, gamma-ray luminosities, variability, and archival radio luminosities--and their correlations are presented and discussed for the different blazar classes. We compare the 1LAC results with predictions regarding the gamma-ray AGN populations, and we comment on the power of the sample to address the question of the blazar sequence.

The Chandra Source Catalog

Abstract: The Chandra Source Catalog (CSC) is a general purpose virtual X-ray astrophysics facility that provides access to a carefully selected set of generally useful quantities for individual X-ray sources, and is designed to satisfy the needs of a broad-based group of scientists, including those who may be less familiar with astronomical data analysis in the X-ray regime. The first release of the CSC includes information about 94,676 distinct X-ray sources detected in a subset of public ACIS imaging observations from roughly the first eight years of the Chandra mission. This release of the catalog includes point and compact sources with observed spatial extents <~ 30''. The catalog (1) provides access to the best estimates of the X-ray source properties for detected sources, with good scientific fidelity, and directly supports scientific analysis using the individual source data; (2) facilitates analysis of a wide range of statistical properties for classes of X-ray sources; and (3) provides efficient access to calibrated observational data and ancillary data products for individual X-ray sources, so that users can perform detailed further analysis using existing tools. The catalog includes real X-ray sources detected with flux estimates that are at least 3 times their estimated 1 sigma uncertainties in at least one energy band, while maintaining the number of spurious sources at a level of <~ 1 false source per field for a 100 ks observation. For each detected source, the CSC provides commonly tabulated quantities, including source position, extent, multi-band fluxes, hardness ratios, and variability statistics, derived from the observations in which the source is detected. In addition to these traditional catalog elements, for each X-ray source the CSC includes an extensive set of file-based data products that can be manipulated interactively.

Magnetic acceleration of relativistic jets

Abstract: This is a brief review of the recent developments in the theory of magnetic acceleration of relativistic jets. We attempt to explain the key results of this complex theory using basic physical arguments and simple calculations. The main focus is on the standard model, which describes steady-state axisymmetric ideal MHD flows. We argue that this model is over-restrictive and discuss various alternatives.

An XMM-Newton spectral survey of 12 micron selected galaxies. I. X-ray data

Abstract: We present an X-ray spectral analysis of 126 galaxies of the 12 micron galaxy sample. We pay particular attention to Compton thick AGN with the help of new spectral fitting models that we have produced, which are based on Monte-Carlo simulations of X-ray radiative transfer, using both a spherical and torus geometry, and taking into account Compton scattering and Fe fluorescence. We use this data to show that with a torus geometry, unobscured sight lines can achieve a maximum EW of the Fe K\alpha line of ~150 eV, originally shown by Ghisellini, Haardt & Matt (1994). In order for this to be exceeded, the line of sight must be obscured with N_H>10^23 cm^-2, as we show for one case, NGC 3690. We also calculate flux suppression factors from the simulated data, the main conclusion from which is that for N_H>10^25 cm^-2, the X-ray flux is suppressed by a factor of >10 in all X-ray bands and at all redshifts, revealing the biases present against these extremely heavily obscured systems inherent in all X-ray surveys. Furthermore, we confirm previous results from Murphy & Yaqoob (2009) that show that the reflection fraction determined from slab geometries is underestimated with respect to toroidal geometries. For the 12 micron selected galaxies, we investigate the distribution of X-ray power-law indices, finding that the mean $ =1.90_{-0.07}^{+0.05}$ and $\sigma_\Gamma = 0.31_{-0.05}^{+0.05}$) is consistent with previous works, and that the distribution of \Gamma for obscured and unobscured sources is consistent with the source populations being the same, in general support of unification schemes. We determine a Compton thick fraction for the X-ray AGN in our sample to be 18+/-5% which is higher than the hard X-ray (>10 keV) selected samples. Finally we find that the obscured fraction for our sample is a strong function of X-ray luminosity, peaking at L_X~10^42-43 ergs s^-1.

Detonations in Sub-Chandrasekhar Mass C+O White Dwarfs

Abstract: Explosions of sub-Chandrasekhar-mass white dwarfs are one alternative to the standard Chandrasekhar-mass model of Type Ia supernovae. They are interesting since binary systems with sub-Chandrasekhar-mass primary white dwarfs should be common and this scenario would suggest a simple physical parameter which determines the explosion brightness, namely the mass of the exploding white dwarf. Here we perform one-dimensional hydrodynamical simulations, associated post-processing nucleosynthesis and multi-wavelength radiation transport calculations for pure detonations of carbon-oxygen white dwarfs. The light curves and spectra we obtain from these simulations are in good agreement with observed properties of Type Ia supernovae. In particular, for white dwarf masses from 0.97 - 1.15 Msun we obtain 56Ni masses between 0.3 and 0.8 Msun, sufficient to capture almost the complete range of Type Ia supernova brightnesses. Our optical light curve rise times, peak colours and decline timescales display trends which are generally consistent with observed characteristics although the range of B-band decline timescales displayed by our current set of models is somewhat too narrow. In agreement with observations, the maximum light spectra of the models show clear features associated with intermediate mass elements and reproduce the sense of the observed correlation between explosion luminosity and the ratio of the Si II lines at 6355 and 5972 Angstroms. We therefore suggest that sub-Chandrasekhar mass explosions are a viable model for Type Ia supernovae for any binary evolution scenario leading to explosions in which the optical display is dominated by the material produced in a detonation of the primary white dwarf.

Gamma-ray Light Curves and Variability of Bright Fermi-Detected Blazars

Abstract: This paper presents light curves and the first systematic characterization of variability of the 106 objects in the Fermi Large Area Telescope (LAT) Bright AGN Sample (LBAS). Weekly light curves obtained during the first 11 months of survey (August 04, 2008 - July 04, 2009), are tested for variability, and their properties are quantified through autocorrelation and structure function analysis. For the brightest sources power density spectra (PDS) and fit of the temporal structure of major flares is performed. More than 50% of the sources are variable, where high states do not exceed 1/4 of the total observation range. Variation amplitudes are larger for FSRQs and low/intermediate synchrotron peaked (LSP/ISP) BL Lac objects. Autocorrelation time scales vary from 4 to a dozen of weeks. Variable sources of the sample have 1/(f^{a}) PDS and show two modes: (1) rather constant baseline with sporadic flaring activity characterized by flatter PDS slopes resembling flickering and red-noise with occasional intermittence, and (2) - measured for a few blazars showing strong activity - complex and structured temporal profiles characterized by longer-term memory and steeper PDS slopes typical of a random-walk underlying mechanism. The average PDS slope of the brightest 22 FSRQs and the 6 brightest BL Lacs is 1.5 and 1.7 respectively. The study of temporal profiles of well resolved flares shows that they generally have symmetric profiles and that their total duration vary between 10 and 100 days.

Three years Swift-BAT Survey of AGN: reconciling theory and observations?

Abstract: It is well accepted that unabsorbed as well as absorbed AGN are needed to explain the nature and the shape of the Cosmic X-ray background, even if the fraction of highly absorbed objects (dubbed Compton-thick sources) substantially still escapes detection. We derive and analyze the absorption distribution using a complete sample of AGN detected by Swift-BAT in the first three years of the survey. The fraction of Compton-thick AGN represents only 4.6% of the total AGN population detected by Swift-BAT. However, we show that once corrected for the bias against the detection of very absorbed sources the real intrinsic fraction of Compton-thick AGN is 20$^{+9}_{-6}$%. We proved for the first time (also in the BAT band) that the anti-correlation of the fraction of absorbed AGN and luminosity it tightly connected to the different behavior of the luminosity functions (XLFs) of absorbed and unabsorbed AGN. This points towards a difference between the two subsamples of objects with absorbed AGN being, on average, intrinsically less luminous than unobscured ones. Moreover the XLFs show that the fraction of obscured AGN might also decrease at very low luminosity. This can be successfully interpreted in the framework of a disk cloud outflow scenario as the disappearance of the obscuring region below a critical luminosity. Our results are discussed in the framework of population synthesis models and the origin of the Cosmic X-ray Background.

Fermi Observations of GRB 090510: A Short Hard Gamma-Ray Burst with an Additional, Hard Power-Law Component from 10 keV to GeV Energies

Abstract: We present detailed observations of the bright short-hard gamma-ray burst GRB 090510 made with the Gamma-ray Burst Monitor (GBM) and Large Area Telescope (LAT) on board the Fermi observatory. GRB 090510 is the first burst detected by the LAT that shows strong evidence for a deviation from a Band spectral fitting function during the prompt emission phase. The time-integrated spectrum is fit by the sum of a Band function with $\Epeak = 3.9\pm 0.3$\,MeV, which is the highest yet measured, and a hard power-law component with photon index $-1.62\pm 0.03$ that dominates the emission below $\approx$\,20\,keV and above $\approx$\,100\,MeV. The onset of the high-energy spectral component appears to be delayed by $\sim$\,0.1\,s with respect to the onset of a component well fit with a single Band function. A faint GBM pulse and a LAT photon are detected 0.5\,s before the main pulse. During the prompt phase, the LAT detected a photon with energy $30.5^{+5.8}_{-2.6}$ GeV, the highest ever measured from a short GRB. Observation of this photon sets a minimum bulk outflow Lorentz factor, $\Gamma\ga$\,1200, using simple $\gamma\gamma$ opacity arguments for this GRB at redshift $z = 0.903$ and a variability time scale on the order of tens of ms for the $\approx$\,100\,keV--few MeV flux. Stricter high confidence estimates imply $\Gamma \ga 1000$ and still require that the outflows powering short GRBs are at least as highly relativistic as those of long duration GRBs. Implications of the temporal behavior and power-law shape of the additional component on synchrotron/synchrotron self-Compton (SSC), external-shock synchrotron, and hadronic models are considered.

Broad emission lines for negatively spinning black holes

Abstract: We present an extended scheme for the calculation of the profiles of emission lines from accretion discs around rotating black holes. The scheme includes discs with angular momenta which are parallel and antiparallel with respect to the black hole's angular momentum, as both configurations are assumed to be stable (King et al., 2005). We discuss line shapes for such discs and present a code for modelling observational data with this scheme in X-ray data analysis programs. Based on a Green's function approach, an arbitrary radius dependence of the disc emissivity and arbitrary limb darkening laws can be easily taken into account, while the amount of precomputed data is significantly reduced with respect to other available models.

Teraelectronvolt Astronomy

Abstract: Ground-based gamma-ray astronomy, which provides access to the TeV energy range, is a young and rapidly developing discipline Recent discoveries in this waveband have important consequences for a wide range of topics in astrophysics and astroparticle physics This article is an attempt to review the experimental status of this field and to provide the basic formulae and concepts required to begin the interpretation of TeV observations

Discrepant hardening observed in cosmic-ray elemental spectra

Abstract: The balloon-borne Cosmic Ray Energetics And Mass (CREAM) experiment launched five times from Antarctica has achieved a cumulative flight duration of about 156 days above 99.5% of the atmosphere. The instrument is configured with complementary and redundant particle detectors designed to extend direct measurements of cosmic-ray composition to the highest energies practical with balloon flights. All elements from protons to iron nuclei are separated with excellent charge resolution. Here we report results from the first two flights of ~70 days, which indicate hardening of the elemental spectra above ~200 GeV/nucleon and a spectral difference between the two most abundant species, protons and helium nuclei. These results challenge the view that cosmic-ray spectra are simple power laws below the so-called knee at ~1015 eV. This discrepant hardening may result from a relatively nearby source, or it could represent spectral concavity caused by interactions of cosmic rays with the accelerating shock. Other possible explanations should also be investigated.

Sub-Chandrasekhar Mass Models For Type Ia Supernovae

Abstract: For carbon-oxygen white dwarfs accreting hydrogen or helium at rates in the range ~1-10 x 10^(-8) Msun/y, a variety of explosive outcomes is possible well before the star reaches the Chandrasekhar mass. These outcomes are surveyed for a range of white dwarf masses (0.7 - 1.1 Msun), accretion rates (1 - 7 x 10^(-8) Msun/y), and initial white dwarf temperatures (0.01 and 1 Lsun). The results are particularly sensitive to the convection that goes on during the last few minutes before the explosion. Unless this convection maintains a shallow temperature gradient, and unless the density is sufficiently high, the accreted helium does not detonate. Below a critical helium ignition density, which we estimate to be 5 - 10 x 10^5 g cm^(-3), either helium novae or helium deflagrations result. The hydrodynamics, nucleosynthesis, light curves, and spectra of a representative sample of detonating and deflagrating models are explored. Some can be quite faint indeed, powered at peak for a few days by the decay of 48Cr and 48V. Only the hottest, most massive white dwarfs considered with the smallest helium layers, show reasonable agreement with the light curves and spectra of common Type Ia supernovae. For the other models, especially those involving lighter white dwarfs, the helium shell mass exceeds 0.05 Msun and the mass of the 56Ni that is synthesized exceeds 0.01 Msun. These explosions do not look like ordinary Type Ia supernovae, or any other frequently observed transient.

Simulations of Magnetized Disks Around Black Holes: Effects of Black Hole Spin, Disk Thickness, and Magnetic Field Geometry

Abstract: The standard general relativistic model of a razor-thin accretion disk around a black hole, developed by Novikov & Thorne (NT) in 1973, assumes the shear stress vanishes at the radius of the innermost stable circular orbit (ISCO) and that, outside the ISCO, the shear stress is produced by an effective turbulent viscosity. However, astrophysical accretion disks are not razor-thin, it is uncertain whether the shear stress necessarily vanishes at the ISCO, and the magnetic field, which is thought to drive turbulence in disks, may contain large-scale structures that do not behave like a simple local scalar viscosity. We describe three-dimensional general relativistic magnetohydrodynamic simulations of accretion disks around black holes with a range of spin parameters, and we use the simulations to assess the validity of the NT model. Our fiducial initial magnetic field consists of multiple (alternating polarity) poloidal field loops whose shape is roughly isotropic in the disk in order to match the isotropic turbulence expected in the poloidal plane. For a thin disk with an aspect ratio |h/r| ~ 0.07 around a non-spinning black hole, we find a decrease in the accreted specific angular momentum of 2.9% relative to the NT model and an excess luminosity from inside the ISCO of 3.5%. The deviations in the case of spinning black holes are of the same order. In addition, the deviations decrease with decreasing |h/r|. We therefore conclude that magnetized thin accretion disks in x-ray binaries in the thermal/high-soft spectral state ought to be well-described by the NT model, especially at luminosities below 30% of Eddington where we expect a very small disk thickness |h/r| <~ 0.05. We also discuss how the stress and the luminosity inside the ISCO depend on the assumed initial magnetic field geometry. (abridged)

Fermi LAT Observations of the Supernova Remnant W28 (G6.4-0.1)

Abstract: We present detailed analysis of the two gamma-ray sources,1FGL J18013-2322c and 1FGL J18005-2359c,that have been found toward the supernova remnant(SNR) W28 with the Large Area Telescope(LAT) on board the Fermi Gamma-ray Space Telescope1FGL J18013-2322c is found to be an extended source within the boundary of SNR W28,and to extensively overlap with the TeV gamma-ray source HESS J1801-233,which is associated with a dense molecular cloud interacting with the supernova remnantThe gamma-ray spectrum measured with LAT from 02--100 GeV can be described by a broken power-law function with a break of ~1GeV,and photon indices of 209$\pm$008(stat)$\pm$028(sys) below the break and 274$\pm$006(stat)$\pm$009(sys) above the breakGiven the clear association between HESS J1801-233 and the shocked molecular cloud and a smoothly connected spectrum in the GeV--TeV band,we consider the origin of the gamma-ray emission in both GeV and TeV ranges to be the interaction between particles accelerated in the SNR and the molecular cloudThe decay of neutral pions produced in interactions between accelerated hadrons and dense molecular gas provide a reasonable explanation for the broadband gamma-ray spectrum 1FGL J18005-2359c, located outside the southern boundary of SNR W28, cannot be resolvedAn upper limit on the size of the gamma-ray emission was estimated to be ~16$'$ using events above ~2GeV under the assumption of a circular shape with uniform surface brightness It appears to coincide with the TeV source HESS J1800-240B,which is considered to be associated with a dense molecular cloud that contains the ultra compact HII region W28A2(G589-039)We found no significant gamma-ray emission in the LAT energy band at the positions of TeV sources HESS J1800-230A and HESS J1800-230CThe LAT data for HESS J1800-230A combined with the TeV data points indicate a spectral break between 10GeV and 100GeV

Particle Acceleration in Relativistic Magnetized Collisionless Electron-Ion Shocks

Abstract: We investigate shock structure and particle acceleration in relativistic magnetized collisionless electron-ion shocks by means of 2.5D particle-in-cell simulations with ion-to-electron mass ratios (m_i/m_e) ranging from 16 to 1000. We explore a range of inclination angles between the pre-shock magnetic field and the shock normal. In "subluminal" shocks, where relativistic particles can escape ahead of the shock along the magnetic field lines, ions are efficiently accelerated via a Fermi-like mechanism. The downstream ion spectrum consists of a relativistic Maxwellian and a high-energy power-law tail, which contains ~5% of ions and ~30% of ion energy. Its slope is -2.1. Upstream electrons enter the shock with lower energy than ions, so they are more strongly tied to the field. As a result, only ~1% of the incoming electrons are Fermi-accelerated at the shock before being advected downstream, where they populate a steep power-law tail (with slope -3.5). For "superluminal" shocks, where relativistic particles cannot outrun the shock along the field, the self-generated turbulence is not strong enough to permit efficient Fermi acceleration, and the ion and electron downstream spectra are consistent with thermal distributions. The incoming electrons are heated up to equipartition with ions, due to strong electromagnetic waves emitted by the shock into the upstream. Thus, efficient electron heating (>15% of the upstream ion energy) is the universal property of relativistic electron-ion shocks, but significant nonthermal acceleration of electrons (>2% by number, >10% by energy, with slope flatter than -2.5) is hard to achieve in magnetized flows and requires weakly magnetized shocks (magnetization <1e-3). These findings place important constraints on the models of AGN jets and Gamma Ray Bursts that invoke particle acceleration in relativistic magnetized electron-ion shocks.

Precise Gamma-Ray Timing and Radio Observations of 17 Fermi Gamma-Ray Pulsars

Abstract: We present precise phase-connected pulse timing solutions for 16 gamma-ray-selected pulsars recently discovered using the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope plus one very faint radio pulsar (PSR J1124-5916) that is more effectively timed with the LAT. We describe the analysis techniques including a maximum likelihood method for determining pulse times of arrival from unbinned photon data. A major result of this work is improved position determinations, which are crucial for multi-wavelength follow up. For most of the pulsars, we overlay the timing localizations on X-ray images from Swift and describe the status of X-ray counterpart associations. We report glitches measured in PSRs J0007+7303, J1124-5916, and J1813-1246. We analyze a new 20 ks Chandra ACIS observation of PSR J0633+0632 that reveals an arcminute-scale X-ray nebula extending to the south of the pulsar. We were also able to precisely localize the X-ray point source counterpart to the pulsar and find a spectrum that can be described by an absorbed blackbody or neutron star atmosphere with a hard powerlaw component. Another Chandra ACIS image of PSR J1732-3131 reveals a faint X-ray point source at a location consistent with the timing position of the pulsar. Finally, we present a compilation of new and archival searches for radio pulsations from each of the gamma-ray-selected pulsars as well as a new Parkes radio observation of PSR J1124-5916 to establish the gamma-ray to radio phase offset.

Observation of Supernova Remnant IC443 with the Fermi Large Area Telescope

Abstract: We report observation of the supernova remnant IC443 (G189.1+3.0) with the Fermi Gamma-ray Space Telescope Large Area Telescope (LAT) in the energy band between 200MeV and 50GeV. IC443 is a shell-type supernova remnant with mixed morphology located off the outer Galactic plane where high-energy emission has been detected in the X-ray, GeV and TeV gamma-ray bands. Past observations suggest IC443 has been interacting with surrounding interstellar matter. Proximity between dense shocked molecular clouds and GeV-TeV gamma-ray emission regions detected by EGRET, MAGIC and VERITAS suggests an interpretation that cosmic-ray (CR) particles are accelerated by the SNR. With the high gamma-ray statistics and broad energy coverage provided by the LAT, we accurately characterize the gamma-ray emission produced by the CRs accelerated at IC443. The emission region is extended in the energy band with theta_68 = 0.27 deg +/- 0.01 deg (stat) +/- 0.03 deg (sys) for an assumed 2-dimensional Gaussian profile and overlaps almost completely with the extended source region of VERITAS. Its centroid is displaced significantly from the known pulsar wind nebula (PWN) which suggests the PWN is not the major contributor in the present energy band. The observed spectrum changes its power-law slope continuously and continues smoothly to the MAGIC and VERITAS data points. The combined gamma-ray spectrum (200MeV

A Comprehensive Analysis of Fermi Gamma-Ray Burst Data. I. Spectral Components and Their Possible Physical Origins of LAT/GBM GRBs

Abstract: We present a systematic analysis of the spectral and temporal properties of 17 GRBs co-detected by GBM and LAT on board the Fermi satellite by May 2010. We performed a time-resolved spectral analysis of all the bursts with the finest temporal resolution allowed by statistics, in order to avoid temporal smearing of different spectral components. We found that the time-resolved spectra of 14 out of 17 GRBs are best modeled with the Band function over the entire Fermi spectral range, which may suggest a common origin for emissions detected by LAT and GBM. GRB 090902B and GRB 090510 require the superposition between an MeV component and an extra power law component, with the former having a sharp cutoff above E_p. For GRB 090902B, this MeV component becomes progressively narrower as the time bin gets smaller, and can be fit with a Planck function as the time bin becomes small enough. In general, we speculate that phenomenologically there may be three elemental spectral components : (I) a Band-function component (e.g. in GRB 080916C) that extends in a wide energy range and does not narrow with reducing time bins, which may be of the non-thermal origin; (II) a quasi-thermal component (e.g. in GRB 090902B) with the spectra progressively narrowing with reducing time bins; and (III) another non-thermal power law component extending to high energies. The spectra of different bursts may be decomposed into one or more of these elemental components. We compare this sample with the BATSE sample and investigate some correlations among spectral parameters. We discuss the physical implications of the data analysis results for GRB prompt emission, including jet compositions (matter-dominated vs. Poynting-flux-dominated outflow), emission sites (internal shock, external shock or photosphere), as well as radiation mechanisms (synchrotron, synchrotron self-Compton, or thermal Compton upscattering).

The nature of "dark" gamma-ray bursts

Abstract: Context: Thirteen years after the discovery of the first afterglows, the nature of dark gamma-ray bursts (GRB) still eludes explanation: while each ng-duration GRB typically has an X-ray afterglow, optical/NIR emission is only seen for 40-60% of them. Aim: Here we use the afterglow detection statistics of the systematic follow-up observations performed with GROND since mid-2007 in order to derive the fraction of "dark bursts" according to different methods, and to distinguish between various scenarios for "dark bursts". Method: Observations were performed with the 7-channel "Gamma-Ray Optical and Near-infrared Detector" (GROND) at the 2.2m MPI/ESO telescope. We used the afterglow detection rate in dependence on the delay time between GRB and the first GROND exposure. Results: For long-duration Swift bursts with a detected X-ray afterglow, we achieve a 90% (35/39) detection rate of optical/NIR afterglows whenever our observations started within less than 240\,min after the burst. Complementing our GROND data with Swift/XRT spectra we construct broad-band spectral energy distributions and derive rest-frame extinctions. e detect 25-40% "dark bursts", depending on the definition used. The faint optical afterglow emission of "dark bursts" is mainly due to a combination of two contributing factors: (i) moderate intrinsic extinction at moderate redshifts, and (ii) about 22% of "dark" bursts at redshift $>$5.

Double-detonation sub-Chandrasekhar supernovae: synthetic observables for minimum helium shell mass models

Abstract: Abridged. In the double detonation scenario for Type Ia supernovae (SNe Ia) a detonation initiates in a shell of He-rich material accreted from a companion star by a sub-Chandrasekhar-mass White Dwarf (WD). This shell detonation drives a shock front into the carbon-oxygen (C/O) WD that triggers a secondary detonation in the core. The core detonation results in a complete disruption of the WD. Earlier studies concluded that this scenario has difficulties in accounting for the observed properties of SNe Ia since the explosion ejecta are surrounded by the products of explosive He burning in the shell. Recently, it was proposed that detonations might be possible for much less massive He shells than previously assumed. Moreover, it was shown that even detonations of these minimum He shell masses robustly trigger detonations of the C/O core. Here we present time-dependent multi-wavelength radiative transfer calculations for models with minimum He shell mass and derive synthetic observables for both the optical and {\gamma}-ray spectral regions. These differ strongly from those found in earlier simulations of sub-Chandrasekhar-mass explosions in which more massive He shells were considered. Our models predict light curves which cover both the range of brightnesses and the rise and decline times of observed SNe Ia. However, their colours and spectra do not match the observations. In particular, their B-V colours are generally too red. We show that this discrepancy is mainly due to the composition of the burning products of the He shell of our models which contain significant amounts of Ti and Cr. Using a toy model, we also show that the burning products of the He shell depend crucially on its initial composition. This leads us to conclude that good agreement between sub-Chandrasekhar-mass explosions and observed SNe Ia may still be feasible but further study of the shell properties is required.

Constraints on the Galactic Halo Dark Matter from Fermi-LAT Diffuse Measurements

Abstract: We have performed an analysis of the diffuse gamma-ray emission with the Fermi Large Area Telescope in the Milky Way Halo region searching for a signal from dark matter annihilation or decay. In the absence of a robust dark matter signal, constraints are presented. We consider both gamma rays produced directly in the dark matter annihilation/decay and produced by inverse Compton scattering of the e+e- produced in the annihilation/decay. Conservative limits are derived requiring that the dark matter signal does not exceed the observed diffuse gamma-ray emission. A second set of more stringent limits is derived based on modeling the foreground astrophysical diffuse emission using the GALPROP code. Uncertainties in the height of the diffusive cosmic-ray halo, the distribution of the cosmic-ray sources in the Galaxy, the index of the injection cosmic-ray electron spectrum and the column density of the interstellar gas are taken into account using a profile likelihood formalism, while the parameters governing the cosmic-ray propagation have been derived from fits to local cosmic-ray data. The resulting limits impact the range of particle masses over which dark matter thermal production in the early Universe is possible, and challenge the interpretation of the PAMELA/Fermi-LAT cosmic ray anomalies as annihilation of dark matter.

Global cosmic-ray related luminosity and energy budget of the Milky Way

Abstract: We use the GALPROP code for cosmic-ray (CR) propagation to calculate the broad-band luminosity spectrum of the Milky Way related to CR propagation and interactions in the interstellar medium. This includes gamma-ray emission from the production and subsequent decay of neutral pions, bremsstrahlung, and inverse Compton scattering, and synchrotron radiation. The Galaxy is found to be nearly a CR electron calorimeter, but {\it only} if gamma ray emitting processes are taken into account. Synchrotron radiation alone accounts for only one third of the total electron energy losses with ~10-20% of the total synchrotron emission from secondary CR electrons and positrons. The relationship between far-infrared and radio luminosity that we find from our models is consistent with that found for galaxies in general. The results will be useful for understanding the connection between diffuse emissions from radio through gamma rays in ``normal'' (non-AGN dominated) galaxies, as well as for estimating the broad-band extragalactic diffuse background from these kinds of galaxies.

Insights Into the High-Energy Gamma-ray Emission of Markarian 501 from Extensive Multifrequency Observations in the Fermi Era

Abstract: We report on the gamma-ray activity of the blazar Mrk 501 during the first 480 days of Fermi operation. We find that the average LAT gamma-ray spectrum of Mrk 501 can be well described by a single power-law function with a photon index of 1.78 +/- 0.03. While we observe relatively mild flux variations with the Fermi-LAT (within less than a factor of 2), we detect remarkable spectral variability where the hardest observed spectral index within the LAT energy range is 1.52 +/- 0.14, and the softest one is 2.51 +/- 0.20. These unexpected spectral changes do not correlate with the measured flux variations above 0.3GeV. In this paper, we also present the first results from the 4.5-month-long multifrequency campaign (2009 March 15 - August 1) on Mrk 501, which included the VLBA, Swift, RXTE, MAGIC and VERITAS, the F-GAMMA, GASP-WEBT, and other collaborations and instruments which provided excellent temporal and energy coverage of the source throughout the entire campaign. The average spectral energy distribution of Mrk 501 is well described by the standard one-zone synchrotron self-Compton model. In the framework of this model, we find that the dominant emission region is characterized by a size <~ 0.1 pc (comparable within a factor of few to the size of the partially-resolved VLBA core at 15-43 GHz), and that the total jet power (~10^{44} erg s^{-1}) constitutes only a small fraction (~10^{-3}) of the Eddington luminosity. The energy distribution of the freshly-accelerated radiating electrons required to fit the time-averaged data has a broken power-law form in the energy range 0.3GeV-10TeV, with spectral indices 2.2 and 2.7 below and above the break energy of 20GeV. We argue that such a form is consistent with a scenario in which the bulk of the energy dissipation within the dominant emission zone of Mrk 501 is due to relativistic, proton-mediated shocks.

Modeling the spectral evolution of PWNe inside SNRs

Abstract: We present a new model for the spectral evolution of Pulsar Wind Nebulae inside Supernova Remnants. The model couples the long-term dynamics of these systems, as derived in the 1-D approximation, with a 1-zone description of the spectral evolution of the emitting plasma. Our goal is to provide a simplified theoretical description that can be used as a tool to put constraints on unknown properties of PWN-SNR systems: a piece of work that is preliminary to any more accurate and sophisticated modeling. In the present paper we apply the newly developed model to a few objects of different ages and luminosities. We find that an injection spectrum in the form of a broken-power law gives a satisfactory description of the emission for all the systems we consider. More surprisingly, we also find that the intrinsic spectral break turns out to be at a similar energy for all sources, in spite of the differences mentioned above. We discuss the implications of our findings on the workings of pulsar magnetospheres, pair multiplicity and on the particle acceleration mechanism(s) that might be at work at the pulsar wind termination shock.

Efficient cosmic ray acceleration, hydrodynamics, and Self-consistent Thermal X-ray Emission applied to SNR RX J1713.7-3946

Abstract: We model the broad-band emission from SNR RX J1713.7-3946 including, for the first time, a consistent calculation of thermal X-ray emission together with non-thermal emission in a nonlinear diffusive shock acceleration (DSA) model. Our model tracks the evolution of the SNR including the plasma ionization state between the forward shock and the contact discontinuity. We use a plasma emissivity code to predict the thermal X-ray emission spectrum assuming the initially cold electrons are heated either by Coulomb collisions with the shock heated protons (the slowest possible heating), or come into instant equilibration with the protons. For either electron heating model, electrons reach >10^7 K rapidly and the X-ray line emission near 1 keV is more than 10 times as luminous as the underlying thermal bremsstrahlung continuum. Since recent Suzaku observations show no detectable line emission, this places strong constraints on the unshocked ambient medium density and on the relativistic electron to proton ratio. For the uniform circumstellar medium (CSM) models we consider, the low densities and high relativistic electron to proton ratios required to match the Suzaku X-ray observations definitively rule out pion-decay as the emission process producing GeV-TeV photons. We show that leptonic models, where inverse-Compton scattering against the cosmic background radiation dominates the GeV-TeV emission, produce better fits to the broad-band thermal and non-thermal observations in a uniform CSM.

The Palermo Swift-BAT hard X-ray catalogue III. Results after 54 months of sky survey

Abstract: We present the Second Palermo Swift-BAT hard X-ray catalogue obtained by analysing data acquired in the first 54 months of the Swift mission. Using our software dedicated to the analysis of data from coded mask telescopes, we analysed the BAT survey data in three energy bands (15-30 keV, 15-70 keV, 15-150 keV), obtaining a list of 1256 detections above a significance threshold of 4.8 standard deviations. The identification of the source counterparts is pursued using two strategies: the analysis of field observations of soft X-ray instruments and cross-correlation of our catalogue with source databases.The survey covers 50% of the sky to a 15--150 keV flux limit of 1.0 x 10^-11 erg s^-1 cm^-2 and 9.2 x 10^-12 erg s^-1 cm^-2 for |b| 10 degrees, respectively. The Second Palermo Swift-BAT hard X-ray catalogue includes 1079 (86%) hard X-ray sources with an associated counterpart (26 with a double association and 2 with a triple association) and 177 BAT excesses (14%) that still lack a counterpart. The distribution of the BAT sources among the different object classes consists of 19% Galactic sources, 57% extragalactic sources, and 10% sources with a counterpart at softer energies whose nature has not yet been determined. About half of the BAT associated sources lack a counterpart in the ROSAT catalogues. This suggests that either moderate or strong absorption may be preventing their detection in the ROSAT energy band. The comparison of our BAT catalogue with the Fermi Large Area Telescope First Source Catalogue identifies 59 BAT/Fermi correspondences: 48 blazars, 3 Seyfert galaxies, 1 interacting galaxy, 3 high mass X-ray binaries, and 4 pulsars/supernova remnants. This small number of correspondences indicates that different populations make the sky shine in these two different energy bands.

Cosmogenic Neutrinos: parameter space and detectabilty from PeV to ZeV

Abstract: While propagating from their source to the observer, ultrahigh energy cosmic rays interact with cosmological photon backgrounds and generate to the so-called "cosmogenic neutrinos". Here we study the parameter space of the cosmogenic neutrino flux given recent cosmic ray data and updates on plausible source evolution models. The shape and normalization of the cosmogenic neutrino flux are very sensitive to some of the current unknowns of ultrahigh energy cosmic ray sources and composition. We investigate various chemical compositions and maximum proton acceleration energies E_p,max which are allowed by current observations. We consider different models of source evolution in redshift and three possible scenarios for the Galactic to extragalactic transition. We summarize the parameter space for cosmogenic neutrinos into three regions: an optimistic scenario that is currently being constrained by observations, a plausible range of models in which we base many of our rate estimates, and a pessimistic scenario that will postpone detection for decades to come. We present the implications of these three scenarios for the detection of cosmogenic neutrinos from PeV to ZeV (10^14-21 eV) with the existing and upcoming instruments. In the plausible range of parameters, the narrow flux variability in the EeV energy region assures low but detectable rates for IceCube (0.06-0.2 neutrino per year) and the Pierre Auger Observatory (0.03-0.06 neutrino per year), and detection should happen in the next decade. If EeV neutrinos are detected, PeV information can help select between competing models of cosmic ray composition at the highest energy and the Galactic to extragalactic transition at ankle energies. With improved sensitivity, ZeV neutrino observatories, such as ANITA and JEM-EUSO could explore and place limits on the maximum acceleration energy.