We report results of 0.05"-resolution observations toward the O-type proto-binary system IRAS 16547-4247 with the Atacama Large Millimeter/submillimeter Array (ALMA). We present dynamical and chemical structures of the circumbinary disk, circumstellar disks, outflows and jets, illustrated by multi-wavelength continuum and various molecular lines. In particular, we detect sodium chloride, silicon compounds, and vibrationally-excited water lines as probes of the individual protostellar disks at a scale of 100 au. These are complementary to typical hot-core molecules tracing the circumbinary structures on a 1000-au scale. The H2O line tracing inner-disks has an upper-state energy of Eu/k>3000K, indicating a high temperature of the disks. On the other hand, despite the detected transitions of NaCl, SiO, and SiS not necessarily having high upper-state energies, they are enhanced only in the vicinity of the protostars. We interpret that these molecules are the products of dust destruction, which only happens in the inner disks. This is the second detection of alkali metal halide in protostellar systems after the case of the disk of Orion Source I, and also one of few massive protostellar disks associated with high-energy transition water and silicon compounds. These new results suggest these "hot-disk" lines may be common in innermost disks around massive protostars, and have great potential for future research of massive star formation. We also tentatively find that the twin disks are counter-rotating, which might give a hint of the origin of the massive proto-binary system IRAS 16547-4247.

We introduce new analysis methods for studying the star cluster formation processes in Orion A, especially examining the scenario of a cloud-cloud collision. We utilize the CARMA-NRO Orion survey $^{13}$CO (1-0) data to compare molecular gas to the properties of YSOs from the SDSS III IN-SYNC survey. We show that the increase of $v_{\rm 13CO} - v_{\rm YSO}$ and $\Sigma$ scatter of older YSOs can be signals of cloud-cloud collision. SOFIA-upGREAT 158$\mu$m [CII] archival data toward the northern part of Orion A are also compared to the $^{13}$CO data to test whether the position and velocity offsets between the emission from these two transitions resemble those predicted by a cloud-cloud collision model. We find that the northern part of Orion A, including regions ONC-OMC-1, OMC-2, OMC-3 and OMC-4, shows qualitative agreements with the cloud-cloud collision scenario, while in one of the southern regions, NGC1999, there is no indication of such a process in causing the birth of new stars. On the other hand, another southern cluster, L1641N, shows slight tendencies of cloud-cloud collision. Overall, our results support the cloud-cloud collision process as being an important mechanism for star cluster formation in Orion A.

We present multi-wavelength images observed with SOFIA-FORCAST from ~10 to 40 $\mu$m of 14 protostars, selected as intermediate-mass protostar candidates, as part of the SOFIA Massive (SOMA) Star Formation Survey. We build protostellar spectral energy distributions (SEDs) with the SOFIA observations, together with archival data from Spitzer, Herschel and IRAS. We then fit the SEDs with radiative transfer (RT) models of Zhang & Tan (2018), based on Turbulent Core Accretion theory, to estimate key properties of the protostars. The SEDs generally indicate the validity of these RT models down to intermediate-mass and/or early-stage protostars. The protostars analyzed so far in the SOMA survey span a range of luminosities from ~$10^{2}$ to ~$10^{6} L_{\odot}$, current protostellar masses from ~0.5 to ~30 $M_{\odot}$, and ambient clump mass surface densities, $\Sigma_{\rm cl}$ of 0.1 - 3 g cm$^{-2}$. A wide range of evolutionary states of the individual protostars and of the protocluster environments are also probed. The 19 to 37 $\mu$m spectral index of the sources correlates with outflow cavity opening angle, ratio of this angle to viewing angle, and evolutionary stage. We have also added a sample of ~50 protostellar sources identified from within Infrared Dark Clouds and expected to be at the earliest stages of their evolution. With this global sample, most of the evolutionary stages of high- and intermediate-mass protostars are probed. From the best fitting models of the protostars, there is no evidence of a threshold value of protocluster clump mass surface density being needed to form protostars up to ~25 $M_\odot$. However, to form more massive protostars, there is tentative evidence that $\Sigma_{\rm cl}$ needs to be at least 1 g cm$^{-2}$. We discuss how this is consistent with expectations from core accretion models that include internal feedback from the forming massive star.

We conduct a census of the high-mass protostellar population of the $\sim70,000\:M_\odot$ Infrared Dark Cloud (IRDC) G028.37+00.07, identifying 35 sources based on their $70\:\mu$m emission, as reported in the {\it Herschel} Hi-GAL catalog of Molinari et al. (2016). We perform aperture photometry to construct spectral energy distributions (SEDs), which are then fit with the massive protostar models of Zhang & Tan (2018). We find that the sources span a range of isotropic luminosities from $\sim$20 to 4,500$\:L_\odot$. The most luminous sources are predicted to have current protostellar masses of $m_{*}\sim10\:M_\odot$ forming from cores of mass $M_{c}\sim40$ to $400\:M_\odot$. The least luminous sources in our sample are predicted to be protostars with masses as low as $\sim 0.5\:M_\odot$ forming from cores with $M_{c}\sim10\:M_\odot$, which are the minimum values explored in the protostellar model grid. The detected protostellar population has a total estimated protostellar mass of $M_{*}\sim 100\:M_\odot$. Allowing for completeness corrections, which are constrained by comparison with an ALMA study in part of the cloud, we estimate a star formation efficiency per free-fall time of $\sim3\%$ in the IRDC. Finally, analyzing the spatial distribution of the sources, we find relatively low degrees of central concentration of the protostars. The protostars, including the most massive ones, do not appear to be especially centrally concentrated in the protocluster as defined by the IRDC boundary.

We present a near-infrared (NIR) variability analysis for an 6\arcmin $\times$ 6\arcmin region, which encompasses the massive protocluster G286.21+0.17. The total sample comprises more than 5000 objects, of which 562 show signs of a circumstellar disk based on their infrared colors. The data includes HST observations taken in two epochs separated by 3 years in the F110W and F160W bands. 363 objects (7% of the sample) exhibit NIR variability at a significant level (Stetson index >1.7), and a higher variability fraction (14%) is found for the young stellar objects (YSOs) with disk excesses. We identified 4 high amplitude (>0.6 mag) variables seen in both NIR bands. Follow up and archival observations of the most variable object in this survey (G286.2032+0.1740) reveal a rising light curve over 8 years from 2011 to 2019, with a K band brightening of 3.5 mag. Overall the temporal behavior of G286.2032+0.1740 resembles that of typical FU Ori objects, however its pre-burst luminosity indicates it has a very low mass ($<0.12\:M_\odot$), making it an extreme case of an outburst event that is still ongoing.

We study the gas kinematics and dynamics of the massive protocluster G286.21+0.17 with the Atacama Large Millimeter/submillimeter Array using spectral lines of $C^{18}O$(2-1), $N_2D^+$(3-2), $DCO^+$(3-2) and $DCN$(3-2). On the parsec clump scale, $C^{18}O$ emission appears highly filamentary around the systemic velocity. $N_2D^+$ and $DCO^+$ are more closely associated with the dust continuum. $DCN$ is strongly concentrated towards the protocluster center, where no or only weak detection is seen for $N_2D^+$ and $DCO^+$, possibly due to this region being at a relatively evolved evolutionary stage. Spectra of 76 continuum defined dense cores, typically a few 1000 AU in size, are analysed to measure their centroid velocities and internal velocity dispersions. There are no statistically significant velocity offsets of the cores among the different dense gas tracers. Furthermore, the majority (71\%) of the dense cores have subthermal velocity offsets with respect to their surrounding, lower density $C^{18}O$ emitting gas. Within the uncertainties, the dense cores in G286 show internal kinematics that are consistent with being in virial equilibrium. On clumps scales, the core to core velocity dispersion is also similar to that required for virial equilibrium in the protocluster potential. However, the distribution in velocity of the cores is largely composed of two spatially distinct groups, which indicates that the dense molecular gas has not yet relaxed to virial equilibrium, perhaps due to there being recent/continuous infall into the system.

We use Gaia DR2 to hunt for runaway and walkaway stars from the Orion Nebula Cluster (ONC). We search a large area of the sky extending 45$^\circ$ in radius around the ONC and out to a distance of 1 kpc to find sources that overlapped in angular position with the cluster in the last ~10 Myr, also accounting for the ONC's own proper motion. We construct a main sample of ~17,000 runaway/walkaway candidates selected by this 2D traceback condition. However, most of these are expected to be contaminants, e.g., caused by Galactic streaming motions of stars at different distances. We thus examine six further tests or "flags" to help identify real runaways, namely: (1) possessing young stellar object (YSO) colors and magnitudes based on Gaia optical photometry; (2) having IR excess consistent with YSOs based on 2MASS and WISE photometry; (3) having a high degree of optical variability; (4) having closest approach distances to the ONC core well constrained within the half-mass radius of the cluster; (5) having direction of ejection from the ONC that avoids the main contamination zone due to Galactic streaming; and (6) having a required radial velocity (RV) for 3D overlap with the ONC that is of reasonable magnitude (or, for the small fraction, 7%, of candidate stars with Gaia measured RVs, satisfying 3D traceback). Thirteen sources, not previously noted as Orion members, pass all of these tests, while another twelve are similarly promising, except that they are located in the zone most contaminated by Galactic streaming. Among these 25 ejection candidates, there are ten with measured RVs that pass the most restrictive 3D traceback condition. We present full lists of runaway/walkaway candidate stars, estimate the high-velocity population ejected from the ONC and discuss its implications for star cluster formation theories via comparison with numerical simulations.

In order to obtain a complete census of the stellar and sub-stellar population, down to a few M$_{Jup}$ in the $\sim1$ Myr old Orion Nebula Cluster, we used the infrared channel of the Wide Field Camera 3 of the Hubble Space Telescope with the F139M and F130N filters. These bandpasses correspond to the $1.4 \mu$m H$_2$O absorption feature and an adjacent line-free continuum region. Out of $4,504$ detected sources, $3,352$ (about $75\%$) appear fainter than m$_{130}=14$ (Vega mag) in the F130N filter, a brightness corresponding to the hydrogen-burning limit mass (M$\simeq 0.072 M_\odot$) at $\sim 1$ Myr. Of these, however, only $742$ sources have a negative F130M-139N color index, indicative of the presence of H$_2$O vapor in absorption, and can therefore be classified as bona-fide M and L dwarfs, with effective temperatures T$\lesssim 2850$ K at an assumed $1$ Myr cluster age. On our color-magnitude diagram, this population of sources with H$_2$O absorption appears clearly distinct from the larger background population of highly reddened stars and galaxies with positive F130M-F139N color index, and can be traced down to the sensitivity limit of our survey, m$_{130}\simeq 21.5$, corresponding to a $1$ Myr old $\simeq 3 $M$_{Jup}$, planetary mass object under about 2 magnitudes of visual extinction. Theoretical models of the BT-Settl family predicting substellar isochrones of $1, 2$ and $3$ Myr (down to $\sim 1 $M$_{Jup}$) fail to reproduce the observed H$_2$O color index at M$\lesssim 20 $M$_{Jup}$. We perform a Bayesian analysis to determine extinction, mass and effective temperature of each sub-stellar member of our sample, together with its membership probability.

We investigate the formation of stars within giant molecular clouds (GMCs) evolving in environments of different global magnetic field strength and large-scale dynamics. Building upon a series of magnetohydrodynamic (MHD) simulations of non-colliding and colliding GMCs, we employ density- and magnetically-regulated star formation sub-grid models in clouds which range from moderately magnetically supercritical to near critical. We examine gas and star cluster morphologies, magnetic field strengths and relative orientations, pre-stellar core densities, temperatures, mass-to-flux ratios and velocities, star formation rates and efficiencies over time, spatial clustering of stars, and kinematics of the stars and natal gas. The large scale magnetic criticality of the region greatly affects the overall gas evolution and star formation properties. GMC collisions enhance star formation rates and efficiencies in magnetically supercritical conditions, but may actually inhibit them in the magnetically critical case. This may have implications for star formation in different Galactic environments such as the Galactic Center and the main Galactic disk.

The bubble size distribution of ionized hydrogen regions, which could be derived from the tomographic imaging data of the redshifted 21~cm signal from the epoch of reionization, probes the information about the morphology of \HII\ bubbles during the reionization. However, 21~cm imaging is observationally very challenging even for the upcoming large radio interferometers. Given that these interferometers promise to measure the 21~cm power spectrum accurately, we propose a new method, which is based on the artificial neural networks (ANN), to reconstruct the \HII\ bubble size distribution from the 21~cm power spectrum. We demonstrate that the reconstruction from the 21~cm power spectrum can be almost as accurate as directly measured from the imaging data with the fractional error $\lesssim 10\%$, even with thermal noise at the sensitivity level of the Square Kilometre Array. Nevertheless, systematic errors might arise from approximations made in reionization simulations used for training the ANN. This paper, as the first in a series, exemplifies the possibility of recovering from the 21~cm power spectrum with ANN additional statistics of cosmic reionization that could not otherwise be inferred from the power spectrum analysis directly in the conventional methods.

Massive protostars generate strong radiation feedback, which may help set the mass they achieve by the end of the accretion process. Studying such feedback is therefore crucial for understanding the formation of massive stars. We report the discovery of a photoionized bipolar outflow towards the massive protostar G45.47+0.05 using high-resolution observations at 1.3 mm with the Atacama Large Millimeter/Submillimeter Array (ALMA) and at 7 mm with the Karl G. Jansky Very Large Array (VLA). By modeling the free-free continuum, the ionized outflow is found to be a photoevaporation flow with an electron temperature of 10,000 K and an electron number density of ~1.5x10^7 cm^-3 at the center, launched from a disk of radius of 110 au. H30alpha hydrogen recombination line emission shows strong maser amplification, with G45 being one of very few sources to show such millimeter recombination line masers. The mass of the driving source is estimated to be 30-50 Msun based on the derived ionizing photon rate, or 30-40 Msun based on the H30alpha kinematics. The kinematics of the photoevaporated material is dominated by rotation close to the disk plane, while accelerated to outflowing motion above the disk plane. The mass loss rate of the photoevaporation outflow is estimated to be ~(2-3.5)x10^-5 Msun/yr. We also found hints of a possible jet embedded inside the wide-angle ionized outflow with non-thermal emissions. The possible co-existence of a jet and a massive photoevaporation outflow suggests that, in spite of the strong photoionization feedback, accretion is still on-going.

Magneto-hydrodynamic (MHD) shocks are violent events that inject large amounts of energy in the interstellar medium (ISM) dramatically modifying its physical properties and chemical composition. Indirect evidence for the presence of such shocks has been reported from the especial chemistry detected toward a variety of astrophysical shocked environments. However, the internal physical structure of these shocks remains unresolved since their expected spatial scales are too small to be measured with current instrumentation. Here we report the first detection of a fully spatially resolved, MHD shock toward the Infrared Dark Cloud (IRDC) G034.77-00.55. The shock, probed by Silicon Monoxide (SiO) and observed with the Atacama Large Millimetre/sub-millimetre Array (ALMA), is associated with the collision between the dense molecular gas of the cloud and a molecular gas flow pushed toward the IRDC by the nearby supernova remnant (SNR) W44. The interaction is occurring on sub-parsec spatial scales thanks to the enhanced magnetic field of the SNR, making the dissipation region of the MHD shock large enough to be resolved with ALMA. Our observations suggest that molecular flow-flow collisions can be triggered by stellar feedback, inducing shocked molecular gas densities compatible with those required for massive star formation.

It is important to determine if massive stars form via disc accretion, like their low-mass counterparts. Theory and observation indicate that protostellar jets are a natural consequence of accretion discs and are likely to be crucial for removing angular momentum during the collapse. However, massive protostars are typically rarer, more distant and more dust enshrouded, making observational studies of their jets more challenging. A fundamental question is whether the degree of ionisation in jets is similar across the mass spectrum. Here we determine an ionisation fraction of $\sim5-12\%$ in the jet from the massive protostar G35.20-0.74N, based on spatially coincident infrared and radio emission. This is similar to the values found in jets from lower-mass young stars, implying a unified mechanism of shock ionisation applies in jets across most of the protostellar mass spectrum, up to at least $\sim10$ solar masses.

We perform a sequence of 3D magnetohydrodynamic (MHD) simulations of the outflow-core interaction for a massive protostar forming via collapse of an initial cloud core of $60~{M_\odot}$. This allows us to characterize the properties of disk wind driven outflows from massive protostars, which can allow testing of different massive star formation theories. It also enables us to assess quantitatively the impact of outflow feedback on protostellar core morphology and overall star formation efficiency. We find that the opening angle of the flow increases with increasing protostellar mass, in agreement with a simple semi-analytic model. Once the protostar reaches $\sim24~{M_\odot}$ the outflow's opening angle is so wide that it has blown away most of the envelope, thereby nearly ending its own accretion. We thus find an overall star formation efficiency of $\sim50\%$, similar to that expected from low-mass protostellar cores. Our simulation results therefore indicate that the MHD disk wind outflow is the dominant feedback mechanism for helping to shape the stellar initial mass function from a given prestellar core mass function.

We present ALMA CO(2-1) observations toward a massive infrared dark cloud G28.37+0.07. The ALMA data reveal numerous molecular (CO) outflows with a wide range of sizes throughout the cloud. Sixty-two 1.3 mm continuum cores were identified to be driving molecular outflows. We have determined the position angle in the plane-of-sky of 120 CO outflow lobes and studied their distribution. We find that the distribution of the plane-of-sky outflow position angles peaks at about 100 degree, corresponding to a concentration of outflows with an approximately east-west direction. For most outflows, we have been able to estimate the plane-of-sky angle between the outflow axis and the filament that harbors the protostar that powers the outflow. Statistical tests strongly indicate that the distribution of outflow-filament orientations is consistent with most outflow axes being mostly orthogonal to their parent filament in 3D. Such alignment may result from filament fragmentation or continuous mass transportation from filament to the embedded protostellar core. The latter is suggested by recent numerical studies with moderately strong magnetic fields.

Almost all massive stars have bound stellar companions, existing in binaries or higher-order multiples. While binarity is theorized to be an essential feature of how massive stars form, essentially all information about such properties is derived from observations of already formed stars, whose orbital properties may have evolved since birth. Little is known about binarity during formation stages. Here we report high angular resolution observations of 1.3 mm continuum and H30alpha recombination line emission, which reveal a massive protobinary with apparent separation of 180 au at the center of the massive star-forming region IRAS07299-1651. From the line-of-sight velocity difference of 9.5 km/s of the two protostars, the binary is estimated to have a minimum total mass of 18 solar masses, consistent with several other metrics, and maximum period of 570 years, assuming a circular orbit. The H30alpha line from the primary protostar shows kinematics consistent with rotation along a ring of radius of 12 au. The observations indicate that disk fragmentation at several hundred au may have formed the binary, and much smaller disks are feeding the individual protostars.

We present multi-wavelength images observed with SOFIA-FORCAST from $\sim$10 to 40 $\mu$m of seven high luminosity massive protostars, as part of the SOFIA Massive (SOMA) Star Formation Survey. Source morphologies at these wavelengths appear to be influenced by outflow cavities and extinction from dense gas surrounding the protostars. Using these images, we build spectral energy distributions (SEDs) of the protostars, also including archival data from Spitzer, Herschel and other facilities. Radiative transfer (RT) models of Zhang & Tan (2018), based on Turbulent Core Accretion theory, are then fit to the SEDs to estimate key properties of the protostars. Considering the best five models fit to each source, the protostars have masses $m_{*} \sim 12-64 \: M_{\odot}$ accreting at rates of $\dot{m}_{*} \sim 10^{-4}-10^{-3} \: M_{\odot} \: \rm yr^{-1}$ inside cores of initial masses $M_{c} \sim 100-500 \: M_{\odot}$ embedded in clumps with mass surface densities $\Sigma_{\rm cl} \sim 0.1-3 \: \rm g \: cm^{-2}$ and span a luminosity range of $10^{4} -10^{6} \: L_{\odot}$. Compared with the first eight protostars in Paper I, the sources analyzed here are more luminous, and thus likely to be more massive protostars. They are often in a clustered environment or have a companion protostar relatively nearby. From the range of parameter space of the models, we do not see any evidence that $\Sigma_{\rm cl}$ needs to be high to form these massive stars. For most sources the RT models provide reasonable fits to the SEDs, though the cold clump material often influences the long wavelength fitting. However, for sources in very clustered environments, the model SEDs may not be such a good description of the data, indicating potential limitations of the models for these regions.

We report molecular line observations of the massive protostellar source G339.88-1.26 with the Atacama Large Millimeter/Submillimeter Array. The observations reveal a highly collimated SiO jet extending from the 1.3 mm continuum source, which connects to a slightly wider but still highly collimated CO outflow. Rotational features perpendicular to the outflow axis are detected in various molecular emissions, including SiO, SO2, H2S, CH3OH, and H2CO emissions. Based on their spatial distributions and kinematics, we find that they trace different parts of the envelope-disk system. The SiO emission traces the disk and inner envelope in addition to the jet. The CH3OH and H2CO emissions mostly trace the infalling-rotating envelope, and are enhanced around the transition region between envelope and disk, i.e., the centrifugal barrier. The SO2 and H2S emissions are enhanced around the centrifugal barrier, and also trace the outer part of the disk. Envelope kinematics are consistent with rotating-infalling motion, while those of the disk are consistent with Keplerian rotation. The radius and velocity of the centrifugal barrier are estimated to be about 530 au and 6 km/s, leading to a central mass of about 11 solar masses, consistent with estimates based on spectral energy distribution fitting. These results indicate that an ordered transition from an infalling-rotating envelope to a Keplerian disk through a centrifugal barrier, accompanied by changes of types of molecular line emissions, is a valid description of this massive protostellar source. This implies that at least some massive stars form in a similar way as low-mass stars via Core Accretion.

We study centimeter continuum emission of eight high- and intermediate-mass protostars that are part of the SOFIA Massive (SOMA) Star Formation Survey, thus building extended spectral energy distributions (SEDs) from the radio to the infrared. We discuss the morphology seen in the centimeter continuum images, which are mostly derived from archival VLA data, and the relation to infrared morphology. We use the SEDs to test new models of high-mass star formation including radiative and disk-wind feedback and associated free-free and dust continuum emission (Tanaka, Tan, & Zhang 2016). We show that interferometric data of the centimeter continuum flux densities provide additional, stringent tests of the models by constraining the ionizing luminosity of the source and help to break degeneracies encountered when modeling the infrared-only SEDs, especially for the protostellar mass. Our derived parameters are consistent with physical parameters estimated by other methods such as dynamical protostellar masses. We find a few examples of additional stellar sources in the vicinity of the high-mass protostars, which may be low-mass young stellar objects. However, the stellar multiplicity of the regions, at least as traced by radio continuum emission, appears to be relatively low.

To probe the initial conditions for high-mass star and cluster formation, we investigate the properties of dense filaments within the infrared dark cloud G035.39-00.33 (IRDC G035.39) in a combined Very Large Array (VLA) and the Green Bank Telescope (GBT) mosaic tracing the NH3 (1,1) and (2,2) emission down to 0.08 pc scales. Using agglomerative hierarchical clustering on multiple line-of-sight velocity component fitting results, we identify seven extended velocity-coherent components in our data, likely representing spatially coherent physical structures, some exhibiting complex gas motions. The velocity gradient magnitude distribution peaks at its mode of 0.35 km/s/pc and has a long tail extending into higher values of 1.5 - 2 km/s/pc, and is generally consistent with those found toward the same cloud in other molecular tracers and with the values found towards nearby low-mass dense cloud cores at the same scales. Contrary to observational and theoretical expectations, we find the non-thermal ammonia line widths to be systematically narrower (by about 20%) than those of N2H+ (1-0) line transition observed with similar resolution. If the observed ordered velocity gradients represent the core envelope solid-body rotation, we estimate the specific angular momentum to be about 2 x 10^21 cm^2/s, similar to the low-mass star-forming cores. Together with the previous finding of subsonic motions in G035.39, our results demonstrate high levels of similarity between kinematics of a high-mass star-forming IRDC and the low-mass star formation regime.

We calculate cosmic distributions in space and time of the formation sites of the first, "Pop III.1" stars, exploring a model in which these are the progenitors of all supermassive black holes (SMBHs), seen in the centers of most large galaxies. Pop III.1 stars are defined to form from primordial composition gas in dark matter minihalos with $\sim10^6\:M_\odot$ that are isolated from neighboring astrophysical sources by a given isolation distance, $d_{\rm{iso}}$. We assume Pop III.1 sources are seeds of SMBHs, based on protostellar support by dark matter annihilation heating that allows them to accrete a large fraction of their minihalo gas, i.e., $\sim10^5\:M_\odot$. Exploring $d_{\rm{iso}}$ from $10 - 100\:\rm{kpc}$ (proper distances), we predict the redshift evolution of Pop III.1 source and SMBH remnant number densities. The local, $z=0$ density of SMBHs constrains $d_{\rm{iso}}\lesssim 100\:\rm{kpc}$ (i.e., $3\:\rm{Mpc}$ comoving distance at $z\simeq30$). In our simulated ($\sim60\:\rm{Mpc}$)$^3$ comoving volume, Pop III.1 stars start forming just after $z=40$. Their formation is largely complete by $z\simeq25$ to $20$ for $d_{\rm{iso}}=100$ to $50\:\rm{kpc}$. We follow source evolution to $z=10$, by which point most SMBHs reside in halos with $\gtrsim10^8\:M_\odot$. Over this period, there is relatively limited merging of SMBHs for these values of $d_{\rm{iso}}$. We also predict SMBH clustering properties at $z=10$: feedback suppression of neighboring sources leads to relatively flat angular correlation functions.

In this paper, we address two issues related to primordial disk evolution in three clusters (NGC 1333, IC 348, and Orion A) observed by the INfrared Spectra of Young Nebulous Clusters (IN-SYNC) project. First, in each cluster, averaged over the spread of age, we investigate how disk lifetime is dependent on stellar mass. The general relation in IC 348 and Orion A is that primordial disks around intermediate mass stars (2--5$M_{\odot}$) evolve faster than those around loss mass stars (0.1--1$M_{\odot}$), which is consistent with previous results. However, considering only low mass stars, we do not find a significant dependence of disk frequency on stellar mass. These results can help to better constrain theories on gas giant planet formation timescales. Secondly, in the Orion A molecular cloud, in the mass range of 0.35--0.7$M_{\odot}$, we provide the most robust evidence to date for disk evolution within a single cluster exhibiting modest age spread. By using surface gravity as an age indicator and employing 4.5 $\mu m$ excess as a primordial disk diagnostic, we observe a trend of decreasing disk frequency for older stars. The detection of intra-cluster disk evolution in NGC 1333 and IC 348 is tentative, since the slight decrease of disk frequency for older stars is a less than 1-$\sigma$ effect.

Orion BN has a large proper motion and radial velocity with respect to the gas and other stars in the region where it is presumed to have formed. Multiple dynamical interaction scenarios have been proposed to explain this motion. In one case BN is thought to have interacted with stars in the Trapezium cluster, while in another it is thought to have interacted with source I while deeply embedded in molecular gas. If there is dense gas that has been retained in close proximity to BN, it may be evidence that the latter scenario is favored. We observed BN at high spectral resolution in three windows near 6 $\mu$m using SOFIA/EXES targeting the $\nu_2$ vibrational band of H$_2$O. Absorption from only three transitions of H$_2$O is detected, and through kinematic analysis is associated with cool, dense foreground gas, not BN itself. We find no evidence for H$_2$O absorption or emission at the systemic velocity of BN.

We investigate the dynamical evolution of star clusters during their formation, assuming that they are born from a turbulent starless clump of a given mass that is embedded within a parent self-gravitating molecular cloud characterized by a particular mass surface density. In contrast to the standard practice of most $N$-body studies, we do not assume that all stars are formed at once. Rather, we explore the effects of different star formation rates on the global structure and evolution of young embedded star clusters, also considering various primordial binary fractions and mass segregation levels. Our fiducial clumps studied in this paper have initial masses of $M_{\rm cl} = 3000\,M_\odot$, are embedded in ambient cloud environments of $\Sigma_{\rm cloud} = 0.1$ and 1 g cm$^{-2}$, and gradually form stars with an overall efficiency of 50% until the gas is exhausted. We investigate star formation efficiencies per free-fall time in the range $\epsilon_{\rm ff}=0.01$ to 1, and also compare to the instantaneous case ($\epsilon_{\rm ff}=\infty$) of Paper I. We show that most of the interesting dynamical processes that determine the future of the cluster, happen during the early formation phase. In particular, the ejected stellar population is sensitive to the duration of star cluster formation: for example, clusters with longer formation times produce more runaway stars, since these clusters remain in a dense state for longer, thus favouring occurrence of dynamical ejections. We also show that the presence of radial age gradients in star clusters depends sensitively on the star formation efficiency per free fall time, with observed values being matched best by our slowest forming clusters with $\epsilon_{\rm ff}\lesssim0.03$.

We present high resolution (0.2", 1000 AU) 1.3 mm ALMA observations of massive infrared dark cloud clump, G028.37+00.07-C1, thought to harbor the early stages of massive star formation. Using $\rm N_2D^+$(3-2) we resolve the previously identified C1-S core, separating the bulk of its emission from two nearby protostellar sources. C1-S is thus identified as a massive ($\sim50\:M_\odot$), compact ($\sim0.1\:$pc diameter) starless core, e.g., with no signs of outflow activity. Being highly deuterated, this is a promising candidate for a pre-stellar core on the verge of collapse. An analysis of its dynamical state indicates a sub-virial velocity dispersion compared to a trans-Alfv\'enic turbulent core model. However, virial equilibrium could be achieved with sub-Alfv\'enic conditions involving $\sim2\:$mG magnetic field strengths.

The large population of Earth to super-Earth sized planets found very close to their host stars has motivated consideration of $in$ $situ$ formation models. In particular, Inside-Out Planet Formation is a scenario in which planets coalesce sequentially in the disk, at the local gas pressure maximum near the inner boundary of the dead zone. The pressure maximum arises from a decline in viscosity, going from the active innermost disk (where thermal ionization of alkalis yields high viscosities via the magneto-rotational instability (MRI)) to the adjacent dead zone (where the MRI is quenched). Previous studies of the pressure maximum, based on $\alpha$-disk models, have assumed ad hoc values for the viscosity parameter $\alpha$ in the active zone, ignoring the detailed physics of the MRI. Here we explicitly couple the MRI criteria to the $\alpha$-disk equations, to find steady-state (constant accretion rate) solutions for the disk structure. We consider the effects of both Ohmic and ambipolar resistivities, and find solutions for a range of disk accretion rates ($\dot{M}$ = $10^{-10}$ - $10^{-8}$ ${\rm M}_{\odot}$/yr), stellar masses ($M_{\ast}$ = 0.1 - 1 ${\rm M}_{\odot}$), and fiducial values of the $non$-MRI $\alpha$-viscosity in the dead zone ($\alpha_{\rm {DZ}} = 10^{-5}$ - $10^{-3}$). We find that: (1) A midplane pressure maximum forms radially $outside$ the inner boundary of the dead zone; (2) Hall resistivity dominates near the midplane in the inner disk, which may explain why close-in planets do $not$ form in $\sim$50% of systems; (3) X-ray ionization can be competitive with thermal ionization in the inner disk, because of the low surface density there in steady-state; and (4) our inner disk solutions are viscously unstable to surface density perturbations.

We study the core mass function (CMF) within 32 dense clumps in seven infrared dark clouds (IRDCs) with the Atacama Large Millimeter/submillimeter Array (ALMA) via 1.3~mm continuum emission at a resolution of $\sim$1". We have identified 107 cores with the dendrogram algorithm, with a median radius of about 0.02 pc. Their masses range from 0.261 to 178 $M_{\odot}$. After applying completeness corrections, we fit the combined IRDC CMF with a power law of the form $d N / d\:{\rm log} M \propto M^{-\alpha}$ and derive an index of $\alpha\simeq0.86\pm0.11$ for $M \geq 0.79\:M_\odot$ and $\alpha\simeq0.70\pm0.13$ for $M\geq 1.26\:M_\odot$, which is a significantly more top-heavy distribution than the Salpeter stellar initial mass function (IMF) index of 1.35. We also make a direct comparison of these IRDC clump CMF results to those measured in the more evolved protocluster G286 derived with similar methods, which have $\alpha\simeq1.29\pm0.19$ and $1.08\pm0.27$ in these mass ranges, respectively. These results provide a hint that, especially for the $M\geq 1.26\:M_\odot$ range where completeness corrections are modest, the CMF in high pressure, early-stage environments of IRDC clumps may be top-heavy compared to that in the more evolved, global environment of the G286 protoclusters. However, larger samples of cores probing these different environments are needed to better establish the robustness of this potential CMF variation.

We conduct a theoretical study of the formation of massive stars over a wide range of metallicities from 1e-5 to 1Zsun and evaluate the star formation efficiencies (SFEs) from prestellar cloud cores taking into account multiple feedback processes. Unlike for simple spherical accretion, in the case of disk accretion feedback processes do not set upper limits on stellar masses. At solar metallicity, launching of magneto-centrifugally-driven outflows is the dominant feedback process to set SFEs, while radiation pressure, which has been regarded to be pivotal, has only minor contribution even in the formation of over-100Msun stars. Photoevaporation becomes significant in over-20Msun star formation at low metallicities of <1e-2Zsun, where dust absorption of ionizing photons is inefficient. We conclude that if initial prestellar core properties are similar, then massive stars are rarer in extremely metal-poor environments of 1e-5 - 1e-3Zsun. Our results give new insight into the high-mass end of the initial mass function and its potential variation with galactic and cosmological environments.

The formation of the most massive O-type stars is poorly understood. We present a case study of a young massive clump from the ATLASGAL survey, G328.2551-0.5321. It exhibits a bolometric luminosity of 1.3$\times$10$^4$ L$_{\odot}$ corresponding to a current protostellar mass of $\sim$11 and 16 M$_{\odot}$. We analyze high angular-resolution observations with ALMA at $\sim$0.17" corresponding a physical scale of $\sim$400 au in dust continuum and molecular lines. The dust continuum emission reveals a single high-mass protostellar envelope and shows evidence for a marginally resolved continuum source. We detect a rotational line of CH$_3$OH within its $v_{\rm t}$=1 torsionally excited state revealing two bright peaks of emission spatially offset from the dust continuum peak, and exhibiting a distinct velocity component $\pm$4.5 km s$^{-1}$ offset compared to the source $v_{\rm lsr}$. Local thermodynamic equilibrium analysis suggests N(CH$_3$OH)=1.2$-$2$\times$10$^{19}$ cm$^{-2}$, and kinetic temperatures of 160$-$200 K at the position of these peaks. Their velocity shifts correspond well to the expected Keplerian velocity around a central object with 15M$_{\odot}$ consistent with the mass estimate based on the source's bolometric luminosity. We propose a picture where the CH$_3$OH emission peaks trace the accretion shocks around the centrifugal barrier, pinpointing the interaction region between the collapsing envelope and an accretion disk. Because the HC$_3$N $v_{\rm 7}$=1e ($J$=38-37) line shows compact emission, and a velocity pattern consistent with models of Keplerian rotation, we suggest that this could be a new tracer for compact accretion disks around high-mass protostars. The estimated physical properties of the accretion disk suggest a specific angular momentum several times larger than typically observed towards low-mass protostars.

Fragmentation of massive dense molecular clouds is the starting point in the formation of rich clusters and massive stars. Theory and numerical simulations indicate that the population of the fragments (number, mass, diameter, separation) resulting from the gravitational collapse of such clumps is probably regulated by the balance between the magnetic field and the other competitors of self-gravity, in particular turbulence and protostellar feedback. We have observed 11 massive, dense and young star-forming clumps with the Atacama Large Millimeter Array (ALMA) in the thermal dust continuum emission at $\sim 1$~mm with an angular resolution of 0.25 arcseconds with the aim of determining their population of fragments. We find fragments on sub-arcsecond scales in 8 out of the 11 sources. The ALMA images indicate two different fragmentation modes: a dominant fragment surrounded by companions with much smaller mass and size, and many ($\geq 8$) fragments with a gradual change in masses and sizes. On average, the largest number of fragments is found towards the warmer and more massive clumps. Also, the warmer clumps tend to form fragments with larger mass and size. To understand the role of the different physical parameters to regulate the final population of the fragments, we have simulated the collapse of a massive clump of $100$ and $300$ M$_{\odot}$ having different magnetic support. The simulations indicate that: (1) fragmentation is inhibited when the initial turbulence is low, independent of the other physical parameters. (2) a filamentary distribution of the fragments is favoured in a highly magnetised clump. We conclude that the clumps that show many fragments distributed in a filamentary-like structure are likely characterised by a strong magnetic field, while the others are possible also in a weaker magnetic field.

We explore scenarios for the dynamical ejection of stars BN and x from source I in the Kleinmann-Low nebula of the Orion Nebula Cluster (ONC), which is important for being the closest region of massive star formation. This ejection would cause source I to become a close binary or a merger product of two stars. We thus consider binary-binary encounters as the mechanism to produce this event. By running a large suite of $N$-body simulations, we find that it is nearly impossible to match the observations when using the commonly adopted masses for the participants, especially a source I mass of $7\:{\rm{M}}_\odot$. The only way to recreate the event is if source I is more massive, i.e., $\sim20\:{\rm{M}}_\odot$. However, even in this case, the likelihood of reproducing the observed system is low. We discuss the implications of these results for understanding this important star-forming region.

Combining stellar atmospheric parameters, such as effective temperature, surface gravity, and metallicity, with barycentric radial velocity data provides insight into the chemo-dynamics of the MilkyWay and our local Galactic environment. We analyze 3075 stars with spectroscopic data from the Sloan Digital Sky Survey III (SDSS-III) MARVELS radial velocity survey and present atmospheric parameters for 2343 dwarf stars using the spectral indices method, a modified version of the equivalent width method. We present barycentric radial velocities for 2610 stars with uncertainties of $\sim$0.3 km s$^{-1}$. We determine stellar ages using two independent methods and calculate ages for 2343 stars with a maximum-likelihood isochronal age-dating method and for 1694 stars with a Bayesian age-dating method. Using previously published parallax data we compute Galactic orbits and space velocities for 1995 stars to explore stellar populations based on kinematic and age parameters. This study combines good ages and exquisite velocities to explore local chemo-kinematics of the Milky Way, which complements many of the recent studies of giant stars with the APOGEE survey, and we find our results to be in agreement with current chemo-dynamical models of the Milky Way. Particularly, we find from our metallicity distributions and velocity-age relations of a kinematically-defined thin disk that the metal rich end has stars of all ages, even after we clean the sample of highly eccentric stars, suggesting that radial migration plays a key role in the metallicity scatter of the thin disk. All stellar parameters and kinematic data derived in this work are catalogued and published online in machine-readable form.

In this paper, we report our investigation of pulsar scintillation phenomena by monitoring PSR B0355$+$54 at 2.25 GHz for three successive months using \emph{Kunming 40-m radio telescope}. We have measured the dynamic spectrum, the two-dimensional correlation function, and the secondary spectrum. In those observations with high signal-to-noise ratio ($S/N\ge100$), we have detected the scintillation arcs, which are rarely observable using such a small telescope. The sub-microsecond scale width of the scintillation arc indicates that the transverse scale of structures on scattering screen is as compact as AU size. Our monitoring has also shown that both the scintillation bandwidth, timescale, and arc curvature of PSR B0355$+$54 were varying temporally. The plausible explanation would need to invoke multiple-scattering-screen or multiple-scattering-structure scenario that different screens or ray paths dominate the scintillation process at different epochs.

We present SOFIA-upGREAT observations of [CII] emission of Infrared Dark Cloud (IRDC) G035.39-00.33, designed to trace its atomic gas envelope and thus test models of the origins of such clouds. Several velocity components of [CII] emission are detected, tracing structures that are at a wide range of distances in the Galactic plane. We find a main component that is likely associated with the IRDC and its immediate surroundings. This strongest emission component has a velocity similar to that of the $^{13}$CO(2-1) emission of the IRDC, but offset by $\sim3\:{\rm km\:s}^{-1}$ and with a larger velocity width of $\sim9\:{\rm km\:s}^{-1}$. The spatial distribution of the [CII] emission of this component is also offset predominantly to one side of the dense filamentary structure of the IRDC. The CII column density is estimated to be of the order of $\sim10^{17}-10^{18}\,{\rm cm}^{-2}$. We compare these results to the [CII] emission from numerical simulations of magnetized, dense gas filaments formed from giant molecular cloud (GMC) collisions, finding similar spatial and kinematic offsets. These observations and modeling of [CII] add further to the evidence that IRDC G035.39-00.33 has been formed by a process of GMC-GMC collision, which may thus be an important mechanism for initiating star cluster formation.

High-mass star forming regions are typically thought to be dominated by supersonic motions. We present combined Very Large Array and Green Bank Telescope (VLA+GBT) observations of NH$_3$ (1,1) and (2,2) in the infrared dark cloud (IRDC) G035.39-00.33, tracing cold and dense gas down to scales of 0.07 pc. We find that, in contrast to previous similar studies of IRDCs, more than a third of the fitted ammonia spectra show subsonic non-thermal motions (mean line width of 0.71 $\mathrm{km~s^{-1}}$), and the sonic Mach number distribution peaks around $\mathcal{M} = 1$. As possible observational and instrumental biases would only broaden the line profiles, our results provide strong upper limits to the actual value of $\mathcal{M}$, further strengthening our findings of narrow line widths. This finding calls for a reevaluation of the role of turbulent dissipation and subsonic regions in massive-star and cluster formation. Based on our findings in G035.39, we further speculate that the coarser spectral resolution used in the previous VLA NH$_3$ studies may have inhibited the detection of subsonic turbulence in IRDCs. The reduced turbulent support suggests that dynamically important magnetic fields of the 1 mG order would be required to support against possible gravitational collapse. Our results offer valuable input into the theories and simulations that aim to recreate the initial conditions of high-mass star and cluster formation.

Stars are born from dense cores in molecular clouds. Observationally, it is crucial to capture the formation of cores in order to understand the necessary conditions and rate of the star formation process. The {\it Atacama Large Mm/sub-mm Array} (ALMA) is extremely powerful for identifying dense gas structures, including cores, at mm wavelengths via their dust continuum emission. Here we use ALMA to carry out a survey of dense gas and cores in the central region of the massive ($\sim10^5\:M_\odot$) Infrared Dark Cloud (IRDC) G28.37+0.07. The observation consists of a mosaic of 86 pointings of the 12m-array and produces an unprecedented view of the densest structures of this IRDC. In this first paper about this data set, we focus on a comparison between the 1.3 mm continuum emission and a mid-infrared (MIR) extinction map of the IRDC. This allows estimation of the "dense gas" detection probability function (DPF), i.e., as a function of the local mass surface density, $\Sigma$, for various choices of thresholds of mm continuum emission to define "dense gas". We then estimate the dense gas mass fraction, $f_{\rm dg}$, in the central region of the IRDC and, via extrapolation with the DPF and the known $\Sigma$ probability distribution function, to the larger-scale surrounding regions, finding values of about 5\% to 15\% for the fiducial choice of threshold. We argue that this observed dense gas is a good tracer of the protostellar core population and, in this context, estimate a star formation efficiency per free-fall time in the central IRDC region of $\epsilon_{\rm ff}\sim$10\%, with approximately a factor of two systematic uncertainties.

Systems with tightly-packed inner planets (STIPs) are very common. Chatterjee & Tan proposed Inside-Out Planet Formation (IOPF), an in situ formation theory, to explain these planets. IOPF involves sequential planet formation from pebble-rich rings that are fed from the outer disk and trapped at the pressure maximum associated with the dead zone inner boundary (DZIB). Planet masses are set by their ability to open a gap and cause the DZIB to retreat outwards. We present models for the disk density and temperature structures that are relevant to the conditions of IOPF. For a wide range of DZIB conditions, we evaluate the gap opening masses of planets in these disks that are expected to lead to truncation of pebble accretion onto the forming planet. We then consider the evolution of dust and pebbles in the disk, estimating that pebbles typically grow to sizes of a few cm during their radial drift from several tens of AU to the inner, $\lesssim1\:$AU-scale disk. A large fraction of the accretion flux of solids is expected to be in such pebbles. This allows us to estimate the timescales for individual planet formation and entire planetary system formation in the IOPF scenario. We find that to produce realistic STIPs within reasonable timescales similar to disk lifetimes requires disk accretion rates of $\sim10^{-9}\:M_\odot\:{\rm yr}^{-1}$ and relatively low viscosity conditions in the DZIB region, i.e., Shakura-Sunyaev parameter of $\alpha\sim10^{-4}$.

Infrared dark clouds (IRDCs) are thought to be potential hosts of the elusive early phases of high-mass star formation. Here we conduct an in-depth kinematic analysis of one such IRDC, G034.43+00.24 (Cloud F), using high sensitivity and high spectral resolution IRAM-30m N$_2$H$^+$ ($1-0$) and C$^{18}$O ($1-0$) observations. To disentangle the complex velocity structure within this cloud we use Gaussian decomposition and hierarchical clustering algorithms. We find that four distinct coherent velocity components are present within Cloud F. The properties of these components are compared to those found in a similar IRDC, G035.39-00.33 (Cloud H). We find that the components in both clouds have: high densities (inferred by their identification in N$_2$H$^+$), trans-to-supersonic non-thermal velocity dispersions with Mach numbers of $\sim$ $1.5-4$, a separation in velocity of $\sim$3 km s$^{-1}$, and a mean red-shift of $\sim$ 0.3 km s$^{-1}$ between the N$_2$H$^+$ (dense gas) and C$^{18}$O emission (envelope gas). The latter of these could suggest that these clouds share a common formation scenario. We investigate the kinematics of the larger-scale Cloud F structures, using lower-density-tracing $^{13}$CO ($1-0$) observations. A good correspondence is found between the components identified in the IRAM-30m observations and the most prominent component in the$^{13}$CO data. We find that the IRDC Cloud F is only a small part of a much larger structure, which appears to be an inter-arm filament of the Milky Way.

We study the core mass function (CMF) of the massive protocluster G286.21+0.17 with the Atacama Large Millimeter/submillimeter Array via 1.3~mm continuum emission at a resolution of 1.0\arcsec\ (2500~au). We have mapped a field of 5.3\arcmin$\times$5.3\arcmin\ centered on the protocluster clump. We measure the CMF in the central region, exploring various core detection algorithms, which give source numbers ranging from 60 to 125, depending on parameter selection. We estimate completeness corrections due to imperfect flux recovery and core identification via artificial core insertion experiments. For masses $M\gtrsim1\:M_\odot$, the fiducial dendrogram-identified CMF can be fit with a power law of the form ${\rm{d}}N/{\rm{d}}{\rm{log}}M\propto{M}^{-\alpha}$ with $\alpha \simeq1.24\pm0.17$, slightly shallower than, but still consistent with, the index of the Salpeter stellar initial mass function of 1.35. Clumpfind-identified CMFs are significantly shallower with $\alpha\simeq0.64\pm0.13$. While raw CMFs show a peak near $1\:M_\odot$, completeness-corrected CMFs are consistent with a single power law extending down to $\sim 0.5\:M_\odot$, with only a tentative indication of a shallowing of the slope around $\sim1\:M_\odot$. We discuss the implications of these results for star and star cluster formation theories.

We present a continuum radiative transfer model grid for fitting observed spectral energy distributions (SEDs) of massive protostars. The model grid is based on the paradigm of core accretion theory for massive star formation with pre-assembled gravitationally-bound cores as initial conditions. In particular, following the Turbulent Core Model, initial core properties are set primarily by their mass and the pressure of their ambient clump. We then model the evolution of the protostar and its surround structures in a self-consistent way. The model grid contains about 9000 SEDs with 4 free parameters: initial core mass, the mean surface density of the environment, the protostellar mass, and the inclination. The model grid is used to fit observed SEDs via chi^2 minimization, with the foreground extinction additionally estimated. We demonstrate the fitting process and results using the example of massive protostar G35.20- 0.74. Compared with other SED model grids currently used for massive star formation studies, in our model grid, the properties of the protostar and its surrounding structures are more physically connected, which reduces the dimensionality of the parameter spaces and the total number of models. This excludes possible fitting of models that are physically unrealistic or that are not internally self-consistent in the context of the Turbulent Core Model. Thus, this model grid serves not only as a fitting tool to estimate properties of massive protostars, but also as a test of core accretion theory. The SED model grid is publicly released with this paper.

Infrared-Dark Clouds (IRDCs) are cold, dense regions of high (optical and infrared) extinction, believed to be the birthplace of high-mass stars and stellar clusters. The physical mechanisms leading to the formation of these IRDCs are not completely understood and it is thus important to study their molecular gas kinematics and chemical content to search for any signature of the IRDCs formation process. Using the 30m-diameter antenna at the Instituto de Radioastronom\'ia Milim\'etrica, we have obtained emission maps of dense gas tracers (H$^{13}$CO$^{+}$ and HN$^{13}$C) and typical shock tracers (SiO and CH$_3$OH) toward three IRDCs, G028.37+00.07, G034.43+00.24 and G034.77-00.55 (clouds C, F and G, respectively). We have studied the molecular gas kinematics in these clouds and, consistent with previous works toward other IRDCs, the clouds show complex gas kinematics with several velocity-coherent sub-structures separated in velocity space by a few km s$^{-1}$. Correlated with these complex kinematic structures, widespread (parsec-scale) emission of SiO and CH$_3$OH is present in all the three clouds. For clouds C and F, known to be actively forming stars, widespread SiO and CH$_3$OH is likely associated with on-going star formation activity. However, for cloud G, which lacks either 8 $\mu$m or 24 $\mu$m sources and 4.5 $\mu$m H$_2$ shock-excited emission, the detected widespread SiO and CH$_3$OH emission may have originated in a large-scale shock interaction, although a scenario involving a population of low-mass stars driving molecular outflows cannot be fully ruled out.

We present a series of adaptive mesh refinement (AMR) hydrodynamic simulations of flat rotation curve galactic gas disks with a detailed treatment of the interstellar medium (ISM) physics of the atomic to molecular phase transition under the influence of diffuse FUV radiation fields and cosmic ray backgrounds. We explore the effects of different FUV intensities, including a model with a radial gradient designed to mimic the Milky Way. The effects of cosmic rays, including radial gradients in their heating and ionization rates, are also explored. The final simulations in this series achieve $4\:$pc resolution across the $\sim20\:$kpc global disk diameter, with heating and cooling followed down to temperatures of $\sim10\:$K. The disks are evolved for $300\:$Myr, which is enough time for the ISM to achieve a quasi-statistical equilibrium. In particular, the mass fraction of molecular gas stabilizes by $\sim$200 Myr. Additional global ISM properties are analysed. Giant molecular clouds (GMCs) are also identified and the statistical properties of their populations examined. GMCs are tracked as the disks evolve. GMC collisions, which may be a means of triggering star cluster formation, are counted and the rates compared with analytic models. Relatively frequent GMC collision rates are seen in these simulations and their implications for understanding GMC properties, including the driving of internal turbulence, are discussed.

We study the occurrence of spectroscopic binaries in young star-forming regions using the INfrared Spectroscopy of Young Nebulous Clusters(IN-SYNC) survey, carried out in SDSS-III with the APOGEE spectrograph. Multi-epoch observations of thousands of low-mass stars in Orion A, NGC 2264, NGC 1333, IC 348, and the Pleiades have been carried out, yielding H-band spectra with a nominal resolution of R=22,500 for sources with H $\le$ 12 mag. Radial velocity precisions of $\sim$0.3 $km\:s^{-1}$ were achieved, which we use to identify radial velocity variations indicative of undetected companions. We use Monte Carlo simulations to assess the types of spectroscopic binaries to which we are sensitive, finding sensitivity to binaries with orbital periods $< 10^{4}$ d, for stars with $2500 {\rm K} \le T_\mathrm{eff} \le 6000 {\rm K}$ and $\it{v} \sin \it{i}$ $\le$ 100 $km\:s^{-1}$. Using Bayesian inference, we find evidence for a decline in the spectroscopic binary fraction, by a factor of 3-4 from the age of our pre-main-sequence sample to the Pleiades age . The significance of this decline is weakened if spot-induced radial-velocity jitter is strong in the sample, and is only marginally significant when comparing any one of the pre-main-sequence clusters against the Pleiades. However, the same decline in both sense and magnitude is found for each of the five pre-main-sequence clusters, and the decline reaches statistical significance of greater than 95% confidence when considering the pre-main-sequence clusters jointly. Our results suggest that dynamical processes disrupt the widest spectroscopic binaries ($P_{\rm orb} \approx 10^3 - 10^4$ d) as clusters age, indicating that this occurs early in the stars' evolution, while they still reside within their nascent clusters.

I review massive star formation in our Galaxy, focusing on initial conditions in Infrared Dark Clouds (IRDCs), including the search for massive pre-stellar cores (PSCs), and modeling of later stages of massive protostars, i.e., hot molecular cores (HMCs). I highlight how developments in astrochemistry, coupled with rapidly improving theoretical/computational and observational capabilities are helping to improve our understanding of the complex process of massive star formation.

We investigate collisions between giant molecular clouds (GMCs) as potential generators of their internal turbulence. Using magnetohydrodynamic (MHD) simulations of self-gravitating, magnetized, turbulent, GMCs, we compare kinematic and dynamic properties of dense gas structures formed when such clouds collide compared to those that form in non-colliding clouds as self-gravity overwhelms decaying turbulence. We explore the nature of turbulence in these structures via distribution functions of density, velocity dispersions, virial parameters, and momentum injection. We find that the dense clumps formed from GMC collisions have higher effective Mach number, greater overall velocity dispersions, sustain near-virial equilibrium states for longer times, and are the conduit for injection of turbulent momentum into high density gas at high rates.

We present calculations of molecular, atomic and ionic line emission from simulations of giant molecular cloud (GMC) collisions. We post-process snapshots of the magneto-hydrodynamical simulations presented in an earlier paper in this series by Wu et al. (2017) of colliding and non-colliding GMCs. Using photodissociation region (PDR) chemistry and radiative transfer we calculate the level populations and emission properties of $^{12}$CO $J=1-0$, [CI] $^3{\rm P}_1\rightarrow{^3{\rm P}}_0$ at $609\,\mu$m, [CII] $158\,\mu$m and [OI] $^3{\rm P}_1\rightarrow{^3{\rm P}}_0$ transition at $63\,\mu$m. From integrated intensity emission maps and position-velocity diagrams, we find that fine-structure lines, particularly the [CII] $158\,\mu$m, can be used as a diagnostic tracer for cloud-cloud collision activity. These results hold even in more evolved systems in which the collision signature in molecular lines has been diminished.

In this series of papers, we model the formation and evolution of the photoionized region and its observational signatures during massive star formation. Here we focus on the early break out of the photoionized region into the outflow cavity. Using results of 3-D magnetohydrodynamic-outflow simulations and protostellar evolution calculations, we perform post-processing radiative-transfer. The photoionized region first appears at a protostellar mass of 10Msun in our fiducial model, and is confined to within 10-100AU by the dense inner outflow, similar to some observed very small hypercompact HII regions. Since the ionizing luminosity of the massive protostar increases dramatically as Kelvin-Helmholz (KH) contraction proceeds, the photoionized region breaks out to the entire outflow region in <10,000yr. Accordingly, the radio free-free emission brightens significantly in this stage. In our fiducial model, the radio luminosity at 10 GHz changes from 0.1 mJy kpc2 at m=11Msun to 100 mJy kpc2 at 16Msun, while the infrared luminosity increases by less than a factor of two. The radio spectral index also changes in the break-out phase from the optically thick value of 2 to the partially optically thin value of 0.6. Additionally, we demonstrate that short-timescale variation in free-free flux would be induced by an accretion burst. The outflow density is enhanced in the accretion burst phase, which leads to a smaller ionized region and weaker free-free emission. The radio luminosity may decrease by one order of magnitude during such bursts, while the infrared luminosity is much less affected, since internal protostellar luminosity dominates over accretion luminosity after KH contraction starts. Such variability may be observable on timescales as short 10-100 yr, if accretion bursts are driven by disk instabilities.

The early evolution during massive star cluster formation is still uncertain. Observing embedded clusters at their earliest stages of formation can provide insight into the spatial and temporal distribution of the stars and thus probe different star cluster formation models. We present near-infrared imaging of an 8'*13'(5.4pc*8.7pc) region around the massive infalling clump G286.21+0.17(also known as BYF73). The stellar content across the field is determined and photometry is derived in order to { obtain} stellar parameters for the cluster members. We find evidence for some sub-structure (on scales less than a pc diameter) within the region with apparently at least three different sub-clusters associated with the molecular clump based on differences in extinction and disk fractions. At the center of the clump we identify a deeply embedded sub-cluster. Near-infrared excess is detected for 39-44% in the two sub-clusters associated with molecular material and 27% for the exposed cluster. Using the disk excess as a proxy for age this suggests the clusters are very young. The current total stellar mass is estimated to be at least 200 Msun. The molecular core hosts a rich population of pre-main sequence stars. There is evidence for multiple events of star formation both in terms of the spatial distribution within the star forming region and possibly from the disk frequency.

We investigate the role of ambipolar diffusion (AD) in collisions between magnetized giant molecular clouds (GMCs), which may be an important mechanism for triggering star cluster formation. Three dimensional simulations of GMC collisions are performed using a version of the Enzo magnetohydrodynamics code that has been extended to include AD. The resistivities are calculated using the 31-species chemical model of Wu et al. (2015). We find that in the weak-field, $10\:{\rm \mu G}$ case, AD has only a modest effect on the dynamical evolution during the collision. However, for the stronger-field, $30\:{\rm \mu G}$ case involving near-critical clouds, AD results in formation of dense cores in regions where collapse is otherwise inhibited. The overall efficiency of formation of cores with $n_{\rm H}\geq10^{6}\:{\rm cm}^{-3}$ in these simulations is increases from about 0.2% to 2% once AD is included, comparable to observed values in star-forming GMCs. The gas around these cores typically has relatively slow infall at speeds that are a modest fraction of the free-fall speed.

We have found many Earth-sized worlds but we have no way of determining if their surfaces are Earth-like. This makes it impossible to quantitatively compare habitability, and pretending we can risks damaging the field.

The kinematics and dynamics of young stellar populations enable us to test theories of star formation. With this aim, we continue our analysis of the SDSS-III/APOGEE IN-SYNC survey, a high resolution near infrared spectroscopic survey of young clusters. We focus on the Orion A star-forming region, for which IN-SYNC obtained spectra of $\sim2700$ stars. In Paper IV we used these data to study the young stellar population. Here we study the kinematic properties through radial velocities ($v_r$). The young stellar population remains kinematically associated with the molecular gas, following a $\sim10\:{\rm{km\:s}}^{-1}$ gradient along filament. However, near the center of the region, the $v_r$ distribution is slightly blueshifted and asymmetric; we suggest that this population, which is older, is slightly in foreground. We find evidence for kinematic subclustering, detecting statistically significant groupings of co-located stars with coherent motions. These are mostly in the lower-density regions of the cloud, while the ONC radial velocities are smoothly distributed, consistent with it being an older, more dynamically evolved cluster. The velocity dispersion $\sigma_v$ varies along the filament. The ONC appears virialized, or just slightly supervirial, consistent with an old dynamical age. Here there is also some evidence for on-going expansion, from a $v_r$--extinction correlation. In the southern filament, $\sigma_v$ is $\sim2$--$3$ times larger than virial in the L1641N region, where we infer a superposition along the line of sight of stellar sub-populations, detached from the gas. On the contrary, $\sigma_v$ decreases towards L1641S, where the population is again in agreement with a virial state.

Aims. Infrared dark clouds represent the earliest stages of high-mass star formation. Detailed observations of their physical conditions on all physical scales are required to improve our understanding of their role in fueling star formation. Methods. We investigate the large-scale structure of the IRDC G035.39-00.33, probing the dense gas with the classical ammonia thermometer. This allows us to put reliable constraints on the temperature of the extended, pc-scale dense gas reservoir and to probe the magnitude of its non-thermal motions. Available far-infrared observations can be used in tandem with the observed ammonia emission to estimate the total gas mass contained in G035.39-00.33. Results. We identify a main velocity component as a prominent filament, manifested as an ammonia emission intensity ridge spanning more than 6 pc, consistent with the previous studies on the Northern part of the cloud. A number of additional line-of-sight components are found, and a large scale, linear velocity gradient of ~0.2 km s$^{-1}$ pc$^{-1}$ is found along the ridge of the IRDC. In contrast to the dust temperature map, an ammonia-derived kinetic temperature map, presented for the entirety of the cloud, reveals local temperature enhancements towards the massive protostellar cores. We show that without properly accounting for the line of sight contamination, the dust temperature is 2-3 K larger than the gas temperature measured with NH$_3$. Conclusions. While both the large scale kinematics and temperature structure are consistent with that of starless dark filaments, the kinetic gas temperature profile on smaller scales is suggestive of tracing the heating mechanism coincident with the locations of massive protostellar cores.

We present an overview and first results of the Stratospheric Observatory For Infrared Astronomy Massive (SOMA) Star Formation Survey, which is using the FORCAST instrument to image massive protostars from $\sim10$--$40\:\rm{\mu}\rm{m}$. These wavelengths trace thermal emission from warm dust, which in Core Accretion models mainly emerges from the inner regions of protostellar outflow cavities. Dust in dense core envelopes also imprints characteristic extinction patterns at these wavelengths, causing intensity peaks to shift along the outflow axis and profiles to become more symmetric at longer wavelengths. We present observational results for the first eight protostars in the survey, i.e., multiwavelength images, including some ancillary ground-based MIR observations and archival {\it{Spitzer}} and {\it{Herschel}} data. These images generally show extended MIR/FIR emission along directions consistent with those of known outflows and with shorter wavelength peak flux positions displaced from the protostar along the blueshifted, near-facing sides, thus confirming qualitative predictions of Core Accretion models. We then compile spectral energy distributions and use these to derive protostellar properties by fitting theoretical radiative transfer models. Zhang and Tan models, based on the Turbulent Core Model of McKee and Tan, imply the sources have protostellar masses $m_*\sim10$--50$\:M_\odot$ accreting at $\sim10^{-4}$--$10^{-3}\:M_\odot\:{\rm{yr}}^{-1}$ inside cores of initial masses $M_c\sim30$--500$\:M_\odot$ embedded in clumps with mass surface densities $\Sigma_{\rm{cl}}\sim0.1$--3$\:{\rm{g\:cm}^{-2}}$. Fitting Robitaille et al. models typically leads to slightly higher protostellar masses, but with disk accretion rates $\sim100\times$ smaller. We discuss reasons for these differences and overall implications of these first survey results for massive star formation theories.

We present radial velocity measurements for 70 high confidence, and 34 potential binary systems in fields containing the Perseus Molecular Cloud, Pleiades, NGC 2264, and the Orion A star forming region. 18 of these systems have been previously identified as binaries in the literature. Candidate double-lined spectroscopic binaries (SB2s) are identified by analyzing the cross-correlation functions (CCFs) computed during the reduction of each APOGEE spectrum. We identify sources whose CCFs are well fit as the sum of two Lorentzians as likely binaries, and provide an initial characterization of the system based on the radial velocities indicated by that dual fit. For systems observed over several epochs, we present mass ratios and systemic velocities; for two systems with observations on eight or more epochs, and which meet our criteria for robust orbital coverage, we derive initial orbital parameters. The distribution of mass ratios for multi-epoch sources in our sample peaks at q=1, but with a significant tail toward lower q values. Tables reporting radial velocities, systemic velocities, and mass ratios are provided online. We discuss future improvements to the radial velocity extraction method we employ, as well as limitations imposed by the number of epochs currently available in the APOGEE database. The Appendix contains brief notes from the literature on each system in the sample, and more extensive notes for select sources of interest.

Star formation from the interstellar medium of galactic disks is a basic process controlling the evolution of galaxies. Understanding the star formation rate in a local patch of a disk with a given gas mass is thus an important challenge for theoretical models. Here we simulate a kiloparsec region of a disk, following the evolution of self-gravitating molecular clouds down to subparsec scales, as they form stars that then inject feedback energy by dissociating and ionizing UV photons and supernova explosions. We assess the relative importance of each feedback mechanism. We find that $\rm H_2$-dissociating feedback results in the largest absolute reduction in star formation compared to the run with no feedback. Subsequently adding photoionization feedback produces a more modest reduction. Our fiducial models that combine all three feedback mechanisms yield, without fine-tuning, star formation rates that are in excellent agreement with observations, with $\rm H_2$-dissociating photons playing a crucial role. Models that only include supernova feedback---a common method in galaxy evolution simulations---settle to similar star formation rates, but with very different temperature and chemical states of the gas, and with very different spatial distributions of young stars.

We study giant molecular cloud (GMC) collisions and their ability to trigger star cluster formation. We further develop our three dimensional magnetized, turbulent, colliding GMC simulations by implementing star formation sub-grid models. Two such models are explored: (1) "Density-Regulated," i.e., fixed efficiency per free-fall time above a set density threshold; (2) "Magnetically-Regulated," i.e., fixed efficiency per free-fall time in regions that are magnetically supercritical. Variations of parameters associated with these models are also explored. In the non-colliding simulations, the overall level of star formation is sensitive to model parameter choices that relate to effective density thresholds. In the GMC collision simulations, the final star formation rates and efficiencies are relatively independent of these parameters. Between non-colliding and colliding cases, we compare the morphologies of the resulting star clusters, properties of star-forming gas, time evolution of the star formation rate (SFR), spatial clustering of the stars, and resulting kinematics of the stars in comparison to the natal gas. We find that typical collisions, by creating larger amounts of dense gas, trigger earlier and enhanced star formation, resulting in 10 times higher SFRs and efficiencies. The star clusters formed from GMC collisions show greater spatial sub-structure and more disturbed kinematics.

We have measured astrometry for members of the Orion Nebula Cluster with images obtained in 2015 with the Wide Field Camera 3 on board the Hubble Space Telescope. By comparing those data to previous measurements with NICMOS on Hubble in 1998, we have discovered that a star in the Kleinmann-Low Nebula, source x from Lonsdale et al. (1982), is moving with an unusually high proper motion of 29 mas/yr, which corresponds to 55 km/s at the distance of Orion. Previous radio observations have found that three other stars in the Kleinmann-Low Nebula (BN and sources I and n) have high proper motions (5-14 mas/yr) and were near a single location ~540 years ago, and thus may have been members of a multiple system that dynamically decayed. The proper motion of source x is consistent with ejection from that same location 540 years ago, which provides strong evidence that the dynamical decay did occur and that the runaway star BN originated in the Kleinmann-Low Nebula rather than the nearby Trapezium cluster. However, our constraint on the motion of source n is significantly smaller than the most recent radio measurement, which indicates that it did not participate in the event that ejected the other three stars.

The early evolution of massive cluster progenitors is poorly understood. We investigate the fragmentation properties from 0.3 pc to 0.06 pc scales of a homogenous sample of infrared-quiet massive clumps within 4.5 kpc selected from the ATLASGAL survey. Using the ALMA 7m array we detect compact dust continuum emission towards all targets, and find that fragmentation, at these scales, is limited. The mass distribution of the fragments uncovers a large fraction of cores above 40 $M_\odot$, corresponding to massive dense cores (MDCs) with masses up to ~400 $M_\odot$. 77 % of the clumps contain at most 3 MDCs per clump, and we also reveal single clumps/MDCs. The most massive cores are formed within the more massive clumps, and a high concentration of mass on small scales reveals a high core formation efficiency. The mass of MDCs highly exceeds the local thermal Jeans-mass, and observational evidence is lacking for a sufficiently high level of turbulence or strong enough magnetic fields to keep the most massive MDCs in equilibrium. If already collapsing, the observed fragmentation properties with a high core formation efficiency are consistent with the collapse setting in at parsec scales.

We investigate the formation and early evolution of star clusters assuming that they form from a turbulent starless clump of given mass bounded inside a parent self-gravitating molecular cloud characterized by a particular mass surface density. As a first step we assume instantaneous star cluster formation and gas expulsion. We draw our initial conditions from observed properties of starless clumps. We follow the early evolution of the clusters up to 20 Myr, investigating effects of different star formation efficiencies, primordial binary fractions and eccentricities and primordial mass segregation levels. We investigate clumps with initial masses of $M_{\rm cl}=3000\:{\rm M}_\odot$ embedded in ambient cloud environments with mass surface densities, $\Sigma_{\rm cloud}=0.1$ and $1\:{\rm g\:cm^{-2}}$. We show that these models of fast star cluster formation result, in the fiducial case, in clusters that expand rapidly, even considering only the bound members. Clusters formed from higher $\Sigma_{\rm cloud}$ environments tend to expand more quickly, so are soon larger than clusters born from lower $\Sigma_{\rm cloud}$ conditions. To form a young cluster of a given age, stellar mass and mass surface density, these models need to assume a parent molecular clump that is many times denser, which is unrealistic compared to observed systems. We also show that in these models the initial binary properties are only slightly modified by interactions, meaning that binary properties, e.g., at 20 Myr, are very similar to those at birth. With this study we set up the basis of future work where we will investigate more realistic models of star formation compared to this instantaneous, baseline case.

We study feedback during massive star formation using semi-analytic methods, considering the effects of disk winds, radiation pressure, photoevaporation and stellar winds, while following protostellar evolution in collapsing massive gas cores. We find that disk winds are the dominant feedback mechanism setting star formation efficiencies (SFEs) from initial cores of ~0.3-0.5. However, radiation pressure is also significant to widen the outflow cavity causing reductions of SFE compared to the disk-wind only case, especially for >100Msun star formation at clump mass surface densities Sigma<0.3g/cm2. Photoevaporation is of relatively minor importance due to dust attenuation of ionizing photons. Stellar winds have even smaller effects during the accretion stage. For core masses Mc~10-1000Msun and Sigma~0.1-3g/cm2, we find the overall SFE to be 0.31(Rc/0.1pc)^{-0.39}, potentially a useful sub-grid star-formation model in simulations that can resolve pre-stellar core radii, Rc=0.057(Mc/60Msun)^{1/2}(Sigma/g/cm2)^{-1/2}pc. The decline of SFE with Mc is gradual with no evidence for a maximum stellar-mass set by feedback processes up to stellar masses of ~300Msun. We thus conclude that the observed truncation of the high-mass end of the IMF is shaped mostly by the pre-stellar core mass function or internal stellar processes. To form massive stars with the observed maximum masses of ~150-300Msun, initial core masses need to be >500-1000Msun. We also apply our feedback model to zero-metallicity primordial star formation, showing that, in the absence of dust, photoevaporation staunches accretion at ~50Msun. Our model implies radiative feedback is most significant at metallicities ~10^{-2}Zsun, since both radiation pressure and photoevaporation are effective in this regime.

We investigate giant molecular cloud (GMCs) collisions and their ability to induce gravitational instability and thus star formation. This mechanism may be a major driver of star formation activity in galactic disks. We carry out a series of three dimensional, magnetohydrodynamics (MHD), adaptive mesh refinement (AMR) simulations to study how cloud collisions trigger formation of dense filaments and clumps. Heating and cooling functions are implemented based on photo-dissociation region (PDR) models that span the atomic to molecular transition and can return detailed diagnostic information. The clouds are initialized with supersonic turbulence and a range of magnetic field strengths and orientations. Collisions at various velocities and impact parameters are investigated. Comparing and contrasting colliding and non-colliding cases, we characterize morphologies of dense gas, magnetic field structure, cloud kinematic signatures, and cloud dynamics. We present key observational diagnostics of cloud collisions, especially: relative orientations between magnetic fields and density structures, like filaments; $^{13}$CO($J$=2-1), $^{13}$CO($J$=3-2), and $^{12}$CO($J$=8-7) integrated intensity maps and spectra; and cloud virial parameters. We compare these results to observed Galactic clouds.

We carry out an ALMA $\rm N_2D^+$(3-2) and 1.3~mm continuum survey towards 32 high mass surface density regions in seven Infrared Dark Clouds with the aim of finding massive starless cores, which may be the initial conditions for the formation of massive stars. Cores showing strong $\rm N_2D^+$(3-2) emission are expected to be highly deuterated and indicative of early, potentially pre-stellar stages of star formation. We also present maps of these regions in ancillary line tracers, including C$^{18}$O(2-1), DCN(3-2) and DCO$^+$(3-2). Over 100 $\rm N_2D^+$ cores are identified with our newly developed core-finding algorithm based on connected structures in position-velocity space. The most massive core has $\sim70\:M_\odot$ (potentially $\sim170\:M_\odot$) and so may be representative of the initial conditions for massive star formation. The existence and dynamical properties of such cores constrain massive star formation theories. We measure the line widths and thus velocity dispersion of six of the cores with strongest $\rm N_2D^+$(3-2) line emission, finding results that are generally consistent with virial equilibrium of pressure confined cores.

High levels of deuterium fraction in N$_2$H$^+$ are observed in some pre-stellar cores. Single-zone chemical models find that the timescale required to reach observed values ($D_{\rm frac}^{{\rm N}_2{\rm H}^+} \equiv {\rm N}_2{\rm D}^+/{\rm N}_2{\rm H}^+ \gtrsim 0.1$) is longer than the free-fall time, possibly ten times longer. Here, we explore the deuteration of turbulent, magnetized cores with 3D magnetohydrodynamics simulations. We use an approximate chemical model to follow the growth in abundances of N$_2$H$^+$ and N$_2$D$^+$. We then examine the dynamics of the core using each tracer for comparison to observations. We find that the velocity dispersion of the core as traced by N$_2$D$^+$ appears slightly sub-virial compared to predictions of the Turbulent Core Model of McKee & Tan, except at late times just before the onset of protostar formation. By varying the initial mass surface density, the magnetic energy, the chemical age, and the ortho-to-para ratio of H$_2$, we also determine the physical and temporal properties required for high deuteration. We find that low initial ortho-to-para ratios ($\lesssim 0.01$) and/or multiple free-fall times ($\gtrsim 3$) of prior chemical evolution are necessary to reach the observed values of deuterium fraction in pre-stellar cores.

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.

High-mass stars shape the interstellar medium in galaxies, and yet, largely because the initial conditions are poorly constrained, we do not know how they form. One possibility is that high-mass stars and star clusters form at the junction of filamentary networks, referred to as 'hubs'. In this letter we present the complex anatomy of a protocluster hub within an Infrared Dark Cloud (IRDC), G035.39-00.33, believed to be in an early phase of its evolution. We use high-angular resolution ($\{\theta_{\rm maj}, \theta_{\rm min}\}=\{1.4{\rm arcsec}, 0.8{\rm arcsec}\}\sim\{0.02{\rm pc}, 0.01{\rm pc}\}$) and high-sensitivity ($0.2$mJy beam$^{-1}$; $\sim0.2$M$_{\odot}$) 1.07 mm dust continuum observations from the Atacama Large Millimeter Array (ALMA) to identify a network of narrow, $0.028\,\pm\,0.005$pc wide, filamentary structures. These are a factor of $\gtrsim3$ narrower than the proposed 'quasi-universal' $\sim0.1$pc width of interstellar filaments. Additionally, 28 compact objects are reported, spanning a mass range $0.3{\rm M_{\odot}}<M_{\rm c}<10.4{\rm M_{\odot}}$. This indicates that at least some low-mass objects are forming coevally with more massive counterparts. Comparing to the popular 'bead-on-a-string' analogy, the protocluster hub is poorly represented by a monolithic clump embedded within a single filament. Instead, it comprises multiple intra-hub filaments, each of which retains its integrity as an independent structure and possesses its own embedded core population.

We study the probability distribution function (PDF) of mass surface densities, $\Sigma$, of infrared dark cloud (IRDC) G028.37+00.07 and its surrounding giant molecular cloud. This PDF constrains the physical processes, such as turbulence, magnetic fields and self-gravity, that are expected to be controlling cloud structure and star formation activity. The chosen IRDC is of particular interest since it has almost 100,000 solar masses within a radius of 8 parsecs, making it one of the most massive, dense molecular structures known and is thus a potential site for the formation of a "super star cluster." We study $\Sigma$ in two ways. First, we use a combination of NIR and MIR extinction maps that are able to probe the bulk of the cloud structure up to $\Sigma\sim1\:{\rm{g\:cm}^{-2}}\:$($A_V\simeq200$ mag). Second, we study the FIR and sub-mm dust continuum emission from the cloud utilizing Herschel PACS and SPIRE images and paying careful attention to the effects of foreground and background contamination. We find that the PDFs from both methods, applied over a $\sim20^\prime$(30 pc)-sized region that contains $\simeq1.5\times10^5\:M_\odot$ and encloses a minimum closed contour with $\Sigma\simeq0.013\:{\rm{g\:cm}^{-2}}\:$($A_V\simeq3$ mag), shows a log-normal shape with the peak measured at $\Sigma\simeq0.021\:{\rm{g\:cm}^{-2}}\:$($A_V\simeq4.7$ mag). There is tentative evidence for the presence of a high-$\Sigma$ power law tail that contains from $\sim3\%$ to 8\% of the mass of the cloud material. We discuss the implications of these results for the physical processes occurring in this cloud.

Massive stars, multiple stellar systems and clusters are born from the gravitational collapse of massive dense gaseous clumps, and the way these systems form strongly depends on how the parent clump fragments into cores during collapse. Numerical simulations show that magnetic fields may be the key ingredient in regulating fragmentation. Here we present ALMA observations at ~0.25'' resolution of the thermal dust continuum emission at ~278 GHz towards a turbulent, dense, and massive clump, IRAS 16061-5048c1, in a very early evolutionary stage. The ALMA image shows that the clump has fragmented into many cores along a filamentary structure. We find that the number, the total mass and the spatial distribution of the fragments are consistent with fragmentation dominated by a strong magnetic field. Our observations support the theoretical prediction that the magnetic field plays a dominant role in the fragmentation process of massive turbulent clump.

We present 3.7 arcsec (~0.05 pc) resolution 3.2 mm dust continuum observations from the IRAM PdBI, with the aim of studying the structure and fragmentation of the filamentary Infrared Dark Cloud G035.39-00.33. The continuum emission is segmented into a series of 13 quasi-regularly spaced (~0.18pc) cores, following the major axis of the IRDC. We compare the spatial distribution of the cores with that predicted by theoretical work describing the fragmentation of hydrodynamic fluid cylinders, finding a significant (a factor of ~8) discrepancy between the two. Our observations are consistent with the picture emerging from kinematic studies of molecular clouds suggesting that the cores are harboured within a complex network of independent sub-filaments. This result emphasizes the importance of considering the underlying physical structure, and potentially, dynamically important magnetic fields, in any fragmentation analysis. The identified cores exhibit a range in (peak) beam-averaged column density ($3.6{\rm x}10^{23}{\rm cm}^{-2}<N_{H,c}<8.0{\rm x}10^{23}{\rm cm}^{-2}$), mass ($8.1M_{\odot}<M_{c}<26.1M_{\odot}$), and number density ($6.1{\rm x}10^{5}{\rm cm}^{-3}<n_{H, c, eq}<14.7{\rm x}10^{5}{\rm cm}^{-3}$). Two of these cores, dark in the mid-infrared, centrally-concentrated, monolithic (with no traceable substructure at our PdBI resolution), and with estimated masses of the order ~20-25$M_{\odot}$, are good candidates for the progenitors of intermediate-to-high-mass stars. Virial parameters span a range $0.2<\alpha_{\rm vir}<1.3$. Without additional support, possibly from dynamically important magnetic fields with strengths of the order 230$\mu$G<B<670$\mu$G, the cores are susceptible to gravitational collapse. These results may imply a multilayered fragmentation process, which incorporates the formation of sub-filaments, embedded cores, and the possibility of further fragmentation.

The Thirty Meter Telescope (TMT) first light instrument IRIS (Infrared Imaging Spectrograph) will complete its preliminary design phase in 2016. The IRIS instrument design includes a near-infrared (0.85 - 2.4 micron) integral field spectrograph (IFS) and imager that are able to conduct simultaneous diffraction-limited observations behind the advanced adaptive optics system NFIRAOS. The IRIS science cases have continued to be developed and new science studies have been investigated to aid in technical performance and design requirements. In this development phase, the IRIS science team has paid particular attention to the selection of filters, gratings, sensitivities of the entire system, and science cases that will benefit from the parallel mode of the IFS and imaging camera. We present new science cases for IRIS using the latest end-to-end data simulator on the following topics: Solar System bodies, the Galactic center, active galactic nuclei (AGN), and distant gravitationally-lensed galaxies. We then briefly discuss the necessity of an advanced data management system and data reduction pipeline.

Giant molecular clouds contain supersonic turbulence that can locally heat small fractions of gas to over 100 K. We run shock models for low-velocity, C-type shocks propagating into gas with densities between 10^3 and 10^5 cm^(-3) and find that CO lines are the most important cooling lines. Comparison to photodissociation region (PDR) models indicates that mid-J CO lines (J = 8-7 and higher) should be dominated by emission from shocked gas. In Papers I and II we presented CO J = 3-2, 8-7, and 9-8 observations toward four primarily quiescent clumps within infrared dark clouds. Here we fit PDR models to the combined spectral line energy distributions and show that the PDR models that best fit the low-J CO emission underpredict the mid-J CO emission by orders of magnitude, strongly hinting at a hot gas component within these clumps. The low-J CO data clearly show that the integrated intensities of both the CO J = 8-7 and 9-8 lines are anomalously high, such that the line ratio can be used to characterize the hot gas component. Shock models are reasonably consistent with the observed mid-J CO emission, with models with densities near 10^(4.5) cm^(-3) providing the best agreement. Where this mid-J CO is detected, the mean volume filling factor of the hot gas is 0.1%. Much of the observed mid-J CO emission, however, is also associated with known protostars and may be due to protostellar feedback.

We study deuterium fractionation in two massive starless/early-stage cores C1-N and C1-S in Infrared Dark Cloud (IRDC) G028.37+00.07, first identified by Tan et al. (2013) with ALMA. Line emission from multiple transitions of $\rm N_2H^+$ and $\rm N_2D^+$ were observed with the ALMA, CARMA, SMA, JCMT, NRO 45m and IRAM 30m telescopes. By simultaneously fitting the spectra, we estimate the excitation conditions and deuterium fraction, $D_{\rm frac}^{\rm N_2H^+} \equiv [\rm N_2D^+]/[N_2H^+]$, with values of $D_{\rm frac}^{\rm N_2H^+} \simeq 0.2$--$0.7$, several orders of magnitude above the cosmic [D]/[H] ratio. Additional observations of o-H$_2$D$^+$ are also presented that help constrain the ortho-to-para ratio of $\rm H_2$, which is a key quantity affecting the degree of deuteration. We then present chemodynamical modeling of the two cores, exploring especially the implications for the collapse rate relative to free-fall, $\alpha_{\rm ff}$. In order to reach the high level of observed deuteration of $\rm N_2H^+$, we find that the most likely evolutionary history of the cores involves collapse at a relatively slow rate, $\lesssim1/10$th of free-fall.

We present ALMA follow-up observations of two massive, early-stage core candidates, C1-N & C1-S, in Infrared Dark Cloud (IRDC) G028.37+00.07, which were previously identified by their N2D+(3-2) emission and show high levels of deuteration of this species. The cores are also dark at far infrared wavelengths up to ~100 microns. We detect 12CO(2-1) from a narrow, highly-collimated bipolar outflow that is being launched from near the center of the C1-S core, which is also the location of the peak 1.3mm dust continuum emission. This protostar, C1-Sa, has associated dense gas traced by C18O(2-1) and DCN(3-2), from which we estimate it has a radial velocity that is near the center of the range exhibited by the C1-S massive core. A second outflow-driving source is also detected within the projected boundary of C1-S, but appears to be at a different radial velocity. After considering properties of the outflows, we conclude C1-Sa is a promising candidate for an early-stage massive protostar and as such it shows that these early phases of massive star formation can involve highly ordered outflow, and thus accretion, processes, similar to models developed to explain low-mass protostars.

Infrared Dark Clouds (IRDCs) are cold, dense regions that are usually found within Giant Molecular Clouds (GMCs). Ongoing star formation within IRDCs is typically still deeply embedded within the surrounding molecular gas. Characterising the properties of relatively quiescent IRDCs may therefore help us to understand the earliest phases of the star formation process. Studies of local molecular clouds have revealed that deuterated species are enhanced in the earliest phases of star formation. In this paper we test this towards IRDC G035.39-00.33. We present an 80 arcsec by 140 arcsec map of the J=2-1 transition of N2D+, obtained with the IRAM-30m telescope. We find that N2D+ is widespread throughout G035.39-00.33. Complementary observations of N2H+(1-0) are used to estimate the deuterium fraction, N(N2D+)/N(N2H+). We report a mean deuterium fraction of 0.04+-0.01, with a maximum of 0.09+-0.02. The mean deuterium fraction is ~3 orders of magnitude greater than the interstellar [D]/[H] ratio. High angular resolution observations are required to exclude beam dilution effects of compact deuterated cores. Using chemical modelling, we find that the average observed values of the deuterium fraction are in agreement with an equilibrium deuterium fraction, given the general properties of the cloud. This implies that the IRDC is at least ~3Myr old, which is ~8 times longer than the mean free-fall time of the observed deuterated region.

We present an evolutionary sequence of models of the photoionized disk-wind outflow around forming massive stars based on the Core Accretion model. The outflow is expected to be the first structure to be ionized by the protostar and can confine the expansion of the HII region, especially in lateral directions in the plane of the accretion disk. The ionizing luminosity increases as Kelvin-Helmholz contraction proceeds, and the HII region is formed when the stellar mass reaches ~10-20Msun depending on the initial cloud core properties. Although some part of outer disk surface remains neutral due to shielding by the inner disk and the disk wind, almost the whole of the outflow is ionized in 1e3-1e4 yr after initial HII region formation. Having calculated the extent and temperature structure of the HII region within the immediate protostellar environment, we then make predictions for the strength of its free-free continuum and recombination line emission. The free-free radio emission from the ionized outflow has a flux density of ~(20-200)x(nu/10GHz)^p mJy for a source at a distance of 1 kpc with a spectral index p~0.4-0.7, and the apparent size is typically ~500AU at 10GHz. The H40alpha line profile has a width of about 100km/s. These properties of our model are consistent with observed radio winds and jets around forming massive protostars.

Infrared dark clouds are kinematically complex molecular structures in the interstellar medium that can host sites of massive star formation. We present 4 square arcminute maps of the 12CO, 13CO, and C18O J = 3 to 2 lines from selected locations within the C and F (G028.37+00.07 and G034.43+00.24) infrared dark clouds (IRDCs), as well as single pointing observations of the 13CO and C18O J = 2 to 1 lines towards three cores within these clouds. We derive CO gas temperatures throughout the maps and find that CO is significantly frozen out within these IRDCs. We find that the CO depletion tends to be the highest near column density peaks, with maximum depletion factors between 5 and 9 in IRDC F and between 16 and 31 in IRDC C. We also detect multiple velocity components and complex kinematic structure in both IRDCs. Therefore, the kinematics of IRDCs seem to point to dynamically evolving structures yielding dense cores with considerable depletion factors.

The origin of massive stars is a fundamental open issue in modern astrophysics. Pre-ALMA interferometric studies reveal precursors to early B to late O type stars with collapsing envelopes of 15-20 M$_\odot$ on 1000-3000 AU size-scales. To search for more massive envelopes we selected the most massive nearby young clumps from the ATLASGAL survey to study their protostellar content with ALMA. Our first results using the intermediate scales revealed by the ALMA ACA array providing 3-5" angular resolution, corresponding to $\sim$0.05-0.1 pc size-scales, reveals a sample of compact objects. These massive dense cores are on average two-times more massive than previous studies of similar types of objects. We expect that once the full survey is completed, it will provide a comprehensive view on the origin of the most massive stars.

We present the results of the SDSS APOGEE INfrared Spectroscopy of Young Nebulous Clusters program (IN-SYNC) survey of the Orion A molecular cloud. This survey obtained high resolution near infrared (NIR) spectroscopy of about 2700 young pre-main sequence stars throughout the region, acquired across five distinct fields spanning 6deg field of view (FOV). With these spectra, we have measured accurate stellar parameters (T_eff, log g, v sin i) and extinctions, and placed the sources in the Hertzsprung-Russel Diagram (HRD). We have also extracted radial velocities for the kinematic characterization of the population. We compare our measurements with literature results for a sub-sample of targets in order to assess the performances and accuracy of the survey. Source extinction shows evidence for dust grains that are larger than those in the diffuse interstellar medium (ISM): we estimate an average R_V=5.5 in the region. Importantly, we find a clear correlation between HRD inferred ages and spectroscopic surface-gravity inferred ages. This clearly indicates a real spread of stellar radii at fixed temperature, and together with additional correlations with extinction and with disk presence, strongly suggests a real spread of ages large than a few Myr. Focussing on the young population around NGC1980 iota Ori, which has previously been suggested to be a separate, foreground, older cluster, we confirm its older (5Myr) age and low A_V, but considering that its radial velocity distribution is indistinguishable from the Orion A's population, we suggest that NGC1980 is part of Orion A's star formation activity. Based on their stellar parameters and kinematic properties, we identify 383 new candidate members of Orion A, most of which are diskless sources in areas of the region poorly studied by previous works.

I review theoretical models of star formation and how they apply across the stellar mass spectrum. Several distinct theories are under active study for massive star formation, especially Turbulent Core Accretion, Competitive Accretion and Protostellar Mergers, leading to distinct observational predictions. These include the types of initial conditions, the structure of infall envelopes, disks and outflows, and the relation of massive star formation to star cluster formation. Even for Core Accretion models, there are several major uncertainties related to the timescale of collapse, the relative importance of different processes for preventing fragmentation in massive cores, and the nature of disks and outflows. I end by discussing some recent observational results that are helping to improve our understanding of these processes.

I review (1) Physics of Star Formation & Open Questions; (2) Structure & Dynamics of Star-Forming Clouds & Young Clusters; (3) Star Formation Rates: Observations & Theoretical Implications.

IRAS 22134+5834 was observed in the centimeter with (E)VLA, 3~mm with CARMA, 2~mm with PdBI, and 1.3~mm with SMA, to study the continuum emission as well as the molecular lines, that trace different physical conditions of the gas to study the influence of massive YSOs on nearby starless cores, and the possible implications in the clustered star formation process. The multi-wavelength centimeter continuum observations revealed two radio sources within the cluster, VLA1 and VLA2. VLA1 is considered to be an optically thin UCHII region with a size of 0.01~pc and sits at the edge of the near-infrared (NIR) cluster. The flux of ionizing photons of the VLA1 corresponds to a B1 ZAMS star. VLA2 is associated with an infrared point source and has a negative spectral index. We resolved six millimeter continuum cores at 2~mm, MM2 is associated with the UCHII region VLA1, and other dense cores are distributed around the UCH{\sc ii} region. Two high-mass starless clumps (HMSC), HMSC-E (east) and HMSC-W (west), are detected around the NIR cluster with N$_2$H$^+$(1--0) and NH$_3$ emission, and show different physical and chemical properties. Two N$_2$D$^+$ cores are detected on an NH$_3$ filament close to the UCHII region, with a projected separation of $\sim$8000~AU at the assumed distance of 2.6~kpc. The kinematic properties of the molecular line emission confirm the expansion of the UCHII region and that the molecular cloud around the near infrared (NIR) cluster is also expanding. Our multi-wavelength study has revealed different generations of star formation in IRAS 22134+5834. The formed intermediate- to massive stars show strong impact on nearby starless clumps. We propose that while the stellar wind from the UCHII region and the NIR cluster drives the large scale bubble, the starless clumps and HMPOs formed at the edge of the cluster.

The Kepler mission has discovered more than 4000 exoplanet candidates. Many are in systems with tightly packed inner planets. Inside-Out Planet Formation (IOPF) has been proposed to explain these systems. It involves sequential in situ planet formation at the local pressure maximum of a retreating dead zone inner boundary (DZIB). Pebbles accumulate at this pressure trap, which builds up a ring, and then a planet. The planet is expected to grow until it opens a gap, which helps to both truncate pebble accretion and induce DZIB retreat that sets the location of formation of the next planet. This simple scenario may be modified if the planet migrates significantly from its formation location. Thus planet-disk interactions play a crucial role in the IOPF scenario. We present numerical simulations that first assess migration of planets of various masses that are forming at the DZIB of an active accretion disk, where the effective viscosity rapidly increases in the radially inward direction. We find that the disk's torques on the planet tend to trap the planet at a location very close to the initial pressure maximum where it formed. We then study gap opening by these planets to assess at what mass a significant gap is created. Finally we present a simple model for DZIB retreat due to penetration of X-rays from the star to the disk midplane. Overall, these simulations help to quantify both the mass scale of first,"Vulcan," planet formation and the orbital separation to the location of second planet formation.

To understand massive star formation requires study of its initial conditions. Two massive starless core candidates, C1-N & C1-S, have been detected in IRDC G028.37+00.07 in $\rm N_2D^+$(3-2) with $ALMA$. From their line widths, either the cores are subvirial and are thus young structures on the verge of near free-fall collapse, or they are threaded by $\sim1$ mG $B$-fields that help support them in near virial equilibrium and potentially have older ages. We modeled the deuteration rate of $\rm N_2H^+$ to constrain collapse rates of the cores. First, to measure their current deuterium fraction, $D_{\rm frac}^{\rm N_2H^+}$ $\equiv [\rm N_2D^+]/[N_2H^+]$, we observed multiple transitions of $\rm N_2H^+$ and $\rm N_2D^+$ with $CARMA$, $SMA$, $JCMT$, $NRO~45m$ and $IRAM~30m$, to complement the $ALMA$ data. For both cores we derived $D_{\rm frac}^{\rm N_2H^+}\sim0.3$, several orders of magnitude above the cosmic [D]/[H] ratio. We then carried out chemodynamical modeling, exploring how collapse rate relative to free-fall, $\alpha_{\rm ff}$, affects the level of $D_{\rm frac}^{\rm N_2H^+}$ that is achieved from a given initial condition. To reach the observed $D_{\rm frac}^{\rm N_2H^+}$, most models require slow collapse with $\alpha_{\rm ff}\sim0.1$, i.e., $\sim1/10$th of free-fall. This makes it more likely that the cores have been able to reach a near virial equilibrium state and we predict that strong $B$-fields will eventually be detected. The methods developed here will be useful for measurement of the pre-stellar core mass function.

The Kepler-discovered Systems with Tightly-packed Inner Planets (STIPs), typically with several planets of Earth to super-Earth masses on well-aligned, sub-AU orbits may host the most common type of planets, including habitable planets, in the Galaxy. They pose a great challenge for planet formation theories, which fall into two broad classes: (1) formation further out followed by inward migration; (2) formation in situ, in the very inner regions of the protoplanetary disk. We review the pros and cons of these classes, before focusing on a new theory of sequential in situ formation from the inside-out via creation of successive gravitationally unstable rings fed from a continuous stream of small (~cm-m size) "pebbles," drifting inward via gas drag. Pebbles first collect at the pressure trap associated with the transition from a magnetorotational instability (MRI)-inactive ("dead zone") region to an inner, MRI-active zone. A pebble ring builds up that begins to dominate the local mass surface density of the disk and spawns a planet. The planet continues to grow, most likely by pebble accretion, until it becomes massive enough to isolate itself from the accretion flow via gap opening. This reduces the local gas density near the planet, leading to enhanced ionization and a retreat of the dead zone inner boundary. The process repeats with a new pebble ring gathering at the new pressure maximum associated with this boundary. We discuss the theory's predictions for planetary masses, relative mass scalings with orbital radius, and minimum orbital separations, and their comparison with observed systems. Finally, we discuss open questions, including potential causes of diversity of planetary system architectures, i.e., STIPs versus Solar System analogs.

We present spectroscopic tests of MIR to FIR extinction laws in IRDC G028.36+00.07, a potential site of massive star and star cluster formation. Lim & Tan (2014) developed methods of FIR extinction mapping of this source using ${\it Spitzer}$-MIPS ${\rm 24\mu m}$ and ${\it Herschel}$-PACS ${\rm 70\mu m}$ images, and by comparing to MIR ${\it Spitzer}$-IRAC $3$--${\rm 8\mu m}$ extinction maps, found tentative evidence for grain growth in the highest mass surface density regions. Here we present results of spectroscopic infrared extinction (SIREX) mapping using ${\it Spitzer}$-IRS (14 to ${\rm 38\mu m}$) data of the same IRDC. These methods allow us to first measure the SED of the diffuse Galactic ISM that is in the foreground of the IRDC. We then carry out our primary investigation of measuring the MIR to FIR opacity law and searching for potential variations as a function of mass surface density within the IRDC. We find relatively flat, featureless MIR-FIR opacity laws that lack the $\sim{\rm 12\mu m}$ and $\sim{\rm 35\mu m}$ features associated with the thick water ice mantle models of Ossenkopf & Henning (1994). Their thin ice mantle models and the coagulating aggregate dust models of Ormel et al. (2011) are a generally better match to the observed opacity laws. We also find evidence for generally flatter MIR to FIR extinction laws as mass surface density increases, strengthening the evidence for grain and ice mantle growth in higher density regions.

We utilize magnetohydrodynamic (MHD) simulations to develop a numerical model for GMC-GMC collisions between nearly magnetically critical clouds. The goal is to determine if, and under what circumstances, cloud collisions can cause pre-existing magnetically subcritical clumps to become supercritical and undergo gravitational collapse. We first develop and implement new photodissociation region (PDR) based heating and cooling functions that span the atomic to molecular transition, creating a multiphase ISM and allowing modeling of non-equilibrium temperature structures. Then in 2D and with ideal MHD, we explore a wide parameter space of magnetic field strength, magnetic field geometry, collision velocity, and impact parameter, and compare isolated versus colliding clouds. We find factors of ~2-3 increase in mean clump density from typical collisions, with strong dependence on collision velocity and magnetic field strength, but ultimately limited by flux-freezing in 2D geometries. For geometries enabling flow along magnetic field lines, greater degrees of collapse are seen. We discuss observational diagnostics of cloud collisions, focussing on 13CO(J=2-1), 13CO(J=3-2), and 12CO(J=8-7) integrated intensity maps and spectra, which we synthesize from our simulation outputs. We find the ratio of J=8-7 to lower-J emission is a powerful diagnostic probe of GMC collisions.

We study Giant Molecular Cloud (GMC) environments surrounding 10 Infrared Dark Clouds (IRDCs), using $^{13}$CO(1-0) emission from the Galactic Ring Survey. We measure physical properties of these IRDCs/GMCs on a range of scales extending to radii, R, of 30 pc. By comparing different methods for defining cloud boundaries and for deriving mass surface densities and velocity dispersions, we settle on a preferred "CE,$\tau$,G" method of "Connected Extraction" in position-velocity space plus Gaussian fitting to opacity-corrected line profiles for velocity dispersion and mass estimation. We examine how cloud definition affects measurements of the magnitude and direction of line-of-sight velocity gradients and velocity dispersions, including associated dependencies on size scale. CE,$\tau$,G-defined GMCs show velocity dispersion versus size relations $\sigma\propto{s}^{1/2}$, which are consistent with the large-scale gradients being caused by turbulence. However, IRDCs have velocity dispersions that are moderately enhanced above those predicted by this scaling relation. We examine the dynamical state of the clouds finding mean virial parameters $\bar{\alpha}_{\rm{vir}}\simeq 1.0$ for GMCs and 1.6 for IRDCs, broadly consistent with models of magnetized virialized pressure-confined polytropic clouds, but potentially indicating that IRDCs have more disturbed kinematics. CE,$\tau$,G-defined clouds exhibit a tight correlation of $\sigma/R^{1/2}\propto\Sigma^n$, with $n\simeq0.7$ for GMCs and 1.3 for IRDCs (c.f., a value of 0.5 expected for a population of virialized clouds). We conclude that while GMCs show evidence for virialization over a range of scales, IRDCs may be moderately super virial. Alternatively, IRDCs could be virialized but have systematically different $^{13}$CO gas phase abundances, i.e., due to freeze-out, affecting mass estimations.

The TMT Detailed Science Case describes the transformational science that the Thirty Meter Telescope will enable. Planned to begin science operations in 2024, TMT will open up opportunities for revolutionary discoveries in essentially every field of astronomy, astrophysics and cosmology, seeing much fainter objects much more clearly than existing telescopes. Per this capability, TMT's science agenda fills all of space and time, from nearby comets and asteroids, to exoplanets, to the most distant galaxies, and all the way back to the very first sources of light in the Universe. More than 150 astronomers from within the TMT partnership and beyond offered input in compiling the new 2015 Detailed Science Case. The contributing astronomers represent the entire TMT partnership, including the California Institute of Technology (Caltech), the Indian Institute of Astrophysics (IIA), the National Astronomical Observatories of the Chinese Academy of Sciences (NAOC), the National Astronomical Observatory of Japan (NAOJ), the University of California, the Association of Canadian Universities for Research in Astronomy (ACURA) and US associate partner, the Association of Universities for Research in Astronomy (AURA).

Most field stars will have encountered the highest stellar density and hence the largest number of interactions in their birth environment. Yet the stellar dynamics during this crucial phase are poorly understood. Here we analyze the radial velocities measured for 152 out of 380 observed stars in the 2-6 Myr old star cluster IC 348 as part of the SDSS-III APOGEE. The radial velocity distribution of these stars is fitted with one or two Gaussians, convolved with the measurement uncertainties including binary orbital motions. Including a second Gaussian improves the fit; the high-velocity outliers that are best fit by this second component may either (1) be contaminants from the nearby Perseus OB2 association, (2) be a halo of ejected or dispersing stars from IC 348, or (3) reflect that IC 348 has not relaxed to a Gaussian velocity distribution. We measure a velocity dispersion for IC 348 of $0.72 \pm 0.07$ km s$^{-1}$ (or $0.64 \pm 0.08$ km s$^{-1}$ if two Gaussians are fitted), which implies a supervirial state, unless the gas contributes more to the gravitational potential than expected. No evidence is found for a dependence of this velocity dispersion on distance from the cluster center or stellar mass. We also find that stars with lower extinction (in the front of the cloud) tend to be redshifted compared with stars with somewhat higher extinction (towards the back of the cloud). This data suggests that the stars in IC 348 are converging along the line of sight. We show that this correlation between radial velocity and extinction is unlikely to be spuriously caused by the small cluster rotation of $0.024 \pm 0.013$ km s$^{-1}$ arcmin$^{-1}$ or by correlations between the radial velocities of neighboring stars. This signature, if confirmed, will be the first detection of line-of-sight convergence in a star cluster(...)

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

Infrared dark clouds (IRDCs) are dense, molecular structures in the interstellar medium that can harbour sites of high-mass star formation. IRDCs contain supersonic turbulence, which is expected to generate shocks that locally heat pockets of gas within the clouds. We present observations of the CO J = 8-7, 9-8, and 10-9 transitions, taken with the Herschel Space Observatory, towards four dense, starless clumps within IRDCs (C1 in G028.37+00.07, F1 and F2 in G034.43+0007, and G2 in G034.77-0.55). We detect the CO J = 8-7 and 9-8 transitions towards three of the clumps (C1, F1, and F2) at intensity levels greater than expected from photodissociation region (PDR) models. The average ratio of the 8-7 to 9-8 lines is also found to be between 1.6 and 2.6 in the three clumps with detections, significantly smaller than expected from PDR models. These low line ratios and large line intensities strongly suggest that the C1, F1, and F2 clumps contain a hot gas component not accounted for by standard PDR models. Such a hot gas component could be generated by turbulence dissipating in low velocity shocks.

We present hydrodynamic simulations of self-gravitating dense gas in a galactic disk, exploring scales ranging from 1 kpc down to $\sim 0.1$~pc. Our primary goal is to understand how dense filaments form in Giant Molecular Clouds (GMCs). These structures, often observed as Infrared Dark Clouds (IRDCs) in the Galactic plane, are thought to be the precursors to massive stars and star clusters, so their formation may be the rate limiting step controlling global star formation rates in galactic systems as described by the Kennicutt-Schmidt relation. Our study follows on from Van Loo et al. (2013, Paper I), which carried out simulations to 0.5~pc resolution and examined global aspects of the formation of dense gas clumps and the resulting star formation rate. Here, using our higher resolution, we examine the detailed structural, kinematic and dynamical properties of dense filaments and clumps, including mass surface density ($\Sigma$) probability distribution functions, filament mass per unit length and its dispersion, lateral $\Sigma$ profiles, filament fragmentation, filament velocity gradients and infall, and degree of filament and clump virialization. Where possible, these properties are compared to observations of IRDCs. By many metrics, especially too large mass fractions of high $\Sigma>1\:{\rm g\:cm^{-2}}$ material, too high mass per unit length dispersion due to dense clump formation, too high velocity gradients and too high velocity dispersion for a given mass per unit length, the simulated filaments differ from observed IRDCs. We thus conclude that IRDCs do not form from global fast collapse of GMCs. Rather, we expect IRDC formation and collapse is slowed significantly by the influence of dynamically important magnetic fields, which may thus play a crucial role in regulating galactic star formation rates.

The deuterium fraction [N$_2$D$^+$]/[N$_2$H$^+$], may provide information about the ages of dense, cold gas structures, important to compare with dynamical models of cloud core formation and evolution. Here we introduce a complete chemical network with species containing up to three atoms, with the exception of the Oxygen chemistry, where reactions involving H$_3$O$^+$ and its deuterated forms have been added, significantly improving the consistency with comprehensive chemical networks. Deuterium chemistry and spin states of H$_2$ and H$_3^+$ isotopologues are included in this primarily gas-phase chemical model. We investigate dependence of deuterium chemistry on model parameters: density ($n_{\rm H}$), temperature, cosmic ray ionization rate, and gas-phase depletion factor of heavy elements ($f_{\rm D}$). We also explore the effects of time-dependent freeze-out of gas-phase species and dynamical evolution of density at various rates relative to free-fall collapse. For a broad range of model parameters, the timescales to reach large values of $D_{\rm frac}^{\rm N_2H^+} \gtrsim 0.1$, observed in some low- and high-mass starless cores, are relatively long compared to the local free-fall timescale. These conclusions are unaffected by introducing time-dependent freeze-out and considering models with evolving density, unless the initial $f_{\rm D} \gtrsim$ 10. For fiducial model parameters, achieving $D_{\rm frac}^{\rm N_2H^+} \gtrsim 0.1$ requires collapse to be proceeding at rates at least several times slower than that of free-fall collapse, perhaps indicating a dynamically important role for magnetic fields in the support of starless cores and thus the regulation of star formation.

Even today in our Galaxy, stars form from gas cores in a variety of environments, which may affect the properties of resulting star and planetary systems. Here we study the role of pressure, parameterized via ambient clump mass surface density, on protostellar evolution and appearance, focussing on low-mass, Sun-like stars and considering a range of conditions from relatively low pressure filaments in Taurus, to intermediate pressures of cluster-forming clumps like the Orion Nebula Cluster (ONC), to very high pressures that may be found in the densest Infrared Dark Clouds (IRDCs) or in the Galactic Center (GC). We present unified analytic and numerical models for collapse of prestellar cores, accretion disks, protostellar evolution and bipolar outflows, coupled to radiative transfer (RT) calculations and a simple astrochemical model to predict CO gas phase abundances. Prestellar cores in high pressure environments are smaller and denser and thus collapse with higher accretion rates and efficiencies, resulting in higher luminosity protostars with more powerful outflows. The protostellar envelope is heated to warmer temperatures, affecting infrared morphologies (and thus classification) and astrochemical processes like CO depletion to dust grain ice mantles (and thus CO morphologies). These results have general implications for star and planet formation, especially via their effect on astrochemical and dust grain evolution during infall to and through protostellar accretion disks.

The regulation of galactic-scale star formation rates (SFRs) is a basic problem for theories of galaxy formation and evolution: which processes are responsible for making observed star formation rates so inefficient compared to maximal rates of gas content divided by dynamical timescale? Here we study the effect of magnetic fields of different strengths on the evolution of giant molecular clouds (GMCs) within a kiloparsec patch of a disk galaxy and resolving scales down to $\simeq0.5\:{\rm{pc}}$. Including an empirically motivated prescription for star formation from dense gas ($n_{\rm{H}}>10^5\:{\rm{cm}^{-3}}$) at an efficiency of 2\% per local free-fall time, we derive the amount of suppression of star formation by magnetic fields compared to the nonmagnetized case. We find GMC fragmentation, dense clump formation and SFR can be significantly affected by the inclusion of magnetic fields, especially in our strongest investigated $B$-field case of $80\:{\rm{\mu}}$G. However, our chosen kpc-scale region, extracted from a global galaxy simulation, happens to contain a starbursting cloud complex that is only modestly affected by these magnetic fields and likely requires internal star formation feedback to regulate its SFR.

The compact multi-transiting systems discovered by Kepler challenge traditional planet formation theories. These fall into two broad classes: (1) formation further out followed by migration; (2) formation in situ from a disk of gas and planetesimals. In the former, an abundance of resonant chains is expected, which the Kepler data do not support. In the latter, required disk mass surface densities may be too high. A recently proposed mechanism hypothesizes that planets form in situ at the pressure trap associated with the dead-zone inner boundary (DZIB) where radially drifting "pebbles" accumulate. This scenario predicts planet masses ($M_p$) are set by the gap-opening process that then leads to DZIB retreat, followed by sequential, inside-out planet formation (IOPF). For typical disk accretion rates, IOPF predictions for $M_p$, $M_p$ versus orbital radius $r$, and planet-planet separations are consistent with observed systems. Here we investigate the IOPF prediction for how the masses, $M_{p,1}$, of the innermost ("Vulcan") planets vary with $r$. We show that for fiducial parameters, $M_{p,1}\simeq5.0(r/{\rm{0.1\:AU}})\:M_\oplus$, independent of the disk's accretion rate at time of planet formation. Then, using Monte Carlo sampling of a population of these innermost planets, we test this predicted scaling against observed planet properties, allowing for intrinsic dispersions in planetary densities and Kepler's observational biases. These effects lead to a slightly shallower relation $M_{p,1}\propto{r}^{0.9\pm0.2}$, which is consistent with $M_{p,1}\propto{r}^{0.7\pm0.2}$ of the observed Vulcans. The normalization of the relation constrains the gap-opening process, favoring relatively low viscosities in the inner dead zone.

The initial velocity dispersion of newborn stars is a major unconstrained aspect of star formation theory. Using near-infrared spectra obtained with the APOGEE spectrograph, we show that the velocity dispersion of young (1-2 Myr) stars in NGC 1333 is 0.92+/-0.12 km/s after correcting for measurement uncertainties and the effect of binaries. This velocity dispersion is consistent with the virial velocity of the region and the diffuse gas velocity dispersion, but significantly larger than the velocity dispersion of the dense, star-forming cores, which have a sub-virial velocity dispersion of 0.5 km/s. Since the NGC 1333 cluster is dynamically young and deeply embedded, this measurement provides a strong constraint on the initial velocity dispersion of newly-formed stars. We propose that the difference in velocity dispersion between stars and dense cores may be due to the influence of a 70 micro-Gauss magnetic field acting on the dense cores, or be the signature of a cluster with initial sub-structure undergoing global collapse.

High-mass Stars are cosmic engines known to dominate the energetics in the Milky Way and other galaxies. However, their formation is still not well understood. Massive, cold, dense clouds, often appearing as Infrared Dark Clouds (IRDCs), are the nurseries of massive stars. No measurements of magnetic fields in IRDCs in a state prior to the onset of high-mass star formation (HMSF) have previously been available, and prevailing HMSF theories do not consider strong magnetic fields. Here, we report observations of magnetic fields in two of the most massive IRDCs in the Milky Way. We show that IRDCs G11.11-0.12 and G0.253+0.016 are strongly magnetized and that the strong magnetic field is as important as turbulence and gravity for HMSF. The main dense filament in G11.11-0.12 is perpendicular to the magnetic field, while the lower density filament merging onto the main filament is parallel to the magnetic field. The implied magnetic field is strong enough to suppress fragmentation sufficiently to allow HMSF. Other mechanisms reducing fragmentation, such as the entrapment of heating from young stars via high mass surface densities, are not required to facilitate HMSF.

An ever growing number of observational and theoretical evidence suggests that the deuterated fraction (column density ratio between a species containing D and its hydrogenated counterpart, Dfrac) is an evolutionary indicator both in the low- and the high-mass star formation process. However, the role of surface chemistry in these studies has not been quantified from an observational point of view. In order to compare how the deuterated fractions of species formed only in the gas and partially or uniquely on grain surfaces evolve with time, we observed rotational transitions of CH3OH, 13CH3OH, CH2DOH, CH3OD at 3 and 1.3~mm, and of NH2D at 3~mm with the IRAM-30m telescope, and the inversion transitions (1,1) and (2,2) of NH3 with the GBT, towards most of the cores already observed by Fontani et al.~(2011, 2014) in N2H+, N2D+, HNC, DNC. NH2D is detected in all but two cores, regardless of the evolutionary stage. Dfrac(NH3) is on average above 0.1, and does not change significantly from the earliest to the most evolved phases, although the highest average value is found in the protostellar phase (~0.3). Few lines of CH2DOH and CH3OD are clearly detected, and only towards protostellar cores or externally heated starless cores. This work clearly confirms an expected different evolutionary trend of the species formed exclusively in the gas (N2D+ and N2H+) and those formed partially (NH2D and NH3) or totally (CH2DOH and CH3OH) on grain mantles. The study also reinforces the idea that Dfrac(N2H+) is the best tracer of massive starless cores, while high values of Dfrac(CH3OH) seem rather good tracers of the early protostellar phases, at which the evaporation/sputtering of the grain mantles is most efficient.

Inside-Out Planet Formation (IOPF; Chatterjee & Tan 2014, hereafter CT14) is a scenario for sequential in situ planet formation at the pressure traps of retreating dead zone inner boundaries (DZIBs) motivated to explain the many systems with tightly packed inner planets (STIPs) discovered by Kepler. The scenario involves build-up of a pebble-dominated protoplanetary ring, supplied by radial drift of pebbles from the outer disk. It may also involve further build-up of planetary masses to gap-opening scales via continued pebble accretion. Here we study radial drift & growth of pebbles delivered to the DZIB in fiducial IOPF disk models.

We study the degree of angular substructure in the stellar position distribution of young members of Galactic star-forming regions, looking for correlations with distance from cluster center, surface number density of stars, and local dynamical age. To this end we adopt the catalog of members in 18 young ($\sim$1-3 Myr) clusters from the Massive Young Star-Forming Complex Study in Infrared and X-ray (MYStIX) Survey and the statistical analysis of the Angular Dispersion Parameter, $\delta_{\rm ADP}$. We find statistically significant correlation between $\delta_{\rm ADP}$ and physical projected distance from the center of the clusters, with the centers appearing smoother than the outskirts, consistent with more rapid dynamical processing on local dynamical, free-fall or orbital timescales. Similarly, smoother distributions are seen in regions of higher surface density, or older dynamical ages. These results indicate that dynamical processing that erases substructure is already well-advanced in young, sometimes still-forming, clusters. Such observations of the dissipation of substructure have the potential to constrain theoretical models of the dynamical evolution of young and forming clusters.

The enormous radiative and mechanical luminosities of massive stars impact a vast range of scales and processes, from the reionization of the universe, to the evolution of galaxies, to the regulation of the interstellar medium, to the formation of star clusters, and even to the formation of planets around stars in such clusters. Two main classes of massive star formation theory are under active study, Core Accretion and Competitive Accretion. In Core Accretion, the initial conditions are self-gravitating, centrally concentrated cores that condense with a range of masses from the surrounding, fragmenting clump environment. They then undergo relatively ordered collapse via a central disk to form a single star or a small-N multiple. In this case, the pre-stellar core mass function has a similar form to the stellar initial mass function. In Competitive Accretion, the material that forms a massive star is drawn more chaotically from a wider region of the clump without passing through a phase of being in a massive, coherent core. In this case, massive star formation must proceed hand in hand with star cluster formation. If stellar densities become very high near the cluster center, then collisions between stars may also help to form the most massive stars. We review recent theoretical and observational progress towards understanding massive star formation, considering physical and chemical processes, comparisons with low and intermediate-mass stars, and connections to star cluster formation.

The spatial morphology and dynamical status of a young, still-forming stellar cluster provide valuable clues on the conditions during the star formation event and the processes that regulated it. We analyze the Orion Nebula Cluster (ONC), utilizing the latest censuses of its stellar content and membership estimates over a large wavelength range. We determine the center of mass of the ONC, and study the radial dependence of angular substructure. The core appears rounder and smoother than the outskirts, consistent with a higher degree of dynamical processing. At larger distances the departure from circular symmetry is mostly driven by the elongation of the system, with very little additional substructure, indicating a somewhat evolved spatial morphology or an expanding halo. We determine the mass density profile of the cluster, which is well fitted by a power law that is slightly steeper than a singular isothermal sphere. Together with the ISM density, estimated from average stellar extinction, the mass content of the ONC is insufficient by a factor $\sim 1.8$ to reproduce the observed velocity dispersion from virialized motions, in agreement with previous assessments that the ONC is moderately supervirial. This may indicate recent gas dispersal. Based on the latest estimates for the age spread in the system and our density profiles, we find that, at the half-mass radius, 90% of the stellar population formed within $\sim 5$-$8$ free-fall times ($t_{\rm ff}$). This implies a star formation efficiency per $t_{\rm ff}$ of $\epsilon_{\rm ff}\sim 0.04$-$0.07$, i.e., relatively slow and inefficient star formation rates during star cluster formation.

Over two years 8,859 high-resolution H-band spectra of 3,493 young (1 - 10 Myr) stars were gathered by the multi-object spectrograph of the APOGEE project as part of the IN-SYNC ancillary program of that SDSS-III survey. Here we present the forward modeling approach used to derive effective temperatures, surface gravities, radial velocities, rotational velocities, and H-band veiling from these near-infrared spectra. We discuss in detail the statistical and systematic uncertainties in these stellar parameters. In addition we present accurate extinctions by measuring the E(J-H) of these young stars with respect to the single-star photometric locus in the Pleiades. Finally we identify an intrinsic stellar radius spread of about 25% for late-type stars in IC 348 using three (nearly) independent measures of stellar radius, namely the extinction-corrected J-band magnitude, the surface gravity and the $R \sin i$ from the rotational velocities and literature rotation periods. We exclude that this spread is caused by uncertainties in the stellar parameters by showing that the three estimators of stellar radius are correlated, so that brighter stars tend to have lower surface gravities and larger $R \sin i$ than fainter stars at the same effective temperature. Tables providing the spectral and photometric parameters for the Pleiades and IC 348 have been provided online.

IRIS (InfraRed Imaging Spectrograph) is a first light near-infrared diffraction limited imager and integral field spectrograph being designed for the future Thirty Meter Telescope (TMT). IRIS is optimized to perform astronomical studies across a significant fraction of cosmic time, from our Solar System to distant newly formed galaxies (Barton et al. [1]). We present a selection of the innovative science cases that are unique to IRIS in the era of upcoming space and ground-based telescopes. We focus on integral field spectroscopy of directly imaged exoplanet atmospheres, probing fundamental physics in the Galactic Center, measuring 10^4 to 10^10 Msun supermassive black hole masses, resolved spectroscopy of young star-forming galaxies (1 < z < 5) and first light galaxies (6 < z < 12), and resolved spectroscopy of strong gravitational lensed sources to measure dark matter substructure. For each of these science cases we use the IRIS simulator (Wright et al. [2], Do et al. [3]) to explore IRIS capabilities. To highlight the unique IRIS capabilities, we also update the point and resolved source sensitivities for the integral field spectrograph (IFS) in all five broadband filters (Z, Y, J, H, K) for the finest spatial scale of 0.004" per spaxel. We briefly discuss future development plans for the data reduction pipeline and quicklook software for the IRIS instrument suite.

We present radiation transfer (RT) simulations of evolutionary sequences of massive protostars forming from massive dense cores in environments of high surface densities. The protostellar evolution is calculated with a detailed multi-zone model, with the accretion rate regulated by feedback from an evolving disk-wind outflow cavity. Disk and envelope evolutions are calculated self-consistently. In this framework, an evolutionary track is determined by three environmental initial conditions: the initial core mass M_c, the mean surface density of the ambient star-forming clump Sigma_cl, and the rotational-to-gravitational energy ratio of the initial core, beta_c. Evolutionary sequences with various M_c, Sigma_cl, beta_c are constructed. We find that in a fiducial model with M_c=60Msun, Sigma_cl=1 g/cm^2 and beta_c=0.02, the final star formation efficiency >~0.43. For each evolutionary track, RT simulations are performed at selected stages, with temperature profiles, SEDs, and images produced. At a given stage the envelope temperature is highly dependent on Sigma_cl, but only weakly dependent on M_c. The SED and MIR images depend sensitively on the evolving outflow cavity, which gradually wides as the protostar grows. The fluxes at <~100 microns increase dramatically, and the far-IR peaks move to shorter wavelengths. We find that, despite scatter caused by different M_c, Sigma_cl, beta, and inclinations, sources at a given evolutionary stage appear in similar regions on color-color diagrams, especially when using colors at >~ 70 microns, where the scatter due to the inclination is minimized, implying that such diagrams can be useful diagnostic tools of evolutionary stages of massive protostars. We discuss how intensity profiles along or perpendicular to the outflow axis are affected by environmental conditions and source evolution.

We use observed radial profiles of mass surface densities of total, $\Sigma_g$, & molecular, $\Sigma_{\rm H2}$, gas, rotation velocity & star formation rate (SFR) surface density, $\Sigma_{\rm sfr}$, of the molecular-rich ($\Sigma_{\rm H2}\ge\Sigma_{\rm HI}/2$) regions of 16 nearby disk galaxies to test several star formation laws: a Kennicutt-Schmidt law, $\Sigma_{\rm sfr}=A_g\Sigma_{g,2}^{1.5}$; a Constant Molecular law, $\Sigma_{\rm sfr}=A_{\rm H2}\Sigma_{\rm H2,2}$; the turbulence-regulated laws of Krumholz & McKee (KM05) and Krumholz et al. (KMT09), a Gas-$\Omega$ law, $\Sigma_{\rm sfr}=B_\Omega\Sigma_g\Omega$; and a shear-driven GMC Collision law, $\Sigma_{\rm sfr}=B_{\rm CC}\Sigma_g\Omega(1-0.7\beta)$, where $\beta\equiv d {\rm ln} v_{\rm circ}/d {\rm ln} r$. If allowed one free normalization parameter for each galaxy, these laws predict the SFR with rms errors of factors of 1.4 - 1.8. If a single normalization parameter is used by each law for the entire galaxy sample, then rms errors range from factors of 1.5 - 2.1. Although the Constant Molecular law gives the smallest errors, the improvement over KMT, Kennicutt-Schmidt & GMC Collision laws is not especially significant, particularly given the different observational inputs that the laws utilize and the scope of included physics, which ranges from empirical relations to detailed treatment of interstellar medium processes. We next search for variation of star formation law parameters with local & global galactic dynamical properties of disk shear rate (related to $\beta$), rotation speed & presence of a bar. We demonstrate with high significance that higher shear rates enhance star formation efficiency per local orbital time. Such a trend is expected if GMC collisions play an important role in star formation, while an opposite trend would be expected if development of disk gravitational instabilities is the controlling physics.

Infrared Dark Clouds (IRDCs) are unique laboratories to study the initial conditions of high-mass star and star cluster formation. We present high-sensitivity and high-angular resolution IRAM PdBI observations of N2H+ (1-0) towards IRDC G035.39-00.33. It is found that G035.39-00.33 is a highly complex environment, consisting of several mildly supersonic filaments (sigma_NT/c_s ~1.5), separated in velocity by <1 km s^-1 . Where multiple spectral components are evident, moment analysis overestimates the non-thermal contribution to the line-width by a factor ~2. Large-scale velocity gradients evident in previous single-dish maps may be explained by the presence of substructure now evident in the interferometric maps. Whilst global velocity gradients are small (<0.7 km s^-1 pc^-1), there is evidence for dynamic processes on local scales (~1.5-2.5 km s^-1 pc^-1 ). Systematic trends in velocity gradient are observed towards several continuum peaks. This suggests that the kinematics are influenced by dense (and in some cases, starless) cores. These trends are interpreted as either infalling material, with accretion rates ~(7 \pm 4)x10^-5 M_sun yr^-1 , or expanding shells with momentum ~24 \pm 12 M_sun km s^-1 . These observations highlight the importance of high-sensitivity and high-spectral resolution data in disentangling the complex kinematic and physical structure of massive star forming regions.

We use deep $8\:\mu m$ Spitzer-IRAC imaging of a massive Infrared Dark Cloud (IRDC) G028.37+00.07 to construct a Mid-Infrared (MIR) extinction map that probes mass surface densities up to $\Sigma\:\sim 1\:\rm{g~cm^{-2}}$ ($A_V\sim200\:$mag), amongst the highest values yet probed by extinction mapping. Merging with a NIR extinction map of the region, creates a high dynamic range map that reveals structures down to $A_V\sim1\:$mag. We utilize the map to: (1) Measure a cloud mass $\sim7\times10^4\:M_\odot$ within a radius of $\sim8\:$pc. $^{13}$CO kinematics indicate that the cloud is gravitationally bound. It thus has the potential to form one of the most massive young star clusters known in the Galaxy. (2) Characterize the structures of 16 massive cores within the IRDC, finding they can be fit by singular polytropic spheres with $\rho\propto{r}^{-k_\rho}$ and $k_\rho=1.3\pm0.3$. They have $\overline{\Sigma}\simeq0.1-0.4\:\rm{g~cm^{-2}}$ --- relatively low values that, along with their measured cold temperatures, suggest magnetic fields, rather than accretion-powered radiative heating, are important for controlling fragmentation of these cores. (3) Determine the $\Sigma$ (equivalently column density or $A_V$) probability distribution function (PDF) for a region that is near complete for $A_V>3\:$mag. The PDF is well fit by a single log-normal with mean $\overline{A}_V\simeq9\:$mag, high compared to other known clouds. It does not exhibit a separate high-end power law tail, which has been claimed to indicate the importance of self-gravity. However, we suggest that the PDF does result from a self-similar, self-gravitating hierarchy of structure being present over a wide range of scales in the cloud.

Some theories of dense molecular cloud formation involve dynamical environments driven by converging atomic flows or collisions between preexisting molecular clouds. The determination of the dynamics and physical conditions of the gas in clouds at the early stages of their evolution is essential to establish the dynamical imprints of such collisions, and to infer the processes involved in their formation. We present multi-transition 13CO and C18O maps toward the IRDC G035.39-00.33, believed to be at the earliest stages of evolution. The 13CO and C18O gas is distributed in three filaments (Filaments 1, 2 and 3), where the most massive cores are preferentially found at the intersecting regions between them. The filaments have a similar kinematic structure with smooth velocity gradients of ~0.4-0.8 km s-1 pc-1. Several scenarios are proposed to explain these gradients, including cloud rotation, gas accretion along the filaments, global gravitational collapse, and unresolved sub-filament structures. These results are complemented by HCO+, HNC, H13CO+ and HN13C single-pointing data to search for gas infall signatures. The 13CO and C18O gas motions are supersonic across G035.39-00.33, with the emission showing broader linewidths toward the edges of the IRDC. This could be due to energy dissipation at the densest regions in the cloud. The average H2 densities are ~5000-7000 cm-3, with Filaments 2 and 3 being denser and more massive than Filament 1. The C18O data unveils three regions with high CO depletion factors (f_D~5-12), similar to those found in massive starless cores.

Progress in understanding star formation requires detailed observational constraints on the initial conditions, i.e. dense clumps and cores in giant molecular clouds that are on the verge of gravitational instability. Such structures have been studied by their extinction of Near-Infrared (NIR) and, more recently, Mid-Infrared (MIR) background light. It has been somewhat more of a surprise to find that there are regions that appear as dark shadows at Far-Infrared (FIR) wavelengths as long as $\sim$100$\mu m$. Here we develop analysis methods of FIR images from Spitzer-MIPS and Herschel-PACS that allow quantitative measurements of cloud mass surface density, $\Sigma$. The method builds upon that developed for MIR extinction mapping (MIREX) (Butler and Tan 2012), in particular involving a search for independent saturated, i.e. very opaque, regions that allow measurement of the foreground intensity. We focus on three massive starless core/clumps in IRDC G028.37+00.07, deriving mass surface density maps from 3.5 to 70$\mu m$. A by-product of this analysis is measurement of the spectral energy distribution of the diffuse foreground emission. The lower opacity at 70$\mu m$ allows us to probe to higher $\Sigma$ values, up to $\sim1\:{\rm{g\:cm}^{-2}}$ in the densest parts of the core/clumps. Comparison of the $\Sigma$ maps at different wavelengths constrains the shape of the MIR-FIR dust opacity law in IRDCs. We find it is most consistent with the thick ice mantle models of Ossenkopf and Henning (1994). There is tentative evidence for grain ice mantle growth as one goes from lower to higher $\Sigma$ regions.

How do stars that are more massive than the Sun form, and thus how is the stellar initial mass function (IMF) established? Such intermediate- and high-mass stars may be born from relatively massive pre-stellar gas cores, which are more massive than the thermal Jeans mass. The Turbulent Core Accretion model invokes such cores as being in approximate virial equilibrium and in approximate pressure equilibrium with their surrounding clump medium. Their internal pressure is provided by a combination of turbulence and magnetic fields. Alternatively, the Competitive Accretion model requires strongly sub-virial initial conditions that then lead to extensive fragmentation to the thermal Jeans scale, with intermediate- and high-mass stars later forming by competitive Bondi-Hoyle accretion. To test these models, we have identified four prime examples of massive (~100Msun) clumps from mid-infrared extinction mapping of infrared dark clouds (IRDCs). Fontani et al. found high deuteration fractions of N2H+ in these objects, which are consistent with them being starless. Here we present ALMA observations of these four clumps that probe the N2D+(3-2) line at 2.3" resolution. We find six N2D+ cores and determine their dynamical state. Their observed velocity dispersions and sizes are broadly consistent with the predictions of the Turbulent Core model of self-gravitating, magnetized (with Alfven Mach number m_A~1) and virialized cores that are bounded by the high pressures of their surrounding clumps. However, in the most massive cores, with masses up to ~60Msun, our results suggest that moderately enhanced magnetic fields (so that m_A~0.3) may be needed for the structures to be in virial and pressure equilibrium. Magnetically regulated core formation may thus be important in controlling the formation of massive cores, inhibiting their fragmentation, and thus helping to establish the stellar IMF.

The compact multi-transiting planet systems discovered by Kepler challenge planet formation theories. Formation in situ from disks with radial mass surface density, $\Sigma$, profiles similar to the minimum mass solar nebula (MMSN) but boosted in normalization by factors $\gtrsim 10$ has been suggested. We propose that a more natural way to create these planets in the inner disk is formation sequentially from the inside-out via creation of successive gravitationally unstable rings fed from a continuous stream of small (~cm--m size) "pebbles", drifting inwards via gas drag. Pebbles collect at the pressure maximum associated with the transition from a magneto-rotational instability (MRI)-inactive ("dead zone") region to an inner MRI-active zone. A pebble ring builds up until it either becomes gravitationally unstable to form an $\sim 1\ M_\oplus$ planet directly or induces gradual planet formation via core accretion. The planet may undergo Type I migration into the active region, allowing a new pebble ring and planet to form behind it. Alternatively if migration is inefficient, the planet may continue to accrete from the disk until it becomes massive enough to isolate itself from the accretion flow. A variety of densities may result depending on the relative importance of residual gas accretion as the planet approaches its isolation mass. The process can repeat with a new pebble ring gathering at the new pressure maximum associated with the retreating dead zone boundary. Our simple analytical model for this scenario of inside-out planet formation yields planetary masses, relative mass scalings with orbital radius, and minimum orbital separations consistent with those seen by Kepler. It provides an explanation of how massive planets can form with tightly-packed and well-aligned system architectures, starting from typical protoplanetary disk properties.

(Abridged) The Census of High- and Medium-mass Protostars (CHaMP) is the first large-scale (280 degree<l<300 degree, -4 degree<b<2 degree), unbiased, sub-parsec resolution survey of Galactic molecular clumps and their embedded stars. Barnes et al. (2011) presented the source catalog of ~300 clumps based on HCO+(1-0) emission, used to estimate masses M. Here we use archival mid-infrared to mm continuum data to construct spectral energy distributions. Fitting two-temperature grey-body models, we derive bolometric luminosities, L. We find the clumps have 10Lsun<L<1E6.5Lsun and 0.1<L/M<1E3, consistent with theoretical expectations of a clump population that spans a range of instantaneous star formation efficiencies from 0 to ~50%. We thus expect L/M to be a useful, strongly-varying indicator of clump evolution during the star cluster formation process. We find correlations of the ratio of warm to cold component fluxes and of cold component temperature with L/M. We also find a near linear relation between L/M and Spitzer-IRAC specific intensity (surface brightness), which may thus also be useful as a star formation efficiency indicator. The lower bound of the clump $L/M$ distribution suggests the star formation efficiency per free-fall time is epsilon<0.2. We do not find strong correlations of L/M with mass surface density, velocity dispersion or virial parameter. We find a linear relation between L and L_{HCO+(1-0}}, although with large scatter for any given individual clump. Fitting together with extragalactic systems, the linear relation still holds, extending over 10 orders of magnitude in luminosity. The complete nature of the CHaMP survey over a several kiloparsec-scale region allows us to derive a measurement at an intermediate scale bridging those of individual clumps and whole galaxies.

We report new near-infrared and mm-wave observational data on a selection of massive Galactic molecular clumps (part of the CHaMP sample) and their associated young star clusters. The clumps show, for the first time in a "dense gas tracer", a significant correlation between HCO+ line emission from cold molecular gas and Br{\gamma} line emission of associated nebulae. This correlation arises in the HCO+ line's brightness, not its linewidth. In contrast, the correlation between the N2H+ line emission and Br{\gamma} is weak or absent. The HCO+/N2H+ line ratio also varies widely from clump to clump: bright HCO+ emission tends to be more closely associated with Br{\gamma} nebulosity, while bright N2H+ emission tends to avoid areas that are bright in Br{\gamma}. Both molecular species show correlations of weak significance with infrared H2 v=1-0 and v=2-1 line emission, in or near the clumps. The H2 emission line ratio is consistent with fluorescent excitation in most of the clumps, although thermal excitation is seen in a few clumps. We interpret these trends as evidence for evolution in the gas conditions due to the effects of ongoing star formation in the clumps, in particular, the importance of UV radiation from massive YSOs as the driving agent that heats the molecular gas and alters its chemistry. This suggests that some traditional dense gas tracers of molecular clouds do not sample a homogeneous population of clumps, i.e., that the HCO+ brightness in particular is directly related to the heating and disruption of cold gas by massive young stars, whereas the N2H+ better samples gas not yet affected by this process. We therefore suggest that the HCO+-N2H+-Br{\gamma} relationship is a useful diagnostic of a molecular clump's progress in forming massive stars.

We aim at characterising dense cores in the clustered environments associated with massive star-forming regions. For this, we present an uniform analysis of VLA NH3(1,1) and (2,2) observations towards a sample of 15 massive star-forming regions, where we identify a total of 73 cores, classify them as protostellar, quiescent starless, or perturbed starless, and derive some physical properties. The average sizes and ammonia column densities are 0.06 pc and 10^15 cm^-2, respectively, with no significant differences between the starless and protostellar cores, while the linewidth and rotational temperature of quiescent starless cores are smaller, 1.0 km/s and 16 K, than those of protostellar (1.8 km/s, 21 K), and perturbed starless (1.4 km/s, 19 K) cores. Such linewidths and temperatures for these quiescent starless cores in the surroundings of massive stars are still significantly larger than the typical values measured in starless cores of low-mass star-forming regions, implying an important non-thermal component. We confirm at high angular resolutions the correlations previously found with single-dish telescopes between the linewidth, the temperature of the cores, and the bolometric luminosity. In addition, we find a correlation between the temperature of each core and the incident flux from the most massive star in the cluster, suggesting that the large temperatures measured in the starless cores of our sample could be due to heating from the nearby massive star. A simple virial equilibrium analysis seems to suggest a scenario of a self-similar, self-graviting, turbulent, virialised hierarchy of structures from clumps (0.1-10 pc) to cores (0.05 pc). A closer inspection of the dynamical state taking into account external pressure effects, reveal that relatively strong magnetic field support may be needed to stabilise the cores, or that they are unstable and thus on the verge of collapse.

We present 30 and 40 micron imaging of the massive protostar G35.20-0.74 with SOFIA-FORCAST. The high surface density of the natal core around the protostar leads to high extinction, even at these relatively long wavelengths, causing the observed flux to be dominated by that emerging from the near-facing outflow cavity. However, emission from the far-facing cavity is still clearly detected. We combine these results with fluxes from the near-infrared to mm to construct a spectral energy distribution (SED). For isotropic emission the bolometric luminosity would be 3.3x10^4 Lsun. We perform radiative transfer modeling of a protostar forming by ordered, symmetric collapse from a massive core bounded by a clump with high mass surface density, Sigma_cl. To fit the SED requires protostellar masses ~20-34 Msun depending on the outflow cavity opening angle (35 - 50 degrees), and Sigma_cl ~ 0.4-1 g cm-2. After accounting for the foreground extinction and the flashlight effect, the true bolometric luminosity is ~ (0.7-2.2)x10^5 Lsun. One of these models also has excellent agreement with the observed intensity profiles along the outflow axis at 10, 18, 31 and 37 microns. Overall our results support a model of massive star formation involving the relatively ordered, symmetric collapse of a massive, dense core and the launching bipolar outflows that clear low density cavities. Thus a unified model may apply for the formation of both low and high mass stars.

(Abridged) We present radiation transfer simulations of a massive (8 Msun) protostar forming from a massive (Mc=60 Msun) protostellar core, extending the model developed by Zhang & Tan (2011). The two principal improvements are (1) developing a model for the density and velocity structure of a disk wind that fills the bipolar outflow cavities; and (2) solving for the radially varying accretion rate in the disk due to a supply of mass and angular momentum from the infall envelope and their loss to the disk wind. One consequence of the launching of the disk wind is a reduction in the amount of accretion power that is radiated by the disk. For the transition from dusty to dust-free conditions where gas opacities dominate, we now implement a gradual change as a more realistic approximation of dust destruction. We study how the above effects, especially the outflow, influence the SEDs and the images of the protostar. Dust in the outflow cavity significantly affects the SEDs at most viewing angles. It further attenuates the short-wavelength flux from the protostar, controlling how the accretion disk may be viewed, and contributes a significant part of the near- and mid-IR fluxes. These fluxes warm the disk, boosting the mid- and far-IR emission. We find that for near face-on views, the SED from the near-IR to about 60 micron is very flat, which may be used to identify such systems. We show that the near-facing outflow cavity and its walls are still the most significant features in images up to 70 micron, dominating the mid-IR emission and determining its morphology. The thermal emission from the dusty outflow itself dominates the flux at ~20 micron. The detailed distribution of the dust in the outflow affects the morphology, for example, even though the outflow cavity is wide, at 10 to 20 micron, the dust in the disk wind can make the outflow appear narrower than in the near-IR bands.

All stars are born in molecular clouds, and most in giant molecular clouds (GMCs), which thus set the star formation activity of galaxies. We first review their observed properties, including measures of mass surface density, Sigma, and thus mass, M. We discuss cloud dynamics, concluding most GMCs are gravitationally bound. Star formation is highly clustered within GMCs, but overall is very inefficient. We compare properties of star-forming clumps with those of young stellar clusters (YSCs). The high central densities of YSCs may result via dynamical evolution of already-formed stars during and after star cluster formation. We discuss theoretical models of GMC evolution, especially addressing how turbulence is maintained, and emphasizing the importance of GMC collisions. We describe how feedback limits total star formation efficiency, epsilon, in clumps. A turbulent and clumpy medium allows higher epsilon, permitting formation of bound clusters even when escape speeds are less than the ionized gas sound speed.

We present hydrodynamic simulations of the evolution of self-gravitating dense gas on scales of 1 kiloparsec down to < parsec in a galactic disk, designed to study dense clump formation from giant molecular clouds (GMCs). These structures are expected to be the precursors to star clusters and this process may be the rate limiting step controling star formation rates in galactic systems as described by the Kennicutt-Schmidt relation. We follow the thermal evolution of the gas down to ~5K using extinction-dependent heating and cooling functions. We do not yet include magnetic fields or localized stellar feedback, so the evolution of the GMCs and clumps is determined solely by self-gravity balanced by thermal and turbulent pressure support and the large scale galactic shear. While cloud structures and densities change significantly during the simulation, GMC virial parameters remain mostly above unity for time scales exceeding the free-fall time of GMCs indicating that energy from galactic shear and large-scale cloud motions continuously cascades down to and within the GMCs. We implement star formation at a slow, inefficient rate of 2% per local free-fall time, but even this yields global star formation rates that are about two orders of magnitude larger than the observed Kennicutt-Schmidt relation due to over-production of dense gas clumps. We expect a combination of magnetic support and localized stellar feedback is required to inhibit dense clump formation to ~1% of the rate that results from the nonmagnetic, zero-feedback limit.

We present simulations of the evolution of self-gravitating dense gas on kiloparsec-size scales in a galactic disk, designed to study dense clump formation from giant molecular clouds (GMCs). These dense clumps are expected to be the precursors to star clusters and this process may be the rate limiting step controling star formation rates in galactic systems as described by the Kennicutt-Schmidt relation. The evolution of these simulated GMCs and clumps is determined by self-gravity balanced by turbulent pressure support and the large scale galactic shear. While the cloud structures and densities significantly change during their evolution, they remain roughly in virial equilibrium for time scales exceeding the free-fall time of GMCs, indicating that energy from the galactic shear continuously cascades down. We implement star formation at a slow, inefficient rate of 2% per local free-fall time, but this yields global star formation rates that are more than ~two orders of magnitude larger than the observed Kennicutt-Schmidt relation due to the over-production of dense clump gas. To explain this discrepancy, we anticipate magnetic fields to provide additional support. Low-resolution simulations indeed show that the magnetic field reduces the star formation rate.

Measuring the mass distribution of infrared dark clouds (IRDCs) over the wide dynamic range of their column densities is a fundamental obstacle in determining the initial conditions of high-mass star formation and star cluster formation. We present a new technique to derive high-dynamic-range, arcsecond-scale resolution column density data for IRDCs and demonstrate the potential of such data in measuring the density variance - sonic Mach number relation in molecular clouds. We combine near-infrared data from the UKIDSS/Galactic Plane Survey with mid-infrared data from the Spitzer/GLIMPSE survey to derive dust extinction maps for a sample of ten IRDCs. We then examine the linewidths of the IRDCs using 13CO line emission data from the FCRAO/Galactic Ring Survey and derive a column density - sonic Mach number relation for them. For comparison, we also examine the relation in a sample of nearby molecular clouds. The presented column density mapping technique provides a very capable, temperature independent tool for mapping IRDCs over the column density range equivalent to A_V=1-100 mag at a resolution of 2". Using the data provided by the technique, we present the first direct measurement of the relationship between the column density dispersion, \sigma_{N/<N>}, and sonic Mach number, M_s, in molecular clouds. We detect correlation between the variables with about 3-sigma confidence. We derive the relation \sigma_{N/<N>} = (0.047 \pm 0.016) Ms, which is suggestive of the correlation coefficient between the volume density and sonic Mach number, \sigma_{\rho/<\rho>} = (0.20^{+0.37}_{-0.22}) Ms, in which the quoted uncertainties indicate the 3-sigma range. When coupled with the results of recent numerical works, the existence of the correlation supports the picture of weak correlation between the magnetic field strength and density in molecular clouds (i.e., B ~ \rho^{0.5}).

Infrared Dark Clouds (IRDCs) host the initial conditions under which massive stars and stellar clusters form. We have obtained high sensitivity and high spectral resolution observations with the IRAM 30m antenna, which allowed us to perform detailed analysis of the kinematics within one IRDC, G035.39-00.33. We focus on the 1-0 and 3-2 transitions of N2H+, C18O (1-0), and make comparison with SiO (2-1) observations and extinction mapping. Three interacting filaments of gas are found. We report large-scale velocity coherence throughout the cloud, evidenced through small velocity gradients and relatively narrow line widths. This suggests that the merging of these filaments is somewhat "gentle", possibly regulated by magnetic fields. This merging of filaments may be responsible for the weak parsec-scale SiO emission detected by Jimenez-Serra et al. 2010, via grain mantle vaporization. A systematic velocity shift between the N2H+ (1-0) and C18O (1-0) gas throughout the cloud of 0.18 +/- 0.04 kms^{-1} is also found, consistent with a scenario of collisions between filaments which is still ongoing. The N2H+ (1-0) is extended throughout the IRDC and it does not only trace dense cores, as found in nearby low-mass star-forming regions. The average H2 number density across the IRDC is ~ 5 x 10^4 cm^{-3}, at least one order of magnitude larger than in nearby molecular clouds where low-mass stars are forming. A temperature gradient perpendicular to the filament is found. From our study, we conclude that G035.39-00.33 (clearly seen in the extinction map and in N2H+) has been formed via the collision between two relatively quiescent filaments with average densities of ~ 5 x 10^3 cm^{-3}, moving with relative velocities of ~ 5 kms^{-1}. The accumulation of material at the merging points started > 1 Myr ago and it is still ongoing.

The initial conditions of massive star and star cluster formation are expected to be cold, dense and high column density regions of the interstellar medium, which can reveal themselves via near, mid and even far-infrared absorption as Infrared Dark Clouds (IRDCs). Elucidating the dynamical state of IRDCs thus constrains theoretical models of these complex processes. In particular, it is important to assess whether IRDCs have reached virial equilibrium, where the internal pressure balances that due to the self-gravitating weight of the cloud plus the pressure of the external environmental. We study this question for the filamentary IRDC G035.39-00.33 by deriving mass from combined NIR & MIR extinction maps and velocity dispersion from C18O (1-0) & (2-1) line emission. In contrast to our previous moderately super-virial results based on 13CO emission and MIR-only extinction mapping, with improved mass measurements we now find that the filament is consistent with being in virial equilibrium, at least in its central parsec-wide region where ~1000 M_Sun snakes along several parsecs. This equilibrium state does not require large-scale net support or confinement by magnetic fields.

The Orion Nebula Cluster (ONC) is the nearest region of massive star formation and thus a crucial testing ground for theoretical models. Of particular interest amongst the ONC's ~1000 members are: \theta^1 Ori C, the most massive binary in the cluster with stars of masses 38 and 9 MSun (Kraus et al. 2009); the Becklin-Neugebauer (BN) object, a 30 km/s runaway star of ~8 MSun (Tan 2004); and the Kleinmann-Low (KL) nebula protostar, a highly-obscured, ~15 MSun object still accreting gas while also driving a powerful, apparently "explosive" outflow (Allen & Burton 1993). The unusual behavior of BN and KL is much debated: How did BN acquire its high velocity? How is this related to massive star formation in the KL nebula? Here we report the results of a systematic survey using ~ 10^7 numerical experiments of gravitational interactions of the \theta^1C and BN stars. We show that dynamical ejection of BN from this triple system at its observed velocity leaves behind a binary with total energy and eccentricity matching those observed for \theta^1C. Five other observed properties of \theta^C are also consistent with it having ejected BN and altogether we estimate there is only a <~ 10^{-5} probability that \theta^1C has these properties by chance. We conclude that BN was dynamically ejected from the \theta^1C system about 4,500 years ago. BN has then plowed through the KL massive-star-forming core within the last 1,000 years causing its recently-enhanced accretion and outflow activity.

(abridged) We develop the mid-infrared extinction (MIREX) mapping technique of Butler & Tan (2009, Paper I), presenting a new method to correct for the Galactic foreground emission based on observed saturation in independent cores. Using Spitzer GLIMPSE 8 micron images, this allows us to accurately probe mass surface densities, Sigma, up to ~0.5g/cm^2 with 2" resolution. We then characterize the structure of 42 massive starless and early-stage IRDC cores and their surrounding clumps, measuring Sigma_cl(r) from the core/clump centers. We first assess the properties of the core/clump at a scale where the total enclosed mass as projected on the sky is M_cl=60Msun. We find these objects have a mean radius of R_cl~0.1pc, mean Sigma_cl=0.3g/cm^2 and, if fit by a power law density profile rho_cl ~ r^{-k_{rho,cl}}, a mean value of k_{rho,cl}=1.1. If we assume a core is embedded in each clump and subtract the surrounding clump envelope to derive the core properties, we find a mean core density power law index of k_{rho,c} = 1.6. We repeat this analysis as a function of radius and derive the best-fitting power law plus uniform clump envelope model for each of the 42 core/clumps. The cores have typical masses of M_c~100Msun and mean Sigma_c~0.1g/cm^2, and are embedded in clumps with comparable mass surface densities. We conclude massive starless cores exist and are well-described by singular polytropic spheres. Their relatively low values of Sigma and the fact that they are IR dark may imply that their fragmentation is inhibited by magnetic fields rather than radiative heating. Comparing to massive star-forming cores, there is tentative evidence for an evolution towards higher densities and steeper density profiles as star formation proceeds.

Light transport in graded index media follows a curved trajectory determined by the Fermat's principle. Besides the effect of variation of the refractive index on the transport of radiative intensity, the curved ray trajectory will induce geometrical effects on the transport of polarization ellipse. This paper presents a complete derivation of vector radiative transfer equation for polarized radiation transport in absorption, emission and scattering graded index media. The derivation is based on the analysis of the conserved quantities for polarized light transport along curved trajectory and a novel approach. The obtained transfer equation can be considered as a generalization of the classic vector radiative transfer equation that is only valid for uniform refractive index media. Several variant forms of the transport equation are also presented, which include the form for Stokes parameters defined with a fixed reference and the Eulerian forms in the ray coordinate and in several common orthogonal coordinate systems.

New laboratory measurements using an Electron Beam Ion Trap (EBIT) and an x-ray microcalorimeter are presented for the n=3 to n=2 Fe XVII emission lines in the 15 {\AA} to 17 {\AA} range, along with new theoretical predictions for a variety of electron energy distributions. This work improves upon our earlier work on these lines by providing measurements at more electron impact energies (seven values from 846 to 1185 eV), performing an in situ determination of the x-ray window transmission, taking steps to minimize the ion impurity concentrations, correcting the electron energies for space charge shifts, and estimating the residual electron energy uncertainties. The results for the 3C/3D and 3s/3C line ratios are generally in agreement with the closest theory to within 10%, and in agreement with previous measurements from an independent group to within 20%. Better consistency between the two experimental groups is obtained at the lowest electron energies by using theory to interpolate, taking into account the significantly different electron energy distributions. Evidence for resonance collision effects in the spectra is discussed. Renormalized values for the absolute cross sections of the 3C and 3D lines are obtained by combining previously published results, and shown to be in agreement with the predictions of converged R-matrix theory. This work establishes consistency between results from independent laboratories and improves the reliability of these lines for astrophysical diagnostics. Factors that should be taken into account for accurate diagnostics are discussed, including electron energy distribution, polarization, absorption/scattering, and line blends.

Infrared Dark Clouds (IRDCs) are cold, high mass surface density and high density structures, likely to be representative of the initial conditions for massive star and star cluster formation. CO emission from IRDCs has the potential to be useful for tracing their dynamics, but may be affected by depleted gas phase abundances due to freeze-out onto dust grains. Here we analyze C18O J=1-0 and J=2-1 emission line data, taken with the IRAM 30m telescope, of the highly filamentary IRDC G035.39.-0033. We derive the excitation temperature as a function of position and velocity, with typical values of ~7K, and thus derive total mass surface densities, Sigma_C18O, assuming standard gas phase abundances and accounting for optical depth in the line, which can reach values of ~1. The mass surface densities reach values of ~0.07 g/cm^2. We compare these results to the mass surface densities derived from mid-infrared (MIR) extinction mapping, Sigma_SMF, by Butler & Tan, which are expected to be insensitive to the dust temperatures in the cloud. With a significance of >10sigma, we find Sigma_C18O/Sigma_SMF decreases by about a factor of 5 as Sigma increases from ~0.02 to ~0.2 g/cm^2, which we interpret as evidence for CO depletion. Several hundred solar masses are being affected, making this one of the most massive clouds in which CO depletion has been observed directly. We present a map of the depletion factor in the filament and discuss implications for the formation of the IRDC.

Thermal images of cold dust in the Central Molecular Zone of the Milky Way, obtained with the far-infrared cameras on-board the Herschel satellite, reveal a 3x10^7 solar masses ring of dense and cold clouds orbiting the Galactic Center. Using a simple toy-model, an elliptical shape having semi-major axes of 100 and 60 parsecs is deduced. The major axis of this 100-pc ring is inclined by about 40 degrees with respect to the plane-of-the-sky and is oriented perpendicular to the major axes of the Galactic Bar. The 100-pc ring appears to trace the system of stable x_2 orbits predicted for the barred Galactic potential. Sgr A* is displaced with respect to the geometrical center of symmetry of the ring. The ring is twisted and its morphology suggests a flattening-ratio of 2 for the Galactic potential, which is in good agreement with the bulge flattening ratio derived from the 2MASS data.

Theory predicts, and observations confirm, that the column density ratio of a molecule containing D to its counterpart containing H can be used as an evolutionary tracer in the low-mass star formation process. Since it remains unclear if the high-mass star formation process is a scaled-up version of the low-mass one, we investigated whether the relation between deuteration and evolution can be applied to the high-mass regime. With the IRAM-30m telescope, we observed rotational transitions of N2D+ and N2H+ and derived the deuterated fraction in 27 cores within massive star-forming regions understood to represent different evolutionary stages of the massive-star formation process. Results. Our results clearly indicate that the abundance of N2D+ is higher at the pre-stellar/cluster stage, then drops during the formation of the protostellar object(s) as in the low-mass regime, remaining relatively constant during the ultra-compact HII region phase. The objects with the highest fractional abundance of N2D+ are starless cores with properties very similar to typical pre-stellar cores of lower mass. The abundance of N2D+ is lower in objects with higher gas temperatures as in the low-mass case but does not seem to depend on gas turbulence. Our results indicate that the N2D+-to-N2H+ column density ratio can be used as an evolutionary indicator in both low- and high-mass star formation, and that the physical conditions influencing the abundance of deuterated species likely evolve similarly during the processes that lead to the formation of both low- and high-mass stars.

Radiative transfer calculations of massive star formation are presented. These are based on the Turbulent Core Model of McKee & Tan and self-consistently included a hydrostatic core, an inside-out expansion wave, a zone of free-falling rotating collapse, wide-angle dust-free outflow cavities, an active accretion disk, and a massive protostar. For the first time for such models, an optically thick inner gas disk extends inside the dust destruction front. This is important to conserve the accretion energy naturally and for its shielding effect on the outer region of the disk and envelope. The simulation of radiation transfer is performed with the Monte Carlo code of Whitney, yielding spectral energy distributions (SEDs) for the model series, from the simplest spherical model to the fiducial one, with the above components each added step-by-step. Images are also presented in different wavebands of various telescope cameras, including Spitzer IRAC and MIPS, SOFIA FORCAST and Herschel PACS and SPIRE. The existence of the optically thick inner disk produces higher optical wavelength fluxes but reduces near- and mid-IR emission. The presence of outflow cavities, the inclination angle to the line of sight, and the thickness of the disk all affect the SEDs and images significantly. For the high mass surface density cores considered here, the mid-IR emission can be dominated by the outflow cavity walls, as has been suggested by De Buizer. The effect of varying the pressure of the environment bounding the surface of the massive core is also studied. With lower surface pressures, the core is larger, has lower extinction and accretion rates, and the observed mid-IR flux from the disk can then be relatively high even though the accretion luminosity is lower. In this case the silicate absorption feature becomes prominent, in contrast to higher density cores forming under higher pressures.

The Census of High- and Medium-mass Protostars (CHaMP) is the first large-scale, unbiased, uniform mapping survey at sub-parsec scale resolution of 90 GHz line emission from massive molecular clumps in the Milky Way. We present the first Mopra (ATNF) maps of the CHaMP survey region (300{\deg}>l>280{\deg}) in the HCO+ J=1-0 line, which is usually thought to trace gas at densities up to 10^11 m-3. In this paper we introduce the survey and its strategy, describe the observational and data reduction procedures, and give a complete catalogue of moment maps of the HCO+ J=1-0 emission from the ensemble of 301 massive molecular clumps. From these maps we also derive the physical parameters of the clumps, using standard molecular spectral-line analysis techniques. This analysis yields the following range of properties: integrated line intensity 1-30 K km s-1, peak line brightness 1-7 K, linewidth 1-10 km s-1, integrated line luminosity 0.5-200 K km s-1 pc^2, FWHM size 0.2-2.5 pc, mean projected axial ratio 2, optical depth 0.08-2, total surface density 30-3000 M{\sun} pc-2, number density 0.2-30 x 10^9 m-3, mass 15-8000 M{\sun}, virial parameter 1-55, and total gas pressure 0.3-700 pPa. We find that the CHaMP clumps do not obey a Larson-type size-linewidth relation. Among the clumps, there exists a large population of subthermally excited, weakly-emitting (but easily detectable) dense molecular clumps, confirming the prediction of Narayanan et al. (2008). These weakly-emitting clumps comprise 95% of all massive clumps by number, and 87% of the molecular mass, in this portion of the Galaxy; their properties are distinct from the brighter massive star-forming regions that are more typically studied. If the clumps evolve by slow contraction, the 95% of fainter clumps may represent a long-lived stage of pressure-confined, gravitationally stable massive clump evolution, while the CHaMP ... (abridged)

Stars and more particularly massive stars, have a drastic impact on galaxy evolution. Yet the conditions in which they form and collapse are still not fully understood. In particular, the influence of the magnetic field on the collapse of massive clumps is relatively unexplored, it is thus of great relevance in the context of the formation of massive stars to investigate its impact. We perform high resolution, MHD simulations of the collapse of hundred solar masses, turbulent and magnetized clouds, using the adaptive mesh refinement code RAMSES. We compute various quantities such as mass distribution, magnetic field and angular momentum within the collapsing core and study the episodic outflows and the fragmentation that occurs during the collapse. The magnetic field has a drastic impact on the cloud evolution. We find that magnetic braking is able to substantially reduce the angular momentum in the inner part of the collapsing cloud. Fast and episodic outflows are being launched with typical velocities of the order of 3-5 km s$^{-1}$ although the highest velocities can be as high as 30-40 km s$^{-1}$. The fragmentation in several objects, is reduced in substantially magnetized clouds with respect to hydrodynamical ones by a factor of the order of 1.5-2. We conclude that magnetic fields have a significant impact on the evolution of massive clumps. In combination with radiation, magnetic fields largely determine the outcome of massive core collapse. We stress that numerical convergence of MHD collapse is a challenging issue. In particular, numerical diffusion appears to be important at high density therefore possibly leading to an over-estimation of the number of fragments.

The dense, cold gas of Infrared Dark Clouds (IRDCs) is thought to be representative of the initial conditions of massive star and star cluster formation. We analyze 13CO(J=1-0) line emission data from the Galactic Ring Survey of Jackson et al. for two filamentary IRDCs, comparing the mass surface densities derived from 13CO, Sigma_13CO, with those derived from mid-infrared small median filter extinction mapping, Sigma_SMF, by Butler & Tan. After accounting for molecular envelopes around the filaments, we find approximately linear relations between Sigma_CO and Sigma_SMF, i.e. an approximately constant ratio Sigma_CO/Sigma_SMF in the clouds. There is a variation of about a factor of two between the two clouds. We find evidence for a modest decrease of Sigma_CO/Sigma_SMF with increasing mass surface density, which may be due to a systematic decrease in temperature, increase in importance of high 13CO opacity cores, increase in dust opacity, or decrease in 13CO abundance due to depletion in regions of higher column density. We perform ellipsoidal and filamentary virial analyses of the clouds, finding that the surface pressure terms are dynamically important and that globally the filaments may not yet have reached virial equilibrium. Some local regions along the filaments appear to be close to virial equilibrium, although still with dynamically important surface pressures, and these appear to be sites where star formation is most active.

We review theoretical models of Population III.1 star formation, focusing on the protostellar feedback processes that are expected to terminate accretion and thus set the mass of these stars. We discuss how dark matter annihilation may modify this standard feedback scenario. Then, under the assumption that dark matter annihilation is unimportant, we predict the mass of stars forming in 12 cosmological minihalos produced in independent numerical simulations. This allows us to make a simple estimate of the Pop III.1 initial mass function and how it may evolve with redshift.

The "Extreme starbursts in the local universe" workshop was held at the Insituto de Astrofisica de Andalucia in Granada, Spain on 21-25 June 2010. Bearing in mind the advent of a new generation of facilities such as JWST, Herschel, ALMA, eVLA and eMerlin, the aim of the workshop was to bring together observers and theorists to review the latest results. The purpose of the workshop was to address the following issues: what are the main modes of triggering extreme starbursts in the local Universe? How efficiently are stars formed in extreme starbursts? What are the star formation histories of local starburst galaxies? How well do the theoretical simulations model the observations? What can we learn about starbursts in the distant Universe through studies of their local counterparts? How important is the role of extreme starbursts in the hierarchical assembly of galaxies? How are extreme starbursts related to the triggering of AGN in the nuclei of galaxies? Overall, 41 talks and 4 posters with their corresponding 10 minutes short talks were presented during the workshop. In addition, the workshop was designed with emphasis on discussions, and therefore, there were 6 discussion sessions of up to one hour during the workshop. Here is presented a summary of the purposes of the workshop as well as a compilation of the abstracts corresponding to each of the presentations. The summary and conclusions of the workshop along with a description of the future prospects by Sylvain Veilleux can be found in the last section of this document. A photo of the assistants is included.

We use observations of the radial profiles of the mass surface density of total, Sigma_g, and molecular, Sigma_H2, gas, rotation velocity and star formation rate surface density, Sigma_sfr, of the molecular dominated regions of 12 disk galaxies from Leroy et al. to test several star formation laws: a "Kennicutt-Schmidt power law", Sigma_sfr=A_g Sigma_{g,2}^{1.5}$; a "Constant molecular law", Sigma_sfr = A_H2 Sigma_{H2,2}; the "Turbulence-regulated laws" of Krumholz & McKee (KM) and Krumholz, McKee & Tumlinson (KMT), a "Gas-Omega law", Sigma_sfr = B_Omega Sigma_g Omega; and a shear-driven "GMC collisions law", Sigma_sfr = B_CC Sigma_g Omega (1 - 0.7beta), where beta is d ln v_circ / d ln r. We find the constant molecular law, KMT turbulence law and GMC collision law are the most accurate, with an rms error of a factor of 1.5 if the normalization constants are allowed to vary between galaxies. Of these three laws, the GMC collision law does not require a change in physics to account for the full range of star formation activity seen from normal galaxies to circumnuclear starbursts. A single global GMC collision law with B_CC=8.0x10^{-3}, i.e. a gas consumption time of 20 orbital times for beta=0, yields an rms error of a factor of 1.8.

Source I in the Orion-KL nebula is believed to be the nearest example of a massive star still in the main accretion phase. It is thus one of the best cases for studying the properties of massive protostars to constrain high-mass star formation theories. Near-infrared radiation from source I escapes through the cavity opened by the OMC1 outflow and is scattered by dust towards our line of sight. The reflected spectrum offers a unique possibility of observing the emission from the innermost regions of the system and probing the nature of source I and its immediate surroundings. We obtained moderately high spectral-resolution (R~9000) observations of the near-infrared diffuse emission in several locations around source I/Orion-KL. We observed a widespread rich absorption line spectrum that we compare with cool stellar photospheres and protostellar accretion disk models. The spectrum is broadly similar to strongly veiled, cool, low-gravity stellar photospheres in the range Teff~3500-4500 K, luminosity class I-III. An exact match explaining all features has not been found, and a plausible explanation is that a range of different temperatures contribute to the observed absorption spectrum. The 1D velocity dispersions implied by the absorption spectra, sigma~30 km/s, can be explained by the emission from a disk around a massive, mstar~10 Msun, protostar that is accreting at a high rate, mdot~3x10^{-3} Msun/yr. Our observations suggest that the near-infrared reflection spectrum observed in the Orion-KL region is produced close to source I and scattered to our line of sight in the OMC1 outflow cavity. The spectrum allows us to exclude source I being a very large, massive protostar rotating at breakup speed. We suggest that the absorption spectrum is produced in a disk surrounding a ~10 Msun protostar, accreting from its disk at a high rate of a few 10^{-3} Msun/yr.

We present high-sensitivity 2'x4' maps of the J=2-1 rotational lines of SiO, CO, 13CO and C18O, observed toward the filamentary Infrared Dark Cloud (IRDC) G035.39-00.33. Single-pointing spectra of the SiO J=2-1 and J=3-2 lines toward several regions in the filament, are also reported. The SiO images reveal that SiO is widespread along the IRDC (size >2 pc), showing two different components: one bright and compact arising from three condensations (N, E and S), and the other weak and extended along the filament. While the first component shows broad lines (linewidths of ~4-7 kms-1) in both SiO J=2-1 and SiO J=3-2, the second one is only detected in SiO J=2-1 and has narrow lines (~0.8 kms-1). The maps of CO and its isotopologues show that low-density filaments are intersecting the IRDC and appear to merge toward the densest portion of the cloud. This resembles the molecular structures predicted by flow-driven, shock-induced and magnetically-regulated cloud formation models. As in outflows associated with low-mass star formation, the excitation temperatures and fractional abundances of SiO toward N, E and S, increase with velocity from ~6 to 40 K, and from ~1E-10 to >1E-8 respectively, over a velocity range of ~7 kms-1. Since 8 micron sources, 24 micron sources and/or extended 4.5 micron emission are detected in N, E and S, broad SiO is likely produced in outflows associated with high-mass protostars. The excitation temperatures and fractional abundances of the narrow SiO lines, however, are very low (~9 K and ~1E-11, respectively), and consistent with the processing of interstellar grains by the passage of a shock with vs~12 kms-1. This emission could be generated i) by a large-scale shock, perhaps remnant of the IRDC formation process; ii) by decelerated or recently processed gas in large-scale outflows driven by 8 micron and 24 micron sources; or iii) by an undetected and widespread population of lower mass protostars. High-angular resolution observations are needed to disentangle between these three scenarios.

We investigate the formation and evolution of giant molecular clouds (GMCs) in a Milky-Way-like disk galaxy with a flat rotation curve. We perform a series of 3D adaptive mesh refinement (AMR) numerical simulations that follow both the global evolution on scales of ~20kpc and resolve down to scales ~<10pc with a multiphase atomic interstellar medium (ISM). In this first study, we omit star formation and feedback, and focus on the processes of gravitational instability and cloud collisions and interactions. We define clouds as regions with n_H>=100cm^-3 and track the evolution of individual clouds as they orbit through the galaxy from their birth to their eventual destruction via merger or via destructive collision with another cloud. After ~140Myr a large fraction of the gas in the disk has fragmented into clouds with masses ~10^6 Msun and a mass spectrum similar to that of Galactic GMCs. The disk settles into a quasi steady state in which gravitational scattering of clouds keeps the disk near the threshold of global gravitational instability. The cloud collision time is found to be a small fraction, ~1/5, of the orbital time, and this is an efficient mechanism to inject turbulence into the clouds. This helps to keep clouds only moderately gravitationally bound, with virial parameters of order unity. Many other observed GMC properties, such as mass surface density, angular momentum, velocity dispersion, and vertical distribution, can be accounted for in this simple model with no stellar feedback.

(Abridged) We use 8 micron Spitzer GLIMPSE images to make extinction maps of 10 IRDCs, selected to be relatively nearby and massive. The extinction mapping technique requires modeling the IR background intensity behind the cloud, which is achieved by correcting for foreground emission and then interpolating from the surrounding regions. The correction for foreground emission can be quite large, thus restricting the utility of this technique to relatively nearby clouds. We investigate three methods for the interpolation, finding systematic differences at about the 10% level, which, for fiducial dust models, corresponds to a mass surface density Sigma = 0.013 g cm^-2, above which we conclude this extinction mapping technique attains validity. We examine the probability distribution function of Sigma in IRDCs. From a qualitative comparison with numerical simulations of astrophysical turbulence, many clouds appear to have relatively narrow distributions suggesting relatively low (<5) Mach numbers and/or dynamically strong magnetic fields. Given cloud kinematic distances, we derive cloud masses. Rathborne, Jackson & Simon identified cores within the clouds and measured their masses via mm dust emission. For 43 cores, we compare these mass estimates with those derived from our extinction mapping, finding good agreement: typically factors of <~2 difference for individual cores and an average systematic offset of <~10% for the adopted fiducial assumptions of each method. We find tentative evidence for a systematic variation of these mass ratios as a function of core density, which is consistent with models of ice mantle formation on dust grains and subsequent grain growth by coagulation, and/or with a temperature decrease in the densest cores.

We have surveyed molecular line emission from Orion BN/KL from 42.3 to 43.6 GHz with the Green Bank Telescope. Sixty-seven lines were identified and ascribed to 13 different molecular species. The spectrum at 7 mm is dominated by SiO, SO2, CH3OCH3, and C2H5CN. Five transitions have been detected from the SiO isotopologues 28SiO, 29SiO, and 30SiO. We report here for the first time the spectra of the 29SiO and 30SiO v=0 J=1-0 emission in Orion BN/KL, and we show that they have double-peaked profiles with velocity extents similar to the main isotopologue. The main motivation for the survey was the search of high-velocity (100-1000 km/s) outflows in the BN/KL region as traced by SiO Doppler components. Some of the unidentified lines in principle could be high-velocity SiO features, but without imaging data their location cannot be established. Wings of emission are present in the v=0 28SiO, 29SiO and 30SiO profiles, and we suggest that the v=0 emission from the three isotopologues might trace a moderately high-velocity (~30-50 km/s) component of the flows around the high-mass protostar Source I in the Orion BN/KL region. We also confirm the 7 mm detection of a complex oxygen-bearing species, acetone (CH3COCH3), which has been recently observed towards the hot core at 3 mm, and we have found further indications of the presence of long cyanopolyynes (HC5N and HC7N) in the quiescent cold gas of the extended ridge.

We show that the proper motion of the Becklin-Neugebauer (BN) object is consistent with its dynamical ejection from the theta^1 Ori C binary, contrary to recent claims by Gomez et al. Continued radio observations of BN and future precise astrometric observations of theta^1 Ori C with SIM and the Orion Nebula Cluster with GAIA can constrain the properties of this ejection event, with implications for theories of how the nearest example of massive star formation is proceeding.

I discuss current theoretical expectations of how primordial, Pop III.1 stars form. Lack of direct observational constraints makes this a challenging task. In particular predicting the mass of these stars requires solving a series of problems, which all affect, perhaps drastically, the final outcome. While there is general agreement on the initial conditions, H_2-cooled gas at the center of dark matter minihalos, the subsequent evolution is more uncertain. In particular, I describe the potential effects of dark matter annihilation heating, fragmentation within the minihalo, magnetic field amplification, and protostellar ionizing feedback. After these considerations, one expects that the first stars are massive >~100Msun, with dark matter annihilation heating having the potential to raise this scale by large factors. Higher accretion rates in later-forming minihalos may cause the Pop III.1 initial mass function to evolve to higher masses.

Supermassive black holes in giant elliptical galaxies are remarkably faint given their expected accretion rates. This motivates models of radiatively inefficient accretion, due to either ion-electron thermal decoupling, generation of outflows that inhibit accretion, or settling of gas to a gravitationally unstable disk that forms stars in preference to feeding the black hole. The latter model predicts the presence of cold molecular gas in a thin disk around the black hole. Here we report Submillimeter Array observations of the nucleus of the giant elliptical galaxy M87 that probe 230 GHz continuum and CO(J=2--1) line emission. Continuum emission is detected from the nucleus and several knots in the jet, including one that has been undergoing flaring behavior. We estimate a conservative upper limit on the mass of molecular gas within ~100pc and +-400km/s line of sight velocity of the central black hole of ~8x10^6Msun, which includes an allowance for possible systematic errors associated with subtraction of the continuum. Ignoring such errors, we have a 3 sigma sensitivity to about 3x10^6Msun. In fact, the continuum-subtracted spectrum shows weak emission features extending up to 4 sigma above the RMS dispersion of the line-free channels. These may be artifacts of the continuum subtraction process. Alternatively, if they are interpreted as CO emission, then the implied molecular gas mass is ~5x10^6Msun spread out over a velocity range of 700km/s. These constraints on molecular gas mass are close to the predictions of the model of self-gravitating, star-forming accretion disks fed by Bondi accretion (Tan & Blackman 2005).

We investigate the effects of weakly-interacting massive particle (WIMP) dark matter annihilation on the formation of Population III.1 stars, which are theorized to form from the collapse of gas cores at the centers of dark matter minihalos. We consider the relative importance of cooling due to baryonic radiative processes and heating due to WIMP annihilation. We analyze the dark matter and gas profiles of several halos formed in cosmological-scale numerical simulations. The heating rate depends sensitively on the dark matter density profile, which we approximate with a power law rho_chi ~ r^{-alpha_chi}, in the numerically unresolved inner regions of the halo. If we assume a self-similar structure so that alpha_chi ~= 1.5 as measured on the resolved scales ~1pc, then for a fiducial WIMP mass of 100GeV, the heating rate is typically much smaller (<10^{-3}) than the cooling rate for densities up to n_H=10^{17}cm^{-3}. In one case, where alpha_chi=1.65, the heating rate becomes similar to the cooling rate by a density of n_H=10^{15}cm^{-3}. The dark matter density profile is expected to steepen in the central baryon-dominated region <~1pc due to adiabatic contraction, and we observe this effect (though with relatively low resolution) in our numerical models. From these we estimate alpha_chi~=2.0. The heating now dominates cooling above n_H~=10^{14}cm^{-3}, in agreement with the previous study of Spolyar, Freese & Gondolo. We expect this leads to the formation of an equilibrium structure with a baryonic and dark matter density distribution exhibiting a flattened central core. Examining such equilibria, we find total luminosities due to WIMP annihilation are relatively constant and ~10^3 L_sun, set by the radiative luminosity of the baryonic core. We discuss the implications for Pop III.1 star formation... (abridged)

We consider the radiative feedback processes that operate during the formation of the first stars, including the photodissociation of H_2, Ly-alpha radiation pressure, formation and expansion of an HII region, and disk photoevaporation. These processes may inhibit continued accretion once the stellar mass has reached a critical value, and we evaluate this mass separately for each process. Photodissociation of H_2 in the local dark matter minihalo occurs relatively early in the growth of the protostar, but we argue this does not affect subsequent accretion since by this time the depth of the potential is large enough for accretion to be mediated by atomic cooling. However, neighboring starless minihalos can be affected. Ionization creates an HII region in the infalling envelope above and below the accretion disk. Ly-alpha radiation pressure acting at the boundary of the HII region is effective at reversing infall from narrow polar directions when the star reaches ~20-30Msun, but cannot prevent infall from other directions. Expansion of the HII region beyond the gravitational escape radius for ionized gas occurs at masses ~50-100Msun, depending on the accretion rate and angular momentum of the inflow. However, again, accretion from the equatorial regions can continue since the neutral accretion disk has a finite thickness and shields a substantial fraction of the accretion envelope from direct ionizing flux. At higher stellar masses, ~140Msun in the fiducial case, the combination of declining accretion rates and increasing photoevaporation-driven mass loss from the disk act to effectively halt the increase in the protostellar mass. We identify this process as the mechanism that terminates the growth of Population III stars... (abridged)

Until recently, simulations that modeled entire galaxies were restricted to an isothermal or fixed 2- or 3-phase interstellar medium (ISM). This obscured the full role of the ISM in shaping the observed galactic-scale star formation relations. In particular the Kennicutt relation suggests that star formation rates depend in a simple way on global galactic quantities, such as mean gas mass surface density and dynamical time. Contrary to this, observations of nearby star-forming regions, including images from the Spitzer telescope, show that all the way down to ``local-scales'' ($\sim$ parsecs) star formation is a highly clustered process with the gas existing at a wide range of densities and multiple phases, and seemingly decoupled from the larger-scale Galaxy. Many different physical processes appear to be influencing the star formation rate, including heating and cooling of the gas, ``turbulent'' energy injection from outflows, winds and supernova blastwaves, ionization, and magnetic field support. Here we present results from global simulations of disk galaxy star formation that include a fully multiphase ISM and compare to simpler isothermal models. We also investigate the effect of background heating. We discuss the process of giant molecular cloud formation as one of the possible links between the galactic and local scales of star formation.

The initial cluster mass function (ICMF) is a fundamental property of star formation in galaxies. To gauge its universality, we measure and compare the ICMFs in irregular and spiral galaxies. Our sample of irregular galaxies is based on 13 nearby galaxies selected from a volume-limited sample from the fifth data release of the Sloan Digital Sky Survey (SDSS). The extinctions, ages, and masses were determined by comparing their u'g'i'z' magnitudes to those generated from starburst models. Completeness corrections were performed using Monte Carlo simulations in which artificial clusters were inserted into each galaxy. We analyzed three nearby spiral galaxies with SDSS data in exactly the same way to derive their ICMF based on a similar number of young, massive clusters as the irregular galaxy ICMF. We find that the ICMFs of irregular and spiral galaxies for masses >3x10^4 M_sun are statistically indistinguishable. For clusters more massive than 3x10^4 M_sun, the ICMF of the irregular galaxies is reasonably well fit by a power law dN(M)/dM ~ M^-a_M with a_M = 1.88 +/- 0.09. Similar results were obtained for the ICMF of the spiral galaxy sample but with a_M = 1.75 +/- 0.06. We discuss the implications of our result for theories of star cluster formation: the shape of the ICMF appears to be independent of metallicity and galactic shear rate.

I discuss different theories of massive star formation: formation from massive cores, competitive Bondi-Hoyle accretion, and protostellar collisions. I summarize basic features of the Turbulent Core Model (TCM). I then introduce the Orion Kleinmann-Low (KL) region, embedded in the Orion Nebula Cluster (ONC) and one of the nearest regions of massive star formation. The KL region contains three principal radio sources, known as "I", "n" and "BN". BN is known to be a runaway star, almost certainly set in motion by dynamical ejection within the ONC from a multiple system of massive stars, that would leave behind a recoiling, hard, massive, probably eccentric binary. I review the debate about whether this binary is Theta^1C, the most massive star in the ONC, or source "I", and argue that it is most likely to be Theta^1C, since this is now known be a recoiling, hard, massive, eccentric binary, with properties that satisfy the energy and momentum constraints implied by BN's motion. Source "n" is a relatively low-mass protostar with extended radio emission suggestive of a bipolar outflow. Source "I", located near the center of the main gas concentration in the region, the Orion Hot Core, is the likely location of a massive protostar that is powering the KL region, and I discuss how its basic properties are consistent with predictions from the TCM. In this scenario, the radio emission from source "I" is the base of a bipolar outflow that is ionized by the massive protostar and should be elongated along the axis of the outflow.

We describe how star formation is expected to proceed in the early metal-free Universe, focusing on the very first generations of stars. We then discuss how the star formation process may change as the effects of metallicity, external radiative feedback, and magnetic and turbulent support of the gas become more important. The very first stars (Pop III.1) have relatively simple initial conditions set by cosmology and the cooling properties of primordial gas. We describe the evolution of these stars as they grow in mass by accretion from their surrounding gas cores and how the accretion process is affected and eventually terminated by radiative feedback processes, especially HII region expansion and disk photoevaporation. The ability of the protostar and its disk to generate dynamically important magnetic fields is reviewed and their effects discussed. Pop III.1 star formation is likely to produce massive (~100-200Msun) stars that then influence their surroundings via ionization, stellar winds, and supernovae. These processes heat, ionize and metal-enrich the gas, thus altering the initial conditions for the next generation of star formation. Stars formed from gas that has been altered significantly by radiative and/or mechanical feedback, but not by metal enrichment (Pop III.2) are expected to have significantly smaller masses than Pop III.1 stars because of more efficient cooling from enhanced HD production. Stars formed from gas that is metal-enriched to levels that affect the dynamics of the collapse (the first Pop II stars) are also expected to have relatively low masses. We briefly compare the above star formation scenarios to what is known about present-day star formation.

Supermassive black hole accretion and star formation appear intimately connected. I review the observational and theoretical evidence for this statement. I then discuss how focussed studies of two systems, our Galactic Center and the nucleus of M87, can help to improve our understanding of these processes.

I review the status of massive star formation theories: accretion from collapsing, massive, turbulent cores; competitive accretion; and stellar collisions. I conclude the observational and theoretical evidence favors the first of these models. I then discuss: the initial conditions of star cluster formation as traced by infrared dark clouds; the cluster formation timescale; and comparison of the initial cluster mass function in different galactic environments.

It has been known for more than 30 years that star formation in giant molecular clouds (GMCs) is slow, in the sense that only ~1% of the gas forms stars every free-fall time. This result is entirely independent of any particular model of molecular cloud lifetime or evolution. Here we survey observational data on higher density objects in the interstellar medium, including infrared dark clouds and dense molecular clumps, to determine if these objects form stars slowly like GMCs, or rapidly, converting a significant fraction of their mass into stars in one free-fall time. We find no evidence for a transition from slow to rapid star formation in structures covering three orders of magnitude in density. This has important implications for models of star formation, since competing models make differing predictions for the characteristic density at which star formation should transition from slow to rapid. The data are inconsistent with models that predict that star clusters form rapidly and in free-fall collapse. Magnetic- and turbulence-regulated star formation models can reproduce the observations qualitatively, and the turbulence-regulated star formation model of Krumholz & McKee quantitatively reproduces the infrared-HCN luminosity correlation recently reported by Gao & Solomon. Slow star formation also implies that the process of star cluster formation cannot be one of global collapse, but must instead proceed over many free-fall times. This suggests that turbulence in star-forming clumps must be driven, and that the competitive accretion mechanism does not operate in typical cluster-forming molecular clumps.

Aims. To understand the roles of infall and protostellar evolution on the envelopes of massive young stellar objects (YSOs). Methods. The chemical evolution of gas and dust is traced, including infall and realistic source evolution. The temperatures are determined self-consistently. Both ad/desorption of ices using recent laboratory temperature-programmed-desorption measurements are included. Results. The observed water abundance jump near 100 K is reproduced by an evaporation front which moves outward as the luminosity increases. Ion-molecule reactions produce water below 100 K. The age of the source is constrained to t \~ 8 +/- 4 x 10^4 yrs since YSO formation. It is shown that the chemical age-dating of hot cores at ~ few x 10^3 - 10^4 yr and the disappearance of hot cores on a timescale of ~ 10^5 yr is a natural consequence of infall in a dynamic envelope and protostellar evolution. Dynamical structures of ~ 350AU such as disks should contain most of the complex second generation species. The assumed order of desorption kinetics does not affect these results.

After a brief introduction to the NIST EBIT facility, we present the results of three different types of experiments that have been carried out there recently: EUV and visible spectroscopy in support of the microelectronics industry, laboratory astrophysics using an x-ray microcalorimeter, and charge exchange studies using extracted beams of highly charged ions.

We argue that rich star clusters take at least several local dynamical times to form, and so are quasi-equilibrium structures during their assembly. Observations supporting this conclusion include morphologies of star-forming clumps, momentum flux of protostellar outflows from forming clusters, age spreads of stars in the Orion Nebula Cluster (ONC) and other clusters, and the age of a dynamical ejection event from the ONC. We show that these long formation timescales are consistent with the expected star formation rate in turbulent gas, as recently evaluated by Krumholz & McKee. Finally, we discuss the implications of these timescales for star formation efficiencies, the disruption of gas by stellar feedback, mass segregation of stars, and the longevity of turbulence in molecular clumps.

Massive stars have a profound influence on the Universe, but their formation remains poorly understood. We review the current status of observational and theoretical research in this field, describing the various stages of an evolutionary sequence that begins with cold, massive gas cores and ends with the dispersal and ionization of gas by the newly-formed star. The physical processes in massive star formation are described and related to their observational manifestations. Feedback processes and the relation of massive stars to star cluster formation are also discussed. We identify key observational and theoretical questions that future studies should address.

The luminosities of the centers of nearby elliptical galaxies are very low compared to models of thin disc accretion to their black holes at the Bondi rate, typically a few hundredths to a few tenths of a solar mass per year. This has motivated models of inefficiently-radiated accretion that invoke weak electron-ion thermal coupling, and/or inhibited accretion rates due to convection or outflows. Here we point out that even if such processes are operating, a significant fraction of the accreting gas is prevented from reaching the central black hole because it condenses into stars in a gravitationally unstable disc. Star formation occurs inside the Bondi radius (typically ~100pc in giant ellipticals), but still relatively far from the black hole in terms of Schwarzschild radii. Star formation depletes and heats the gas disc, eventually leading to a marginally stable, but much reduced, accretion flow to the black hole. We predict the presence of cold (~100K), dusty gas discs, containing clustered H-alpha emission and occasional type II supernovae, both resulting from the presence of massive stars. Star formation accounts for several features of the M87 system: a thin disc, traced by H-alpha emission, is observed on scales of about 100pc, with features reminiscent of spiral arms and dust lanes; the star formation rate inferred from the intensity of H-alpha emission is consistent with the Bondi accretion rate of the system. Star formation may therefore help suppress accretion onto the central engines of massive ellipticals. We also discuss some implications for the fueling of the Galactic center and quasars.

I review some important questions in the field of massive star formation: What are the initial conditions for proto star clusters and how do they arise? What are the initial conditions for individual massive star formation within star clusters? How do massive protostars accumulate their mass? I compare the Turbulent Core Model (McKee & Tan 2003) to several nearby regions, including Orion KL. Here I also discuss the origin of BN's high proper motion.

In the present-day universe, it appears that most, and perhaps all, massive stars are born in star clusters. It also appears that all star clusters contain stars drawn from an approximately universal initial mass function, so that almost all rich young star clusters contain massive stars. In this review I discuss the physical processes associated with both massive star formation and with star cluster formation. First I summarize the observed properties of star-forming gas clumps, then address the following questions. How do these clumps emerge from giant molecular clouds? In these clustered environments, how do individual stars form and gain mass? Can a forming star cluster be treated as an equilibrium system or is this process too rapid for equilibrium to be established? How does feedback affect the formation process?

We consider whether the variability properties of Gamma-Ray Bursts (GRBs) that produce bright optical and longer wavelength transient afterglows (A-GRBs) are the same as a larger, inclusive sample of bright, long-duration GRBs, selected only by their gamma-ray emission. This sample may include a significant population of physically distinct ``dark'' or ``faint-afterglow'' GRBs with different variability properties, or may be composed of a single population, some of which lack afterglows only because of observational selection effects. We argue that the structure function is the most appropriate method for measuring the variability of bursts because of their transient and aperiodic nature. We define a simple statistic: the ratio of the integrated structure function from 0.1 to 1 s compared to that from 0.1 to 10 s, as measured in the observer frame. To avoid instrumental effects we restrict our analysis to GRBs with BATSE data. Comparing 10 A-GRBs to a ``main'' sample of about 500 bursts, we find there is a probability of only 0.03 of the samples being drawn from the same population, with the A-GRBs tending to have relatively less power on sub-second timescales. We conclude that this result is tentative evidence for variations in the properties of GRB progenitors that affect both the gamma-ray and afterglow properties of long-duration GRBs. In addition, our method of analyzing variability identifies a characteristic timescale of ~1 s, below which variability is suppressed, and finds a trend of increased short timescale variability at higher gamma-ray energies. The long-duration GRBs that we identify as having the most sub-second timescale variability, may be relatively bright examples of short-duration GRBs.

I discuss two questions about the origin of star clusters: How long does the process take? What is the mode of individual star formation? I argue that observations of Galactic star-forming regions, particularly the Orion Nebula Cluster (ONC), indicate that cluster formation often takes several Myr, which is many local dynamical timescales. Individual stars and binaries, including massive stars, appear to form from the collapse of gas cores.

We attempt to explain the properties of the Becklin-Neugebauer (BN) object as a runaway B star, as originally proposed by Plambeck et al. (1995). This is one of the best-studied bright infrared sources, located in the Orion Nebula Cluster -- an important testing ground for massive star formation theories. From radio observations of BN's proper motion, we trace its trajectory back to Trapezium star theta^1C, the most massive (45 Msun) in the cluster and a relatively tight (17 AU) visual binary with a B star secondary. This origin would be the most recent known runaway B star ejection event, occurring only \~4000 yr ago and providing a unique test of models of ejection from multiple systems of massive stars. Although highly obscured, we can constrain BN's mass (~7 Msun) from both its bolometric luminosity and the recoil of theta^1C. Interaction of a runaway B star with dense ambient gas should produce a compact wind bow shock. We suggest that X-ray emission from this shocked gas may have been seen by Chandra: the offset from the radio position is ~300 AU in the direction of BN's motion. Given this model, we constrain the ambient density, wind mass-loss rate and wind velocity. BN made closest approach to the massive protostar, source ``I'', 500 yr ago. This may have triggered enhanced accretion and thus outflow, consistent with previous interpretations of the outflow being a recent (~10^3 yr) "explosive" event.

It is well known that the energy input from massive stars dominates the thermal and mechanical heating of typical regions in the interstellar medium of galaxies. These effects are amplified tremendously in the immediate environment of young massive star clusters, which may contain thousands of O and B stars. We present models of star cluster formation that attempt to account for the interplay between feedback and self-gravity in forming clusters, with application to systems ranging from the Orion Nebula Cluster, the Galactic Center clusters like the Arches, and super star clusters.

We propose that supermassive stars may form in quasar accretion disks, and we discuss possible observational consequences. The structure and stability of very massive stars are reviewed. Because of high accretion rates, quasar disks are massive and the fringes of their optically luminous parts are prone to fragmentation. Starting from a few hundred solar masses, a dominant fragment will grow to the isolation mass, which is a significant fraction of the disk mass, more quickly than the fragment contracts onto the stellar main sequence. A gap will form in the disk and the star will migrate inward on the accretion timescale, which is comparable to the star's main sequence lifetime. By interrupting the gas supply to the inner disk, the gap may temporarily dim and redden the quasar. The final stages of stellar migration will be a strong source of low-frequency gravitational waves.

We consider the absorption and scattering of X-rays observed from the Galactic center. One objective is to characterize the intrinsic X-ray emission from the central black hole, SgrA*, in its quiescent and flaring states. We correct the fluxes observed by the Chandra and XMM telescopes for absorption and scattering, using realistic models for the dust grain size distribution, gas and dust abundances and spatial distributions. However, there are still large uncertainties. Since much of the intervening dust is close to SgrA*, the scattered halo of X-ray photons is very concentrated: its intensity can dominate the Chandra PSF inside 1", and so affects estimates of the point source flux. It also broadens the radial intensity profiles of Galactic center sources - observations of this effect can help constrain models of the spatial distribution of the dust. We estimate that the combined scattering halos from observed Galactic center sources within 10" of SgrA* contribute up to ~10% of the observed diffuse emission in this region. Dust-scattered photons suffer a time delay relative to the photons that arrive directly. For dust that is 100 pc towards us from the Galactic center, this delay is about 1000 s at angles of 1", and 100 ks at 10". We illustrate how the evolution of the scattering halo following X-ray flares from SgrA* or other sources can also help to constrain the dust's line of sight distribution. We discuss the implications of X-ray scattering halos for the intensity of diffuse emission that has been reported within a few arcseconds of SgrA*: in the most extreme, yet viable, model we consider, ~1/3 of it is due to dust scattering of an unresolved source. The remainder results from an extended source of emission.

Magnetically mediated disk outflows are a leading paradigm to explain winds and jets in a variety of astrophysical sources, but where do the fields come from? Since accretion of mean magnetic flux may be disfavored in a thin turbulent disk, and only fields generated with sufficiently large scale can escape before being shredded by turbulence, in situ field production is desirable. Nonlinear helical inverse dynamo theory can provide the desired fields for coronae and outflows. We discuss the implications for contemporary protostellar disks, where the MRI (magneto-rotational instability) can drive turbulence in the inner regions, and primordial protostellar disks, where gravitational instability drives the turbulence. We emphasize that helical dynamos are compatible with the magneto-rotational instability, and clarify the relationship between the two.

Are magnetic fields important in primordial star formation? Assuming that star formation occurs via an accretion disk that is turbulent, initially because of local gravitational instability, we calculate the disk structure for realistic accretion rates. We predict that local gravitational viscosity is able to drive accretion, without the disk fragmenting. We then estimate the rate of dynamo amplification of seed magnetic field. Turbulence in a stratified disk can be helical, with different signs of the helicity in each hemisphere. This provides a key ingredient for production of global scale magnetic fields whose sign of flux is sustained over many orbit times. The resulting fields can drive collimated protostellar outflows that reduce the star formation efficiency from the initial gas cloud, especially once the protostar has contracted to the main sequence, at ~100 Msun. We estimate that the outflows are powerful enough to eject some material from the host dark matter halo and to initiate relatively strong magnetization of the local intergalactic medium. Close to the protostar, the outflow acts to shield the disk and equatorial regions from radiative feedback, such as ionizing photons, and this may enable accretion up to relatively large stellar masses. We conclude that magnetic fields cannot be ignored from models of primordial star formation.

We present a theoretical model for primordial star formation. First we describe the structure of the initial gas cores as virialized, quasi-hydrostatic objects in accord with recent high resolution numerical studies. The accretion rate can then be related to characteristic densities and temperatures that are set by the cooling properties of molecular hydrogen. We allow for rotation of the gas core, assuming angular momentum conservation inside the sonic point of the flow. In the typical case, most mass then reaches the star via an accretion disk. The structure of the inner region of this disk is described with the standard theory of viscous disks, but with allowance for the substantial energies absorbed in ionizing and dissociating the gas. The size of the protostar and its luminosity depend upon the accretion rate, the energetics of the accreting gas, and the ability of the radiation to escape from the stellar accretion shock. We combine these models for the infall rate, inner disk structure, and protostellar evolution to predict the radiation field that is the basis for radiative feedback processes acting against infall (Paper II). For realistic initial angular momenta, the photosphere of the protostar is much smaller and hotter than in the spherical case, leading to stronger radiative feedback at earlier stages in the evolution. In particular, once the star is older than its Kelvin-Helmholtz time, contraction towards the main sequence causes a rapid increase in ionizing and far-ultraviolet luminosity at masses ~30Msun in the fiducial case. Since the cores out of which the first stars formed were much more massive than 30Msun and since feedback is dynamically unimportant at lower masses, we conclude that the first stars should have had masses >~30Msun.

I compare theoretical models of massive star formation with observations of the Orion Hot Core, which harbors one of the closest massive protostars. Although this region is complicated, many of its features (size, luminosity, accretion disk, HII region, outflow) may be understood by starting with a simple model in which the star forms from a massive gas core that is a coherent entity in approximate pressure balance with its surroundings. The dominant contribution to the pressure is from turbulent motions and magnetic fields. The collapse can be perturbed by interactions with other stars in the forming cluster, which may induce sporadic enhancements of the accretion and outflow rate.

If massive stars form by disk accretion, then bipolar outflows should be generated as in the case of low-mass star formation. High accretion rates lead to high outflow rates and make the wind density very large near the protostar. We therefore predict that massive protostars have very small, jet-like HII regions confined by their outflows; we identify these HII regions with "hypercompact" (<0.01 pc) HII regions. Their lifetime is approximately equal to the accretion timescale (~10^5 yr), much longer than the sound-crossing and dynamical timescales of the ionized region. We present an analytic description of the density distribution of the outflow, relate the overall mass loss rate to the accretion rate, and normalize to values appropriate to massive protostars. For a given ionizing luminosity, we calculate the extent of the HII region and its radio spectrum. A detailed comparison is made with observations of radio source "I" in the Orion Hot Core. The spectra and morphologies of many other sources are broadly consistent with this model. We argue that confinement of ionization allows disk and equatorial accretion to continue up to high masses, potentially overcoming a major difficulty of standard accretion models.

We use published ROSAT observations of the X-ray Nova V1974 Cygni 1992 to test a model for interstellar dust, consisting of a mixture of carbonaceous grains and silicate grains. The time-dependent X-ray emission from the nova is modelled, and X-ray scattering is calculated for a dust mixture with a realistic size distribution. Model results are compared with the scattered X-ray halos measured by ROSAT at 9 epochs. The observed X-ray halos are consistent with the halos calculated for the size distribution of Weingartner & Draine which reproduces the Milky Way extinction with R_V=3.1, provided that the reddening E(B-V) approx 0.20, consistent with E(B-V) = 0.19 inferred from the late-time Balmer decrement. The time delay of the scattered halo relative to the direct flux from the nova is clearly detected. Models with smoothly-distributed dust give good overall agreement with the observed scattering halo, but tend to produce somewhat more scattering than observed at 200-300", and insufficient scattering at 50-100". While an additional population of large grains can increase the scattered intensity at 50-100", this could also be achieved by having ~30% of the dust in a cloud at a distance from us equal to ~95% of the distance to the nova. Such a model also improves agreement with the data at larger angles, and illustrates the sensitivity of X-ray scattering halos to the location of the dust. The observations therefore do not require a population of micron-sized dust grains.

How massive were the first stars? This question is of fundamental importance for galaxy formation and cosmic reionization. Here we consider how protostellar feedback can limit the mass of a forming star. For this we must understand the rate at which primordial protostars accrete, how they and their feedback output evolve, and how this feedback interacts with the infalling matter. We describe the accretion rate with an ``isentropic accretion'' model: the rate is initially very large (~0.03 M_sun/yr when m_* =1 M_sun) and declines as m_*^{-3/7}. Protostellar evolution is treated with a model that tracks the total energy of the star. A key difference compared to previous studies is allowance for rotation of the infalling envelope. This leads to photospheric conditions at the star and dramatic differences in the feedback. Two feedback mechanisms are considered: HII region breakout and radiation pressure from Lyman-alpha and FUV photons. Radiation pressure appears to be the dominant mechanism for suppressing infall, becoming dynamically important around 20 M_sun.

Observations indicate that massive stars form in regions of very high surface density, ~1 g cm^-2. Clusters containing massive stars and globular clusters have a comparable column density. The total pressure in clouds of such a column density is P/k~10^8-10^9 K cm^-3, far greater than that in the diffuse ISM or the average in GMCs. Observations show that massive star-forming regions are supersonically turbulent, and we show that the molecular cores out of which individual massive stars form are as well. The protostellar accretion rate in such a core is approximately equal to the instantaneous mass of the star divided by the free-fall time of the gas that is accreting onto the star (Stahler, Shu, & Taam 1980). The star-formation time in this Turbulent Core model for massive star formation is several mean free-fall timesscales of the core, but is about equal to that of the region in which the core is embedded. The typical time for a massive star to form is about 10^5 yr and the accretion rate is high enough to overcome radiation pressure due to the luminosity of the star. For the typical case we consider, in which the cores out of which the stars form have a density structure varying as r^{-1.5}, the protostellar accretion rate grows linearly with time. We calculate the evolution of the radius of a protostar and determine the accretion luminosity. At the high accretion rates that are typical in regions of massive star formation, protostars join the main sequence at about 20 solar masses. We apply these results to predict the properties of protostars thought to be powering several observed hot molecular cores, including the Orion hot core and W3(H2O). In the Appendixes, we discuss the pressure in molecular clouds and we argue that ``logatropic'' models for molecular clouds are incompatible with observation.

I review the arguments motivating models for massive star formation via stellar collisions. I then describe how the standard accretion scenario, involving the collapse of a quasi-hydrostatic gas core, can produce high-mass stars in the pressurized regions of forming star clusters. I argue that the observational evidence, particularly in the Orion hot core, favors the standard accretion paradigm.

Massive stars (with mass m_* > 8 solar masses) are fundamental to the evolution of galaxies, because they produce heavy elements, inject energy into the interstellar medium, and possibly regulate the star formation rate. The individual star formation time, t_*f, determines the accretion rate of the star; the value of the former quantity is currently uncertain by many orders of magnitude, leading to other astrophysical questions. For example, the variation of t_*f with stellar mass dictates whether massive stars can form simultaneously with low-mass stars in clusters. Here we show that t_*f is determined by conditions in the star's natal cloud, and is typically ~10^5 yr. The corresponding mass accretion rate depends on the pressure within the cloud - which we relate to the gas surface density - and on both the instantaneous and final stellar masses. Characteristic accretion rates are sufficient to overcome radiation pressure from ~100 solar mass protostars, while simultaneously driving intense bipolar gas outflows. The weak dependence of t_*f on the final mass of the star allows high- and low-mass star formation to occur nearly simultaneously in clusters.

We model the formation of high-mass stars, specifying the accretion rate in terms of the instantaneous and final mass of the star, the ambient pressure of the star-forming region and the form of polytropic pressure support of the pre-stellar gas core. The high pressures typical of Galactic regions of massive star formation allow a 100 solar mass star to form in ~10^5 years with a final accretion rate ~10^-3 solar masses per year. By modeling protostellar evolution we predict the properties of several nearby massive protostars. We model cluster formation by applying this theory to many stars. We use the observed intensity of outflows from protoclusters to estimate the star formation rate, finding that clusters take at least several free-fall times to form; for a cluster similar to the Orion Nebula Cluster, we predict a formation timescale of ~1x10^6 years.

Supernova (SN) 1998bw and gamma-ray burst (GRB) 980425 offer the first direct evidence that supernovae are the progenitors of some GRBs. However, this burst was unusually dim, smooth and soft compared to other bursts with known afterglows. Whether it should be considered a prototype for cosmological GRBs depends largely on whether the supernova explosion and burst were asymmetrical or can be modeled as spherical. We address this question by treating the acceleration of the supernova shock in the outermost layers of the stellar envelope, the transition to relativistic flow, and the subsequent expansion (and further acceleration) of the ejecta into the surrounding medium. We find that GRB 980425 could plausibly have been produced by a collision between the relativistic ejecta from SN 1998bw and the star's pre-supernova wind; the model requires no significant asymmetry. This event therefore belongs to a dim subclass of GRBs and is not a prototype for jet-like cosmological GRBs.

The spatiotemporal coincidence of supernova (SN) 1998bw and gamma-ray burst (GRB) 980425 and this supernova's unusual optical and radio properties have prompted many theoretical models that produce GRBs from supernovae. We review the salient features of our simple, spherical model in which an energetic supernova explosion shock accelerates a small fraction of the progenitor's stellar envelope to mildly relativistic velocities. This material carries sufficient energy to produce a weak GRB and a bright radio supernova through an external shock against a dense stellar wind.

We investigate the acceleration of shock waves to relativistic velocities in the outer layers of exploding stars. By concentrating the explosion energy in the outermost ejecta, such trans-relativistic blast waves can serve as the progenitors of gamma-ray bursts (GRBs); in particular, the ``baryon-loading'' problem that plagues many models of GRBs is circumvented. We present physically motivated and numerically validated analytic expressions to describe trans-relativistic blast waves in supernovae. We find that relativistic ejecta are enhanced in more centrally condensed envelopes, e.g., for radiative envelopes, when the luminosity approaches the Eddington limit. We present convenient formulae for estimating the relativistic ejecta from a given progenitor. We apply our analytic and numerical methods to a model of SN 1998bw, finding significantly enhanced relativistic ejecta compared to previous studies. We propose that GRB 980425 is associated with SN 1998bw and may have resulted from an approximately spherical explosion producing ~10^-6 M_sun of mildly relativistic ejecta with mean Lorentz factor ~2, which then interacted with a dense circumstellar wind with mass loss rate ~few x 10^-4 M_sun/yr. A highly asymmetric explosion is not required. An extreme model of ``hypernova'' explosions in massive stars is able to account for the energetics and relativistic ejecta velocities required by many of the observed cosmological GRBs. We present simplified models of explosions resulting from accretion-induced collapse of white dwarfs and phase transitions of neutron stars. While we find increased energies in relativistic ejecta compared to previous studies, these explosions are unlikely to be observed at cosmological distances with current detectors. (abridged)

We examine the effects of ionization, radiation pressure and main sequence winds from massive stars on self-gravitating, clumpy molecular clouds, thereby modeling the formation and pre-supernova feedback of massive star clusters. We find the process of ``turbulent mass loading'' is effective in confining HII regions. Extrapolating typical Galactic high-mass star forming regions to greater initial gas cloud masses and assuming steady star formation rates, we determine the timescales for cloud disruption. We find that a dense n_c ~ 2 x 10^5 cm^-3 cloud with initial mass M_c ~ 4 x 10^5 M_sun is able to form ~2 x 10^5 M_sun of stars (50% efficiency) before feedback disperses the gas after ~3 Myr. This mass and age are typical of young, optically visible super star clusters (SSCs). The high efficiency permits the creation of a bound stellar system.

We invoke star formation triggered by cloud-cloud collisions to explain global star formation rates of disk galaxies and circumnuclear starbursts. Previous theories based on the growth rate of gravitational perturbations ignore the dynamically important presence of magnetic fields. Theories based on triggering by spiral density waves fail to explain star formation in systems without such waves. Furthermore, observations suggest gas and stellar disk instabilities are decoupled. Following Gammie, Ostriker & Jog (1991), the cloud collision rate is set by the shear velocity of encounters with initial impact parameters of a few tidal radii, due to differential rotation in the disk. This, together with the effective confinement of cloud orbits to a two dimensional plane, enhances the collision rate above that for particles in a three dimensional box. We predict Sigma_{SFR}(R) proportional to Sigma_{gas} Omega (1-0.7 beta). For constant circular velocity (beta = 0), this is in agreement with recent observations (Kennicutt 1998). We predict a B-band Tully-Fisher relation: L_{B} proportional to v_{circ}^{7/3}, also consistent with observations. As additional tests, we predict enhanced star formation in regions with relatively high shear rates, and lower star formation efficiencies in clouds of higher mass.

We consider the effect of lensing by galaxy clusters on the angular distribution of submillimeter wavelength objects. While lensing does not change the total flux and number counts of submillimeter sources, it can affect the number counts and fluxes of flux-limited samples. Therefore imposing a flux cut on point sources not only reduces the overall Poisson noise, but imprints the correlations between lensing clusters on the unresolved flux distribution. Using a simple model, we quantify the lensing anisotropy induced in flux-limited samples and compare this to Poisson noise. We find that while the level of induced anisotropies on the scale of the cluster angular correlation length is comparable to Poisson noise for a slowly evolving cluster model, it is negligible for more realistic models of cluster evolution. Thus the removal of point sources is not expected to induce measurable structure in the microwave or far-infrared backgrounds.

We present a simple model of the infra-red universe, based as much as possible on local observations. We model the luminosity and number evolution of disk and starburst galaxies, including the effects of dust, gas and spectral evolution. Although simple, our approach is able to reproduce observations of galaxy number counts and the infra-red and sub-millimeter extra-galactic backgrounds. It provides a useful probe of galaxy formation and evolution out to high redshift. The model demonstrates the significant role of the starburst population and predicts high star formation rates at redshifts 3 to 4, consistent with recent extinction-corrected observations of Lyman break galaxies. Starbursting galaxies are predicted to dominate the current SCUBA surveys. Their star formation is driven predominantly by strong tidal interactions and mergers of galaxies. This leads to the creation of spheroidal stellar systems, which may act as the seeds for disk formation as gas infalls. We predict the present-day baryonic mass in bulges and halos is comparable to that in disks. From observations of the extra-galactic background, the model predicts that the vast majority of star formation in the Universe occurs at z<5.