The radio source 3C 264, hosted by the giant elliptical galaxy NGC 3862, was observed with VERITAS between February 2017 and May 2019. These deep observations resulted in the discovery of very-high-energy (VHE; E $>100$ GeV) $\gamma$-ray emission from this active galaxy. An analysis of $\sim$57 hours of quality-selected live time yields a detection at the position of the source, corresponding to a statistical significance of 7.8 standard deviations above background. The observed VHE flux is variable on monthly time scales, with an elevated flux seen in 2018 observations. The VHE emission during this elevated state is well-characterized by a power-law spectrum with a photon index $\Gamma = 2.20 \pm 0.27$ and flux F($>315$ GeV) = ($7.6\pm 1.2_{\mathrm stat} \pm 2.3_{\mathrm syst})\times 10^{-13}$ cm$^{-2}$ s$^{-1}$, or approximately 0.7% of the Crab Nebula flux above the same threshold. 3C 264 ($z = 0.0217$) is the most distant radio galaxy detected at VHE, and the elevated state is thought to be similar to that of the famously outbursting jet in M 87. Consequently, extensive contemporaneous multi-wavelength data were acquired in 2018 at the time of the VHE high state. An analysis of these data, including VLBA, VLA, HST, Chandra and Swift observations in addition to the VERITAS data, is presented, along with a discussion of the resulting spectral energy distribution.

We report on variability and correlation studies using multiwavelength observations of the blazar Mrk 421 during the month of February, 2010 when an extraordinary flare reaching a level of $\sim$27~Crab Units above 1~TeV was measured in very-high-energy (VHE) $\gamma$-rays with the VERITAS observatory. This is the highest flux state for Mrk 421 ever observed in VHE $\gamma$-rays. Data are analyzed from a coordinated campaign across multiple instruments including VHE $\gamma$-ray (VERITAS, MAGIC), high-energy (HE) $\gamma$-ray (Fermi-LAT), X-ray (Swift}, RXTE, MAXI), optical (including the GASP-WEBT collaboration and polarization data) and radio (Mets\"ahovi, OVRO, UMRAO). Light curves are produced spanning multiple days before and after the peak of the VHE flare, including over several flare `decline' epochs. The main flare statistics allow 2-minute time bins to be constructed in both the VHE and optical bands enabling a cross-correlation analysis that shows evidence for an optical lag of $\sim$25-55 minutes, the first time-lagged correlation between these bands reported on such short timescales. Limits on the Doppler factor ($\delta \gtrsim 33$) and the size of the emission region ($ \delta^{-1}R_B \lesssim 3.8\times 10^{13}\,\,\mbox{cm}$) are obtained from the fast variability observed by VERITAS during the main flare. Analysis of 10-minute-binned VHE and X-ray data over the decline epochs shows an extraordinary range of behavior in the flux-flux relationship: from linear to quadratic to lack of correlation to anti-correlation. Taken together, these detailed observations of an unprecedented flare seen in Mrk 421 are difficult to explain by the classic single-zone synchrotron self-Compton model.

We report on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on 15 minute time bins, and over three X-ray bands (3-7 keV, 7-30 keV and 30-80 keV) and three very-high-energy (>0.1 TeV, hereafter VHE) gamma-ray bands (0.2-0.4 TeV, 0.4-0.8 TeV and >0.8 TeV). We detected substantial flux variations on multi-hour and sub-hour timescales in all the X-ray and VHE gamma-ray bands. The characteristics of the sub-hour flux variations are essentially energy-independent, while the multi-hour flux variations can have a strong dependence on the energy of the X-ray and the VHE gamma rays. The three VHE bands and the three X-ray bands are positively correlated with no time-lag, but the strength and the characteristics of the correlation changes substantially over time and across energy bands. Our findings favour multi-zone scenarios for explaining the achromatic/chromatic variability of the fast/slow components of the light curves, as well as the changes in the flux-flux correlation on day-long timescales. We interpret these results within a magnetic reconnection scenario, where the multi-hour flux variations are dominated by the combined emission from various plasmoids of different sizes and velocities, while the sub-hour flux variations are dominated by the emission from a single small plasmoid moving across the magnetic reconnection layer.

We report the results of decade-long (2008-2018) $\gamma$-ray to 1 GHz radio monitoring of the blazar 3C 279, including GASP/WEBT, $\it{Fermi}$ and $\it{Swift}$ data, as well as polarimetric and spectroscopic data. The X-ray and $\gamma$-ray light curves correlate well, with no delay > 3 hours, implying general co-spatiality of the emission regions. The $\gamma$-ray-optical flux-flux relation changes with activity state, ranging from a linear to a more complex dependence. The behaviour of the Stokes parameters at optical and radio wavelengths, including 43 GHz VLBA images, supports either a predominantly helical magnetic field or motion of the radiating plasma along a spiral path. Apparent speeds of emission knots range from 10 to 37c, with the highest values requiring bulk Lorentz factors close to those needed to explain $\gamma$-ray variability on very short time scales. The Mg II emission line flux in the `blue' and `red' wings correlates with the optical synchrotron continuum flux density, possibly providing a variable source of seed photons for inverse Compton scattering. In the radio bands we find progressive delays of the most prominent light curve maxima with decreasing frequency, as expected from the frequency dependence of the $\tau=1$ surface of synchrotron self-absorption. The global maximum in the 86 GHz light curve becomes less prominent at lower frequencies, while a local maximum, appearing in 2014, strengthens toward decreasing frequencies, becoming pronounced at $\sim5$ GHz. These tendencies suggest different Doppler boosting of stratified radio-emitting zones in the jet.

We present a multiwavelength study of the flat-spectrum radio quasar CTA 102 during 2013-2017. We use radio-to-optical data obtained by the Whole Earth Blazar Telescope, 15 GHz data from the Owens Valley Radio Observatory, 91 and 103 GHz data from the Atacama Large Millimeter Array, near-infrared data from the Rapid Eye Monitor telescope, as well as data from the Swift (optical-UV and X-rays) and Fermi ($\gamma$ rays) satellites to study flux and spectral variability and the correlation between flux changes at different wavelengths. Unprecedented $\gamma$-ray flaring activity was observed during 2016 November-2017 February, with four major outbursts. A peak flux of (2158 $\pm$ 63)$\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$, corresponding to a luminosity of (2.2 $\pm$ 0.1)$\times$10$^{50}$ erg s$^{-1}$, was reached on 2016 December 28. These four $\gamma$-ray outbursts have corresponding events in the near-infrared, optical, and UV bands, with the peaks observed at the same time. A general agreement between X-ray and $\gamma$-ray activity is found. The $\gamma$-ray flux variations show a general, strong correlation with the optical ones with no time lag between the two bands and a comparable variability amplitude. This $\gamma$-ray/optical relationship is in agreement with the geometrical model that has successfully explained the low-energy flux and spectral behaviour, suggesting that the long-term flux variations are mainly due to changes in the Doppler factor produced by variations of the viewing angle of the emitting regions. The difference in behaviour between radio and higher energy emission would be ascribed to different viewing angles of the jet regions producing their emission.

A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of $\sim$0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was $\sim$3 CU, and the peak of the high-energy spectral component was found to be at $\sim$2 TeV. This study reports very hard X-ray spectra, and the hardest VHE spectra measured to date for Mrk 501. The fractional variability was found to increase with energy, with the highest variability occurring at VHE, and a significant correlation between the X-ray and VHE bands. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency- peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The one-zone synchrotron self-Compton (SSC) scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behaviour seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays.

We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS, Whipple 10-m, Fermi-LAT, RXTE, Swift, GASP-WEBT and VLBA. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found in the acquired data set. The higher variability in the very high energy (>100 GeV, VHE) gamma-ray emission and the lack of correlation with the X-ray emission indicate that the highest-energy electrons that are responsible for the VHE gamma-rays do not make a dominant contribution to the ~1 keV emission. Alternatively, there could be a very variable component contributing to the VHE gamma-ray emission in addition to that coming from the synchrotron self-Compton (SSC) scenarios. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters. We find that there is some degeneracy in both the one-zone and the two-zone SSC scenarios that were probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The SSC model grid-scan shows that the flaring activity around 2009 May 22 cannot be modeled adequately with a one-zone SSC scenario, while it can be suitably described within a two-independent-zone SSC scenario. The observation of an electric vector polarization angle rotation coincident with the gamma-ray flare from 2009 May 1 resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties.

Here we report on the results of the WEBT photo-polarimetric campaign targeting the blazar S5~0716+71, organized in March 2014 to monitor the source simultaneously in BVRI and near IR filters. The campaign resulted in an unprecedented dataset spanning $\sim 110$\,h of nearly continuous, multi-band observations, including two sets of densely sampled polarimetric data mainly in R filter. During the campaign, the source displayed pronounced variability with peak-to-peak variations of about $30\%$ and "bluer-when-brighter" spectral evolution, consisting of a day-timescale modulation with superimposed hourlong microflares characterized by $\sim 0.1$\,mag flux changes. We performed an in-depth search for quasi-periodicities in the source light curve; hints for the presence of oscillations on timescales of $\sim 3$\,h and $\sim 5$\,h do not represent highly significant departures from a pure red-noise power spectrum. We observed that, at a certain configuration of the optical polarization angle relative to the positional angle of the innermost radio jet in the source, changes in the polarization degree led the total flux variability by about 2\,h; meanwhile, when the relative configuration of the polarization and jet angles altered, no such lag could be noted. The microflaring events, when analyzed as separate pulse emission components, were found to be characterized by a very high polarization degree ($> 30\%$) and polarization angles which differed substantially from the polarization angle of the underlying background component, or from the radio jet positional angle. We discuss the results in the general context of blazar emission and energy dissipation models.

Over the past few years, several occasions of large, continuous rotations of the electric vector position angle (EVPA) of linearly polarized optical emission from blazars have been reported. These events are often coincident with high energy gamma-ray flares and they have attracted considerable attention, as they could allow one to probe the magnetic field structure in the gamma-ray emitting region of the jet. The flat-spectrum radio quasar 3C279 is one of the most prominent examples showing this behaviour. Our goal is to study the observed EVPA rotations and to distinguish between a stochastic and a deterministic origin of the polarization variability. We have combined multiple data sets of R-band photometry and optical polarimetry measurements of 3C279, yielding exceptionally well-sampled flux density and polarization curves that cover a period of 2008-2012. Several large EVPA rotations are identified in the data. We introduce a quantitative measure for the EVPA curve smoothness, which is then used to test a set of simple random walk polarization variability models against the data. 3C279 shows different polarization variation characteristics during an optical low-flux state and a flaring state. The polarization variation during the flaring state, especially the smooth approx. 360 degrees rotation of the EVPA in mid-2011, is not consistent with the tested stochastic processes. We conclude that during the two different optical flux states, two different processes govern the polarization variation, possibly a stochastic process during the low-brightness state and a deterministic process during the flaring activity.

After several years of quiescence, the blazar CTA 102 underwent an exceptional outburst in 2012 September-October. The flare was tracked from gamma-ray to near-infrared frequencies, including Fermi and Swift data as well as photometric and polarimetric data from several observatories. An intensive GASP-WEBT collaboration campaign in optical and NIR bands, with an addition of previously unpublished archival data and extension through fall 2015, allows comparison of this outburst with the previous activity period of this blazar in 2004-2005. We find remarkable similarity between the optical and gamma-ray behaviour of CTA 102 during the outburst, with a time lag between the two light curves of ~1 hour, indicative of co-spatiality of the optical and gamma-ray emission regions. The relation between the gamma-ray and optical fluxes is consistent with the SSC mechanism, with a quadratic dependence of the SSC gamma-ray flux on the synchrotron optical flux evident in the post-outburst stage. However, the gamma-ray/optical relationship is linear during the outburst; we attribute this to changes in the Doppler factor. A strong harder-when-brighter spectral dependence is seen both the in gamma-ray and optical non-thermal emission. This hardening can be explained by convexity of the UV-NIR spectrum that moves to higher frequencies owing to an increased Doppler shift as the viewing angle decreases during the outburst stage. The overall pattern of Stokes parameter variations agrees with a model of a radiating blob or shock wave that moves along a helical path down the jet.

We study the multi-band variability and correlations of the TeV blazar Mrk 421 on year time scales, which can bring additional insight on the processes responsible for its broadband emission. We observed Mrk 421 in the very high energy (VHE) gamma-ray range with the Cherenkov telescope MAGIC-I from March 2007 to June 2009 for a total of 96 hours of effective time after quality cuts. The VHE flux variability is quantified with several methods, including the Bayesian Block algorithm, which is applied to data from Cherenkov telescopes for the first time. The 2.3 year long MAGIC light curve is complemented with data from the Swift/BAT and RXTE/ASM satellites and the KVA, GASP-WEBT, OVRO, and Mets\"ahovi telescopes from February 2007 to July 2009, allowing for an excellent characterisation of the multi-band variability and correlations over year time scales. Mrk 421 was found in different gamma-ray emission states during the 2.3 year long observation period. Flares and different levels of variability in the gamma-ray light curve could be identified with the Bayesian Block algorithm. The same behaviour of a quiet and active emission was found in the X-ray light curves measured by Swift/BAT and the RXTE/ASM, with a direct correlation in time. The behaviour of the optical light curve of GASP-WEBT and the radio light curves by OVRO and Mets\"ahovi are different as they show no coincident features with the higher energetic light curves and a less variable emission. The fractional variability is overall increasing with energy. The comparable variability in the X-ray and VHE bands and their direct correlation during both high- and low-activity periods spanning many months show that the electron populations radiating the X-ray and gamma-ray photons are either the same, as expected in the Synchrotron-Self-Compton mechanism, or at least strongly correlated, as expected in electromagnetic cascades.

We present coordinated multiwavelength observations of the bright, nearby BL Lac object Mrk 421 taken in 2013 January-March, involving GASP-WEBT, Swift, NuSTAR, Fermi-LAT, MAGIC, VERITAS, and other collaborations and instruments, providing data from radio to very-high-energy (VHE) gamma-ray bands. NuSTAR yielded previously unattainable sensitivity in the 3-79 keV range, revealing that the spectrum softens when the source is dimmer until the X-ray spectral shape saturates into a steep power law with a photon index of approximately 3, with no evidence for an exponential cutoff or additional hard components up to about 80 keV. For the first time, we observed both the synchrotron and the inverse-Compton peaks of the spectral energy distribution (SED) simultaneously shifted to frequencies below the typical quiescent state by an order of magnitude. The fractional variability as a function of photon energy shows a double-bump structure which relates to the two bumps of the broadband SED. In each bump, the variability increases with energy which, in the framework of the synchrotron self-Compton model, implies that the electrons with higher energies are more variable. The measured multi-band variability, the significant X-ray-to-VHE correlation down to some of the lowest fluxes ever observed in both bands, the lack of correlation between optical/UV and X-ray flux, the low degree of polarization and its significant (random) variations, the short estimated electron cooling time, and the significantly longer variability timescale observed in the NuSTAR light curves point toward in-situ electron acceleration, and suggest that there are multiple compact regions contributing to the broadband emission of Mrk 421 during low-activity states.

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

The occurrence of low-amplitude flux variations in blazars on hourly timescales, commonly known as microvariability, is still a widely debated subject in high-energy astrophysics. Several competing scenarios have been proposed to explain such occurrences, including various jet plasma instabilities leading to the formation of shocks, magnetic reconnection sites, and turbulence. In this letter we present the results of our detailed investigation of a prominent, five-hour-long optical microflare detected during recent WEBT campaign in 2014, March 2-6 targeting the blazar 0716+714. After separating the flaring component from the underlying base emission continuum of the blazar, we find that the microflare is highly polarized, with the polarization degree $\sim (40-60)\%$$\pm (2-10)\%$, and the electric vector position angle $\sim (10 - 20)$deg$\pm (1-8)$deg slightly misaligned with respect to the position angle of the radio jet. The microflare evolution in the $(Q,\,U)$ Stokes parameter space exhibits a looping behavior with a counter-clockwise rotation, meaning polarization degree decreasing with the flux (but higher in the flux decaying phase), and approximately stable polarization angle. The overall very high polarization degree of the flare, its symmetric flux rise and decay profiles, and also its structured evolution in the $Q-U$ plane, all imply that the observed flux variation corresponds to a single emission region characterized by a highly ordered magnetic field. As discussed in the paper, a small-scale but strong shock propagating within the outflow, and compressing a disordered magnetic field component, provides a natural, though not unique, interpretation of our findings.

A flare from the TeV blazar Mrk 421, occurring in March 2010, was observed for 13 consecutive days from radio to very high energy (VHE, E > 100 GeV) gamma-rays with MAGIC, VERITAS, Whipple, FermiLAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We model the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (SSC) models, investigate the physical parameters, and evaluate whether the observed broadband SED variability can be associated to variations in the relativistic particle population. Flux variability was remarkable in the X-ray and VHE bands while it was minor or not significant in the other bands. The one-zone SSC model can describe reasonably well the SED of each day for the 13 consecutive days. This flaring activity is also very well described by a two-zone SSC model, where one zone is responsible for the quiescent emission while the other smaller zone, which is spatially separated from the first one, contributes to the daily-variable emission occurring in X-rays and VHE gamma-rays. Both the one-zone SSC and the two-zone SSC models can describe the daily SEDs via the variation of only four or five model parameters, under the hypothesis that the variability is associated mostly to the underlying particle population. This shows that the particle acceleration and cooling mechanism producing the radiating particles could be the main one responsible for the broadband SED variations during the flaring episodes in blazars. The two-zone SSC model provides a better agreement to the observed SED at the narrow peaks of the low- and high-energy bumps during the highest activity, although the reported one-zone SSC model could be further improved by the variation of the parameters related to the emitting region itself ($\delta$, $B$ and $R$), in addition to the parameters related to the particle population.

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.

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.

Quasar 3C 279 is known to exhibit episodes of optical polarization angle rotation. We present new, well-sampled optical polarization data for 3C 279 and introduce a method to distinguish between random and deterministic electric vector position angle (EVPA) variations. We observe EVPA rotations in both directions with different amplitudes and find that the EVPA variation shows characteristics of both random and deterministic cases. Our analysis indicates that the EVPA variation is likely dominated by a random process in the low brightness state of the jet and by a deterministic process in the flaring state.

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

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.

Since 2005, the blazar 3C 454.3 has shown remarkable flaring activity at all frequencies, and during the last four years it has exhibited more than one gamma-ray flare per year, becoming the most active gamma-ray blazar in the sky. We present for the first time the multi-wavelength AGILE, SWIFT, INTEGRAL, and GASP-WEBT data collected in order to explain the extraordinary gamma-ray flare of 3C 454.3 which occurred in November 2010. On 2010 November 20 (MJD 55520), 3C 454.3 reached a peak flux (E>100 MeV) of F_gamma(p) = (6.8+-1.0)E-5 ph/cm2/s on a time scale of about 12 hours, more than a factor of 6 higher than the flux of the brightest steady gamma-ray source, the Vela pulsar, and more than a factor of 3 brighter than its previous super-flare on 2009 December 2-3. The multi-wavelength data make a thorough study of the present event possible: the comparison with the previous outbursts indicates a close similarity to the one that occurred in 2009. By comparing the broadband emission before, during, and after the gamma-ray flare, we find that the radio, optical and X-ray emission varies within a factor 2-3, whereas the gamma-ray flux by a factor of 10. This remarkable behavior is modeled by an external Compton component driven by a substantial local enhancement of soft seed photons.

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.

The GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT) monitored BL Lacertae in 2008-2009 at radio, near-IR, and optical frequencies. During this period, high-energy observations were performed by XMM-Newton, Swift, and Fermi. We analyse these data with particular attention to the calibration of Swift UV data, and apply a helical jet model to interpret the source broad-band variability. The GASP-WEBT observations show an optical flare in 2008 February-March, and oscillations of several tenths of mag on a few-day time scale afterwards. The radio flux is only mildly variable. The UV data from both XMM-Newton and Swift seem to confirm a UV excess that is likely caused by thermal emission from the accretion disc. The X-ray data from XMM-Newton indicate a strongly concave spectrum, as well as moderate flux variability on an hour time scale. The Swift X-ray data reveal fast (interday) flux changes, not correlated with those observed at lower energies. We compare the spectral energy distribution (SED) corresponding to the 2008 low-brightness state, which was characterised by a synchrotron dominance, to the 1997 outburst state, where the inverse-Compton emission was prevailing. A fit with an inhomogeneous helical jet model suggests that two synchrotron components are at work with their self inverse-Compton emission. Most likely, they represent the radiation from two distinct emitting regions in the jet. We show that the difference between the source SEDs in 2008 and 1997 can be explained in terms of pure geometrical variations. The outburst state occurred when the jet-emitting regions were better aligned with the line of sight, producing an increase of the Doppler beaming factor. Our analysis demonstrates that the jet geometry can play an extremely important role in the BL Lacertae flux and spectral variability.

During the month of December, 2009 the blazar 3C 454.3 became the brightest gamma-ray source in the sky, reaching a peak flux F ~2000E-8 ph/cm2/s for E > 100 MeV. Starting in November, 2009 intensive multifrequency campaigns monitored the 3C 454 gamma-ray outburst. Here we report the results of a 2-month campaign involving AGILE, INTEGRAL, Swift/XRT, Swift/BAT, RossiXTE for the high-energy observations, and Swift/UVOT, KANATA, GRT, REM for the near-IR/optical/UV data. The GASP/WEBT provided radio and additional optical data. We detected a long-term active emission phase lasting ~1 month at all wavelengths: in the gamma-ray band, peak emission was reached on December 2-3, 2009. Remarkably, this gamma-ray super-flare was not accompanied by correspondingly intense emission in the optical/UV band that reached a level substantially lower than the previous observations in 2007-2008. The lack of strong simultaneous optical brightening during the super-flare and the determination of the broad-band spectral evolution severely constrain the theoretical modelling. We find that the pre- and post-flare broad-band behavior can be explained by a one-zone model involving SSC plus external Compton emission from an accretion disk and a broad-line region. However, the spectra of the Dec. 2-3, 2009 super-flare and of the secondary peak emission on Dec. 9, 2009 cannot be satisfactorily modelled by a simple one-zone model. An additional particle component is most likely active during these states.

Since September 2005, the Whipple 10m Gamma-ray Telescope has been operated primarily as a blazar monitor. The five Northern Hemisphere blazars that have already been detected at the Whipple Observatory, Markarian 421, H1426+428, Markarian 501, 1ES 1959+650 and 1ES 2344+514, are monitored routinely each night that they are visible. We report on the Markarian 421 observations taken from November 2005 to June 2006 in the gamma-ray, X-ray, optical and radio bands. During this time, Markarian 421 was found to be variable at all wavelengths probed. Both the variability and the correlations among different energy regimes are studied in detail here. A tentative correlation, with large spread, was measured between the X-ray and gamma-ray bands, while no clear correlation was evident among the other energy bands. In addition to this, the well-sampled spectral energy distribution of Markarian 421 (1101+384) is presented for three different activity levels. The observations of the other blazar targets will be reported separately.

In September 2005, the observing program of the Whipple 10 m gamma-ray telescope was redefined to be dedicated almost exclusively to AGN monitoring. Since then the five Northern Hemisphere blazars that had already been detected at Whipple are monitored routinely each night that they are visible. Thanks to the efforts of a large number of multiwavelength collaborators, the first year of this program has been very successful. We report here on the analysis of Markarian 421 observations taken from November, 2005 to May, 2006 in the gamma-ray, X-ray, optical and radio bands.

We report on a multiwavelength campaign on the TeV gamma-ray blazar Markarian (Mrk) 421 performed during December 2002 and January 2003. These target of opportunity observations were initiated by the detection of X-ray and TeV gamma-ray flares with the All Sky Monitor (ASM) on board the Rossi X-ray Timing Explorer (RXTE) and the 10 m Whipple gamma-ray telescope.The campaign included observational coverage in the radio (University of Michigan Radio Astronomy Observatory), optical (Boltwood, La Palma KVA 0.6m, WIYN 0.9m), X-ray (RXTE pointed telescopes), and TeV gamma-ray (Whipple and HEGRA) bands. At TeV energies, the observations revealed several flares at intermediate flux levels, peaking between 1 and 1.5 times the flux from the Crab Nebula. While the time averaged spectrum can be fitted with a single power law of photon index Gamma =2.8, we find some evidence for spectral variability. Confirming earlier results, the campaign reveals a rather loose correlation between the X-ray and TeV gamma-ray fluxes. In one case, a very strong X-ray flare is not accompanied by a comparable TeV gamma-ray flare. Although the source flux was variable in the optical and radio bands, the sparse sampling of the optical and radio light curves does not allow us to study the correlation properties in detail. We present a simple analysis of the data with a synchrotron-self Compton model, emphasizing that models with very high Doppler factors and low magnetic fields can describe the data.

Following the detection of strong TeV gamma-ray flares from the BL Lac object 1ES 1959+650 with the Whipple 10 m Cherenkov telescope on May 16 and 17, 2002, we performed intensive Target of Opportunity (ToO) radio, optical, X-ray and TeV gamma-ray observations from May 18, 2002 to August 14, 2002. Observations with the X-ray telescope RXTE and the Whipple and HEGRA gamma-ray telescopes revealed several strong flares, enabling us to sensitively test the X-ray/gamma-ray flux correlation properties. Although the X-ray and gamma-ray fluxes seemed to be correlated in general, we found an ``orphan'' gamma-ray flare that was not accompanied by an X-ray flare. After describing in detail the radio (UMRAO, VLA), optical (Boltwood, Abastumani), X-ray (RXTE) and gamma-ray (Whipple, HEGRA) light curves and Spectral Energy Distributions (SEDs) we present initial modeling of the SED with a simple Synchrotron Self-Compton (SSC) model. With the addition of another TeV blazar with good broadband data, we consider the set of all TeV blazars to begin to look for a connection of the jet properties to the properties of the central accreting black hole thought to drive the jet. Remarkably, the temporal and spectral X-ray and gamma-ray emission characteristics of TeV blazars are very similar, even though the masses estimates of their central black holes differ by up to one order of magnitude.

New optical and radio data on the BL Lacertae object AO 0235+16 have been collected in the last four years by a wide international collaboration, which confirm the intense activity of this source. The optical data also include the results of the Whole Earth Blazar Telescope (WEBT) first-light campaign organized in November 1997. The optical spectrum is observed to basically steepen when the source gets fainter. We have investigated the existence of typical variability time scales and of possible correlations between the optical and radio emissions by means of visual inspection, Discrete Correlation Function analysis, and Discrete Fourier Transform technique. The major radio outbursts are found to repeat quasi-regularly with a periodicity of about 5.7 years; this period is also in agreement with the occurrence of some of the major optical outbursts, but not all of them.

We report the discovery of X-ray emission from a radio knot at a projected distance of 25" from the nucleus of the Seyfert galaxy, 3C 120. The data were obtained with the ROSAT High Resolution Imager (HRI). Optical upper limits for the knot preclude a simple power law extension of the radio spectrum and we calculate some of the physical parameters for thermal bremsstrahlung and synchrotron self-Compton models. We conclude that no simple model is consistent with the data but if the knot contains small regions with flat spectra, these could produce the observed X-rays (via synchrotron emission) without being detected at other wavebands.