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.

The XENON1T experiment uses a time projection chamber (TPC) with liquid Xenon to search for Weakly Interacting Massive Particles (WIMPs), a proposed Dark Matter particle, via direct detection. As this experiment relies on capturing rare events, the focus is on achieving a high recall of WIMP events. Hence the ability to distinguish between WIMP and the background is extremely important. To accomplish this, we suggest using Convolutional Neural Networks (CNNs); a Machine Learning procedure mainly used in image recognition tasks. To explore this technique we use XENON collaboration open-source software to simulate the TPC graphical output of Dark Matter signals and main backgrounds. A CNN turns out to be a suitable tool for this purpose, as it can identify features in the images that differentiate the two types of events without the need to manipulate or remove data in order to focus on a particular region of the detector. We find that the CNN can distinguish between the dominant background events (ER) and 500 GeV WIMP events with an accuracy of 87.0\% and a recall of 88.2\%.

We have conducted a multiwavelength survey of 42 radio loud narrow-1ine Seyfert 1 galaxies (RLNLS1s), selected by searching among all the known sources of this type and omitting those with steep radio spectra. We analyse data from radio frequencies to X-rays, and supplement these with information available from online catalogs and the literature in order to cover the full electromagnetic spectrum. This is the largest known multiwavelength survey for this type of source. We detected 90% of the sources in X-rays and found 17% at gamma rays. Extreme variability at high energies was also found, down to timescales as short as hours. In some sources, dramatic spectral and flux changes suggest interplay between a relativistic jet and the accretion disk. The estimated masses of the central black holes are in the range $\sim 10^{6-8}M_{\odot}$, smaller than those of blazars, while the accretion luminosities span a range from $\sim 0.01$ to $\sim 0.49$ times the Eddington limit, similar to those of quasars. The distribution of the calculated jet power spans a range from $\sim 10^{42.6}$ to $\sim 10^{45.6}$ erg s$^{-1}$, generally lower than quasars and BL Lac objects, but partially overlapping with the latter. Once normalised by the mass of the central black holes, the jet power of the three types of active galactic nuclei are consistent with each other, indicating the scalability of the jet. Despite the observational differences, the central engine of RLNLS1s is apparently quite similar to that of blazars. The historical difficulties in finding radio-loud narrow-line Seyfert 1 galaxies might be due to their low power and to intermittent jet activity.

We apply classical machine vision and machine deep learning methods to prototype signal classifiers for the search for extraterrestrial intelligence. Our novel approach uses two-dimensional spectrograms of measured and simulated radio signals bearing the imprint of a technological origin. The studies are performed using archived narrow-band signal data captured from real-time SETI observations with the Allen Telescope Array and a set of digitally simulated signals designed to mimic real observed signals. By treating the 2D spectrogram as an image, we show that high quality parametric and non-parametric classifiers based on automated visual analysis can achieve high levels of discrimination and accuracy, as well as low false-positive rates. The (real) archived data were subjected to numerous feature-extraction algorithms based on the vertical and horizontal image moments and Huff transforms to simulate feature rotation. The most successful algorithm used a two-step process where the image was first filtered with a rotation, scale and shift-invariant affine transform followed by a simple correlation with a previously defined set of labeled prototype examples. The real data often contained multiple signals and signal ghosts, so we performed our non-parametric evaluation using a simpler and more controlled dataset produced by simulation of complex-valued voltage data with properties similar to the observed prototypes. The most successful non-parametric classifier employed a wide residual (convolutional) neural network based on pre-existing classifiers in current use for object detection in ordinary photographs. These results are relevant to a wide variety of research domains that already employ spectrogram analysis from time-domain astronomy to observations of earthquakes to animal vocalization analysis.

Note: This is a revised version of the paper that _corrects_a_calculation_error in translating observed Jansky units to EIRP in Watts. Mistakes are labeled below. Motivated by the hypothesis that Oumuamua could conceivably be an interstellar probe, we used the Allen Telescope Array to search for radio transmissions that would indicate a non-natural origin for this object. Observations were made at radio frequencies between 1 and 10 GHz using the Array's correlator receiver with a channel bandwidth of 100 kHz. In frequency regions not corrupted by man-made interference, we find no signal flux with frequency-dependent lower limits of 0.01 Jy at 1 GHz and 0.1 Jy at 7 GHz. For a putative isotropic object, these limits correspond to transmitter powers of (was mistakenly 30 mW) 10 W and (was mistakenly 300 mW) 100 W, respectively. In frequency ranges that are heavily utilized for satellite communications, our sensitivity to weak signals is badly impinged, but we can still place an upper limit of (was mistakenly 10 W) 3 kW for a transmitter on the asteroid. For comparison and validation should a transmitter be discovered, contemporaneous measurements were made on the solar system asteroids 2017 UZ and 2017 WC with comparable sensitivities. Because they are closer to Earth, we place upper limits on transmitter power to be 0.1 and 0.001 times the limits for Oumuamua, respectively.

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 describe a new approach and algorithm for the detection of artificial signals and their classification in the search for extraterrestrial intelligence (SETI). The characteristics of radio signals observed during SETI research are often most apparent when those signals are represented as spectrograms. Additionally, many observed signals tend to share the same characteristics, allowing for sorting of the signals into different classes. For this work, complex-valued time-series data were simulated to produce a corpus of 140,000 signals from seven different signal classes. A wide residual neural network was then trained to classify these signal types using the gray-scale 2D spectrogram representation of those signals. An average $F_1$ score of 95.11\% was attained when tested on previously unobserved simulated signals. We also report on the performance of the model across a range of signal amplitudes.

We report the first results of a survey on 74 narrow-line Seyfert 1 galaxies (NLS1s) carried out in 2015 with the Karl G. Jansky Very Large Array (JVLA) at 5 GHz in A-configuration. So far, this is the largest survey aimed to image the radio continuum of NLS1s. We produced radio maps in order to compare the general properties of three different samples of objects: radio-quiet NLS1s (RQNLS1s), steep-spectrum radio-loud NLS1s (S-NLS1s), and flat-spectrum radio-loud NLS1s (F-NLS1s). We find that the three classes correspond to different radio morphologies, with F-NLS1s being more compact, and RQNLS1s often showing diffuse emission on kpc scales. We also find that F-NLS1s might be low-luminosity and possibly young blazars, and that S-NLS1s are part of the parent population of F-NLS1s. Dedicated studies to RQNLS1s are needed in order to fully understand their role in the unification pictures.

We present a photometric detection of the first brightness dips of the unique variable star KIC 8462852 since the end of the Kepler space mission in 2013 May. Our regular photometric surveillance started in October 2015, and a sequence of dipping began in 2017 May continuing on through the end of 2017, when the star was no longer visible from Earth. We distinguish four main 1-2.5% dips, named "Elsie," "Celeste," "Skara Brae," and "Angkor", which persist on timescales from several days to weeks. Our main results so far are: (i) there are no apparent changes of the stellar spectrum or polarization during the dips; (ii) the multiband photometry of the dips shows differential reddening favoring non-grey extinction. Therefore, our data are inconsistent with dip models that invoke optically thick material, but rather they are in-line with predictions for an occulter consisting primarily of ordinary dust, where much of the material must be optically thin with a size scale <<1um, and may also be consistent with models invoking variations intrinsic to the stellar photosphere. Notably, our data do not place constraints on the color of the longer-term "secular" dimming, which may be caused by independent processes, or probe different regimes of a single process.

We present the first results of a radio survey of 79 narrow-line Seyfert 1 (NLS1) carried out with the Karl G. Jansky Very Large Array (JVLA) at 5 GHz in A configuration aimed at studying the radio properties of these sources. We report the detection of extended emission in one object: Mrk 783. This is intriguing, since the radio-loudness parameter R of this object is close to the threshold between radio-quiet and radio-loud active galactic nuclei (AGN). The galaxy is one of the few NLS1 showing such an extended emission at z < 0.1. The radio emission is divided in a compact core component and an extended component, observed on both sides of the nucleus and extending from 14 kpc south-east to 12 kpc north-west. There is no sign of a collimated jet, and the shape of the extended component is similar to those of some Seyfert galaxies. The properties of the emission are compatible with a relic produced by the intermittent activity cycle of the AGN.

Observations at 37 GHz, performed at Mets\"ahovi Radio Observatory, are presented for a sample of 78 radio-loud and radio-quiet narrow-line Seyfert 1 (NLS1) galaxies, together with additional lower and higher frequency radio data from RATAN-600, Owens Valley Radio Observatory, and the Planck satellite. Most of the data have been gathered between February 2012 and April 2015 but for some sources even longer lightcurves exist. The detection rate at 37 GHz is around 19%, comparable to other populations of active galactic nuclei presumed to be faint at radio frequencies, such as BL Lac objects. Variability and spectral indices are determined for sources with enough detections. Based on the radio data, many NLS1 galaxies show a blazar-like radio spectra exhibiting significant variability. The spectra at a given time are often inverted or convex. The source of the high-frequency radio emission in NLS1 galaxies, detected at 37 GHz, is most probably a relativistic jet rather than star formation. Jets in NLS1 galaxies are therefore expected to be a much more common phenomenon than earlier assumed.

Blazars are known to show periods of quiescence followed by outbursts visible throughout the electromagnetic spectrum. We present a novel maximum likelihood approach to capture this bimodal behavior by examining blazar radio variability in the flux-density domain. We separate quiescent and flaring components of a source's light curve by modeling its flux-density distribution as a series of "off" and "on" states. Our modeling allows us to extract information regarding the flaring ratio, duty cycle, and the modulation index in the "off"-state, in the "on"-state, as well as throughout the monitoring period of each blazar. We apply our method to a flux-density-limited subsample from the Owens Valley Radio observatory's 15 GHz blazar monitoring program, and explore differences in the variability characteristics between BL Lacs and FSRQs as well as between $\gamma$-ray detected and non-detected sources. We find that: (1) BL Lacs are more variable and have relatively larger outbursts than the FSRQs, (2) unclassified blazar candidates in our sample show similar variability characteristics as the FSRQs, and (3) $\gamma$-ray detected differ from the $\gamma$-ray non-detected sources in all their variability properties, suggesting a link between the production of $\gamma$-rays and the mechanism responsible for the radio variability. Finally, we fit distributions for blazar flaring ratios, duty cycles, and on- and off- modulation indices that can be used in population studies of variability-dependent blazar properties.

The radio light-curve of J1415+1320 (PKS 1413+135) shows time-symmetric and recurring U-shaped features across the cm-wave and mm-wave bands. The symmetry of these features points to lensing by an intervening object as the cause. U-shaped events in radio light curves in the cm-wave band have previously been attributed to Extreme scattering events (ESE). ESEs are thought to be the result of lensing by compact plasma structures in the Galactic interstellar medium, but the precise nature of these plasma structures remains unknown. Since the strength of a plasma lens evolves with wavelength $\lambda$ as $\lambda^2$, the presence of correlated variations at over a wide wavelength range casts doubt on the canonical ESE interpretation for J1415+1320. In this paper, we critically examine the evidence for plasma lensing in J1415+1320. We compute limits on the lensing strength, and the associated free-free opacity of the putative plasma lenses. We compare the observed and model ESE light curves, and also derive a lower limit on the lens distance based on the effects of parallax due to the Earth's orbit around the Sun. We conclude that plasma lensing is not a viable interpretation for J1415+1320's light curves and that symmetric U-shaped features in the radio light curves of extragalactic sources do not present {\em prima facie} evidence for ESEs. The methodology presented here is generic enough to be applicable to any plasma lensing candidate.

We report the discovery of a rare new form of long-term radio variability in the light-curves of active galaxies (AG) --- Symmetric Achromatic Variability (SAV) --- a pair of opposed and strongly skewed peaks in the radio flux density observed over a broad frequency range. We propose that SAV arises through gravitational milli-lensing when relativistically moving features in AG jets move through gravitational lensing caustics created by $10^3-10^6 \;{\rm M}_{\odot}$ subhalo condensates or black holes located within intervening galaxies. The lower end of this mass range has been inaccessible with previous gravitational lensing techniques. This new interpretation of some AG variability can easily be tested and if it passes these tests, will enable a new and powerful probe of cosmological matter distribution on these intermediate mass scales, as well as provide, for the first time, micro-arcsecond resolution of the nuclei of AG --- a factor of 30--100 greater resolution than is possible with ground-based millimeter VLBI.

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.

The MAGIC (Major Atmospheric Gamma-ray Imaging Cherenkov) telescopes observed the BL Lac object H1722+119 (redshift unknown) for six consecutive nights between 2013 May 17 and 22, for a total of 12.5 h. The observations were triggered by high activity in the optical band measured by the KVA (Kungliga Vetenskapsakademien) telescope. The source was for the first time detected in the very high energy (VHE, $E > 100$ GeV) $\gamma$-ray band with a statistical significance of 5.9 $\sigma$. The integral flux above 150 GeV is estimated to be $(2.0\pm 0.5)$ per cent of the Crab Nebula flux. We used contemporaneous high energy (HE, 100 MeV $ < E < 100$ GeV) $\gamma$-ray observations from Fermi-LAT (Large Area Telescope) to estimate the redshift of the source. Within the framework of the current extragalactic background light models, we estimate the redshift to be $z = 0.34 \pm 0.15$. Additionally, we used contemporaneous X-ray to radio data collected by the instruments on board the Swift satellite, the KVA, and the OVRO (Owens Valley Radio Observatory) telescope to study multifrequency characteristics of the source. We found no significant temporal variability of the flux in the HE and VHE bands. The flux in the optical and radio wavebands, on the other hand, did vary with different patterns. The spectral energy distribution (SED) of H1722+119 shows surprising behaviour in the $\sim 3\times10^{14} - 10^{18}$ Hz frequency range. It can be modelled using an inhomogeneous helical jet synchrotron self-Compton model.

After a long low-activity period, a gamma-ray flare from the narrow-line Seyfert 1 PKS 1502+036 (z=0.4089) was detected by the Large Area Telescope (LAT) on board Fermi in 2015. On 2015 December 20 the source reached a daily peak flux, in the 0.1-300 GeV band, of (93 $\pm$ 19) $\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$, attaining a flux of (237 $\pm$ 71) $\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$ on 3-hr time-scales, which corresponds to an isotropic luminosity of (7.3 $\pm$ 2.1) $\times$10$^{47}$ erg/s. The gamma-ray flare was not accompanied by significant spectral changes. We report on multi-wavelength radio-to-gamma-ray observations of PKS 1502+036 during 2008 August-2016 March by Fermi-LAT, Swift, XMM-Newton, Catalina Real-Time Transient Survey, and the Owens Valley Radio Observatory (OVRO). An increase in activity was observed on 2015 December 22 by Swift in optical, UV, and X-rays. The OVRO 15 GHz light curve reached the highest flux density observed from this source on 2016 January 12, indicating a delay of about three weeks between the gamma-ray and 15 GHz emission peaks. This suggests that the gamma-ray emitting region is located beyond the broad line region. We compared the spectral energy distribution (SED) of an average activity state with that of the flaring state. The two SED, with the high-energy bump modelled as an external Compton component with seed photons from a dust torus, could be fitted by changing the electron distribution parameters as well as the magnetic field. The fit of the disc emission during the average state constrains the black hole mass to values lower than 10$^8$ solar masses. The SED, high-energy emission mechanisms, and gamma-ray properties of the source resemble those of a flat spectrum radio quasar.

To fully exploit the scientific potential of the Fermi mission, we initiated the F-GAMMA program. Between 2007 and 2015 it was the prime provider of complementary multi-frequency monitoring in the radio regime. We quantify the radio variability of gamma-ray blazars. We investigate its dependence on source class and examine whether the radio variability is related to the gamma-ray loudness. Finally, we assess the validity of a putative correlation between the two bands. The F-GAMMA monitored monthly a sample of about 60 sources at up to twelve radio frequencies between 2.64 and 228.39 GHz. We perform a time series analysis on the first 2.5-year dataset to obtain variability parameters. A maximum likelihood analysis is used to assess the significance of a correlation between radio and gamma-ray fluxes. We present light curves and spectra (coherent within ten days) obtained with the Effelsberg 100-m and IRAM 30-m telescopes. All sources are variable across all frequency bands with amplitudes increasing with frequency up to rest frame frequencies of around 60 - 80 GHz as expected by shock-in-jet models. Compared to FSRQs, BL Lacs show systematically lower variability amplitudes, brightness temperatures and Doppler factors at lower frequencies, while the difference vanishes towards higher ones. The time scales appear similar for the two classes. The distribution of spectral indices appears flatter or more inverted at higher frequencies for BL Lacs. Evolving synchrotron self-absorbed components can naturally account for the observed spectral variability. We find that the Fermi-detected sources show larger variability amplitudes as well as brightness temperatures and Doppler factors, than non-detected ones. Flux densities at 86.2 and 142.3 GHz correlate with 1 GeV fluxes at a significance level better than 3sigma, implying that gamma rays are produced very close to the mm-band emission region.

We report radio SETI observations on a large number of known exoplanets and other nearby star systems using the Allen Telescope Array (ATA). Observations were made over about 19000 hours from May 2009 to Dec 2015. This search focused on narrow-band radio signals from a set totaling 9293 stars, including 2015 exoplanet stars and Kepler objects of interest and an additional 65 whose planets may be close to their Habitable Zone. The ATA observations were made using multiple synthesized beams and an anticoincidence filter to help identify terrestrial radio interference. Stars were observed over frequencies from 1- 9 GHz in multiple bands that avoid strong terrestrial communication frequencies. Data were processed in near-real time for narrow-band (0.7- 100 Hz) continuous and pulsed signals, with transmitter/receiver relative accelerations from -0.3 to 0.3 m/s^2. A total of 1.9 x 10^8 unique signals requiring immediate follow-up were detected in observations covering more than 8 x 10^6 star-MHz. We detected no persistent signals from extraterrestrial technology exceeding our frequency-dependent sensitivity threshold of 180 - 310 x 10^-26 W / m^2.

We compare the variability properties of very high energy gamma-ray emitting BL Lac objects in the optical and radio bands. We use the variability information to distinguish multiple emission components in the jet, to be used as a guidance for spectral energy distribution modelling. Our sample includes 32 objects in the Northern sky that have data for at least 2 years in both bands. We use optical R-band data from the Tuorla blazar monitoring program and 15 GHz radio data from the Owens Valley Radio Observatory blazer monitoring program. We estimate the variability amplitudes using the intrinsic modulation index, and study the time-domain connection by cross-correlating the optical and radio light curves assuming power law power spectral density. Our sample objects are in general more variable in the optical than radio. We find correlated flares in about half of the objects, and correlated long-term trends in more than 40% of the objects. In these objects we estimate that at least 10%-50% of the optical emission originates in the same emission region as the radio, while the other half is due to faster variations not seen in the radio. This implies that simple single-zone spectral energy distribution models are not adequate for many of these objects.

1ES 1011+496 $(z=0.212)$ was discovered in very high energy (VHE, E >100 GeV) $\gamma$-rays with MAGIC in 2007. The absence of simultaneous data at lower energies led to a rather incomplete characterization of the broadband spectral energy distribution (SED). We study the source properties and the emission mechanisms, probing whether a simple one-zone synchrotron-self-Compton (SSC) scenario is able to explain the observed broadband spectrum. We analyzed VHE to radio data from 2011 and 2012 collected by MAGIC, $Fermi$-LAT, $Swift$, KVA, OVRO, and Mets\"ahovi in addition to optical polarimetry data and radio maps from the Liverpool Telescope and MOJAVE. The VHE spectrum was fit with a simple power law with a photon index of $3.69\pm0.22$ and a flux above 150 GeV of $(1.46\pm0.16)\times10^{-11}$ ph cm$^{-2}$ s$^{-1}$. 1ES 1011+496 was found to be in a generally quiescent state at all observed wavelengths, showing only moderate variability from radio to X-rays. A low degree of polarization of less than 10% was measured in optical, while some bright features polarized up to 60% were observed in the radio jet. A similar trend in the rotation of the electric vector position angle was found in optical and radio. The radio maps indicated a superluminal motion of $1.8\pm0.4\,c$, which is the highest speed statistically significantly measured so far in a high-frequency-peaked BL Lac. For the first time, the high-energy bump in the broadband SED of 1ES 1011+496 could be fully characterized from 0.1 GeV to 1 TeV which permitted a more reliable interpretation within the one-zone SSC scenario. The polarimetry data suggest that at least part of the optical emission has its origin in some of the bright radio features, while the low polarization in optical might be due to the contribution of parts of the radio jet with different orientations of the magnetic field to the optical emission.

We report on a search for the presence of signals from extraterrestrial intelligence in the direction of the star system KIC 8462852. Observations were made at radio frequencies between 1-10 GHz using the Allen Telescope Array. No narrowband radio signals were found at a level of 180-300 Jy in a 1 Hz channel, or medium band signals above 10 Jy in a 100 kHz channel.

We present 1625 new 15 GHz (2 cm) VLBA images of 295 jets associated with active galactic nuclei (AGNs) from the MOJAVE and 2 cm VLBA surveys, spanning observations between 1994 Aug 31 and 2013 Aug 20. For 274 AGNs with at least 5 VLBA epochs, we have analyzed the kinematics of 961 individual bright features in their parsec-scale jets. A total of 122 of these jets have not been previously analyzed by the MOJAVE program. In the case of 451 jet features that had at least 10 epochs, we also examined their kinematics for possible accelerations. At least half of the well-sampled features have non-radial and/or accelerating trajectories, indicating that non-ballistic motion is common in AGN jets. Since it is impossible to extrapolate any accelerations that occurred before our monitoring period, we could only determine reliable ejection dates for about 24% of those features that had significant proper motions. The distribution of maximum apparent jet speeds in all 295 AGNs measured by our program to date is peaked below 5c, with very few jets with apparent speeds above 30c. The fastest speed in our survey is about 50c, measured in the jet of the quasar PKS 0805-07, and is indicative of a maximum jet Lorentz factor of about 50 in the parent population. The Fermi LAT-detected gamma-ray AGNs in our sample have, on average, higher jet speeds than non LAT-detected AGNs, indicating a strong correlation between pc-scale jet speed and gamma-ray Doppler boosting factor. We have identified 11 moderate-redshift (z<0.35) AGNs with fast apparent speeds (>10c) that are strong candidates for future TeV gamma-ray detection. Of the five gamma-ray loud narrow-lined Seyfert I AGNs in our sample, three show highly superluminal jet motions, while the others have sub-luminal speeds. (abridged)

Narrow-line Seyfert 1 galaxies (NLS1s) are active galactic nuclei (AGN) recently identified as a new class of $\gamma$-ray sources. The high energy emission is explained by the presence of a relativistic jet observed at small angles, just like in the case of blazars. When the latter are observed at larger angles they appear as radio-galaxies, but an analogue parent population for beamed NLS1s has not yet been determined. In this work we analyze this problem by studying the physical properties of three different samples of parent sources candidates: steep-spectrum radio-loud NLS1s, radio-quiet NLS1s, and disk-hosted radio-galaxies, along with compact steep-spectrum sources. In our approach, we first derived black hole mass and Eddington ratio from the optical spectra, then we investigated the interaction between the jet and the narrow-line region from the [O III] $\lambda\lambda$4959,5007 lines. Finally, the radio luminosity function allowed us to compare their jet luminosity and hence determine the relations between the samples.

We report about recent updates of broad-band properties of radio-loud narrow-line Seyfert 1 galaxies.

Narrow-line Seyfert 1 galaxies (NLS1) are active galactic nuclei (AGN) known to have small masses of the central black hole and high accretion rates. NLS1s are generally radio-quiet, but a small part of them (about 7\%) are radio-loud. The recent discovery of powerful relativistic jets in radio-loud NLS1s (RLNLS1s), emitting at high-energy $\gamma$-rays, opened intriguing questions. The observed luminosity of the jet is generally weak, smaller than blazars, although when rescaled for the mass of the central black hole, it becomes of the same order of magnitude of the latter. The weak luminosity, and hence observed flux, resulted in a small number of known RLNLS1. From a recent survey of RLNLS1s, it was found that only 8 out of 42 sources had radio flux density at 1.4 GHz greater than 100 mJy, while 21 out of 42 had flux density smaller than 10 mJy. In addition, given the strong variability at all wavelengths, with present-day facilities RLNLS1s can often only be detected during high activity periods. The Square Kilometer Array (SKA), with its superior sensitivity, will break this limit, allowing us to unveil a relatively unknown population of jetted AGN. We present the results of a study aimed at evaluating the scenario that could emerge after the advent of SKA.

We present results on multifrequency Very Long Baseline Array (VLBA) monitoring observations of the double-image gravitationally lensed blazar JVAS B0218+357. Multi-epoch observations started less than one month after the gamma-ray flare detected in 2012 by the Large Area Telescope on board Fermi, and spanned a 2-month interval. The radio light curves did not reveal any significant flux density variability, suggesting that no clear correlation between the high energy and low-energy emission is present. This behaviour was confirmed also by the long-term Owens Valley Radio Observatory monitoring data at 15 GHz. The milliarcsecond-scale resolution provided by the VLBA observations allowed us to resolve the two images of the lensed blazar, which have a core-jet structure. No significant morphological variation is found by the analysis of the multi-epoch data, suggesting that the region responsible for the gamma-ray variability is located in the core of the AGN, which is opaque up to the highest observing frequency of 22 GHz.

Outbursts from gamma-ray quasars provide insights on the relativistic jets of active galactic nuclei and constraints on the diffuse radiation fields that fill the Universe. The detection of significant emission above 100 GeV from a distant quasar would show that some of the radiated gamma rays escape pair-production interactions with low-energy photons, be it the extragalactic background light (EBL), or the radiation near the supermassive black hole lying at the jet's base. VERITAS detected gamma-ray emission up to 200 GeV from PKS 1441+25 (z=0.939) during April 2015, a period of high activity across all wavelengths. This observation of PKS 1441+25 suggests that the emission region is located thousands of Schwarzschild radii away from the black hole. The gamma-ray detection also sets a stringent upper limit on the near-ultraviolet to near-infrared EBL intensity, suggesting that galaxy surveys have resolved most, if not all, of the sources of the EBL at these wavelengths.

To explore the hypothesis that KIC 8462852's aperiodic dimming is caused by artificial megastructures in orbit (Wright et al. 2015), rather than a natural cause such as cometary fragments in a highly elliptical orbit (Marengo et al. 2015), we searched for electromagnetic signals from KIC 8462852 indicative of extraterrestrial intelligence. The primary observations were in the visible optical regime using the Boquete Optical SETI Observatory in Panama. In addition, as a preparatory exercise for the possible future detection of a candidate signal (Heidmann 1991), three of six observing runs simultaneously searched radio frequencies at the Allen Telescope Array in California. No periodic optical signals greater than 67 photons/m2 within a time frame of 25 ns were seen. This limit corresponds to isotropic optical pulses of 8E22 joules. If, however, any inhabitants of KIC 8462852 were targeting our solar system (Shostak & Villard 2004), the required energy would be reduced greatly. The limits on narrowband radio signals were 180 - 300 Jy Hz at 1 and 8 GHz, respectively, corresponding to a transmitter with an effective isotropic radiated power of 4E15 W (and 7E15 W) at the distance of KIC 8462852. While these powers requirements are high, even modest targeting could - just as for optical signals - lower these numbers substantially.

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 invariance of the speed of light implies a series of consequences related to our perception of simultaneity and of time itself. Whilst these consequences are experimentally well studied for subluminal speeds, the kinematics of superluminal motion lack direct evidence. Using high temporal resolution imaging techniques, we demonstrate that if a source approaches an observer at superluminal speeds, the temporal ordering of events is inverted and its image appears to propagate backwards. If the source changes its speed, crossing the interface between sub- and super-luminal propagation, we observe image pair annihilation and creation. These results show that it is not possible to unambiguously determine the kinematics of an event from imaging and time-resolved measurements alone.

The aim of this work is to understand the nature of the parent population of beamed narrow-line Seyfert 1 galaxies (NLS1s), by studying the physical properties of three parent candidates samples: steep-spectrum radio-loud NLS1s, radio-quiet NLS1s and disk-hosted radio-galaxies. In particular, we focused on the black hole mass and Eddington ratio distribution and on the interactions between the jet and the narrow-line region.

We report on a multiwavelength survey of a sample of 42 flat-spectrum radio-loud narrow-line Seyfert 1 galaxies (RLNLS1s). This is the largest known sample of this type of active galactic nucleus (AGN) to date. We found that 17% of sources were detected at high-energy gamma rays (E>100 MeV), and 90% at X-rays (0.3-10 keV). The masses of the central black holes are in the range $\sim 10^{6-8}M_{\odot}$, smaller than the values of blazars. The disk luminosities are about 1-49% of the Eddington value, with one outlier at 0.3%, comparable with the luminosities observed in flat-spectrum radio quasars (FSRQs). The jet powers are $\sim 10^{42-46}$ erg s$^{-1}$, comparable with BL Lac Objects, yet relatively smaller than FSRQs. However, once renormalized by the mass of the central black hole, the jet powers of RLNLS1s, BL Lacs, and FSRQs are consistent each other, indicating the scalability of the jets. We found episodes of extreme variability at high energies on time scales of hours. In some cases, dramatic spectral and flux changes are interpreted as the interplay between the relativistic jet and the accretion disk. We conclude that, despite the distinct observational properties, the central engines of RLNLS1s are similar to those of blazars.

We use the complete MOJAVE 1.5 Jy sample of active galactic nuclei (AGN) to examine the gamma-ray detection statistics of the brightest radio-loud blazars in the northern sky. We find that 23% of these AGN were not detected above 0.1 GeV by the Fermi LAT during the 4-year 3FGL catalog period partly because of an instrumental selection effect, and partly due to their lower Doppler boosting factors. Blazars with synchrotron peaks in their spectral energy distributions located below $10^{13.4}$ Hz also tend to have high-energy peaks that lie below the 0.1 GeV threshold of the LAT, and are thus less likely to be detected by Fermi. The non-detected AGN in the 1.5 Jy sample also have significantly lower 15 GHz radio modulation indices and apparent jet speeds, indicating that they have lower than average Doppler factors. Since the effective amount of relativistic Doppler boosting is enhanced in gamma-rays (particularly in the case of external inverse-Compton scattering), this makes them less likely to appear in the 3FGL catalog. Based on their observed properties, we have identified several bright radio-selected blazars that are strong candidates for future detection by Fermi.

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.

The high-frequency-peaked BL Lac (HBL) 1ES 0806+524 (z = 0.138) was discovered in VHE $\gamma$ rays in 2008. Until now, the broad-band spectrum of 1ES 0806+524 has been only poorly characterized, in particular at high energies. We analysed multiwavelength observations from $\gamma$ rays to radio performed from 2011 January to March, which were triggered by the high activity detected at optical frequencies. These observations constitute the most precise determination of the broad-band emission of 1ES 0806+524 to date. The stereoscopic MAGIC observations yielded a $\gamma$-ray signal above 250 GeV of $(3.7 \pm 0.7)$ per cent of the Crab Nebula flux with a statistical significance of 9.9 $\sigma$. The multiwavelength observations showed significant variability in essentially all energy bands, including a VHE $\gamma$-ray flare that lasted less than one night, which provided unprecedented evidence for short-term variability in 1ES 0806+524. The spectrum of this flare is well described by a power law with a photon index of $2.97 \pm 0.29$ between $\sim$150 GeV and 1 TeV and an integral flux of $(9.3 \pm 1.9)$ per cent of the Crab Nebula flux above 250 GeV. The spectrum during the non-flaring VHE activity is compatible with the only available VHE observation performed in 2008 with VERITAS when the source was in a low optical state. The broad-band spectral energy distribution can be described with a one-zone Synchrotron Self Compton model with parameters typical for HBLs, indicating that 1ES 0806+524 is not substantially different from the HBLs previously detected.

Flat-spectrum radio-loud Narrow-Line Seyfert 1 galaxies (NLS1s) are a recently discovered class of $\gamma$-ray emitting Active Galactic Nuclei (AGN), that exhibit some blazar-like properties which are explained with the presence of a relativistic jet viewed at small angles. When blazars are observed at larger angles they appear as radio-galaxies, and we expect to observe an analogue parent population for beamed NLS1s. However, the number of known NLS1s with the jet viewed at large angles is not enough. Therefore, we tried to understand the origin of this deficit. Current hypotheses about the nature of parent sources are steep-spectrum radio-loud NLS1s, radio-quiet NLS1s and disk-hosted radio-galaxies. To test these hypotheses we built three samples of candidate sources plus a control sample, and calculated their black hole mass and Eddington ratio using their optical spectra. We then performed a Kolmogorov-Smirnov statistical test to investigate the compatibility of our different samples with a beamed population. Our results indicate that, when the inclination angle increases, a beamed source appears as a steep-spectrum radio-loud NLS1, or possibly even as a disk-hosted radio-galaxy with low black hole mass and high Eddington ratio. Further investigations, involving larger complete samples and observations at radio frequency, are needed to understand the incidence of disk-hosted radio-galaxies in the parent population, and to assess whether radio-quiet NLS1s can play a role, as well.

PG 1553+113 is a very-high-energy (VHE, $E>100\,\mathrm{GeV}$) $\gamma$-ray emitter classified as a BL Lac object. Its redshift is constrained by intergalactic absorption lines in the range $0.4<z<0.58$. The MAGIC telescopes have monitored the source's activity since 2005. In early 2012, PG 1553+113 was found in a high-state, and later, in April of the same year, the source reached its highest VHE flux state detected so far. Simultaneous observations carried out in X-rays during 2012 April show similar flaring behaviour. In contrast, the $\gamma$-ray flux at $E<100\,\mathrm{GeV}$ observed by Fermi-LAT is compatible with steady emission. In this paper, a detailed study of the flaring state is presented. The VHE spectrum shows clear curvature, being well fitted either by a power law with an exponential cut-off or by a log-parabola. A simple power-law fit hypothesis for the observed shape of the PG 1553+113 VHE $\gamma$-ray spectrum is rejected with a high significance (fit probability P=2.6 $\times 10^{-6}$). The observed curvature is compatible with the extragalactic background light (EBL) imprint predicted by current generation EBL models assuming a redshift $z\sim0.4$. New constraints on the redshift are derived from the VHE spectrum. These constraints are compatible with previous limits and suggest that the source is most likely located around the optical lower limit, $z=0.4$, based on the detection of Ly$\alpha$ absorption. Finally, we find that the synchrotron self-Compton (SSC) model gives a satisfactory description of the observed multi-wavelength spectral energy distribution during the flare.

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.

In 2012 Markarian 421 underwent the largest flare ever observed in this blazar at radio frequencies. In the present study, we start exploring this unique event and compare it to a less extreme event in 2013. We use 15 GHz radio data obtained with the Owens Valley Radio Observatory 40-m telescope, 95 GHz millimeter data from the Combined Array for Research in Millimeter-Wave Astronomy, and GeV gamma-ray data from the Fermi Gamma-ray Space Telescope. The radio light curves during the flaring periods in 2012 and 2013 have very different appearances, both in shape and peak flux density. Assuming that the radio and gamma-ray flares are physically connected, we attempt to model the most prominent sub-flares of the 2012 and 2013 activity periods by using the simplest possible theoretical framework. We first fit a one-zone synchrotron self-Compton (SSC) model to the less extreme 2013 flare and estimate parameters describing the emission region. We then model the major gamma-ray and radio flares of 2012 using the same framework. The 2012 gamma-ray flare shows two distinct spikes of similar amplitude, so we examine scenarios associating the radio flare with each spike in turn. In the first scenario, we cannot explain the sharp radio flare with a simple SSC model, but we can accommodate this by adding plausible time variations to the Doppler beaming factor. In the second scenario, a varying Doppler factor is not needed, but the SSC model parameters require fine tuning. Both alternatives indicate that the sharp radio flare, if physically connected to the preceding gamma-ray flares, can be reproduced only for a very specific choice of parameters.

We study the kinematics of ridge lines on the pc-scale jet of the active galactic nucleus BL Lac. We show that the ridge lines display transverse patterns that move superluminally downstream, and that the moving patterns are analogous to waves on a whip. Their apparent speeds $\beta_\mathrm{app}$ (units of $c$) range from 3.9 to 13.5, corresponding to $\beta_\mathrm{wave}^\mathrm{gal}= 0.981 - 0.998$ in the galaxy frame. We show that the magnetic field in the jet is well-ordered with a strong transverse component, and assume that it is helical and that the transverse patterns are Alfv\'en waves propagating downstream on the longitudinal component of the magnetic field. The wave-induced transverse speed of the jet is non-relativistic ($\beta_\mathrm{tr}^\mathrm{gal} \lesssim 0.09$). In 2010 the wave activity subsided and the jet then displayed a mild wiggle that had a complex oscillatory behaviour. The Alfv\'en waves appear to be excited by changes in the position angle of the recollimation shock, in analogy to exciting a wave on a whip by shaking the handle. A simple model of the system with plasma sound speed $\beta_\mathrm{s}=0.3$ and apparent speed of a slow MHD wave $\beta_\mathrm{app,S}=4$ yields Lorentz factor of the beam $\Gamma_\mathrm{beam} \sim 4.5$, pitch angle of the helix (in the beam frame) $\alpha\sim 67^\circ$, Alfv\'en speed $\beta_\mathrm{A}\sim 0.64$, and magnetosonic Mach number $M_\mathrm{ms}\sim 4.7$. This describes a plasma in which the magnetic field is dominant and in a rather tight helix, and Alfv\'en waves are responsible for the moving transverse patterns.

We have discovered kiloparsec-scale extended radio emission in three narrow-line Seyfert 1 galaxies (NLS1s) in sub-arcsecond resolution 9 GHz images from the Karl G. Jansky Very Large Array (VLA). We find all sources show two-sided, mildly core-dominated jet structures with diffuse lobes dominated by termination hotspots. These span 20-70 kpc with morphologies reminiscent of FR II radio galaxies, while the extended radio luminosities are intermediate between FR I and FR II sources. In two cases the structure is linear, while a $45^{\circ}$ bend is apparent in the third. Very Long Baseline Array images at 7.6 GHz reveal parsec-scale jet structures, in two cases with extended structure aligned with the inner regions of the kiloparsec-scale jets. Based on this alignment, the ratio of the radio core luminosity to the optical luminosity, the jet/counter-jet intensity and extension length ratios, and moderate core brightness temperatures ($\lesssim10^{10}$ K), we conclude these jets are mildly relativistic ($\beta\lesssim0.3$, $\delta\sim1$-$1.5$) and aligned at moderately small angles to the line of sight (10-15$^{\circ}$). The derived kinematic ages of $\sim10^6$-$10^7$ y are much younger than radio galaxies but comparable to other NLS1s. Our results increase the number of radio-loud NLS1s with known kiloparsec-scale extensions from seven to ten and suggest that such extended emission may be common, at least among the brightest of these sources.

To study the effect of supermassive black holes (SMBHs) on their host galaxies it is important to study the hosts when the SMBH is near its peak activity. A method to investigate the host galaxies of high luminosity quasars is to obtain optical spectra at positions offset from the nucleus where the relative contribution of the quasar and host are comparable. However, at these extended radii the galaxy surface brightness is often low (20-22 mag per arcsec$^{2}$) and the resulting spectrum might have such low S/N that it hinders analysis with standard stellar population modeling techniques. To address this problem we have developed a method that can recover galaxy star formation histories (SFHs) from rest frame optical spectra with S/N $\sim$ 5~\AA$^{-1}$. This method uses the statistical technique diffusion k-means to tailor the stellar population modeling basis set. Our diffusion k-means minimal basis set, composed of 4 broad age bins, is successful in recovering a range of galaxy SFHs. Additionally, using an analytic prescription for seeing conditions, we are able to simultaneously model scattered quasar light and the SFH of quasar host galaxies (QHGs). We use synthetic data to compare results of our novel method with previous techniques. We also present the modeling results on a previously published QHG and show that galaxy properties recovered from a diffusion k-means basis set are less sensitive to noise added to this quasar host galaxy spectrum. Our new method has a clear advantage in recovering information from QHGs and could also be applied to the analysis of other low S/N galaxy spectra such as those typically obtained for high redshift objects or integral field spectroscopic surveys.

Several narrow-line Seyfert 1 galaxies (NLS1s) have now been detected in gamma rays, providing firm evidence that at least some of this class of active galactic nuclei (AGN) produce relativistic jets. The presence of jets in NLS1s is surprising, as these sources are typified by comparatively small black hole masses and near- or super-Eddington accretion rates. This challenges the current understanding of the conditions necessary for jet production. Comparing the properties of the jets in NLS1s with those in more familiar jetted systems is thus essential to improve jet production models. We present early results from our campaign to monitor the kinematics and polarization of the parsec-scale jets in a sample of 15 NLS1s through multifrequency observations with the Very Long Baseline Array. These observations are complemented by fast-cadence 15 GHz monitoring with the Owens Valley Radio Observatory 40m telescope and optical spectroscopic monitoring with with the 2m class telescope at the Guillermo Haro Astrophysics Observatory in Cananea, Mexico.

A fundamental challenge for wide-field imaging surveys is obtaining follow-up spectroscopic observations: there are > $10^9$ photometrically cataloged sources, yet modern spectroscopic surveys are limited to ~few x $10^6$ targets. As we approach the Large Synoptic Survey Telescope (LSST) era, new algorithmic solutions are required to cope with the data deluge. Here we report the development of a machine-learning framework capable of inferring fundamental stellar parameters (Teff, log g, and [Fe/H]) using photometric-brightness variations and color alone. A training set is constructed from a systematic spectroscopic survey of variables with Hectospec/MMT. In sum, the training set includes ~9000 spectra, for which stellar parameters are measured using the SEGUE Stellar Parameters Pipeline (SSPP). We employed the random forest algorithm to perform a non-parametric regression that predicts Teff, log g, and [Fe/H] from photometric time-domain observations. Our final, optimized model produces a cross-validated root-mean-square error (RMSE) of 165 K, 0.39 dex, and 0.33 dex for Teff, log g, and [Fe/H], respectively. Examining the subset of sources for which the SSPP measurements are most reliable, the RMSE reduces to 125 K, 0.37 dex, and 0.27 dex, respectively, comparable to what is achievable via low-resolution spectroscopy. For variable stars this represents a ~12-20% improvement in RMSE relative to models trained with single-epoch photometric colors. As an application of our method, we estimate stellar parameters for ~54,000 known variables. We argue that this method may convert photometric time-domain surveys into pseudo-spectrographic engines, enabling the construction of extremely detailed maps of the Milky Way, its structure, and history.

Since mid-2007 we have carried out a dedicated long-term monitoring programme at 15 GHz using the Owens Valley Radio Observatory 40 meter telescope. One of the main goals of this programme is to study the relation between the radio and gamma-ray emission in blazars and to use it as a tool to locate the site of high energy emission. Using this large sample of objects we are able to characterize the radio variability, and study the significance of correlations between the radio and gamma-ray bands. We find that the radio variability of many sources can be described using a simple power law power spectral density, and that when taking into account the red-noise characteristics of the light curves, cases with significant correlation are rare. We note that while significant correlations are found in few individual objects, radio variations are most often delayed with respect to the gamma-ray variations. This suggests that the gamma-ray emission originates upstream of the radio emission. Because strong flares in most known gamma-ray-loud blazars are infrequent, longer light curves are required to settle the issue of the strength of radio-gamma cross-correlations and establish confidently possible delays between the two. For this reason continuous multiwavelength monitoring over a longer time period is essential for statistical tests of jet emission models.

We report on the acceleration properties of 329 features in 95 blazar jets from the MOJAVE VLBA program. Nearly half the features and three-quarters of the jets show significant changes in speed and/or direction. In general, apparent speed changes are distinctly larger than changes in direction, indicating that changes in the Lorentz factors of jet features dominate the observed speed changes rather than bends along the line of sight. Observed accelerations tend to increase the speed of features near the jet base, $\lesssim 10-20$ parsecs projected, and decrease their speed at longer distances. The range of apparent speeds at fixed distance in an individual jet can span a factor of a few, indicating that shock properties and geometry may influence the apparent motions; however, we suggest that the broad trend of jet features increasing their speed near the origin is due to an overall acceleration of the jet flow out to de-projected distances of order $10^2$ parsecs, beyond which the flow begins to decelerate or remains nearly constant in speed. We estimate intrinsic rates of change of the Lorentz factors in the galaxy frame of order $\dot{\Gamma}/\Gamma \simeq 10^{-3}$ to $10^{-2}$ per year which can lead to total Lorentz factor changes of a factor of a few on the length scales observed here. Finally, we also find evidence for jet collimation at projected distances of $\lesssim 10$ parsecs in the form of the non-radial motion and bending accelerations that tend to better align features with the inner jet.

We report on multifrequency observations performed during 2012 December-2013 August of the first narrow-line Seyfert 1 galaxy detected in gamma rays, PMN J0948+0022 ($z$ = 0.5846). A gamma-ray flare was observed by the Large Area Telescope on board Fermi during 2012 December-2013 January, reaching a daily peak flux in the 0.1-100 GeV energy range of (155 $\pm$ 31) $\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$ on 2013 January 1, corresponding to an apparent isotropic luminosity of about 1.5$\times$10$^{48}$ erg s$^{-1}$. The gamma-ray flaring period triggered Swift and VERITAS observations in addition to radio and optical monitoring by OVRO, MOJAVE, and CRTS. A strong flare was observed in optical, UV, and X-rays on 2012 December 30, quasi-simultaneously to the gamma-ray flare, reaching a record flux for this source from optical to gamma rays. VERITAS observations at very high energy (E > 100 GeV) during 2013 January 6-17 resulted in an upper limit of F (> 0.2 TeV) < 4.0$\times$10$^{-12}$ ph cm$^{-2}$ s$^{-1}$. We compared the spectral energy distribution (SED) of the flaring state in 2013 January with that of an intermediate state observed in 2011. The two SEDs, modelled as synchrotron emission and an external Compton scattering of seed photons from a dust torus, can be modelled by changing both the electron distribution parameters and the magnetic field.

We aim to characterize the broadband emission from 2FGL J2001.1+4352, which has been associated with the unknown-redshift blazar MG4 J200112+4352. Based on its gamma-ray spectral properties, it was identified as a potential very high energy (VHE; E > 100 GeV) gamma-ray emitter. The source was observed with MAGIC first in 2009 and later in 2010 within a multi-instrument observation campaign. The MAGIC observations yielded 14.8 hours of good quality stereoscopic data. The object was monitored at radio, optical and gamma-ray energies during the years 2010 and 2011. The source, named MAGIC J2001+439, is detected for the first time at VHE with MAGIC at a statistical significance of 6.3 {\sigma} (E > 70 GeV) during a 1.3-hour long observation on 2010 July 16. The multi-instrument observations show variability in all energy bands with the highest amplitude of variability in the X-ray and VHE bands. We also organized deep imaging optical observations with the Nordic Optical Telescope in 2013 to determine the source redshift. We determine for the first time the redshift of this BL Lac object through the measurement of its host galaxy during low blazar activity. Using the observational evidence that the luminosities of BL Lac host galaxies are confined to a relatively narrow range, we obtain z = 0.18 +/- 0.04. Additionally, we use the Fermi-LAT and MAGIC gamma-ray spectra to provide an independent redshift estimation, z = 0.17 +/- 0.10. Using the former (more accurate) redshift value, we adequately describe the broadband emission with a one-zone SSC model for different activity states and interpret the few-day timescale variability as produced by changes in the high-energy component of the electron energy distribution.

We present a practical implementation of a Monte Carlo method to estimate the significance of cross-correlations in unevenly sampled time series of data, whose statistical properties are modeled with a simple power-law power spectral density. This implementation builds on published methods, we introduce a number of improvements in the normalization of the cross-correlation function estimate and a bootstrap method for estimating the significance of the cross-correlations. A closely related matter is the estimation of a model for the light curves, which is critical for the significance estimates. We present a graphical and quantitative demonstration that uses simulations to show how common it is to get high cross-correlations for unrelated light curves with steep power spectral densities. This demonstration highlights the dangers of interpreting them as signs of a physical connection. We show that by using interpolation and the Hanning sampling window function we are able to reduce the effects of red-noise leakage and to recover steep simple power-law power spectral densities. We also introduce the use of a Neyman construction for the estimation of the errors in the power-law index of the power spectral density. This method provides a consistent way to estimate the significance of cross-correlations in unevenly sampled time series of data.

In order to determine the location of the gamma-ray emission site in blazars, we investigate the time-domain relationship between their radio and gamma-ray emission. Light-curves for the brightest detected blazars from the first 3 years of the mission of the Fermi Gamma-ray Space Telescope are cross-correlated with 4 years of 15GHz observations from the OVRO 40-m monitoring program. The large sample and long light-curve duration enable us to carry out a statistically robust analysis of the significance of the cross-correlations, which is investigated using Monte Carlo simulations including the uneven sampling and noise properties of the light-curves. Modeling the light-curves as red noise processes with power-law power spectral densities, we find that only one of 41 sources with high quality data in both bands shows correlations with significance larger than 3-sigma (AO 0235+164), with only two more larger than even 2.25-sigma (PKS 1502+106 and B2 2308+34). Additionally, we find correlated variability in Mrk 421 when including a strong flare that occurred in July-September 2012. These results demonstrate very clearly the difficulty of measuring statistically robust multiwavelength correlations and the care needed when comparing light-curves even when many years of data are used. This should be a caution. In all four sources the radio variations lag the gamma-ray variations, suggesting that the gamma-ray emission originates upstream of the radio emission. Continuous simultaneous monitoring over a longer time period is required to obtain high significance levels in cross-correlations between gamma-ray and radio variability in most blazars.

We study the multifrequency emission and spectral properties of the quasar 3C 279. We observed 3C 279 in very high energy (VHE, E>100GeV) gamma rays, with the MAGIC telescopes during 2011, for the first time in stereoscopic mode. We combine these measurements with observations at other energy bands: in high energy (HE, E>100MeV) gamma rays from Fermi-LAT, in X-rays from RXTE, in the optical from the KVA telescope and in the radio at 43GHz, 37GHz and 15GHz from the VLBA, Mets\"ahovi and OVRO radio telescopes and optical polarisation measurements from the KVA and Liverpool telescopes. During the MAGIC observations (February to April 2011) 3C 279 was in a low state in optical, X-ray and gamma rays. The MAGIC observations did not yield a significant detection. These upper limits are in agreement with the extrapolation of the HE gamma-ray spectrum, corrected for extragalactic background light absorption, from Fermi-LAT. The second part of the MAGIC observations in 2011 was triggered by a high activity state in the optical and gamma-ray bands. During the optical outburst the optical electric vector position angle rotatated of about 180 degrees. There was no simultaneous rotation of the 43GHz radio polarisation angle. No VHE gamma rays were detected by MAGIC, and the derived upper limits suggest the presence of a spectral break or curvature between the Fermi-LAT and MAGIC bands. The combined upper limits are the strongest derived to date for the source at VHE and below the level of the previously detected flux by a factor 2. Radiation models that include synchrotron and inverse Compton emissions match the optical to gamma-ray data, assuming an emission component inside the broad line region (BLR) responsible for the high-energy emission and one outside the BLR and the infrared torus causing optical and low-energy emission. We interpreted the optical polarisation with a bent trajectory model.

Among more than fifty blazars detected in very high energy (VHE, E>100GeV) gamma-rays, only three belong to the subclass of Flat Spectrum Radio Quasars (FSRQs). MAGIC observed FSRQ PKS 1510-089 in February-April 2012 during a high activity state in the high energy (HE, E>100 MeV) gamma-ray band observed by AGILE and Fermi. MAGIC observations result in the detection of a source with significance of 6.0 sigma. In agreement with the previous VHE observations of the source, we find no statistically significant variability during the MAGIC observations in daily, weekly or monthly time scales. The other two known VHE FSRQs have shown daily scale to sub-hour variability. We study the multifrequency behaviour of the source at the epoch of MAGIC observation, collecting quasi-simultaneous data at radio and optical (GASP-WEBT and F-Gamma collaborations, REM, Steward, Perkins, Liverpool, OVRO and VLBA telescopes), X-ray (Swift satellite) and HE gamma-ray frequencies. The gamma-ray SED combining AGILE, Fermi and MAGIC data joins smoothly and shows no hint of a break. The multifrequency light curves suggest a common origin for the millimeter radio and HE gamma-ray emission and the HE gamma-ray flaring starts when the new component is ejected from the 43GHz VLBA core. The quasi-simultaneous multifrequency SED is modelled with a one-zone inverse Compton model. We study two different origins of the seed photons for the inverse Compton scattering, namely the infra-red torus and a slow sheath surrounding the jet around the VLBA core. Both models fit the data well. However, the fast HE gamma-ray variability requires that within the modelled large emitting region, there must exist more compact regions. We suggest that these observed signatures would be most naturally explained by a turbulent plasma flowing at a relativistic speed down the jet and crossing a standing conical shock.

We present a study of the very high energy (VHE; E > 100 GeV) gamma-ray emission of the blazar PKS 1424+240 observed with the MAGIC telescopes. The primary aim of this paper is the multiwavelength spectral characterization and modeling of this blazar, which is made particularly interesting by the recent discovery of a lower limit of its redshift of z > 0.6 and makes it a promising candidate to be the most distant VHE source. The source has been observed with the MAGIC telescopes in VHE gamma rays for a total observation time of ~33.6 h from 2009 to 2011. The source was marginally detected in VHE gamma rays during 2009 and 2010, and later, the detection was confirmed during an optical outburst in 2011. The combined significance of the stacked sample is ~7.2 sigma. The differential spectra measured during the different campaigns can be described by steep power laws with the indices ranging from 3.5 +/- 1.2 to 5.0 +/- 1.7. The MAGIC spectra corrected for the absorption due to the extragalactic background light connect smoothly, within systematic errors, with the mean spectrum in 2009-2011 observed at lower energies by the Fermi-LAT. The absorption-corrected MAGIC spectrum is flat with no apparent turn down up to 400 GeV. The multiwavelength light curve shows increasing flux in radio and optical bands that could point to a common origin from the same region of the jet. The large separation between the two peaks of the constructed non-simultaneous spectral energy distribution also requires an extremely high Doppler factor if an one zone synchrotron self-Compton model is applied. We find that a two-component synchrotron self-Compton model describes the spectral energy distribution of the source well, if the source is located at z~0.6.

Next-generation synoptic photometric surveys will yield unprecedented (for the astronomical community) volumes of data and the processes of discovery and rare-object identification are, by necessity, becoming more autonomous. Such autonomous searches can be used to find objects of interest applicable to a wide range of outstanding problems in astronomy, and in this paper we present the methods and results of a largely autonomous search for highly eccentric detached eclipsing binary systems in the Machine-learned ASAS Classification Catalog. 106 detached eclipsing binaries with eccentricities greater than 0.1 are presented, most of which are identified here for the first time. We also present new radial-velocity curves and absolute parameters for 6 of those systems with the long-term goal of increasing the number of highly eccentric systems with orbital solutions, thereby facilitating further studies of the tidal circularization process in binary stars.

Parsec-scale VLBA images of BL Lac at 15 GHz show that the jet contains a permanent quasi-stationary emission feature 0.26 mas (0.34 pc projected) from the core, along with numerous moving features. In projection, the tracks of the moving features cluster around an axis at position angle -166.6 deg that connects the core with the standing feature. The moving features appear to emanate from the standing feature in a manner strikingly similar to the results of numerical 2-D relativistic magneto-hydrodynamic (RMHD) simulations in which moving shocks are generated at a recollimation shock. Because of this, and the close analogy to the jet feature HST-1 in M87, we identify the standing feature in BL Lac as a recollimation shock. We assume that the magnetic field dominates the dynamics in the jet, and that the field is predominantly toroidal. From this we suggest that the moving features are compressions established by slow and fast mode magneto-acoustic MHD waves. We illustrate the situation with a simple model in which the slowest moving feature is a slow-mode wave, and the fastest feature is a fast-mode wave. In the model the beam has Lorentz factor about 3.5 in the frame of the host galaxy, and the fast mode wave has Lorentz factor about 1.6 in the frame of the beam. This gives a maximum apparent speed for the moving features 10c. In this model the Lorentz factor of the pattern in the galaxy frame is approximately 3 times larger than that of the beam itself.

Previously unremarkable, the extragalactic radio source GB 1310+487 showed a gamma-ray flare on 2009 November 18, reaching a daily flux of ~10^-6 photons/cm^2/s at energies E>100 MeV and becoming one of the brightest GeV sources for about two weeks. Its optical spectrum is not typical for a blazar, instead, it resembles those of narrow emission-line galaxies. We investigate changes of the object's radio-to-GeV spectral energy distribution (SED) during and after the prominent GeV flare with the aim to determine the nature of the object and constrain the origin of the variable high-energy emission. The data collected by the Fermi and AGILE satellites at gamma-ray energies, Swift at X-ray and ultraviolet, Kanata, NOT, and Keck telescopes at optical, OAGH and WISE at infrared, and IRAM 30m, OVRO 40m, Effelsberg 100m, RATAN-600, and VLBA at radio, are analysed together to trace the SED evolution on timescales of months. The gamma-ray/radio-loud narrow-line active galactic nucleus (AGN) is located at redshift z=0.638. It is shining through an unrelated foreground galaxy at z=0.500. The AGN light is likely amplified by a factor of a few because of gravitational lensing. The AGN SED shows a two-humped structure typical of blazars and gamma-ray-loud NLSy1 galaxies, with the high-energy (inverse-Compton) emission dominating by more than an order of magnitude over the low-energy (synchrotron) emission during gamma-ray flares. The difference between the two SED humps is smaller during the low-activity state. Fermi observations reveal a strong correlation between the gamma-ray flux and spectral index, with the hardest spectrum observed during the brightest gamma-ray state. If the gamma-ray flux is a mixture of synchrotron self-Compton (SSC) and external Compton (EC) emission, the observed GeV spectral variability may result from varying relative contributions of these two emission components.

We use photometry from the recent AllWISE Data Release of the Wide-field In-frared Survey Explorer (WISE) of 129 calibration stars, combined with prior distances obtained from the established $M_V-$[Fe/H] relation and Hubble Space Telescope trigonometric parallax, to derive mid-infrared period--luminosity relations for RR Lyrae pulsating variable stars. We derive relations in the W1, W2, and W3 wavebands (3.4, 4.6, and 12 {\mu}m, respectively), and for each of the two main RR Lyrae sub-types (RRab and RRc). We report an error on the period--luminosity relation slope for RRab stars of 0.2. We also fit posterior distances for the calibration catalog and find a median fractional distance error of 0.8 per cent.

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

We use optical data from the Palomar Transient Factory (PTF) and the Catalina Real-Time Transient Survey (CRTS) to study the variability of gamma-ray detected and non-detected objects in a large population of active galactic nuclei (AGN) selected from the Candidate Gamma-Ray Blazar Survey and Fermi Gamma-Ray Space Telescope catalogs. Our samples include 714 sources with PTF data and 1244 sources with CRTS data. We calculate the intrinsic modulation index to quantify the optical variability amplitude in these samples. We find the gamma-ray detected objects to be more variable than the non-detected ones. The flat spectrum radio quasars (FSRQs) are more variable than the BL Lac objects in our sample, but the significance of the difference depends on the sample used. When dividing the objects based on their synchrotron peak frequency, we find the low synchrotron peaked (LSP) objects to be significantly more variable than the high synchrotron peaked (HSP) ones, explaining the difference between the FSRQs and BL Lacs. This could be due to the LSPs being observed near their electron energy peak, while in the HSPs the emission is caused by lower energy electrons, which cool more slowly. We also find a significant correlation between the optical and gamma-ray fluxes that is stronger in the HSP BL Lacs than in the FSRQs. The FSRQs in our sample are also more Compton dominated than the HSP BL Lacs. These findings are consistent with models where the gamma-ray emission of HSP objects is produced by the synchrotron self-Compton mechanism, while the LSP objects need an additional external Compton component that increases the scatter in the flux-flux correlation.

We report on radio-to-gamma-ray observations during 2011 May-September of PMN J0948+0022, the first narrow-line Seyfert 1 (NLSy1) galaxy detected in gamma-rays by Fermi-LAT. Strong variability was observed in gamma-rays, with two flaring periods peaking on 2011 June 20 and July 28. The variability observed in optical and near-infrared seems to have no counterpart in gamma-rays. This different behaviour could be related to a bending and inhomogeneous jet or a turbulent extreme multi-cell scenario. The radio spectra showed a variability pattern typical of relativistic jets. The XMM spectrum shows that the emission from the jet dominates above 2 keV, while a soft X-ray excess is evident in the low-energy part of the X-ray spectrum. Models where the soft emission is partly produced by blurred reflection or Comptonisation of the thermal disc emission provide good fits to the data. The X-ray spectral slope is similar to that found in radio-quiet NLSy1, suggesting that a standard accretion disc is present, as expected from the high accretion rate. Except for the soft X-ray excess, unusual in jet-dominated AGNs, PMN J0948+0022 shows all characteristics of the blazar class.

We present a group of 3740 previously identified RR Lyrae variables well-observed with the Wide-field Infrared Survey Explorer (WISE). We explore how the shape of the generic RR Lyrae mid-infrared light curve evolves in period-space, comparing light curves in mid-infrared and optical bands. We find that optical light curves exhibit high amplitudes and a large spectrum of light curve shapes, while mid-infrared light curves have low amplitudes and uniform light curve shapes. From the period-space analysis, we hope to improve the classification methods of RR Lyrae variables and enable reliable discovery of these pulsators in the WISE catalog and future mid-infrared surveys such as the James Webb Space Telescope (JWST). We provide mid-infrared templates for typical RR Lyrae stars and demonstrate how these templates can be applied to improve estimates of mid-infrared RR Lyrae mean magnitude, which is used for distance measurement. This method of template fitting is particularly beneficial for improving observational efficiency. For example, using light curves with observational noise of 0.05 mag, we obtain the same level of accuracy in mean magnitude estimates for light curves randomly sampled at 12 data points with template fitting as with light curves randomly sampled at 20 data points with harmonic modelling.

We present results from four years of twice-weekly 15 GHz radio monitoring of about 1500 blazars with the Owens Valley Radio Observatory 40 m telescope. Using the intrinsic modulation index to measure variability amplitude, we find that, with $>\!6\sigma$ significance, the radio variability of radio-selected gamma-ray-loud blazars is stronger than that of gamma-ray-quiet blazars. Our extended data set also includes at least 21 months of data for all AGN with `clean' associations in the Fermi Large Area Telescope First AGN catalogue, 1LAC. With these additional data we examine the radio variability properties of a gamma-ray-selected blazar sample. Within this sample, we find no evidence for a connection between radio variability amplitude and optical classification. In contrast, for our radio-selected sample we find that the BL Lac object subpopulation is more variable than the flat spectrum radio quasar (FSRQ) subpopulation. Radio variability is found to correlate with the synchrotron peak frequency, with low- and intermediate-synchrotron-peaked blazars varying less than high-synchrotron-peaked ones. We find evidence for a significant negative correlation between redshift and radio variability among bright FSRQs.

During the LIGO and Virgo joint science runs in 2009-2010, gravitational wave (GW) data from three interferometer detectors were analyzed within minutes to select GW candidate events and infer their apparent sky positions. Target coordinates were transmitted to several telescopes for follow-up observations aimed at the detection of an associated optical transient. Images were obtained for eight such GW candidates. We present the methods used to analyze the image data as well as the transient search results. No optical transient was identified with a convincing association with any of these candidates, and none of the GW triggers showed strong evidence for being astrophysical in nature. We compare the sensitivities of these observations to several model light curves from possible sources of interest, and discuss prospects for future joint GW-optical observations of this type.

Fermi has provided the largest sample of gamma-ray selected blazars to date. In this work we use a uniformly selected set of 211 BL Lacertae (BL Lac) objects detected by it Fermi during its first year of operation. We have obtained redshift constraints for 206 out of the 211 BL Lacs in our sample making it the largest and most complete sample of BL Lacs available in the literature. We use this sample to determine the luminosity function of BL Lacs and its evolution with cosmic time. We find that for most BL Lac classes, the evolution is positive with a space density peaking at modest redshift (z~1.2). The low-luminosity, high-synchrotron peaked (HSP) BL Lacs are an exception, showing strong negative evolution, with number density increasing for z$\lesssim$0.5. Since this rise corresponds to a drop-off in the density of flat-spectrum radio quasars (FSRQs), a possible interpretation is that these HSPs represent an accretion-starved end-state of an earlier merger-driven gas-rich phase. We additionally find that the known BL Lac correlation between luminosity and photon spectral index persists after correction for the substantial observational selection effects with implications for the so called `blazar sequence'. Finally, estimating the beaming corrections to the luminosity function, we find that BL Lacs have an average Lorentz factor of $\gamma=6.1^{+1.1}_{-0.8}$, and that most are seen within 10$^{\circ}$ of the jet axis.

In September 2012, the high-synchrotron-peaked (HSP) blazar Markarian 421 underwent a rapid wideband radio flare, reaching nearly twice the brightest level observed in the centimeter band in over three decades of monitoring. In response to this event we carried out a five epoch centimeter- to millimeter-band multifrequency Very Long Baseline Array (VLBA) campaign to investigate the aftermath of this emission event. Rapid radio variations are unprecedented in this object and are surprising in an HSP BL Lac object. In this flare, the 15 GHz flux density increased with an exponential doubling time of about 9 days, then faded to its prior level at a similar rate. This is comparable with the fastest large-amplitude centimeter-band radio variability observed in any blazar. Similar flux density increases were detected up to millimeter bands. This radio flare followed about two months after a similarly unprecedented GeV gamma-ray flare (reaching a daily E>100 MeV flux of (1.2 +/- 0.7)x10^(-6) ph cm^(-2) s^(-1)) reported by the Fermi Large Area Telescope (LAT) collaboration, with a simultaneous tentative TeV detection by ARGO-YBJ. A cross-correlation analysis of long-term 15 GHz and LAT gamma-ray light curves finds a statistically significant correlation with the radio lagging ~40 days behind, suggesting that the gamma-ray emission originates upstream of the radio emission. Preliminary results from our VLBA observations show brightening in the unresolved core region and no evidence for apparent superluminal motions or substantial flux variations downstream.

We present the discovery and extensive early-time observations of the Type Ic supernova (SN) PTF12gzk. Our finely sampled light curves show a rise of 0.8mag within 2.5hr. Power-law fits [f(t)\sim(t-t_0)^n] to these data constrain the explosion date to within one day. We cannot rule out the expected quadratic fireball model, but higher values of n are possible as well for larger areas in the fit parameter space. Our bolometric light curve and a dense spectral sequence are used to estimate the physical parameters of the exploding star and of the explosion. We show that the photometric evolution of PTF12gzk is slower than that of most SNe Ic, and its high ejecta velocities (~30,000km/s four days after explosion) are closer to the observed velocities of broad-lined SNe Ic associated with gamma-ray bursts (GRBs) than to the observed velocities in normal Type Ic SNe. The high velocities are sustained through the SN early evolution, and are similar to those of GRB-SNe when the SN reach peak magnitude. By comparison with the spectroscopically similar SN 2004aw, we suggest that the observed properties of PTF12gzk indicate an initial progenitor mass of 25-35 solar mass and a large (5-10E51 erg) kinetic energy, close to the regime of GRB-SN properties. The host-galaxy characteristics are consistent with GRB-SN hosts, and not with normal SN Ic hosts as well, yet this SN does not show the broad lines over extended periods of time that are typical of broad-line Type Ic SNe.

We have developed a new Bayesian Markov Chain Monte Carlo method to deconvolve light curves of blazars into individual flares, including proper estimation of the fit errors. We use the method to fit 15\,GHz light curves obtained within the OVRO 40-m blazar monitoring program where a large number of AGN have been monitored since 2008 in support of the Fermi Gamma-Ray Space Telescope mission. The time scales obtained from the fitted models are used to calculate the variability brightness temperature of the sources. Additionally, we have calculated brightness temperatures of a sample of these objects using Very Long Baseline Array data from the MOJAVE survey. Combining these two data sets enables us to study the intrinsic brightness temperature distribution in these blazars at 15 GHz. Our preliminary results indicate that the mean intrinsic brightness temperature in a sample of 14 sources is near the equipartition brightness temperature of $\sim 10^{11}$K.

The OVRO 40-m telescope has been monitoring the 15 GHz radio flux density of over 1200 blazars since 2008. The 15 GHz light curve of the flat spectrum radio quasar J1359+4011 shows a strong and persistent quasi-periodic oscillation. The time-scale of the oscillation varies between 120 and 150 days over a $\sim4$ year time span. We interpret this as the active galactic nucleus mass-scaled analog of low-frequency quasi-periodic oscillations from Galactic microquasars, or as evidence of modulation of the accretion flow by thermal instabilites in the "inner" accretion disc.

The narrow-line Seyfert 1 galaxy SBS 0846+513 was first detected by the Large Area Telescope (LAT) on-board Fermi in 2011 June-July when it underwent a period of flaring activity. Since then, as Fermi continues to accumulate data on this source, its flux has been monitored on a daily basis. Two further gamma-ray flaring episodes from SBS 0846+513 were observed in 2012 May and August, reaching a daily peak flux integrated above 100 MeV of (50+/-12)x10^-8 ph/cm^2/s, and (73+/-14)x10^-8 ph/cm^2/s on May 24 and August 7, respectively. Three outbursts were detected at 15 GHz by the Owens Valley Radio Observatory 40-m telescope in 2012 May, 2012 October, and 2013 January, suggesting a complex connection with the gamma-ray activity. The most likely scenario suggests that the 2012 May gamma-ray flare may not be directly related to the radio activity observed over the same period, while the two gamma-ray flaring episodes may be related to the radio activity observed at 15 GHz in 2012 October and 2013 January. The gamma-ray flare in 2012 May triggered Swift observations that confirmed that SBS 0846+513 was also exhibiting high activity in the optical, UV and X-ray bands, thus providing a firm identification between the gamma-ray source and the lower-energy counterpart. We compared the spectral energy distribution (SED) of the flaring state in 2012 May with that of a quiescent state. The two SEDs, modelled as an external Compton component of seed photons from a dust torus, could be fitted by changing the electron distribution parameters as well as the magnetic field. No significant evidence of thermal emission from the accretion disc has been observed. Interestingly, in the 5 GHz radio luminosity vs. synchrotron peak frequency plot SBS 0846+513 seems to lie in the flat spectrum radio quasar part of the so-called `blazar sequence'.

We describe the parsec-scale kinematics of 200 AGN jets based on 15 GHz VLBA data obtained between 1994 Aug 31 and 2011 May 1. We present new VLBA 15 GHz images of these and 59 additional AGN from the MOJAVE and 2 cm Survey programs. Nearly all of the 60 most heavily observed jets show significant changes in their innermost position angle over a 12 to 16 year interval, ranging from 10 deg to 150 deg on the sky, corresponding to intrinsic variations of ~0.5 deg to ~2 deg. The BL Lac jets show smaller variations than quasars. Roughly half of the heavily observed jets show systematic position angle trends with time, and 20 show indications of oscillatory behavior. The time spans of the data sets are too short compared to the fitted periods (5 to 12 y), however, to reliably establish periodicity. The rapid changes and large jumps in position angle seen in many cases suggest that the superluminal AGN jet features occupy only a portion of the entire jet cross section, and may be energized portions of thin instability structures within the jet. We have derived vector proper motions for 887 moving features in 200 jets having at least five VLBA epochs. For 557 well-sampled features, there are sufficient data to additionally study possible accelerations. We find that the moving features are generally non-ballistic, with 70% of the well-sampled features showing either significant accelerations or non-radial motions. Inward motions are rare (2% of all features), are slow (<0.1 mas per y), are more prevalent in BL Lac jets, and are typically found within 1 mas of the unresolved core feature. There is a general trend of increasing apparent speed with distance down the jet for both radio galaxies and BL Lac objects. In most jets, the speeds of the features cluster around a characteristic value, yet there is a considerable dispersion in the distribution. (abridged)

[Abridged] Context. To construct and interpret the spectral energy distribution (SED) of BL Lacertae objects, simultaneous broad-band observations are mandatory. Aims. We present the results of a dedicated multi-wavelength study of the high-frequency peaked BL Lacertae (HBL) object and known TeV emitter 1ES 2344+514 by means of a pre-organised campaign. Methods. The observations were conducted during simultaneous visibility windows of MAGIC and AGILE in late 2008. The measurements were complemented by Mets\"ahovi, RATAN-600, KVA+Tuorla, Swift and VLBA pointings. Additional coverage was provided by the ongoing long-term F-GAMMA and MOJAVE programs, the OVRO 40-m and CrAO telescopes as well as the Fermi satellite. The obtained SEDs are modelled using a one-zone as well as a self-consistent two-zone synchrotron self-Compton model. Results. 1ES 2344+514 was found at very low flux states in both X-rays and very high energy gamma rays. Variability was detected in the low frequency radio and X-ray bands only, where for the latter a small flare was observed. The X-ray flare was possibly caused by shock acceleration characterised by similar cooling and acceleration time scales. MOJAVE VLBA monitoring reveals a static jet whose components are stable over time scales of eleven years, contrary to previous findings. There appears to be no significant correlation between the 15 GHz and R-band monitoring light curves. The observations presented here constitute the first multi-wavelength campaign on 1ES 2344+514 from radio to VHE energies and one of the few simultaneous SEDs during low activity states. The quasi-simultaneous Fermi-LAT data poses some challenges for SED modelling, but in general the SEDs are described well by both applied models. The resulting parameters are typical for TeV emitting HBLs. Consequently it remains unclear whether a so-called quiescent state was found in this campaign.

We report on the first detection of the theoretically-predicted rare phenomenon of multiple parsec-scale imaging of an active galactic nucleus induced by refractive effects due to localized foreground electron density enhancements, e.g., in an AU-scale plasma lens(es) in the ionized component of the Galactic interstellar medium. We detected multiple imaging in the low galactic latitude (b=-2 deg) quasar 2023+335 from the 15.4 GHz MOJAVE observations when the source was undergoing an ESE. While the parsec-scale jet of the source normally extends along PA -20 deg, in the 28 May 2009 and 23 July 2009 images a highly significant multi-component pattern of secondary images is stretched out nearly along the constant galactic latitude line with a local PA 40 deg, indicating that the direction of relative motion of the plasma lens is close to orbital. Weaker but still detectable imaging patterns at similar position angles are sporadically manifest at several other epochs. Modeling the ESE that occurred in early 2009 and lasted ~0.14 yr, we determined that the foreground screen has a double-lens structure, with proper motion (~6.8 mas/yr), and angular size (~0.27 mas). We also found that the angular separation between the two brightest sub-images roughly follows a wavelength-squared dependence expected from plasma scattering. Furthermore, by analyzing archival non-simultaneous VLBA observations covering a wide frequency range from 1.4 to 86 GHz, we found that the scattered angular size of the VLBI core follows a $\nu^{-1.89}$ dependence, implying the presence of a turbulent, refractive dominated scattering screen that has a confined structure or is truncated transverse to the line of sight toward 2023+335.

We report on multifrequency observations of the gamma-ray emitting narrow-line Seyfert 1 galaxy PKS 1502+036 performed from radio to gamma-rays during 2008 August-2012 November by Fermi-LAT, Swift (XRT and UVOT), OVRO, VLBA, and VLA. No significant variability has been observed in gamma-rays, with 0.1-100 GeV flux that ranged between (3-7)x10^-8 ph/cm^2/s using 3-month time bins. The photon index of the LAT spectrum (Gamma=2.60+/-0.06) and the apparent isotropic gamma-ray luminosity, L(0.1-100 GeV)= 7.8x10^45 erg/s, over 51 months are typical of a flat spectrum radio quasar. The radio spectral variability and the one-sided structure, in addition to the observed gamma-ray luminosity, suggest a relativistic jet with a high Doppler factor. In contrast to SBS 0846+513, the VLBA at 15 GHz did not observe superluminal motion for PKS 1502+036. Despite having the optical characteristics typical of a narrow-line Seyfert 1 galaxy, radio and gamma-ray properties of PKS 1502+036 are found to be similar to those of a blazar at the low end of the black hole mass distribution for blazars. This is in agreement with what has been found in the case of the other gamma-ray emitting narrow-line Seyfert 1 SBS 0846+513.

We present broadband observations and analysis of Swift gamma-ray burst (GRB) 120119A. Our early-time afterglow detections began under 15 s after the burst in the host frame (redshift z = 1.73), and they yield constraints on the burst energetics and local environment. Late-time afterglow observations of the burst show evidence for a moderate column of dust (A_V ~ 1.1 mag) similar to, but statistically distinct from, dust seen along Small Magellanic Cloud sightlines. Deep late-time observations reveal a dusty, rapidly star-forming host galaxy. Most notably, our early-time observations exhibit a significant red-to-blue colour change in the first ~200 s after the trigger at levels heretofore unseen in GRB afterglows. This colour change, which is coincident with the final phases of the prompt emission, is a hallmark prediction of the photodestruction of dust in GRB afterglows. We test whether dust-destruction signatures are significantly distinct from other sources of colour change, namely a change in the intrinsic spectral index {\beta}. We find that a time-varying power-law spectrum alone cannot adequately describe the observed colour change, and allowing for dust destruction (via a time-varying A_V) significantly improves the fit. While not definitively ruling out other possibilities, this event provides the best support yet for the direct detection of dust destruction in the local environment of a GRB.

[Abridged] We study the spectral properties of intermediate mass galaxies as a function of colour and morphology. We use Galaxy Zoo to define three morphological classes of galaxies, namely early-types (ellipticals), late-type (disk-dominated) face-on spirals and early-type (bulge-dominated) face-on spirals. We classify these galaxies as blue or red according to their SDSS g-r colour and use the spectral fitting code VESPA to calculate time-resolved star-formation histories, metallicity and total starlight dust extinction from their SDSS fibre spectra. We find that red late-type spirals show less star-formation in the last 500 Myr than blue late-type spirals by up to a factor of three, but share similar star-formation histories at earlier times. This decline in recent star-formation explains their redder colour: their chemical and dust content are the same. We postulate that red late-type spirals are recent descendants of blue late-type spirals, with their star-formation curtailed in the last 500 Myrs. The red late-type spirals are however still forming stars approximately 17 times faster than red ellipticals over the same period. Red early-type spirals lie between red late-type spirals and red ellipticals in terms of recent-to-intermediate star-formation and dust content. Therefore, it is plausible that these galaxies represent an evolutionary link between these two populations. They are more likely to evolve directly into red ellipticals than red late-type spirals. Blue ellipticals show similar star-formation histories as blue spirals (regardless of type), except they have formed less stars in the last 100 Myrs. However, blue ellipticals have different dust content, which peaks at lower extinction values than all spiral galaxies.

Blazars are powerful, variable emitters from radio to gamma-ray wavelengths. Even though the general picture of synchrotron emission at low energies and inverse Compton at the high energy end is well established, many important aspects of these remarkable objects are still not well understood. For example, even the location of the gamma-ray emission region is still not clearly established, with some theories locating it close to the black hole/accretion disk while others place it at parsec scales in the radio jet. Since mid-2007 we have carried out a large scale monitoring program at 15 GHz using the OVRO 40 m telescope. We are currently observing about 1700 blazars twice per week. The sample includes all the Fermi-LAT detected blazars north of declination -20 degrees. Here, we study the existence of correlated variability between these two bands for 86 sources bright enough to be detected weekly by LAT. The existence of correlated variability can be interpreted as an indication of a related spatial locations for the radio and gamma-ray emission, making the evaluation of its statistical significance a key goal of our program. A study of the statistical significance of these cross-correlations is presented along with a discussion of the Monte Carlo simulations used to evaluate them. More information about the conditions on the radio emission zone can be obtained through polarization monitoring which tells us about the configuration of the magnetic fields in this region. To study radio polarization variability we are building KuPol, a radio polarization receiver for the 12 to 18 GHz band that will replace the current total power receiver at the OVRO 40 meter telescope.

The Q/U Imaging ExperimenT (QUIET) has observed the cosmic microwave background (CMB) at 43 and 95GHz. The 43-GHz results have been published in QUIET Collaboration et al. (2011), and here we report the measurement of CMB polarization power spectra using the 95-GHz data. This data set comprises 5337 hours of observations recorded by an array of 84 polarized coherent receivers with a total array sensitivity of 87 uK sqrt(s). Four low-foreground fields were observed, covering a total of ~1000 square degrees with an effective angular resolution of 12.8', allowing for constraints on primordial gravitational waves and high-signal-to-noise measurements of the E-modes across three acoustic peaks. The data reduction was performed using two independent analysis pipelines, one based on a pseudo-Cl (PCL) cross-correlation approach, and the other on a maximum-likelihood (ML) approach. All data selection criteria and filters were modified until a predefined set of null tests had been satisfied before inspecting any non-null power spectrum. The results derived by the two pipelines are in good agreement. We characterize the EE, EB and BB power spectra between l=25 and 975 and find that the EE spectrum is consistent with LCDM, while the BB power spectrum is consistent with zero. Based on these measurements, we constrain the tensor-to-scalar ratio to r=1.1+0.9-0.8 (r<2.8 at 95% C.L.) as derived by the ML pipeline, and r=1.2+0.9-0.8 (r<2.7 at 95% C.L.) as derived by the PCL pipeline. In one of the fields, we find a correlation with the dust component of the Planck Sky Model, though the corresponding excess power is small compared to statistical errors. Finally, we derive limits on all known systematic errors, and demonstrate that these correspond to a tensor-to-scalar ratio smaller than r=0.01, the lowest level yet reported in the literature.

We present the results of a series of radio, optical, X-ray and gamma-ray observations of the BL Lac object S50716+714 carried out between April 2007 and January 2011. The multi-frequency observations were obtained using several ground and space based facilities. The intense optical monitoring of the source reveals faster repetitive variations superimposed on a long-term variability trend at a time scale of ~350 days. Episodes of fast variability recur on time scales of ~ 60-70 days. The intense and simultaneous activity at optical and gamma-ray frequencies favors the SSC mechanism for the production of the high-energy emission. Two major low-peaking radio flares were observed during this high optical/gamma-ray activity period. The radio flares are characterized by a rising and a decaying stage and are in agreement with the formation of a shock and its evolution. We found that the evolution of the radio flares requires a geometrical variation in addition to intrinsic variations of the source. Different estimates yield a robust and self-consistent lower limits of \delta > 20 and equipartition magnetic field B_eq > 0.36 G. Causality arguments constrain the size of emission region \theta < 0.004 mas. We found a significant correlation between flux variations at radio frequencies with those at optical and gamma-rays. The optical/GeV flux variations lead the radio variability by ~65 days. The longer time delays between low-peaking radio outbursts and optical flares imply that optical flares are the precursors of radio ones. An orphan X-ray flare challenges the simple, one-zone emission models, rendering them too simple. Here we also describe the spectral energy distribution modeling of the source from simultaneous data taken through different activity periods.

We report on spectroscopic observations covering most of the 475 BL Lacs in the 2nd Fermi LAT catalog of AGN. Including archival measurements (correcting several erroneous literature values) we now have spectroscopic redshifts for 44% of the BL Lacs. We establish firm lower redshift limits via intervening absorption systems and statistical lower limits via searches for host galaxies for an additional 51% of the sample leaving only 5% of the BL Lacs unconstrained. The new redshifts raise the median spectroscopic z from 0.23 to 0.33 and include redshifts as large as z=2.471. Spectroscopic redshift minima from intervening absorbers have ~ z= 0.70, showing a substantial fraction at large z and arguing against strong negative evolution. We find that detected BL Lac hosts are bright ellipticals with black hole masses M_\bullet ~ 10^{8.5-9}, substantially larger than the mean of optical AGN and LAT Flat Spectrum Radio Quasar samples. A slow increase in M_\bullet with z may be due to selection bias. We find that the power-law dominance of the optical spectrum extends to extreme values, but this does not strongly correlate with the gamma-ray properties, suggesting that strong beaming is the primary cause of the range in continuum dominance.

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

We investigate the radio and gamma-ray variability of the flat spectrum radio quasar PKS 1510-089 in the time range between 2010 November and 2012 January. In this period the source showed an intense activity, with two major gamma-ray flares detected in 2011 July and October. During the latter episode both the gamma-ray and the radio flux density reached their historical peak. Multiwavelength analysis shows a rotation of about 380 deg of the optical polarization angle close in time with the rapid and strong gamma-ray flare in 2011 July. An enhancement of the optical emission and an increase of the fractional polarization both in the optical and in radio bands is observed about three weeks later, close in time with another gamma-ray outburst. On the other hand, after 2011 September a huge radio outburst has been detected, first in the millimeter regime followed with some time delay at centimeter down to decimeter wavelengths. This radio flare is characterized by a rising and a decaying stage, in agreement with the formation of a shock and its evolution, as a consequence of expansion and radiative cooling. If the gamma-ray flare observed in 2011 October is related to this radio outburst, then this strongly indicates that the region responsible for the gamma-ray variability is not within the broad line, but a few parsecs downstream along the jet.

We present more than three years of observations at different frequencies, from radio to high-energy gamma-rays, of the Narrow-Line Seyfert 1 (NLS1) Galaxy PMN J0948+0022 (z=0.585). This source is the first NLS1 detected at energies above 100 MeV and therefore can be considered the prototype of this emerging new class of gamma-ray emitting active galactic nuclei (AGN). The observations performed from 2008 August 1 to 2011 December 31 confirmed that PMN J0948+0022 generates a powerful relativistic jet, able to develop an isotropic luminosity at gamma-rays of the order of 10^48 erg s^-1, at the level of powerful quasars. The evolution of the radiation emission of this source in 2009 and 2010 followed the canonical expectations of relativistic jets, with correlated multiwavelength variability (gamma-rays followed by radio emission after a few months), but it was difficult to retrieve a similar pattern in the light curves of 2011. The comparison of gamma-ray spectra before and including 2011 data suggested that there was a softening of the high-energy spectral slope. We selected five specific epochs to be studied by modelling the broad-band spectrum, characterised by an outburst at gamma-rays or very low/high flux at other wavelengths. The observed variability can largely be explained either by changes in the injected power, the bulk Lorentz factor of the jet or the electron spectrum. The characteristic time scale of doubling/halving flux ranges from a few days to a few months, depending on the frequency and the sampling rate. The shortest doubling time scale at gamma-rays is 2.3+-0.5 days. These small values underline the need of highly-sampled multiwavelength campaigns to better understand the physics of these sources.

Modern time-domain surveys continuously monitor large swaths of the sky to look for astronomical variability. Astrophysical discovery in such data sets is complicated by the fact that detections of real transient and variable sources are highly outnumbered by bogus detections caused by imperfect subtractions, atmospheric effects and detector artefacts. In this work we present a machine learning (ML) framework for discovery of variability in time-domain imaging surveys. Our ML methods provide probabilistic statements, in near real time, about the degree to which each newly observed source is astrophysically relevant source of variable brightness. We provide details about each of the analysis steps involved, including compilation of the training and testing sets, construction of descriptive image-based and contextual features, and optimization of the feature subset and model tuning parameters. Using a validation set of nearly 30,000 objects from the Palomar Transient Factory, we demonstrate a missed detection rate of at most 7.7% at our chosen false-positive rate of 1% for an optimized ML classifier of 23 features, selected to avoid feature correlation and over-fitting from an initial library of 42 attributes. Importantly, we show that our classification methodology is insensitive to mis-labelled training data up to a contamination of nearly 10%, making it easier to compile sufficient training sets for accurate performance in future surveys. This ML framework, if so adopted, should enable the maximization of scientific gain from future synoptic survey and enable fast follow-up decisions on the vast amounts of streaming data produced by such experiments.

We present the results of a machine-learning (ML) based search for new R Coronae Borealis (RCB) stars and DY Persei-like stars (DYPers) in the Galaxy using cataloged light curves from the All-Sky Automated Survey (ASAS) Catalog of Variable Stars (ACVS). RCB stars - a rare class of hydrogen-deficient carbon-rich supergiants - are of great interest owing to the insights they can provide on the late stages of stellar evolution. DYPers are possibly the low-temperature, low-luminosity analogs to the RCB phenomenon, though additional examples are needed to fully establish this connection. While RCB stars and DYPers are traditionally identified by epochs of extreme dimming that occur without regularity, the ML search framework more fully captures the richness and diversity of their photometric behavior. We demonstrate that our ML method can use newly discovered RCB stars to identify additional candidates within the same data set. Our search yields 15 candidates that we consider likely RCB stars/DYPers: new spectroscopic observations confirm that four of these candidates are RCB stars and four are DYPers. Our discovery of four new DYPers increases the number of known Galactic DYPers from two to six; noteworthy is that one of the new DYPers has a measured parallax and is m ~ 7 mag, making it the brightest known DYPer to date. Future observations of these new DYPers should prove instrumental in establishing the RCB connection. We consider these results, derived from a machine-learned probabilistic classification catalog, as an important proof-of-concept for the efficient discovery of rare sources with time-domain surveys.

The Q/U Imaging ExperimenT (QUIET) is designed to measure polarization in the Cosmic Microwave Background, targeting the imprint of inflationary gravitational waves at large angular scales (~ 1 degree). Between 2008 October and 2010 December, two independent receiver arrays were deployed sequentially on a 1.4 m side-fed Dragonian telescope. The polarimeters which form the focal planes use a highly compact design based on High Electron Mobility Transistors (HEMTs) that provides simultaneous measurements of the Stokes parameters Q, U, and I in a single module. The 17-element Q-band polarimeter array, with a central frequency of 43.1 GHz, has the best sensitivity (69 uK sqrt(s)) and the lowest instrumental systematic errors ever achieved in this band, contributing to the tensor-to-scalar ratio at r < 0.1. The 84-element W-band polarimeter array has a sensitivity of 87 uK sqrt(s) at a central frequency of 94.5 GHz. It has the lowest systematic errors to date, contributing at r < 0.01. The two arrays together cover multipoles in the range l= 25-975. These are the largest HEMT-based arrays deployed to date. This article describes the design, calibration, performance of, and sources of systematic error for the instrument.

This paper presents both the result of a search for fossil systems (FSs) within the XMM Cluster Survey and the Sloan Digital Sky Survey and the results of a study of the stellar mass assembly and stellar populations of their fossil galaxies. In total, 17 groups and clusters are identified at z < 0.25 with large magnitude gaps between the first and fourth brightest galaxies. All the information necessary to classify these systems as fossils is provided. For both groups and clusters, the total and fractional luminosity of the brightest galaxy is positively correlated with the magnitude gap. The brightest galaxies in FSs (called fossil galaxies) have stellar populations and star formation histories which are similar to normal brightest cluster galaxies (BCGs). However, at fixed group/cluster mass, the stellar masses of the fossil galaxies are larger compared to normal BCGs, a fact that holds true over a wide range of group/cluster masses. Moreover, the fossil galaxies are found to contain a significant fraction of the total optical luminosity of the group/cluster within 0.5R200, as much as 85%, compared to the non-fossils, which can have as little as 10%. Our results suggest that FSs formed early and in the highest density regions of the universe and that fossil galaxies represent the end products of galaxy mergers in groups and clusters. The online FS catalog can be found at http://www.astro.ljmu.ac.uk/~xcs/Harrison2012/XCSFSCat.html.

The Q/U Imaging ExperimenT (QUIET) employs coherent receivers at 43GHz and 95GHz, operating on the Chajnantor plateau in the Atacama Desert in Chile, to measure the anisotropy in the polarization of the CMB. QUIET primarily targets the B modes from primordial gravitational waves. The combination of these frequencies gives sensitivity to foreground contributions from diffuse Galactic synchrotron radiation. Between 2008 October and 2010 December, >10,000hours of data were collected, first with the 19-element 43GHz array (3458hours) and then with the 90-element 95GHz array. Each array observes the same four fields, selected for low foregrounds, together covering ~1000deg^2. This paper reports initial results from the 43GHz receiver which has an array sensitivity to CMB fluctuations of 69uK sqrt(s). The data were extensively studied with a large suite of null tests before the power spectra, determined with two independent pipelines, were examined. Analysis choices, including data selection, were modified until the null tests passed. Cross correlating maps with different telescope pointings is used to eliminate a bias. This paper reports the EE, BB and EB power spectra in the multipole range ell=25-475. With the exception of the lowest multipole bin for one of the fields, where a polarized foreground, consistent with Galactic synchrotron radiation, is detected with 3sigma significance, the E-mode spectrum is consistent with the LCDM model, confirming the only previous detection of the first acoustic peak. The B-mode spectrum is consistent with zero, leading to a measurement of the tensor-to-scalar ratio of r=0.35+1.06-0.87. The combination of a new time-stream double-demodulation technique, Mizuguchi-Dragone optics, natural sky rotation, and frequent boresight rotation leads to the lowest level of systematic contamination in the B-mode power so far reported, below the level of r=0.1

We report Fermi-LAT observations of the radio-loud AGN SBS 0846+513 (z=0.5835), optically classified as a Narrow-Line Seyfert 1 galaxy, together with new and archival radio-to-X-ray data. The source was not active at gamma-ray energies during the first two years of Fermi operation. A significant increase in activity was observed during 2010 October-2011 August. In particular a strong gamma-ray flare was observed in 2011 June reaching an isotropic gamma-ray luminosity (0.1-300 GeV) of 1.0x10^48 erg/s, comparable to that of the brightest flat spectrum radio quasars, and showing spectral evolution in gamma rays. An apparent superluminal velocity of (8.2+/-1.5)c in the jet was inferred from 2011-2012 VLBA images, suggesting the presence of a highly relativistic jet. Both the power released by this object during the flaring activity and the apparent superluminal velocity are strong indications of the presence of a relativistic jet as powerful as those of blazars. In addition, variability and spectral properties in radio and gamma-ray bands indicate blazar-like behaviour, suggesting that, except for some distinct optical characteristics, SBS 0846+513 could be considered as a young blazar at the low end of the blazar's black hole mass distribution.

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

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

We present simultaneous Planck, Swift, Fermi, and ground-based data for 105 blazars belonging to three samples with flux limits in the soft X-ray, hard X-ray, and gamma-ray bands. Our unique data set has allowed us to demonstrate that the selection method strongly influences the results, producing biases that cannot be ignored. Almost all the BL Lac objects have been detected by Fermi-LAT, whereas ~40% of the flat-spectrum radio quasars (FSRQs) in the radio, soft X-ray, and hard X-ray selected samples are still below the gamma-ray detection limit even after integrating 27 months of Fermi-LAT data. The radio to sub-mm spectral slope of blazars is quite flat up to ~70GHz, above which it steepens to <\alpha>~-0.65. BL Lacs have significantly flatter spectra than FSRQs at higher frequencies. The distribution of the rest-frame synchrotron peak frequency (\nupS) in the SED of FSRQs is the same in all the blazar samples with <\nupS>=10^13.1 Hz, while the mean inverse-Compton peak frequency, <\nupIC>, ranges from 10^21 to 10^22 Hz. The distributions of \nupS and of \nupIC of BL Lacs are much broader and are shifted to higher energies than those of FSRQs and strongly depend on the selection method. The Compton dominance of blazars ranges from ~0.2 to ~100, with only FSRQs reaching values >3. Its distribution is broad and depends strongly on the selection method, with gamma-ray selected blazars peaking at ~7 or more, and radio-selected blazars at values ~1, thus implying that the assumption that the blazar power is dominated by high-energy emission is a selection effect. Simple SSC models cannot explain the SEDs of most of the gamma-ray detected blazars in all samples. The SED of the blazars that were not detected by Fermi-LAT may instead be consistent with SSC emission. Our data challenge the correlation between bolometric luminosity and \nupS predicted by the blazar sequence.

The parsec-scale radio properties of blazars detected by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope have been investigated using observations with the Very Long Baseline Array (VLBA). Comparisons between LAT and non-LAT detected samples were made using contemporaneous data. In total, 232 sources were used in the LAT-detected sample. This very large, radio flux-limited sample of active galactic nuclei (AGN) provides insights into the mechanism that produces strong gamma-ray emission. It has been found that LAT-detected BL Lac objects are very similar to the non-LAT BL Lac objects in most properties, although LAT BL Lac objects may have longer jets. The LAT flat spectrum radio quasars (FSRQs) are significantly different from non-LAT FSRQs and are likely extreme members of the FSRQ population. Contemporaneous observations showed a strong correlation, whereas no correlation is found using archival radio data. Most of the differences between the LAT and non-LAT populations are related to the cores of the sources, indicating that the gamma-ray emission may originate near the base of the jets (i.e., within a few pc of the central engine). There is some indication that LAT-detected sources may have larger jet opening angles than the non-LAT sources. Strong core polarization is significantly more common among the LAT sources, suggesting that gamma-ray emission is related to strong, uniform magnetic fields at the base of the jets of the blazars. Observations of sources in two epochs indicate that core fractional polarization was higher when the objects were detected by the LAT. Included in our sample are several non-blazar AGN such as 3C84, M82, and NGC 6251.

A large sample of known and likely gamma-ray blazars has been monitored twice per week since late 2007 at 15 GHz with the Owens Valley Radio Observatory (OVRO) 40-meter Telescope. The sample contains about 1700 sources, including the initial sample of 1158 sources above declination -20 degrees from the Candidate Gamma-Ray Blazar Survey (CGRaBS) plus all the blazars associated with Fermi-LAT detections as released in the Fermi AGN catalogs. Using statistical likelihood analyses, we compare the variability amplitude for various sub-populations within our sample. These include comparisons of gamma-ray-loud versus quiet objects, BL Lac objects versus flat-spectrum radio quasars, and a study of the variability amplitude trend with redshift. To learn about the location of the gamma-ray emission region we study the significance of peaks in the radio/gamma-ray cross-correlation using Monte Carlo simulations. First results for 52 sources with data from both the high-confidence Fermi Large Area Telescope Bright AGN Sample and the first 2 years of our monitoring program are presented. We find that assuming a power spectral density with power law slope of -2 at 15 GHz and -1.5 at gamma-ray energies, 7 of our objects show cross-correlations at the 3sigma level. We are now studying the physical significance of these correlations by further exploring the range of power law slopes that are consistent with the data. An extension of this to a larger sample and longer light curves is underway and preliminary results are presented. We also describe KuPol, the new digital Ku-band receiver being constructed for the 40-meter telescope. This new receiver will provide total intensity and linear polarization measurements over the 12-18 GHz band, with 16 MHz spectral resolution. The polarization data will provide important clues about the magnetic field configuration in the radio emission region.

With growing data volumes from synoptic surveys, astronomers must become more abstracted from the discovery and introspection processes. Given the scarcity of follow-up resources, there is a particularly sharp onus on the frameworks that replace these human roles to provide accurate and well-calibrated probabilistic classification catalogs. Such catalogs inform the subsequent follow-up, allowing consumers to optimize the selection of specific sources for further study and permitting rigorous treatment of purities and efficiencies for population studies. Here, we describe a process to produce a probabilistic classification catalog of variability with machine learning from a multi-epoch photometric survey. In addition to producing accurate classifications, we show how to estimate calibrated class probabilities, and motivate the importance of probability calibration. We also introduce a methodology for feature-based anomaly detection, which allows discovery of objects in the survey that do not fit within the predefined class taxonomy. Finally, we apply these methods to sources observed by the All Sky Automated Survey (ASAS), and unveil the Machine-learned ASAS Classification Catalog (MACC), which is a 28-class probabilistic classification catalog of 50,124 ASAS sources. We estimate that MACC achieves a sub-20% classification error rate, and demonstrate that the class posterior probabilities are reasonably calibrated. MACC classifications compare favorably to the classifications of several previous domain-specific ASAS papers and to the ASAS Catalog of Variable Stars, which had classified only 24% of those sources into one of 12 science classes. The MACC is publicly available at http://www.bigmacc.info.

Whether a correlation exists between the radio and gamma-ray flux densities of blazars is a long-standing question, and one that is difficult to answer confidently because of various observational biases which may either dilute or apparently enhance any intrinsic correlation between radio and gamma-ray luminosities. We introduce a novel method of data randomization to evaluate quantitatively the effect of these biases and to assess the intrinsic significance of an apparent correlation between radio and gamma-ray flux densities of blazars. The novelty of the method lies in a combination of data randomization in luminosity space (to ensure that the randomized data are intrinsically, and not just apparently, uncorrelated) and significance assessment in flux space (to explicitly avoid Malmquist bias and automatically account for the limited dynamical range in both frequencies). The method is applicable even to small samples that are not selected with strict statistical criteria. For larger samples we describe a variation of the method in which the sample is split in redshift bins, and the randomization is applied in each bin individually; this variation is designed to yield the equivalent to luminosity-function sampling of the underlying population in the limit of very large, statistically complete samples. We show that for a smaller number of redshift bins, the method yields a worse significance, and in this way it is conservative in that it does not assign a stronger, artificially enhanced significance. We demonstrate how our test performs as a function of number of sources, strength of correlation, and number of redshift bins used, and we show that while our test is robust against common-distance biases and associated false positives for uncorrelated data, it retains the power of other methods in rejecting the null hypothesis of no correlation for correlated data.

Parameter estimation in astrophysics often requires the use of complex physical models. In this paper we study the problem of estimating the parameters that describe star formation history (SFH) in galaxies. Here, high-dimensional spectral data from galaxies are appropriately modeled as linear combinations of physical components, called simple stellar populations (SSPs), plus some nonlinear distortions. Theoretical data for each SSP is produced for a fixed parameter vector via computer modeling. Though the parameters that define each SSP are continuous, optimizing the signal model over a large set of SSPs on a fine parameter grid is computationally infeasible and inefficient. The goal of this study is to estimate the set of parameters that describes the SFH of each galaxy. These target parameters, such as the average ages and chemical compositions of the galaxy's stellar populations, are derived from the SSP parameters and the component weights in the signal model. Here, we introduce a principled approach of choosing a small basis of SSP prototypes for SFH parameter estimation. The basic idea is to quantize the vector space and effective support of the model components. In addition to greater computational efficiency, we achieve better estimates of the SFH target parameters. In simulations, our proposed quantization method obtains a substantial improvement in estimating the target parameters over the common method of employing a parameter grid. Sparse coding techniques are not appropriate for this problem without proper constraints, while constrained sparse coding methods perform poorly for parameter estimation because their objective is signal reconstruction, not estimation of the target parameters.

Efficient and automated classification of periodic variable stars is becoming increasingly important as the scale of astronomical surveys grows. Several recent papers have used methods from machine learning and statistics to construct classifiers on databases of labeled, multi--epoch sources with the intention of using these classifiers to automatically infer the classes of unlabeled sources from new surveys. However, the same source observed with two different synoptic surveys will generally yield different derived metrics (features) from the light curve. Since such features are used in classifiers, this survey-dependent mismatch in feature space will typically lead to degraded classifier performance. In this paper we show how and why feature distributions change using OGLE and \textit{Hipparcos} light curves. To overcome survey systematics, we apply a method, \textit{noisification}, which attempts to empirically match distributions of features between the labeled sources used to construct the classifier and the unlabeled sources we wish to classify. Results from simulated and real--world light curves show that noisification can significantly improve classifier performance. In a three--class problem using light curves from \textit{Hipparcos} and OGLE, noisification reduces the classifier error rate from 27.0% to 7.0%. We recommend that noisification be used for upcoming surveys such as Gaia and LSST and describe some of the promises and challenges of applying noisification to these surveys.

A Bayesian approach to calibrating period-luminosity (PL) relations has substantial benefits over generic least-squares fits. In particular, the Bayesian approach takes into account the full prior distribution of the model parameters, such as the a priori distances, and refits these parameters as part of the process of settling on the most highly-constrained final fit. Additionally, the Bayesian approach can naturally ingest data from multiple wavebands and simultaneously fit the parameters of PL relations for each waveband in a procedure that constrains the parameter posterior distributions so as to minimize the scatter of the final fits appropriately in all wavebands. Here we describe the generalized approach to Bayesian model fitting and then specialize to a detailed description of applying Bayesian linear model fitting to the mid-infrared PL relations of RR Lyrae variable stars. For this example application we quantify the improvement afforded by using a Bayesian model fit. We also compare distances previously predicted in our example application to recently published parallax distances measured with the Hubble Space Telescope and find their agreement to be a vindication of our methodology. Our intent with this article is to spread awareness of the benefits and applicability of this Bayesian approach and encourage future PL relation investigations to consider employing this powerful analysis method.

The Fermi Large Area Telescope (LAT) reported the first definitive GeV detections of the binaries LS I +61\degree 303 and LS 5039 in the first year after its launch in June, 2008. These detections were unambiguous as a consequence of the reduced positional uncertainty and the detection of modulated gamma-ray emission on the corresponding orbital periods. An analysis of new data from the LAT, comprising 30 months of observations, identifies a change in the gamma-ray behavior of LS I +61\degree 303. An increase in flux is detected in March 2009 and a steady decline in the orbital flux modulation is observed. Significant emission up to 30GeV is detected by the LAT; prior datasets led to upper limits only. Contemporaneous TeV observations no longer detected the source, or found it -in one orbit- close to periastron, far from the phases at which the source previously appeared at TeV energies. The detailed numerical simulations and models that exist within the literature do not predict or explain many of these features now observed at GeV and TeV energies. New ideas and models are needed to fully explain and understand this behavior. A detailed phase-resolved analysis of the spectral characterization of LS I +61\degree 303 in the GeV regime ascribes a power law with an exponential cutoff spectrum along each analyzed portion of the system's orbit. The on-source exposure of LS 5039 is also substantially increased with respect to our prior publication. In this case, whereas the general gamma-ray properties remain consistent, the increased statistics of the current dataset allows for a deeper investigation of its orbital and spectral evolution.

With frequent flaring activity of its relativistic jets, Cygnus X-3 is one of the most active microquasars and is the only Galactic black hole candidate with confirmed high energy Gamma-ray emission, thanks to detections by Fermi/LAT and AGILE. In 2011, Cygnus X-3 was observed to transit to a soft X-ray state, which is known to be associated with high-energy Gamma-ray emission. We present the results of a multi-wavelength campaign covering a quenched state, when radio emission from Cygnus X-3 is at its weakest and the X-ray spectrum is very soft. A giant (~ 20 Jy) optically thin radio flare marks the end of the quenched state, accompanied by rising non-thermal hard X-rays. Fermi/LAT observations (E >100 MeV) reveal renewed Gamma-ray activity associated with this giant radio flare, suggesting a common origin for all non-thermal components. In addition, current observations unambiguously show that the Gamma-ray emission is not exclusively related to the rare giant radio flares. A 3-week period of Gamma-ray emission is also detected when Cygnus X-3 was weakly flaring in radio, right before transition to the radio quenched state. No Gamma rays are observed during the ~ one-month long quenched state, when the radio flux is weakest. Our results suggest transitions into and out of the ultrasoft X-ray (radio quenched) state trigger Gamma-ray emission, implying a connection to the accretion process, and also that the Gamma-ray activity is related to the level of radio flux (and possibly shock formation), strengthening the connection to the relativistic jets.

We report on optical spectroscopy of 165 Flat Spectrum Radio Quasars (FSRQs) in the Fermi 1LAC sample, which have helped allow a nearly complete study of this population. Fermi FSRQ show significant evidence for non-thermal emission even in the optical; the degree depends on the gamma-ray hardness. They also have smaller virial estimates of hole mass than the optical quasar sample. This appears to be largely due to a preferred (axial) view of the gamma-ray FSRQ and non-isotropic (H/R ~ 0.4) distribution of broad-line velocities. Even after correction for this bias, the Fermi FSRQ show higher mean Eddington ratios than the optical population. A comparison of optical spectral properties with Owens Valley Radio Observatory radio flare activity shows no strong correlation.

As the number of observed Gamma-Ray Bursts (GRBs) continues to grow, follow-up resources need to be used more efficiently in order to maximize science output from limited telescope time. As such, it is becoming increasingly important to rapidly identify bursts of interest as soon as possible after the event, before the afterglows fade beyond detectability. Studying the most distant (highest redshift) events, for instance, remains a primary goal for many in the field. Here we present our Random forest Automated Triage Estimator for GRB redshifts (RATE GRB-z) for rapid identification of high-redshift candidates using early-time metrics from the three telescopes onboard Swift. While the basic RATE methodology is generalizable to a number of resource allocation problems, here we demonstrate its utility for telescope-constrained follow-up efforts with the primary goal to identify and study high-z GRBs. For each new GRB, RATE GRB-z provides a recommendation - based on the available telescope time - of whether the event warrants additional follow-up resources. We train RATE GRB-z using a set consisting of 135 Swift bursts with known redshifts, only 18 of which are z > 4. Cross-validated performance metrics on this training data suggest that ~56% of high-z bursts can be captured from following up the top 20% of the ranked candidates, and ~84% of high-z bursts are identified after following up the top ~40% of candidates. We further use the method to rank 200+ Swift bursts with unknown redshifts according to their likelihood of being high-z.

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.