In 2010/3, the Large Area Telescope on board Fermi revealed a transient gamma-ray source, positionally coincident with the optical nova in the symbiotic binary, V407Cyg. This event marked the first discovery of gamma-ray emission from a nova. We aimed to obtain resolved radio imaging of the material involved in the nova event; to determine the ejecta geometry and advance velocity directly in the image plane; to constrain the physical conditions of the system. We observed the source with the EVN and the VLBA over 16 epochs, between 20 days and 6 months after the optical discovery. The source is initially very dim but it later shows a substantial increase in brightness and a resolved shell-like structure 40 to 90 days after the optical event. The shell has a projected elliptical shape and is asymmetric in brightness and spectral index, being brighter and characterised by a rising spectrum at the S-E edge. We determine a projected velocity of ~3500 km/s in the initial phase, and ~2100 km/s between day 20 and 91. We also found an emitting feature about 350 mas (940 AU) to the N-W, advancing at a projected velocity of ~700 km/s along the polar axis of the binary. The total flux density in the VLBI images is significantly lower than that previously reported at similar epochs and over much wider angular scales with the VLA. Optical spectra demonstrated that in 2010 we were viewing V407Cyg along the equatorial plane and from behind the Mira. Our radio observations image the bipolar flow of the ejecta perpendicular to the orbital plane, where deceleration is much lower than through the equatorial plane probed by the truncated profile of optical emission lines. The separated polar knot at 350 mas and the bipolar flow strictly resemble the similar arrangement seen in Hen 2-104. The observed ~700 km/s expansion constrains the launch-date of the polar knot around 2004. [Abridged]

With the Expanded Long Wavelength Array (ELWA) and pulsar binning techniques, we searched for off-pulse emission from PSR B0950+08 at 76 MHz. Previous studies suggest that off-pulse emission can be due to pulsar wind nebulae (PWNe) in younger pulsars. Other studies, such as that done by Basu et al. (2012), propose that in older pulsars this emission extends to some radius that is on the order of the light cylinder radius, and is magnetospheric in origin. Through imaging analysis we conclude that this older pulsar with a spin-down age of 17 Myr has a surrounding PWN, which is unexpected since as a pulsar ages its PWN spectrum is thought to shift from being synchrotron to inverse-Compton-scattering dominated. At 76 MHz, the average flux density of the off-pulse emission is 0.59 $\pm$ 0.16 Jy. The off-pulse emission from B0950+08 is $\sim$ 110 $\pm$ 17 arcseconds (0.14 $\pm$ 0.02 pc) in size, extending well-beyond the light cylinder diameter and ruling out a magnetospheric origin. Using data from our observation and the surveys VLSSr, TGSS, NVSS, FIRST, and VLASS, we have found that the spectral index for B0950+08 is about -1.36 $\pm$ 0.20, while the PWN's spectral index is steeper than -1.85 $\pm$ 0.45.

The Very Large Array Sky Survey (VLASS) is a synoptic, all-sky radio sky survey with a unique combination of high angular resolution ($\approx$2.5"), sensitivity (a 1$\sigma$ goal of 70 $\mu$Jy/beam in the coadded data), full linear Stokes polarimetry, time domain coverage, and wide bandwidth (2-4 GHz). The first observations began in September 2017, and observing for the survey will finish in 2024. VLASS will use approximately 5500 hours of time on the Karl G. Jansky Very Large Array (VLA) to cover the whole sky visible to the VLA (Declination $>-40^{\circ}$), a total of 33,885 deg$^2$. The data will be taken in three epochs to allow the discovery of variable and transient radio sources. The survey is designed to engage radio astronomy experts, multi-wavelength astronomers, and citizen scientists alike. By utilizing an "on the fly" interferometry mode, the observing overheads are much reduced compared to a conventional pointed survey. In this paper, we present the science case and observational strategy for the survey, and also results from early survey observations.

We present the results of the first transient survey from the Owens Valley Radio Observatory Long Wavelength Array (OVRO-LWA) using 31 hr of data, in which we place the most constraining limits on the instantaneous transient surface density at timescales of 13 s to a few minutes and at frequencies below 100 MHz. The OVRO-LWA is a dipole array that images the entire viewable hemisphere with 58 MHz of bandwidth from 27 to 84 MHz at 13 s cadence. No transients are detected above a 6.5$\sigma$ flux density limit of 10.5 Jy, implying an upper limit to the transient surface density of $2.5\times10^{-8}$ deg$^{-2}$ at the shortest timescales probed, which is orders of magnitude deeper than has been achieved at sub-100 MHz frequencies and comparable flux densities to date. The nondetection of transients in the OVRO-LWA survey, particularly at minutes-long timescales, allows us to place further constraints on the rate of the potential population of transients uncovered by Stewart et al. (2016). From their transient rate, we expect a detection of $8.4^{+31.8}_{-8.0}$ events, and the probability of our null detection is $1.9^{+644}_{-1.9}\times10^{-3}$, ruling out a transient rate $>1.4\times10^{-4} \text{days}^{-1} \text{deg}^{-2}$ with 95% confidence at a flux density limit of 18.1 Jy, under the assumption of a flat spectrum and wide bandwidth. We discuss the implications of our nondetection for this population and further constraints that can be made on the source spectral index, intrinsic emission bandwidth, and resulting luminosity distribution.

A successful ground array Radio Frequency (RF)-only self-trigger on 10 high-energy cosmic ray events is demonstrated with 256 dual-polarization antennas of the Owens Valley Radio Observatory Long Wavelength Array (OVRO-LWA). This RF-only capability is predicated on novel techniques for Radio Frequency Interference (RFI) identification and mitigation with an analysis efficiency of 45\% for shower-driven events with a Signal-to-noise ratio $\gtrsim$ 5 against the galactic background noise power of individual antennas. This technique enables more efficient detection of cosmic rays over a wider range of zenith angles than possible via triggers from in-situ particle detectors and can be easily adapted by neutrino experiments relying on RF-only detection. This paper discusses the system design, RFI characterization and mitigation techniques, and initial results from 10 cosmic ray events identified within a 40-hour observing window. A design for a future optimized commensal cosmic-ray detector for the OVRO-LWA is presented, as well as recommendations for developing a similar capability for other experiments -- these designs either reduce data-rate or increase sensitivity by an order of magnitude for many configurations of radio instruments.

The 21\,cm transition of neutral hydrogen is opening an observational window into the cosmic dawn of the universe---the epoch of first star formation. We use 28\,hr of data from the Owens Valley Radio Observatory Long Wavelength Array (OVRO-LWA) to place upper limits on the spatial power spectrum of 21\,cm emission at $z \approx 18.4$ ($\Delta_{21} \lesssim 10^4\,\text{mK}$), and within the absorption feature reported by the EDGES experiment (Bowman et al. 2018). In the process we demonstrate the first application of the double Karhunen-Lo\`{e}ve transform for foreground filtering, and diagnose the systematic errors that are currently limiting the measurement. We also provide an updated model for the angular power spectrum of low-frequency foreground emission measured from the northern hemisphere, which can be used to refine sensitivity forecasts for next-generation experiments.

We have carried out VLA imaging and a Fermi timing analysis of the 115 ms gamma-ray and radio pulsar PSR J0002+6216. We found that the pulsar lies at the apex of a narrowly collimated cometary-like 7 arcmin tail of non-thermal radio emission which we identify as a bow-shock pulsar wind nebula. The tail of the nebula points back toward the geometric center of the supernova remnant CTB 1 (G116.9+0.2) 28 arcmin away, at a position angle $\theta_\mu=113^\circ$. We measure a proper motion with 2.9$\sigma$ significance from a Fermi timing analysis giving $\mu$=115$\pm$33 mas yr$^{-1}$ and $\theta_\mu=121^\circ\pm{13}^\circ$, corresponding to a large transverse pulsar velocity of 1100 km s$^{-1}$ at a distance of 2 kpc. This proper motion is of the right magnitude and direction to support the claim that PSR J0002+6216 was born from the same supernova that produced CTB 1. We explore the implications for pulsar birth periods, asymmetric supernova explosions, and mechanisms for pulsar natal kick velocities.

Astrometry, the measurement of positions and motions of the stars, is one of the oldest disciplines in Astronomy, extending back at least as far as Hipparchus' discovery of the precession of Earth's axes in 190 BCE by comparing his catalog with those of his predecessors. Astrometry is fundamental to Astronomy, and critical to many aspects of Astrophysics and Geodesy. In order to understand our planet's and solar system's context within their surroundings, we must be able to to define, quantify, study, refine, and maintain an inertial frame of reference relative to which all positions and motions can be unambiguously and self-consistently described. It is only by using this inertial reference frame that we are able to disentangle our observations of the motions of celestial objects from our own complex path around our star, and its path through the galaxy, and the local group. Every aspect of each area outlined in the call for scientific frontiers in astronomy in the era of the 2020-2030 timeframe will depend on the quality of the inertial reference frame. In this white paper, we propose support for development of radio Very Long Baseline Interferometry (VLBI) capabilities, including the Next Generation Very Large Array (ngVLA), a radio astronomy observatory that will not only support development of a next generation reference frame of unprecedented accuracy, but that will also serve as a highly capable astronomical instrument in its own right. Much like its predecessors, the Very Long Baseline Array (VLBA) and other VLBI telescopes, the proposed ngVLA will provide the foundation for the next three decades for the fundamental reference frame, benefitting astronomy, astrophysics, and geodesy alike.

Multimessenger observations of the binary neutron star merger GW170817 have enabled the discovery of a diverse array of electromagnetic counterparts to compact binary mergers, including an unambiguous kilonova, a short gamma-ray burst, and a late-time radio jet. Beyond these counterparts, compact binary mergers are additionally predicted to be accompanied by prompt low-frequency radio emission. The successful observation of a prompt radio counterpart would be immensely valuable, but is made difficult by the short delay between the gravitational-wave and prompt electromagnetic signals as well as the poor localization of gravitational-wave sources. Here, we present the first search for prompt radio emission accompanying a gravitational-wave event, targeting the binary black hole merger GW170104 detected by the Advanced LIGO and Virgo gravitational-wave observatories during their second (O2) observing run. Using the Owens Valley Radio Observatory Long Wavelength Array (OVRO-LWA), we search a $\sim900\,\mathrm{deg}^2$ region for transient radio emission within approximately one hour of GW170104, obtaining an upper limit of $2.5\times10^{41}\,\mathrm{erg}\,\mathrm{s}^{-1}$ on its equivalent isotropic luminosity between 27-84 MHz. We additionally discuss plans to target binary neutron star mergers in Advanced LIGO and Virgo's upcoming O3 observing run.

The Large-Aperture Experiment to Detect the Dark Age (LEDA) was designed to detect the predicted O(100)mK sky-averaged absorption of the Cosmic Microwave Background by Hydrogen in the neutral pre- and intergalactic medium just after the cosmological Dark Age. The spectral signature would be associated with emergence of a diffuse Ly$\alpha$ background from starlight during 'Cosmic Dawn'. Recently, Bowman et al. (2018) have reported detection of this predicted absorption feature, with an unexpectedly large amplitude of 530 mK, centered at 78 MHz. Verification of this result by an independent experiment, such as LEDA, is pressing. In this paper, we detail design and characterization of the LEDA radiometer systems, and a first-generation pipeline that instantiates a signal path model. Sited at the Owens Valley Radio Observatory Long Wavelength Array, LEDA systems include the station correlator, five well-separated redundant dual polarization radiometers and backend electronics. The radiometers deliver a 30-85MHz band (16<z<34) and operate as part of the larger interferometric array, for purposes ultimately of in situ calibration. Here, we report on the LEDA system design, calibration approach, and progress in characterization as of January 2016. The LEDA systems are currently being modified to improve performance near 78 MHz in order to verify the purported absorption feature.

We have conducted the most sensitive low frequency (below 100 MHz) search to date for prompt, low-frequency radio emission associated with short-duration gamma-ray bursts (GRBs), using the Owens Valley Radio Observatory Long Wavelength Array (OVRO-LWA). The OVRO-LWA's nearly full-hemisphere field-of-view ($\sim20$,$000$ square degrees) allows us to search for low-frequency (sub-$100$ MHz) counterparts for a large sample of the subset of GRB events for which prompt radio emission has been predicted. Following the detection of short GRB 170112A by Swift, we used all-sky OVRO-LWA images spanning one hour prior to and two hours following the GRB event to search for a transient source coincident with the position of GRB 170112A. We detect no transient source, with our most constraining $1\sigma$ flux density limit of $650~\text{mJy}$ for frequencies spanning $27~\text{MHz}-84~\text{MHz}$. We place constraints on a number of models predicting prompt, low-frequency radio emission accompanying short GRBs and their potential binary neutron star merger progenitors, and place an upper limit of $L_\text{radio}/L_\gamma \lesssim 7\times10^{-16}$ on the fraction of energy released in the prompt radio emission. These observations serve as a pilot effort for a program targeting a wider sample of both short and long GRBs with the OVRO-LWA, including bursts with confirmed redshift measurements which are critical to placing the most constraining limits on prompt radio emission models, as well as a program for the follow-up of gravitational wave compact binary coalescence events detected by advanced LIGO and Virgo.

A host of new low-frequency radio telescopes seek to measure the 21-cm transition of neutral hydrogen from the early universe. These telescopes have the potential to directly probe star and galaxy formation at redshifts $20 \gtrsim z \gtrsim 7$, but are limited by the dynamic range they can achieve against foreground sources of low-frequency radio emission. Consequently, there is a growing demand for modern, high-fidelity maps of the sky at frequencies below 200 MHz for use in foreground modeling and removal. We describe a new widefield imaging technique for drift-scanning interferometers, Tikhonov-regularized $m$-mode analysis imaging. This technique constructs images of the entire sky in a single synthesis imaging step with exact treatment of widefield effects. We describe how the CLEAN algorithm can be adapted to deconvolve maps generated by $m$-mode analysis imaging. We demonstrate Tikhonov-regularized $m$-mode analysis imaging using the Owens Valley Long Wavelength Array (OVRO-LWA) by generating 8 new maps of the sky north of $\delta=-30^\circ$ with 15 arcmin angular resolution, at frequencies evenly spaced between 36.528 MHz and 73.152 MHz, and $\sim$800 mJy/beam thermal noise. These maps are a 10-fold improvement in angular resolution over existing full-sky maps at comparable frequencies, which have angular resolutions $\ge 2^\circ$. Each map is constructed exclusively from interferometric observations and does not represent the globally averaged sky brightness. Future improvements will incorporate total power radiometry, improved thermal noise, and improved angular resolution -- due to the planned expansion of the OVRO-LWA to 2.6 km baselines. These maps serve as a first step on the path to the use of more sophisticated foreground filters in 21-cm cosmology incorporating the measured angular and frequency structure of all foreground contaminants.

We present a community study exploring the low frequency (5 - 800 MHz) options and opportunities for the ngVLA project and its infrastructure. We describe a Next Generation LOw Band Observatory (ngLOBO) that will provide access to the low frequency sky in a commensal fashion, operating independently from the ngVLA, but leveraging common infrastructure. This approach provides continuous coverage through an aperture array (called ngLOBO-Low) below 150 MHz and by accessing the primary focus of the ngVLA antennas (called ngLOBO-High) above 150 MHz. ngLOBO preconditions include a) non-interference and b) low relative cost (<5%) with respect to ngVLA. ngLOBO has three primary scientific missions: (1) Radio Large Synoptic Survey Telescope (Radio-LSST): one naturally wide beam, commensal with ngVLA, will conduct a continuous synoptic survey of large swaths of the sky for both slow and fast transients; (2) This same commensal beam will provide complementary low frequency images of all ngVLA targets and their environment {\it when such data enhances their value}. (3) Independent beams from the ngLOBO-Low aperture array will conduct research in astrophysics, Earth science and space weather applications, engaging new communities and attracting independent resources. If ngVLA operates down to 2 GHz or lower, ngLOBO data will enhance ngVLA calibration and dynamic scheduling. Finally, non-variable field sources outside the ngVLA field of view can be harvested for serendipitous science, e.g. population studies for thermal and non-thermal continuum sources. The ngVLA will be a superb, high frequency instrument; ngLOBO will allow it to participate in the worldwide renaissance in low frequency science as well.

Radio astronomy observatories with high throughput back end instruments require real-time data processing. While computing hardware continues to advance rapidly, development of real-time processing pipelines remains difficult and time-consuming, which can limit scientific productivity. Motivated by this, we have developed Bifrost: an open-source software framework for rapid pipeline development. Bifrost combines a high-level Python interface with highly efficient reconfigurable data transport and a library of computing blocks for CPU and GPU processing. The framework is generalizable, but initially it emphasizes the needs of high-throughput radio astronomy pipelines, such as the ability to process data buffers as if they were continuous streams, the capacity to partition processing into distinct data sequences (e.g., separate observations), and the ability to extract specific intervals from buffered data. Computing blocks in the library are designed for applications such as interferometry, pulsar dedispersion and timing, and transient search pipelines. We describe the design and implementation of the Bifrost framework and demonstrate its use as the backbone in the correlation and beamforming back end of the Long Wavelength Array station in the Sevilleta National Wildlife Refuge, NM.

We present a survey of the radio sky accessible from the first station of the Long Wavelength Array (LWA1). Images are presented at nine frequencies between 35 and 80 MHz with spatial resolutions ranging from $4.7^\circ$ to $2.0^\circ$, respectively. The maps cover the sky north of a declination of $-40^\circ$ and represent the most modern systematic survey of the diffuse Galactic emission within this frequency range. We also combine our survey with other low frequency sky maps to create an updated model of the low frequency sky. Due to the low frequencies probed by our survey, the updated model better accounts for the effects of free-free absorption from Galactic ionized Hydrogen. A longer term motivation behind this survey is to understand the foreground emission that obscures the redshifted 21 cm transition of neutral hydrogen from the cosmic dark ages ($z>10$) and, at higher frequencies, the epoch of reionization ($z>6$).

The third Fermi Large Area Telescope (LAT) $\gamma$-ray source catalog (3FGL) contains over 1000 objects for which there is no known counterpart at other wavelengths. The physical origin of the $\gamma$-ray emission of those objects is unknown. Such objects are commonly referred to as unassociated and mostly do not exhibit significant $\gamma$-ray flux variability. We performed a survey of all unassociated $\gamma$-ray sources found in 3FGL using the Australia Telescope Compact Array and Very Large Array in the range of 4.0-10.0 GHz. We found 2097 radio candidates for association with $\gamma$-ray sources. The follow-up with very long baseline interferometry for a subset of those candidates yielded 142 new AGN associations with $\gamma$-ray sources, provided alternative associations for 7 objects, and improved positions for another 144 known associations to the milliarcsecond level of accuracy. In addition, for 245 unassociated $\gamma$-ray sources we did not find a single compact radio source above 2 mJy within 3$\sigma$ of their $\gamma$-ray localization. A significant fraction of these empty fields, 39%, are located away from the galactic plane. We also found 36 extended radio sources that are candidates for association with a corresponding $\gamma$-ray object, 19 of which are most likely supernova remnants or HII regions, whereas 17 could be radio galaxies.

Blazars are highly variable active galactic nuclei which emit radiation at all wavelengths from radio to gamma-rays. Polarized radiation from blazars is one key piece of evidence for synchrotron radiation at low energies and it also varies dramatically. The polarization of blazars is of interest for understanding the origin, confinement, and propagation of jets. However, even though numerous measurements have been performed, the mechanisms behind jet creation, composition and variability are still debated. We performed simultaneous gamma-ray and optical photopolarimetry observations of 45 blazars between Jul. 2008 and Dec. 2014 to investigate the mechanisms of variability and search for a basic relation between the several subclasses of blazars. We identify a correlation between the maximum degree of optical linear polarization and the gamma-ray luminosity or the ratio of gamma-ray to optical fluxes. Since the maximum polarization degree depends on the condition of the magnetic field (chaotic or ordered), this result implies a systematic difference in the intrinsic alignment of magnetic fields in pc-scale relativistic jets between different types blazars (FSRQs vs. BL Lacs), and consequently between different types of radio galaxies (FR Is vs. FR IIs).

Rotating Radio Transients (RRATs) are a subclass of pulsars first identified in 2006 that are detected only in searches for single pulses and not through their time averaged emission. Here, we present the results of observations of 19 RRATs using the first station of the Long Wavelength Array (LWA1) at frequencies between 30 MHz and 88 MHz. The RRATs observed here were first detected in higher frequency pulsar surveys. Of the 19 RRATs observed, 2 sources were detected and their dispersion measures, periods, pulse profiles, and flux densities are reported and compared to previous higher frequency measurements. We find a low detection rate (11%), which could be a combination of the lower sensitivity of LWA1 compared to the higher frequency telescopes, and the result of scattering by the interstellar medium or a spectral turnover.

Utilizing the all-sky imaging capabilities of the LWA1 radio telescope along with a host of all-sky optical cameras, we have now observed 44 optical meteor counterparts to radio afterglows. Combining these observations we have determined the geographic positions of all 44 afterglows. Comparing the number of radio detections as a function of altitude above sea level to the number of expected bright meteors we find a strong altitudinal dependence characterized by a cutoff below $\sim$ 90 km, below which no radio emission occurs, despite the fact that many of the observed optical meteors penetrated well below this altitude. This cutoff suggests that wave damping from electron collisions is an important factor for the evolution of radio afterglows, which agrees with the hypothesis that the emission is the result of electron plasma wave emission.

We present a large survey of giant pulses from the Crab Pulsar as observed with the first station of the Long Wavelength Array. Automated methods for detecting giant pulses at low frequencies where scattering becomes prevalent are also explored. More than 1400 pulses were detected across four frequency bands between 20 - 84 MHz over a seven month period beginning in 2013, with additional followup observations in late 2014 and early 2015. A handful of these pulses were detected simultaneously across all four frequency bands. We examine pulse characteristics, including pulse broadening and power law indices for amplitude distributions. We find that the flux density increases toward shorter wavelengths, consistent with a spectral turnover at 100 MHz. Our observations uniquely span multiple scattering epochs, manifesting as a notable trend in the number of detections per observation. These results are characteristic of the variable interface between the synchrotron nebula and the surrounding interstellar medium.

Using the narrowband all-sky imager mode of the LWA1 we have now detected 30 transients at 25.6 MHz, 1 at 34 MHz, and 93 at 38.0 MHz. While we have only optically confirmed that 37 of these events are radio afterglows from meteors, evidence suggests that most, if not all, are. Using the beam-forming mode of the LWA1 we have also captured the broadband spectra between 22.0 and 55.0 MHz of four events. We compare the smooth, spectral components of these four events and fit the frequency dependent flux density to a power law, and find that the spectral index is time variable, with the spectrum steepening over time for each meteor afterglow. Using these spectral indices along with the narrow band flux density measurements of the 123 events at 25.6 and 38 MHz, we predict the expected flux densities and rates for meteor afterglows potentially observable by other low frequency radio telescopes.

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

The birth of the first luminous sources and the ensuing epoch of reionization are best studied via the redshifted 21-cm emission line, the signature of the first two imprinting the last. In this work we present a fully-Bayesian method, \textsc{hibayes}, for extracting the faint, global (sky-averaged) 21-cm signal from the much brighter foreground emission. We show that a simplified (but plausible), Gaussian model of the 21-cm emission from the Cosmic Dawn epoch ($15 \lesssim z \lesssim 30$), parameterized by an amplitude $A_{\rm HI}$, a frequency peak $\nu_{\rm HI}$ and a width $\sigma_{\rm HI}$, can be extracted even in the presence of a structured foreground frequency spectrum (parameterized as a $7^{\rm th}$-order polynomial), provided sufficient signal-to-noise (400~hours of observation with a single dipole). We apply our method to an early, 19-minute long observation from the Large aperture Experiment to detect the Dark Ages, constraining the 21-cm signal amplitude and width to be $-890 < A_{\rm HI} < 0$ mK and $\sigma_{\rm HI} > 6.5$ MHz (corresponding to $\Delta z > 1.9$ at redshift $z \simeq 20$) respectively at the 95-per-cent confidence level in the range $13.2 < z < 27.4$ ($100 > \nu > 50$ MHz).

We present dynamic spectra from the LWA1 telescope of two large meteors (fireballs) observed to emit between 37 and 54 MHz. These spectra show the first ever recorded broadband measurements of this newly discovered VHF emission. The spectra show that the emission is smooth and steep, getting very bright at lower frequencies. We suggest that this signal is possibly emission of Langmuir waves and that these waves could be excited by a weak electron beam within the trail. The spectra of one fireball displays broadband temporal frequency sweeps. We suggest that these sweeps are evidence of individual expanding clumps of emitting plasma. While some of these proposed clumps may have formed at the very beginning of the fireball event, others must have formed seconds after the initial event.

The third catalog of active galactic nuclei (AGNs) detected by the Fermi-LAT (3LAC) is presented. It is based on the third Fermi-LAT catalog (3FGL) of sources detected between 100 MeV and 300 GeV with a Test Statistic (TS) greater than 25, between 2008 August 4 and 2012 July 31. The 3LAC includes 1591 AGNs located at high Galactic latitudes (|b|>10{\deg}), a 71% increase over the second catalog based on 2 years of data. There are 28 duplicate associations, thus 1563 of the 2192 high-latitude gamma-ray sources of the 3FGL catalog are AGNs. Most of them (98%) are blazars. About half of the newly detected blazars are of unknown type, i.e., they lack spectroscopic information of sufficient quality to determine the strength of their emission lines. Based on their gamma-ray spectral properties, these sources are evenly split between flat-spectrum radio quasars (FSRQs) and BL~Lacs. The most abundant detected BL~Lacs are of the high-synchrotron-peaked (HSP) type. About 50% of the BL~Lacs have no measured redshifts. A few new rare outliers (HSP-FSRQs and high-luminosity HSP BL~Lacs) are reported. The general properties of the 3LAC sample confirm previous findings from earlier catalogs. The fraction of 3LAC blazars in the total population of blazars listed in BZCAT remains non-negligible even at the faint ends of the BZCAT-blazar radio, optical and X-ray flux distributions, which is a clue that even the faintest known blazars could eventually shine in gamma rays at LAT-detection levels. The energy-flux distributions of the different blazar populations are in good agreement with extrapolation from earlier catalogs.

We present initial pulsar results from the first station of the Long Wavelength Array (LWA1) obtained during the commissioning period of LWA1 and early science results. We present detections of periodic emission from 44 previously known pulsars, including 3 millisecond pulsars (MSPs). The effects of the interstellar medium on pulsar emission are significantly enhanced at the low frequencies of the LWA1 band (10--88 MHz), making LWA1 a very sensitive instrument for characterizing changes in dispersion measures (DM) and other effects from the interstellar medium. Pulsars also often have significant evolution in their pulse profile at low frequency and a break in their spectral index. We report DM measurements for 44 pulsars, mean flux density measurements for 36 pulsars, and multi-frequency component spacing and widths for 15 pulsars with more than one profile component. For 27 pulsars, we report spectral index measurements within our frequency range. We also introduce the LWA1 Pulsar Data Archive, which stores reduced data products from LWA1 pulsar observations. Reduced data products for the observations presented here can be found on the archive. Reduced data products from future LWA1 pulsar observations will also be made available through the archive.

We present a description of the Prototype All-Sky Imager (PASI), a backend correlator and imager of the first station of the Long Wavelength Array (LWA1). PASI cross-correlates a live stream of 260 dual-polarization dipole antennas of the LWA1, creates all-sky images, and uploads them to the LWA-TV website in near real-time. PASI has recorded over 13,000 hours of all-sky images at frequencies between 10 and 88 MHz creating opportunities for new research and discoveries. We also report rate density and pulse energy density limits on transients at 38, 52, and 74 MHz, for pulse widths of 5 s. We limit transients at those frequencies with pulse energy densities of $>2.7\times 10^{-23}$, $>1.1\times 10^{-23}$, and $>2.8\times 10^{-23}$ J m$^{-2}$ Hz$^{-1}$ to have rate densities $<1.2\times10^{-4}$, $<5.6\times10^{-4}$, and $<7.2\times10^{-4}$ yr$^{-1}$ deg$^{-2}$

We report the detection of giant pulse emission from PSR B0950+08 in 24 hours of observations made at 39.4 MHz, with a bandwidth of 16 MHz, using the first station of the Long Wavelength Array, LWA1. We detected 119 giant pulses from PSR B0950+08 (at its dispersion measure), which we define as having SNRs at least 10 times larger than for the mean pulse in our data set. These 119 pulses are 0.035% of the total number of pulse periods in the 24 hours of observations. The rate of giant pulses is about 5.0 per hour. The cumulative distribution of pulse strength $S$ is a steep power law, $N(>S)\propto S^{-4.7}$, but much less steep than would be expected if we were observing the tail of a Gaussian distribution of normal pulses. We detected no other transient pulses in a dispersion measure range from 1 to 90 pc cm$^{-3}$, in the beam tracking PSR B0950+08. The giant pulses have a narrower temporal width than the mean pulse (17.8 ms, on average, vs. 30.5 ms). The pulse widths are consistent with a previously observed weak dependence on observing frequency, which may be indicative of a deviation from a Kolmogorov spectrum of electron density irregularities along the line of sight. The rate and strength of these giant pulses is less than has been observed at $\sim$100 MHz. Additionally, the mean (normal) pulse flux density we observed is less than at $\sim$100 MHz. These results suggest this pulsar is weaker and produces less frequent giant pulses at 39 MHz than at 100 MHz.

We present the results of an all-sky radio survey between 5 and 9 GHz of the fields surrounding all unassociated gamma-ray objects listed in the Fermi Large Area Telescope Second Source Catalog (2FGL). The goal of these observations is to find all new gamma-ray AGN associations with radio sources >10 mJy at 8 GHz. We observed with the Very Large Array and the Australia Telescope Compact Array the areas around unassociated sources, providing localizations of weak radio point sources found in 2FGL fields at arcmin scales. Then we followed-up a subset of those with the Very Long Baseline and the Long Baseline Arrays to confirm detections of radio emission on parsec-scales. We quantified association probabilities based on known statistics of source counts and assuming a uniform distribution of background sources. In total we found 865 radio sources at arcsec scales as candidates for association and detected 95 of 170 selected for follow-up observations at milliarcsecond resolution. Based on this we obtained firm associations for 76 previously unknown gamma-ray AGN. Comparison of these new AGN associations with the predictions from using the WISE color-color diagram shows that half of the associations are missed. We found that 129 out of 588 observed gamma-ray sources at arcmin scales not a single radio continuum source was detected above our sensitivity limit within the 3-sigma gamma-ray localization. These "empty" fields were found to be particularly concentrated at low Galactic latitudes. The nature of these Galactic gamma-ray emitters is not yet determined.

A "large-N" correlator that makes use of Field Programmable Gate Arrays and Graphics Processing Units has been deployed as the digital signal processing system for the Long Wavelength Array station at Owens Valley Radio Observatory (LWA-OV), to enable the Large Aperture Experiment to Detect the Dark Ages (LEDA). The system samples a ~100MHz baseband and processes signals from 512 antennas (256 dual polarization) over a ~58MHz instantaneous sub-band, achieving 16.8Tops/s and 0.236 Tbit/s throughput in a 9kW envelope and single rack footprint. The output data rate is 260MB/s for 9 second time averaging of cross-power and 1 second averaging of total-power data. At deployment, the LWA-OV correlator was the largest in production in terms of N and is the third largest in terms of complex multiply accumulations, after the Very Large Array and Atacama Large Millimeter Array. The correlator's comparatively fast development time and low cost establish a practical foundation for the scalability of a modular, heterogeneous, computing architecture.

New observations of Jupiter's decametric radio emissions have been made with the Long Wavelength Array Station 1 (LWA1) which is capable of making high quality observations as low as 11 MHz. Full Stokes parameters were determined for bandwidths of 16 MHz. Here we present the first LWA1 results for the study of six Io-related events at temporal resolutions as fine as 0.25 ms. LWA1 data show excellent spectral detail in Jovian DAM such as simultaneous left hand circular (LHC) and right hand circular (RHC) polarized Io-related arcs and source envelopes, modulation lane features, S-bursts structures, narrow band N-events, and interactions between S-bursts and N-events. The sensitivity of the LWA1 combined with the low radio frequency interference environment allow us to trace the start of the LHC Io-C source region to much earlier CMLIII than typically found in the literature. We find the Io-C starts as early as CMLIII = 230 degrees at frequencies near 11 MHz. This early start of the Io-C emission may be valuable for refining models of the emission mechanism. We also detect modulation lane structures that appear continuous across LHC and RHC emissions, suggesting that both polarizations may originate from the same hemisphere of Jupiter. We present a study of rare S-bursts detected during an Io-D event and show drift rates are consistent with those from other Io-related sources. Finally, S-N burst events are seen in high spectral and temporal resolution and our data strongly support the co-spatial origins of these events.

We present the findings from the Prototype All-Sky Imager (PASI), a backend correlator of the first station of the Long Wavelength Array (LWA1), which has recorded over 11,000 hours of all-sky images at frequencies between 25 and 75 MHz. In a search of this data for radio transients, we have found 49 long (10s of seconds) duration transients. Ten of these transients correlate both spatially and temporally with large meteors (fireballs), and their signatures suggest that fireballs emit a previously undiscovered low frequency, non-thermal pulse. This emission provides a new probe into the physics of meteors and identifies a new form of naturally occurring radio transient foreground.

As a backend to the first station of the Long Wavelength Array (LWA1) the Prototype All Sky Imager (PASI) has been imaging the sky $>$ -26$^{\circ}$ declination during 34 Gamma Ray Bursts (GRBs) between January 2012 and May 2013. Using this data we were able to put the most stringent limits to date on prompt low frequency emission from GRBs. While our limits depend on the zenith angle of the observed GRB, we estimate a 1$\sigma$ RMS sensitivity of 68, 65 and 70 Jy for 5 second integrations at 37.9, 52.0, and 74.0 MHz at zenith. These limits are relevant for pulses $\geq$ 5 s and are limited by dispersion smearing. For pulses of length 5 s we are limited to dispersion measures ($DM$s) $\leq$ 220, 570, and 1,600 pc cm$^{-3}$ for the frequencies above. For pulses lasting longer than 5s, the $DM$ limits increase linearly with the duration of the pulse. We also report two interesting transients, which are, as of yet, of unknown origin, and are not coincident with any known GRBs. For general transients, we give rate density limits of $\leq$ $7.5\times10^{-3}$, $2.9\times10^{-2}$, and $1.4\times10^{-2}$ yr$^{-1}$ deg$^{-2}$ with pulse energy densities $>1.3\times 10^{-22}$, $1.1\times 10^{-22}$, and $1.4\times 10^{-22}$ J m$^{-2}$ Hz$^{-1}$ and pulse widths of 5 s at the frequencies given above.

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.

Radio astronomical imaging arrays comprising large numbers of antennas, O(10^2-10^3) have posed a signal processing challenge because of the required O(N^2) cross correlation of signals from each antenna and requisite signal routing. This motivated the implementation of a Packetized Correlator architecture that applies Field Programmable Gate Arrays (FPGAs) to the O(N) "F-stage" transforming time domain to frequency domain data, and Graphics Processing Units (GPUs) to the O(N^2) "X-stage" performing an outer product among spectra for each antenna. The design is readily scalable to at least O(10^3) antennas. Fringes, visibility amplitudes and sky image results obtained during field testing are presented.

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

We present flux density measurements and pulse profiles for the millisecond pulsar PSR J2145-0750 spanning 37 to 81 MHz using data obtained from the first station of the Long Wavelength Array. These measurements represent the lowest frequency detection of pulsed emission from a millisecond pulsar to date. We find that the pulse profile is similar to that observed at 102 MHz. We also find that the flux density spectrum between ~40 MHz to 5 GHz is suggestive of a break and may be better fit by a model that includes spectral curvature with a rollover around 730 MHz rather than a single power law.

We report results of the first phase of observations with the Australia Telescope Compact Array (ATCA) at 5 and 9 GHz of the fields around 411 gamma-ray sources with declinations < +10 deg detected by Fermi but marked as unassociated in the 2FGL catalogue. We have detected 424 sources with flux densities in a range of 2 mJy to 6 Jy that lie within the 99 per cent localisation uncertainty of 283 gamma-ray sources. Of these, 146 objects were detected in both the 5 and 9 GHz bands. We found 84 sources in our sample with a spectral index flatter than -0.5. The majority of detected sources are weaker than 100 mJy and for this reason were not found in previous surveys. Approximately 1/3 of our sample, 128 objects, have the probability of being associated by more than 10 times than the probability of being a background source found in the vicinity of a gamma-ray object by chance. We present the catalogue of positions of these sources, estimates of their flux densities and spectral indices where available.

We report on multiwavelength observations of the blazar PKS 0537-441 (z = 0.896) obtained from microwaves through gamma rays by SMA, REM, ATOM, Swift and Fermi during 2008 August-2010 April. Strong variability has been observed in gamma rays, with two major flaring episodes (2009 July and 2010 March) and a harder-when-brighter behaviour, quite common for FSRQs and low-synchrotron-peaked BL Lacs, in 2010 March. In the same way the SED of the source cannot be modelled by a simple synchrotron self-Compton model, as opposed to many BL Lacs, but the addition of an external Compton component of seed photons from a dust torus is needed. The 230 GHz light curve showed an increase simultaneous with the gamma-ray one, indicating co-spatiality of the mm and gamma-ray emission region likely at large distance from the central engine. The low, average, and high activity SED of the source could be fit changing only the electron distribution parameters, but two breaks in the electron distribution are necessary. The ensuing extra spectral break, located at NIR-optical frequencies, together with that in gamma rays seem to indicate a common origin, most likely due to an intrinsic feature in the underlying electron distribution. An overall correlation between the gamma-ray band with the R-band and K-band has been observed with no significant time lag. On the other hand, when inspecting the light curves on short time scales some differences are evident. In particular, flaring activity has been detected in NIR and optical bands with no evident gamma-ray counterparts in 2009 September and November. Moderate variability has been observed in X-rays with no correlation between flux and photon index. An increase of the detected X-ray flux with no counter part at the other wavelengths has been observed in 2008 October, suggesting once more a complex correlation between the emission at different energy bands.

LWA1 is a new radio telescope operating in the frequency range 10-88 MHz, located in central New Mexico. The telescope consists of 258 pairs of dipole-type antennas whose outputs are individually digitized and formed into beams. Simultaneously, signals from all dipoles can be recorded using one of the instrument's "all dipoles" modes, facilitating all-sky imaging. Notable features of the instrument include high intrinsic sensitivity (about 6 kJy zenith system equivalent flux density), large instantaneous bandwidth (up to 78 MHz), and 4 independently-steerable beams utilizing digital "true time delay" beamforming. This paper summarizes the design of LWA1 and its performance as determined in commissioning experiments. We describe the method currently in use for array calibration, and report on measurements of sensitivity and beamwidth.

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.

The parsec-scale radio properties of 232 active galactic nuclei (AGNs), most of which are blazars, detected by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope have been observed contemporaneously by the Very Long Baseline Array (VLBA) at 5 GHz. Data from both the first 11 months (1FGL) and the first 2 years (2FGL) of the Fermi mission were used to investigate these sources' gamma-ray properties. We use the ratio of the gamma-ray to radio luminosity as a measure of gamma-ray loudness. We investigate the relationship of several radio properties to gamma-ray loudness and to the synchrotron peak frequency. There is a tentative correlation between gamma-ray loudness and synchrotron peak frequency for BL Lac objects in both 1FGL and 2FGL, and for flat-spectrum radio quasars (FSRQs) in 2FGL. We find that the apparent opening angle tentatively correlates with gamma-ray loudness for FSRQs, but only when we use the 2FGL data. We also find that the total VLBA flux density correlates with the synchrotron peak frequency for BL Lac objects and FSRQs. The core brightness temperature also correlates with synchrotron peak frequency, but only for the BL Lac objects. The low-synchrotron peaked (LSP) BL Lac object sample shows indications of contamination by FSRQs which happen to have undetectable emission lines. There is evidence that the LSP BL Lac objects are more strongly beamed than the rest of the BL Lac object population.

The first station of the Long Wavelength Array (LWA1) was completed in April 2011 and is currently performing observations resulting from its first call for proposals in addition to a continuing program of commissioning and characterization observations. The instrument consists of 258 dual-polarization dipoles, which are digitized and combined into beams. Four independently-steerable dual-polarization beams are available, each with two "tunings" of 16 MHz bandwidth that can be independently tuned to any frequency between 10 MHz and 88 MHz. The system equivalent flux density for zenith pointing is ~3 kJy and is approximately independent of frequency; this corresponds to a sensitivity of ~5 Jy/beam (5sigma, 1 s); making it one of the most sensitive meter-wavelength radio telescopes. LWA1 also has two "transient buffer" modes which allow coherent recording from all dipoles simultaneously, providing instantaneous all-sky field of view. LWA1 provides versatile and unique new capabilities for Galactic science, pulsar science, solar and planetary science, space weather, cosmology, and searches for astrophysical transients. Results from LWA1 will detect or tightly constrain the presence of hot Jupiters within 50 parsecs of Earth. LWA1 will provide excellent resolution in frequency and in time to examine phenomena such as solar bursts, and pulsars over a 4:1 frequency range that includes the poorly understood turnover and steep-spectrum regimes. Observations to date have proven LWA1's potential for pulsar observing, and just a few seconds with the completed 256-dipole LWA1 provide the most sensitive images of the sky at 23 MHz obtained yet. We are operating LWA1 as an open skies radio observatory, offering ~2000 beam-hours per year to the general community.

We report the development of a semi-automatic pipeline for the calibration of 86 GHz full-polarization observations performed with the Global Millimeter-VLBI array (GMVA) and describe the calibration strategy followed in the data reduction. Our calibration pipeline involves non-standard procedures, since VLBI polarimetry at frequencies above 43 GHz is not yet well established. We also present, for the first time, a full-polarization global-VLBI image at 86 GHz (source 3C 345), as an example of the final product of our calibration pipeline, and discuss the effect of instrumental limitations on the fidelity of the polarization images. Our calibration strategy is not exclusive for the GMVA, and could be applied on other VLBI arrays at millimeter wavelengths. The use of this pipeline will allow GMVA observers to get fully-calibrated datasets shortly after the data correlation.

Aims: 3C345 was recently identified as a gamma-ray emitter, based on the first 20 months of Fermi-LAT data and optical monitoring. In this paper, a connection between the gamma-ray and optical variability of 3C345 and properties of its parsec-scale radio emission is investigated. Methods: The Fermi-LAT data of 3C345, covering an energy range of 0.1-300 GeV, were combined with 32 Very Long Baseline Array observations of the object made at 43.2 GHz in the period of January 2008 - March 2010. Results: The VLBA data reveal morphology and kinematics of the flow on scales of up to ~5 milliarcseconds (mas; deprojected linear distances of 380 parsecs). The brightness temperature, T_b(r), measured along the jet first decreases with distance proportional to r^-(0.95 +/-0.69) and later exhibits a break at ~0.3 mas, with T_b(r) proportional to r^-(4.11 +/-0.85) at larger separations. Variations of the gamma-ray, optical and parsec-scale radio emission show a similar long-term trend persistent during the entire VLBA monitoring period. The gamma-ray and optical variations on shorter time scales are related to structural changes in the jet on scales of ~0.3 mas (~23 parsecs, deprojected), with the gamma-ray and optical flares possibly related to the evolution of four distinct superluminal components identified in the flow. Conclusions: The observations indicate that both the quiescent and flaring components of the gamma-ray emission are produced in a region of the jet of ~23 pc in extent. This region may mark the Compton-loss dominated zone of the flow and its large extent may favor the synchrotron self-Compton mechanism for gamma-ray production in the relativistic jet of the quasar 3C345.

For more than 15 years, since the days of the Energetic Gamma-Ray Experiment Telescope (EGRET) on board the Compton Gamma-Ray Observatory (CGRO; 1991-2000), it has remained an open question why the prominent blazar 3C 345 was not reliably detected at gamma-ray energies <=20 MeV. Recently a bright gamma-ray source (0FGL J1641.4+3939/1FGL J1642.5+3947), potentially associated with 3C 345, was detected by the Large Area Telescope (LAT) on Fermi. Multiwavelength observations from radio bands to X-rays (mainly GASP-WEBT and Swift) of possible counterparts (3C 345, NRAO 512, B3 1640+396) were combined with 20 months of Fermi-LAT monitoring data (August 2008 - April 2010) to associate and identify the dominating gamma-ray emitting counterpart of 1FGL J1642.5+3947. The source 3C 345 is identified as the main contributor for this gamma-ray emitting region. However, after November 2009 (15 months), a significant excess of photons from the nearby quasar NRAO 512 started to contribute and thereafter was detected with increasing gamma-ray activity, possibly adding flux to 1FGL J1642.5+3947. For the same time period and during the summer of 2010, an increase of radio, optical and X-ray activity of NRAO 512 was observed. No gamma-ray emission from B3 1640+396 was detected.

We report on the discovery and observations of the extremely luminous optical transient CSS100217:102913+404220 (CSS100217 hereafter). Spectroscopic observations show this transient was coincident with a galaxy at redshift z=0.147, and reached an apparent magnitude of V ~ 16.3. After correcting for foreground Galactic extinction we determine the absolute magnitude to be M_V =-22.7 approximately 45 days after maximum light. Based on our unfiltered optical photometry the peak optical emission was L = 1.3 x 10^45 erg s^-1, and over a period of 287 rest-frame days had an integrated bolometric luminosity of 1.2 x 10^52 erg. Analysis of the pre-outburst SDSS spectrum of the source shows features consistent with a Narrow-line Seyfert1 (NLS1) galaxy. High-resolution HST and Keck followup observations show the event occurred within 150pc of nucleus of the galaxy, suggesting a possible link to the active nuclear region. However, the rapid outburst along with photometric and spectroscopic evolution are much more consistent with a luminous supernova. Line diagnostics suggest that the host galaxy is undergoing significant star formation. We use extensive follow-up of the event along with archival CSS and SDSS data to investigate the three most likely sources of such an event; 1) an extremely luminous supernova; 2) the tidal disruption of a star by the massive nuclear black hole; 3) variability of the central AGN. We find that CSS100217 was likely an extremely luminous type IIn supernova that occurred within range of the narrow-line region of an AGN. We discuss how similar events may have been missed in past supernova surveys because of confusion with AGN activity.

In this paper we report on the two-year-long Fermi-LAT observation of the peculiar blazar 4C +21.35 (PKS 1222+216). This source was in a quiescent state from the start of science operations of the Fermi Gamma-ray Space Telescope in 2008 August until 2009 September, and then became more active, with gradually increasing flux and some moderately-bright flares. In 2010 April and June, 4C +21.35 underwent a very strong GeV outburst composed of several major flares characterized by rise and decay timescales of the order of a day. During the outburst, the GeV spectra of 4C +21.35 displayed a broken power-law form with spectral breaks observed near 1-3 GeV photon energies. We demonstrate that, at least during the major flares, the jet in 4C +21.35 carried a total kinetic luminosity comparable to the total accretion power available to feed the outflow. We also discuss the origin of the break observed in the flaring spectra of 4C +21.35. We show that, in principle, a model involving annihilation of the GeV photons on the He II Lyman recombination continuum and line emission of "broad line region" clouds may account for such. However, we also discuss the additional constraint provided by the detection of 4C +21.35 at 0.07-0.4 TeV energies by the MAGIC telescope, which coincided with one of the GeV flares of the source. We argue that there are reasons to believe that the $lesssim$,TeV emission of 4C +21.35 (as well as the GeV emission of the source, if co-spatial), is not likely to be produced inside the broad line region zone of highest ionization ($sim 10^{17}$,cm from the nucleus), but instead originates further away from the active center, namely around the characteristic scale of the hot dusty torus surrounding the 4C +21.35 nucleus ($sim 10^{19}$,cm).

3C 345 is one of the archetypical active galactic nuclei, showing structural and flux variability on parsec scales near a compact unresolved radio core. During the last 2 years, the source has been undergoing a period of high activity visible in the broad spectral range, from radio through high-energy bands. We have been monitoring parsec-scale radio emission in 3C 345 during this period at monthly intervals, using the VLBA at 15, 24, and 43 GHz. Our radio observations are compared with gamma-ray emission detected by Fermi-LAT in the region including 3C 345 (1FGL J1642.5+3947). Three distinct gamma-ray events observed in this region are associated with the propagation of relativistic plasma condensations inside the radio jet of 3C 345. We report on evidence for the gamma-rays to be produced in a region of the jet of up to 40 pc (de-projected) in extent. This suggests the synchrotron self-Compton process as the most likely mechanism for production of gamma-rays in the source.

Single-zone synchrotron self-Compton and external Compton models are widely used to explain broad-band Spectral Energy Distributions (SEDs) of blazars from infrared to gamma-rays. These models bear obvious similarities to the homogeneous synchrotron cloud model which is often applied to explain radio emission from individual components of parsec-scale radio jets. The parsec-scale core, typically the brightest and most compact feature of blazar radio jet, could be the source of high-energy emission. We report on ongoing work to test this hypothesis by deriving the physical properties of parsec-scale radio emitting regions of twenty bright Fermi blazars using dedicated 5-43 GHz VLBA observations and comparing these parameters to results of SED modeling.

The 16th magnitude quasar 3C 345 (redshift z=0.5928) shows structural and emission variability on parsec scales around a compact unresolved radio core. For the last three decades it has been closely monitored with very long baseline interferometry (VLBI), yielding a wealth of information about the physics of relativistic outflows and dynamics of the central regions in AGN. We present here preliminary results for the long-term jet evolution, based on the 15 GHz monitoring data collected by the MOJAVE survey and various other groups over the last ~14 years and combined with data from earlier VLBI observations of 3C 345 which started in 1979. We discuss the trajectories, kinematics, and flux density evolution of enhanced emission regions embedded in the jet and present evidence for geometrical (e.g. precession) and physical (e.g. relativistic shocks and plasma instability) factors determining the morphology and dynamics of relativistic flows on parsec scales.

SN 2001em is a peculiar supernova, originally classified as Type Ib/c. About two years after the SN it was detected in the radio, showing a rising radio flux with an optically thin spectral slope, and it also displayed a large X-ray luminosity (~10^{41} erg/s). Thus it was suspected to harbor a decelerating (by then, mildly) relativistic jet pointing away from us. About 3 years after its discovery the optical spectrum of SN 2001em showed a broad H-alpha line, and it was therefore reclassified as Type IIn. Here we constrain its proper motion and expansion velocity by analyzing four epochs of VLBI observations, extending out to 5.4 years after the SN. The supernova is still unresolved 5.4 years after the explosion. For the proper motion we obtain (23,000 +/- 30,000) km/s while our 2-sigma upper limit on the expansion velocity is 6000 km/s. These limits are somewhat tighter than those derived by Bietenholz & Bartel, and confirm their conclusion that late time emission from SN 2001em, a few years after the explosion, is not driven by a relativistic jet. VLA observations of the radio flux density, at 8.46 GHz, show a decay as t^{-1.23 +/- 0.40} starting ~2.7 years after the SN. Collectively, the observations suggest interaction of the SN ejecta with a very dense circumstellar medium, though the implied opacity constraints still present a challenge.