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 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 the first results from a reverberation-mapping campaign undertaken during the first half of 2012, with additional data on one AGN (NGC 3227) from a 2014 campaign. Our main goals are (1) to determine the black hole masses from continuum-Hbeta reverberation signatures, and (2) to look for velocity-dependent time delays that might be indicators of the gross kinematics of the broad-line region. We successfully measure Hbeta time delays and black hole masses for five AGNs, four of which have previous reverberation mass measurements. The values measured here are in agreement with earlier estimates, though there is some intrinsic scatter beyond the formal measurement errors. We observe velocity dependent Hbeta lags in each case, and find that the patterns have changed in the intervening five years for three AGNs that were also observed in 2007.

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.

The size, shape and degree of emptiness of void interiors sheds light on the details of galaxy formation. A particularly interesting question is whether void interiors are completely empty or contain a dwarf population. However the nearby voids that are most conducive for dwarf searches have large angular diameters, on the order of a steradian, making it difficult to redshift-map a statistically significant portion of their volume to the magnitude limit of dwarf galaxies. As part of addressing this problem, we investigate here the usefulness of number counts in establishing the best locations to search inside nearby (d < 300 Mpc) galaxy voids, utilizing Wolf plots of log(n < m) vs. m as the basic diagnostic. To illustrate expected signatures, we consider the signature of three void profiles, "cut out", "built up", and "universal profile" carved into Monte-Carlo Schechter function models. We then investigate the signatures of voids in the Millennium Run dark matter simulation and the Sloan Digital Sky Survey. We find in all of these the evidence for cut-out and built-up voids is most discernible when the void diameter is similar to the distance to its center. However the density distribution of the universal profile that is characteristic of actual voids is essentially undetectable at any distance. A useful corollary of this finding is that galaxy counts are a reliable measure of survey completeness and stellar contamination even when sampling through significant voids.

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.

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 have measured Johnson BV and Cousins RI magnitudes for comparison stars near five TeV blazars. We compare our values with published values, spanning 25 years in some cases, to identify those stars that are most likely proven stable. To avoid zero-point offsets mimicking long-term variability, we based our analysis on the standard deviation between measurements after a mean offset between data sets was removed. We found most stars to be stable at the 0.04 magnitude level. We confirm two stars as variable and identify two others as possibly being variable. In each of the five fields there are at least two stars, and typically many more, that show no evidence of variability.

Observatories constructed solely for photometric monitoring make it possible to understand the temporal nature of objects over time scales that historically have been difficult to achieve. We report on one such observatory, the Remote Observatory for Variable Object Research (ROVOR) which was constructed to enable both long-term and rapid cadence observations of brighter objects. ROVOR is an 0.4m optical telescope located in central Utah and commissioned for scientific observations in 2008. Principle research has been monitoring blazars, x-ray binaries, AGN, and an occasional gamma-ray burst afterglow. We describe the observatory, the control system, and its unique roof.

The Seyfert 1 galaxy Zw 229-015 is among the brightest active galaxies being monitored by the Kepler mission. In order to determine the black hole mass in Zw 229-015 from H-beta reverberation mapping, we have carried out nightly observations with the Kast Spectrograph at the Lick 3m telescope during the dark runs from June through December 2010, obtaining 54 spectroscopic observations in total. We have also obtained nightly V-band imaging with the Katzman Automatic Imaging Telescope at Lick Observatory and with the 0.9m telescope at the Brigham Young University West Mountain Observatory over the same period. We detect strong variability in the source, which exhibited more than a factor of 2 change in broad H-beta flux. From cross-correlation measurements, we find that the H-beta light curve has a rest-frame lag of 3.86(+0.69,-0.90) days with respect to the V-band continuum variations. We also measure reverberation lags for H-alpha and H-gamma and find an upper limit to the H-delta lag. Combining the H-beta lag measurement with a broad H-beta width of sigma = 1590+/-47 km/s measured from the root-mean-square variability spectrum, we obtain a virial estimate of M_BH = 1.00(-0.24,+0.19)*10^7 solar masses for the black hole in Zw 229-015. As a Kepler target, Zw 229-015 will eventually have one of the highest-quality optical light curves ever measured for any active galaxy, and the black hole mass determined from reverberation mapping will serve as a benchmark for testing relationships between black hole mass and continuum variability characteristics in active galactic nuclei.

The KPNO International Spectroscopic Survey (KISS) is an objective-prism survey designed to detect extragalactic emission-line objects. It combines many of the features of previous slitless spectroscopic surveys with the advantages of modern CCD detectors, and is the first purely digital objective-prism survey for emission-line galaxies (ELGs). Here we present the third list of ELG candidates selected from our red spectral data, which cover the wavelength range 6400 to 7200 A. In most cases, the detected emission line is H-alpha. The current survey list covers the region of the NOAO Deep Wide-Field Survey (NDWFS). This survey covers two fields; the first is 3 x 3 degrees square and located at RA = 14h 30m, DEC = 34.5 deg (B1950), the second is 2.3 x 4.0 degrees and centered at RA = 2h 7.5m, DEC = -4.75 deg. A total area of 19.65 deg^2 is covered by the KISS data. A total of 261 candidate emission-line objects have been selected for inclusion in the survey list (13.3 per deg^2). We tabulate accurate coordinates and photometry for each source, as well as estimates of the redshift, emission-line flux and line equivalent width based on measurements of the digital objective-prism spectra. The properties of the KISS ELGs are examined using the available observational data. When combined with the wealth of multi-wavelength data already available for the NDWFS fields, the current list of KISS ELGs should provide a valuable tool for studying star-formation and nuclear activity in galaxies in the local universe.

The KPNO International Spectroscopic Survey (KISS) is an objective-prism survey for extragalactic emission-line objects. It combines many of the features of previous slitless spectroscopic surveys with the advantages of modern CCD detectors, and is the first purely digital objective-prism survey for emission-line galaxies. Here we present the second list of emission-line galaxy candidates selected from our red spectral data, which cover the wavelength range 6400 to 7200 A. In most cases, the detected emission line is H-alpha. The current survey list covers a 1.6-degree-wide strip located at Dec(1950) = 43d 30' and spans the RA range 11h 55m to 16h 15m. The survey strip runs through the center of the Bootes Void, and has enough depth to adequately sample the far side of the void. An area of 65.8 sq. deg. is covered. A total of 1029 candidate emission-line objects have been selected for inclusion in the survey list (15.6 per sq. deg.). We tabulate accurate coordinates and photometry for each source, as well as estimates of the redshift and emission-line flux and equivalent width based on measurements of the digital objective-prism spectra. The properties of the KISS emission-line galaxies are examined using the available observational data. Although the current survey covers only a modest fraction of the total volume of the Bootes Void, we catalog at least twelve objects that appear to be located within the void. Only one of these objects has been recognized previously as a void galaxy.

The KPNO International Spectroscopic Survey (KISS) is an objective-prism survey for extragalactic emission-line objects. It combines many of the features of previous slitless spectroscopic surveys with the advantages of modern CCD detectors, and is the first purely digital objective-prism survey for emission-line galaxies. Here we present the first list of emission-line galaxy candidates selected from our blue spectral data, which cover the wavelength range 4800 -- 5500 A. In most cases, the detected emission line is [O III]5007. The current survey list covers a one-degree-wide strip located at Dec = 29 deg 30' (B1950.0) and spanning the right ascension range 8h 30m to 17h 0m. An area of 116.6 deg^2 is covered. A total of 223 candidate emission-line objects have been selected for inclusion in the survey list (1.91 deg^-2). We tabulate accurate coordinates and photometry for each source, as well as estimates of the redshift, emission-line flux, and equivalent width based on measurements of the digital objective-prism spectra. The median apparent magnitude of the sample is B = 18.2, and galaxies with redshifts approaching z = 0.1 are detected. The properties of the KISS emission-line galaxies are examined using the available observational data, and compared to previous surveys carried out with photographic plates as well as with the H-alpha-selected portion of KISS.

The KPNO International Spectroscopic Survey (KISS) is a new objective-prism survey for extragalactic emission-line objects. It combines many of the features of previous slitless spectroscopic surveys with the advantages of modern CCD detectors, and is the first purely digital objective-prism survey for emission-line galaxies. Here we present the first list of emission-line galaxy candidates selected from our red spectral data, which cover the spectral range 6400 to 7200 Angstroms. In most cases, the detected emission line is H-alpha. The current survey list covers a one-degree-wide strip located at Dec (1950) = 29.5 degrees and spanning the RA range 12h 15m to 17h 0m. An area of 62.2 sq. degrees is covered. A total of 1128 candidate emission-line objects have been selected for inclusion in the survey list (18.1 per sq. degree). We tabulate accurate coordinates and photometry for each source, as well as estimates of the redshift and emission-line flux and equivalent width based on measurements of the digital objective-prism spectra. The properties of the KISS emission-line galaxies are examined using the available observational data.

The KPNO International Spectroscopic Survey (KISS) is a new objective-prism survey for extragalactic emission-line objects. It combines many of the features of previous slitless spectroscopic surveys that were carried out with Schmidt telescopes using photographic plates with the advantages of modern CCD detectors. It is the first purely digital objective-prism survey, and extends previous photographic surveys to substantially fainter flux limits. In this, the first paper in the series, we give an overview of the survey technique, describe our data processing procedures, and present examples of the types of objects found by KISS. Our first H-alpha-selected survey list detects objects at the rate of 18.1 per square degree, which is 181 times higher than the surface density of the Markarian survey.Since the sample is line-selected, there is an imposed redshift limit of z ~ 0.095 due to the filter employed for the objective-prism observations. We evaluate the quality of the observed parameters derived from the survey data, which include accurate astrometry, photometry, redshifts, and line fluxes. Finally, we describe some of the many applications the KISS database will have for addressing specific questions in extragalactic astronomy. Subsequent papers in this series will present our survey lists of emission-line galaxy candidates.