The accretion-powered X-ray pulsar GX 301-2 was observed with the balloon-borne X-Calibur hard X-ray polarimeter during late December 2018, with contiguous observations by the NICER X-ray telescope, the Swift X-ray Telescope and Burst Alert Telescope, and the Fermi Gamma-ray Burst Monitor spanning several months. The observations detected the pulsar in a rare apastron flaring state coinciding with a significant spin-up of the pulsar discovered with the Fermi GBM. The X-Calibur, NICER, and Swift observations reveal a pulse profile strongly dominated by one main peak, and the NICER and Swift data show strong variation of the profile from pulse to pulse. The X-Calibur observations constrain for the first time the linear polarization of the 15-35 keV emission from a highly magnetized accreting neutron star, indicating a polarization degree of (27+38-27)% (90% confidence limit) averaged over all pulse phases. We discuss the spin-up and the X-ray spectral and polarimetric results in the context of theoretical predictions. We conclude with a discussion of the scientific potential of future observations of highly magnetized neutron stars with the more sensitive follow-up mission XL-Calibur.

It is commonly assumed that in black hole accretion disks the angular momenta of the disk and the black hole are aligned. However, for a significant fraction of stellar mass black holes and supermassive black holes, the momenta may not be aligned. In such systems, the interplay of disk viscosity and general relativistic frame dragging can cause the disk to warp or break into two (or more) distinct planes; this is called the Bardeen-Petterson effect. We have developed a general relativistic ray-tracing code to find the energy spectra and polarization of warped accretion disks, accounting for the emission from the disk and for photons reflecting one or multiple times off the warped accretion disk segments. We find that polarization angle can be used to give a lower limit on the misalignment angle when a previous measurement of the jet, which is thought be aligned with the black hole angular momentum, can be spatially resolved.

Since the birth of X-ray Astronomy, spectacular advances have been seen in the imaging, spectroscopic and timing studies of the hot and violent X-ray Universe, and further leaps forward are expected in the future. On the other hand, polarimetry is very much lagging behind: after the measurements of the Crab Nebula and Scorpius X-1, obtained by OSO-8 in the 70s, no more observations have been performed in the classical X-ray band, even if some interesting results have been obtained in hard X-rays and in soft gamma-rays. The NASA/ASI mission IXPE, scheduled for the launch in 2021, is going to provide for the first time imaging X-ray polarimetry in the 2-8 keV band thanks to its photoelectric polarimeter, coupled with ~25'' angular resolution X-ray mirrors. Its orders of magnitude improvement in sensitivity with respect to the OSO-8 Bragg polarimeter implies scientifically meaningful polarimetric measurements for at least the brightest specimens of most classes of X-ray sources. In 2027, the Chinese-led mission eXTP should also be launched. In addition to timing and spectroscopic instruments, eXTP will have on board photoelectric polarimeters very similar to those of IXPE, but with a total effective area 2-3 times larger. Building on IXPE results, eXTP will increase the number of sources for which significant polarimetric measurements could be obtained. However, further progresses, such as exploring a broader energy range, considering a larger effective area, improving the angular resolution, and performing wide-field polarization measurements, are needed to reach a mature phase for X-ray polarimetry. In the first part of this White Paper we will discuss a few scientific cases in which a next generation X-ray Polarimetry mission can provide significant advances. In the second part, a possible concept for a medium-class Next Generation X-ray Polarimetry (NGXP) mission will be sketched.

The X-ray Polarization Probe (XPP) is a second generation X-ray polarimeter following up on the Imaging X-ray Polarimetry Explorer (IXPE). The XPP will offer true broadband polarimetery over the wide 0.2-60 keV bandpass in addition to imaging polarimetry from 2-8 keV. The extended energy bandpass and improvements in sensitivity will enable the simultaneous measurement of the polarization of several emission components. These measurements will give qualitatively new information about how compact objects work, and will probe fundamental physics, i.e. strong-field quantum electrodynamics and strong gravity.

In this paper, we present the analysis and results of a direct measurement of the cosmic-ray proton spectrum with the CALET instrument onboard the International Space Station, including the detailed assessment of systematic uncertainties. The observation period used in this analysis is from October 13, 2015 to August 31, 2018 (1054 days). We have achieved the very wide energy range necessary to carry out measurements of the spectrum from 50 GeV to 10 TeV covering, for the first time in space, with a single instrument the whole energy interval previously investigated in most cases in separate subranges by magnetic spectrometers (BESS-TeV, PAMELA, and AMS-02) and calorimetric instruments (ATIC, CREAM, and NUCLEON). The observed spectrum is consistent with AMS-02 but extends to nearly an order of magnitude higher energy, showing a very smooth transition of the power-law spectral index from -2.81 +- 0.03 (50--500 GeV) neglecting solar modulation effects (or -2.87 +- 0.06 including solar modulation effects in the lower energy region) to -2.56 +- 0.04 (1--10 TeV), thereby confirming the existence of spectral hardening and providing evidence of a deviation from a single power law by more than 3 sigma.

This white paper highlights compact object and fundamental physics science opportunities afforded by high-throughput broadband (0.1-60 keV) X-ray polarization observations. X-ray polarimetry gives new observables with geometric information about stellar remnants which are many orders of magnitude too small for direct imaging. The X-ray polarimetric data also reveal details about the emission mechanisms and the structure of the magnetic fields in and around the most extreme objects in the Universe. Whereas the Imaging X-ray Polarimetry Explorer (IXPE) to be launched in 2021 will obtain first results for bright objects, a follow-up mission could be one order of magnitude more sensitive and would be able to use a broader bandpass to perform physics type experiments for representative samples of sources.

Current X-ray observations and simulations show that gravitational lensing can be used to infer the structure near the event horizons of black holes, constrain the dynamics and evolution of black-hole accretion and outflows, test general relativity in the strong-gravity regime and place constraints on the evolution of dark matter in the lensing galaxies. These science goals currently cannot be achieved in a statistically large sample of z = 0.5 - 5 lensed quasars due to the limited capabilities of current X-ray telescopes and the relatively low number (~200) of known lensed quasars. The latter limitation will be resolved with the multi-band and wide-field photometric optical survey of LSST that is expected to lead to the discovery of > 4,000 additional gravitationally lensed systems. As we show in this white paper, these science goals can be reached with an X-ray telescope having a spatial resolution of <0.5arcsec to resolve the lensed images and a collecting area of >0.5 m^2 at 1 keV.

The CALorimetric Electron Telescope (CALET) is a high-energy astroparticle physics space experiment installed on the International Space Station (ISS), developed and operated by Japan in collaboration with Italy and the United States. The CALET mission goals include the investigation of possible nearby sources of high-energy electrons, of the details of galactic particle acceleration and propagation, and of potential signatures of dark matter. CALET measures the cosmic-ray electron + positron flux up to 20 TeV, gamma-rays up to 10 TeV, and nuclei with Z=1 to 40 up to 1,000 TeV for the more abundant elements during a long-term observation aboard the ISS. Starting science operation in mid-October 2015, CALET performed continuous observation without major interruption with close to 20 million triggered events over 10 GeV per month. Based on the data taken during the first two-years, we present an overview of CALET observations: uses w/o major interruption 1) Electron + positron energy spectrum, 2) Nuclei analysis, 3) Gamma-ray observation including a characterization of on-orbit performance. Results of the electromagnetic counterpart search for LIGO/Virgo gravitational wave events are discussed as well.

The Chandra observations of several gravitationally lensed quasars show evidence for flux and spectral variability of the X-ray emission that is uncorrelated between images and is thought to result from the microlensing by stars in the lensing galaxy. We report here on the most detailed modeling of such systems to date, including simulations of the emission of the Fe K-alpha fluorescent radiation from the accretion disk with a general relativistic ray tracing code, the use of realistic microlensing magnification maps derived from inverse ray shooting calculations, and the simulation of the line detection biases. We use lensing and black hole parameters appropriate for the quadruply lensed quasar RX J1131-1231, and compare the simulated results with the observational results. The simulations cannot fully reproduce the distribution of the detected line energies indicating that some of the assumptions underlying the simulations are not correct, or that the simulations are missing some important physics. We conclude by discussing several possible explanations.

X-ray studies of stellar mass black holes in X-ray binaries and mass-accreting supermassive black holes in Active Galactic Nuclei have achieved a high degree of maturity and have delivered detailed information about the astrophysical sources and the physics of black hole accretion. In this article, I review recent progress made towards using the X-ray observations for testing the "Kerr hypothesis" that the background spacetimes of all astrophysical quasi-stationary black holes are described by the Kerr metric. Although the observations have indeed revealed clear evidence for relativistic effects in strong-field gravity, quantitative tests of the Kerr hypothesis still struggle with theoretical and practical difficulties. In this article, I describe several recently introduced test metrics and review the status of constraining the background spacetimes of mass accreting stellar mass and supermassive black holes with these test metrics. The main conclusion of the discussion is that astrophysical uncertainties are large compared to the rather small observational differences between the Kerr and non-Kerr metrics precluding quantitative constraints on deviations from the Kerr metric at this point in time. I conclude with discussing future progress enabled by more detailed numerical simulations and by future X-ray spectroscopy, timing, polarimetry, and interferometry missions.

We present results on searches for gamma-ray counterparts of the LIGO/Virgo gravitational-wave events using CALorimetric Electron Telescope ({\sl CALET}) observations. The main instrument of {\sl CALET}, CALorimeter (CAL), observes gamma-rays from $\sim1$ GeV up to 10 TeV with a field of view of nearly 2 sr. In addition, the {\sl CALET} gamma-ray burst monitor (CGBM) views $\sim$3 sr and $\sim2\pi$ sr of the sky in the 7 keV -- 1 MeV and the 40 keV -- 20 MeV bands, respectively, by using two different crystal scintillators. The {\sl CALET} observations on the International Space Station started in October 2015, and here we report analyses of events associated with the following gravitational wave events: GW151226, GW170104, GW170608, GW170814 and GW170817. Although only upper limits on gamma-ray emission are obtained, they correspond to a luminosity of $10^{49}\sim10^{53}$ erg s$^{-1}$ in the GeV energy band depending on the distance and the assumed time duration of each event, which is approximately the order of luminosity of typical short gamma-ray bursts. This implies there will be a favorable opportunity to detect high-energy gamma-ray emission in further observations if additional gravitational wave events with favorable geometry will occur within our field-of-view. We also show the sensitivity of {\sl CALET} for gamma-ray transient events which is the order of $10^{-7}$~erg\,cm$^{-2}$\,s$^{-1}$ for an observation of 100~s duration.

Extended results on the cosmic-ray electron + positron spectrum from 11 GeV to 4.8 TeV are presented based on observations with the Calorimetric Electron Telescope (CALET) on the International Space Station utilizing the data up to November 2017. The analysis uses the full detector acceptance at high energies, approximately doubling the statistics compared to the previous result. CALET is an all-calorimetric instrument with a total thickness of 30 $X_0$ at normal incidence and fine imaging capability, designed to achieve large proton rejection and excellent energy resolution well into the TeV energy region. The observed energy spectrum in the region below 1 TeV shows good agreement with Alpha Magnetic Spectrometer (AMS-02) data. In the energy region below $\sim$300 GeV, CALET's spectral index is found to be consistent with the AMS-02, Fermi Large Area Telescope (Fermi-LAT) and Dark Matter Particle Explorer (DAMPE), while from 300 to 600 GeV the spectrum is significantly softer than the spectra from the latter two experiments. The absolute flux of CALET is consistent with other experiments at around a few tens of GeV. However, it is lower than those of DAMPE and Fermi-LAT with the difference increasing up to several hundred GeV. The observed energy spectrum above $\sim$1 TeV suggests a flux suppression consistent within the errors with the results of DAMPE, while CALET does not observe any significant evidence for a narrow spectral feature in the energy region around 1.4 TeV. Our measured all-electron flux, including statistical errors and a detailed breakdown of the systematic errors, is tabulated in the Supplemental Material in order to allow more refined spectral analyses based on our data.

The CALorimetric Electron Telescope (CALET), launched for installation on the International Space Station (ISS) in August, 2015, has been accumulating scientific data since October, 2015. CALET is intended to perform long-duration observations of high-energy cosmic rays onboard the ISS. CALET directly measures the cosmic-ray electron spectrum in the energy range of 1 GeV to 20 TeV with a 2% energy resolution above 30 GeV. In addition, the instrument can measure the spectrum of gamma rays well into the TeV range, and the spectra of protons and nuclei up to a PeV. In order to operate the CALET onboard ISS, JAXA Ground Support Equipment (JAXA-GSE) and the Waseda CALET Operations Center (WCOC) have been established. Scientific operations using CALET are planned at WCOC, taking into account orbital variations of geomagnetic rigidity cutoff. Scheduled command sequences are used to control the CALET observation modes on orbit. Calibration data acquisition by, for example, recording pedestal and penetrating particle events, a low-energy electron trigger mode operating at high geomagnetic latitude, a low-energy gamma-ray trigger mode operating at low geomagnetic latitude, and an ultra heavy trigger mode, are scheduled around the ISS orbit while maintaining maximum exposure to high-energy electrons and other high-energy shower events by always having the high-energy trigger mode active. The WCOC also prepares and distributes CALET flight data to collaborators in Italy and the United States. As of August 31, 2017, the total observation time is 689 days with a live time fraction of the total time of approximately 84%. Nearly 450 million events are collected with a high-energy (E>10 GeV) trigger. By combining all operation modes with the excellent-quality on-orbit data collected thus far, it is expected that a five-year observation period will provide a wealth of new and interesting results.

First results of a cosmic-ray electron + positron spectrum, from 10 GeV to 3 TeV, is presented based upon observations with the CALET instrument on the ISS starting in October, 2015. Nearly a half million electron + positron events are included in the analysis. CALET is an all-calorimetric instrument with total vertical thickness of 30 $X_0$ and a fine imaging capability designed to achieve a large proton rejection and excellent energy resolution well into the TeV energy region. The observed energy spectrum over 30 GeV can be fit with a single power law with a spectral index of -3.152 $\pm$ 0.016 (stat.+ syst.). Possible structure observed above 100 GeV requires further investigation with increased statistics and refined data analysis.

In August 2015, the CALorimetric Electron Telescope (CALET), designed for long exposure observations of high energy cosmic rays, docked with the International Space Station (ISS) and shortly thereafter began tocollect data. CALET will measure the cosmic ray electron spectrum over the energy range of 1 GeV to 20 TeV with a very high resolution of 2% above 100 GeV, based on a dedicated instrument incorporating an exceptionally thick 30 radiation-length calorimeter with both total absorption and imaging (TASC and IMC) units. Each TASC readout channel must be carefully calibrated over the extremely wide dynamic range of CALET that spans six orders of magnitude in order to obtain a degree of calibration accuracy matching the resolution of energy measurements. These calibrations consist of calculating the conversion factors between ADC units and energy deposits, ensuring linearity over each gain range, and providing a seamless transition between neighboring gain ranges. This paper describes these calibration methods in detail, along with the resulting data and associated accuracies. The results presented in this paper show that a sufficient accuracy was achieved for the calibrations of each channel in order to obtain a suitable resolution over the entire dynamic range of the electron spectrum measurement.

Hard X-rays observed in Active Galactic Nuclei (AGNs) are thought to originate from the Comptonization of the optical/UV accretion disk photons in a hot corona. Polarization studies of these photons can help to constrain the corona geometry and the plasma properties. We have developed a ray-tracing code that simulates the Comptonization of accretion disk photons in coronae of arbitrary shape, and use it here to study the polarization of the X-ray emission from wedge and spherical coronae. We study the predicted polarization signatures for the fully relativistic and various approximate treatments of the elemental Compton scattering processes. We furthermore use the code to evaluate the impact of non-thermal electrons and cyclo-synchrotron photons on the polarization properties. Finally, we model the NuSTAR observations of the Seyfert I galaxy Mrk 335 and predict the associated polarization signal. Our studies show that X-ray polarimetry missions such as NASA's Imaging X-ray Polarimetry Explorer (IXPE) and the X-ray Imaging Polarimetry Explorer (XIPE) proposed to ESA will provide valuable new information about the physical properties of the plasma close to the event horizon of AGN black holes.

Observations of the fluorescent Fe K-alpha emission line from the inner accretion flows of stellar mass black holes in X-ray binaries and supermassive black holes in Active Galactic Nuclei have become an important tool to study the magnitude and inclination of the black hole spin, and the structure of the accretion flow close to the event horizon of the black hole. Modeling spectral, timing, and soon also X-ray polarimetric observations of the Fe K-alpha emission requires to calculate the specific intensity in the rest frame of the emitting plasma. We revisit the derivation of the equation used for calculating the illumination of the accretion disk by the corona. We present an alternative derivation leading to a simpler equation, and discuss the relation to the previously published results.

Infrared, optical and ultraviolet spectropolarimetric observations have proven to be ideal tools for the study of the hidden nuclei of type-2 active galactic nuclei (AGN) and for constraining the composition and morphology of the sub-parsec scale emission components. In this paper, we extend the analysis to the polarization of the X-rays from type-2 AGN. Combining two radiative transfer codes, we performed the first simulations of photons originating in the gravity dominated vicinity of the black hole and scattering in structures all the way out to the parsec-scale torus and polar winds. We demonstrate that, when strong gravity effects are accounted for, the X-ray polarimetric signal of Seyfert-2s carries as much information about the central AGN components as spectropolarimetric observations of Seyfert-1s. The spectropolarimetric measurements can constrain the spin of the central supermassive black hole even in edge-on AGN, the hydrogen column density along the observer's line-of-sight, and the composition of the polar outflows. However, the polarization state of the continuum source is washed out by multiple scattering, and should not be measurable unless the initial polarization is exceptionally strong. Finally, we estimate that modern X-ray polarimeters, either based on the photo-electric effect or on Compton scattering, will require long observational times on the order of a couple of mega-seconds to be able to properly measure the polarization of type-2 AGN.

The analysis of the Chandra X-ray observations of the gravitationally lensed quasar RX J1131-1231 revealed the detection of multiple and energy-variable spectral peaks. The spectral variability is thought to result from the microlensing of the Fe K-alpha emission, selectively amplifying the emission from certain regions of the accretion disk with certain effective frequency shifts of the Fe K-alpha line emission. In this paper, we combine detailed simulations of the emission of Fe K-alpha photons from the accretion disk of a Kerr black hole with calculations of the effect of gravitational microlensing on the observed energy spectra. The simulations show that microlensing can indeed produce multiply peaked energy spectra. We explore the dependence of the spectral characteristics on black hole spin, accretion disk inclination, corona height, and microlensing amplification factor, and show that the measurements can be used to constrain these parameters. We find that the range of observed spectral peak energies of QSO RX J1131-1231 can only be reproduced for black hole inclinations exceeding 70 degree and for lamppost corona heights of less than 30 gravitational radii above the black hole. We conclude by emphasizing the scientific potential of studies of the microlensed Fe K$\alpha$ quasar emission and the need for more detailed modeling that explores how the results change for more realistic accretion disk and corona geometries and microlensing magnification patterns. A full analysis should furthermore model the signal-to-noise ratio of the observations and the resulting detection biases.

Lorentz invariance is the fundamental symmetry of Einstein's theory of special relativity, and has been tested to great level of detail. However, theories of quantum gravity at the Planck scale indicate that Lorentz symmetry may be broken at that scale motivating further tests. While the Planck energy is currently unreachable by experiment, tiny residual effects at attainable energies can become measurable when photons propagate over sufficiently large distances. The Standard-Model Extension (SME) is an effective field theory approach to describe low-energy effects of quantum gravity theories. Lorentz and CPT symmetry violating effects are introduced by adding additional terms to the Standard Model Lagrangian. These terms can be ordered by the mass dimension of the corresponding operator, and the leading terms of interest have dimension d = 5. Effects of these operators are a linear variation of the speed of light with photon energy, and a rotation of the linear polarization of photons quadratic in photon energy, as well as anisotropy. We analyzed optical polarization data from 72 AGN and GRBs and derived the first set of limits on all 16 coefficients of mass dimension d = 5 of the SME photon sector. Our constraints imply a lower limit on the energy scale of quantum gravity of $10^6$ times the Planck energy, severly limiting the phase space for any theory that predicts a rotation of the photon polarization quadratic in energy.

We present a promising new technique, the g-distribution method, for measuring the inclination angle (i), the innermost stable circular orbit (ISCO), and the spin of a supermassive black hole. The g-distribution method uses measurements of the energy shifts in the relativistic iron line emitted by the accretion disk of a supermassive black hole due to microlensing by stars in a foreground galaxy relative to the g-distribution shifts predicted from microlensing caustic calculations. We apply the method to the gravitationally lensed quasars RX J1131-1231 (z_s=0.658, z_l=0.295), QJ 0158-4325 (z_s=1.294, z_l=0.317), and SDSS 1004+4112 (z_s=1.734, z_l=0.68). For RX J1131-1231 our initial results indicate that r_ISCO<8.5 gravitational radii (r_g) and i > 76 degrees. We detect two shifted Fe lines, in several observations, as predicted in our numerical simulations of caustic crossings. The current DeltaE-distribution of RX J1131-1231 is sparsely sampled but further X-ray monitoring of RX J1131-1231 and other lensed quasars will provide improved constraints on the inclination angles, ISCO radii and spins of the black holes of distant quasars.

X-ray polarimetry has seen a growing interest in recent years. Improvements in detector technology and focusing X-ray optics now enable sensitive astrophysical X-ray polarization measurements. These measurements will provide new insights into the processes at work in accreting black holes, the emission of X-rays from neutron stars and magnetars, and the structure of AGN jets. X-Calibur is a balloon-borne hard X-ray scattering polarimeter. An X-ray mirror with a focal length of 8 m focuses X-rays onto the detector, which consists of a plastic scattering element surrounded by Cadmium-Zinc-Telluride detectors, which absorb and record the scattered X-rays. Since X-rays preferentially scatter perpendicular to their polarization direction, the polarization properties of an X-ray beam can be inferred from the azimuthal distribution of scattered X-rays. A close alignment of the X-ray focal spot with the center of the detector is required in order to reduce systematic uncertainties and to maintain a high photon detection efficiency. This places stringent requirements on the mechanical and thermal stability of the telescope structure. During the flight on a stratospheric balloon, X-Calibur makes use of the Wallops Arc-Second Pointer (WASP) to point the telescope at astrophysical sources. In this paper, we describe the design, construction, and test of the telescope structure, as well as its performance during a 25-hour flight from Ft. Sumner, New Mexico. The carbon fiber-aluminum composite structure met the requirements set by X-Calibur and its design can easily be adapted for other types of experiments, such as X-ray imaging or spectroscopic telescopes.

We present upper limits in the hard X-ray and gamma-ray bands at the time of the LIGO gravitational-wave event GW 151226 derived from the CALorimetric Electron Telescope (CALET) observation. The main instrument of CALET, CALorimeter (CAL), observes gamma-rays from ~1 GeV up to 10 TeV with a field of view of ~2 sr. The CALET gamma-ray burst monitor (CGBM) views ~3 sr and ~2pi sr of the sky in the 7 keV - 1 MeV and the 40 keV - 20 MeV bands, respectively, by using two different scintillator-based instruments. The CGBM covered 32.5% and 49.1% of the GW 151226 sky localization probability in the 7 keV - 1 MeV and 40 keV - 20 MeV bands respectively. We place a 90% upper limit of 2 x 10^{-7} erg cm-2 s-1 in the 1 - 100 GeV band where CAL reaches 15% of the integrated LIGO probability (~1.1 sr). The CGBM 7 sigma upper limits are 1.0 x 10^{-6} erg cm-2 s-1 (7-500 keV) and 1.8 x 10^{-6} erg cm-2 s-1 (50-1000 keV) for one second exposure. Those upper limits correspond to the luminosity of 3-5 x 10^{49} erg s-1 which is significantly lower than typical short GRBs.

Lorentz invariance, the fundamental symmetry of Einstein's theory of Special Relativity, has been established and tested by many classical and modern experiments. However, many theories that unify the Standard Model of particle physics and General Relativity predict a violation of Lorentz invariance at the Planck scale. While this energy range cannot be reached by current experiments, minute deviations from Lorentz symmetry may be present at lower energies. Astrophysical experiments are very suitable to search for these deviations, since their effects accumulate as photons travel across large distances. In this paper, we describe astrophysical methods that we used to constrain the photon dispersion and vacuum birefringence.

The BL Lacertae object 1ES 1440+122 was observed in the energy range from 85 GeV to 30 TeV by the VERITAS array of imaging atmospheric Cherenkov telescopes. The observations, taken between 2008 May and 2010 June and totalling 53 hours, resulted in the discovery of $\gamma$-ray emission from the blazar, which has a redshift $z$=0.163. 1ES 1440+122 is detected at a statistical significance of 5.5 standard deviations above the background with an integral flux of (2.8$\pm0.7_{\mathrm{stat}}\pm0.8_{\mathrm{sys}}$) $\times$ 10$^{-12}$ cm$^{-2}$ s$^{-1}$ (1.2\% of the Crab Nebula's flux) above 200 GeV. The measured spectrum is described well by a power law from 0.2 TeV to 1.3 TeV with a photon index of 3.1 $\pm$ 0.4$_{\mathrm{stat}}$ $\pm$ 0.2$_{\mathrm{sys}}$. Quasi-simultaneous multi-wavelength data from the Fermi Large Area Telescope (0.3--300 GeV) and the Swift X-ray Telescope (0.2--10 keV) are additionally used to model the properties of the emission region. A synchrotron self-Compton model produces a good representation of the multi-wavelength data. Adding an external-Compton or a hadronic component also adequately describes the data.

The very high energy (VHE; E $>$ 100 GeV) blazar Markarian 501 was observed between April 17 and May 5 (MJD 54938--54956), 2009, as part of an extensive multi-wavelength campaign from radio to VHE. Strong VHE $\gamma$-ray activity was detected on May 1st with Whipple and VERITAS, when the flux (E $>$ 400 GeV) increased to 10 times the pre-flare baseline flux ($3.9{\times 10^{-11}}~{\rm ph~cm^{-2}~s^{-1}}$), reaching five times the flux of the Crab Nebula. This coincided with a decrease in the optical polarization and a rotation of the polarization angle by 15$^{\circ}$. This VHE flare showed a fast flux variation with an increase of a factor $\sim$4 in 25 minutes, and a falling time of $\sim$50 minutes. We present the observations of the quiescent state previous to the flare and of the high state after the flare, focusing on the flux and spectral variability from Whipple, VERITAS, Fermi-LAT, RXTE, and Swift combined with optical and radio data.

Observations with RXTE (Rossi X-ray Timing Explorer) revealed the presence of High Frequency Quasi-Periodic Oscillations (HFQPOs) of the X-ray flux from several accreting stellar mass Black Holes. HFQPOs (and their counterparts at lower frequencies) may allow us to study general relativity in the strong gravity regime. However, the observational evidence today does not yet allow us to distinguish between different HFQPO models. In this paper we use a general relativistic ray-tracing code to investigate X-ray timing-spectroscopy and polarization properties of HFQPOs in the orbiting Hotspot model. We study observational signatures for the particular case of the 166 Hz quasi-periodic oscillation (QPO) in the galactic binary GRS 1915+105. We conclude with a discussion of the observability of spectral signatures with a timing-spectroscopy experiment like the LOFT (Large Observatory for X-ray Timing) and polarization signatures with space-borne X-ray polarimeters such as IXPE (Imaging X-ray Polarimetry Explorer), PolSTAR (Polarization Spectroscopic Telescope Array), PRAXyS (Polarimetry of Relativistic X-ray Sources), or XIPE (X-ray Imaging Polarimetry Explorer). A high count-rate mission like LOFT would make it possible to get a QPO phase for each photon, enabling the study of the QPO-phase-resolved spectral shape and the correlation between this and the flux level. Owing to the short periods of the HFQPOs, first-generation X-ray polarimeters would not be able to assign a QPO phase to each photon. The study of QPO-phase-resolved polarization energy spectra would thus require simultaneous observations with a first-generation X-ray polarimeter and a LOFT-type mission.

Despite its success in the weak gravity regime, General Relativity (GR) has yet to be verified in the regime of strong gravity. In this paper, we present the results of detailed ray tracing simulations aiming at clarifying if the combined information from X-ray spectroscopy, timing, and polarization observations of stellar mass and supermassive black holes can be used to test GR's no-hair theorem. The latter states that stationary astrophysical black holes are described by the Kerr-family of metrics with the black hole mass and spin being the only free parameters. We use four "non-Kerr metrics", some phenomenological in nature and others motivated by alternative theories of gravity, and study the observational signatures of deviations from the Kerr metric. Particular attention is given to the case when all the metrics are set to give the same Innermost Stable Circular Orbit (ISCO) in quasi-Boyer Lindquist coordinates. We give a detailed discussion of similarities and differences of the observational signatures predicted for BHs in the Kerr metric and the non-Kerr metrics. We emphasize that even though some regions of the parameter space are nearly degenerate even when combining the information from all observational channels, X-ray observations of very rapidly spinning black holes can be used to exclude large regions of the parameter space of the alternative metrics. Although it proves difficult to distinguish between the Kerr and non-Kerr metrics for some portions of the parameter space, the observations of very rapidly spinning black holes like Cyg X-1 can be used to rule out large regions for several black hole metrics.

This paper describes the Polarization Spectroscopic Telescope Array (PolSTAR), a mission proposed to NASA's 2014 Small Explorer (SMEX) announcement of opportunity. PolSTAR measures the linear polarization of 3-50 keV (requirement; goal: 2.5-70 keV) X-rays probing the behavior of matter, radiation and the very fabric of spacetime under the extreme conditions close to the event horizons of black holes, as well as in and around magnetars and neutron stars. The PolSTAR design is based on the technology developed for the Nuclear Spectroscopic Telescope Array (NuSTAR) mission launched in June 2012. In particular, it uses the same X-ray optics, extendable telescope boom, optical bench, and CdZnTe detectors as NuSTAR. The mission has the sensitivity to measure ~1% linear polarization fractions for X-ray sources with fluxes down to ~5 mCrab. This paper describes the PolSTAR design as well as the science drivers and the potential science return.

While Lorentz invariance, the fundamental symmetry of Einstein's theory of General Relativity, has been tested to a great level of detail, Grand Unified Theories that combine gravity with the other three fundamental forces may result in a violation of Lorentz symmetry at the Planck scale. These energies are unattainable experimentally. However, minute deviations from Lorentz invariance may still be present at much lower energies. These deviations can accumulate over large distances, making astrophysical measurements the most sensitive tests of Lorentz symmetry. One effect of Lorentz invariance violation is an energy dependent photon dispersion of the vacuum resulting in differences of the light travel time from distant objects. The Standard-Model Extension (SME) is an effective theory to describe the low-energy behavior of a more fundamental Grand Unified Theory, including Lorentz and CPT violating terms. In the SME the Lorentz violating operators can in part be classified by their mass-dimension d, with the lowest order being d=5. However, measurements of photon polarization have constrained operators with d=5 setting lower limits on the energy at which they become dominant well beyond the Planck scale. On the other hand, these operators also violate CPT, and thus d=6 could be the leading order. In this paper we present constraints on all 25 real coefficients describing anisotropic non-birefringent Lorentz invariance violation at mass dimension d=6 in the SME. We used Fermi-LAT observations of 25 active galactic nuclei to constrain photon dispersion and combined our results with previously published limits in order to simultaneously constrain all 25 coefficients. This represents the first set of constraints on these coefficients of a mass-dimension d=6, whereas previous measurements were only able to constrain linear combinations of all 25 coefficients.

The VERITAS array of Cherenkov telescopes has carried out a deep observational program on the nearby dwarf spheroidal galaxy Segue 1. We report on the results of nearly 48 hours of good quality selected data, taken between January 2010 and May 2011. No significant $\gamma$-ray emission is detected at the nominal position of Segue 1, and upper limits on the integrated flux are derived. According to recent studies, Segue 1 is the most dark matter-dominated dwarf spheroidal galaxy currently known. We derive stringent bounds on various annihilating and decaying dark matter particle models. The upper limits on the velocity-weighted annihilation cross-section are $\mathrm{<\sigma v >^{95% CL} \lesssim 10^{-23} cm^{3} s^{-1}}$, improving our limits from previous observations of dwarf spheroidal galaxies by at least a factor of two for dark matter particle masses $\mathrm{m_{\chi}\gtrsim 300 GeV}$. The lower limits on the decay lifetime are at the level of $\mathrm{\tau^{95% CL} \gtrsim 10^{24} s}$. Finally, we address the interpretation of the cosmic ray lepton anomalies measured by ATIC and PAMELA in terms of dark matter annihilation, and show that the VERITAS observations of Segue 1 disfavor such a scenario.

It has long been recognised that quasi-periodic oscillations (QPOs) in the X-ray light curves of accreting black hole and neutron star binaries have the potential to be powerful diagnostics of strong field gravity. However, this potential cannot be fulfilled without a working theoretical model, which has remained elusive. Perhaps the most promising model associates the QPO with Lense-Thirring precession of the inner accretion flow, with the changes in viewing angle and Doppler boosting modulating the flux over the course of a precession cycle. Here, we consider the polarization signature of a precessing inner accretion flow. We use simple assumptions about the Comptonization process generating the emitted spectrum and take all relativistic effects into account, parallel transporting polarization vectors towards the observer along null geodesics in the Kerr metric. We find that both the degree of linear polarization and the polarization angle should be modulated on the QPO frequency. We calculate the predicted absolute rms variability amplitude of the polarization degree and angle for a specific model geometry. We find that it should be possible to detect these modulations for a reasonable fraction of parameter space with a future X-ray polarimeter such as NASA's Polarization Spectroscopic Telescope Array (PolSTAR: the satellite incarnation of the recent balloon experiment X-Calibur).

The relativistic jet in M87 offers a unique opportunity for understanding the detailed jet structure and emission processes due to its proximity. In particular, the peculiar jet region HST-1 at ~1 arcsecond (or 80 pc, projected) from the nucleus has attracted a great deal of interest in the last decade because of its superluminal motion and broadband radio-to-X-ray outbursts, which may be further connected to the gamma-ray productions up to TeV energies. Over the last five years, we have been doing an intensive monitoring of HST-1 with EVN at 5GHz in order to examine the detailed structural evolution and its possible connection to high-energy activities. While this program already yielded interesting results in terms of the detailed mas-scale structure, proper motion measurements and structural variations, the recent HST-1 brightness is continuously decreasing at this frequency. To counter this, we have shifted our monitoring frequency to 1.7GHz from October 2013. This strategy successfully recovered the fainter emission that was missed in the last 5GHz session. Moreover, we again discovered the sudden emergence of a new component at the upstream edge of HST-1, demonstrating that the use of EVN 1.7GHz is indeed powerful to probe the current weak nature of HST-1. Here we report early results from the 1.7GHz monitoring as well as further progress on the long-term kinematic study.

X-ray polarimetry promises to deliver unique information about the geometry of the inner accretion flow of astrophysical black holes and the nature of matter and electromagnetism in and around neutron stars. In this paper, we discuss the possibility to use Stokes parameters - a commonly used tool in radio, infrared, and optical polarimetry - to analyze the data from X-ray polarimeters such as scattering polarimeters and photoelectric effect polarimeters, which measure the linear polarization of the detected X-rays. Based on the azimuthal scattering angle (in the case of a scattering polarimeter) or the azimuthal component of the angle of the electron ejection (in the case of a photoelectric effect polarimeter), the Stokes parameters can be calculated for each event recorded in the detector. Owing to the additive nature of Stokes parameters, the analysis reduces to adding the Stokes parameters of the individual events and subtracting the Stokes parameters characterizing the background (if present). The main strength of this kind of analysis is that the errors on the Stokes parameters can be computed easily and are well behaved - in stark contrast of the errors on the polarization fraction and polarization direction. We demonstrate the power of the Stokes analysis by deriving several useful formulae, e.g. the expected error on the polarization fraction and polarization direction for a detection of $N_S$ signal and $N_{BG}$ background events, the optimal observation times of the signal and background regions in the presence of non-negligible background contamination of the signal, and the minimum detectable polarization (MDP) that can be achieved when following this prescription.

The nearby radio galaxy M87 offers a unique opportunity for exploring the connection between gamma-ray production and jet formation at an unprecedented linear resolution. However, the origin and location of the gamma-rays in this source is still elusive. Based on previous radio/TeV correlation events, the unresolved jet base (radio core) and the peculiar knot HST-1 at 120 pc from the nucleus are proposed as candidate site(s) of gamma-ray production. Here we report our intensive, high-resolution radio monitoring observations of the M87 jet with the VLBI Exploration of Radio Astrometry (VERA) and the European VLBI Network (EVN) from February 2011 to October 2012, together with contemporaneous high-energy gamma-ray light curves obtained by the Fermi Large Area Telescope. During this period, an elevated level of the M87 flux is reported at TeV with VERITAS. We detected a remarkable flux increase in the radio core with VERA at 22/43 GHz coincident with the VHE activity. Meanwhile, HST-1 remained quiescent in terms of its flux density and structure in the radio band. These results strongly suggest that the TeV gamma-ray activity in 2012 originates in the jet base within 0.03 pc (projected) from the central black hole.

We performed a 4.5-month multi-instrument campaign (from radio to VHE gamma rays) on Mrk421 between January 2009 and June 2009, which included VLBA, F-GAMMA, GASP-WEBT, Swift, RXTE, Fermi-LAT, MAGIC, and Whipple, among other instruments and collaborations. Mrk421 was found in its typical (non-flaring) activity state, with a VHE flux of about half that of the Crab Nebula, yet the light curves show significant variability at all wavelengths, the highest variability being in the X-rays. We determined the power spectral densities (PSD) at most wavelengths and found that all PSDs can be described by power-laws without a break, and with indices consistent with pink/red-noise behavior. We observed a harder-when-brighter behavior in the X-ray spectra and measured a positive correlation between VHE and X-ray fluxes with zero time lag. Such characteristics have been reported many times during flaring activity, but here they are reported for the first time in the non-flaring state. We also observed an overall anti-correlation between optical/UV and X-rays extending over the duration of the campaign. The harder-when-brighter behavior in the X-ray spectra and the measured positive X-ray/VHE correlation during the 2009 multi-wavelength campaign suggests that the physical processes dominating the emission during non-flaring states have similarities with those occurring during flaring activity. In particular, this observation supports leptonic scenarios as being responsible for the emission of Mrk421 during non-flaring activity. Such a temporally extended X-ray/VHE correlation is not driven by any single flaring event, and hence is difficult to explain within the standard hadronic scenarios. The highest variability is observed in the X-ray band, which, within the one-zone synchrotron self-Compton scenario, indicates that the electron energy distribution is most variable at the highest energies.

Cadmium Zinc Telluride and Cadmium Telluride are the detector materials of choice for the detection of X-rays in the X-ray energy band E >= 5keV with excellent spatial and spectral resolution and without cryogenic cooling. Owing to recent breakthroughs in grazing incidence mirror technology, next-generation hard X-ray telescopes will achieve angular resolution between 5 and 10 arc seconds - about an order of magnitude better than that of the NuSTAR hard X-ray telescope. As a consequence, the next generation of X-ray telescopes will require pixelated X-ray detectors with pixels on a grid with a lattice constant of <= 250um. Additional detector requirements include a low energy threshold of less than 5keV and an energy resolution of less than one keV. The science drivers for a high angular-resolution X-ray mission include studies and measurements of black hole spins, the cosmic evolution of super-massive black holes, active galactic nuclei feedback, and the behaviour of matter at very high densities. In this contribution, we report on our R&D studies with the goal to optimise small-pixel Cadmium Zinc Telluride and Cadmium Telluride detectors.

The nearby radio galaxy M87 offers a unique opportunity for exploring the connection between gamma-ray production and jet formation at an unprecedented linear resolution. However, the origin and location of the gamma-rays in this source is still elusive. Based on previous radio/TeV correlation events, the unresolved jet base (radio core) and the peculiar knot HST-1 at >120 pc from the nucleus are proposed as candidate site(s) of gamma-ray production. Here we report our intensive, high-resolution radio monitoring observations of the M87 jet with the VLBI Exploration of Radio Astrometry (VERA) and the European VLBI Network (EVN) from February 2011 to October 2012. During this period, an elevated level of the M87 flux is reported at TeV with VERITAS. We detected a remarkable flux increase in the radio core with VERA at 22/43 GHz coincident with the VHE activity. Meanwhile, HST-1 remained quiescent in terms of its flux density and structure at radio. These results strongly suggest that the TeV gamma-ray activity in 2012 originates in the jet base within 0.03 pc (projected) from the central supermassive black hole.

X-ray polarimetry promises to give qualitatively new information about high-energy astrophysical sources, such as binary black hole systems, micro-quasars, active galactic nuclei, neutron stars, and gamma-ray bursts. We designed, built and tested a X-ray polarimeter, X-Calibur, to be used in the focal plane of the balloon-borne InFOCuS grazing incidence X-ray telescope. X-Calibur combines a low-Z scatterer with a CZT detector assembly to measure the polarization of 20-80keV X-rays making use of the fact that polarized photons scatter preferentially perpendicular to the electric field orientation. X-Calibur achieves a high detection efficiency of ~80%. The X-Calibur detector assembly is completed, tested, and fully calibrated. The response to a polarized X-ray beam was measured successfully at the Cornell High Energy Synchrotron Source. This paper describes the design, calibration and performance of the X-Calibur polarimeter. In principle, a similar space-borne scattering polarimeter could operate over the broader 2-100keV energy band.

We summarize broadband observations of the TeV-emitting blazar 1ES 1959+650, including optical R-band observations by the robotic telescopes Super-LOTIS and iTelescope, UV observations by Swift UVOT, X-ray observations by the Swift X-ray Telescope (XRT), high-energy gamma-ray observations with the Fermi Large Area Telescope (LAT) and very-high-energy (VHE) gamma-ray observations by VERITAS above 315 GeV, all taken between 17 April 2012 and 1 June 2012 (MJD 56034 and 56079). The contemporaneous variability of the broadband spectral energy distribution is explored in the context of a simple synchrotron self Compton (SSC) model. In the SSC emission scenario, we find that the parameters required to represent the high state are significantly different than those in the low state. Motivated by possible evidence of gas in the vicinity of the blazar, we also investigate a reflected-emission model to describe the observed variability pattern. This model assumes that the non-thermal emission from the jet is reflected by a nearby cloud of gas, allowing the reflected emission to re-enter the blob and produce an elevated gamma-ray state with no simultaneous elevated synchrotron flux. The model applied here, although not required to explain the observed variability pattern, represents one possible scenario which can describe the observations. As applied to an elevated VHE state of 66% of the Crab Nebula flux, observed on a single night during the observation period, the reflected-emission scenario does not support a purely leptonic non-thermal emission mechanism. The reflected emission model does, however, predict a reflected photon field with sufficient energy to enable elevated gamma-ray emission via pion production with protons of energies between 10 and 100 TeV.

X-ray polarimetry has great scientific potential and new experiments, such as X-Calibur, PoGOLite, XIPE, and GEMS, will not only be orders of magnitude more sensitive than previous missions, but also provide the capability to measure polarization over a wide energy range. However, the measured spectra depend on the collection area, detector responses, and, in case of balloon-borne experiments, the absorption of X-rays in the atmosphere, all of which are energy dependent. Combined with the typically steep source spectra, this leads to significant biases that need to be taken into account to correctly reconstruct energy-resolved polarization properties. In this paper, we present a method based on an iterative unfolding algorithm that makes it possible to simultaneously reconstruct the energy spectrum and the polarization properties as a function of true photon energy. We apply the method to a simulated X-Calibur data set and show that it is able to recover both the energy spectrum and the energy-dependent polarization fraction.

We present results from VERITAS observations of the BL Lac object PG 1553+113 spanning the years 2010, 2011, and 2012. The time-averaged spectrum, measured between 160 and 560\,GeV, is well described by a power law with a spectral index of $4.33 \pm 0.09$. The time-averaged integral flux above $200\,$GeV measured for this period was $(1.69 \pm 0.06) \times 10^{-11} \, \mathrm{ph} \, \mathrm{cm}^{-2} \, \mathrm{s}^{-1}$, corresponding to 6.9\% of the Crab Nebula flux. We also present the combined $\gamma$-ray spectrum from the Fermi Large Area Telescope and VERITAS covering an energy range from 100~MeV to 560~GeV. The data are well fit by a power law with an exponential cutoff at $\rm {101.9 \pm 3.2 \, \mathrm{GeV}} $. The origin of the cutoff could be intrinsic to PG~1553+113 or be due to the $\gamma$-ray opacity of our universe through pair production off the extragalactic background light (EBL). Given lower limits to the redshift of $\rm z \negthinspace > \negthinspace 0.395$ based on optical/UV observations of PG~1553+113, the cutoff would be dominated by EBL absorption. Conversely, the small statistical uncertainties of the VERITAS energy spectrum have allowed us to provide a robust upper limit on the redshift of PG 1553+113 of $z \negthinspace \leq \negthinspace 0.62$. A strongly-elevated mean flux of $(2.50 \pm 0.14) \times 10^{-11} \, \mathrm{ph} \, \mathrm{cm}^{-2} \, \mathrm{s}^{-1}$ (10.3\% of the Crab Nebula flux) was observed during 2012, with the daily flux reaching as high as $(4.44 \pm 0.71) \times 10^{-11} \, \mathrm{ph} \, \mathrm{cm}^{-2} \, \mathrm{s}^{-1}$ (18.3\% of the Crab Nebula flux) on MJD 56048. The light curve measured during the 2012 observing season is marginally inconsistent with a steady flux, giving a $\chi^2$ probability for a steady flux of 0.03\%.

Mrk 501 is one of the brightest blazars at TeV energies and has been extensively studied since its first VHE detection in 1996. Our goal is to characterize in detail the source gamma-ray emission, together with the radio-to-X-ray emission, during the non-flaring (low) activity, which is less often studied than the occasional flaring (high) activity. We organized a multiwavelength (MW) campaign on Mrk 501 between March and May 2008. This multi-instrument effort included the most sensitive VHE gamma-ray instruments in the northern hemisphere, namely the imaging atmospheric Cherenkov telescopes MAGIC and VERITAS, as well as Swift, RXTE, the F-GAMMA, GASP-WEBT, and other collaborations and instruments. Mrk 501 was found to be in a low state of activity during the campaign, with a VHE flux in the range of 10%-20% of the Crab nebula flux. Nevertheless, significant flux variations were detected with various instruments, with a trend of increasing variability with energy. The broadband spectral energy distribution during the two different emission states of the campaign can be adequately described within the homogeneous one-zone synchrotron self-Compton model, with the (slightly) higher state described by an increase in the electron number density. This agrees with previous studies of the broadband emission of this source during flaring and non-flaring states. We report for the first time a tentative X-ray-to-VHE correlation during a low VHE activity. Although marginally significant, this positive correlation between X-ray and VHE, which has been reported many times during flaring activity, suggests that the mechanisms that dominate the X-ray/VHE emission during non-flaring-activity are not substantially different from those that are responsible for the emission during flaring activity.

Prompt emission from the very fluent and nearby (z=0.34) gamma-ray burst GRB 130427A was detected by several orbiting telescopes and by ground-based, wide-field-of-view optical transient monitors. Apart from the intensity and proximity of this GRB, it is exceptional due to the extremely long-lived high-energy (100 MeV to 100 GeV) gamma-ray emission, which was detected by the Large Area Telescope on the Fermi Gamma-ray Space Telescope for ~70 ks after the initial burst. The persistent, hard-spectrum, high-energy emission suggests that the highest-energy gamma rays may have been produced via synchrotron self-Compton processes though there is also evidence that the high-energy emission may instead be an extension of the synchrotron spectrum. VERITAS, a ground-based imaging atmospheric Cherenkov telescope array, began follow-up observations of GRB 130427A ~71 ks (~20 hr) after the onset of the burst. The GRB was not detected with VERITAS; however, the high elevation of the observations, coupled with the low redshift of the GRB, make VERITAS a very sensitive probe of the emission from GRB 130427A for E > 100 GeV. The non-detection and consequent upper limit derived place constraints on the synchrotron self-Compton model of high-energy gamma-ray emission from this burst.

We report our intensive radio monitoring observations of the jet in M87 with the VLBI Exploration of Radio Astrometry (VERA) and the European VLBI Network (EVN) from February 2011 to October 2012, together with contemporaneous high-energy gamma-ray light curves obtained by the Fermi-LAT. During this period, an elevated level of the M87 flux is reported at VHE gamma rays. We detected a remarkable increase of the radio flux density from the unresolved jet base (radio core) with VERA at 22 and 43GHz coincident with the VHE activity. Meanwhile, we confirmed with EVN at 5GHz that HST-1 (an alternative gamma-ray production candidate site) remained quiescent in terms of its flux density and structure. These results in the radio bands strongly suggest that the VHE gamma-ray activity in 2012 originates in the jet base within 0.03pc or 56 Schwarzschild radii from the central supermassive black hole. We further conducted VERA astrometry for the M87 core during the flaring period, and detected core shifts between 22 and 43GHz. We also discovered a clear frequency-dependent evolution of the radio core flare at 43, 22 and 5GHz; the radio flux density increased more rapidly at higher frequencies with a larger amplitude, and the light curves clearly showed a time-lag between the peaks at 22 and 43GHz. This indicates that a new radio-emitting component was created near the black hole in the period of the VHE event, and then propagated outward with progressively decreasing synchrotron opacity. By combining these results, we estimated an apparent speed of the newborn component, and derived a sub-luminal speed of less than ~0.2c. This value is significantly slower than the super-luminal (~1.1c) features that appeared from the core during the prominent VHE flaring event in 2008, suggesting that the stronger VHE activity can be associated with the production of the higher Lorentz factor jet.

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 $\Gamma=$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 $\tau=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.

TeV J2032+4130 was the first unidentified source discovered at very high energies (VHE; E $>$ 100 GeV), with no obvious counterpart in any other wavelength. It is also the first extended source to be observed in VHE gamma rays. Following its discovery, intensive observational campaigns have been carried out in all wavelengths in order to understand the nature of the object, which have met with limited success. We report here on a deep observation of TeV J2032+4130, based on 48.2 hours of data taken from 2009 to 2012 by the VERITAS (Very Energetic Radiation Imaging Telescope Array System) experiment. The source is detected at 8.7 standard deviations ($\sigma$) and is found to be extended and asymmetric with a width of 9.5$^{\prime}$$\pm$1.2$^{\prime}$ along the major axis and 4.0$^{\prime}$$\pm$0.5$^{\prime}$ along the minor axis. The spectrum is well described by a differential power law with an index of 2.10 $\pm$ 0.14$_{stat}$ $\pm$ 0.21$_{sys}$ and a normalization of (9.5 $\pm$ 1.6$_{stat}$ $\pm$ 2.2$_{sys}$) $\times$ 10$^{-13}$TeV$^{-1}$ cm$^{-2}$ s$^{-1}$ at 1 TeV. We interpret these results in the context of multiwavelength scenarios which particularly favor the pulsar wind nebula (PWN) interpretation.

The high-frequency-peaked BL Lacertae object 1ES 0229+200 is a relatively distant (z = 0.1396), hard-spectrum (Gamma ~ 2.5), very-high-energy-emitting (E > 100 GeV) gamma-ray blazar. Very-high-energy measurements of this active galactic nucleus have been used to place constraints on the intensity of the extragalactic background light and the intergalactic magnetic field. A multi-wavelength study of this object centered around very-high-energy observations by VERITAS is presented. This study obtained, over a period of three years, an 11.7 standard deviation detection and an average integral flux F(E>300 GeV) = (23.3 +- 2.8_stat +- 5.8_sys) x 10^-9 photons m^-2 s^-1, or 1.7% of the Crab Nebula's flux (assuming the Crab Nebula spectrum measured by H.E.S.S). Supporting observations from Swift and RXTE are analyzed. The Swift observations are combined with previously published Fermi observations and the very-high-energy measurements to produce an overall spectral energy distribution which is then modeled assuming one-zone synchrotron-self-Compton emission. The chi^2 probability of the TeV flux being constant is 1.6%. This, when considered in combination with measured variability in the X-ray band, and the demonstrated variability of many TeV blazars, suggests that the use of blazars such as 1ES 0229+200 for intergalactic magnetic field studies may not be straightforward and challenges models that attribute hard TeV spectra to secondary gamma-ray production along the line of sight.

HESS J0632+057 is the only gamma-ray binary known so far whose position in the sky allows observations with ground-based observatories both in the northern and southern hemispheres. Here we report on long-term observations of HESS J0632+057 conducted with the VERITAS and H.E.S.S. Cherenkov Telescopes and the X-ray Satellite Swift, spanning a time range from 2004 to 2012 and covering most of the system's orbit. The VHE emission is found to be variable, and is correlated with that at X-ray energies. An orbital period of $315 ^{+6}_{-4}$ days is derived from the X-ray data set, which is compatible with previous results, $P = (321 \pm 5$) days. The VHE light curve shows a distinct maximum at orbital phases close to 0.3, or about 100 days after periastron passage, which coincides with the periodic enhancement of the X-ray emission. Furthermore, the analysis of the TeV data shows for the first time a statistically significant ($> 6.5 \sigma$) detection at orbital phases 0.6--0.9. The obtained gamma-ray and X-ray light curves and the correlation of the source emission at these two energy bands are discussed in the context of the recent ephemeris obtained for the system. Our results are compared to those reported for other gamma-ray binaries.

We report results from TeV gamma-ray observations of the microquasar Cygnus X-3. The observations were made with the Very Energetic Radiation Imaging Telescope Array System (VERITAS) over a time period from 2007 June 11 to 2011 November 28. VERITAS is most sensitive to gamma rays at energies between 85 GeV to 30 TeV. The effective exposure time amounts to a total of about 44 hours, with the observations covering six distinct radio/X-ray states of the object. No significant TeV gamma-ray emission was detected in any of the states, nor with all observations combined. The lack of a positive signal, especially in the states where GeV gamma rays were detected, places constraints on TeV gamma-ray production in Cygnus X-3. We discuss the implications of the results.

The variability of the blazar Markarian 421 in TeV gamma rays over a 14-year time period has been explored with the Whipple 10 m telescope. It is shown that the dynamic range of its flux variations is large and similar to that in X-rays. A correlation between the X-ray and TeV energy bands is observed during some bright flares and when the complete data sets are binned on long timescales. The main database consists of 878.4 hours of observation with the Whipple telescope, spread over 783 nights. The peak energy response of the telescope was 400 GeV with 20% uncertainty. This is the largest database of any TeV-emitting active galactic nucleus (AGN) and hence was used to explore the variability profile of Markarian 421. The time-averaged flux from Markarian 421 over this period was 0.446$\pm$0.008 Crab flux units. The flux exceeded 10 Crab flux units on three separate occasions. For the 2000-2001 season the average flux reached 1.86 Crab units, while in the 1996-1997 season the average flux was only 0.23 Crab units.

Report of the CF6 Working Group at Snowmass 2013. Topics addressed include ultra-high energy cosmic rays, neutrinos, gamma rays, baryogenesis, and experiments probing the fundamental nature of spacetime.

We report on VERITAS observations of the BL Lac object B2 1215+30 between 2008 and 2012. During this period, the source was detected at very high energies (VHE; E > 100 GeV) by VERITAS with a significance of $8.9\sigma$ and showed clear variability on time scales larger than months. In 2011, the source was found to be in a relatively bright state and a power-law fit to the differential photon spectrum yields a spectral index of $3.6 \pm 0.4_{\mathrm{stat}} \pm 0.3_{\mathrm{syst}}$ with an integral flux above 200 GeV of $(8.0 \pm 0.9_{\mathrm{stat}} \pm 3.2_{\mathrm{syst}}) \times 10^{-12}\, \mathrm{cm}^{-2} \mathrm{s}^{-1}$. No short term variability could be detected during the bright state in 2011. Multi-wavelength data were obtained contemporaneous with the VERITAS observations in 2011 and cover optical (Super-LOTIS, MDM, Swift-UVOT), X-ray (Swift-XRT), and gamma-ray (Fermi-LAT) frequencies. These were used to construct the spectral energy distribution (SED) of B2 1215+30. A one-zone leptonic model is used to model the blazar emission and the results are compared to those of MAGIC from early 2011 and other VERITAS-detected blazars. The SED can be well reproduced with model parameters typical for VHE-detected BL Lacs.

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.

Compilation of papers contributed by the VERITAS Collaboration to the 33rd International Cosmic Ray Conference, held 2-9 July, 2013, in Rio de Janeiro, Brazil.

We report the detection of a new TeV gamma-ray source, VER J0521+211, based on observations made with the VERITAS imaging atmospheric Cherenkov telescope array. These observations were motivated by the discovery of a cluster of >30GeV photons in the first year of Fermi-LAT observations. VER J0521+211 is relatively bright at TeV energies, with a mean photon flux of 1.93 +/- 0.13_stat +/- 0.78_sys 10^-11 cm-2 s-1 above 0.2 TeV during the period of the VERITAS observations. The source is strongly variable on a daily timescale across all wavebands, from optical to TeV, with a peak flux corresponding to ~0.3 times the steady Crab Nebula flux at TeV energies. Follow-up observations in the optical and X-ray bands classify the newly-discovered TeV source as a BL Lac-type blazar with uncertain redshift, although recent measurements suggest z=0.108. VER J0521+211 exhibits all the defining properties of blazars in radio, optical, X-ray, and gamma-ray wavelengths.

X-ray and gamma-ray observations of astrophysical objects at cosmological distances can be used to probe the energy dependence of the speed of light with high accuracy and to search for violations of Lorentz invariance and CPT symmetry at the Planck energy scale. In this conference contribution, we discuss these searches in the theoretical framework of the Standard-Model Extension. We present new limits on the dispersion relation governed by operators of mass dimension d=5 and d=6, and we discuss avenues for future progress.

We report on the VERITAS observations of the high-frequency peaked BL Lac object 1ES 1959+650 in the period 2007-2011. This source is detected at TeV energies by VERITAS at 16.4 standard deviation (sigma) significance in 7.6 hours of observation in a low flux state. A multiwavelength spectral energy distribution (SED) is constructed from contemporaneous data from VERITAS, Fermi-LAT, RXTE PCA, and Swift UVOT. Swift XRT data is not included in the SED due to a lack of simultaneous observations with VERITAS. In contrast to the orphan gamma-ray flare exhibited by this source in 2002, the X-ray flux of the source is found to vary by an order of magnitude, while other energy regimes exhibit less variable emission. A quasi-equilibrium synchrotron self-Compton model with an additional external radiation field is used to describe three SEDs corresponding to the lowest, highest, and average X-ray states. The variation in the X-ray spectrum is modeled by changing the electron injection spectral index, with minor adjustments of the kinetic luminosity in electrons. This scenario produces small-scale flux variability of order >~2 in the HE (E>1 MeV) and VHE (E>100 GeV) gamma-ray regimes, which is corroborated by the Fermi-LAT, VERITAS, and Whipple 10m telescope light curves.

We report the discovery of an unidentified, extended source of very-high-energy (VHE) gamma-ray emission, VER J2019+407, within the radio shell of the supernova remnant SNR G78.2+2.1, using 21.4 hours of data taken by the VERITAS gamma-ray observatory in 2009. These data confirm the preliminary indications of gamma-ray emission previously seen in a two-year (2007-2009) blind survey of the Cygnus region by VERITAS. VER J2019+407, which is detected at a post-trials significance of 7.5 standard deviations in the 2009 data, is localized to the northwestern rim of the remnant in a region of enhanced radio and X-ray emission. It has an intrinsic extent of 0.23^{\circ} \pm 0.03^{\circ} (stat)+0.04^{\circ}_{-0.02}^{\circ}(sys) and its spectrum is well-characterized by a differential power law (dN/dE = N_0 \times (E/TeV)^{-\Gamma}) with a photon index of {\Gamma} = 2.37 \pm 0.14 (stat) \pm 0.20 (sys) and a flux normalization of N0 = 1.5 \pm 0.2 (stat) \pm 0.4(sys) \times 10^-12 ph TeV^{-1} cm^{-2} s^{-1}. This yields an integral flux of 5.2 \pm 0.8 (stat) \pm 1.4 (sys) \times 10^-12 ph cm^{-2} s^{-1} above 320 GeV, corresponding to 3.7% of the Crab Nebula flux. We consider the relationship of the TeV gamma-ray emission with the GeV gamma-ray emission seen from SNR G78.2+2.1 as well as that seen from a nearby cocoon of freshly accelerated cosmic rays. Multiple scenarios are considered as possible origins for the TeV gamma-ray emission, including hadronic particle acceleration at the supernova remnant shock.

We discuss the prospects of doing tests of Lorentz invariance with gamma-rays observed with present and future ground based gamma-ray observatories.

Active Galactic Nuclei (hereafter AGN) produce powerful outflows which offer excellent conditions for efficient particle acceleration in internal and external shocks, turbulence, and magnetic reconnection events. The jets as well as particle accelerating regions close to the supermassive black holes (hereafter SMBH) at the intersection of plasma inflows and outflows, can produce readily detectable very high energy gamma-ray emission. As of now, more than 45 AGN including 41 blazars and 4 radiogalaxies have been detected by the present ground-based gamma-ray telescopes, which represents more than one third of the cosmic sources detected so far in the VHE gamma-ray regime. The future Cherenkov Telescope Array (CTA) should boost the sample of AGN detected in the VHE range by about one order of magnitude, shedding new light on AGN population studies, and AGN classification and unification schemes. CTA will be a unique tool to scrutinize the extreme high-energy tail of accelerated particles in SMBH environments, to revisit the central engines and their associated relativistic jets, and to study the particle acceleration and emission mechanisms, particularly exploring the missing link between accretion physics, SMBH magnetospheres and jet formation. Monitoring of distant AGN will be an extremely rewarding observing program which will inform us about the inner workings and evolution of AGN. Furthermore these AGN are bright beacons of gamma-rays which will allow us to constrain the extragalactic infrared and optical backgrounds as well as the intergalactic magnetic field, and will enable tests of quantum gravity and other "exotic" phenomena.

In this document, we describe the scientific potential of blazar observations with a X-ray polarimetry mission like GEMS (Gravity and Extreme Magnetism SMEX). We describe five blazar science investigations that such a mission would enable: (i) the structure and the role of magnetic fields in AGN jets, (ii) analysis of the polarization of the synchrotron X-ray emission from AGN jets, (iii) discrimination between synchrotron self-Compton and external Compton models for blazars with inverse Compton emission in the X-ray band, (iv) a precision study of the polarization properties of the X-ray emission from Cen-A, (v) tests of Lorentz Invariance based on X-ray polarimetric observations of blazars. We conclude with a discussion of a straw man observation program and recommended accompanying multiwavelength observations.

We examine the expected X-ray polarization properties of neutron-star X-ray sources of various types, e.g., accretion and rotation powered pulsars, magnetars, and low-mass X-ray binaries. We summarize the model calculations leading to these expected properties. We describe how a comparison of these with their observed properties, as inferred from GEMS data, will probe the essential dynamical, electromagnetic, plasma, and emission processes in neutron-star binaries, discriminate between models of these processes, and constrain model parameters. An exciting goal is the first observational demonstration in this context of the existence of vacuum resonance, a fundamental quantum electrodynamical phenomenon first described in the 1930s.

In this paper we give a brief review of the astrophysics of active galactic nuclei (AGN). After a general introduction motivating the study of AGNs, we discuss our present understanding of the inner workings of the central engines, most likely accreting black holes with masses between a million and ten billion solar masses. We highlight recent results concerning the jets (collimated outflows) of AGNs derived from X-ray observations (Chandra) of kpc-scale jets and gamma-ray observations of AGNs (Fermi, Cherenkov telescopes) with jets closely aligned with the lines of sight (blazars), and discuss the interpretation of these observations. Subsequently, we summarize our knowledge about the cosmic history of AGN formation and evolution. We conclude with a description of upcoming observational opportunities.

We present here a summary of the scientific goals behind the Gravity and Extreme Magnetism SMEX (GEMS) X-ray polarimetry mission's black hole (BH) observing program. The primary targets can be divided into two classes: stellar-mass galactic BHs in accreting binaries, and super-massive BHs in the centers of active galactic nuclei (AGN). The stellar-mass BHs can in turn be divided into various X-ray spectral states: thermal-dominant (disk), hard (radio jet), and steep power-law (hot corona). These different spectral states are thought to be generated by different accretion geometries and emission mechanisms. X-ray polarization is an ideal tool for probing the geometry around these BHs and revealing the specific properties of the accreting gas.

We report the discovery of TeV gamma-ray emission coincident with the shell-type radio supernova remnant (SNR) CTA 1 using the VERITAS gamma-ray observatory. The source, VER J0006+729, was detected as a 6.5 standard deviation excess over background and shows an extended morphology, approximated by a two-dimensional Gaussian of semi-major (semi-minor) axis 0.30 degree (0.24 degree) and a centroid 5' from the Fermi gamma-ray pulsar PSR J0007+7303 and its X-ray pulsar wind nebula (PWN). The photon spectrum is well described by a power-law dN/dE = N_0 (E/3 TeV)^(-\Gamma), with a differential spectral index of \Gamma = 2.2 +- 0.2_stat +- 0.3_sys, and normalization N_0 = (9.1 +- 1.3_stat +- 1.7_sys) x 10^(-14) cm^(-2) s^(-1) TeV^(-1). The integral flux, F_\gamma = 4.0 x 10^(-12) erg cm^(-2) s^(-1) above 1 TeV, corresponds to 0.2% of the pulsar spin-down power at 1.4 kpc. The energetics, co-location with the SNR, and the relatively small extent of the TeV emission strongly argue for the PWN origin of the TeV photons. We consider the origin of the TeV emission in CTA 1.

We report on the optimization of the hard X-ray polarimeter X-Calibur for a high-altitude balloon-flight in the focal plane of the InFOC{\mu}S X-ray telescope from Fort Sumner (NM) in Fall 2013. X-Calibur combines a low-Z scintillator slab to Compton-scatter photons with a high-Z Cadmium Zinc Telluride (CZT) detector assembly to photo-absorb the scattered photons. The detector makes use of the fact that polarized photons Compton scatter preferentially perpendicular to the electric field orientation. X-Calibur achieves a high detection efficiency of order unity and reaches a sensitivity close to the best theoretically possible. In this paper, we discuss the optimization of the design of the instrument based on Monte Carlo simulations of polarized and unpolarized X-ray beams and of the most important background components. We calculate the sensitivity of the polarimeter for the upcoming balloon flight from Fort Sumner and for additional longer balloon flights with higher throughput mirrors. We conclude by emphasizing that Compton polarimeters on satellite borne missions can be used down to energies of a few keV.

We report on the detection of a very rapid TeV gamma-ray flare from BL Lacertae on 2011 June 28 with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). The flaring activity was observed during a 34.6-minute exposure, when the integral flux above 200 GeV reached $(3.4\pm0.6) \times 10^{-6} \;\text{photons}\;\text{m}^{-2}\text{s}^{-1}$, roughly 125% of the Crab Nebula flux measured by VERITAS. The light curve indicates that the observations missed the rising phase of the flare but covered a significant portion of the decaying phase. The exponential decay time was determined to be $13\pm4$ minutes, making it one of the most rapid gamma-ray flares seen from a TeV blazar. The gamma-ray spectrum of BL Lacertae during the flare was soft, with a photon index of $3.6\pm 0.4$, which is in agreement with the measurement made previously by MAGIC in a lower flaring state. Contemporaneous radio observations of the source with the Very Long Baseline Array (VLBA) revealed the emergence of a new, superluminal component from the core around the time of the TeV gamma-ray flare, accompanied by changes in the optical polarization angle. Changes in flux also appear to have occurred at optical, UV, and GeV gamma-ray wavelengths at the time of the flare, although they are difficult to quantify precisely due to sparse coverage. A strong flare was seen at radio wavelengths roughly four months later, which might be related to the gamma-ray flaring activities. We discuss the implications of these multiwavelength results.

We report on VERITAS very-high-energy (VHE; E>100 GeV) observations of six blazars selected from the Fermi Large Area Telescope First Source Catalog (1FGL). The gamma-ray emission from 1FGL sources was extrapolated up to the VHE band, taking gamma-ray absorption by the extragalactic background light into account. This allowed the selection of six bright, hard-spectrum blazars that were good candidate TeV emitters. Spectroscopic redshift measurements were attempted with the Keck Telescope for the targets without Sloan Digital Sky Survey (SDSS) spectroscopic data. No VHE emission is detected during the observations of the six sources described here. Corresponding TeV upper limits are presented, along with contemporaneous Fermi observations and non-concurrent Swift UVOT and XRT data. The blazar broadband spectral energy distributions (SEDs) are assembled and modeled with a single-zone synchrotron self-Compton model. The SED built for each of the six blazars show a synchrotron peak bordering between the intermediate- and high-spectrum-peak classifications, with four of the six resulting in particle-dominated emission regions.

We present the results of a joint observational campaign between the Green Bank radio telescope and the VERITAS gamma-ray telescope, which searched for a correlation between the emission of very-high-energy (VHE) gamma rays ($E_{\gamma} >$ 150 GeV) and Giant Radio Pulses (GRPs) from the Crab pulsar at 8.9 GHz. A total of 15366 GRPs were recorded during 11.6 hours of simultaneous observations, which were made across four nights in December 2008 and in November and December 2009. We searched for an enhancement of the pulsed gamma-ray emission within time windows placed around the arrival time of the GRP events. In total, 8 different time windows with durations ranging from 0.033 ms to 72 s were positioned at three different locations relative to the GRP to search for enhanced gamma-ray emission which lagged, led, or was concurrent with, the GRP event. Further, we performed separate searches on main pulse GRPs and interpulse GRPs and on the most energetic GRPs in our data sample. No significant enhancement of pulsed VHE emission was found in any of the preformed searches. We set upper limits of 5-10 times the average VHE flux of the Crab pulsar on the flux simultaneous with interpulse GRPs on single-rotation-period time scales. On $\sim$8-second time scales around interpulse GRPs, we set an upper limit of 2-3 times the average VHE flux. Within the framework of recent models for pulsed VHE emission from the Crab pulsar, the expected VHE-GRP emission correlations are below the derived limits.

Observations of radio halos and relics in galaxy clusters indicate efficient electron acceleration. Protons should likewise be accelerated, suggesting that clusters may also be sources of very high-energy (VHE; E>100 GeV) gamma-ray emission. We report here on VHE gamma-ray observations of the Coma galaxy cluster with the VERITAS array of imaging Cherenkov telescopes, with complementing Fermi-LAT observations at GeV energies. No significant gamma-ray emission from the Coma cluster was detected. Integral flux upper limits at the 99% confidence level were measured to be on the order of (2-5)*10^-8\ ph. m^-2 s^-1 (VERITAS, >220 GeV} and ~2*10^-6 ph. m^-2 s^-1 (Fermi, 1-3 GeV), respectively. We use the gamma-ray upper limits to constrain CRs and magnetic fields in Coma. Using an analytical approach, the CR-to-thermal pressure ratio is constrained to be < 16% from VERITAS data and < 1.7% from Fermi data (averaged within the virial radius). These upper limits are starting to constrain the CR physics in self-consistent cosmological cluster simulations and cap the maximum CR acceleration efficiency at structure formation shocks to be <50%. Assuming that the radio-emitting electrons of the Coma halo result from hadronic CR interactions, the observations imply a lower limit on the central magnetic field in Coma of (2 - 5.5) muG, depending on the radial magnetic-field profile and on the gamma-ray spectral index. Since these values are below those inferred by Faraday rotation measurements in Coma (for most of the parameter space), this {renders} the hadronic model a very plausible explanation of the Coma radio halo. Finally, since galaxy clusters are dark-matter (DM) dominated, the VERITAS upper limits have been used to place constraints on the thermally-averaged product of the total self-annihilation cross section and the relative velocity of the DM particles, <\sigma v>. (abr.)

We present observations of the BL Lac object 1ES 0414+009 in the >200 GeV gamma-ray band by the VERITAS array of Cherenkov telescopes. 1ES 0414+009 was observed by VERITAS between January 2008 and February 2011, resulting in 56.2 hours of good quality pointed observations. These observations resulted in a detection of 822 events from the source corresponding to a statistical significance of 6.4 standard deviations (6.4 sigma) above the background. The source flux, showing no evidence for variability, is measured as 5.2 +/- 1.1_stat +/- 2.6_sys * 10^-12 photons cm^-2 s^-1 above 200 GeV, equivalent to approximately 2% of the Crab Nebula flux above this energy. The differential photon spectrum from 230 GeV to 850 GeV is well fit by a power law with an photon index of Gamma 3.4 +/- 0.5_stat +/- 0.3_sys and a flux normalization of 1.6 +/- 0.3_stat +/- 0.8_sys * 10^-11 photons cm^-2 s^-1 at 300 GeV. We also present multiwavelength results taken in the optical (MDM), X-ray (Swift-XRT), and GeV (Fermi-LAT) bands and use these results to construct a broadband spectral energy distribution (SED). Modeling of this SED indicates that homogenous one-zone leptonic scenarios are not adequate to describe emission from the system, with a lepto-hadronic model providing a better fit to the data.

Although General Relativity (GR) has been tested extensively in the weak gravity regime, similar tests in the strong gravity regime are still missing. In this paper we explore the possibility to use X-ray spectropolarimetric observations of black holes in X-ray binaries to distinguish between the Kerr metric and the phenomenological metrics introduced by Johannsen and Psaltis (2011) (which are not vacuum solutions of Einstein's equation) and thus to test the no-hair theorem of GR. To this end, we have developed a numerical code that calculates the radial brightness profiles of accretion disks and parallel transports the wave vector and polarization vector of photons through the Kerr and non-GR spacetimes. We used the code to predict the observational appearance of GR and non-GR accreting black hole systems. We find that the predicted energy spectra and energy dependent polarization degree and polarization direction do depend strongly on the underlying spacetime. However, for large regions of the parameter space, the GR and non-GR metrics lead to very similar observational signatures, making it difficult to observationally distinguish between the two types of models.

We report on very high energy (E > 100 GeV) gamma-ray observations of V407 Cygni, a symbiotic binary that underwent a nova outburst producing 0.1-10 GeV gamma rays during 2010 March 10-26. Observations were made with the Very Energetic Radiation Imaging Telescope Array System during 2010 March 19-26 at relatively large zenith angles, due to the position of V407 Cyg. An improved reconstruction technique for large zenith angle observations is presented and used to analyze the data. We do not detect V407 Cygni and place a differential upper limit on the flux at 1.6 TeV of 2.3 \times 10^(-12) erg cm^(-2) s^(-1) (at the 95% confidence level). When considered jointly with data from Fermi-LAT, this result places limits on the acceleration of very high energy particles in the nova.

The current generation of Imaging Atmospheric Cherenkov Telescopes (IACTs), including the H.E.S.S., MAGIC, and VERITAS telescope arrays, have made substantial contributions to our knowledge about the structure and composition of the highly relativistic jets from Active Galactic Nuclei (AGNs). In this paper, we discuss some of the outstanding scientific questions and give a qualitative overview of AGN related science topics which will be explored with the next-generation Cherenkov Telescope Array (CTA). CTA is expected to further constrain the structure and make-up of jets, and thus, to constrain models of jet formation, acceleration, and collimation. Furthermore, being the brightest well-established extragalactic sources of TeV {\gamma}- rays, AGNs can be used to probe the EBL, intergalactic magnetic fields, and the validity of the Lorentz Invariance principle at high photon energies.

In this white paper, we discuss the concept of a next-generation X-ray mission called BEST (Black hole Evolution and Space Time). The mission concept uses a 3000 square centimeter effective area mirror (at 6 keV) to achieve unprecedented sensitivities for hard X-ray imaging spectrometry (5-70 keV) and for broadband X-ray polarimetry (2-70 keV). BEST can make substantial contributions to our understanding of the inner workings of accreting black holes, our knowledge about the fabric of extremely curved spacetime, and the evolution of supermassive black holes. BEST will allow for time resolved studies of accretion disks. With a more than seven times larger mirror area and a seven times wider bandpass than GEMS, BEST will take X-ray polarimetry to a new level: it will probe the time variability of the X-ray polarization from stellar mass and supermassive black holes, and it will measure the polarization properties in 30 independent energy bins. These capabilities will allow BEST to conduct tests of accretion disk models and the underlying spacetimes. With three times larger mirror area and ten times better angular resolution than NuSTAR, BEST will be able to make deep field observations with a more than 15 times better sensitivity than NuSTAR. The mission will be able to trace the evolution of obscured and unobscured black holes in the redshift range from zero to six, covering the most important epoch of supermassive black hole growth. The hard X-ray sensitivity of BEST will enable a deep census of non-thermal particle populations. BEST will give us insights into AGN feedback by measuring the particle luminosity injected by AGNs into the interstellar medium (ISM) of their hosts, and will map the emission from particles accelerated at large scale structure shocks. Finally, BEST has the potential to constrain the equation of state of neutron stars (NS).

In this review we confront the current theoretical understanding of particle acceleration at relativistic outflows with recent observational results on various source classes thought to involve such outflows, e.g. gamma-ray bursts, active galactic nuclei, and pulsar wind nebulae. We highlight the possible contributions of these sources to ultra-high-energy cosmic rays.

We report on the discovery of high-energy (HE; E > 0.1 GeV) and very high-energy (VHE; E > 100 GeV) gamma-ray emission from the high-frequency-peaked BL Lac object RBS 0413. VERITAS, a ground-based gamma-ray observatory, detected VHE gamma rays from RBS 0413 with a statistical significance of 5.5 standard deviations (sigma) and a gamma-ray flux of (1.5 \pm 0.6stat \pm 0.7syst) \times 10^(-8) photons m^(-2) s^(-1) (\sim 1% of the Crab Nebula flux) above 250 GeV. The observed spectrum can be described by a power law with a photon index of 3.18 \pm 0.68stat \pm 0.30syst. Contemporaneous observations with the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope detected HE gamma rays from RBS 0413 with a statistical significance of more than 9 sigma, a power-law photon index of 1.57 \pm 0.12stat +0.11sys -0.12sys and a gamma-ray flux between 300 MeV and 300 GeV of (1.64 \pm 0.43stat +0.31sys -0.22sys) \times 10^(-5) photons m^(-2) s^(-1). We present the results from Fermi-LAT and VERITAS, including a spectral energy distribution modeling of the gamma-ray, quasi-simultaneous X-ray (Swift-XRT), ultraviolet (Swift-UVOT) and R-band optical (MDM) data. We find that, if conditions close to equipartition are required, both the combined synchrotron self-Compton/external-Compton and the lepto-hadronic models are preferred over a pure synchrotron self-Compton model.

Abridged: The giant radio galaxy M 87 with its proximity, famous jet, and very massive black hole provides a unique opportunity to investigate the origin of very high energy (VHE; E>100 GeV) gamma-ray emission generated in relativistic outflows and the surroundings of super-massive black holes. M 87 has been established as a VHE gamma-ray emitter since 2006. The VHE gamma-ray emission displays strong variability on timescales as short as a day. In this paper, results from a joint VHE monitoring campaign on M 87 by the MAGIC and VERITAS instruments in 2010 are reported. During the campaign, a flare at VHE was detected triggering further observations at VHE (H.E.S.S.), X-rays (Chandra), and radio (43 GHz VLBA). The excellent sampling of the VHE gamma-ray light curve enables one to derive a precise temporal characterization of the flare: the single, isolated flare is well described by a two-sided exponential function with significantly different flux rise and decay times. While the overall variability pattern of the 2010 flare appears somewhat different from that of previous VHE flares in 2005 and 2008, they share very similar timescales (~day), peak fluxes (Phi(>0.35 TeV) ~= (1-3) x 10^-11 ph cm^-2 s^-1), and VHE spectra. 43 GHz VLBA radio observations of the inner jet regions indicate no enhanced flux in 2010 in contrast to observations in 2008, where an increase of the radio flux of the innermost core regions coincided with a VHE flare. On the other hand, Chandra X-ray observations taken ~3 days after the peak of the VHE gamma-ray emission reveal an enhanced flux from the core. The long-term (2001-2010) multi-wavelength light curve of M 87, spanning from radio to VHE and including data from HST, LT, VLA and EVN, is used to further investigate the origin of the VHE gamma-ray emission. No unique, common MWL signature of the three VHE flares has been identified.

Aims: We aim to constrain the structural variations within the HST-1 region downstream of the radio jet of M87, in general as well as in connection to the episodes of activity at very high energy (VHE). Methods: We analyzed and compared 26 VLBI observations of the M87 jet, obtained between 2006 and 2011 with the Very Long Baseline Array (VLBA) at 1.7 GHz and the European VLBI Network (EVN) at 5 GHz. Results: HST-1 is detected at all epochs; we model-fitted its complex structure with two or more components, the two outermost of which display a significant proper motion with a superluminal velocity around ~4c. The motion of a third feature that is detected upstream is more difficult to characterize. The overall position angle of HST-1 has changed during the time of our observations from -65deg to -90deg, while the structure has moved by over 80 mas downstream. Our results on the component evolution suggest that structural changes at the upstream edge of HST-1 can be related to the VHE events.

VERITAS has been monitoring the very-high-energy (VHE; >100GeV) gamma-ray activity of the radio galaxy M87 since 2007. During 2008, flaring activity on a timescale of a few days was observed with a peak flux of (0.70 +- 0.16) X 10^{-11} cm^{-2} s^{-1} at energies above 350GeV. In 2010 April, VERITAS detected a flare from M87 with peak flux of (2.71 +- 0.68) X 10^{-11} cm^{-2} s^{-1} for E>350GeV. The source was observed for six consecutive nights during the flare, resulting in a total of 21 hr of good quality data. The most rapid flux variation occurred on the trailing edge of the flare with an exponential flux decay time of 0.90^{+0.22}_{-0.15} days. The shortest detected exponential rise time is three times as long, at 2.87^{+1.65}_{-0.99} days. The quality of the data sample is such that spectral analysis can be performed for three periods: rising flux, peak flux, and falling flux. The spectra obtained are consistent with power-law forms. The spectral index at the peak of the flare is equal to 2.19 +- 0.07. There is some indication that the spectrum is softer in the falling phase of the flare than the peak phase, with a confidence level corresponding to 3.6 standard deviations. We discuss the implications of these results for the acceleration and cooling rates of VHE electrons in M87 and the constraints they provide on the physical size of the emitting region.

We present the results of sixteen Swift-triggered GRB follow-up observations taken with the VERITAS telescope array from January, 2007 to June, 2009. The median energy threshold and response time of these observations was 260 GeV and 320 s, respectively. Observations had an average duration of 90 minutes. Each burst is analyzed independently in two modes: over the whole duration of the observations and again over a shorter time scale determined by the maximum VERITAS sensitivity to a burst with a t^-1.5 time profile. This temporal model is characteristic of GRB afterglows with high-energy, long-lived emission that have been detected by the Large Area Telescope (LAT) on-board the Fermi satellite. No significant VHE gamma-ray emission was detected and upper limits above the VERITAS threshold energy are calculated. The VERITAS upper limits are corrected for gamma-ray extinction by the extragalactic background light (EBL) and interpreted in the context of the keV emission detected by Swift. For some bursts the VHE emission must have less power than the keV emission, placing constraints on inverse Compton models of VHE emission.

Compilation of papers contributed by the VERITAS Collaboration to the 32nd International Cosmic Ray Conference, held 11-18 August 2011 in Beijing, China.

We report on the VERITAS discovery of very-high-energy (VHE) gamma- ray emission above 200 GeV from the high-frequency-peaked BL Lac object RXJ0648.7+1516 (GBJ0648+1516), associated with 1FGLJ0648.8+1516. The photon spectrum above 200 GeV is fit by a power law dN/dE = F0(E/E0)-{\Gamma} with a photon index {\Gamma} of 4.4 {\pm} 0.8stat {\pm}0.3syst and a flux normalization F0 of (2.3 {\pm}0.5stat {\pm} 1.2sys) {\times} 10-11 TeV-1cm-2s-1 with E0 = 300 GeV. No VHE vari- ability is detected during VERITAS observations of RXJ0648.7+1516 between 2010 March 4 and April 15. Following the VHE discovery, the optical identifica- tion and spectroscopic redshift were obtained using the Shane 3-m Telescope at the Lick Observatory, showing the unidentified object to be a BL Lac type with a redshift of z = 0.179. Broadband multiwavelength observations contemporaneous with the VERITAS exposure period can be used to sub-classify the blazar as a high-frequency-peaked BL Lac (HBL) object, including data from the MDM ob- servatory, Swift-UVOT and XRT, and continuous monitoring at photon energies above 1 GeV from the Fermi Large Area Telescope (LAT). We find that in the absence of undetected, high-energy rapid variability, the one-zone synchrotron self-Compton model (SSC) overproduces the high-energy gamma-ray emission measured by the Fermi-LAT over 2.3 years. The SED can be parameterized sat- isfactorily with an external-Compton or lepto-hadronic model, which have two and six additional free parameters, respectively, compared to the one-zone SSC model.

The VERITAS telescope array has been operating smoothly since 2007, and has detected gamma-ray emission above 100 GeV from 40 astrophysical sources. These include blazars, pulsar wind nebulae, supernova remnants, gamma-ray binary systems, a starburst galaxy, a radio galaxy, the Crab pulsar, and gamma-ray sources whose origin remains unidentified. In 2009, the array was reconfigured, greatly improving the sensitivity. We summarize the current status of the observatory, describe some of the scientific highlights since 2009, and outline plans for the future.

X-ray polarimetry promises to give new information about high-energy astrophysical sources, such as binary black hole systems, micro-quasars, active galactic nuclei, and gamma-ray bursts. We designed, built and tested a hard X-ray polarimeter X-Calibur to be used in the focal plane of the InFOCuS grazing incidence hard X-ray telescope. X-Calibur combines a low-Z Compton scatterer with a CZT detector assembly to measure the polarization of 10-80 keV X-rays making use of the fact that polarized photons Compton scatter preferentially perpendicular to the electric field orientation. X-Calibur achieves a high detection efficiency of order unity.

X-ray polarimetry promises to give qualitatively new information about high-energy astrophysical sources, such as binary black hole systems, micro-quasars, active galactic nuclei, and gamma-ray bursts. We designed, built and tested a hard X-ray polarimeter X-Calibur to be used in the focal plane of the InFOCuS grazing incidence hard X-ray telescope. X-Calibur combines a low-Z Compton scatterer with a CZT detector assembly to measure the polarization of 10-80 keV X-rays making use of the fact that polarized photons Compton scatter preferentially perpendicular to the electric field orientation. X-Calibur achieves a high detection efficiency of order unity.

NASA's Small Explorer Mission GEMS (Gravity and Extreme Magnetism SMEX), scheduled for launch in 2014, will have the sensitivity to detect and measure the linear polarization properties of the 0.5 keV and 2-10 keV X-ray emission of a considerable number of galactic and extragalactic sources. The prospect of sensitive X-ray polarimetry justifies a closer look at the polarization properties of the basic emission mechanisms. In this paper, we present analytical and numerical calculations of the linear polarization properties of inverse Compton scattered radiation. We describe a generally applicable formalism that can be used to numerically compute the polarization properties in the Thomson and Klein-Nishina regimes. We use the code to perform for the first time a detailed comparison of numerical results and the earlier analytical results derived by Bonometto et al. (1970) for scatterings in the Thomson regime. Furthermore, we use the numerical formalism to scrutinize the polarization properties of synchrotron self-Compton emission, and of inverse Compton radiation emitted in the Klein-Nishina regime. We conclude with a discussion of the scientific potential of future GEMS observations of blazars. The GEMS mission will be able to confirm the synchrotron origin of the low-energy emission component from high frequency peaked BL Lac objects. Furthermore, the observations have the potential to decide between a synchrotron self-Compton and external-Compton origin of the high-energy emission component from flat spectrum radio quasars and low frequency peaked BL Lac objects.

We report the detection of pulsed gamma rays from the Crab pulsar at energies above 100 Gigaelectronvolts (GeV) with the VERITAS array of atmospheric Cherenkov telescopes. The detection cannot be explained on the basis of current pulsar models. The photon spectrum of pulsed emission between 100 Megaelectronvolts (MeV) and 400 GeV is described by a broken power law that is statistically preferred over a power law with an exponential cutoff. It is unlikely that the observation can be explained by invoking curvature radiation as the origin of the observed gamma rays above 100 GeV. Our findings require that these gamma rays be produced more than 10 stellar radii from the neutron star.

We successfully carried out the first high-altitude balloon flight of a wide-field hard X-ray coded-aperture telescope ProtoEXIST1, which was launched from the Columbia Scientific Balloon Facility at Ft. Sumner, New Mexico on October 9, 2009. ProtoEXIST1 is the first implementation of an advanced CdZnTe (CZT) imaging detector in our ongoing program to establish the technology required for next generation wide-field hard X-ray telescopes. The CZT detector plane in ProtoEXIST1 consists of an 8 x 8 array of closely tiled 2 cm x 2 cm x 0.5 cm thick pixellated CZT crystals, each with 8 x 8 pixels, covering a 256 cm^2 active area with 2.5 mm pixels. A tungsten mask, mounted at 90 cm above the detector provides shadowgrams of X-ray sources in the 30 - 600 keV band for imaging, allowing a fully coded field of view of 9 Deg x 9 Deg with an angular resolution of 20 arcmin. To reduce the background radiation, the detector is surrounded by semi-graded (Pb/Sn/Cu) passive shields on the four sides all the way to the mask. On the back side, a 26 cm x 26 cm x 2 cm CsI(Na) active shield provides signals to tag charged particle induced events as well as >~ 100 keV background photons from below. The flight duration was only about 7.5 hours due to strong winds (60 knots) at float altitude (38-39 km). Throughout the flight, the CZT detector performed excellently. The telescope observed Cyg X-1 for ~ 1 hour at the end of the flight. Despite a few problems with the pointing and aspect systems that caused the telescope to track about 6.4 deg off the target, the analysis of the Cyg X-1 data revealed an X-ray source at 7.2 sigma in the 30-100 keV energy band at the expected location from the optical images taken by the onboard daytime star camera. The success of this first flight is very encouraging for the future development of the advanced CZT imaging detectors (ProtoEXIST2, with 0.6 mm pixels).

Giant X-ray outbursts, with luminosities of about $ 10^{37}$ erg s$^{-1}$, are observed roughly every 5 years from the nearby Be/pulsar binary 1A 0535+262. In this article, we present observations of the source with VERITAS at very-high energies (VHE; E$>$100 GeV) triggered by the X-ray outburst in December 2009. The observations started shortly after the onset of the outburst, and they provided comprehensive coverage of the episode, as well as the 111-day binary orbit. No VHE emission is evident at any time. We also examined data from the contemporaneous observations of 1A 0535+262 with the Fermi/LAT at high energy photons (HE; E$>$0.1 GeV) and failed to detect the source at GeV energies. The X-ray continua measured with the Swift/XRT and the RXTE/PCA can be well described by the combination of blackbody and Comptonized emission from thermal electrons. Therefore, the gamma-ray and X-ray observations suggest the absence of a significant population of non-thermal particles in the system. This distinguishes 1A~0535+262 from those Be X-ray binaries (such as PSR B1259--63 and LS I +61$^{\circ}$303) that have been detected at GeV--TeV energies. We discuss the implications of the results on theoretical models.

We report on very-high-energy ($>$100 GeV) gamma-ray observations of Swift J164449.3+573451, an unusual transient object first detected by the {\it Swift} Observatory and later detected by multiple radio, optical and X-ray observatories. A total exposure of 28 hours was obtained on Swift J164449.3+573451 with VERITAS during 2011 March 28 -- April 15. We do not detect the source and place a differential upper limit on the emission at 500 GeV during these observations of $1.4 \times 10^{-12}$ erg cm$^{-2}$ s$^{-1}$ (99% confidence level). We also present time-resolved upper limits and use a flux limit averaged over the X-ray flaring period to constrain various emission scenarios that can accommodate both the radio-through-X-ray emission detected from the source and the lack of detection by VERITAS.

Multiwavelength observations of the high-frequency-peaked blazar 1ES2344+514 were performed from 2007 October to 2008 January. The campaign represents the first contemporaneous data on the object at very high energy (VHE, E >100 GeV) {\gamma}-ray, X-ray, and UV energies. Observations with VERITAS in VHE {\gamma}-rays yield a strong detection of 20 {\sigma} with 633 excess events in a total exposure of 18.1 hours live-time. A strong VHE {\gamma}-ray flare on 2007 December 7 is measured at F(>300 GeV) = (6.76 \pm 0.62) \times 10-11 ph cm-2 s-1, corresponding to 48% of the Crab Nebula flux. Excluding this flaring episode, nightly variability at lower fluxes is observed with a time-averaged mean of F(>300 GeV) = (1.06 \pm 0.09) \times 10-11 ph cm-2 s-1 (7.6% of the Crab Nebula flux). The differential photon spectrum between 390 GeV and 8.3 TeV for the time-averaged observations excluding 2007 December 7 is well described by a power law with a photon index of {\Gamma} = 2.78 \pm 0.09stat \pm 0.15syst. Over the full period of VERITAS observations contemporaneous X-ray and UV data were taken with Swift and RXTE. The measured 2-10 keV flux ranged by a factor of ~7 during the campaign. On 2007 December 8 the highest ever observed X-ray flux from 1ES 2344+514 was measured by Swift XRT at a flux of F(2-10 keV) = (6.28 \pm 0.31) \times 10-11 erg cm-2 s-1. Evidence for a correlation between the X-ray flux and VHE {\gamma}-ray flux on nightly time-scales is indicated with a Pearson correlation coefficient of r = 0.60 \pm 0.11. Contemporaneous spectral energy distributions (SEDs) of 1ES 2344+514 are presented for two distinct flux states. A one-zone synchrotron self-Compton (SSC) model describes both SEDs using parameters consistent with previous SSC modeling of 1ES 2344+514 from non-contemporaneous observations.

We report on TeV gamma-ray observations of the blazar Mrk 421 (redshift of 0.031) with the VERITAS observatory and the Whipple 10m Cherenkov telescope. The excellent sensitivity of VERITAS allowed us to sample the TeV gamma-ray fluxes and energy spectra with unprecedented accuracy where Mrk 421 was detected in each of the pointings. A total of 47.3 hrs of VERITAS and 96 hrs of Whipple 10m data were acquired between January 2006 and June 2008. We present the results of a study of the TeV gamma-ray energy spectra as a function of time, and for different flux levels. On May 2nd and 3rd, 2008, bright TeV gamma-ray flares were detected with fluxes reaching the level of 10 Crab. The TeV gamma-ray data were complemented with radio, optical, and X-ray observations, with flux variability found in all bands except for the radio waveband. The combination of the RXTE and Swift X-ray data reveal spectral hardening with increasing flux levels, often correlated with an increase of the source activity in TeV gamma-rays. Contemporaneous spectral energy distributions were generated for 18 nights, each of which are reasonably described by a one-zone SSC model.

We present the results of observations of the TeV binary LS I +61 303 with the VERITAS telescope array between 2008 and 2010, at energies above 300 GeV. In the past, both ground-based gamma-ray telescopes VERITAS and MAGIC have reported detections of TeV emission near the apastron phases of the binary orbit. The observations presented here show no strong evidence for TeV emission during these orbital phases; however, during observations taken in late 2010, significant emission was detected from the source close to the phase of superior conjunction (much closer to periastron passage) at a 5.6 standard deviation (5.6 sigma) post-trials significance. In total, between October 2008 and December 2010 a total exposure of 64.5 hours was accumulated with VERITAS on LS I +61 303, resulting in an excess at the 3.3 sigma significance level for constant emission over the entire integrated dataset. The flux upper limits derived for emission during the previously reliably active TeV phases (i.e. close to apastron) are less than 5% of the Crab Nebula flux in the same energy range. This result stands in apparent contrast to previous observations by both MAGIC and VERITAS which detected the source during these phases at >10% of the Crab Nebula flux. During the two year span of observations, a large amount of X-ray data were also accrued on LS I +61 303 by the Swift X-ray Telescope (XRT) and the Rossi X-ray Timing Explorer Timing (RXTE) Proportional Counter Array (PCA). We find no evidence for a correlation between emission in the X-ray and TeV regimes during 20 directly overlapping observations. We also comment on data obtained contemporaneously by the Fermi Large Area Telescope (LAT).

We report the discovery of TeV gamma-ray emission from the Type Ia supernova remnant (SNR) G120.1+1.4, known as Tycho's supernova remnant. Observations performed in the period 2008-2010 with the VERITAS ground-based gamma-ray observatory reveal weak emission coming from the direction of the remnant, compatible with a point source located at $00^{\rm h} \ 25^{\rm m} \ 27.0^{\rm s},\ +64^{\circ} \ 10^{\prime} \ 50^{\prime\prime}$ (J2000). The TeV photon spectrum measured by VERITAS can be described with a power-law $dN/dE = C(E/3.42\;\textrm{TeV})^{-\Gamma}$ with $\Gamma = 1.95 \pm 0.51_{stat} \pm 0.30_{sys}$ and $C = (1.55 \pm 0.43_{stat} \pm 0.47_{sys}) \times 10^{-14}$ cm$^{-2}$s$^{-1}$TeV$^{-1}$. The integral flux above 1 TeV corresponds to $\sim 0.9%$ percent of the steady Crab Nebula emission above the same energy, making it one of the weakest sources yet detected in TeV gamma rays. We present both leptonic and hadronic models which can describe the data. The lowest magnetic field allowed in these models is $\sim 80 \mu$G, which may be interpreted as evidence for magnetic field amplification.

X-ray polarimetry has the potential to make key-contributions to our understanding of galactic compact objects like binary black hole systems and neutron stars, and extragalactic objects like active galactic nuclei, blazars, and Gamma Ray Bursts. Furthermore, several particle astrophysics topics can be addressed including uniquely sensitive tests of Lorentz invariance. In the energy range from 10 keV to several MeV, Compton polarimeters achieve the best performance. In this paper we evaluate the benefit that comes from using the azimuthal and polar angles of the Compton scattered photons in the analysis, rather than using the azimuthal scattering angles alone. We study the case of an ideal Compton polarimeter and show that a Maximum Likelihood analysis which uses the two scattering angles lowers the Minimum Detectable Polarization (MDP) by ~20% compared to a standard analysis based on the azimuthal scattering angles alone. The accuracies with which the polarization fraction and the polarization direction can be measured improve by a similar amount. We conclude by discussing potential applications of Maximum Likelihood analysis methods for various polarimeter experiments.

CdZnTe (CZT) is a wide bandgap II-VI semiconductor developed for the spectroscopic detection of X-rays and {\gamma}-rays at room temperature. The Swift Burst Alert Telescope is using an 5240 cm2 array of 2 mm thick CZT detectors for the detection of 15-150 keV X-rays from Gamma-Ray Bursts. We report on the systematic tests of thicker (\geq 0.5 cm) CZT detectors with volumes between 2 cm3 and 4 cm3 which are potential detector choices for a number of future X-ray telescopes that operate in the 10 keV to a few MeV energy range. The detectors contacted in our laboratory achieve Full Width Half Maximum energy resolutions of 2.7 keV (4.5%) at 59 keV, 3 keV (2.5%) at 122 keV and 4 keV (0.6%) at 662 keV. The 59 keV and 122 keV energy resolutions are among the world-best results for \geq 0.5 cm thick CZT detectors. We use the data set to study trends of how the energy resolution depends on the detector thickness and on the pixel pitch. Unfortunately, we do not find clear trends, indicating that even for the extremely good energy resolutions reported here, the achievable energy resolutions are largely determined by the properties of individual crystals. Somewhat surprisingly, we achieve the reported results without applying a correction of the anode signals for the depth of the interaction. Measuring the interaction depths thus does not seem to be a pre-requisite for achieving sub-1% energy resolutions at 662 keV.

We report on Monte Carlo studies of the hard X-ray polarimeter X-Calibur. The polarimeter will be used in the focal plane of a grazing incidence hard X-ray telescope. It combines a low-Z Compton scatterer with a high-Z Cadmium Zinc Telluride (CZT) detector assembly to measure the polarization of 10 keV - 80 keV X-rays. X-Calibur makes use of the fact that polarized photons Compton scatter preferentially perpendicular to the electric field orientation. In contrast of competing designs, which use only a small fraction of the incoming X-rays, X-Calibur achieves a high detection efficiency of order unity. In this contributions, we discuss a Monte Carlo study which compares X-Calibur's polarimeteric performance achieved using different scattering materials (Scintillator, Be, LiH, Li), and calculate the sensitivity of X-Calibur when used with different balloon-borne and space-borne mirror assemblies.

The very high energy (VHE; E > 100 GeV) blazar Markarian 501 has a well-studied history of extreme spectral variability and is an excellent laboratory for studying the physical processes within the jets of active galactic nuclei. However, there are few detailed multiwavelength studies of Markarian 501 during its quiescent state, due to its low luminosity. A short-term multiwavelength study of Markarian 501 was coordinated in March 2009, focusing around a multi-day observation with the Suzaku X-ray satellite and including {\gamma}-ray data from VERITAS, MAGIC, and the Fermi Gamma-ray Space Telescope with the goal of providing a well-sampled multiwavelength baseline measurement of Markarian 501 in the quiescent state. The results of these quiescent-state observations are compared to the historically extreme outburst of April 16, 1997, with the goal of examining variability of the spectral energy distribution between the two states. The derived broadband spectral energy distribution shows the characteristic double-peaked profile. We find that the X-ray peak shifts by over two orders of magnitude in photon energy between the two flux states while the VHE peak varies little. The limited shift in the VHE peak can be explained by the transition to the Klein-Nishina regime. Synchrotron self-Compton models are matched to the data and the implied Klein-Nishina effects are explored.

The intermediate-frequency peaked BL Lacertae (IBL) object 3C 66A is detected during 2007 - 2008 in VHE (very high energy: E > 100 GeV) gamma-rays with the VERITAS stereoscopic array of imaging atmospheric Cherenkov telescopes. An excess of 1791 events is detected, corresponding to a significance of 21.2 standard deviations (sigma), in these observations (32.8 hours live time). The observed integral flux above 200 GeV is 6% of the Crab Nebula's flux and shows evidence for variability on the time-scale of days. The measured energy spectrum is characterized by a soft power law with photon index Gamma = 4.1 +- 0.4_stat +- 0.6_sys. The radio galaxy 3C 66B is excluded as a possible source of the VHE emission.

X-ray polarimetry promises to give qualitatively new information about high-energy sources. Examples of interesting source classes are binary black hole systems, rotation and accretion powered neutron stars, Microquasars, Active Galactic Nuclei and Gamma-Ray Bursts. Furthermore, X-ray polarimetry affords the possibility for testing fundamental physics, e.g. to observe signatures of light bending in the strong gravitational field of a black hole, to detect third order Quantum Electrodynamic effects in the magnetosphere of Magnetars, and to perform sensitive tests of Lorentz Invariance. In this paper we discuss scientific drivers of hard (>10 keV) X-ray polarimetry emphasizing how observations in the hard band can complement observations at lower energies (0.1 - 10 keV). Subsequently, we describe four different technical realizations of hard X-ray polarimeters suitable for small to medium sized space borne missions, and study their performance in the signal-dominated case based on Monte Carlo simulations. We end with confronting the instrument requirements for accomplishing the science goals with the capabilities of the four polarimeters.

One of the most intriguing open questions of today's astrophysics is the jet physical properties and the location and the mechanisms for the production of MeV, GeV, and TeV gamma-rays in AGN jets. M87 is a privileged laboratory for a detailed study of the properties of jets, owing to its proximity, its massive black hole, and its conspicuous emission at radio wavelengths and above. We started on November 2009 a monitoring program with the e-EVN at 5 GHz. We present here results of these multi-epoch observations and discuss the two episodes of activity at energy E>100 GeV that occured in this period. One of these observations was obtained at the same day of the first high energy flare. We added to our results literature data obtained with the VLBI and VLA. A clear change in the proper motion velocity of HST-1 is present at the epoch ~ 2005.5. In the time range 1998 -- 2005.5 the apparent velocity is subluminal, and superluminal (~ 2.7c) after 2005.5.

The hard X-ray sky now being studied by INTEGRAL and Swift and soon by NuSTAR is rich with energetic phenomena and highly variable non-thermal phenomena on a broad range of timescales. The High Energy Telescope (HET) on the proposed Energetic X-ray Imaging Survey Telescope (EXIST) mission will repeatedly survey the full sky for rare and luminous hard X-ray phenomena at unprecedented sensitivities. It will detect and localize (<20", at 5 sigma threshold) X-ray sources quickly for immediate followup identification by two other onboard telescopes - the Soft X-ray imager (SXI) and Optical/Infrared Telescope (IRT). The large array (4.5 m^2) of imaging (0.6 mm pixel) CZT detectors in the HET, a coded-aperture telescope, will provide unprecedented high sensitivity (~0.06 mCrab Full Sky in a 2 year continuous scanning survey) in the 5 - 600 keV band. The large field of view (90 deg x 70 deg) and zenith scanning with alternating-orbital nodding motion planned for the first 2 years of the mission will enable nearly continuous monitoring of the full sky. A 3y followup pointed mission phase provides deep UV-Optical-IR-Soft X-ray and Hard X-ray imaging and spectroscopy for thousands of sources discovered in the Survey. We review the HET design concept and report the recent progress of the CZT detector development, which is underway through a series of balloon-borne wide-field hard X-ray telescope experiments, ProtoEXIST. We carried out a successful flight of the first generation of fine pixel large area CZT detectors (ProtoEXIST1) on Oct 9, 2009. We also summarize our future plan (ProtoEXIST2 & 3) for the technology development needed for the HET.

Indirect dark matter searches with ground-based gamma-ray observatories provide an alternative for identifying the particle nature of dark matter that is complementary to that of direct search or accelerator production experiments. We present the results of observations of the dwarf spheroidal galaxies Draco, Ursa Minor, Bootes 1, and Willman 1 conducted by VERITAS. These galaxies are nearby dark matter dominated objects located at a typical distance of several tens of kiloparsecs for which there are good measurements of the dark matter density profile from stellar velocity measurements. Since the conventional astrophysical background of very high energy gamma rays from these objects appears to be negligible, they are good targets to search for the secondary gamma-ray photons produced by interacting or decaying dark matter particles. No significant gamma-ray flux above 200 GeV was detected from these four dwarf galaxies for a typical exposure of ~20 hours. The 95% confidence upper limits on the integral gamma-ray flux are in the range 0.4-2.2x10^-12 photons cm^-2s^-1. We interpret this limiting flux in the context of pair annihilation of weakly interacting massive particles and derive constraints on the thermally averaged product of the total self-annihilation cross section and the relative velocity of the WIMPs. The limits are obtained under conservative assumptions regarding the dark matter distribution in dwarf galaxies and are approximately three orders of magnitude above the generic theoretical prediction for WIMPs in the minimal supersymmetric standard model framework. However significant uncertainty exists in the dark matter distribution as well as the neutralino cross sections which under favorable assumptions could further lower the limits.

We report on the findings of a 364 ksec observation of the BL LAC object Mrk 421 with the X-ray observatory Suzaku. The analysis in this paper uses fluxes and hardness ratios in the broad energy range from 0.5 keV to 30 keV. During the course of the observation, the 0.5 keV - 30 keV flux decreased by a factor of $\sim$2 and was accompanied by several large flares occurring on timescales of a few hours. We find that fitting a broken power model to spectra from isolated epochs during the observation describes the data well. Different flares exhibit different spectral and hardness ratio evolutions. The cumulative observational evidence indicates that the particle acceleration mechanism in the Mrk 421 jet produces electron energy distributions with a modest range of spectral indices and maximum energies. We argue that the short-timescale X-ray spectral variability in the flares can be attributed mostly to intrinsic changes in the acceleration process, dominating other influences such as fluctuations in the Doppler beaming factor, or radiative cooling in or outside the acceleration zone.

We report the discovery of very high energy gamma-ray emission from the direction of the SNR G54.1+0.3 using the VERITAS ground-based gamma-ray observatory. The TeV signal has an overall significance of 6.8$\sigma$ and appears point-like given the 5$^{arcminute}$ resolution of the instrument. The integral flux above 1 TeV is 2.5% of the Crab Nebula flux and significant emission is measured between 250 GeV and 4 TeV, well described by a power-law energy spectrum dN/dE $\sim$ E$^{-\Gamma}$ with a photon index $\Gamma= 2.39\pm0.23_{stat}\pm0.30_{sys}$. We find no evidence of time variability among observations spanning almost two years. Based on the location, the morphology, the measured spectrum, the lack of variability and a comparison with similar systems previously detected in the TeV band, the most likely counterpart of this new VHE gamma-ray source is the PWN in the SNR G54.1+0.3. The measured X-ray to VHE gamma-ray luminosity ratio is the lowest among all the nebulae supposedly driven by young rotation-powered pulsars, which could indicate a particle-dominated PWN.

One of the most intriguing open questions of today's astrophysics is the one concerning the location and the mechanisms for the production of MeV, GeV, and TeV gamma-rays in AGN jets. M87 is a privileged laboratory for a detailed study of the properties of jets, owing to its proximity, its massive black hole, and its conspicuous emission at radio wavelengths and above. We started on November 2009 a monitoring program with the e-EVN at 5 GHz, during which two episodes of activity at energy E > 100 GeV have occured. We present here results of these multi-epoch observations. The inner jet and HST-1 are both detected and resolved in our datasets. One of these observations was obtained at the same day of the first high energy flare. A clear change in the proper motion velocity of HST-1 is present at the epoch ~2005.5. In the time range 2003 -- 2005.5 the apparent velocity is subluminal, and superluminal (~ 2.7c) after 2005.5.

We report on the design and estimated performance of a balloon-borne hard X-ray polarimeter called HX-POL. The experiment uses a combination of Si and Cadmium Zinc Telluride detectors to measure the polarization of 50 keV-400 keV X-rays from cosmic sources through the dependence of the angular distribution of Compton scattered photons on the polarization direction. On a one-day balloon flight, HX-POL would allow us to measure the polarization of bright Crab-like sources for polarization degrees well below 10%. On a longer (15-30 day) flight from Australia or Antarctica, HX-POL would be be able to measure the polarization of bright galactic X-ray sources down to polarization degrees of a few percent. Hard X-ray polarization measurements provide unique venues for the study of particle acceleration processes by compact objects and relativistic outflows. In this paper, we discuss the overall instrument design and performance. Furthermore, we present results from laboratory tests of the Si and CZT detectors.

M87 is a nearby radio galaxy that is detected at energies ranging from radio to VHE gamma-rays. Its proximity and its jet, misaligned from our line-of-sight, enable detailed morphological studies and extensive modeling at radio, optical, and X-ray energies. Flaring activity was observed at all energies, and multi-wavelength correlations would help clarify the origin of the VHE emission. In this paper, we describe a detailed temporal and spectral analysis of the VERITAS VHE gamma-ray observations of M87 in 2008 and 2009. In the 2008 observing season, VERITAS detected an excess with a statistical significance of 7.2 sigma from M87 during a joint multi-wavelength monitoring campaign conducted by three major VHE experiments along with the Chandra X-ray Observatory. In February 2008, VERITAS observed a VHE flare from M87 occurring over a 4-day timespan. The peak nightly flux above 250GeV was 7.7% of the Crab Nebula flux. M87 was marginally detected before this 4-day flare period, and was not detected afterwards. Spectral analysis of the VERITAS observations showed no significant change in the photon index between the flare and pre-flare states. Shortly after the VHE flare seen by VERITAS, the Chandra X-ray Observatory detected the flux from the core of M87 at a historical maximum, while the flux from the nearby knot HST-1 remained quiescent. Acciari et al. (2009) presented the 2008 contemporaneous VHE gamma-ray, Chandra X-ray, and VLBA radio observations which suggest the core as the most likely source of VHE emission, in contrast to the 2005 VHE flare that was simultaneous with an X-ray flare in the HST-1 knot. In 2009, VERITAS continued its monitoring of M87 and marginally detected a 4.2 sigma excess corresponding to a flux of ~1% of the Crab Nebula. No VHE flaring activity was observed in 2009.