This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

We report ALMA and VLA continuum observations that potentially identify the four youngest protostars in the Orion Molecular Clouds taken as part of the Orion VANDAM program. These are distinguished by bright, extended, irregular emission at 0.87 mm and 8 mm and are optically thick at 0.87 mm. These structures are distinct from the disk or point-like morphologies seen toward the other Orion protostars. The 0.87 mm emission implies temperatures of 41-170 K, requiring internal heating. The bright 8 mm emission implies masses of 0.5 to 1.2 M_sun assuming standard dust opacity models. One source has a Class 0 companion, while another exhibits substructure indicating a companion-candidate. Three compact outflows are detected, two of which may be driven by companions, with dynamical times of ~300 to ~400 years. The slowest outflow may be driven by a first hydrostatic core. These protostars appear to trace an early phase when the centers of collapsing fragments become optically thick to their own radiation and compression raises the gas temperature. This phase is thought to accompany the formation of hydrostatic cores. A key question is whether these structures are evolving on free fall times of ~100 years, or whether they are evolving on Kelvin-Helmholtz times of several thousand years. The number of these sources imply a lifetime of ~6000 years, in closer agreement with the Kelvin-Helmholtz time. In this case, rotational and/or magnetic support could be slowing the collapse.

We have conducted a survey of 328 protostars in the Orion molecular clouds with ALMA at 0.87 mm at a resolution of $\sim$0.1" (40 au), including observations with the VLA at 9 mm toward 148 protostars at a resolution of $\sim$0.08" (32 au). This is the largest multi-wavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 mm and 9 mm to measure the dust disk radii and masses toward the Class 0, Class I, and Flat Spectrum protostars, characterizing the evolution of these disk properties in the protostellar phase. The mean dust disk radii for the Class 0, Class I, and Flat Spectrum protostars are 44.9$^{+5.8}_{-3.4}$, 37.0$^{+4.9}_{-3.0}$, and 28.5$^{+3.7}_{-2.3}$ au, respectively, and the mean protostellar dust disk masses are 25.9$^{+7.7}_{-4.0}$, 14.9$^{+3.8}_{-2.2}$, 11.6$^{+3.5}_{-1.9}$ Earth masses, respectively. The decrease in dust disk masses is expected from disk evolution and accretion, but the decrease in disk radii may point to the initial conditions of star formation not leading to the systematic growth of disk radii or that radial drift is keeping the dust disk sizes small. At least 146 protostellar disks (35% out of 379 detected 0.87 mm continuum sources plus 42 non-detections) have disk radii greater than 50 au in our sample. These properties are not found to vary significantly between different regions within Orion. The protostellar dust disk mass distributions are systematically larger than that of Class II disks by a factor of $>$4, providing evidence that the cores of giant planets may need to at least begin their formation during the protostellar phase.

We present ALMA (0.87~mm) and VLA (9~mm) observations toward OMC2-FIR4 and OMC2-FIR3 within the Orion integral-shaped filament that are thought to be the nearest regions of intermediate mass star formation. We characterize the continuum sources within these regions on $\sim$40~AU (0\farcs1) scales and associated molecular line emission at a factor of $\sim$30 better resolution than previous observations at similar wavelengths. We identify six compact continuum sources within OMC2-FIR4, four in OMC2-FIR3, and one additional source just outside OMC2-FIR4. This continuum emission is tracing the inner envelope and/or disk emission on less than 100~AU scales. HOPS-108 is the only protostar in OMC2-FIR4 that exhibits emission from high-excitation transitions of complex organic molecules (e.g., methanol and other lines) coincident with the continuum emission. HOPS-370 in OMC2-FIR3 with L~$\sim$~360~\lsun, also exhibits emission from high-excitation methanol and other lines. The methanol emission toward these two protostars is indicative of temperatures high enough to thermally evaporate methanol from icy dust grains; overall these protostars have characteristics similar to hot corinos. We do not identify a clear outflow from HOPS-108 in \twco, but find evidence of interaction between the outflow/jet from HOPS-370 and the OMC2-FIR4 region. The multitude of observational constraints indicate that HOPS-108 is likely a low to intermediate-mass protostar in its main mass accretion phase and it is the most luminous protostar in OMC2-FIR4. The high resolution data presented here are essential for disentangling the embedded protostars from their surrounding dusty environments and characterizing them.

We perform simulations to test the effects of a moving gas filament on a young star cluster (i.e. the "Slingshot" Model). We model Orion Nebula Cluster-like clusters as Plummer spheres and the Integral Shaped Filament gas as a cylindrical potential. We observe that in a static filament, an initially spherical cluster evolves naturally into an elongated distribution of stars. For sinusoidal moving filaments, we observe different remnants, and classify them into 4 categories.%: 3 different objects and one transition object. "Healthy" clusters, where almost all the stars stay inside the filament and the cluster; "destroyed" clusters are the opposite case, with almost no particles in the filament or near the centre of density of the clusters; "ejected" clusters, where a large fraction of stars are close to the centre of density of the stars , but almost none of them in the filament; and "transition" clusters, where roughly the same number of particles is ejected from the cluster and from the filament. An {{Orion Nebula Cluster-like}} cluster might stay inside the filament or be ejected, but it will not be destroyed.

We present analysis of the gas kinematics of the Integral Shaped Filament (ISF) in Orion~A using four different molecular lines, $^{12}$CO (1-0), $^{13}$CO (1-0), NH$_3$ (1,1), and N$_2$H$^+$ (1-0). We describe our method to visualize the position-velocity (PV) structure using the intensity-weighted line velocity centroid, which enables us to identify structures that were previously muddled or invisible. We observe a north to south velocity gradient in all tracers that terminates in a velocity peak near the center of the Orion Nebula Cluster (ONC), consistent with the previously reported "wave-like" properties of the ISF. We extract the velocity dispersion profiles and compare the non-thermal line widths to the gas gravitational potential. We find supersonic Mach number profiles, yet the line widths are consistent with the gas being deeply gravitationally bound. We report the presence of two $^{12}$CO velocity components along the northern half of the ISF; if interpreted as circular rotation, the angular velocity is $\omega=1.4\,{\rm Myr}^{-1}$. On small scales we report the detection of N$_2$H$^+$ and NH$_3$ "twisting and turning" structures, with short associated timescales that give the impression of a torsional wave. Neither the nature of these structures nor their relation to the larger scale wave is presently understood.

We have investigated the formation and kinematics of sub-mm continuum cores in the Orion A molecular cloud. A comparison between sub-mm continuum and near infrared extinction shows a continuum core detection threshold of $A_V\sim$ 5-10 mag. The threshold is similar to the star formation extinction threshold of $A_V\sim$ 7 mag proposed by recent work, suggesting a universal star formation extinction threshold among clouds within 500 pc to the Sun. A comparison between the Orion A cloud and a massive infrared dark cloud G28.37+0.07 indicates that Orion A produces more dense gas within the extinction range 15 mag $\lesssim A_V \lesssim$ 60 mag. Using data from the CARMA-NRO Orion Survey, we find that dense cores in the integral-shaped filament (ISF) show sub-sonic core-to-envelope velocity dispersion that is significantly less than the local envelope line dispersion, similar to what has been found in nearby clouds. Dynamical analysis indicates that the cores are bound to the ISF. An oscillatory core-to-envelope motion is detected along the ISF. Its origin is to be further explored.

We use APEX mapping observations of 13CO, and C18O (2-1) to investigate the internal gas kinematics of the filamentary cloud G350.54+0.69, composed of the two distinct filaments G350.5-N and G350.5-S. G350.54+0.69 as a whole is supersonic and gravitationally bound. We find a large-scale periodic velocity oscillation along the entire G350.5-N filament with a wavelength of ~1.3 pc and an amplitude of ~0.12 km/s. Comparing with gravitational-instability induced core formation models, we conjecture that this periodic velocity oscillation could be driven by a combination of longitudinal gravitational instability and a large-scale periodic physical oscillation along the filament. The latter may be an example of an MHD transverse wave. This hypothesis can be tested with Zeeman and dust polarization measurements.

The Integral Shaped Filament (ISF) is the nearest molecular cloud with rapid star formation, including massive stars, and it is therefore a star-formation laboratory. We use Gaia parallaxes, to show that the distances to young Class II stars ('disks') projected along the spine of this filament are related to the gas radial velocity by $$ v = -{D\over\tau} + K;\qquad \tau = 4\,{\rm Myr}, $$ where $K$ is a constant. This implies that the ISF is a standing wave, which is consistent with the Stutz & Gould (2016) 'Slingshot' prediction. The $\tau=4\,{\rm Myr}$ timescale is consistent with the 'Slingshot' picture that the Orion Nebula Cluster (ONC) is the third cluster to be violently split off from the Orion A cloud (following NGC 1981 and NGC 1987) at few-Myr intervals due to gravito-magnetic oscillations. We also present preliminary evidence that the truncation of the ISF is now taking place $16^\prime$ south of the ONC and is mediated by a torsional wave that is propagating south with a characteristic timescale $\tau_{\rm torsion} = 0.5\,{\rm Myr}$, i.e. eight times shorter. The relation between these two wave phenomena is not presently understood.

We analyze Herschel Space Observatory observations of 104 young stellar objects with protoplanetary disks in the ~1.5 Myr star-forming region Lynds 1641 (L1641) within the Orion A Molecular Cloud. We present spectral energy distributions from the optical to the far-infrared including new photometry from the Herschel Photodetector Array Camera and Spectrometer (PACS) at 70 microns. Our sample, taken as part of the Herschel Orion Protostar Survey, contains 24 transitional disks, eight of which we identify for the first time in this work. We analyze the full disks with irradiated accretion disk models to infer dust settling properties. Using forward modeling to reproduce the observed nKS-[70] index for the full disk sample, we find the observed disk indices are consistent with models that have depletion of dust in the upper layers of the disk relative to the mid plane, indicating significant dust settling. We perform the same analysis on full disks in Taurus with Herschel data and find that Taurus is slightly more evolved, although both samples show signs of dust settling. These results add to the growing literature that significant dust evolution can occur in disks by ~1.5 Myr.

We use the POLARIS radiative transport code to generate predictions of the two main observables directly sensitive to the magnetic field morphology and strength in filaments: dust polarization and gas Zeeman line splitting. We simulate generic gas filaments with power-law density profiles assuming two density-field strength dependencies, six different filament inclinations, and nine distinct magnetic field morphologies, including helical, toroidal, and warped magnetic field geometries. We present idealized spatially resolved dust polarization and Zeeman-derived field strengths and directions maps. Under the assumption that dust grains are aligned by radiative torques (RATs), dust polarization traces the projected plane-of-the-sky magnetic field morphology. Zeeman line splitting delivers simultaneously the intensity-weighted line-of-sight field strength and direction. We show that linear dust polarization alone is unable to uniquely constrain the 3D field morphology. We demonstrate that these ambiguities are ameliorated or resolved with the addition of the Zeeman directional information. Thus, observations of both the dust polarization and Zeeman splitting together provide the most promising means for obtaining constraints of the 3D magnetic field configuration. We find that the Zeeman-derived field strengths are at least a factor of a few below the input field strengths due to line-of-sight averaging through the filament density gradient. Future observations of both dust polarization and Zeeman splitting are essential for gaining insights into the role of magnetic fields in star and cluster forming filaments.

We investigate the global properties of the straight and isolated filamentary cloud G350.54+0.69 using Herschel continuum and APEX molecular line data. The overall straight morphology is similar to two other well studied nearby filaments (Musca and Taurus-B211/3) while the isolated nature of G350.54+0.69 appears similar to Musca. G350.54+0.69 is composed of two distinct filaments with a length ~5.9 pc for G350.5-N (~2.3 pc for G350.5-S), a total mass of ~810 $M_{\odot}$ (~ 110 $M_{\odot}$), and a mean temperature of ~ 18.2 K (~17.7 K). We identify 9 dense and gravitationally bound cores in the whole cloud G350.54+0.69. The separations between cores and the line mass of the whole cloud appear to follow the predictions of the "sausage" instability theory, which suggests that G350.54+0.69 could have undergone radial collapse and fragmentation. The presence of young protostars is consistent with this hypothesis. The line masses of the two filaments (~120 $M_{\odot}$ pc$^{-1}$ for G350.5-N, and ~45 $M_{\odot}$ pc$^{-1}$ for G350.5-S), mass-size distributions of the dense cores, and low-mass protostars collectively suggest that G350.54+0.69 is a site of ongoing low-mass star formation. Based on the above evidence, we place G350.54+0.69 in an intermediate evolutionary state between Musca and Taurus-B211/3. We suggest that investigations into straight (and isolated) versus those distributed inside molecular clouds may provide important clues into filament formation and evolution.

Models aiming to explain the formation of massive black hole seeds, and in particular the direct collapse scenario, face substantial difficulties. These are rooted in rather ad hoc and fine-tuned initial conditions, such as the simultaneous requirements of extremely low metallicities and strong radiation backgrounds. Here we explore a modification of such scenarios where a massive primordial star cluster is initially produced. Subsequent stellar collisions give rise to the formation of massive (10^4 - 10^5 solar mass) objects. Our calculations demonstrate that the interplay between stellar dynamics, gas accretion and protostellar evolution is particularly relevant. Gas accretion onto the protostars enhances their radii, resulting in an enhanced collisional cross section. We show that the fraction of collisions can increase from 0.1-1% of the initial population to about 10% when compared to gas-free models or models of protostellar clusters in the local Universe. We conclude that very massive objects can form in spite of initial fragmentation, making the first massive protostellar clusters viable candidate birth places for observed supermassive black holes.

The nature of the far-IR source HOPS 108 has been a matter of debate in the last years. Previous radio observations detected a 3.6 cm source (VLA 12), coincident with HOPS 108, that was interpreted as a radio jet from this protostar. We present new multi-wavelength (0.7-5 cm), multi-configuration VLA observations as well as archive data (3.6 cm) that reveal VLA 12 as three knots of non-thermal emission, with HOPS 108 close to the central knot. We show that these knots have not been ejected by HOPS 108. We propose that the VLA 12 knots are actually part of a radio jet driven by VLA 11 (HOPS 370), a strong nearby source clearly elongated in the direction of the knots. The position of HOPS 108 in the path of the VLA 11-VLA 12 jet suggests an appealing new scenario: the triggered formation of HOPS 108 by the interaction of the jet with the surrounding medium.

We present intensity-corrected Herschel maps at 100 um, 160 um, 250 um, 350 um, and 500 um for 56 isolated low-mass clouds. We determine the zero-point corrections for Herschel PACS and SPIRE maps from the Herschel Science Archive (HSA) using Planck data. Since these HSA maps are small, we cannot correct them using typical methods. Here, we introduce a technique to measure the zero-point corrections for small Herschel maps. We use radial profiles to identify offsets between the observed HSA intensities and the expected intensities from Planck. Most clouds have reliable offset measurements with this technique. In addition, we find that roughly half of the clouds have underestimated HSA-SPIRE intensities in their outer envelopes relative to Planck, even though the HSA-SPIRE maps were previously zero-point corrected. Using our technique, we produce corrected Herschel intensity maps for all 56 clouds and determine their line-of-sight average dust temperatures and optical depths from modified black body fits. The clouds have typical temperatures of ~ 14-20 K and optical depths of ~ 1e-5 - 1e-3. Across the whole sample, we find an anti-correlation between temperature and optical depth. We also find lower temperatures than what was measured in previous Herschel studies, which subtracted out a background level from their intensity maps to circumvent the zero-point correction. Accurate Herschel observations of clouds are key to obtain accurate density and temperature profiles. To make such future analyses possible, intensity-corrected maps for all 56 clouds are publicly available in the electronic version.

The Orion nebula is a prime example of a massive star-forming region in our galaxy. Observations have shown that gravitational and magnetic energy are comparable in its integral shaped filament (ISF) on a scale of ~1 pc, and that the population of pre-main sequence stars appears dynamically heated compared to the protostars. These results have been attributed to a slingshot mechanism resulting from the oscillation of the filament (Stutz & Gould 2016). In this paper, we show that radially contracting filaments naturally evolve toward a state where gravitational, magnetic, and rotational energy are comparable. While the contraction of the filament will preferentially amplify the axial component of the magnetic field, the presence of rotation leads to a helical field structure. We show how magnetic tension can give rise to a filament oscillation, and estimate a typical timescale of 0.7 million years for the motion of the filament to the position of maximum displacement, consistent with the characteristic timescale of the ejected stars. Furthermore, the presence of helical magnetic fields is expected to give rise to magneto-hydrodynamical instabilities. We show here that the presence of a magnetic field significantly enhances the overall instability, which operates on a characteristic scale of about 1 pc. We expect the physics discussed here to be generally relevant in massive star forming regions, and encourage further investigations in the future.

We characterize the stellar and gas volume density, potential, and gravitational field profiles in the central $\sim$ 0.5 pc of the Orion Nebula Cluster (ONC), the nearest embedded star cluster (or rather, proto-cluster) hosting massive star formation available for detailed observational scrutiny. We find that the stellar volume density is well characterized by a Plummer profile $\rho_{stars}(r) = 5755\,{\rm M}_{\odot}\,{\rm pc}^{-3}\,(1+(r/a)^2)^{-5/2}$, where $a = 0.36$ pc. The gas density follows a cylindrical power law $\rho_{gas}(R) = 25.9\,{\rm M}_{\odot}\,{\rm pc}^{-3}\,(R/{\rm pc})^{-1.775}$. The stellar density profile dominates over the gas density profile inside $r\,\sim\,1$ pc. The gravitational field is gas-dominated at all radii, but the contribution to the total field by the stars is nearly equal to that of the gas at $r\,\sim\,a$. This fact alone demonstrates that the proto-cluster cannot be considered a gas-free system or a virialized system dominated by its own gravity. The stellar proto-cluster core is dynamically young, with an age of $\sim$ 2-3 Myr, a 1D velocity dispersion of $\sigma_{\rm obs} = 2.6$ km s$^{-1}$, and a crossing time of $\sim$ 0.55 Myr. This timescale is almost identical to the gas filament oscillation timescale estimated recently by Stutz & Gould (2016). This provides strong evidence that the proto-cluster structure is regulated by the gas filament. The proto-cluster structure may be set by tidal forces due to the oscillating filamentary gas potential. Such forces could naturally suppress low density stellar structures on scales $\gtrsim\,a$. The analysis presented here leads to a new suggestion that clusters form by an analog of the "slingshot mechanism" previously proposed for stars.

We present the GIGAS (Gigahertz Identification of Giant Air Shower) microwave radio sensor arrays of the EASIER project (Extensive Air Shower Identification with Electron Radiometers), deployed at the site of the Pierre Auger cosmic ray observatory. The aim of these novel arrays is to probe the intensity of the molecular bremsstrahlung radiation expected from the development of the extensive air showers produced by the interaction of ultra high energy cosmic rays in the atmosphere. In the designed setup, the sensors are embedded within the surface detector array of the Pierre Auger observatory allowing us to use the particle signals at ground level to trigger the radio system. A series of seven, then 61 sensors have been deployed in the C-band, followed by a new series of 14 higher sensitivity ones in the C-band and the L-band. The design, the operation, the calibration and the sensitivity to extensive air showers of these arrays are described in this paper.

Density profiles of isolated cores derived from thermal dust continuum emission rely on models of dust properties, such as mass opacity, which are poorly constrained. With complementary measures from near-infrared extinction maps, we can assess the reliability of commonly-used dust models. In this work, we compare Herschel-derived maps of the optical depth with equivalent maps derived from CFHT WIRCAM near-infrared observations for three isolated cores: CB68, L429, and L1552. We assess the dust opacities provided from four models: OH1a, OH5a, Orm1, and Orm4. Although the consistency of the models differs between the three sources, the results suggest that the optical properties of dust in the envelopes of the cores are best described by either silicate and bare graphite grains (e.g., Orm1) or carbonaceous grains with some coagulation and either thin or no ice mantles (e.g., OH5a). None of the models, however, individually produced the most consistent optical depth maps for every source. The results suggest that either the dust in the cores is not well described by any one dust property model, the application of the dust models cannot be extended beyond the very center of the cores, or more complex SED fitting functions are necessary.

Deuterium fractionation is considered as an important process to infer the chemical ages of prestellar cores in filaments. We present here the first magneto-hydrodynamical simulations including a chemical network to study deuterium fractionation in magnetized and turbulent filaments and their substructures. The filaments typically show widespread deuterium fractionation with average values $\gtrsim0.01$. For individual cores of similar age, we observe the deuteration fraction to increase with time, but also to be independent of their average properties such as density, virial or mass-to-magnetic flux ratio. We further find a correlation of the deuteration fraction with core mass, average H$_2$ density and virial parameter only at late evolutionary stages of the filament and attribute this to the lifetime of the individual cores. Specifically, chemically old cores reveal higher deuteration fractions. Within the radial profiles of selected cores, we notice differences in the structure of the deuteration fraction or surface density, which we can attribute to their different turbulent properties. High deuteration fractions of the order $0.01-0.1$ may be reached within approximately $200$~kyrs, corresponding to two free-fall times, as defined for cylindrical systems, of the filaments

Observations of the Orion-A integral shaped filament (ISF) have shown indications of an oscillatory motion of the gas filament. This evidence is based on both the wave-like morphology of the filament as well as the kinematics of the gas and stars, where the characteristic velocities of the stars require a dynamical heating mechanism. As proposed by Stutz and Gould (2016), such a heating mechanism (the "Slingshot") may be the result of an oscillating gas filament in a gas-dominated (as opposed to stellar-mass dominated) system. Here we test this hypothesis with the first stellar-dynamical simulations in which the stars are subjected to the influence of an oscillating cylindrical potential. The accelerating, cylindrical background potential is populated with a narrow distribution of stars. By coupling the potential to N-body dynamics, we are able to measure the influence of the potential on the stellar distribution. The simulations provide evidence that the slingshot mechanism can successfully reproduce several stringent observational constraints. These include the stellar spread (both in projected position and in velocity) around the filament, the symmetry in these distributions, and a bulk motion of the stars with respect to the filament. Using simple considerations we show that star-star interactions are incapable of reproducing these spreads on their own when properly accounting for the gas potential. Thus, properly accounting for the gas potential is essential for understanding the dynamical evolution of star forming filamentary systems in the era of Gaia.

Herschel has revolutionized our ability to measure column densities (N$_{\rm H}$) and temperatures (T) of molecular clouds thanks to its far infrared multiwavelength coverage. However, the lack of a well defined background intensity level in the Herschel data limits the accuracy of the N$_{\rm H}$ and T maps. We provide a method that corrects the missing Herschel background intensity levels using the Planck model for foreground Galactic thermal dust emission. We present a Fourier method that combines the publicly available Planck model on large angular scales with the Herschel images on smaller angular scales. We apply our method to two regions spanning a range of Galactic environments: Perseus and the Galactic plane region around $l = 11\deg$ (HiGal--11). We post-process the combined dust continuum emission images to generate column density and temperature maps. We compare these to previously adopted constant--offset corrections. We find significant differences ($\gtrsim$20\%) over significant ($\sim$15\%) areas of the maps, at low column densities ($N_{\rm H}\lesssim10^{22}$\,cm$^{-2}$) and relatively high temperatures ($T\gtrsim20$\,K). We also apply our method to synthetic observations of a simulated molecular cloud to validate our method. Our method successfully corrects the Herschel images, including both the constant--offset intensity level and the scale-dependent background variations measured by Planck. Our method improves the previous constant--offset corrections, which did not account for variations in the background emission levels.

The Herschel Orion Protostar Survey obtained well-sampled 1.2 - 870 micron spectral energy distributions (SEDs) of over 300 protostars in the Orion molecular clouds, home to most of the young stellar objects (YSOs) in the nearest 500 pc. We plot the bolometric luminosities and temperatures for 330 Orion YSOs, 315 of which have bolometric temperatures characteristic of protostars. The histogram of bolometric temperature is roughly flat; 29% of the protostars are in Class 0. The median luminosity decreases by a factor of four with increasing bolometric temperature; consequently, the Class 0 protostars are systematically brighter than the Class I protostars, with a median luminosity of 2.3 L_sun as opposed to 0.87 L_sun. At a given bolometric temperature, the scatter in luminosities is three orders of magnitude. Using fits to the SEDs, we analyze how the luminosities corrected for inclination and foreground reddening relate to the mass in the inner 2500 AU of the best-fit model envelopes. The histogram of envelope mass is roughly flat, while the median corrected luminosity peaks at 15 L_sun for young envelopes and falls to 1.7 L_sun for late-stage protostars with remnant envelopes. The spread in luminosity at each envelope mass is three orders of magnitude. Envelope masses that decline exponentially with time explain the flat mass histogram and the decrease in luminosity, while the formation of a range of stellar masses explains the dispersion in luminosity.

We carried out multiwavelength (0.7-5 cm), multiepoch (1994-2015) Very Large Array (VLA) observations toward the region enclosing the bright far-IR sources FIR 3 (HOPS 370) and FIR 4 (HOPS 108) in OMC-2. We report the detection of 10 radio sources, seven of them identified as young stellar objects. We image a well-collimated radio jet with a thermal free-free core (VLA 11) associated with the Class I intermediate-mass protostar HOPS 370. The jet presents several knots (VLA 12N, 12C, 12S) of non-thermal radio emission (likely synchrotron from shock-accelerated relativistic electrons) at distances of ~7,500-12,500 au from the protostar, in a region where other shock tracers have been previously identified. These knots are moving away from the HOPS 370 protostar at ~ 100 km/s. The Class 0 protostar HOPS 108, which itself is detected as an independent, kinematically decoupled radio source, falls in the path of these non-thermal radio knots. These results favor the previously proposed scenario where the formation of HOPS 108 has been triggered by the impact of the HOPS 370 outflow with a dense clump. However, HOPS 108 presents a large proper motion velocity of ~ 30 km/s, similar to that of other runaway stars in Orion, whose origin would be puzzling within this scenario. Alternatively, an apparent proper motion could result because of changes in the position of the centroid of the source due to blending with nearby extended emission, variations in the source shape, and /or opacity effects.

We study the fragmentation of the nearest high line-mass filament, the integral shaped filament (ISF, line-mass $\sim$ 400 M$_\odot$ pc$^{-1}$) in the Orion A molecular cloud. We have observed a 1.6 pc long section of the ISF with the Atacama Large Millimetre/submillimeter Array (ALMA) at 3 mm continuum emission, at a resolution of $\sim$3" (1 200 AU). We identify from the region 43 dense cores with masses about a solar mass. 60% of the ALMA cores are protostellar and 40\% are starless. The nearest neighbour separations of the cores do not show a preferred fragmentation scale; the frequency of short separations increases down to 1 200 AU. We apply a two-point correlation analysis on the dense core separations and show that the ALMA cores are significantly grouped at separations below $\sim$17 000 AU and strongly grouped below $\sim$6 000 AU. The protostellar and starless cores are grouped differently: only the starless cores group strongly below $\sim$6 000 AU. In addition, the spatial distribution of the cores indicates periodic grouping of the cores into groups of $\sim$30 000 AU in size, separated by $\sim$50 000 AU. The groups coincide with dust column density peaks detected by Herschel. These results show hierarchical, two-mode fragmentation in which the maternal filament periodically fragments into groups of dense cores. Critically, our results indicate that the fragmentation models for lower line-mass filaments ($\sim$ 16 M$_\odot$ pc$^{-1}$) fail to capture the observed properties of the ISF. We also find that the protostars identified with Spitzer and Herschel in the ISF are grouped at separations below $\sim$17 000 AU. In contrast, young stars with disks do not show significant grouping. This suggests that the grouping of dense cores is partially retained over the protostar lifetime, but not over the lifetime of stars with disks.

Ultrahigh energy cosmic ray air showers probe particle physics at energies beyond the reach of accelerators. Here we introduce a new method to test hadronic interaction models without relying on the absolute energy calibration, and apply it to events with primary energy 6-16 EeV (E_CM = 110-170 TeV), whose longitudinal development and lateral distribution were simultaneously measured by the Pierre Auger Observatory. The average hadronic shower is 1.33 +- 0.16 (1.61 +- 0.21) times larger than predicted using the leading LHC-tuned models EPOS-LHC (QGSJetII-04), with a corresponding excess of muons.

We present the first detection of a jet in the far-IR [O I] lines from an intermediate mass protostar. We have carried out a Herschel/PACS spectral mapping study in the [O I] lines of OMC-2 FIR 3 and FIR 4, two of the most luminous protostars in Orion outside of the Orion Nebula. The spatial morphology of the fine structure line emission reveals the presence of an extended photodissociation region (PDR) and a narrow, but intense jet connecting the two protostars. The jet seen in [O I] emission is spatially aligned with the Spitzer/IRAC 4.5 micron jet and the CO (6-5) molecular outflow centered on FIR 3. The mass loss rate derived from the total [O I] 63 micron line luminosity of the jet is 7.7 x 10^-6 M_sun/yr, more than an order of magnitude higher than that measured for typical low mass class 0 protostars. The implied accretion luminosity is significantly higher than the observed bolometric luminosity of FIR 4, indicating that the [O I] jet is unlikely to be associated with FIR 4. We argue that the peak line emission seen toward FIR 4 originates in the terminal shock produced by the jet driven by FIR 3. The higher mass-loss rate that we find for FIR 3 is consistent with the idea that intermediate mass protostars drive more powerful jets than their low-mass counterparts. Our results also call into question the nature of FIR 4.

We investigate the evolution of far-IR CO emission from protostars observed with Herschel/PACS for 50 sources from the combined sample of HOPS and DIGIT Herschel key programs. From the uniformly sampled spectral energy distributions, we computed $L_{\rm{bol}}$, $T_{\rm{bol}}$ and $L_{\rm {bol}}/L_{\rm {smm}}$ for these sources to search for correlations between far-IR CO emission and protostellar properties. We find a strong and tight correlation between far-IR CO luminosity ($L^{\rm fir}_{\rm CO}$) and the bolometric luminosity ($L_{\rm{bol}}$) of the protostars with $L^{\rm fir}_{\rm CO}$ $\propto$ $L_{\rm{bol}}^{0.7}$. We, however, do not find a strong correlation between $L^{\rm fir}_{\rm CO}$ and protostellar evolutionary indicators, $T_{\rm{bol}}$ and $L_{\rm {bol}}/L_{\rm {smm}}$. FIR CO emission from protostars traces the currently shocked gas by jets/outflows, and $L^{\rm fir}_{\rm CO}$ is proportional to the instantaneous mass loss rate, $\dot{M}_{\rm{out}}$. The correlation between $L^{\rm fir}_{\rm CO}$ and $L_{\rm{bol}}$ is indicative of instantaneous $\dot{M}_{\rm{out}}$ tracking instantaneous $\dot{M}_{\rm{acc}}$. The lack of correlation between $L^{\rm fir}_{\rm CO}$ and evolutionary indicators $T_{\rm{bol}}$ and $L_{\rm {bol}}/L_{\rm {smm}}$ suggests that $\dot{M}_{\rm{out}}$ and, therefore, $\dot{M}_{\rm{acc}}$ do not show any clear evolutionary trend. These results are consistent with mass accretion/ejection in protostars being episodic. Taken together with the previous finding that the time-averaged mass ejection/accretion rate declines during the protostellar phase (e.g., Bontemps et al. 1996), our results suggest that the instantaneous accretion/ejection rate of protostars is highly time variable and episodic, but the amplitude and/or frequency of this variability decreases with time such that the time averaged accretion/ejection rate declines with system age.

We present CARMA CO (J=1-0) observations and Herschel PACS spectroscopy, characterizing the outflow properties toward extremely young and deeply embedded protostars in the Orion molecular clouds. The sample comprises a subset of the Orion protostars known as the PACS Bright Red Sources (PBRS) (Stutz et al. 2013). We observed 14 PBRS with CARMA and 8 of these 14 with Herschel, acquiring full spectral scans from 55 micron to 200 micron. Outflows are detected in CO (J=1-0) from 8 of 14 PBRS, with two additional tentative detections; outflows are also detected from the outbursting protostar HOPS 223 (V2775 Ori) and the Class I protostar HOPS 68. The outflows have a range of morphologies, some are spatially compact, <10000 AU in extent, while others extend beyond the primary beam. The outflow velocities and morphologies are consistent with being dominated by intermediate inclination angles (80 deg > i > 20 deg). This confirms the interpretation of the very red 24 micron to 70 micron colors of the PBRS as a signpost of high envelope densities, with only one (possibly two) cases of the red colors resulting from edge-on inclinations. We detect high-J (J_up > 13) CO lines and/or H_2O lines from 5 of 8 PBRS and only for those with detected CO outflows. The far-infrared CO rotation temperatures of the detected PBRS are marginally colder (~230 K) than those observed for most protostars (~300 K), and only one of these 5 PBRS has detected [OI] 63 micron emission. The high envelope densities could be obscuring some [OI] emission and cause a ~20 K reduction to the CO rotation temperatures.

We measure the energy emitted by extensive air showers in the form of radio emission in the frequency range from 30 to 80 MHz. Exploiting the accurate energy scale of the Pierre Auger Observatory, we obtain a radiation energy of 15.8 \pm 0.7 (stat) \pm 6.7 (sys) MeV for cosmic rays with an energy of 1 EeV arriving perpendicularly to a geomagnetic field of 0.24 G, scaling quadratically with the cosmic-ray energy. A comparison with predictions from state-of-the-art first-principle calculations shows agreement with our measurement. The radiation energy provides direct access to the calorimetric energy in the electromagnetic cascade of extensive air showers. Comparison with our result thus allows the direct calibration of any cosmic-ray radio detector against the well-established energy scale of the Pierre Auger Observatory.

The Auger Engineering Radio Array (AERA) is part of the Pierre Auger Observatory and is used to detect the radio emission of cosmic-ray air showers. These observations are compared to the data of the surface detector stations of the Observatory, which provide well-calibrated information on the cosmic-ray energies and arrival directions. The response of the radio stations in the 30 to 80 MHz regime has been thoroughly calibrated to enable the reconstruction of the incoming electric field. For the latter, the energy deposit per area is determined from the radio pulses at each observer position and is interpolated using a two-dimensional function that takes into account signal asymmetries due to interference between the geomagnetic and charge-excess emission components. The spatial integral over the signal distribution gives a direct measurement of the energy transferred from the primary cosmic ray into radio emission in the AERA frequency range. We measure 15.8 MeV of radiation energy for a 1 EeV air shower arriving perpendicularly to the geomagnetic field. This radiation energy -- corrected for geometrical effects -- is used as a cosmic-ray energy estimator. Performing an absolute energy calibration against the surface-detector information, we observe that this radio-energy estimator scales quadratically with the cosmic-ray energy as expected for coherent emission. We find an energy resolution of the radio reconstruction of 22% for the data set and 17% for a high-quality subset containing only events with at least five radio stations with signal.

The azimuthal asymmetry in the risetime of signals in Auger surface detector stations is a source of information on shower development. The azimuthal asymmetry is due to a combination of the longitudinal evolution of the shower and geometrical effects related to the angles of incidence of the particles into the detectors. The magnitude of the effect depends upon the zenith angle and state of development of the shower and thus provides a novel observable, $(\sec \theta)_\mathrm{max}$, sensitive to the mass composition of cosmic rays above $3 \times 10^{18}$ eV. By comparing measurements with predictions from shower simulations, we find for both of our adopted models of hadronic physics (QGSJETII-04 and EPOS-LHC) an indication that the mean cosmic-ray mass increases slowly with energy, as has been inferred from other studies. However, the mass estimates are dependent on the shower model and on the range of distance from the shower core selected. Thus the method has uncovered further deficiencies in our understanding of shower modelling that must be resolved before the mass composition can be inferred from $(\sec \theta)_\mathrm{max}$.

The Pierre Auger Observatory has begun a major Upgrade of its already impressive capabilities, with an emphasis on improved mass composition determination using the surface detectors of the Observatory. Known as AugerPrime, the upgrade will include new 4 m$^2$ plastic scintillator detectors on top of all 1660 water-Cherenkov detectors, updated and more flexible surface detector electronics, a large array of buried muon detectors, and an extended duty cycle for operations of the fluorescence detectors. This Preliminary Design Report was produced by the Collaboration in April 2015 as an internal document and information for funding agencies. It outlines the scientific and technical case for AugerPrime. We now release it to the public via the arXiv server. We invite you to review the large number of fundamental results already achieved by the Observatory and our plans for the future.

We present key results from the Herschel Orion Protostar Survey (HOPS): spectral energy distributions (SEDs) and model fits of 330 young stellar objects, predominantly protostars, in the Orion molecular clouds. This is the largest sample of protostars studied in a single, nearby star-formation complex. With near-infrared photometry from 2MASS, mid- and far-infrared data from Spitzer and Herschel, and sub-millimeter photometry from APEX, our SEDs cover 1.2-870 $\mu$m and sample the peak of the protostellar envelope emission at ~100 $\mu$m. Using mid-IR spectral indices and bolometric temperatures, we classify our sample into 92 Class 0 protostars, 125 Class I protostars, 102 flat-spectrum sources, and 11 Class II pre-main-sequence stars. We implement a simple protostellar model (including a disk in an infalling envelope with outflow cavities) to generate a grid of 30400 model SEDs and use it to determine the best-fit model parameters for each protostar. We argue that far-IR data are essential for accurate constraints on protostellar envelope properties. We find that most protostars, and in particular the flat-spectrum sources, are well-fit. The median envelope density and median inclination angle decrease from Class 0 to Class I to flat-spectrum protostars, despite the broad range in best-fit parameters in each of the three categories. We also discuss degeneracies in our model parameters. Our results confirm that the different protostellar classes generally correspond to an evolutionary sequence with a decreasing envelope infall rate, but the inclination angle also plays a role in the appearance, and thus interpretation, of the SEDs.

By comparing 3 constituents of Orion A (gas, protostars, and pre-main-sequence stars), both morhologically and kinematically, we derive the following. The gas surface density near the integral-shaped filament (ISF) is well represented by a power law, Sigma(b)=72 Msun/pc^2(b/pc)^{-5/8} for our entire range, 0.05<b/pc<8.5, of distance from the filament ridge. Essentially all protostars lie on the ISF or other filament ridges, while almost all pre-main-sequence stars do not. Combined with the fact that protostars move <1 kms relative to the filaments while stars move several times faster, this implies that protostellar accretion is terminated by a slingshot ejection from the filaments. The ISF is the 3rd in a series of star bursts that are progressively moving south, with separations of a few Myr in time and 3 pc in space. This, combined with the filament's observed undulations (spatial and velocity), suggests that repeated propagation of transverse waves thru the filament is progressively digesting the material that formerly connected Orion A and B into stars in discrete episodes. We construct an axially symmetric gas density profile rho(r)=16 Msun/pc^3(r/pc)^{-13/8}. The model implies that the observed magnetic fields are supercritical on scales of the observed undulations, suggesting that the filament's transverse waves are magnetically induced. Because the magnetic fields are subcritical on scales of the filament on larger scales, the system as a whole is relatively stable and long lived. Protostellar ejection occurs because the gas accelerates away from the protostars, not the other way around. The model also implies that the ISF is kinematically young, which is consistent with other lines of evidence. The southern filament has a broken power law, which matches the ISF profile for 2.5<b/pc<8.5, but is shallower closer in. It is also kinematically older than the ISF.

We present new measurements of the dust emissivity index, beta, for the high-mass, star-forming OMC 2/3 filament. We combine 160-500 um data from Herschel with long-wavelength observations at 2 mm and fit the spectral energy distributions across a ~ 2 pc long, continuous section of OMC 2/3 at 15000 AU (0.08 pc) resolution. With these data, we measure beta and reconstruct simultaneously the filtered-out large-scale emission at 2 mm. We implement both variable and fixed values of beta, finding that beta = 1.7 - 1.8 provides the best fit across most of OMC 2/3. These beta values are consistent with a similar analysis carried out with filtered Herschel data. Thus, we show that beta values derived from spatial filtered emission maps agree well with those values from unfiltered data at the same resolution. Our results contradict the very low beta values (~ 0.9) previously measured in OMC 2/3 between 1.2 mm and 3.3 mm data, which we attribute to elevated fluxes in the 3.3 mm observations. Therefore, we find no evidence or rapid, extensive dust grain growth in OMC 2/3. Future studies with Herschel data and complementary ground-based long-wavelength data can apply our technique to obtain robust determinations of beta in nearby cold molecular clouds.

The dust extinction curve is a critical component of many observational programs and an important diagnostic of the physics of the interstellar medium. Here we present new measurements of the dust extinction curve and its variation towards tens of thousands of stars, a hundred-fold larger sample than in existing detailed studies. We use data from the APOGEE spectroscopic survey in combination with ten-band photometry from Pan-STARRS1, 2MASS, and WISE. We find that the extinction curve in the optical through infrared is well characterized by a one-parameter family of curves described by R(V). The extinction curve is more uniform than suggested in past works, with sigma(R(V)) = 0.18, and with less than one percent of sight lines having R(V) > 4. Our data and analysis have revealed two new aspects of Galactic extinction: first, we find significant, wide-area variations in R(V) throughout the Galactic plane. These variations are on scales much larger than individual molecular clouds, indicating that R(V) variations must trace much more than just grain growth in dense molecular environments. Indeed, we find no correlation between R(V) and dust column density up to E(B-V) ~ 2. Second, we discover a strong relationship between R(V) and the far-infrared dust emissivity.

This paper presents the results of different searches for correlations between very high-energy neutrino candidates detected by IceCube and the highest-energy cosmic rays measured by the Pierre Auger Observatory and the Telescope Array. We first consider samples of cascade neutrino events and of high-energy neutrino-induced muon tracks, which provided evidence for a neutrino flux of astrophysical origin, and study their cross-correlation with the ultrahigh-energy cosmic ray (UHECR) samples as a function of angular separation. We also study their possible directional correlations using a likelihood method stacking the neutrino arrival directions and adopting different assumptions on the size of the UHECR magnetic deflections. Finally, we perform another likelihood analysis stacking the UHECR directions and using a sample of through-going muon tracks optimized for neutrino point-source searches with sub-degree angular resolution. No indications of correlations at discovery level are obtained for any of the searches performed. The smallest of the p-values comes from the search for correlation between UHECRs with IceCube high-energy cascades, a result that should continue to be monitored.

A measurement of the cosmic-ray spectrum for energies exceeding $4{\times}10^{18}$ eV is presented, which is based on the analysis of showers with zenith angles greater than $60^{\circ}$ detected with the Pierre Auger Observatory between 1 January 2004 and 31 December 2013. The measured spectrum confirms a flux suppression at the highest energies. Above $5.3{\times}10^{18}$ eV, the "ankle", the flux can be described by a power law $E^{-\gamma}$ with index $\gamma=2.70 \pm 0.02 \,\text{(stat)} \pm 0.1\,\text{(sys)}$ followed by a smooth suppression region. For the energy ($E_\text{s}$) at which the spectral flux has fallen to one-half of its extrapolated value in the absence of suppression, we find $E_\text{s}=(5.12\pm0.25\,\text{(stat)}^{+1.0}_{-1.2}\,\text{(sys)}){\times}10^{19}$ eV.

Neutrinos in the cosmic ray flux with energies near 1 EeV and above are detectable with the Surface Detector array of the Pierre Auger Observatory. We report here on searches through Auger data from 1 January 2004 until 20 June 2013. No neutrino candidates were found, yielding a limit to the diffuse flux of ultra-high energy neutrinos that challenges the Waxman-Bahcall bound predictions. Neutrino identification is attempted using the broad time-structure of the signals expected in the SD stations, and is efficiently done for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for "Earth-skimming" neutrino interactions in the case of tau neutrinos. In this paper the searches for downward-going neutrinos in the zenith angle bins $60^\circ-75^\circ$ and $75^\circ-90^\circ$ as well as for upward-going neutrinos, are combined to give a single limit. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high energy neutrinos with an $E^{-2}$ spectrum in the energy range $1.0 \times 10^{17}$ eV - $2.5 \times 10^{19}$ eV is $E_\nu^2 dN_\nu/dE_\nu < 6.4 \times 10^{-9}~ {\rm GeV~ cm^{-2}~ s^{-1}~ sr^{-1}}$.

We have surveyed 84 Class 0, Class I, and flat-spectrum protostars in mid-infrared [Si II], [Fe II] and [S I] line emission, and 11 of these in far-infrared [O I] emission. We use the results to derive their mass outflow rates. Thereby we observe a strong correlation of mass outflow rates with bolometric luminosity, and with the inferred mass accretion rates of the central objects, which continues through the Class 0 range the trend observed in Class II young stellar objects. Along this trend from large to small mass-flow rates, the different classes of young stellar objects lie in the sequence Class 0 -- Class I/flat-spectrum -- Class II, indicating that the trend is an evolutionary sequence in which mass outflow and accretion rates decrease together with increasing age, while maintaining rough proportionality. The survey results include two which are key tests of magnetocentrifugal outflow-acceleration mechanisms: the distribution of the outflow/accretion branching ratio b, and limits on the distribution of outflow speeds. Neither rule out any of the three leading outflow-acceleration, angular-momentum-ejection mechanisms, but they provide some evidence that disk winds and accretion-powered stellar winds (APSWs) operate in many protostars. An upper edge observed in the branching-ratio distribution is consistent with the upper bound of b = 0.6 found in models of APSWs, and a large fraction (0.31) of the sample have branching ratio sufficiently small that only disk winds, launched on scales as large as several AU, have been demonstrated to account for them.

We have conducted three searches for correlations between ultra-high energy cosmic rays detected by the Telescope Array and the Pierre Auger Observatory, and high-energy neutrino candidate events from IceCube. Two cross-correlation analyses with UHECRs are done: one with 39 cascades from the IceCube `high-energy starting events' sample and the other with 16 high-energy `track events'. The angular separation between the arrival directions of neutrinos and UHECRs is scanned over. The same events are also used in a separate search using a maximum likelihood approach, after the neutrino arrival directions are stacked. To estimate the significance we assume UHECR magnetic deflections to be inversely proportional to their energy, with values $3^\circ$, $6^\circ$ and $9^\circ$ at 100 EeV to allow for the uncertainties on the magnetic field strength and UHECR charge. A similar analysis is performed on stacked UHECR arrival directions and the IceCube sample of through-going muon track events which were optimized for neutrino point-source searches.

Joint contributions of the Pierre Auger Collaboration and the Telescope Array Collaboration to the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands.

We present observations of HDCO and H2CO emission toward a sample of 15 Class 0 protostars in the Orion A and B clouds. Of these, eleven protostars are Herschel-identified PACS Bright Red Sources (PBRS) and four are previously identified protostars. Our observations revealed the chemical properties of the PBRS envelope for the first time. The column densities of HDCO and H2CO are derived from single dish observations at an angular resolution of ~20 arcsec (~8400 AU). The degree of deuteration in H2CO ([HDCO]/[H2CO]) was estimated to range from 0.03 to 0.31. The deuterium fractionation of most PBRS (70%) is similar to that of the non-PBRS sources. Three PBRS (30%) exhibit high deuterium fractionation, larger than 0.15. The large variation of the deuterium fractionation of H2CO in the whole PBRS sample may reflect the diversity in the initial conditions of star forming cores. There is no clear correlation between the [HDCO]/[H2CO] ratio and the evolutionary sequence of protostars.

Contributions of the Pierre Auger Collaboration to the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands

We present the first study of the relationship between the column density distribution of molecular clouds within nearby Galactic spiral arms and their evolutionary status as measured from their stellar content. We analyze a sample of 195 molecular clouds located at distances below 5.5 kpc, identified from the ATLASGAL 870 micron data. We define three evolutionary classes within this sample: starless clumps, star-forming clouds with associated young stellar objects, and clouds associated with HII regions. We find that the N(H2) probability density functions (N-PDFs) of these three classes of objects are clearly different: the N-PDFs of starless clumps are narrowest and close to log-normal in shape, while star-forming clouds and HII regions exhibit a power-law shape over a wide range of column densities and log-normal-like components only at low column densities. We use the N-PDFs to estimate the evolutionary time-scales of the three classes of objects based on a simple analytic model from literature. Finally, we show that the integral of the N-PDFs, the dense gas mass fraction, depends on the total mass of the regions as measured by ATLASGAL: more massive clouds contain greater relative amounts of dense gas across all evolutionary classes.

We compare the structure of star-forming molecular clouds in different regions of Orion A to determine how the column density probability distribution function (N-PDF) varies with environmental conditions such as the fraction of young protostars. A correlation between the N-PDF slope and Class 0 protostar fraction has been previously observed in a low-mass star-formation region (Perseus) by Sadavoy; here we test if a similar correlation is observed in a high-mass star-forming region. We use Herschel data to derive a column density map of Orion A. We use the Herschel Orion Protostar Survey catalog for accurate identification and classification of the Orion A young stellar object (YSO) content, including the short-lived Class 0 protostars (with a $\sim$ 0.14 Myr lifetime). We divide Orion A into eight independent 13.5 pc$^2$ regions; in each region we fit the N-PDF distribution with a power-law, and we measure the fraction of Class 0 protostars. We use a maximum likelihood method to measure the N-PDF power-law index without binning. We find that the Class 0 fraction is higher in regions with flatter column density distributions. We test the effects of incompleteness, YSO misclassification, resolution, and pixel-scale. We show that these effects cannot account for the observed trend. Our observations demonstrate an association between the slope of the power-law N-PDF and the Class 0 fractions within Orion A. Various interpretations are discussed including timescales based on the Class 0 protostar fraction assuming a constant star-formation rate. The observed relation suggests that the N-PDF can be related to an "evolutionary state" of the gas. If universal, such a relation permits an evaluation of the evolutionary state from the N-PDF power-law index at much greater distances than those accesible with protostar counts. (abridged)

We report the dramatic mid-infrared brightening between 2004 and 2006 of HOPS 383, a deeply embedded protostar adjacent to NGC 1977 in Orion. By 2008, the source became a factor of 35 brighter at 24 microns with a brightness increase also apparent at 4.5 microns. The outburst is also detected in the submillimeter by comparing APEX/SABOCA to SCUBA data, and a scattered-light nebula appeared in NEWFIRM K_s imaging. The post-outburst spectral energy distribution indicates a Class 0 source with a dense envelope and a luminosity between 6 and 14 L_sun. Post-outburst time-series mid- and far-infrared photometry shows no long-term fading and variability at the 18% level between 2009 and 2012. HOPS 383 is the first outbursting Class 0 object discovered, pointing to the importance of episodic accretion at early stages in the star formation process. Its dramatic rise and lack of fading over a six-year period hint that it may be similar to FU Ori outbursts, although the luminosity appears to be significantly smaller than the canonical luminosities of such objects.

We present CARMA 2.9 mm dust continuum emission observations of a sample of 14 Herschel-detected Class 0 protostars in the Orion A and B molecular clouds, drawn from the PACS Bright Red Sources (PBRS) sample (Stutz et al.). These objects are characterized by very red 24 \micron\ to 70 \micron\ colors and prominent submillimeter emission, suggesting that they are very young Class 0 protostars embedded in dense envelopes. We detect all of the PBRS in 2.9 mm continuum emission and emission from 4 protostars and 1 starless core in the fields toward the PBRS; we also report 1 new PBRS source. The ratio of 2.9 mm luminosity to bolometric luminosity is higher by a factor of $\sim$5 on average, compared to other well-studied protostars in the Perseus and Ophiuchus clouds. The 2.9 mm visibility amplitudes for 6 of the 14 PBRS are very flat as a function of uv-distance, with more than 50\% of the source emission arising from radii $<$ 1500 AU. These flat visibility amplitudes are most consistent with spherically symmetric envelope density profiles with $\rho$~$\propto$~R$^{-2.5}$. Alternatively, there could be a massive unresolved structure like a disk or a high-density inner envelope departing from a smooth power-law. The large amount of mass on scales $<$ 1500 AU (implying high average central densities) leads us to suggest that that the PBRS with flat visibility amplitude profiles are the youngest PBRS and may be undergoing a brief phase of high mass infall/accretion and are possibly among the youngest Class 0 protostars. The PBRS with more rapidly declining visibility amplitudes still have large envelope masses, but could be slightly more evolved.

Context: The initial conditions for the gravitational collapse of molecular cloud cores and the subsequent birth of stars are still not well constrained. The characteristic cold temperatures (about 10 K) in such regions require observations at sub-millimetre and longer wavelengths. The Herschel Space Observatory and complementary ground-based observations presented in this paper have the unprecedented potential to reveal the structure and kinematics of a prototypical core region at the onset of stellar birth. Aims: This paper aims to determine the density, temperature, and velocity structure of the star-forming Bok globule CB 17. This isolated region is known to host (at least) two sources at different evolutionary stages: a dense core, SMM1, and a Class I protostar, IRS. Methods: We modeled the cold dust emission maps from 100 micron to 1.2 mm with both a modified blackbody technique to determine the optical depth-weighted line-of-sight temperature and column density and a ray-tracing technique to determine the core temperature and volume density structure. Furthermore, we analysed the kinematics of CB17 using the high-density gas tracer N2H+. Results: From the ray-tracing analysis, we find a temperature in the centre of SMM1 of 10.6 K, a flat density profile with radius 9500 au, and a central volume density of n(H) = 2.3x10^5 cm-3. The velocity structure of the N2H+ observations reveal global rotation with a velocity gradient of 4.3 km/s/pc. Superposed on this rotation signature we find a more complex velocity field, which may be indicative of differential motions within the dense core. Conclusions: SMM is a core in an early evolutionary stage at the verge of being bound, but the question of whether it is a starless or a protostellar core remains unanswered.

Spherical harmonic moments are well-suited for capturing anisotropy at any scale in the flux of cosmic rays. An unambiguous measurement of the full set of spherical harmonic coefficients requires full-sky coverage. This can be achieved by combining data from observatories located in both the northern and southern hemispheres. To this end, a joint analysis using data recorded at the Telescope Array and the Pierre Auger Observatory above $10^{19}$ eV is presented in this work. The resulting multipolar expansion of the flux of cosmic rays allows us to perform a series of anisotropy searches, and in particular to report on the angular power spectrum of cosmic rays above $10^{19}$ eV. No significant deviation from isotropic expectations is found throughout the analyses performed. Upper limits on the amplitudes of the dipole and quadrupole moments are derived as a function of the direction in the sky, varying between 7% and 13% for the dipole and between 7% and 10% for a symmetric quadrupole.

Recent observational and theoretical investigations have emphasised the importance of filamentary networks within molecular clouds as sites of star formation. Since such environments are more complex than those of isolated cores, it is essential to understand how the observed line profiles from collapsing cores with non-spherical geometry are affected by filaments. In this study, we investigate line profile asymmetries by performing radiative transfer calculations on hydrodynamic models of three collapsing cores that are embedded in filaments. We compare the results to those that are expected for isolated cores. We model the five lowest rotational transition line (J = 1-0, 2-1, 3-2, 4-3, and 5-4) of both optically thick (HCN, HCO$^+$) as well as optically thin (N$_2$H$^+$, H$^{13}$CO$^+$) molecules using constant abundance laws. We find that less than 50% of simulated (1-0) transition lines show blue infall asymmetries due to obscuration by the surrounding filament. However, the fraction of collapsing cores that have a blue asymmetric emission line profile rises to 90% when observed in the (4-3) transition. Since the densest gas towards the collapsing core can excite higher rotational states, upper level transitions are more likely to produce blue asymmetric emission profiles. We conclude that even in irregular, embedded cores one can trace infalling gas motions with blue asymmetric line profiles of optically thick lines by observing higher transitions. The best tracer of collapse motions of our sample is the (4-3) transition of HCN, but the (3-2) and (5-4) transitions of both HCN and HCO$^+$ are also good tracers.