Phosphorus is a key ingredient in terrestrial biochemistry, but is rarely observed in the molecular ISM and therefore little is known about how it is inherited during the star and planet formation sequence. We present observations of the phosphorus-bearing molecules PO and PN towards the Class I low-mass protostar B1-a using the IRAM 30m telescope, representing the second detection of phosphorus carriers in a Solar-type star forming region. The P/H abundance contained in PO and PN is ~10$^{-10}$-10$^{-9}$ depending on the assumed source size, accounting for just 0.05-0.5% of the solar phosphorus abundance and implying significant sequestration of phosphorus in refractory material. Based on a comparison of the PO and PN line profiles with the shock tracers SiO, SO$_2$, and CH$_3$OH, the phosphorus molecule emission seems to originate from shocked gas and is likely associated with a protostellar outflow. We find a PO/PN column density ratio of ~1-3, which is consistent with the values measured in the shocked outflow of the low-mass protostar L1157, the massive star-forming regions W51 and W3(OH), and the galactic center GMC G+0.693-0.03. This narrow range of PO/PN ratios across sources with a range of environmental conditions is surprising, and likely encodes information on how phosphorus carriers are stored in grain mantles.

It is not known whether the original carriers of Earth's nitrogen were molecular ices or refractory dust. To investigate this question, we have used data and results of Herschel observations towards two protostellar sources: the high-mass hot core of Orion KL, and the low-mass protostar IRAS 16293-2422. Towards Orion KL, our analysis of the molecular inventory of Crockett et al. (2014) indicates that HCN is the organic molecule that contains by far the most nitrogen, carrying $74_{-9}^{+5}\%$ of nitrogen-in-organics. Following this evidence, we explore HCN towards IRAS 16293-2422, which we consider a solar analog. Towards IRAS 16293-2422, we have reduced and analyzed Herschel spectra of HCN, and fit these observations against "jump" abundance models of IRAS 16293-2422's protostellar envelope. We find an inner-envelope HCN abundance $X_{\textrm{in}} = 5.9\pm0.7 \times 10^{-8}$ and an outer-envelope HCN abundance $X_{\textrm{out}} = 1.3 \pm 0.1 \times 10^{-9}$. We also find the sublimation temperature of HCN to be $T_{\textrm{jump}} = 71 \pm 3$~K; this measured $T_{\textrm{jump}}$ enables us to predict an HCN binding energy $E_{\textrm{B}}/k = 3840 \pm 140$~K. Based on a comparison of the HCN/H2O ratio in these protostars to N/H2O ratios in comets, we find that HCN (and, by extension, other organics) in these protostars is incapable of providing the total bulk N/H2O in comets. We suggest that refractory dust, not molecular ices, was the bulk provider of nitrogen to comets. However, interstellar dust is not known to have 15N enrichment, while high 15N enrichment is seen in both nitrogen-bearing ices and in cometary nitrogen. This may indicate that these 15N-enriched ices were an important contributor to the nitrogen in planetesimals and likely to the Earth.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

The all-Galaxy CO survey of Dame, Hartmann, & Thaddeus (2001) is by far the most uniform, large-scale Galactic CO survey. Using a dendrogram-based decomposition of this survey, we present a catalog of 1064 massive molecular clouds throughout the Galactic plane. This catalog contains $2.5 \times 10^8$ solar masses, or $25^{+10.7}_{-5.8} \%$ of the Milky Way's estimated H$_2$ mass. We track clouds in some spiral arms through multiple quadrants. The power index of Larson's first law, the size-linewidth relation, is consistent with 0.5 in all regions - possibly due to an observational bias - but clouds in the inner Galaxy systematically have significantly (~ 30%) higher linewidths at a given size, indicating that their linewidths are set in part by Galactic environment. The mass functions of clouds in the inner Galaxy versus the outer Galaxy are both qualitatively and quantitatively distinct. The inner Galaxy mass spectrum is best described by a truncated power-law with a power index of $\gamma=-1.6\pm0.1$ and an upper truncation mass $M_0 = (1.0 \pm 0.2) \times 10^7 M_\odot$, while the outer Galaxy mass spectrum is better described by a non-truncating power law with $\gamma=-2.2\pm0.1$ and an upper mass $M_0 = (1.5 \pm 0.5) \times 10^6 M_\odot$, indicating that the inner Galaxy is able to form and host substantially more massive GMCs than the outer Galaxy. Additionally, we have simulated how the Milky Way would appear in CO from extragalactic perspectives, for comparison with CO maps of other galaxies.

Using the United Kingdom Infrared Telescope on Mauna Kea, we have carried out a new near-infrared J, H, K monitoring survey of almost a square degree of the star-forming Orion Nebula Cluster with observations on 120 nights over three observing seasons, spanning a total of 894 days. We monitored ~15,000 stars down to J=20 using the WFCAM instrument, and have extracted 1203 significantly variable stars from our data. By studying variability in young stellar objects (YSOs) in the H-K, K color-magnitude diagram, we are able to distinguish between physical mechanisms of variability. Many variables show color behavior indicating either dust-extinction or disk/accretion activity, but we find that when monitored for longer periods of time, a number of stars shift between these two variability mechanisms. Further, we show that the intrinsic timescale of disk/accretion variability in young stars is longer than that of dust-extinction variability. We confirm that variability amplitude is statistically correlated with evolutionary class in all bands and colors. Our investigations of these 1203 variables have revealed 73 periodic AA Tau type variables, many large-amplitude and long-period (P > 15 day) YSOs, including three stars showing widely-spaced periodic brightening events consistent with circumbinary disk activity, and four new eclipsing binaries. These phenomena and others indicate the activity of long-term disk/accretion variability processes taking place in young stars. We have made the light curves and associated data for these 1203 variables available online.

We present an analysis of near-infrared time-series photometry in J, H, and K bands for about 100 epochs of a 1 square degree region of the Lynds 1003/1004 dark cloud in the Cygnus OB7 region. Augmented by data from the Wide-field Infrared Survey Explorer (WISE), we identify 96 candidate disk bearing young stellar objects (YSOs) in the region. Of these, 30 are clearly Class I or earlier. Using the Wide-Field imaging CAMera (WFCAM) on the United Kingdom InfraRed Telescope (UKIRT), we were able to obtain photometry over three observing seasons, with photometric uncertainty better than 0.05 mag down to J ~17. We study detailed light curves and color trajectories of ~50 of the YSOs in the monitored field. We investigate the variability and periodicity of the YSOs and find the data are consistent with all YSOs being variable in these wavelengths on time scales of a few years. We divide the variability into four observational classes: 1) stars with periodic variability stable over long timescales, 2) variables which exhibit short-lived cyclic behavior, 3) long duration variables, and 4) stochastic variables. Some YSO variability defies simple classification. We can explain much of the observed variability as being due to dynamic and rotational changes in the disk, including an asymmetric or changing blocking fraction, changes to the inner disk hole size, as well as changes to the accretion rate. Overall, we find that the Class I:Class II ratio of the cluster is consistent with an age of < 1Myr, with at least one individual, wildly varying, source ~ 100,000 yr old. We have also discovered a Class II eclipsing binary system with a period of 17.87 days.

We present a subset of the results of a three season, 124 night, near-infrared monitoring campaign of the dark clouds Lynds 1003 and Lynds 1004 in the Cygnus OB7 star forming region. In this paper, we focus on the field star population. Using three seasons of UKIRT J, H and K band observations spanning 1.5 years, we obtained high-quality photometry on 9,200 stars down to J=17 mag, with photometric uncertainty better than 0.04 mag. After excluding known disk bearing stars we identify 149 variables - 1.6% of the sample. Of these, about 60 are strictly periodic, with periods predominantly < 2 days. We conclude this group is dominated by eclipsing binaries. A few stars have long period signals of between 20 and 60 days. About 25 stars have weak modulated signals, but it was not clear if these were periodic. Some of the stars in this group may be diskless young stellar objects with relatively large variability due to cool star spots. The remaining ~60 stars showed variations which appear to be purely stochastic.

We present the first results from a 124 night J, H, K near-infrared monitoring campaign of the dark cloud L 1003 in Cygnus OB7, an active star-forming region. Using 3 seasons of UKIRT observations spanning 1.5 years, we obtained high-quality photometry on 9,200 stars down to J=17 mag, with photometric uncertainty better than 0.04 mag. On the basis of near-infrared excesses from disks, we identify 30 pre-main sequence stars, including 24 which are newly discovered. We analyze those stars and find the NIR excesses are significantly variable. All 9,200 stars were monitored for photometric variability; among the field star population, about 160 exhibited near-infrared variability (1.7% of the sample). Of the 30 YSOs (young stellar objects), 28 of them (93%) are variable at a significant level. 25 of the 30 YSOs have near-infrared excess consistent with simple disk-plus-star classical T Tauri models. Nine of these (36%) drift in color space over the course of these observations and/or since 2MASS observations such that they cross the boundary defining the NIR excess criteria; effectively, they have a transient near-infrared excess. About half of the YSOs have color-space variations parallel to either the classical T Tauri star locus or a hybrid track which includes the dust reddening trajectory. This indicates that the NIR variability in YSOs that possess accretion disks arises from a combination of variable extinction and changes in the inner accretion disk: either in accretion rate, central hole size and/or the inclination of the inner disk. While some variability may be due to stellar rotation, the level of variability on the individual stars can exceed a magnitude. This is a strong empirical suggestion that protoplanetary disks are quite dynamic and exhibit more complex activity on short timescales than is attributable to rotation alone or captured in static disk models.