Source author record

J. Chiang

J. Chiang appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

53works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

53 published item(s)

preprint2021arXiv

SNIa-Cosmology Analysis Results from Simulated LSST Images: from Difference Imaging to Constraints on Dark Energy

The Vera Rubin Observatory Legacy Survey of Space and Time (LSST) is expected to process ${\sim}10^6$ transient detections per night. For precision measurements of cosmological parameters and rates, it is critical to understand the detection efficiency, magnitude limits, artifact contamination levels, and biases in the selection and photometry. Here we rigorously test the LSST Difference Image Analysis (DIA) pipeline using simulated images from the Rubin Observatory LSST Dark Energy Science Collaboration (DESC) Data Challenge (DC2) simulation for the Wide-Fast-Deep (WFD) survey area. DC2 is the first large-scale (300 deg$^2$) image simulation of a transient survey that includes realistic cadence, variable observing conditions, and CCD image artifacts. We analyze ${\sim}$15 deg$^2$ of DC2 over a 5-year time-span in which artificial point-sources from Type Ia Supernovae (SNIa) light curves have been overlaid onto the images. We measure the detection efficiency as a function of Signal-to-Noise Ratio (SNR) and find a $50\%$ efficiency at $\rm{SNR}=5.8$. The magnitude limits for each filter are: $u=23.66$, $g=24.69$, $r=24.06$, $i=23.45$, $z=22.54$, $y=21.62$ $\rm{mag}$. The artifact contamination is $\sim90\%$ of detections, corresponding to $\sim1000$ artifacts/deg$^2$ in $g$ band, and falling to 300 per deg$^2$ in $y$ band. The photometry has biases $<1\%$ for magnitudes $19.5 < m <23$. Our DIA performance on simulated images is similar to that of the Dark Energy Survey pipeline applied to real images. We also characterize DC2 image properties to produce catalog-level simulations needed for distance bias corrections. We find good agreement between DC2 data and simulations for distributions of SNR, redshift, and fitted light-curve properties. Applying a realistic SNIa-cosmology analysis for redshifts $z<1$, we recover the input cosmology parameters to within statistical uncertainties.

preprint2016arXiv

Development of the Model of Galactic Interstellar Emission for Standard Point-Source Analysis of Fermi Large Area Telescope Data

Most of the celestial gamma rays detected by the Large Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope originate from the interstellar medium when energetic cosmic rays interact with interstellar nucleons and photons. Conventional point and extended source studies rely on the modeling of this diffuse emission for accurate characterization. We describe here the development of the Galactic Interstellar Emission Model (GIEM) that is the standard adopted by the LAT Collaboration and is publicly available. The model is based on a linear combination of maps for interstellar gas column density in Galactocentric annuli and for the inverse Compton emission produced in the Galaxy. We also include in the GIEM large-scale structures like Loop I and the Fermi bubbles. The measured gas emissivity spectra confirm that the cosmic-ray proton density decreases with Galactocentric distance beyond 5 kpc from the Galactic Center. The measurements also suggest a softening of the proton spectrum with Galactocentric distance. We observe that the Fermi bubbles have boundaries with a shape similar to a catenary at latitudes below 20 degrees and we observe an enhanced emission toward their base extending in the North and South Galactic direction and located within 4 degrees of the Galactic Center.

preprint2016arXiv

First NuSTAR observations of the BL Lac - type blazar PKS~2155-304: constraints on the jet content and distribution of radiating particles

We report the first hard X-ray observations with NuSTAR of the BL Lac type blazar PKS 2155-304, augmented with soft X-ray data from XMM-Newton and gamma-ray data from the Fermi Large Area Telescope, obtained in April 2013 when the source was in a very low flux state. A joint NuSTAR and XMM spectrum, covering the energy range 0.5 - 60 keV, is best described by a model consisting of a log-parabola component with curvature beta = 0.3(+0.2,-0.1) and a (local) photon index 3.04 +/- 0.15 at photon energy of 2 keV, and a hard power-law tail with photon index 2.2 +/- 0.4. The hard X-ray tail can be smoothly joined to the quasi-simultaneous gamma-ray spectrum by a synchrotron self-Compton component produced by an electron distribution with index p = 2.2. Assuming that the power-law electron distribution extends down to the minimum electron Lorentz factor gamma_min = 1 and that there is one proton per electron, an unrealistically high total jet power L_p of roughly 10^47 erg/s is inferred. This can be reduced by two orders of magnitude either by considering a significant presence of electron-positron pairs with lepton-to-proton ratio of at least 30, or by introducing an additional, low-energy break in the electron energy distribution at the electron Lorentz factor gamma_br1 of roughly 100. In either case, the jet composition is expected to be strongly matter-dominated.

preprint2016arXiv

Localization and broadband follow-up of the gravitational-wave transient GW150914

A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams.

preprint2016arXiv

Measurement of the high-energy gamma-ray emission from the Moon with the Fermi Large Area Telescope

We have measured the gamma-ray emission spectrum of the Moon using the data collected by the Large Area Telescope onboard the Fermi satellite during its first 7 years of operation, in the energy range from 30 MeV up to a few GeV. We have also studied the time evolution of the flux, finding a correlation with the solar activity. We have developed a full Monte Carlo simulation describing the interactions of cosmic rays with the lunar surface. The results of the present analysis can be explained in the framework of this model, where the production of gamma rays is due to the interactions of cosmic-ray proton and helium nuclei with the surface of the Moon. Finally, we have used our simulation to derive the cosmic-ray proton and helium spectra near Earth from the Moon gamma-ray data.

preprint2016arXiv

Supplement: Localization and broadband follow-up of the gravitational-wave transient GW150914

This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the different bands.

preprint2015arXiv

First NuSTAR Observations of Mrk 501 within a Radio to TeV Multi-Instrument Campaign

We report on simultaneous broadband observations of the TeV-emitting blazar Markarian 501 between 1 April and 10 August 2013, including the first detailed characterization of the synchrotron peak with Swift and NuSTAR. During the campaign, the nearby BL Lac object was observed in both a quiescent and an elevated state. The broadband campaign includes observations with NuSTAR, MAGIC, VERITAS, the Fermi Large Area Telescope (LAT), Swift X-ray Telescope and UV Optical Telescope, various ground-based optical instruments, including the GASP-WEBT program, as well as radio observations by OVRO, Metsähovi and the F-Gamma consortium. Some of the MAGIC observations were affected by a sand layer from the Saharan desert, and had to be corrected using event-by-event corrections derived with a LIDAR (LIght Detection And Ranging) facility. This is the first time that LIDAR information is used to produce a physics result with Cherenkov Telescope data taken during adverse atmospheric conditions, and hence sets a precedent for the current and future ground-based gamma-ray instruments. The NuSTAR instrument provides unprecedented sensitivity in hard X-rays, showing the source to display a spectral energy distribution between 3 and 79 keV consistent with a log-parabolic spectrum and hard X-ray variability on hour timescales. None (of the four extended NuSTAR observations) shows evidence of the onset of inverse-Compton emission at hard X-ray energies. We apply a single-zone equilibrium synchrotron self-Compton model to five simultaneous broadband spectral energy distributions. We find that the synchrotron self-Compton model can reproduce the observed broadband states through a decrease in the magnetic field strength coinciding with an increase in the luminosity and hardness of the relativistic leptons responsible for the high-energy emission.

preprint2015arXiv

Gamma-ray flaring activity from the gravitationally lensed blazar PKS 1830-211 observed by Fermi LAT

The Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope routinely detects the highly dust-absorbed, reddened, and MeV-peaked flat spectrum radio quasar PKS 1830-211 (z=2.507). Its apparent isotropic gamma-ray luminosity (E>100 MeV) averaged over $\sim$ 3 years of observations and peaking on 2010 October 14/15 at 2.9 X 10^{50} erg s^{-1}, makes it among the brightest high-redshift Fermi blazars. No published model with a single lens can account for all of the observed characteristics of this complex system. Based on radio observations, one expects time delayed variability to follow about 25 days after a primary flare, with flux about a factor 1.5 less. Two large gamma-ray flares of PKS 1830-211 have been detected by the LAT in the considered period and no substantial evidence for such a delayed activity was found. This allows us to place a lower limit of about 6 on the gamma rays flux ratio between the two lensed images. Swift XRT observations from a dedicated Target of Opportunity program indicate a hard spectrum and with no significant correlation of X-ray flux with the gamma-ray variability. The spectral energy distribution can be modeled with inverse Compton scattering of thermal photons from the dusty torus. The implications of the LAT data in terms of variability, the lack of evident delayed flare events, and different radio and gamma-ray flux ratios are discussed. Microlensing effects, absorption, size and location of the emitting regions, the complex mass distribution of the system, an energy-dependent inner structure of the source, and flux suppression by the lens galaxy for one image path may be considered as hypotheses for understanding our results.

preprint2015arXiv

Multiwavelength Evidence for Quasi-periodic Modulation in the Gamma-ray Blazar PG 1553+113

We report for the first time a gamma-ray and multi-wavelength nearly-periodic oscillation in an active galactic nucleus. Using the Fermi Large Area Telescope (LAT) we have discovered an apparent quasi-periodicity in the gamma-ray flux (E >100 MeV) from the GeV/TeV BL Lac object PG 1553+113. The marginal significance of the 2.18 +/-0.08 year-period gamma-ray cycle is strengthened by correlated oscillations observed in radio and optical fluxes, through data collected in the OVRO, Tuorla, KAIT, and CSS monitoring programs and Swift UVOT. The optical cycle appearing in ~10 years of data has a similar period, while the 15 GHz oscillation is less regular than seen in the other bands. Further long-term multi-wavelength monitoring of this blazar may discriminate among the possible explanations for this quasi-periodicity.

preprint2015arXiv

NuSTAR and multifrequency study of the two high-redshift blazars S5 0836+710 and PKS 2149-306

The most powerful blazars are the flat spectrum radio quasars whose emission is dominated by a Compton component peaking between a few hundred keV and a few hundred MeV. We selected two bright blazars, PKS 2149-306 at redshift z=2.345 and S5 0836+710 at z=2.172, in order to observe them in the hard X-ray band with the NuSTAR satellite. In this band the Compton component is rapidly rising almost up to the peak of the emission. Simultaneous soft-X-rays and UV-optical observations were performed with the Swift satellite, while near-infrared (NIR) data were obtained with the REM telescope. To study their variability, we repeated these observations for both sources on a timescale of a few months. While no fast variability was detected during a single observation, both sources were found to be variable in the X-ray band, up to 50%, between the two observations, with larger variability at higher energies. No variability was detected in the optical/NIR band. These data together with Fermi-LAT, WISE and other literature data are then used to study the overall spectral energy distributions (SEDs) of these blazars. Although the jet non-thermal emission dominates the SED, it leaves the UV band unhidden, allowing us to detect the thermal emission of the disc and to estimate the mass of the black hole. The non-thermal emission is well reproduced by a one-zone leptonic model. The non-thermal radiative processes are synchrotron, self-Compton and external Compton using seed photons from both the broad-line region (BLR) and the torus. We find that our data are better reproduced if we assume that the location of the dissipation region of the jet, R_diss, is in-between the torus, (at R_torus), and the BLR (R_torus>R_diss>R_BLR). The observed variability is explained by changing a minimum number of model parameters by a very small amount.

preprint2015arXiv

PSR J1906+0722: An Elusive Gamma-ray Pulsar

We report the discovery of PSR J1906+0722, a gamma-ray pulsar detected as part of a blind survey of unidentified Fermi Large Area Telescope (LAT) sources being carried out on the volunteer distributed computing system, Einstein@Home. This newly discovered pulsar previously appeared as the most significant remaining unidentified gamma-ray source without a known association in the second Fermi-LAT source catalog (2FGL) and was among the top ten most significant unassociated sources in the recent third catalog (3FGL). PSR J1906+0722 is a young, energetic, isolated pulsar, with a spin frequency of $8.9$ Hz, a characteristic age of $49$ kyr, and spin-down power $1.0 \times 10^{36}$ erg s$^{-1}$. In 2009 August it suffered one of the largest glitches detected from a gamma-ray pulsar ($Δf / f \approx 4.5\times10^{-6}$). Remaining undetected in dedicated radio follow-up observations, the pulsar is likely radio-quiet. An off-pulse analysis of the gamma-ray flux from the location of PSR J1906+0722 revealed the presence of an additional nearby source, which may be emission from the interaction between a neighboring supernova remnant and a molecular cloud. We discuss possible effects which may have hindered the detection of PSR J1906+0722 in previous searches and describe the methods by which these effects were mitigated in this survey. We also demonstrate the use of advanced timing methods for estimating the positional, spin and glitch parameters of difficult-to-time pulsars such as this.

preprint2015arXiv

Rapid Variability of Blazar 3C 279 during Flaring States in 2013-2014 with Joint Fermi-LAT, NuSTAR, Swift, and Ground-Based Multi-wavelength Observations

We report the results of a multi-band observing campaign on the famous blazar 3C 279 conducted during a phase of increased activity from 2013 December to 2014 April, including first observations of it with NuSTAR. The $γ$-ray emission of the source measured by Fermi-LAT showed multiple distinct flares reaching the highest flux level measured in this object since the beginning of the Fermi mission, with $F(E > 100\,{\rm MeV})$ of $10^{-5}$ photons cm$^{-2}$ s$^{-1}$, and with a flux doubling time scale as short as 2 hours. The $γ$-ray spectrum during one of the flares was very hard, with an index of $Γ_γ= 1.7 \pm 0.1$, which is rarely seen in flat spectrum radio quasars. The lack of concurrent optical variability implies a very high Compton dominance parameter $L_γ/L_{\rm syn} > 300$. Two 1-day NuSTAR observations with accompanying Swift pointings were separated by 2 weeks, probing different levels of source activity. While the 0.5$-$70 keV X-ray spectrum obtained during the first pointing, and fitted jointly with Swift-XRT is well-described by a simple power law, the second joint observation showed an unusual spectral structure: the spectrum softens by $ΔΓ_{\rm X} \simeq 0.4$ at $\sim$4 keV. Modeling the broad-band SED during this flare with the standard synchrotron plus inverse Compton model requires: (1) the location of the $γ$-ray emitting region is comparable with the broad line region radius, (2) a very hard electron energy distribution index $p \simeq 1$, (3) total jet power significantly exceeding the accretion disk luminosity $L_{\rm j}/L_{\rm d} \gtrsim 10$, and (4) extremely low jet magnetization with $L_{\rm B}/L_{\rm j} \lesssim 10^{-4}$. We also find that single-zone models that match the observed $γ$-ray and optical spectra cannot satisfactorily explain the production of X-ray emission.

preprint2015arXiv

Search for Early Gamma-ray Production in Supernovae Located in a Dense Circumstellar Medium with the Fermi LAT

Supernovae (SNe) exploding in a dense circumstellar medium (CSM) are hypothesized to accelerate cosmic rays in collisionless shocks and emit GeV gamma rays and TeV neutrinos on a time scale of several months. We perform the first systematic search for gamma-ray emission in Fermi LAT data in the energy range from 100 MeV to 300 GeV from the ensemble of 147 SNe Type IIn exploding in dense CSM. We search for a gamma-ray excess at each SNe location in a one year time window. In order to enhance a possible weak signal, we simultaneously study the closest and optically brightest sources of our sample in a joint-likelihood analysis in three different time windows (1 year, 6 months and 3 months). For the most promising source of the sample, SN 2010jl (PTF10aaxf), we repeat the analysis with an extended time window lasting 4.5 years. We do not find a significant excess in gamma rays for any individual source nor for the combined sources and provide model-independent flux upper limits for both cases. In addition, we derive limits on the gamma-ray luminosity and the ratio of gamma-ray-to-optical luminosity ratio as a function of the index of the proton injection spectrum assuming a generic gamma-ray production model. Furthermore, we present detailed flux predictions based on multi-wavelength observations and the corresponding flux upper limit at 95% confidence level (CL) for the source SN 2010jl (PTF10aaxf).

preprint2015arXiv

Search for Gamma-Ray Emission from DES Dwarf Spheroidal Galaxy Candidates with Fermi-LAT Data

Due to their proximity, high dark-matter content, and apparent absence of non-thermal processes, Milky Way dwarf spheroidal satellite galaxies (dSphs) are excellent targets for the indirect detection of dark matter. Recently, eight new dSph candidates were discovered using the first year of data from the Dark Energy Survey (DES). We searched for gamma-ray emission coincident with the positions of these new objects in six years of Fermi Large Area Telescope data. We found no significant excesses of gamma-ray emission. Under the assumption that the DES candidates are dSphs with dark matter halo properties similar to the known dSphs, we computed individual and combined limits on the velocity-averaged dark matter annihilation cross section for these new targets. If the estimated dark-matter content of these dSph candidates is confirmed, they will constrain the annihilation cross section to lie below the thermal relic cross section for dark matter particles with masses < 20 GeV annihilating via the b-bbar or tau+tau- channels.

preprint2015arXiv

The 1st Fermi Lat Supernova Remnant Catalog

To uniformly determine the properties of supernova remnants (SNRs) at high energies, we have developed the first systematic survey at energies from 1 to 100 GeV using data from the Fermi Large Area Telescope. Based on the spatial overlap of sources detected at GeV energies with SNRs known from radio surveys, we classify 30 sources as likely GeV SNRs. We also report 14 marginal associations and 245 flux upper limits. A mock catalog in which the positions of known remnants are scrambled in Galactic longitude, allows us to determine an upper limit of 22% on the number of GeV candidates falsely identified as SNRs. We have also developed a method to estimate spectral and spatial systematic errors arising from the diffuse interstellar emission model, a key component of all Galactic Fermi LAT analyses. By studying remnants uniformly in aggregate, we measure the GeV properties common to these objects and provide a crucial context for the detailed modeling of individual SNRs. Combining our GeV results with multiwavelength (MW) data, including radio, X-ray, and TeV, demonstrates the need for improvements to previously sufficient, simple models describing the GeV and radio emission from these objects. We model the GeV and MW emission from SNRs in aggregate to constrain their maximal contribution to observed Galactic cosmic rays.

preprint2015arXiv

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

The third catalog of active galactic nuclei (AGNs) detected by the Fermi-LAT (3LAC) is presented. It is based on the third Fermi-LAT catalog (3FGL) of sources detected between 100 MeV and 300 GeV with a Test Statistic (TS) greater than 25, between 2008 August 4 and 2012 July 31. The 3LAC includes 1591 AGNs located at high Galactic latitudes (|b|>10°), a 71% increase over the second catalog based on 2 years of data. There are 28 duplicate associations, thus 1563 of the 2192 high-latitude gamma-ray sources of the 3FGL catalog are AGNs. Most of them (98%) are blazars. About half of the newly detected blazars are of unknown type, i.e., they lack spectroscopic information of sufficient quality to determine the strength of their emission lines. Based on their gamma-ray spectral properties, these sources are evenly split between flat-spectrum radio quasars (FSRQs) and BL~Lacs. The most abundant detected BL~Lacs are of the high-synchrotron-peaked (HSP) type. About 50% of the BL~Lacs have no measured redshifts. A few new rare outliers (HSP-FSRQs and high-luminosity HSP BL~Lacs) are reported. The general properties of the 3LAC sample confirm previous findings from earlier catalogs. The fraction of 3LAC blazars in the total population of blazars listed in BZCAT remains non-negligible even at the faint ends of the BZCAT-blazar radio, optical and X-ray flux distributions, which is a clue that even the faintest known blazars could eventually shine in gamma rays at LAT-detection levels. The energy-flux distributions of the different blazar populations are in good agreement with extrapolation from earlier catalogs.

preprint2014arXiv

Dark Matter Constraints from Observations of 25 Milky Way Satellite Galaxies with the Fermi Large Area Telescope

The dwarf spheroidal satellite galaxies of the Milky Way are some of the most dark-matter-dominated objects known. Due to their proximity, high dark matter content, and lack of astrophysical backgrounds, dwarf spheroidal galaxies are widely considered to be among the most promising targets for the indirect detection of dark matter via gamma rays. Here we report on gamma-ray observations of 25 Milky Way dwarf spheroidal satellite galaxies based on 4 years of Fermi Large Area Telescope (LAT) data. None of the dwarf galaxies are significantly detected in gamma rays, and we present gamma-ray flux upper limits between 500 MeV and 500 GeV. We determine the dark matter content of 18 dwarf spheroidal galaxies from stellar kinematic data and combine LAT observations of 15 dwarf galaxies to constrain the dark matter annihilation cross section. We set some of the tightest constraints to date on the the annihilation of dark matter particles with masses between 2 GeV and 10 TeV into prototypical Standard Model channels. We find these results to be robust against systematic uncertainties in the LAT instrument performance, diffuse gamma-ray background modeling, and assumed dark matter density profile.

preprint2014arXiv

Deep Broadband Observations of the Distant Gamma-ray Blazar PKS 1424+240

We present deep VERITAS observations of the blazar PKS 1424+240, along with contemporaneous Fermi Large Area Telescope, Swift X-ray Telescope and Swift UV Optical Telescope data between 2009 February 19 and 2013 June 8. This blazar resides at a redshift of $z\ge0.6035$, displaying a significantly attenuated gamma-ray flux above 100 GeV due to photon absorption via pair-production with the extragalactic background light. We present more than 100 hours of VERITAS observations from three years, a multiwavelength light curve and the contemporaneous spectral energy distributions. The source shows a higher flux of (2.1$\pm0.3$)$\times10^{-7}$ ph m$^{-2}$s$^{-1}$ above 120 GeV in 2009 and 2011 as compared to the flux measured in 2013, corresponding to (1.02$\pm0.08$)$\times10^{-7}$ ph m$^{-2}$s$^{-1}$ above 120 GeV. The measured differential very high energy (VHE; $E\ge100$ GeV) spectral indices are $Γ=$3.8$\pm$0.3, 4.3$\pm$0.6 and 4.5$\pm$0.2 in 2009, 2011 and 2013, respectively. No significant spectral change across the observation epochs is detected. We find no evidence for variability at gamma-ray opacities of greater than $τ=2$, where it is postulated that any variability would be small and occur on longer than year timescales if hadronic cosmic-ray interactions with extragalactic photon fields provide a secondary VHE photon flux. The data cannot rule out such variability due to low statistics.

preprint2014arXiv

Detection of significant cm to sub-mm band radio and gamma-ray correlated variability in Fermi bright blazars

The exact location of the gamma-ray emitting region in blazars is still controversial. In order to attack this problem we present first results of a cross-correlation analysis between radio (11 cm to 0.8 mm wavelength, F-GAMMA program) and gamma-ray (0.1-300 GeV) ~ 3.5 year light curves of 54 Fermi-bright blazars. We perform a source stacking analysis and estimate significances and chance correlations using mixed source correlations. Our results reveal: (i) the first highly significant multi-band radio and gamma-ray correlations (radio lagging gamma rays) when averaging over the whole sample, (ii) average time delays (source frame: 76+/-23 to 7+/-9 days), systematically decreasing from cm to mm/sub-mm bands with a frequency dependence tau_r,gamma (nu) ~ nu^-1, in good agreement with jet opacity dominated by synchrotron self-absorption, (iii) a bulk gamma-ray production region typically located within/upstream of the 3 mm core region (tau_3mm,gamma=12+/-8 days), (iv) mean distances between the region of gamma-ray peak emission and the radio "tau=1 photosphere" decreasing from 9.8+/-3.0 pc (11 cm) to 0.9+/-1.1 pc (2 mm) and 1.4+/-0.8 pc (0.8 mm), (v) 3 mm/gamma-ray correlations in 9 individual sources at a significance level where one is expected by chance (probability: 4 times 10^-6), (vi) opacity and "time lag core shift" estimates for quasar 3C 454.3 providing a lower limit for the distance of the bulk gamma-ray production region from the supermassive black hole (SMBH) of ~ 0.8-1.6 pc, i.e. at the outer edge of the Broad Line Region (BLR) or beyond. A 3 mm tau=1 surface at ~ 2-3 pc from the jet-base (i.e. well outside the "canonical BLR") finally suggests that BLR material extends to several pc distances from the SMBH.

preprint2014arXiv

Fermi-LAT Detection of Gravitational Lens Delayed Gamma-ray Flares from Blazar B0218+357

Using data from the Fermi Large Area Telescope (LAT), we report the first clear gamma-ray measurement of a delay between flares from the gravitationally lensed images of a blazar. The delay was detected in B0218+357, a known double-image lensed system, during a period of enhanced gamma-ray activity with peak fluxes consistently observed to reach >20-50 times its previous average flux. An auto-correlation function analysis identified a delay in the gamma-ray data of 11.46 +/- 0.16 days (1 sigma) that is ~1 day greater than previous radio measurements. Considering that it is beyond the capabilities of the LAT to spatially resolve the two images, we nevertheless decomposed individual sequences of superposing gamma-ray flares/delayed emissions. In three such ~8-10 day-long sequences within a ~4-month span, considering confusion due to overlapping flaring emission and flux measurement uncertainties, we found flux ratios consistent with ~1, thus systematically smaller than those from radio observations. During the first, best-defined flare, the delayed emission was detailed with a Fermi pointing, and we observed flux doubling timescales of ~3-6 hrs implying as well extremely compact gamma-ray emitting regions.

preprint2014arXiv

Inferred cosmic-ray spectrum from ${\it Fermi}$-LAT $γ$-ray observations of the Earth's limb

Recent accurate measurements of cosmic-ray (CR) species by ATIC-2, CREAM, and PAMELA reveal an unexpected hardening in the proton and He spectra above a few hundred GeV, a gradual softening of the spectra just below a few hundred GeV, and a harder spectrum of He compared to that of protons. These newly-discovered features may offer a clue to the origin of high-energy CRs. We use the ${\it Fermi}$ Large Area Telescope observations of the $γ$-ray emission from the Earth's limb for an indirect measurement of the local spectrum of CR protons in the energy range $\sim 90~$GeV-$6~$TeV (derived from a photon energy range $15~$GeV-$1~$TeV). Our analysis shows that single power law and broken power law spectra fit the data equally well and yield a proton spectrum with index $2.68 \pm 0.04$ and $2.61 \pm 0.08$ above $\sim 200~$GeV, respectively.

preprint2014arXiv

Multifrequency Studies of the Peculiar Quasar 4C +21.35 During the 2010 Flaring Activity

The discovery of rapidly variable Very High Energy (VHE; E > 100 GeV) gamma-ray emission from 4C +21.35 (PKS 1222+216) by MAGIC on 2010 June 17, triggered by the high activity detected by the Fermi Large Area Telescope (LAT) in high energy (HE; E > 100 MeV) gamma-rays, poses intriguing questions on the location of the gamma-ray emitting region in this flat spectrum radio quasar. We present multifrequency data of 4C +21.35 collected from centimeter to VHE during 2010 to investigate the properties of this source and discuss a possible emission model. The first hint of detection at VHE was observed by MAGIC on 2010 May 3, soon after a gamma-ray flare detected by Fermi-LAT that peaked on April 29. The same emission mechanism may therefore be responsible for both the HE and VHE emission during the 2010 flaring episodes. Two optical peaks were detected on 2010 April 20 and June 30, close in time but not simultaneous with the two gamma-ray peaks, while no clear connection was observed between the X-ray an gamma-ray emission. An increasing flux density was observed in radio and mm bands from the beginning of 2009, in accordance with the increasing gamma-ray activity observed by Fermi-LAT, and peaking on 2011 January 27 in the mm regime (230 GHz). We model the spectral energy distributions (SEDs) of 4C +21.35 for the two periods of the VHE detection and a quiescent state, using a one-zone model with the emission coming from a very compact region outside the broad line region. The three SEDs can be fit with a combination of synchrotron self-Compton and external Compton emission of seed photons from a dust torus, changing only the electron distribution parameters between the epochs. The fit of the optical/UV part of the spectrum for 2010 April 29 seems to favor an inner disk radius of <6 gravitational radii, as one would expect from a prograde-rotating Kerr black hole.

preprint2014arXiv

Search for cosmic-ray induced gamma-ray emission in Galaxy Clusters

Current theories predict relativistic hadronic particle populations in clusters of galaxies in addition to the already observed relativistic leptons. In these scenarios hadronic interactions give rise to neutral pions which decay into $γ$ rays, that are potentially observable with the Large Area Telescope (LAT) on board the Fermi space telescope. We present a joint likelihood analysis searching for spatially extended $γ$-ray emission at the locations of 50 galaxy clusters in 4 years of Fermi-LAT data under the assumption of the universal cosmic-ray model proposed by Pinzke & Pfrommer (2010). We find an excess at a significance of $2.7σ$ which upon closer inspection is however correlated to individual excess emission towards three galaxy clusters: Abell 400, Abell 1367 and Abell 3112. We discuss these cases in detail and conservatively attribute the emission to unmodeled background (for example, radio galaxies within the clusters). Through the combined analysis of 50 clusters we exclude hadronic injection efficiencies in simple hadronic models above 21% and establish limits on the cosmic-ray to thermal pressure ratio within the virial radius, $R_{200}$, to be below 1.2-1.4% depending on the morphological classification. In addition we derive new limits on the $γ$-ray flux from individual clusters in our sample.

preprint2014arXiv

The spectrum of isotropic diffuse gamma-ray emission between 100 MeV and 820 GeV

The γ-ray sky can be decomposed into individually detected sources, diffuse emission attributed to the interactions of Galactic cosmic rays with gas and radiation fields, and a residual all-sky emission component commonly called the isotropic diffuse γ-ray background (IGRB). The IGRB comprises all extragalactic emissions too faint or too diffuse to be resolved in a given survey, as well as any residual Galactic foregrounds that are approximately isotropic. The first IGRB measurement with the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope (Fermi) used 10 months of sky-survey data and considered an energy range between 200 MeV and 100 GeV. Improvements in event selection and characterization of cosmic-ray backgrounds, better understanding of the diffuse Galactic emission, and a longer data accumulation of 50 months, allow for a refinement and extension of the IGRB measurement with the LAT, now covering the energy range from 100 MeV to 820 GeV. The IGRB spectrum shows a significant high-energy cutoff feature, and can be well described over nearly four decades in energy by a power law with exponential cutoff having a spectral index of $2.32\pm0.02$ and a break energy of $(279\pm52)$ GeV using our baseline diffuse Galactic emission model. The total intensity attributed to the IGRB is $(7.2\pm0.6) \times 10^{-6}$ cm$^{-2}$ s$^{-1}$ sr$^{-1}$ above 100 MeV, with an additional $+15$%/$-30$% systematic uncertainty due to the Galactic diffuse foregrounds.

preprint2013arXiv

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

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.

preprint2013arXiv

Constraints on the Galactic Population of TEV Pulsar Wind Nebulae Using Fermi Large Area Telescope Observations

Pulsar wind nebulae (PWNe) have been established as the most populous class of TeV gamma-ray emitters. Since launch, the Fermi Large Area Telescope (LAT)identified five high-energy (100MeV <E< 100 GeV) gamma-ray sources as PWNe, and detected a large number of PWNe candidates, all powered by young and energetic pulsars. The wealth of multi-wavelength data available and the new results provided by Fermi-LAT give us an opportunity to find new PWNe and to explore the radiative processes taking place in known ones. The TeV gamma-ray unidentifiedsources (UNIDs) are the best candidates for finding new PWNe. Using 45 months of Fermi-LAT data for energies above 10 GeV, an analysis was performed near the position of 58TeV PWNe and UNIDs within 5deg of the Galactic Plane to establish new constraints on PWNe properties and find new clues on the nature of UNIDs. Of the 58 sources, 30 were detected, and this work provides their gamma-rayfluxes for energies above 10 GeV. The spectral energy distributions (SED) andupper limits, in the multi-wavelength context, also provide new information on the source nature and can help distinguish between emission scenarios, i.e. between classification as a pulsar candidate or as a PWN candidate. Six new GeV PWNe candidates are described in detail and compared with existing models. A population study of GeV PWNe candidates as a function of the pulsar/PWN system characteristics is presented.

preprint2013arXiv

Detection of the Characteristic Pion-Decay Signature in Supernova Remnants

Cosmic rays are particles (mostly protons) accelerated to relativistic speeds. Despite wide agreement that supernova remnants (SNRs) are the sources of galactic cosmic rays, unequivocal evidence for the acceleration of protons in these objects is still lacking. When accelerated protons encounter interstellar material, they produce neutral pions, which in turn decay into gamma rays. This offers a compelling way to detect the acceleration sites of protons. The identification of pion-decay gamma rays has been difficult because high-energy electrons also produce gamma rays via bremsstrahlung and inverse Compton scattering. We detected the characteristic pion-decay feature in the gamma-ray spectra of two SNRs, IC 443 and W44, with the Fermi Large Area Telescope. This detection provides direct evidence that cosmic-ray protons are accelerated in SNRs.

preprint2013arXiv

First Results from NuSTAR Observations of Mkn 421

Mkn 421 is a nearby active galactic nucleus dominated at all wavelengths by a very broad non-thermal continuum thought to arise from a relativistic jet seen at a small angle to the line of sight. Its spectral energy distribution peaks in the X-ray and TeV gamma-ray bands, where the energy output is dominated by cooling of high-energy electrons in the jet. In order to study the electron distribution and its evolution, we carried out a dedicated multi-wavelength campaign, including extensive observations by the recently launched highly sensitive hard X-ray telescope NuSTAR, between December 2012 and May 2013. Here we present some initial results based on NuSTAR data from January through March 2013, as well as calibration observations conducted in 2012. Although the observations cover some of the faintest hard X-ray flux states ever observed for Mkn 421, the sensitivity is high enough to resolve intra-day spectral variability. We find that in this low state the dominant flux variations are smooth on timescales of hours, with typical intra-hour variations of less than 5%. We do not find evidence for either a cutoff in the hard X-ray spectrum, or a rise towards a high-energy component, but rather that at low flux the spectrum assumes a power law shape with a photon index of approximately 3. The spectrum is found to harden with increasing brightness.

preprint2013arXiv

The First Pulse of the Extremely Bright GRB 130427A: A Test Lab for Synchrotron Shocks

Gamma-ray burst (GRB) 130427A is one of the most energetic GRBs ever observed. The initial pulse up to 2.5 s is possibly the brightest well-isolated pulse observed to date. A fine time resolution spectral analysis shows power-law decays of the peak energy from the onset of the pulse, consistent with models of internal synchrotron shock pulses. However, a strongly correlated power-law behavior is observed between the luminosity and the spectral peak energy that is inconsistent with curvature effects arising in the relativistic outflow. It is difficult for any of the existing models to account for all of the observed spectral and temporal behaviors simultaneously.

preprint2012arXiv

Binary Millisecond Pulsar Discovery via Gamma-Ray Pulsations

Millisecond pulsars, old neutron stars spun-up by accreting matter from a companion star, can reach high rotation rates of hundreds of revolutions per second. Until now, all such "recycled" rotation-powered pulsars have been detected by their spin-modulated radio emission. In a computing-intensive blind search of gamma-ray data from the Fermi Large Area Telescope (with partial constraints from optical data), we detected a 2.5-millisecond pulsar, PSR J1311-3430. This unambiguously explains a formerly unidentified gamma-ray source that had been a decade-long enigma, confirming previous conjectures. The pulsar is in a circular orbit with an orbital period of only 93 minutes, the shortest of any spin-powered pulsar binary ever found.

preprint2012arXiv

Fermi LAT Search for Dark Matter in Gamma-ray Lines and the Inclusive Photon Spectrum

Dark matter particle annihilation or decay can produce monochromatic gamma-ray lines and contribute to the diffuse gamma-ray background. Flux upper limits are presented for gamma-ray spectral lines from 7 to 200 GeV and for the diffuse gamma-ray background from 4.8 GeV to 264 GeV obtained from two years of Fermi Large Area Telescope data integrated over most of the sky. We give cross section upper limits and decay lifetime lower limits for dark matter models that produce gamma-ray lines or contribute to the diffuse spectrum, including models proposed as explanations of the PAMELA and Fermi cosmic-ray data.

preprint2012arXiv

Gamma-ray Activity in the Crab Nebula: The Exceptional Flare of April 2011

The Large Area Telescope on board the Fermi satellite observed a gamma-ray flare in the Crab nebula lasting for approximately nine days in April of 2011. The source, which at optical wavelengths has a size of ~11 ly across, doubled its gamma-ray flux within eight hours. The peak photon flux was (186 +- 6) 10-7 cm-2 s-1 above 100 MeV, which corresponds to a 30-fold increase compared to the average value. During the flare, a new component emerged in the spectral energy distribution, which peaked at an energy of (375 +- 26) MeV at flare maximum. The observations imply that the emission region was likely relativistically beamed toward us and that variations in its motion are responsible for the observed spectral variability.

preprint2012arXiv

Gamma-ray observations of the Orion Molecular Clouds with the Fermi Large Area Telescope

We report on the gamma-ray observations of giant molecular clouds Orion A and B with the Large Area Telescope (LAT) on-board the Fermi Gamma-ray Space Telescope. The gamma-ray emission in the energy band between \sim100 MeV and \sim100 GeV is predicted to trace the gas mass distribution in the clouds through nuclear interactions between the Galactic cosmic rays (CRs) and interstellar gas. The gamma-ray production cross-section for the nuclear interaction is known to \sim10% precision which makes the LAT a powerful tool to measure the gas mass column density distribution of molecular clouds for a known CR intensity. We present here such distributions for Orion A and B, and correlate them with those of the velocity integrated CO intensity (WCO) at a 1° \times1° pixel level. The correlation is found to be linear over a WCO range of ~10 fold when divided in 3 regions, suggesting penetration of nuclear CRs to most of the cloud volumes. The Wco-to-mass conversion factor, Xco, is found to be \sim2.3\times10^20 cm-2(K km s-1)-1 for the high-longitude part of Orion A (l > 212°), \sim1.7 times higher than \sim1.3 \times 10^20 found for the rest of Orion A and B. We interpret the apparent high Xco in the high-longitude region of Orion A in the light of recent works proposing a non-linear relation between H2 and CO densities in the diffuse molecular gas. Wco decreases faster than the H2 column density in the region making the gas "darker" to Wco.

preprint2012arXiv

GeV Observations of Star-forming Galaxies with \textit{Fermi} LAT

Recent detections of the starburst galaxies M82 and NGC 253 by gamma-ray telescopes suggest that galaxies rapidly forming massive stars are more luminous at gamma-ray energies compared to their quiescent relatives. Building upon those results, we examine a sample of 69 dwarf, spiral, and luminous and ultraluminous infrared galaxies at photon energies 0.1-100 GeV using 3 years of data collected by the Large Area Telescope (LAT) on the \textit{Fermi Gamma-ray Space Telescope} (\textit{Fermi}). Measured fluxes from significantly detected sources and flux upper limits for the remaining galaxies are used to explore the physics of cosmic rays in galaxies. We find further evidence for quasi-linear scaling relations between gamma-ray luminosity and both radio continuum luminosity and total infrared luminosity which apply both to quiescent galaxies of the Local Group and low-redshift starburst galaxies (conservative $P$-values $\lesssim0.05$ accounting for statistical and systematic uncertainties). The normalizations of these scaling relations correspond to luminosity ratios of $\log(L_{0.1-100 \rm{GeV}}/L_{1.4 \rm{GHz}}) = 1.7 \pm 0.1_{\rm (statistical)} \pm 0.2_{\rm (dispersion)}$ and $\log(L_{0.1-100 \rm{GeV}}/L_{8-1000 μ\rm{m}}) = -4.3 \pm 0.1_{\rm (statistical)} \pm 0.2_{\rm (dispersion)}$ for a galaxy with a star formation rate of 1 $M_{\odot}$ yr$^{-1}$, assuming a Chabrier initial mass function. Using the relationship between infrared luminosity and gamma-ray luminosity, the collective intensity of unresolved star-forming galaxies at redshifts $0<z<2.5$ above 0.1 GeV is estimated to be 0.4-2.4 $\times 10^{-6}$ ph cm$^{-2}$ s$^{-1}$ sr$^{-1}$ (4-23% of the intensity of the isotropic diffuse component measured with the LAT). We anticipate that $\sim10$ galaxies could be detected by their cosmic-ray induced gamma-ray emission during a 10-year \textit{Fermi} mission.

preprint2012arXiv

Measurement of separate cosmic-ray electron and positron spectra with the Fermi Large Area Telescope

We measured separate cosmic-ray electron and positron spectra with the Fermi Large Area Telescope. Because the instrument does not have an onboard magnet, we distinguish the two species by exploiting the Earth's shadow, which is offset in opposite directions for opposite charges due to the Earth's magnetic field. We estimate and subtract the cosmic-ray proton background using two different methods that produce consistent results. We report the electron-only spectrum, the positron-only spectrum, and the positron fraction between 20 GeV and 200 GeV. We confirm that the fraction rises with energy in the 20-100 GeV range. The three new spectral points between 100 and 200 GeV are consistent with a fraction that is continuing to rise with energy.

preprint2012arXiv

Multi-wavelength observations of blazar AO 0235+164 in the 2008-2009 flaring state

The blazar AO 0235+164 (z = 0.94) has been one of the most active objects observed by Fermi Large Area Telescope (LAT) since its launch in Summer 2008. In addition to the continuous coverage by Fermi, contemporaneous observations were carried out from the radio to γ -ray bands between 2008 September and 2009 February. In this paper, we summarize the rich multi-wavelength data collected during the campaign (including F-GAMMA, GASP- WEBT, Kanata, OVRO, RXTE, SMARTS, Swift, and other instruments), examine the cross-correlation between the light curves measured in the different energy bands, and interpret the resulting spectral energy distributions in the context of well-known blazar emission models. We find that the γ -ray activity is well correlated with a series of near-IR/optical flares, accompanied by an increase in the optical polarization degree. On the other hand, the X-ray light curve shows a distinct 20 day high state of unusually soft spectrum, which does not match the extrapolation of the optical/UV synchrotron spectrum. We tentatively interpret this feature as the bulk Compton emission by cold electrons contained in the jet, which requires an accretion disk corona with an effective covering factor of 19% at a distance of 100 Rg . We model the broadband spectra with a leptonic model with external radiation dominated by the infrared emission from the dusty torus.

preprint2012arXiv

Search for Dark Matter Satellites using the FERMI-LAT

Numerical simulations based on the Lambda-CDM model of cosmology predict a large number of as yet unobserved Galactic dark matter satellites. We report the results of a Large Area Telescope (LAT) search for these satellites via the gamma-ray emission expected from the annihilation of weakly interacting massive particle (WIMP) dark matter. Some dark matter satellites are expected to have hard gamma-ray spectra, finite angular extents, and a lack of counterparts at other wavelengths. We sought to identify LAT sources with these characteristics, focusing on gamma-ray spectra consistent with WIMP annihilation through the $b \bar b$ channel. We found no viable dark matter satellite candidates using one year of data, and we present a framework for interpreting this result in the context of numerical simulations to constrain the velocity-averaged annihilation cross section for a conventional 100 GeV WIMP annihilating through the $b \bar b$ channel.

preprint2012arXiv

The structure and emission model of the relativistic jet in the quasar 3C 279 inferred from radio to high-energy gamma-ray observations in 2008-2010

We present time-resolved broad-band observations of the quasar 3C 279 obtained from multi-wavelength campaigns conducted during the first two years of the Fermi Gamma-ray Space Telescope mission. While investigating the previously reported gamma-ray/optical flare accompanied by a change in optical polarization, we found that the optical emission appears delayed with respect to the gamma-ray emission by about 10 days. X-ray observations reveal a pair of `isolated' flares separated by ~90 days, with only weak gamma-ray/optical counterparts. The spectral structure measured by Spitzer reveals a synchrotron component peaking in the mid-infrared band with a sharp break at the far-infrared band during the gamma-ray flare, while the peak appears in the mm/sub-mm band in the low state. Selected spectral energy distributions are fitted with leptonic models including Comptonization of external radiation produced in a dusty torus or the broad-line region. Adopting the interpretation of the polarization swing involving propagation of the emitting region along a curved trajectory, we can explain the evolution of the broad-band spectra during the gamma-ray flaring event by a shift of its location from ~ 1 pc to ~ 4 pc from the central black hole. On the other hand, if the gamma-ray flare is generated instead at sub-pc distance from the central black hole, the far-infrared break can be explained by synchrotron self-absorption. We also model the low spectral state, dominated by the mm/sub-mm peaking synchrotron component, and suggest that the corresponding inverse-Compton component explains the steady X-ray emission.

preprint2011arXiv

Constraining Dark Matter Models from a Combined Analysis of Milky Way Satellites with the Fermi Large Area Telescope

Satellite galaxies of the Milky Way are among the most promising targets for dark matter searches in gamma rays. We present a search for dark matter consisting of weakly interacting massive particles, applying a joint likelihood analysis to 10 satellite galaxies with 24 months of data of the Fermi Large Area Telescope. No dark matter signal is detected. Including the uncertainty in the dark matter distribution, robust upper limits are placed on dark matter annihilation cross sections. The 95% confidence level upper limits range from about 1e-26 cm^3 s^-1 at 5 GeV to about 5e-23 cm^3 s^-1 at 1 TeV, depending on the dark matter annihilation final state. For the first time, using gamma rays, we are able to rule out models with the most generic cross section (~3e-26 cm^3 s^-1 for a purely s-wave cross section), without assuming additional boost factors.

preprint2011arXiv

Fermi Large Area Telescope Detection of Bright Gamma-ray Outbursts from a Peculiar Quasar 4C +21.35

In this paper we report on the two-year-long Fermi-LAT observation of the peculiar blazar 4C +21.35 (PKS 1222+216). This source was in a quiescent state from the start of science operations of the Fermi Gamma-ray Space Telescope in 2008 August until 2009 September, and then became more active, with gradually increasing flux and some moderately-bright flares. In 2010 April and June, 4C +21.35 underwent a very strong GeV outburst composed of several major flares characterized by rise and decay timescales of the order of a day. During the outburst, the GeV spectra of 4C +21.35 displayed a broken power-law form with spectral breaks observed near 1-3 GeV photon energies. We demonstrate that, at least during the major flares, the jet in 4C +21.35 carried a total kinetic luminosity comparable to the total accretion power available to feed the outflow. We also discuss the origin of the break observed in the flaring spectra of 4C +21.35. We show that, in principle, a model involving annihilation of the GeV photons on the He II Lyman recombination continuum and line emission of "broad line region" clouds may account for such. However, we also discuss the additional constraint provided by the detection of 4C +21.35 at 0.07-0.4 TeV energies by the MAGIC telescope, which coincided with one of the GeV flares of the source. We argue that there are reasons to believe that the $lesssim$,TeV emission of 4C +21.35 (as well as the GeV emission of the source, if co-spatial), is not likely to be produced inside the broad line region zone of highest ionization ($sim 10^{17}$,cm from the nucleus), but instead originates further away from the active center, namely around the characteristic scale of the hot dusty torus surrounding the 4C +21.35 nucleus ($sim 10^{19}$,cm).

preprint2011arXiv

Simultaneous multi-wavelength campaign on PKS 2005-489 in a high state

The high-frequency peaked BL Lac object PKS 2005-489 was the target of a multi-wavelength campaign with simultaneous observations in the TeV gamma-ray (H.E.S.S.), GeV gamma-ray (Fermi/LAT), X-ray (RXTE, Swift), UV (Swift) and optical (ATOM, Swift) bands. This campaign was carried out during a high flux state in the synchrotron regime. The flux in the optical and X-ray bands reached the level of the historical maxima. The hard GeV spectrum observed with Fermi/LAT connects well to the very high energy (VHE, E>100GeV) spectrum measured with H.E.S.S. with a peak energy between ~5 and 500 GeV. Compared to observations with contemporaneous coverage in the VHE and X-ray bands in 2004, the X-ray flux was ~50 times higher during the 2009 campaign while the TeV gamma-ray flux shows marginal variation over the years. The spectral energy distribution during this multi-wavelength campaign was fit by a one zone synchrotron self-Compton model with a well determined cutoff in X-rays. The parameters of a one zone SSC model are inconsistent with variability time scales. The variability behaviour over years with the large changes in synchrotron emission and small changes in the inverse Compton emission does not warrant an interpretation within a one-zone SSC model despite an apparently satisfying fit to the broadband data in 2009.

preprint2011arXiv

The Discovery and Nature of Optical Transient CSS100217:102913+404220

We report on the discovery and observations of the extremely luminous optical transient CSS100217:102913+404220 (CSS100217 hereafter). Spectroscopic observations show this transient was coincident with a galaxy at redshift z=0.147, and reached an apparent magnitude of V ~ 16.3. After correcting for foreground Galactic extinction we determine the absolute magnitude to be M_V =-22.7 approximately 45 days after maximum light. Based on our unfiltered optical photometry the peak optical emission was L = 1.3 x 10^45 erg s^-1, and over a period of 287 rest-frame days had an integrated bolometric luminosity of 1.2 x 10^52 erg. Analysis of the pre-outburst SDSS spectrum of the source shows features consistent with a Narrow-line Seyfert1 (NLS1) galaxy. High-resolution HST and Keck followup observations show the event occurred within 150pc of nucleus of the galaxy, suggesting a possible link to the active nuclear region. However, the rapid outburst along with photometric and spectroscopic evolution are much more consistent with a luminous supernova. Line diagnostics suggest that the host galaxy is undergoing significant star formation. We use extensive follow-up of the event along with archival CSS and SDSS data to investigate the three most likely sources of such an event; 1) an extremely luminous supernova; 2) the tidal disruption of a star by the massive nuclear black hole; 3) variability of the central AGN. We find that CSS100217 was likely an extremely luminous type IIn supernova that occurred within range of the narrow-line region of an AGN. We discuss how similar events may have been missed in past supernova surveys because of confusion with AGN activity.

preprint2010arXiv

Constraints on Cosmological Dark Matter Annihilation from the Fermi-LAT Isotropic Diffuse Gamma-Ray Measurement

The first published Fermi large area telescope (Fermi-LAT) measurement of the isotropic diffuse gamma-ray emission is in good agreement with a single power law, and is not showing any signature of a dominant contribution from dark matter sources in the energy range from 20 to 100 GeV. We use the absolute size and spectral shape of this measured flux to derive cross section limits on three types of generic dark matter candidates: annihilating into quarks, charged leptons and monochromatic photons. Predicted gamma-ray fluxes from annihilating dark matter are strongly affected by the underlying distribution of dark matter, and by using different available results of matter structure formation we assess these uncertainties. We also quantify how the dark matter constraints depend on the assumed conventional backgrounds and on the Universe's transparency to high-energy gamma-rays. In reasonable background and dark matter structure scenarios (but not in all scenarios we consider) it is possible to exclude models proposed to explain the excess of electrons and positrons measured by the Fermi-LAT and PAMELA experiments. Derived limits also start to probe cross sections expected from thermally produced relics (e.g. in minimal supersymmetry models) annihilating predominantly into quarks. For the monochromatic gamma-ray signature, the current measurement constrains only dark matter scenarios with very strong signals.

preprint2010arXiv

Fermi LAT Search for Photon Lines from 30 to 200 GeV and Dark Matter Implications

Dark matter (DM) particle annihilation or decay can produce monochromatic $γ$-rays readily distinguishable from astrophysical sources. $γ$-ray line limits from 30 GeV to 200 GeV obtained from 11 months of Fermi Large Area Space Telescope data from 20-300 GeV are presented using a selection based on requirements for a $γ$-ray line analysis, and integrated over most of the sky. We obtain $γ$-ray line flux upper limits in the range $0.6-4.5\times 10^{-9}\mathrm{cm}^{-2}\mathrm{s}^{-1}$, and give corresponding DM annihilation cross-section and decay lifetime limits. Theoretical implications are briefly discussed.

preprint2010arXiv

Radio and Gamma-Ray Constraints on the Emission Geometry and Birthplace of PSR J2043+2740

We report on the first year of Fermi gamma-ray observations of pulsed high-energy emission from the old PSR J2043+2740. The study of the gamma-ray efficiency of such old pulsars gives us an insight into the evolution of pulsars' ability to emit in gammma rays as they age. The gamma-ray lightcurve of this pulsar above 0.1 GeV is clearly defined by two sharp peaks, 0.353+/-0.035 periods apart. We have combined the gamma-ray profile characteristics of PSR J2043+2740 with the geometrical properties of the pulsar's radio emission, derived from radio polarization data, and constrained the pulsar-beam geometry in the framework of a Two Pole Caustic and an Outer Gap model. The ranges of magnetic inclination and viewing angle were determined to be {alpha,zeta}~{52-57,61-68} for the Two Pole Caustic model, and {alpha,zeta}~{62-73,74-81} and {alpha,zeta}~{72-83,60-75} for the Outer Gap model. Based on this geometry, we assess possible birth locations for this pulsar and derive a likely proper motion, sufficiently high to be measurable with VLBI. At a characteristic age of 1.2 Myr, PSR J2043+2740 is the third oldest of all discovered, non-recycled, gamma-ray pulsars: it is twice as old as the next oldest, PSR J0357+32, and younger only than the recently discovered PSR J1836+5925 and PSR J2055+25, both of which are at least 5 and 10 times less energetic, respectively.

preprint2010arXiv

The Discovery of gamma-Ray Emission From The Blazar RGB J0710+591

The high-frequency-peaked BL Lacertae object RGB J0710+591 was observed in the very high-energy (VHE; E > 100 GeV) wave band by the VERITAS array of atmospheric Cherenkov telescopes. The observations, taken between 2008 December and 2009 March and totaling 22.1 hr, yield the discovery of VHE gamma rays from the source. RGB J0710+591 is detected at a statistical significance of 5.5 standard deviations (5.5σ) above the background, corresponding to an integral flux of (3.9 +/- 0.8) x 10-12 cm-2 s-1 (3% of the Crab Nebula's flux) above 300 GeV. The observed spectrum can be fit by a power law from 0.31 to 4.6 TeV with a photon spectral index of 2.69 +/- 0.26stat +/- 0.20sys. These data are complemented by contemporaneous multiwavelength data from the Fermi Large Area Telescope, the Swift X-ray Telescope, the Swift Ultra-Violet and Optical Telescope, and the Michigan-Dartmouth-MIT observatory. Modeling the broadband spectral energy distribution (SED) with an equilibrium synchrotron self-Compton model yields a good statistical fit to the data. The addition of an external-Compton component to the model does not improve the fit nor brings the system closer to equipartition. The combined Fermi and VERITAS data constrain the properties of the high-energy emission component of the source over 4 orders of magnitude and give measurements of the rising and falling sections of the SED.

preprint2009arXiv

Identification and properties of the photospheric emission in GRB090902B

The Fermi Gamma-ray Space Telescope observed the bright and long GRB090902B, lying at a redshift of z = 1.822. Together the Large Area Telescope (LAT) and the Gamma-ray Burst Monitor (GBM) cover the spectral range from 8 keV to >300 GeV is covered. Here we show that the prompt burst spectrum is consistent with emission from the jet photosphere combined with non-thermal emission described by a single power-law with photon index -1.9. The photosphere gives rise to a strong quasi-blackbody spectrum which is somewhat broader than a single Planck function and has a characteristic temperature of ~290keV. We derive the photospheric radius Rph = (1.1 \pm 0.3) x 10^12 Y^{1/4} cm and the bulk Lorentz factor of the flow, which is found to vary by a factor of two and has a maximal value of Gamma = 750 Y^{1/4}. Here Y is the ratio between the total fireball energy and the energy emitted in the gamma-rays. We find that during the first quarter of the prompt phase the photospheric emission dominates, which explains the delayed onset of the observed flux in the LAT compared to the GBM. We model the photospheric emission with a multi-color blackbody and its shape indicates that the photospheric radius increases at higher latitudes. We interpret the broad band emission as synchrotron emission at R ~ 4x10^15 cm. Our analysis emphasize the importance of having high temporal resolution when performing spectral analysis on GRBs, since there is strong spectral evolution.

preprint2009arXiv

Prospects for GRB Science with the Fermi Large Area Telescope

The LAT instrument on the Fermi mission will reveal the rich spectral and temporal gamma-ray burst phenomena in the > 100 MeV band. The synergy with Fermi's GBM detectors will link these observations to those in the well explored 10-1000 keV range; the addition of the > 100 MeV band observations will resolve theoretical uncertainties about burst emission in both the prompt and afterglow phases. Trigger algorithms will be applied to the LAT data both onboard the spacecraft and on the ground. The sensitivity of these triggers will differ because of the available computing resources onboard and on the ground. Here we present the LAT's burst detection methodologies and the instrument's GRB capabilities.

preprint2009arXiv

Source detection using a 3D sparse representation: application to the Fermi gamma-ray space telescope

The multiscale variance stabilization Transform (MSVST) has recently been proposed for Poisson data denoising. This procedure, which is nonparametric, is based on thresholding wavelet coefficients. We present in this paper an extension of the MSVST to 3D data (in fact 2D-1D data) when the third dimension is not a spatial dimension, but the wavelength, the energy, or the time. We show that the MSVST can be used for detecting and characterizing astrophysical sources of high-energy gamma rays, using realistic simulated observations with the Large Area Telescope (LAT). The LAT was launched in June 2008 on the Fermi Gamma-ray Space Telescope mission. The MSVST algorithm is very fast relative to traditional likelihood model fitting, and permits efficient detection across the time dimension and immediate estimation of spectral properties. Astrophysical sources of gamma rays, especially active galaxies, are typically quite variable, and our current work may lead to a reliable method to quickly characterize the flaring properties of newly-detected sources.

preprint2009arXiv

The Evolution of Swift/BAT blazars and the origin of the MeV background

We use 3 years of data from the Swift/BAT survey to select a complete sample of X-ray blazars above 15 keV. This sample comprises 26 Flat-Spectrum Radio Quasars (FSRQs) and 12 BL Lac objects detected over a redshift range of 0.03<z<4.0. We use this sample to determine, for the first time in the 15--55 keV band, the evolution of blazars. We find that, contrary to the Seyfert-like AGNs detected by BAT, the population of blazars shows strong positive evolution. This evolution is comparable to the evolution of luminous optical QSOs and luminous X-ray selected AGNs. We also find evidence for an epoch-dependence of the evolution as determined previously for radio-quiet AGNs. We interpret both these findings as a strong link between accretion and jet activity. In our sample, the FSRQs evolve strongly, while our best-fit shows that BL Lacs might not evolve at all. The blazar population accounts for 10-20 % (depending on the evolution of the BL Lacs) of the Cosmic X-ray background (CXB) in the 15--55 keV band. We find that FSRQs can explain the entire CXB emission for energies above 500 keV solving the mystery of the generation of the MeV background. The evolution of luminous FSRQs shows a peak in redshift ($z_c$=4.3$\pm0.5$) which is larger than the one observed in QSOs and X-ray selected AGNs. We argue that FSRQs can be used as tracers of massive elliptical galaxies in the early Universe.

preprint2009arXiv

The Spectral Energy Distribution of Fermi bright blazars

(Abridged) We have conducted a detailed investigation of the broad-band spectral properties of the γ-ray selected blazars of the Fermi LAT Bright AGN Sample (LBAS). By combining our accurately estimated Fermi gamma-ray spectra with Swift, radio, infra-red, optical and other hard X-ray/gamma-ray data, collected within three months of the LBAS data taking period, we were able to assemble high-quality and quasi-simultaneous Spectral Energy Distributions (SED) for 48 LBAS blazars.The SED of these gamma-ray sources is similar to that of blazars discovered at other wavelengths, clearly showing, in the usual Log $ν$ - Log $ν$ F$_ν$ representation, the typical broad-band spectral signatures normally attributed to a combination of low-energy synchrotron radiation followed by inverse Compton emission of one or more components. We have used these SEDs to characterize the peak intensity of both the low and the high-energy components. The results have been used to derive empirical relationships that estimate the position of the two peaks from the broad-band colors (i.e. the radio to optical and optical to X-ray spectral slopes) and from the gamma-ray spectral index. Our data show that the synchrotron peak frequency $ν_p^S$ is positioned between 10$^{12.5}$ and 10$^{14.5}$ Hz in broad-lined FSRQs and between $10^{13}$ and $10^{17}$ Hz in featureless BL Lacertae objects.We find that the gamma-ray spectral slope is strongly correlated with the synchrotron peak energy and with the X-ray spectral index, as expected at first order in synchrotron - inverse Compton scenarios. However, simple homogeneous, one-zone, Synchrotron Self Compton (SSC) models cannot explain most of our SEDs, especially in the case of FSRQs and low energy peaked (LBL) BL Lacs. (...)

preprint1998arXiv

Synchrotron and SSC Emission and the Blast-Wave Model of Gamma-Ray Bursts

We investigate the dynamics and radiation from a relativistic blast-wave which decelerates as it sweeps up ambient matter. The bulk kinetic energy of the blast-wave shell is converted into internal energy by the process of accreting external matter. If it takes the form of non-thermal electrons and magnetic fields, then this internal energy will be emitted as synchrotron and synchrotron self-Compton radiation. We perform analytic and numerical calculations for the deceleration and radiative processes and present time-resolved spectra throughout the evolution of the blast-wave. We also examine the dependence of the burst spectra and light curves on various parameters describing the magnetic field and non-thermal electron distributions. We find that for bursts such as GRB~910503, GRB~910601 and GRB~910814, the spectral shapes of the prompt gamma-ray emission at the peaks in $νF_ν$ strongly constrain the magnetic fields in these bursts to be well below ($\la 10^{-2}$) the equipartition values. These calculations are also considered in the context of the afterglow emission from the recently detected gamma-ray burst counterparts.