We report the detection of a potential cosmic radio transient source using the two stations of the Long Wavelength Array. The transient was detected on 18 October 2017 08:47 UTC near the celestial equator while reducing 10,240 hours of archival all-sky images from the LWA1 and LWA-SV stations. The detected transient at 34 MHz has a duration of 15 - 20 seconds and a flux density of 842 +/- 116 Jy at LWA1 and 830 +/- 92 Jy at LWA-SV. The transient source has not repeated, and its nature is not well understood. The Pan-STARRS optical telescope has detected a supernova that occurred on the edge of the position error circle of the transient on the same day.

Radio emission from meteors or meteor radio afterglows (MRAs) were first detected using the all-sky imaging capabilities of the first station of the Long Wavelength Array (LWA1). In this work, we use the recently commissioned LWA Sevilleta (LWA-SV) station along with the LWA1 to carry out co-ordinated observations. The combined all-sky observations with LWA1 and LWA-SV have co-observed 32 MRAs and 21 transmitter reflections from meteors (meteor scatter events) which are believed to be specular reflections from overdense trails. The flux density of the events observed by each station were measured from the all-sky images. Triangulating the angular direction of events from each station gave the physical location and the distance of the event to each station. The luminosity of the events in each station were calculated using the flux distance relation for an isotropic source. The luminosity distribution for MRAs and meteor scatter events observed by each station shows a clear distinction between these two types of events as the ratio of luminosities are closer to unity for MRAs than the meteor scatter events. Furthermore, we find that MRAs follow an isotropic radiation pattern. This suggests, either a complete incoherent emission mechanism or an incoherent addition of coherently emitting small regions within the meteor trail.

Using the first station of the Long Wavelength Array (LWA1), we examine polarized pulsar emission between 25 and 88 MHz. Polarized light from pulsars undergoes Faraday rotation as it passes through the magnetized interstellar medium. Observations from low frequency telescopes are ideal for obtaining precise rotation measures (RMs) because the effect of Faraday rotation is proportional to the square of the observing wavelength. With these RMs, we obtained polarized pulse profiles to see how polarization changes in the 25-88 MHz range. The RMs were also used to derive values for the electron density weighted average Galactic magnetic field along the line of sight. We present rotation measures and polarization profiles of 15 pulsars acquired using data from LWA1. These results provide new insight into low-frequency polarization characteristics and pulsar emission heights, and complement measurements at higher frequencies.

We observed the flare stars AD Leonis, Wolf 424, EQ Pegasi, EV Lacertae, and UV Ceti for nearly 135 hours. These stars were observed between 63 and 83 MHz using the interferometry mode of the Long Wavelength Array. Given that emission from flare stars is typically circularly polarized, we used the condition that any significant detection present in Stokes I must also be present in Stokes V at the same time in order for us to consider it a possible flare. Following this, we made one marginal flare detection for the star EQ Pegasi. This flare had a flux density of 5.91 Jy in Stokes I and 5.13 Jy in Stokes V, corresponding to a brightness temperature $1.75 \times 10^{16}(r/r_*)^{-2}$ K.

The search for the spectral signature of hydrogen from the formation of the first stars, known as Cosmic Dawn or First Light, is an ongoing effort around the world. The signature should present itself as a decrease in the temperature of the 21-cm transition relative to that of the Cosmic Microwave Background and is believed to reside somewhere below 100 MHz. A potential detection was published by the Experiment to Detect the Global EoR Signal (EDGES) collaboration with a profile centered around 78 MHz of both unexpected depth and width (Bowman et. al 2018; arXiv:1810.05912). If validated, this detection will have profound impacts on the current paradigm of structure formation within $\Lambda$CDM cosmology. We present an attempt to detect the spectral signature reported by the EDGES collaboration with the Long Wavelength Array station located on the Sevilleta National Wildlife Refuge in New Mexico, USA (LWA-SV). LWA-SV differs from other instruments in that it is a 256 element antenna array and offers beamforming capabilities that should help with calibration and detection. We report first limits from LWA-SV and look toward future plans to improve these limits.

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.

PSR B1508+55 is known to have a single component profile above 300 MHz. However, when we study it at frequencies below 100 MHz using the first station of the Long Wavelength Array, it shows multiple components. These include the main pulse, a precursor, a postcursor, and a trailing component. The separation of the trailing component from the main peak evolves over the course of a three year study. This evolution is likely an effect of the pulse signal getting refracted off an ionized gas cloud (acting as a lens) leading to what appears to be a trailing component in the profile as the pulsar signal traverses the interstellar medium. Using this interpretation, we identify the location and electron density of the lens affecting the pulse profile.

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.

The recent demonstration of a real-time direct imaging radio interferometry correlator represents a new capability in radio astronomy. However wide field imaging with this method is challenging since wide-field effects and array non-coplanarity degrade image quality if not compensated for. Here we present an alternative direct imaging correlation strategy using a Direct Fourier Transform (DFT), modelled as a linear operator facilitating a matrix multiplication between the DFT matrix and a vector of the electric fields from each antenna. This offers perfect correction for wide field and non-coplanarity effects. When implemented with data from the Long Wavelength Array (LWA), it offers comparable computational performance to previously demonstrated direct imaging techniques, despite having a theoretically higher floating point cost. It also has additional benefits, such as imaging sparse arrays and control over which sky co-ordinates are imaged, allowing variable pixel placement across an image. It is in practice a highly flexible and efficient method of direct radio imaging when implemented on suitable arrays. A functioning Electric Field Direct imaging architecture using the DFT is presented, alongside an exploration of techniques for wide-field imaging similar to those in visibility based imaging, and an explanation of why they do not fit well to imaging directly with the digitized electric field data. The DFT imaging method is demonstrated on real data from the LWA telescope, alongside a detailed performance analysis, as well as an exploration of its applicability to other arrays.

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.

The future of radio astronomy will require instruments with large collecting areas for higher sensitivity, wide fields of view for faster survey speeds, and efficient computing and data rates relative to current capabilities. We describe the first successful deployment of the E-field Parallel Imaging Correlator (EPIC) on the LWA station in Sevilleta, New Mexico, USA (LWA-SV). EPIC is a solution to the computational problem of large interferometers. By gridding and spatially Fourier transforming channelised electric fields from the antennas in real-time, EPIC removes the explicit cross multiplication of all pairs of antenna voltages to synthesize an aperture, reducing the computational scaling from $\mathcal{O}(n_a^2)$ to $\mathcal{O}(n_g \log_2 n_g)$, where $n_a$ is the number of antennas and $n_g$ is the number of grid points. Not only does this save computational costs for dense arrays but it produces very high time resolution images in real time. The GPU-based implementation uses existing LWA-SV hardware and the high performance streaming framework, Bifrost. We examine the practical details of the EPIC deployment and verify the imaging performance by detecting a meteor impact on the atmosphere using continuous all-sky imaging at 50 ms time resolution.

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.

Interstellar scattering causes pulsar profiles to grow asymmetrically, thus affecting the pulsar timing residuals, and is strongest at lower frequencies. Different Interstellar medium models predict different frequency ($\nu$) and dispersion measure (DM) dependencies for the scattering time-scale $\tau_{sc}$. For Gaussian inhomogeneity the expected scaling relation is $\tau_{sc} \propto \nu^{-4}\ DM^{2}$, while for a Kolmogorov distribution of irregularities, the expected relation is $\tau_{sc} \propto \nu^{-4.4}\ DM^{2.2}$. Previous scattering studies show a wide range of scattering index across all ranges of DM. A scattering index below 4 is believed to be either due to limitations of the underlying assumptions of the thin screen model or an anisotropic scattering mechanism. We present a study of scattering for seven nearby pulsars (DM $< 50$ pc cm$^{-3}$) observed at low frequencies ($10-88$ MHz), using the first station of the Long Wavelength Array (LWA1). We examine the scattering spectral index and DM variation over a period of about three years. The results yield insights into the small-scale structure of ISM as well as the applicability of the thin screen model for low DM pulsars.

As telescope facilities become increasingly more capable they also become increasingly complex and require additional resources to operate. This is particularly true for the current and future generations of "software defined telescopes" that can support a variety of observing programs simultaneously, either through commensal observations or through support for multiple pointing centers as in the case of dipole arrays or dishes equipped with phased array feeds. At the same time, many current and future facilities are also distributed over large geographic areas, making monitoring and maintenance more difficult and costly. For these reasons we have developed a new paradigm for telescope operations called the "swarm telescope" that breaks large, single facilities into smaller groups of independent systems that can collaboratively work together to function as a single facility but with much less operational overhead. In this paper we outline the swarm telescope concept and an example of its implementation at the Long Wavelength Array. We also discuss potential advantages of using this approach for other facilities, in particular the Next Generation Very Large Array.

The ionosphere contributes time-varying Faraday Rotation (FR) to radio signals passing through it. Correction for the effect of the ionosphere is important for deriving magnetic field information from FR observations of polarized cosmic radio sources, as well as providing valuable diagnostics of the structure of the ionosphere. In this paper, we evaluate the accuracy of models commonly used to correct for its effects using new observations of pulsars at low frequencies, which provide total rotation measures (RM) at better precision than previously available. We evaluate models of the ionosphere derived from modern digital ionosondes that provide electron density information as a function of height, as well as GPS-derived Total Electron Content (TEC) measurements. We combine these density models with reference global magnetic field models to derive ionospheric RM contributions. We find that the models disagree substantially with each other and seek corrections that may explain the differences in RM prediction. Additionally we compare these models to global TEC models and find that local high-cadence TEC measurements are superior to global models for ionospheric RM correction.

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.

The recent detection of the "cosmic dawn" redshifted 21 cm signal at 78 MHz by the EDGES experiment differs significantly from theoretical predictions. In particular, the absorption trough is roughly a factor of two stronger than the most optimistic theoretical models. The early interpretations of the origin of this discrepancy fall into two categories. The first is that there is increased cooling of the gas due to interactions with dark matter, while the second is that the background radiation field includes a contribution from a component in addition to the cosmic microwave background. In this paper we examine the feasibility of the second idea using new data from the first station of the Long Wavelength Array. The data span 40 to 80 MHz and provide important constraints on the present-day background in a frequency range where there are few surveys with absolute temperature calibration suitable for measuring the strength of the radio monopole. We find support for a strong, diffuse radio background that was suggested by the ARCARDE 2 results in the 3 to 10 GHz range. We find that this background is well modeled by a power law with a spectral index of $-$2.58$\pm$0.05 and a temperature at the rest frame 21 cm frequency of 603$^{+102}_{-92}$ mK.

We summarize the radio synchrotron background workshop that took place July 19-21, 2017 at the University of Richmond. This first scientific meeting dedicated to the topic was convened because current measurements of the diffuse radio monopole reveal a surface brightness that is several times higher than can be straightforwardly explained by known Galactic and extragalactic sources and processes, rendering it by far the least well understood photon background at present. It was the conclusion of a majority of the participants that the radio monopole level is at or near that reported by the ARCADE 2 experiment and inferred from several absolutely calibrated zero level lower frequency radio measurements, and unanimously agreed that the production of this level of surface brightness, if confirmed, represents a major outstanding question in astrophysics. The workshop reached a consensus on the next priorities for investigations of the radio synchrotron background.

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

We present results of the coordinated observing campaign that made the first subarcsecond localization of a Fast Radio Burst, FRB 121102. During this campaign, we made the first simultaneous detection of an FRB burst by multiple telescopes: the VLA at 3 GHz and the Arecibo Observatory at 1.4 GHz. Of the nine bursts detected by the Very Large Array at 3 GHz, four had simultaneous observing coverage at other observatories. We use multi-observatory constraints and modeling of bursts seen only at 3 GHz to confirm earlier results showing that burst spectra are not well modeled by a power law. We find that burst spectra are characterized by a ~500 MHz envelope and apparent radio energy as high as $10^{40}$ erg. We measure significant changes in the apparent dispersion between bursts that can be attributed to frequency-dependent profiles or some other intrinsic burst structure that adds a systematic error to the estimate of DM by up to 1%. We use FRB 121102 as a prototype of the FRB class to estimate a volumetric birth rate of FRB sources $R_{FRB} \approx 5x10^{-5}/N_r$ Mpc$^{-3}$ yr$^{-1}$, where $N_r$ is the number of bursts per source over its lifetime. This rate is broadly consistent with models of FRBs from young pulsars or magnetars born in superluminous supernovae or long gamma-ray bursts, if the typical FRB repeats on the order of thousands of times during its lifetime.

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.

On 2015 February 21, simultaneous observations of Jupiter's decametric radio emission between 10 and 33 MHz were carried out using three powerful low-frequency radio telescopes: Long Wavelength Array Station One (LWA1) in USA; Nan\c{c}ay Decameter Array (NDA) in France; and URAN2 telescope in Ukraine. We measure lag times of short-bursts (S-bursts) for 105-minutes of data over effective baselines up to 8460 km by using cross-correlation analysis of the spectrograms from each instrument. Of particular interest is the measurement of the beaming thickness of S-bursts, testing if either flashlight- or beacon-like beaming is emanating from Jupiter. We find that the lag times for all pairs drift slightly as time elapses, in agreement with expectations from the flashlight-like beaming model. This leads to a new constraint of the minimum beaming thickness of 2.66". Also, we find that most of the analyzed data abound with S-bursts, whose occurrence probability peaks at 17-18 MHz.

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.

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 explore opportunities for multi-messenger astronomy using gravitational waves (GWs) and prompt, transient low-frequency radio emission to study highly energetic astrophysical events. We review the literature on possible sources of correlated emission of gravitational waves and radio transients, highlighting proposed mechanisms that lead to a short-duration, high-flux radio pulse originating from the merger of two neutron stars or from a superconducting cosmic string cusp. We discuss the detection prospects for each of these mechanisms by low-frequency dipole array instruments such as LWA1, LOFAR and MWA. We find that a broad range of models may be tested by searching for radio pulses that, when de-dispersed, are temporally and spatially coincident with a LIGO/Virgo GW trigger within a $\usim 30$ second time window and $\usim 200 \mendash 500 \punits{deg}^{2}$ sky region. We consider various possible observing strategies and discuss their advantages and disadvantages. Uniquely, for low-frequency radio arrays, dispersion can delay the radio pulse until after low-latency GW data analysis has identified and reported an event candidate, enabling a \emph{prompt} radio signal to be captured by a deliberately targeted beam. If neutron star mergers do have detectable prompt radio emissions, a coincident search with the GW detector network and low-frequency radio arrays could increase the LIGO/Virgo effective search volume by up to a factor of $\usim 2$. For some models, we also map the parameter space that may be constrained by non-detections.

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.

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.

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.

We have obtained radial velocity measurements for 51 new globular clusters around the Sombrero galaxy. These measurements were obtained using spectroscopic observations from the AAOmega spectrograph on the Anglo-Australian Telescope and the Hydra spectrograph at WIYN. Combined with our own past measurements and velocity measurements obtained from the literature we have constructed a large database of radial velocities that contains a total of 360 confirmed globular clusters. Previous studies' analyses of the kinematics and mass profile of the Sombrero globular cluster system have been constrained to the inner ~9' (~24 kpc or ~5 effective radii), but our new measurements have increased the radial coverage of the data, allowing us to determine the kinematic properties of M104 out to ~15' (~41 kpc or ~9 effective radii). We use our set of radial velocities to study the GC system kinematics and to determine the mass profile and V-band mass-to-light profile of the galaxy. We find that the V-band mass-to-light ratio increases from 4.5 at the center to a value of 20.9 at 41 kpc (~9 effective radii or 15'), which implies that the dark matter halo extends to the edge of our available data set. We compare our mass profile at 20 kpc (~4 effective radii or ~7.4') to the mass computed from x-ray data and find good agreement. We also use our data to look for rotation in the globular cluster system as a whole, as well as in the red and blue subpopulations. We find no evidence for significant rotation in any of these samples.

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.

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.

LWA1 is a new large radio telescope array operating in the frequency range 10-88 MHz, located in central New Mexico. The telescope consists of 260 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 telescope's "all dipoles" modes, facilitating all-sky imaging. Notable features of the instrument include four independently-steerable beams utilizing digital "true time delay" beamforming, high intrinsic sensitivity (about 6 kJy zenith system equivalent flux density), large instantaneous bandwidth (up to 78 MHz), and large field of view (about 3-10 degrees, depending on frequency and zenith angle of pointing). This paper summarizes the design of LWA1, its performance as determined in commissioning experiments, and results from early science observations demonstrating the capabilities of the instrument.

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 present results from a study of the globular cluster (GC) systems of four spiral and S0 galaxies imaged as part of an ongoing wide-field survey of the GC systems of giant galaxies. The target galaxies -- the SB0 galaxy NGC1023, the SBb galaxy NGC1055, and an isolated pair comprised of the Sbc galaxy NGC7339 and the S0 galaxy NGC7332 -- were observed in BVR filters with the WIYN 3.5-m telescope and Minimosaic camera. For two of the galaxies, we combined the WIYN imaging with previously-published data from the Hubble Space Telescope and the Keck Observatory to help characterize the GC distribution in the central few kiloparsecs. We determine the radial distribution (surface density of GCs versus projected radius) of each galaxy's GC system and use it to calculate the total number of GCs (N_GC). We find N_GC = 490+/-30, 210+/-40, 175+/-15, and 75+/-10 for NGC1023, NGC1055, NGC7332, and NGC7339, respectively. We also calculate the GC specific frequency (N_GC normalized by host galaxy luminosity or mass) and find values typical of those of the other spiral and E/S0 galaxies in the survey. The two lenticular galaxies have sufficient numbers of GC candidates for us to perform statistical tests for bimodality in the GC color distributions. We find evidence at a high confidence level (>95%) for two populations in the B-R distribution of the GC system of NGC1023. We find weaker evidence for bimodality (>81% confidence) in the GC color distribution of NGC7332. Finally, we identify eight GC candidates that may be associated with the Magellanic dwarf galaxy NGC1023A, a satellite of NGC1023.

The Long Wavelength Array Software Library (LSL) is a Python module that provides a collection of utilities to analyze and export data collected at the first station of the Long Wavelength Array, LWA1. Due to the nature of the data format and large-N ($\gtrsim$100 inputs) challenges faced by the LWA, currently available software packages are not suited to process the data. Using tools provided by LSL, observers can read in the raw LWA1 data, synthesize a filter bank, and apply incoherent de-dispersion to the data. The extensible nature of LSL also makes it an ideal tool for building data analysis pipelines and applying the methods to other low frequency arrays.

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.

Neutral hydrogen clouds are found in the Milky Way and Andromeda halo both as large complexes and smaller isolated clouds. Here we present a search for Hi clouds in the halo of M33, the third spiral galaxy of the Local Group. We have used two complementary data sets: a 3^o x 3^o map of the area provided by the Arecibo Legacy Fast ALFA (ALFALFA) survey and deeper pointed observations carried out with the Arecibo telescope in two fields that permit sampling of the north eastern and south-western edges of the HI disc. The total amount of Hi around M33 detected by our survey is $\sim 10^7$ M$_{\odot}$. At least 50% of this mass is made of HI clouds that are related both in space and velocity to the galaxy. We discuss several scenarios for the origin of these clouds focusing on the two most interesting ones: $(a)$ dark-matter dominated gaseous satellites, $(b)$ debris from filaments flowing into M33 from the intergalactic medium or generated by a previous interaction with M31. Both scenarios seem to fit with the observed cloud properties. Some structures are found at anomalous velocities, particularly an extended HI complex previously detected by Thilker et al. (2002). Even though the ALFALFA observations seem to indicate that this cloud is possibly connected to M33 by a faint gas bridge, we cannot firmly establish its extragalactic nature or its relation to M33. Taking into account that the clouds associated with M33 are likely to be highly ionised by the extragalactic UV radiation, we predict that the total gas mass associated with them is > 5 x 10^7 M$_{\odot}$. If the gas is steadily falling towards the M33 disc it can provide the fuel needed to sustain a current star formation rate of 0.5 M$_{\odot}$ yr$^{-1}$.

The initial cluster mass function (ICMF) is a fundamental property of star formation in galaxies. To gauge its universality, we measure and compare the ICMFs in irregular and spiral galaxies. Our sample of irregular galaxies is based on 13 nearby galaxies selected from a volume-limited sample from the fifth data release of the Sloan Digital Sky Survey (SDSS). The extinctions, ages, and masses were determined by comparing their u'g'i'z' magnitudes to those generated from starburst models. Completeness corrections were performed using Monte Carlo simulations in which artificial clusters were inserted into each galaxy. We analyzed three nearby spiral galaxies with SDSS data in exactly the same way to derive their ICMF based on a similar number of young, massive clusters as the irregular galaxy ICMF. We find that the ICMFs of irregular and spiral galaxies for masses >3x10^4 M_sun are statistically indistinguishable. For clusters more massive than 3x10^4 M_sun, the ICMF of the irregular galaxies is reasonably well fit by a power law dN(M)/dM ~ M^-a_M with a_M = 1.88 +/- 0.09. Similar results were obtained for the ICMF of the spiral galaxy sample but with a_M = 1.75 +/- 0.06. We discuss the implications of our result for theories of star cluster formation: the shape of the ICMF appears to be independent of metallicity and galactic shear rate.