We present high cadence multi-frequency radio observations of the long Gamma-Ray Burst (GRB) 190829A, which was detected at photon energies above 100 GeV by the High Energy Stereoscopic System (H.E.S.S.). Observations with the Meer Karoo Array Telescope (MeerKAT, 1.3 GHz), and Arcminute Microkelvin Imager - Large Array (AMI-LA, 15.5 GHz) began one day post-burst and lasted nearly 200 days. We used complementary data from Swift X-Ray Telescope (XRT), which ran to 100 days post-burst. We detected a likely forward shock component with both MeerKAT and XRT up to over 100 days post-burst. Conversely, the AMI-LA light curve appears to be dominated by reverse shock emission until around 70 days post-burst when the afterglow flux drops below the level of the host galaxy. We also present previously unpublished observations of the other H.E.S.S.-detected GRB, GRB 180720B from AMI-LA, which shows likely forward shock emission that fades in less than 10 days. We present a comparison between the radio emission from the three GRBs with detected very high energy (VHE) gamma-ray emission and a sensitivity-limited radio afterglow sample. GRB 190829A has the lowest isotropic radio luminosity of any GRB in our sample, but the distribution of luminosities is otherwise consistent, as expected, with the VHE GRBs being drawn from the same parent distribution as the other radio-detected long GRBs.

ThunderKAT is the image-plane transients programme for MeerKAT. The goal as outlined in 2010, and still today, is to find, identify and understand high-energy astrophysical processes via their radio emission (often in concert with observations at other wavelengths). Through a comprehensive and complementary programme of surveying and monitoring Galactic synchrotron transients (across a range of compact accretors and a range of other explosive phenomena) and exploring distinct populations of extragalactic synchrotron transients (microquasars, supernovae and possibly yet unknown transient phenomena) - both from direct surveys and commensal observations - we will revolutionise our understanding of the dynamic and explosive transient radio sky. As well as performing targeted programmes of our own, we have made agreements with the other MeerKAT large survey projects (LSPs) that we will also search their data for transients. This commensal use of the other surveys, which remains one of our key programme goals in 2016, means that the combined MeerKAT LSPs will produce by far the largest GHz-frequency radio transient programme to date.

The blazar 3C 279 exhibited a giant $\gamma$-ray outburst in 2013 December and 2014 April. Apart from the very fast $\gamma$-ray flux variability, the spectral nature of the flares were also found to vary significantly between these two flaring events. A prominent curvature in the $\gamma$-ray spectrum was noticed in 2014 April flare, on the other hand, the 2013 December displayed an extreme hardening of the spectrum. These observations, thus, put strong constraints on our understanding of the underlying particle acceleration mechanisms.

The blazar 3C 279 exhibited twin $\gamma$-ray flares of similar intensity in 2013 December and 2014 April. In this work, we present a detailed multi-wavelength analysis of the 2013 December flaring event. Multi-frequency observations reveal the uncorrelated variability patterns with X-ray and optical-UV fluxes peaking after the $\gamma$-ray maximum. The broadband spectral energy distribution (SED) at the peak of the $\gamma$-ray activity shows a rising $\gamma$-ray spectrum but a declining optical-UV flux. This observation along with the detection of uncorrelated variability behavior rules out the one-zone leptonic emission scenario. We, therefore, adopt two independent methodologies to explain the SED: a time dependent lepto-hadronic modeling and a two-zone leptonic radiative modeling approach. In the lepto-hadronic modeling, a distribution of electrons and protons subjected to a randomly orientated magnetic field produces synchrotron radiation. Electron synchrotron is used to explain the IR to UV emission while proton synchrotron emission is used to explain the high energy $\gamma$-ray emission. A combination of both electron synchrotron self Compton emission and proton synchrotron emission is used to explain the X-ray spectral break seen during the later stage of the flare. In the two-zone modeling, we assume a large emission region emitting primarily in IR to X-rays and $\gamma$-rays to come primarily from a fast moving compact emission region. We conclude by noting that within a span of 4 months, 3C 279 has shown the dominance of a variety of radiative processes over each other and this reflects the complexity involved in understanding the physical properties of blazar jets in general.

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 well studied blazar Markarian 421 (Mrk 421, $z$=0.031) was the subject of an intensive multi-wavelength campaign when it flared in 2013 April. The recorded X-ray and very high energy (VHE, E$>$100 GeV) $\gamma$-ray fluxes are the highest ever measured from this object. At the peak of the activity, it was monitored by the hard X-ray focusing telescope {\it Nuclear Spectroscopic Telescope Array} ({\it NuSTAR}) and {\it Swift} X-Ray Telescope (XRT). In this work, we present a detailed variability analysis of {\it NuSTAR} and {\it Swift}-XRT observations of Mrk 421 during this flaring episode. We obtained the shortest flux doubling time of 14.01$\pm$5.03 minutes, which is the shortest hard X-ray (3$-$79 keV) variability ever recorded from Mrk 421 and is on the order of the light crossing time of the black hole's event horizon. A pattern of extremely fast variability events superposed on slowly varying flares is found in most of the {\it NuSTAR} observations. We suggest that these peculiar variability patterns may be explained by magnetic energy dissipation and reconnection in a fast moving compact emission region within the jet. Based on the fast variability, we derive a lower limit on the magnetic field strength of $B \ge 0.73 \delta_1^{-2/3} \, \nu_{19}^{1/3}$~G, where $\delta_1$ is the Doppler factor in units of 10, and $\nu_{19}$ is the characteristic X-ray synchrotron frequency in units of $10^{19}$~Hz.

We present the results of extensive multi-band intra-night optical monitoring of BL Lacertae during 2010--2012. BL Lacertae was very active in this period and showed intense variability in almost all wavelengths. We extensively observed it for a total for 38 nights; on 26 of them observations were done quasi-simultaneously in B, V, R and I bands (totaling 113 light curves), with an average sampling interval of around 8 minutes. BL Lacertae showed significant variations on hour-like timescales in a total of 19 nights in different optical bands. We did not find any evidence for periodicities or characteristic variability time-scales in the light curves. The intranight variability amplitude is generally greater at higher frequencies and decreases as the source flux increases. We found spectral variations in BL Lacertae in the sense that the optical spectrum becomes flatter as the flux increases but in several flaring states deviates from the linear trend suggesting different jet components contributing to the emission at different times.

We present a detailed investigation of the flaring activity observed from a BL Lac object, S5 0716+714 , during its brightest ever optical state in the second half of January 2015. Observed almost simultaneously in the optical, X-rays and {\gamma}-rays, a significant change in the degree of optical polarization (PD) and a swing in the position angle (PA) of polarization were recorded. A detection in the TeV (VHE) was also reported by the MAGIC consortium during this flaring episode. Two prominent sub-flares, peaking about 5-days apart, were seen in almost all the energy bands. The multi-wavelength light-curves, spectral energy distribution (SED) and polarization are modeled using the time-dependent code developed by Zhang et al. (2014). This model assumes a straight jet threaded by large scale helical magnetic fields taking into account the light travel time effects, incorporating synchrotron flux and polarization in 3D geometry. The rapid variation in PD and rotation in PA are most likely due to re-connections happening in the emission region in the jet, as suggested by the change in the ratio of toroidal to poloidal components of magnetic field during quiescent and flaring states.

To search for optical variability on a wide range of timescales, we have carried out photometric monitoring of two flat spectrum radio quasars, 3C 454.3 and 3C 279, plus one BL Lac, S5 0716+714, all of which have been exhibiting remarkably high activity and pronounced variability at all wavelengths. CCD magnitudes in B, V, R and I pass-bands were determined for $\sim$ 7000 new optical observations from 114 nights made during 2011 - 2014, with an average length of $\sim$ 4 h each, at seven optical telescopes: four in Bulgaria, one in Greece, and two in India. We measured multiband optical flux and colour variations on diverse timescales. Discrete correlation functions were computed among B, V, R, and I observations, to search for any time delays. We found weak correlations in some cases with no significant time lags. The structure function method was used to estimate any characteristic time-scales of variability. We also investigated the spectral energy distribution of the three blazars using B, V, R, I, J and K pass-band data. We found that the sources almost always follows a bluer-when-brighter trend. We discuss possible physical causes of the observed spectral variability.

We investigate potential $\gamma-\gamma$ absorption of gamma-ray emission from blazars arising from inhomogeneities along the line of sight, beyond the diffuse Extragalactic Background Light (EBL). As plausible sources of excess $\gamma-\gamma$ opacity, we consider (1) foreground galaxies, including cases in which this configuration leads to strong gravitational lensing, (2) individual stars within these foreground galaxies, and (3) individual stars within our own galaxy, which may act as lenses for microlensing events. We found that intervening galaxies close to the line-of-sight are unlikely to lead to significant excess $\gamma-\gamma$ absorption. This opens up the prospect of detecting lensed gamma-ray blazars at energies above 10 GeV with their gamma-ray spectra effectively only affected by the EBL. The most luminous stars located either in intervening galaxy or in our galaxy provides an environment in which these gamma-rays could, in principle, be significantly absorbed. However, despite a large microlensing probability due to stars located in intervening galaxies, gamma-rays avoid absorption by being deflected by the gravitational potentials of such intervening stars to projected distances ("impact parameters"') where the resulting $\gamma-\gamma$ opacities are negligible. Thus, neither of the intervening excess photon fields considered here, provide a substantial source of excess $\gamma-\gamma$ opacity beyond the EBL, even in the case of very close alignments between the background blazar and a foreground star or galaxy.

Context. The international whole earth blazar telescope (WEBT) consortium planned and carried out three days of intensive micro-variability observations of S5 0716+714 from February 22, 2009 to February 25, 2009. This object was chosen due to its bright apparent magnitude range, its high declination, and its very large duty cycle for micro-variations. Aims. We report here on the long continuous optical micro-variability light curve of 0716+714 obtained during the multi-site observing campaign during which the Blazar showed almost constant variability over a 0.5 magnitude range. The resulting light curve is presented here for the first time. Observations from participating observatories were corrected for instrumental differences and combined to construct the overall smoothed light curve. Methods. Thirty-six observatories in sixteen countries participated in this continuous monitoring program and twenty of them submitted data for compilation into a continuous light curve. The light curve was analyzed using several techniques including Fourier transform, Wavelet and noise analysis techniques. Those results led us to model the light curve by attributing the variations to a series of synchrotron pulses. Results. We have interpreted the observed microvariations in this extended light curve in terms of a new model consisting of individual stochastic pulses due to cells in a turbulent jet which are energized by a passing shock and cool by means of synchrotron emission. We obtained an excellent fit to the 72-hour light curve with the synchrotron pulse model.

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.

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 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.

B2013+370 and B2023+336 are two blazars at low-galactic latitude that were previously proposed to be the counterparts for the EGRET unidentified sources, 3EG J2016+3657 and 3EG J2027+3429. Gamma-ray emission associated with the EGRET sources has been detected by the Fermi Gamma-ray Space Telescope, and the two sources, 1FGL J2015.7+3708 and 1FGL J2027.6+3335, have been classified as unidentified in the 1-year catalog. This analysis of the Fermi-LAT data collected during 31 months reveals that the 1FGL sources are spatially compatible with the blazars, and are significantly variable, supporting the hypothesis of extragalactic origin for the gamma-ray emission. The gamma-ray light curves are compared with 15 GHz radio light curves from the 40-m telescope at the Owens Valley Radio Observatory (OVRO). Simultaneous variability is seen in both bands for the two blazar candidates. The study is completed with the X-ray analysis of 1FGL J2015.7+3708 using Swift observations that were triggered in August 2010 by a Fermi-detected flare. The resulting spectral energy distribution shows a two-component structure typical of blazars. We also identify a second source in the field of view of 1FGL J2027.6+3335 with similar characteristics to the known LAT pulsars. This study gives solid evidence favoring blazar counterparts for these two unidentified EGRET and Fermi sources, supporting the hypothesis that a number of unidentified gamma-ray sources at low galactic latitudes are indeed of extragalactic origin.

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.

We report results from an intensive multiwavelength campaign on the intermediate-frequency-peaked BL Lacertae object W Com (z=0.102) during a strong outburst of very high energy gamma-ray emission in June 2008. The very high energy gamma-ray signal was detected by VERITAS on 2008 June 7-8 with a flux F(>200 GeV) = (5.7+-0.6)x10^-11 cm-2s-1, about three times brighter than during the discovery of gamma-ray emission from W Com by VERITAS in 2008 March. The initial detection of this flare by VERITAS at energies above 200 GeV was followed by observations in high energy gamma-rays (AGILE, E>100 MeV), and X-rays (Swift and XMM-Newton), and at UV, and ground-based optical and radio monitoring through the GASP-WEBT consortium and other observatories. Here we describe the multiwavelength data and derive the spectral energy distribution (SED) of the source from contemporaneous data taken throughout the flare.

We present results from a long-term monitoring campaign on the TeV binary LSI +61 303 with VERITAS at energies above 500 GeV, and in the 2-10 keV hard X-ray bands with RXTE and Swift, sampling nine 26.5 day orbital cycles between September 2006 and February 2008. The binary was observed by VERITAS to be variable, with all integrated observations resulting in a detection at the 8.8 sigma (2006/2007) and 7.3 sigma (2007/2008) significance level for emission above 500 GeV. The source was detected during active periods with flux values ranging from 5 to 20% of the Crab Nebula, varying over the course of a single orbital cycle. Additionally, the observations conducted in the 2007-2008 observing season show marginal evidence (at the 3.6 sigma significance level) for TeV emission outside of the apastron passage of the compact object around the Be star. Contemporaneous hard X-ray observations with RXTE and Swift show large variability with flux values typically varying between 0.5 and 3.0*10^-11 ergs cm^-2 s^-1 over a single orbital cycle. The contemporaneous X-ray and TeV data are examined and it is shown that the TeV sampling is not dense enough to detect a correlation between the two bands.

Using the data from the Rossi X-ray Timing Explorer (RXTE), we report the different energy dependence of the variability of the four persistent hard X-ray sources in the low-hard state: Cygnus X-1, GX 339-4, GRS 1758-258 and 1E 1740.7-2942. Cygnus X-1 is found to have a flatter power density spectrum (PDS) shape at higher energies. The other three sources have energy independent PDS shapes. The energy dependence of the overall variability (the integrated rms amplitude) varies from source to source and from observation to observation. 1E~1740.7-2942, for example, has a variability generally increasing with energy while GX 339-4 has a decreasing variability. A general trend is found in the four sources that the integrated rms amplitude anti-correlates with the X-ray flux. We compare these distinct energy dependent behaviors with several emission models. None of the models can fully explain all the features that we have found.

GX 339-4 is an unusual black hole candidate in that it is a persistent source, being detected by X-ray telescopes most of the time, but it also has nova-like flaring states. In 1996 we performed radio, optical, X-ray, and gamma-ray observations of GX 339-4 when it was in a hard state (= soft X-ray low state). Here we present a brief summary of some of the results.