On 2012 August 2, two CMEs (CME-1 and CME-2) erupted from the west limb of the Sun as viewed from Earth, and were observed in images from the white light coronagraphs on the SOHO and STEREO spacecraft. These events were also observed by the Very Large Array (VLA), which was monitoring the Sun at radio wavelengths, allowing time-dependent Faraday rotation observations to be made of both events. We use the white-light imaging and radio data to model the 3-D field geometry of both CMEs, assuming a magnetic flux rope geometry. For CME-2, we also consider 1 au in situ field measurements in the analysis, as this CME hits STEREO-A on August~6, making this the first CME with observational constraints from stereoscopic coronal imaging, radio Faraday rotation, and in situ plasma measurements combined. The imaging and in situ observations of CME-2 provide two clear predictions for the radio data; namely that VLA should observe positive rotation measures (RMs) when the radio line of sight first encounters the CME, and that the sign should reverse to negative within a couple hours. The initial positive RMs are in fact observed. The expected sign reversal is not, but the VLA data unfortunately end too soon to be sure of the significance of this discrepancy. We interpret an RM increase prior to the expected occultation time of the CME as a signature of a sheath region of deflected field ahead of the CME itself.

As part of an on-going study of radio transients in Epoch 1 (2017-2019) of the Very Large Array Sky Survey (VLASS), we have discovered a sample of 0.2<z<3.2 active galactic nuclei (AGN) selected in the optical/infrared that have recently brightened dramatically in the radio. These sources would have previously been classified as radio-quiet based on upper limits from the Faint Images of the Radio Sky at Twenty-centimeters (FIRST; 1993-2011) survey; however, they are now consistent with radio-loud quasars. We present a quasi-simultaneous, multi-band (1-18 GHz) VLA follow-up campaign of our sample of AGN with extreme radio variability. We conclude that the radio properties are most consistent with AGN that have recently launched jets within the past few decades, potentially making them among the youngest radio AGN known.

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.

We report on a search for radio transients at 340 MHz with the Jansky Very Large Array (VLA) Low band Ionospheric and Transient Experiment (VLITE). Between 2015 July 29 and 2015 September 27, operating in commensal mode, VLITE imaged approximately 2800 pointings covering 12,000 deg$^2$ on the sky, sampling timescales ranging from tens of seconds to several hours on a daily basis. In addition, between 2015 February 25 and 2015 May 9, VLITE observed 55 epochs of roughly 2-4 hours each toward the COSMOS field. Using existing radio source catalogs we have searched all the daily VLITE images for transients, while for the COSMOS field we compared individual images and the summed image to search for new sources in repeated observations of the same field. The wide range of timescales makes VLITE sensitive to both coherent and incoherent transient source classes. No new transients are found, allowing us to set stringent upper limits on transients at milliJansky levels and at low frequencies where comparatively few such surveys have been carried out to date. An all-sky isotropic rate of bursting radio transients with similar rates, duration, and intensity as the unusual transient GCRT J1745$-$3009, discovered in wide-field monitoring toward the Galactic center, is ruled out with high confidence. The resulting non-detections allows us to argue that this is a coherent source, whose properties most resemble the growing class of nulling pulsars. We end with a discussion of the future prospects for the detection of transients by VLITE and other experiments.

The US Naval Research Laboratory (NRL) and the National Radio Astronomy Observatory (NRAO) have collaborated to develop, install, and commission a new commensal system on the Karl G. Jansky Very Large Array (VLA). The VLA Low Band Ionospheric and Transient Experiment (VLITE) makes use of dedicated samplers and fibers to tap the signal from 10 VLA low band receivers and correlate those through a real-time DiFX correlator. VLITE allows for the simultaneous use of the VLA to observe primary science using the higher frequencies receivers (1-50 GHz) through the NRAO WIDAR correlator and lower frequencies through the DiFX correlator. VLITE operates during nearly all observing programs and provides 64 MHz of bandwidth centered at 352 MHz. The operation of VLITE requires no additional resources from the VLA system running the primary science and produces an ad-hoc sky survey. The commensal system greatly expands the capabilities of the VLA through value-added PI science, stand-alone astrophysics, the opening of a new window on transient searches, and serendipity. In the first year of operation we have recorded more than 6300 hours spread across the sky. We present an overview of the VLITE program, discuss the sky coverage and depth obtained during the first year of operation, and briefly outline the astrophysics and transients programs.

We use N-body simulations of dark matter haloes in cold dark matter (CDM) and a large set of different warm dark matter (WDM) cosmologies to demonstrate that the spherically averaged density profile of dark matter haloes has a shape that depends on the power spectrum of matter perturbations. Density profiles are steeper in WDM but become shallower at scales less than one percent of the virial radius. Virialization isotropizes the velocity dispersion in the inner regions of the halo but does not erase the memory of the initial conditions in phase space. The location of the observed deviations from CDM in the density profile and in phase space can be directly related to the ratio between the halo mass and the filtering mass and are most evident in small mass haloes, even for a 34 keV thermal relic WDM. The rearrangement of mass within the haloes supports analytic models of halo structure that include angular momentum. We also find evidence of a dependence of the slope of the inner density profile in CDM cosmologies on the halo mass with more massive haloes exhibiting steeper profiles, in agreement with the model predictions and with previous simulation results. Our work complements recent studies of microhaloes near the filtering scale in CDM and strongly argue against a universal shape for the density profile.

We report the detection of far-infrared (FIR) CO rotational emission from nearby active galactic nuclei (AGN) and starburst galaxies, as well as several merging systems and Ultra-Luminous Infrared Galaxies (ULIRGs). Using Herschel-PACS, we have detected transitions in the J$_{upp}$ = 14 - 20 range ($\lambda \sim$ 130 - 185 $\mu$m, $\nu \sim$ 1612 - 2300 GHz) with upper limits on (and in two cases, detections of) CO line fluxes up to J$_{upp}$ = 30. The PACS CO data obtained here provide the first well-sampled FIR extragalactic CO SLEDs for this range, and will be an essential reference for future high redshift studies. We find a large range in the overall SLED shape, even amongst galaxies of similar type, demonstrating the uncertainties in relying solely on high-J CO diagnostics to characterize the excitation source of a galaxy. Combining our data with low-J line intensities taken from the literature, we present a CO ratio-ratio diagram and discuss its potential diagnostic value in distinguishing excitation sources and physical properties of the molecular gas. The position of a galaxy on such a diagram is less a signature of its excitation mechanism, than an indicator of the presence (or absence) of warm, dense molecular gas. We then quantitatively analyze the CO emission from a subset of the detected sources with Large Velocity Gradient (LVG) radiative transfer models to fit the CO SLEDs. Using both single-component and two-component LVG models to fit the kinetic temperature, velocity gradient, number density and column density of the gas, we derive the molecular gas mass and the corresponding CO-to-H$_2$ conversion factor, $\alpha_{CO}$, for each respective source. For the ULIRGs we find $\alpha$ values in the canonical range 0.4 - 5 M$_\odot$/(K kms$^{-1}$pc$^2$), while for the other objects, $\alpha$ varies between 0.2 and 14.} Finally, we compare our best-fit LVG model ..

We report on the Herschel/PACS observations of OH in Mrk 231, with detections in 9 doublets observed within the PACS range, and present radiative transfer models for the outflowing OH. Signatures of outflowing gas are found in up to 6 OH doublets with different excitation requirements. At least two outflowing components are identified, one with OH radiatively excited, and the other with low excitation, presumably spatially extended. Particularly prominent, the blue wing of the absorption detected in the in-ladder 2Pi_{3/2} J=9/2-7/2 OH doublet at 65 um, with E_lower=290 K, indicates that the excited outflowing gas is generated in a compact and warm (circum)nuclear region. Because the excited, outflowing OH gas in Mrk 231 is associated with the warm, far-IR continuum source, it is likely more compact (diameter of 200-300 pc) than that probed by CO and HCN. Nevertheless, its mass-outflow rate per unit of solid angle as inferred from OH is similar to that previously derived from CO, >~70x(2.5x10^{-6}/X_{OH}) Msun yr^{-1} sr^{-1}, where X_{OH} is the OH abundance relative to H nuclei. In spherical symmetry, this would correspond to >~850x(2.5x10^{-6}/X_{OH}) Msun yr^{-1}, though significant collimation is inferred from the line profiles. The momentum flux of the excited component attains ~15 L_{AGN}/c, with an OH column density of (1.5-3)x10^{17} cm^-2 and a mechanical luminosity of ~10^{11} Lsun. The detection of very excited OH peaking at central velocities indicates the presence of a nuclear reservoir of gas rich in OH, plausibly the 130-pc scale circumnuclear torus previously detected in OH megamaser emission, that may be feeding the outflow. An exceptional ^{18}OH enhancement, with OH/^{18}OH<~30 at both central and blueshifted velocities, is likely the result of interstellar-medium processing by recent starburst/SNe activity.

We investigate the claim that the largest subhaloes in high resolution dissipationless cold dark matter (CDM) simulations of the Milky Way are dynamically inconsistent with observations of its most luminous satellites. We find that the inconsistency is largely attributable to the large values of \sigma_8 and n_s adopted in the discrepant simulations producing satellites that form too early and therefore are too dense. We find the tension between observations and simulations adopting parameters consistent with WMAP9 is greatly diminished making the satellites a sensitive test of CDM. We find the Via Lactea II halo to be atypical for haloes in a WMAP3 cosmology, a discrepancy that we attribute to its earlier formation epoch than the mean for its mass. We also explore warm dark matter (WDM) cosmologies for 1--4 keV thermal relics. In 1 keV cosmologies subhaloes have circular velocities at kpc scales ~ 60% lower than their CDM counterparts, but are reduced by only 10% in 4 keV cosmologies. Since relic masses < 2-3 keV are ruled out by constraints from the number of Milky Way satellites and Lyman-\alpha\ forest, WDM has a minor effect in reducing the densities of massive satellites. Given the uncertainties on the mass and formation epoch of the Milky Way, the need for reducing the satellite densities with baryonic effects or WDM is alleviated.

We report the results from a systematic search for molecular (OH-119 um) outflows with Herschel-PACS in a sample of 43 nearby (z < 0.3) galaxy mergers, mostly ultraluminous infrared galaxies (ULIRGs) and QSOs. We find that the character of the OH feature (strength of the absorption relative to the emission) correlates with that of the 9.7-um silicate feature, a measure of obscuration in ULIRGs. Unambiguous evidence for molecular outflows, based on the detection of OH absorption profiles with median velocities more blueshifted than -50 km/sec, is seen in 26 (70%) of the 37 OH-detected targets, suggesting a wide-angle (~145 degrees) outflow geometry. Conversely, unambiguous evidence for molecular inflows, based on the detection of OH absorption profiles with median velocities more redshifted than +50 km/sec, is seen in only 4 objects, suggesting a planar or filamentary geometry for the inflowing gas. Terminal outflow velocities of ~-1000 km/sec are measured in several objects, but median outflow velocities are typically ~-200 km s^{-1}. While the outflow velocities show no statistically significant dependence on the star formation rate, they are distinctly more blueshifted among systems with large AGN fractions and luminosities [log (L_AGN / L_sun) > 11.8 +/- 0.3]. The quasars in these systems play a dominant role in driving the molecular outflows. In contrast, the most AGN dominated systems, where OH is seen purely in emission, show relatively modest OH line widths, despite their large AGN luminosities, perhaps indicating that molecular outflows subside once the quasar has cleared a path through the obscuring material.

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 present a comprehensive analysis of the halo model of cosmological large to small-scale structure statistics in the case of warm dark matter (WDM) structure formation scenarios. We include the effects of WDM on the linear matter power spectrum, halo density profile, halo concentration relation, halo mass function, subhalo density profile, subhalo mass function and biasing of the smooth dark matter component. As expected, we find large differences at the smallest physical scales in the nonlinear matter power spectrum predicted in the halo model between WDM and cold dark matter even for reasonably high-scale WDM particle masses. We find that significant effects are contributed from the alteration of the halo density profile and concentration, as well as the halo mass function. We further find that the effects of WDM on the subhalo population are important but sub-dominant. Clustering effects of the biasing of the smooth component in WDM is not largely significant.

We combine high resolution N-body simulations with deep observations of neutral hydrogen (HI) in nearby galaxy groups in order to explore two well-known theories of HI cloud formation: HI stripping by galaxy interactions and dark matter minihalos with embedded HI gas. This paper presents new data from three galaxy groups, Canes Venatici I, NGC 672, and NGC 45, and assembles data from our previous galaxy group campaign to generate a rich HI cloud archive to compare to our simulated data. We find no HI clouds in the Canes Venatici I, NGC 672, or NGC 45 galaxy groups. We conclude that HI clouds in our detection space are most likely to be generated through recent, strong galaxy interactions. We find no evidence of HI clouds associated with dark matter halos above M_HI = 10^6 M_Sun, within +/- 700 km/s of galaxies, and within 50 kpc projected distance of galaxies.

We have conducted N-body simulations of the growth of Milky Way-sized halos in cold and warm dark matter cosmologies. The number of dark matter satellites in our simulated Milky Ways decreases with decreasing mass of the dark matter particle. Assuming that the number of dark matter satellites exceeds or equals the number of observed satellites of the Milky Way we derive lower limits on the dark matter particle mass. We find with 95% confidence m_s > 13.3 keV for a sterile neutrino produced by the Dodelson and Widrow mechanism, m_s > 8.9 keV for the Shi and Fuller mechanism, m_s > 3.0 keV for the Higgs decay mechanism, and m_{WDM} > 2.3 keV for a thermal dark matter particle. The recent discovery of many new dark matter dominated satellites of the Milky Way in the Sloan Digital Sky Survey allows us to set lower limits comparable to constraints from the complementary methods of Lyman-alpha forest modeling and X-ray observations of the unresolved cosmic X-ray background and of dark matter halos from dwarf galaxy to cluster scales. Future surveys like LSST, DES, PanSTARRS, and SkyMapper have the potential to discover many more satellites and further improve constraints on the dark matter particle mass.

This paper presents a search for radio transients at a frequency of 73.8 MHz (4 m wavelength) using the all-sky imaging capabilities of the Long Wavelength Demonstrator Array (LWDA). The LWDA was a 16-dipole phased array telescope, located on the site of the Very Large Array in New Mexico. The field of view of the individual dipoles was essentially the entire sky, and the number of dipoles was sufficiently small that a simple software correlator could be used to make all-sky images. From 2006 October to 2007 February, we conducted an all-sky transient search program, acquiring a total of 106 hr of data; the time sampling varied, being 5 minutes at the start of the program and improving to 2 minutes by the end of the program. We were able to detect solar flares, and in a special-purpose mode, radio reflections from ionized meteor trails during the 2006 Leonid meteor shower. We detected no transients originating outside of the solar system above a flux density limit of 500 Jy, equivalent to a limit of no more than about 10^{-2} events/yr/deg^2, having a pulse energy density >~ 1.5 x 10^{-20} J/m^2/Hz at 73.8 MHz for pulse widths of about 300 s. This event rate is comparable to that determined from previous all-sky transient searches, but at a lower frequency than most previous all-sky searches. We believe that the LWDA illustrates how an all-sky imaging mode could be a useful operational model for low-frequency instruments such as the Low Frequency Array, the Long Wavelength Array station, the low-frequency component of the Square Kilometre Array, and potentially the Lunar Radio Array.