Massive stars produce so much light that the radiation pressure they exert on the gas and dust around them is stronger than their gravitational attraction, a condition that has long been expected to prevent them from growing by accretion. We present three-dimensional radiation-hydrodynamic simulations of the collapse of a massive prestellar core and find that radiation pressure does not halt accretion. Instead, gravitational and Rayleigh-Taylor instabilities channel gas onto the star system through nonaxisymmetric disks and filaments that self-shield against radiation while allowing radiation to escape through optically thin bubbles. Gravitational instabilities cause the disk to fragment and form a massive companion to the primary star. Radiation pressure does not limit stellar masses, but the instabilities that allow accretion to continue lead to small multiple systems.

Forming stars emit a substantial amount of radiation into their natal environment. We use ORION, an adaptive mesh refinement (AMR) three-dimensional gravito-radiation-hydrodyanics code, to simulate low-mass star formation in a turbulent molecular cloud. We compare the distributions of stellar masses, accretion rates, and temperatures in the cases with and without radiative transfer, and we demonstrate that radiative feedback has a profound effect on accretion, multiplicity, and mass by reducing the number of stars formed and the total rate at which gas turns into stars. We also show that once the star formation reaches a steady state, protostellar radiation is by far the dominant source of energy in the simulation, exceeding viscous dissipation and compressional heating by at least an order of magnitude. Calculations that omit radiative feedback from protstars significantly underestimate the gas temperature and the strength of this effect. Although heating from protostars is mainly confined to the protostellar cores, we find that it is sufficient to suppress disk fragmentation that would otherwise result in very low-mass companions or brown dwarfs. We demonstrate that the mean protostellar accretion rate increases with the final stellar mass so that the star formation time is only a weak function of mass.

We have conducted a survey of 328 protostars in the Orion molecular clouds with the Atacama Large Millimeter/ submillimeter Array at 0.87 mm at a resolution of ∼0 1 (40 au), including observations with the Very Large Array at 9mm toward 148 protostars at a resolution of ∼0 08 (32 au). This is the largest multiwavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 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 3.4 5.8 , -+ 37.0 3.0 4.9 , and -+ 28.5 2.3 3.7 au, respectively, and the mean protostellar dust disk masses are 25.9 -+ 4.0 7.7 , -+ 14.9 2.2 3.8 , -+11.6 1.93.5 Å M , 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% of 379 detected 0.87 mm continuum sources plus 42 nondetections) 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 those 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 characterize the infall rate onto protostellar systems forming in self-gravitating radiation-hydrodynamics simulations. Using two dimensionless parameters to determine the disks' susceptibility to gravitational fragmentation, we infer limits on protostellar system multiplicity and the mechanism of binary formation. We show that these parameters give robust predictions even in the case of marginally resolved protostellar disks. We find that protostellar systems with radiation feedback predominately form binaries via turbulent fragmentation, not disk instability, and predict that turbulent fragmentation is the dominant channel for binary formation for low-mass stars. We clearly demonstrate that systems forming in simulations including radiative feedback have fundamentally different parameters than those in purely hydrodynamics simulations.

We present an overview of the first data release (DR1) and first-look science from the Green Bank Ammonia Survey (GAS). GAS is a Large Program at the Green Bank Telescope to map all Gould Belt star-forming regions with A V 7 mag visible from the northern hemisphere in emission from NH 3 and other key molecular tracers. This first release includes the data for four regions in Gould Belt clouds: B18 in Taurus, NGC 1333 in Perseus, L1688 in Ophiuchus, and Orion A North in Orion. We compare the NH 3 emission to dust continuum emission from Herschel, and find that the two tracers correspond closely. NH 3 is present in over 60 % of lines-of-sight with A V 7 mag in three of the four DR1 regions, in agreement with expectations from previous observations. The sole exception is B18, where NH 3 is detected toward ∼ 40 % of lines-of-sight with A V 7 mag. Moreover, we find that the NH 3 emission is generally extended beyond the typical 0.1 pc length scales of dense cores. We produce maps of the gas kinematics, temperature, and NH 3 column densities through forward modeling of the hyperfine structure of the NH 3 (1,1) and (2,2) lines. We show that the NH 3 velocity dispersion, σ v , and gas kinetic temperature, T K , vary systematically between the regions included in this release, with an increase in both the mean value and spread of σ v and T K with increasing star formation activity. The data presented in this paper are publicly available.

The protostellar luminosity function (PLF) is the present-day luminosity function of the protostars in a region of star formation. It is determined using the protostellar mass function (PMF) in combination with a stellar evolutionary model that provides the luminosity as a function of instantaneous and final stellar mass. As in McKee & Offner (2010), we consider three main accretion models: the Isothermal Sphere model, the Turbulent Core model, and an approximation of the Competitive Accretion model. We also consider the effect of an accretion rate that tapers off linearly in time and an accelerating star formation rate. For each model, we characterize the luminosity distribution using the mean, median, maximum, ratio of the median to the mean, standard deviation of the logarithm of the luminosity, and the fraction of very low luminosity objects. We compare the models with bolometric luminosities observed in local star forming regions and find that models with an approximately constant accretion time, such as the Turbulent Core and Competitive Accretion models, appear to agree better with observation than those with a constant accretion rate, such as the Isothermal Sphere model. We show that observations of the mean protostellar luminosity in these nearby regions of low-mass star formation suggest a mean star formation time of 0.3±0.1 Myr. Such a timescale, together with some accretion that occurs non-radiatively and some that occurs in high-accretion, episodic bursts, resolves the classical "luminosity problem" in low-mass star formation, in which observed protostellar luminosities are significantly less than predicted. An accelerating star formation rate is one possible way of reconciling the observed star formation time and mean luminosity. Future observations will place tighter constraints on the observed luminosities, star formation time, and episodic accretion, enabling better discrimination between star formation models and clarifying the influence of variable accretion on the PLF.

We revisit the problem of low-mass pre-main-sequence stellar evolution and its observational consequences for where stars fall on the Hertzsprung-Russell diagram (HRD). In contrast to most previous work, our models follow stars as they grow from small masses via accretion, and we perform a systematic study of how the stars' HRD evolution is influenced by their initial radius, by the radiative properties of the accretion flow, and by the accretion history, using both simple idealized accretion histories and histories taken from numerical simulations of star cluster formation. We compare our numerical results to both non-accreting isochrones and to the positions of observed stars in the HRD, with a goal of determining whether both the absolute ages and the age dispersions inferred from non-accreting isochrones are reliable. We show that non-accreting isochrones can sometimes overestimate stellar ages for more massive stars (those with effective temperatures above ∼3500 K), thereby explaining why non-accreting isochrones often suggest a systematic age difference between more and less massive stars in the same cluster. However, we also find the only way to produce a similar overestimate for the ages of cooler stars is if these stars grow from ∼0.01 M seed protostars that are an order of magnitude smaller than predicted by current theoretical models, and if the size of the seed protostar correlates systematically with the final stellar mass at the end of accretion. We therefore conclude that, unless both of these conditions are met, inferred ages and age spreads for cool stars are reliable, at least to the extent that the observed bolometric luminosities and temperatures are accurate. Finally, we note that the time dependence of the mass accretion rate has remarkably little effect on low-mass stars' evolution on the HRD, and that such time dependence may be neglected for all stars except those with effective temperatures above ∼4000 K.

We use gas temperature and velocity dispersion data from the Green Bank Ammonia Survey and core masses and sizes from the James Clerk Maxwell Telescope Gould Belt Survey to estimate the virial states of dense cores within the Orion A molecular cloud. Surprisingly, we find that almost none of the dense cores are sufficiently massive to be bound when considering only the balance between selfgravity and the thermal and non-thermal motions present in the dense gas. Including the additional pressure binding imposed by the weight of the ambient molecular cloud material and additional smaller pressure terms, however, suggests that most of the dense cores are pressure confined.