We present the confusion-limited 1.28 GHz MeerKAT DEEP2 image covering one $\approx 68'$ FWHM primary beam area with $7.6''$ FWHM resolution and $0.55 \pm 0.01$ $\mu$Jy/beam rms noise. Its J2000 center position $\alpha=04^h 13^m 26.4^s$, $\delta=-80^\circ 00' 00''$ was selected to minimize artifacts caused by bright sources. We introduce the new 64-element MeerKAT array and describe commissioning observations to measure the primary beam attenuation pattern, estimate telescope pointing errors, and pinpoint $(u,v)$ coordinate errors caused by offsets in frequency or time. We constructed a 1.4 GHz differential source count by combining a power-law count fit to the DEEP2 confusion $P(D)$ distribution from $0.25$ to $10$ $\mu$Jy with counts of individual DEEP2 sources between $10$ $\mu$Jy and $2.5$ mJy. Most sources fainter than $S \sim 100$ $\mu$Jy are distant star-forming galaxies obeying the FIR/radio correlation, and sources stronger than $0.25$ $\mu$Jy account for $\sim93\%$ of the radio background produced by star-forming galaxies. For the first time, the DEEP2 source count has reached the depth needed to reveal the majority of the star formation history of the universe. A pure luminosity evolution of the 1.4 GHz local luminosity function consistent with the Madau & Dickinson (2014) model for the evolution of star-forming galaxies based on UV and infrared data underpredicts our 1.4 GHz source count in the range $-5 \lesssim \log[S(\mathrm{Jy})] \lesssim -4$.

Binospec is a high-throughput, 370 to 1000 nm, imaging spectrograph that addresses two adjacent 8' by 15' fields of view. Binospec was commissioned in late 2017 at the f/5 focus of the 6.5m MMT and is now available to all MMT observers. Here we describe the Binospec software used for observation planning, instrument control, and data reduction. The software and control systems incorporate a high level of automation to minimize observer workload. Instrument configuration and observation sequencing is implemented using a database-driven approach to maximize observatory efficiency. A web-based interface allows users to define observations, monitor status, and retrieve data products.

Galaxies grow their supermassive black holes in concert with their stars, although the relationship between these major galactic components is poorly understood. Observations of the cosmic growth of stars and black holes in galaxies suffer from disjoint samples and the strong effects of dust attenuation. The thermal infrared holds incredible potential for simultaneously measuring both the star formation and black hole accretion rates in large samples of galaxies covering a wide range of physical conditions. Spitzer demonstrated this potential at low redshift, and by observing some of the most luminous galaxies at z~2. JWST will apply these methods to normal galaxies at these epochs, but will not be able to generate large spectroscopic samples or access the thermal infrared at high-redshift. An order of magnitude gap in our wavelength coverage will persist between JWST and ALMA. A large, cold infrared telescope can fill this gap to determine when (in cosmic time), and where (within the cosmic web), stars and black holes co-evolve, by measuring these processes simultaneously in statistically complete and unbiased samples of galaxies to z>8. A next-generation radio interferometer will have the resolution and sensitivity to measure star-formation and nuclear accretion in even the dustiest galaxies. Together, the thermal infrared and radio can uniquely determine how stars and supermassive blackholes co-evolve in galaxies over cosmic time.

We identified 15,658 NVSS radio sources among the 55,288 2MASX galaxies brighter than $k_\mathrm{20fe} = 12.25$ at $\lambda = 2.16\,\mu\mathrm{m}$ and covering the $\Omega =7.016$ sr of sky defined by J2000 $\delta > -40^\circ$ and $\vert b \vert > 20^\circ$. The complete sample of 15,043 galaxies with 1.4 GHz flux densities $S \geq 2.45 \mathrm{~mJy}$ contains a 99.9% spectroscopically complete subsample of 9,517 galaxies with $k_\mathrm{20fe} \leq 11.75$. We used only radio and infrared data to quantitatively distinguish radio sources powered primarily by recent star formation from those powered by active galactic nuclei. The radio sources with $\log[L(\mathrm{W~Hz}^{-1})] > 19.3$ that we used to derive the local spectral luminosity and power-density functions account for $>99$% of the total 1.4~GHz spectral power densities $U_\mathrm{SF} = (1.54 \pm 0.20) \times 10^{19} \mathrm{~W~Hz}^{-1} \mathrm{~Mpc}^{-3}$ and $U_\mathrm{AGN} = (4.23 \pm 0.78) \times 10^{19} \mathrm{~W~Hz}^{-1} \mathrm{~Mpc}^{-3}$ in the universe today, and the spectroscopic subsample is large enough that the quoted errors are dominated cosmic variance. The recent comoving star-formation rate density indicated by $U_\mathrm{SF}$ is $\psi \approx 0.015~ M_\odot \mathrm{~yr}^{-1} \mathrm{~Mpc}^{-3}$.

We use Atacama Large Millimeter Array CO(3-2) observations in conjunction with optical observations from the Hubble Space Telescope to determine the ratio of stellar to gas mass for regions in the Antennae Galaxies. We adopt the term "instantaneous mass ratio" IMR(t) = M$_{stars}$/(M$_{gas}$ +M$_{stars}$), that is equivalent to the star formation efficiency for an idealized system at t = 0. We use two complementary approaches to determining the IMR(t) based on 1) the enclosed stellar and molecular mass within circular apertures centered on optically-identified clusters, and 2) a tessellation algorithm that defines regions based on CO emission. We find that only a small number of clusters appear to have IMR(0) = SFE > 0.2, which suggests that only a small fraction of these clusters will remain bound. The results suggest that by ages of $10^{6.7}$ years, some clusters will have lost all of their associated molecular gas, and by $10^{7.5}$ years this is true for the majority of clusters. There appears to be slight dependence of the IMR(t) on the CO surface brightness, which could support the idea that dense molecular environments are more likely to form bound clusters. However, the IMR(t) appears to have a strong dependence on extinction, which likely traces the evolutionary state of clusters.

The homogeneous, isotropic, and flat $\Lambda$CDM universe favored by observations of the cosmic microwave background can be described using only Euclidean geometry, locally correct Newtonian mechanics, and the basic postulates of special and general relativity. We present simple derivations of the most useful equations connecting astronomical observables (redshift, flux density, angular diameter, brightness, local space density,...) with the corresponding intrinsic properties of distant sources (lookback time, distance, spectral luminosity, linear size, specific intensity, source counts,...). We also present an analytic equation for lookback time that is accurate within 0.1% for all redshifts $z$. The exact equation for comoving distance is an elliptic integral that must be evaluated numerically, but we found a simple approximation with errors $< 0.2$% for all redshifts up to $z \approx 50$.

We present high-precision timing data over time spans of up to 11 years for 45 millisecond pulsars observed as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project, aimed at detecting and characterizing low-frequency gravitational waves. The pulsars were observed with the Arecibo Observatory and/or the Green Bank Telescope at frequencies ranging from 327 MHz to 2.3 GHz. Most pulsars were observed with approximately monthly cadence, with six high--timing-precision pulsars observed weekly, and all were observed at widely separated frequencies at each observing epoch in order to fit for time-variable dispersion delays. We describe our methods for data processing, time-of-arrival (TOA) calculation, and the implementation of a new, automated method for removing outlier TOAs. We fit a timing model for each pulsar that includes spin, astrometric, and, if necessary, binary parameters, in addition to time-variable dispersion delays and parameters that quantify pulse-profile evolution with frequency. The new timing solutions provide three new parallax measurements, two new Shapiro delay measurements, and two new measurements of large orbital-period variations. We fit models that characterize sources of noise for each pulsar. We find that 11 pulsars show significant red noise, with generally smaller spectral indices than typically measured for non-recycled pulsars, possibly suggesting a different origin. Future papers will use these data to constrain or detect the signatures of gravitational-wave signals.

We made two new sensitive (rms noise sigma_n ~ 1 microJy/beam) high resolution (theta = 3.0" and theta = 0.66" FWHM) S--band (2 < nu < 4 GHz) images covering a single JVLA primary beam (FWHM ~ 14') centered on J2000 RA = 10 46, Dec = 59 01 in the Lockman Hole. These images yielded a catalog of 792 radio sources, 97.7 +/- 0.8% of which have infrared counterparts stronger than S ~ 2 microJy at lambda = 4.5 micron. About 91% of the radio sources found in our previously published, comparably sensitive low resolution (theta = 8" FWHM) image covering the same area were also detected at 0.66" resolution, so most radio sources with S_3GHz >~ 5 microJy have angular structure phi <~ 0.66". The ratios of peak brightness in the 0.66" and 3" images have a distribution indicating that most microJy radio sources are quite compact, with a median Gaussian angular diameter <phi> = 0.3" +/- 0.1" FWHM and an rms scatter sigma_phi <~ 0.3" of individual sizes. Most of our microJy radio sources obey the tight far-infrared/radio correlation, indicating that they are powered by star formation. The median effective angular radius enclosing half the light emitted by an exponential disk is <rho_e> ~ <phi>/2.43 ~ 0.12", so the median effective radius of star-forming galaxies at redshifts z~1 is <r_e> ~ 1.0 kpc.

We find two chemically distinct populations separated relatively cleanly in the [Fe/H] - [Mg/Fe] plane, but also distinguished in other chemical planes, among metal-poor stars (primarily with metallicities [Fe/H] $< -0.9$) observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE) and analyzed for Data Release 13 (DR13) of the Sloan Digital Sky Survey. These two stellar populations show the most significant differences in their [X/Fe] ratios for the $\alpha$-elements, C+N, Al, and Ni. In addition to these populations having differing chemistry, the low metallicity high-Mg population (which we denote the HMg population) exhibits a significant net Galactic rotation, whereas the low-Mg population (or LMg population) has halo-like kinematics with little to no net rotation. Based on its properties, the origin of the LMg population is likely as an accreted population of stars. The HMg population shows chemistry (and to an extent kinematics) similar to the thick disk, and is likely associated with $\it in$ $\it situ$ formation. The distinction between the LMg and HMg populations mimics the differences between the populations of low- and high-$\alpha$ halo stars found in previous studies, suggesting that these are samples of the same two populations.

Using the nine-year radio-pulsar timing data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), collected at Arecibo Observatory and the Green Bank Telescope, we have measured the positions, proper motions, and parallaxes for 37 millisecond pulsars. We report twelve significant parallax measurements and distance measurements, and eighteen lower limits on distance. We compare these measurements to distances predicted by the NE2001 interstellar electron density model and find them to be in general agreement. We use measured orbital-decay rates and spin-down rates to confirm two of the parallax distances and to place distance upper limits on other sources; these distance limits agree with the parallax distances with one exception, PSR J1024-0719, which we discuss at length. Using the proper motions of the 37 NANOGrav pulsars in combination with other published measurements, we calculate the velocity dispersion of the millisecond pulsar population in Galactocentric coordinates. We find the radial, azimuthal, and perpendicular dispersions to be 46, 40, and 24 km s-1, respectively, in a model that allows for high-velocity outliers; or 81, 58, and 62 km s-1 for the full population. These velocity dispersions are far smaller than those of the canonical pulsar population, and are similar to older Galactic disk populations. This suggests that millisecond pulsar velocities are largely attributable to their being an old population rather than being artifacts of their birth and evolution as neutron star binary systems. The components of these velocity dispersions follow similar proportions to other Galactic populations, suggesting that our results are not biased by selection effects.

We describe the method devised to reconstruct inclined cosmic-ray air showers with zenith angles greater than $60^\circ$ detected with the surface array of the Pierre Auger Observatory. The measured signals at the ground level are fitted to muon density distributions predicted with atmospheric cascade models to obtain the relative shower size as an overall normalization parameter. The method is evaluated using simulated showers to test its performance. The energy of the cosmic rays is calibrated using a sub-sample of events reconstructed with both the fluorescence and surface array techniques. The reconstruction method described here provides the basis of complementary analyses including an independent measurement of the energy spectrum of ultra-high energy cosmic rays using very inclined events collected by the Pierre Auger Observatory.

A flux of neutrons from an astrophysical source in the Galaxy can be detected in the Pierre Auger Observatory as an excess of cosmic-ray air showers arriving from the direction of the source. To avoid the statistical penalty for making many trials, classes of objects are tested in combinations as nine "target sets", in addition to the search for a neutron flux from the Galactic Center or from the Galactic Plane. Within a target set, each candidate source is weighted in proportion to its electromagnetic flux, its exposure to the Auger Observatory, and its flux attenuation factor due to neutron decay. These searches do not find evidence for a neutron flux from any class of candidate sources. Tabulated results give the combined p-value for each class, with and without the weights, and also the flux upper limit for the most significant candidate source within each class. These limits on fluxes of neutrons significantly constrain models of EeV proton emission from non-transient discrete sources in the Galaxy.

Measurements of air showers made using the hybrid technique developed with the fluorescence and surface detectors of the Pierre Auger Observatory allow a sensitive search for point sources of EeV photons anywhere in the exposed sky. A multivariate analysis reduces the background of hadronic cosmic rays. The search is sensitive to a declination band from -85{\deg} to +20{\deg}, in an energy range from 10^17.3 eV to 10^18.5 eV. No photon point source has been detected. An upper limit on the photon flux has been derived for every direction. The mean value of the energy flux limit that results from this, assuming a photon spectral index of -2, is 0.06 eV cm^-2 s^-1, and no celestial direction exceeds 0.25 eV cm^-2 s^-1. These upper limits constrain scenarios in which EeV cosmic ray protons are emitted by non-transient sources in the Galaxy.

The MMT and Magellan infrared spectrograph (MMIRS) is a cryogenic multiple slit spectrograph operating in the wavelength range 0.9-2.4 micron. MMIRS' refractive optics offer a 6.9 by 6.9 arcmin field of view for imaging with a spatial resolution of 0.2 arcsec per pixel on a HAWAII-2 array. For spectroscopy, MMIRS can be used with long slits up to 6.9 arcmin long, or with custom slit masks having slitlets distributed over a 4 by 6.9 arcmin area. A range of dispersers offer spectral resolutions of 800 to 3000. MMIRS is designed to be used at the f/5 foci of the MMT or Magellan Clay 6.5m telescopes. MMIRS was commissioned in 2009 at the MMT and has been in routine operation at the Magellan Clay Telescope since 2010. MMIRS is being used for a wide range of scientific investigations from exoplanet atmospheres to Ly-alpha emitters.