While magnetic fields are important in contemporary star formation, their role in primordial star formation is unknown. Magnetic fields of order 10^-16 G are produced by the Biermann battery due to the curved shocks and turbulence associated with the infall of gas into the dark matter minihalos that are the sites of formation of the first stars. These fields are rapidly amplified by a small-scale dynamo until they saturate at or near equipartition with the turbulence in the central region of the gas. Analytic results are given for the outcome of the dynamo, including the effect of compression in the collapsing gas. The mass-to-flux ratio in this gas is 2-3 times the critical value, comparable to that in contemporary star formation. Predictions of the outcomes of simulations using smooth particle hydrodynamics (SPH) and grid-based adaptive mesh refinement (AMR) are given. Because the numerical viscosity and resistivity for the standard resolution of 64 cells per Jeans length are several orders of magnitude greater than the physical values, dynamically significant magnetic fields affect a much smaller fraction of the mass in simulations than in reality. An appendix gives an analytic treatment of free-fall collapse, including that in a constant density background. Another appendix presents a new method of estimating the numerical viscosity; results are given for both SPH and grid-based codes.

We construct a new analytic phenomenological model for the extended circumgalactic material (CGM) of $L^*$ galaxies. Our model reproduces the OVII/OVIII absorption observations of the Milky Way (MW) and the OVI measurements reported by the COS-Halos and eCGM surveys. The warm/hot gas is in hydrostatic equilibrium in a MW gravitational potential, and we adopt a barotropic equation of state, resulting in a temperature variation as a function of radius. A pressure component with an adiabatic index of $\gamma=4/3$ is included to approximate the effects of a magnetic field and cosmic rays. We introduce a metallicity gradient motivated by the enrichment of the inner CGM by the Galaxy. We then present our fiducial model for the corona, tuned to reproduce the observed OVI-OVIII column densities, and with a total mass of $M_{\rm gas} \approx 5.5 \times 10^{10}~{\rm M_{\odot}}$ inside $r_{\rm cgm} \approx 280$ kpc. The gas densities in the CGM are low ($n_{\rm H} = 10^{-5} - 3 \times 10^{-4}~{\rm cm^{-3}}$) and its collisional ionization state is modified by the metagalactic radiation field (MGRF). We show that for OVI-bearing warm/hot gas with typical observed column densities $N_{\rm OVI} \sim 3 \times 10^{14}~{\rm cm^{-2}}$ at large ($\gtrsim 100$ kpc) impact parameters from the central galaxies, the ratio of the cooling to dynamical times, $t_{\rm cool}/t_{\rm dyn}$, has a model-independent upper limit of $\lesssim 4$. In our model, $t_{\rm cool}/t_{\rm dyn}$ at large radii is $\sim 2-3$. We present predictions for a wide range of future observations of the warm/hot CGM, from UV/X-ray absorption and emission spectroscopy, to dispersion measure (DM) and Sunyaev-Zeldovich CMB measurements. We provide the model outputs in machine-readable data files, for easy comparison and analysis.

Gravitationally-bound clusters that survive gas removal represent an unusual mode of star formation in the Milky Way and similar spiral galaxies. While forming, they can be distinguished observationally from unbound star formation by their high densities, virialised velocity structures, and star formation histories that accelerate toward the present, but extend multiple free-fall times into the past. In this paper we examine several proposed scenarios for how such structures might form and evolve, and carry out a Bayesian analysis to test these models against observed distributions of protostellar age, counts of young stellar objects relative to gas, and the overall star formation rate of the Milky Way. We show that models in which the acceleration of star formation is due either to a large-scale collapse or a time-dependent increase in star formation efficiency are unable to satisfy the combined set of observational constraints. In contrast, models in which clusters form in a "conveyor belt" mode where gas accretion and star formation occur simultaneously, but the star formation rate per free-fall time is low, can match the observations.

Star clusters stand at the intersection of much of modern astrophysics: the interstellar medium, gravitational dynamics, stellar evolution, and cosmology. Here we review observations and theoretical models for the formation, evolution, and eventual disruption of star clusters. Current literature suggests a picture of this life cycle with several phases: (1) Clusters form in hierarchically-structured, accreting molecular clouds that convert gas into stars at a low rate per dynamical time until feedback disperses the gas. (2) The densest parts of the hierarchy resist gas removal long enough to reach high star formation efficiency, becoming dynamically-relaxed and well-mixed. These remain bound after gas removal. (3) In the first $\sim 100$ Myr after gas removal, clusters disperse moderately fast, through a combination of mass loss and tidal shocks by dense molecular structures in the star-forming environment. (4) After $\sim 100$ Myr, clusters lose mass via two-body relaxation and shocks by giant molecular clouds, processes that preferentially affect low-mass clusters and cause a turnover in the cluster mass function to appear on $\sim 1-10$ Gyr timescales. (5) Even after dispersal, some clusters remain coherent and thus detectable in chemical or action space for multiple galactic orbits. In the next decade a new generation of space- and AO-assisted ground-based telescopes will enable us to test and refine this picture.

For dark matter (DM) particles with masses in the 0.6 - 6 m_p range, we set stringent constraints on the interaction cross-sections for scattering with ordinary baryonic matter. These constraints follow from the recognition that such particles can be captured by - and thermalized within - the Earth, leading to a substantial accumulation and concentration of DM that interact with baryons. Here, we discuss the probability that DM intercepted by the Earth will be captured, the number of DM particles thereby accumulated over Earth's lifetime, the fraction of such particles retained in the face of evaporation, and the density distribution of such particles within the Earth. In the latter context, we note that a previous treatment of the density distribution of DM, presented by Gould and Raffelt and applied subsequently to DM in the Sun, is inconsistent with considerations of hydrostatic equilibrium. Our analysis provides an estimate of the DM particle density at Earth's surface, which may exceed 1.E+14 cm-3 for the mass range under consideration. Based upon our determination of the DM density at Earth's surface, we derive constraints on the scattering cross-sections. These constraints are placed by four considerations: (1) the lifetime of the relativistic proton beam at the Large Hadron collider (LHC); (2) the orbital decay of spacecraft in low Earth orbit (LEO); (3) the vaporization rate of cryogenic liquids in well-insulated storage dewars; and (4) the thermal conductivity of Earth's crust. As an example application of our results, we show that for the scattering cross-sections that were invoked recently in Barkana's original explanation for the anomalously deep 21 cm absorption reported by EDGES, DM particle masses in the 0.6 - 4 m_p range are ruled out.

Many spiral galaxies have extended outer H~I disks and display low levels of star formation, inferred from the far-ultraviolet emission detected by {\it GALEX}, well beyond the optical radius. Here, we investigate the supernova (SN) rate in the outskirts of galaxies, using the largest and most homogeneous set of nearby supernovae (SNe) from the Lick Observatory Supernova Search (LOSS). While SN rates have been measured with respect to various galaxy properties, such as stellar mass and metallicity, their relative frequency in the outskirts versus the inner regions has not yet been studied. Understanding the SN rate as a function of intragalactic environment has many ramifications, including the interpretation of LIGO observations, the formation of massive stars, and the puzzlingly high velocity dispersion of the outer H~I disk. Using data from the LOSS survey, we find that the rate beyond the optical radius of spiral galaxies is $2.5 \pm 0.5$ SNe per millennium, while dwarf galaxies host $ 4.0 \pm 2.2$ SNe per millennium. The rates of core-collapse SNe (that may collapse to form the massive black holes detected by LIGO/Virgo) in the outer disks of spirals is $1.5 \pm 0.15$ SNe per millennium and in dwarf galaxies is $2.6 \pm 1.5$ SNe per millennium. Core-collapse SNe in spiral outskirts contribute $7600 \pm 1700$\,SNe\,~Gpc$^{-3}$\,yr$^{-1}$ to the volumetric rate, and dwarf galaxies have a rate of $31,000 \pm 18,000$\,SNe\,Gpc$^{-3}$\,yr$^{-1}$. The relative ratio of core-collapse to Type Ia SNe is comparable in the inner and outer parts of spirals, and in dwarf galaxies.

Recent papers have found that the inferred slope of the high-mass ($>1.5$ M$_\odot$) IMF for field stars in the solar vicinity has a larger value ($\sim 1.7-2.1$) than the slopes ($\sim 1.2-1.7$; Salpeter= 1.35) inferred from numerous studies of young clusters. We attempt to reconcile this apparent contradiction. Stars mostly form in Giant Molecular Clouds, and the more massive stars ($\gtrsim 3$ M$_\odot$) may have insufficient time before their deaths to uniformly populate the solar circle of the Galaxy. We examine the effect of small sample volumes on the {\it apparent} slope, $\Gamma_{\rm app}$, of the high-mass IMF by modeling the present day mass function (PDMF) over the mass range $1.5-6$ M$_\odot$. Depending on the location of the observer along the solar circle and the size of the sample volume, the apparent slope of the IMF can show a wide variance, with typical values steeper than the underlying universal value $\Gamma$. We show, for example, that the PDMFs observed in a small (radius $\sim 200$ pc) volume randomly placed at the solar circle have a $\sim 15-30$\% likelihood of resulting in $\Gamma_{\rm app} \gtrsim \Gamma+ 0.35$ because of inhomogeneities in the surface densities of more massive stars. If we add the a priori knowledge that the Sun currently lies in an interarm region, where the star formation rate is lower than the average at the solar circle, we find an even higher likelihood ($\sim 50-60\%$ ) of $\Gamma_{\rm app} \gtrsim \Gamma+0.35$, corresponding to $\Gamma_{\rm app} \gtrsim 1.7$ when the underlying $\Gamma= 1.35$.

We present a large suite of simulations of the formation of low-mass star clusters. Our simulations include an extensive set of physical processes -- magnetohydrodynamics, radiative transfer, and protostellar outflows -- and span a wide range of virial parameters and magnetic field strengths. Comparing the outcomes of our simulations to observations, we find that simulations remaining close to virial balance throughout their history produce star formation efficiencies and initial mass function (IMF) peaks that are stable in time and in reasonable agreement with observations. Our results indicate that small-scale dissipation effects near the protostellar surface provide a feedback loop for stabilizing the star formation efficiency. This is true regardless of whether the balance is maintained by input of energy from large scale forcing or by strong magnetic fields that inhibit collapse. In contrast, simulations that leave virial balance and undergo runaway collapse form stars too efficiently and produce an IMF that becomes increasingly top-heavy with time. In all cases we find that the competition between magnetic flux advection toward the protostar and outward advection due to magnetic interchange instabilities, and the competition between turbulent amplification and reconnection close to newly-formed protostars renders the local magnetic field structure insensitive to the strength of the large-scale field, ensuring that radiation is always more important than magnetic support in setting the fragmentation scale and thus the IMF peak mass. The statistics of multiple stellar systems are similarly insensitive to variations in the initial conditions and generally agree with observations within the range of statistical uncertainty.

Star forming molecular clouds are observed to be both highly magnetized and turbulent. Consequently the formation of protostellar disks is largely dependent on the complex interaction between gravity, magnetic fields, and turbulence. Studies of non-turbulent protostellar disk formation with realistic magnetic fields have shown that these fields are efficient in removing angular momentum from the forming disks, preventing their formation. However, once turbulence is included, disks can form in even highly magnetized clouds, although the precise mechanism remains uncertain. Here we present several high resolution simulations of turbulent, realistically magnetized, high-mass molecular clouds with both aligned and random turbulence to study the role that turbulence, misalignment, and magnetic fields have on the formation of protostellar disks. We find that when the turbulence is artificially aligned so that the angular momentum is parallel to the initial uniform field, no rotationally supported disks are formed, regardless of the initial turbulent energy. We conclude that turbulence and the associated misalignment between the angular momentum and the magnetic field are crucial in the formation of protostellar disks in the presence of realistic magnetic fields.

Star formation in a filamentary infrared dark cloud (IRDC) is simulated over a dynamic range of 4.2 pc to 28 au for a period of $3.5\times 10^5$ yr, including magnetic fields and both radiative and outflow feedback from the protostars. At the end of the simulation, the star formation efficiency is 4.3 per cent and the star formation rate per free fall time is $\epsilon_{\rm ff}\simeq 0.04$, within the range of observed values (Krumholz et al. 2012a). The total stellar mass increases as $\sim\,t^2$, whereas the number of protostars increases as $\sim\,t^{1.5}$. We find that the density profile around most of the simulated protostars is $\sim\,\rho\propto r^{-1.5}$, as predicted by Murray & Chang (2015). At the end of the simulation, the protostellar mass function approaches the Chabrier (2005) stellar initial mass function. We infer that the time to form a star of median mass $0.2\,M_\odot$ is about $1.4\times 10^5$~yr from the median mass accretion rate. We find good agreement among the protostellar luminosities observed in the large sample of Dunham et al. (2013), our simulation, and a theoretical estimate, and conclude that the classical protostellar luminosity problem Kenyon et al. (1990) is resolved. The multiplicity of the stellar systems in the simulation agrees to within a factor 2 of observations of Class I young stellar objects; most of the simulated multiple systems are unbound. Bipolar protostellar outflows are launched using a sub-grid model, and extend up to 1 pc from their host star. The mass-velocity relation of the simulated outflows is consistent with both observation and theory.

Star formation in our Galaxy occurs in molecular clouds that are self-gravitating, highly turbulent, and magnetized. We study the conditions under which cloud cores inherit large-scale magnetic field morphologies and how the field is governed by cloud turbulence. We present four moving-mesh simulations of supersonic, turbulent, isothermal, self-gravitating gas with a range of magnetic mean-field strengths characterized by the Alfv\'enic Mach number $\mathcal{M}_{{\rm A}, 0}$, resolving pre-stellar core formation from parsec to a few AU scales. In our simulations with the turbulent kinetic energy density dominating over magnetic pressure ($\mathcal{M}_{{\rm A}, 0}>1$), we find that the collapse is approximately isotropic with $B\propto\rho^{2/3}$, core properties are similar regardless of initial mean-field strength, and the field direction on $100$ AU scales is uncorrelated with the mean field. However, in the case of a dominant large-scale magnetic field ($\mathcal{M}_{{\rm A}, 0}=0.35$), the collapse is anisotropic with $B\propto\rho^{1/2}$. This transition at $\mathcal{M}_{{\rm A}, 0}\sim1$ is not expected to be sharp, but clearly signifies two different paths for magnetic field evolution in star formation. Based on observations of different star forming regions, we conclude that star formation in the interstellar medium may occur in both regimes. Magnetic field correlation with the mean-field extends to smaller scales as $\mathcal{M}_{{\rm A}, 0}$ decreases, making future ALMA observations useful for constraining $\mathcal{M}_{{\rm A}, 0}$ of the interstellar medium.

We construct an analytic phenomenological model for extended warm/hot gaseous coronae of $L_*$ galaxies. We consider UV OVI COS-Halos absorption line data in combination with Milky Way X-ray OVII and OVIII absorption and emission. We fit these data with a single model representing the COS-Halos galaxies and a Galactic corona. Our model is multi-phased, with hot and warm gas components, each with a (turbulent) log-normal distribution of temperatures and densities. The hot gas, traced by the X-ray absorption and emission, is in hydrostatic equilibrium in a Milky Way gravitational potential. The median temperature of the hot gas is $1.5 \times 10^6$~K and the mean hydrogen density is $\sim 5 \times 10^{-5}~{\rm cm^{-3}}$. The warm component as traced by the OVI, is gas that has cooled out of the high density tail of the hot component. The total warm/hot gas mass is high and is $1.2 \times 10^{11}~{\rm M_{\odot}}$. The gas metallicity we require to reproduce the oxygen ion column densities is $0.5$ solar. The warm OVI component has a short cooling time ($\sim 2 \times 10^8$ years), as hinted by observations. The hot component, however, is $\sim 80\%$ of the total gas mass and is relatively long-lived, with $t_{cool} \sim 7 \times 10^{9}$ years. Our model supports suggestions that hot galactic coronae can contain significant amounts of gas. These reservoirs may enable galaxies to continue forming stars steadily for long periods of time and account for "missing baryons" in galaxies in the local universe.

The pressure exerted by massive stars' radiation fields is an important mechanism regulating their formation. Detailed simulation of massive star formation therefore requires an accurate treatment of radiation. However, all published simulations have either used a diffusion approximation of limited validity; have only been able to simulate a single star fixed in space, thereby suppressing potentially-important instabilities; or did not provide adequate resolution at locations where instabilities may develop. To remedy this we have developed a new, highly accurate radiation algorithm that properly treats the absorption of the direct radiation field from stars and the re-emission and processing by interstellar dust. We use our new tool to perform three-dimensional radiation-hydrodynamic simulations of the collapse of massive pre-stellar cores with laminar and turbulent initial conditions and properly resolve regions where we expect instabilities to grow. We find that mass is channeled to the stellar system via gravitational and Rayleigh-Taylor (RT) instabilities, in agreement with previous results using stars capable of moving, but in disagreement with methods where the star is held fixed or with simulations that do not adequately resolve the development of RT instabilities. For laminar initial conditions, proper treatment of the direct radiation field produces later onset of instability, but does not suppress it entirely provided the edges of radiation-dominated bubbles are adequately resolved. Instabilities arise immediately for turbulent pre-stellar cores because the initial turbulence seeds the instabilities. Our results suggest that RT features are significant and should be present around accreting massive stars throughout their formation.

To understand the conditions under which dense, molecular gas is able to form within a galaxy, we post-process a series of three-dimensional galactic-disk-scale simulations with ray-tracing based radiative transfer and chemical network integration to compute the equilibrium chemical and thermal state of the gas. In performing these simulations we vary a number of parameters, such as the ISRF strength, vertical scale height of stellar sources, cosmic ray flux, to gauge the sensitivity of our results to these variations. Self-shielding permits significant molecular hydrogen (H2) abundances in dense filaments around the disk midplane, accounting for approximately ~10-15% of the total gas mass. Significant CO fractions only form in the densest, n>~10^3 cm^-3, gas where a combination of dust, H2, and self-shielding attenuate the FUV background. We additionally compare these ray-tracing based solutions to photochemistry with complementary models where photo-shielding is accounted for with locally computed prescriptions. With some exceptions, these local models for the radiative shielding length perform reasonably well at reproducing the distribution and amount of molecular gas as compared with a detailed, global ray tracing calculation. Specifically, an approach based on the Jeans Length with a T=40K temperature cap performs the best in regards to a number of different quantitative measures based on the H2 and CO abundances.

As star-forming clouds collapse, the gas within them fragments to ever-smaller masses. Naively one might expect this process to continue down to the smallest mass that is able to radiate away its binding energy on a dynamical timescale, the opacity limit for fragmentation, at $\sim 0.01$ $M_\odot$. However, the observed peak of the initial mass function (IMF) lies a factor of $20-30$ higher in mass, suggesting that some other mechanism halts fragmentation before the opacity limit is reached. In this paper we analyse radiation-magnetohydrodynamic simulations of star cluster formation in typical Milky Way environments in order to determine what physical process limits fragmentation in them. We examine the regions in the vicinity of stars that form in the simulations to determine the amounts of mass that are prevented from fragmenting by thermal and magnetic pressure. We show that, on small scales, thermal pressure enhanced by stellar radiation heating is the dominant mechanism limiting the ability of the gas to further fragment. In the brown dwarf mass regime, $\sim 0.01$ $M_\odot$, the typical object that forms in the simulations is surrounded by gas whose mass is several times its own that is unable to escape or fragment, and instead is likely to accrete. This mechanism explains why $\sim 0.01$ $M_\odot$ objects are rare: unless an outside agent intervenes (e.g., a shock strips away the gas around them), they will grow by accreting the warmed gas around them. In contrast, by the time stars grow to masses of $\sim 0.2$ $M_\odot$, the mass of heated gas is only tens of percent of the central star mass, too small to alter its final mass by a large factor. This naturally explains why the IMF peak is at $\sim 0.2$ $M_\odot$.

The surface density and vertical distribution of stars, stellar remnants, and gas in the solar vicinity form important ingredients for understanding the star formation history of the Galaxy as well as for inferring the local density of dark matter by using stellar kinematics to probe the gravitational potential. In this paper we review the literature for these baryonic components, reanalyze data, and provide tables of the surface densities and exponential scale heights of main sequence stars, giants, brown dwarfs, and stellar remnants. We also review three components of gas (H2, HI, and HII), give their surface densities at the solar circle, and discuss their vertical distribution. We find a local total surface density of M dwarfs of 17.3 pm 2.3 Mo/pc^2. Our result for the total local surface density of visible stars, 27.0 pm 2.7 Mo/pc^2, is close to previous estimates due to a cancellation of opposing effects: more mass in M dwarfs, less mass in the others. The total local surface density in white dwarfs is 4.9 pm 0.6 Mo/pc^2; in brown dwarfs, it is ~1.2 Mo/pc^2. We find that the total local surface density of stars and stellar remnants is 33.4 pm 3 Mo/pc^2, somewhat less than previous estimates. We analyze data on 21 cm emission and absorption and obtain good agreement with recent results on the local amount of neutral atomic hydrogen obtained with the Planck satellite. The local surface density of gas is 13.7 pm 1.6 Mo/pc^2. The total baryonic mass surface density that we derive for the solar neighborhood is 47.1 pm 3.4 Mo/pc^2. Combining these results with others' measurements of the total surface density of matter within 1-1.1 kpc of the plane, we find that the local density of dark matter is 0.013 pm 0.003Mo/pc^3.The local density of all matter is 0.097 pm 0.013 Mo/pc^3. We discuss limitations on the properties of a possible thin disk of dark matter.

The intermittent dissipation of interstellar turbulence is an important energy source in the diffuse ISM. Though on average smaller than the heating rates due to cosmic rays and the photoelectric effect on dust grains, the turbulent cascade can channel large amounts of energy into a relatively small fraction of the gas that consequently undergoes significant heating and chemical enrichment. In particular, this mechanism has been proposed as a solution to the long-standing problem of the high abundance of CH+ along diffuse molecular sight lines, which steady-state, low temperature models under-produce by over an order of magnitude. While much work has been done on the structure and chemistry of these small-scale dissipation zones, comparatively little attention has been paid to relating these zones to the properties of the large-scale turbulence. In this paper, we attempt to bridge this gap by estimating the temperature and CH+ column density along diffuse molecular sight-lines by post-processing 3-dimensional MHD turbulence simulations. Assuming reasonable values for the cloud density (30 / cm^3), size (20 pc), and velocity dispersion (2.3 km / s), we find that our computed abundances compare well with CH+ column density observations, as well as with observations of emission lines from rotationally excited H2 molecules.

The most accurate measurements of magnetic fields in star-forming gas are based on the Zeeman observations analyzed by Crutcher et al. (2010). We show that their finding that the 3D magnetic field scales approximately as density$^{0.65}$ can also be obtained from analysis of the observed line-of-sight fields. We present two large-scale AMR MHD simulations of several thousand $M_\odot$ of turbulent, isothermal, self-gravitating gas, one with a strong initial magnetic field (Alfven Mach number $M_{A,0}= 1$) and one with a weak initial field ($M_{A,0}=10$). We construct samples of the 100 most massive clumps in each simulation and show that they exhibit a power-law relation between field strength and density in excellent agreement with the observed one. Our results imply that the average field in molecular clumps in the interstellar medium is $<B_{tot}> \sim 42 n_{H,4}^{0.65} \mu$G. Furthermore, the median value of the ratio of the line-of-sight field to density$^{0.65}$ in the simulations is within a factor of about (1.3, 1.7) of the observed value for the strong and weak field cases, respectively. The median value of the mass-to-flux ratio, normalized to the critical value, is 70% of the line-of-sight value. This is larger than the 50% usually cited for spherical clouds because the actual mass-to-flux ratio depends on the volume-weighted field, whereas the observed one depends on the mass-weighted field. Our results indicate that the typical molecular clump in the ISM is significantly supercritical (~ factor of 3). The results of our strong-field model are in very good quantitative agreement with the observations of Li et al. (2009), which show a strong correlation in field orientation between small and large scales. Because there is a negligible correlation in the weak-field model, we conclude that molecular clouds form from strongly magnetized (although magnetically supercritical) gas.

The turbulent environment from which stars form may lead to misalignment between the stellar spin and the remnant protoplanetary disk. By using hydrodynamic and magnetohydrodynamic simulations, we demonstrate that a wide range of stellar obliquities may be produced as a by-product of forming a star within a turbulent environment. We present a simple semi-analytic model that reveals this connection between the turbulent motions and the orientation of a star and its disk. Our results are consistent with the observed obliquity distribution of hot Jupiters. Migration of misaligned hot Jupiters may, therefore, be due to tidal dissipation in the disk, rather than tidal dissipation of the star-planet interaction.

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.

In regions of star formation, protostars and newborn stars accrete mass from their natal clouds. These clouds are threaded by magnetic fields with a strength characterized by the plasma beta---the ratio of thermal and magnetic pressures. Observations show molecular clouds have beta <= 1, so magnetic fields can play a significant role in the accretion process. We have carried out a numerical study of the effect of large-scale magnetic fields on the rate of accretion onto a uniformly moving point particle from a uniform, non-self-gravitating, isothermal gas. We consider gas moving with sonic Mach numbers of up M ~ 45, magnetic fields that are either parallel, perpendicular, or oriented 45 degrees to the flow, and beta as low as 0.01. Our simulations utilize AMR to obtain high spatial resolution where needed; this also allows the simulation boundaries to be far from the accreting object. Additionally, we show our results are independent of our exact prescription for accreting mass in the sink particle. We give simple expressions for the steady-state accretion rate as a function of beta, M, and field orientation. Using typical molecular clouds values of M ~ 5 and beta ~ 0.04 from the literature, our fits suggest a 0.4 M_Sun star accretes ~ 4*10^{-9} M_Sun/year, almost a factor of two less than accretion rates predicted by hydrodynamic models. This disparity grows to orders of magnitude for stronger fields and lower Mach numbers. We discuss the applicability of these accretion rates versus accretion rates expected from gravitational collapse, and when a steady state is possible. This reduction in Mdot increases the time required to form stars in competitive accretion models, making such models less efficient. In numerical codes, our results should enable accurate subgrid models of sink particles accreting from magnetized media.

We present three Orion simulations of star cluster formation in a 1000 Msun, turbulent molecular cloud clump, including the effects of radiative transfer, protostellar outflows, and magnetic fields. Our simulations all use self-consistent turbulent initial conditions and vary the mean mass-to-flux ratio relative to the critical value over 2, 10, and infinity to gauge the influence of magnetic fields on star cluster formation. We find, in good agreement with previous studies, that magnetic fields of typically observed strengths lower the star formation rate by a factor of 2.4 and reduce the amount of fragmentation by a factor of 2 relative to the zero-field case. We also find that the field increases the characteristic sink particle mass, again by a factor of 2.4. The magnetic field also increases the degree of clustering in our simulations, such that the maximum stellar densities in the strong field case are higher than the others by again a factor of 2. This clustering tends to encourage the formation of multiple systems, which are more common in the rad-MHD runs than the rad-hydro run. The companion frequency in our simulations is consistent with observations of multiplicity in Class I sources, particularly for the strong field case. Finally, we find evidence of primordial mass segregation in our simulations reminiscent of that observed in star clusters like the Orion Nebula Cluster.

We review recent advances in the analytical and numerical modeling of the star formation rate in molecular clouds and discuss the available observational constraints. We focus on molecular clouds as the fundamental star formation sites, rather than on the larger-scale processes that form the clouds and set their properties. Molecular clouds are shaped into a complex filamentary structure by supersonic turbulence, with only a small fraction of the cloud mass channeled into collapsing protostars over a free-fall time of the system. In recent years, the physics of supersonic turbulence has been widely explored with computer simulations, leading to statistical models of this fragmentation process, and to the prediction of the star formation rate as a function of fundamental physical parameters of molecular clouds, such as the virial parameter, the rms Mach number, the compressive fraction of the turbulence driver, and the ratio of gas to magnetic pressure. Infrared space telescopes, as well as ground-based observatories have provided unprecedented probes of the filamentary structure of molecular clouds and the location of forming stars within them.

We present dark-matter minihalo models for the Ultra-Compact High Velocity HI Clouds (UCHVCs) recently discovered in the 21 cm ALFALFA survey. We assume gravitational confinement of 10^4 K HI gas by flat-cored dark-matter subhalos within the Local Group. We show that for flat cores, typical (median) tidally-stripped cosmological subhalos at redshift z=0 have dark-matter masses of ~10^7 M_{sun} within the central 300 pc (independent of total halo mass), consistent with the "Strigari mass scale" observed in low-luminosity dwarf galaxies. Flat-cored subhalos also resolve the mass-discrepancy between simulated and observed satellites around the Milky Way. For the UCHVCs we calculate the photoionization-limited hydrostatic gas profiles for any distance-dependent total observed HI mass and predict the associated (projected) HI half-mass radii, assuming the clouds are embedded in distant (d > 300 kpc) and unstripped subhalos. For a typical UCHVC (0.9 Jy km/s) we predict physical HI half-mass radii of 0.18 to 0.35 kpc (or angular sizes of 0.6 to 2.1 arcmin) for distances ranging from 300 kpc to 2 Mpc. As a consistency check we model the gas rich dwarf galaxy Leo T, for which there is a well-resolved HI column density profile and a known distance (420 kpc). For Leo T we find that a subhalo with M_{300} ~ 8 (+/-0.2) 10^6 M_{sun} best fits the observed HI profile. We derive an upper limit of P_{HIM} < 150 K/cm^3 for the pressure of any enveloping hot IGM gas at the distance of Leo T. Our analysis suggests that some of the UCHVCs may in fact constitute a population of 21-cm-selected but optically-faint dwarf galaxies in the Local Group.

We present a model in which the 22 GHz H$_2$O masers observed in star-forming regions occur behind shocks propagating in dense regions (preshock density $n_0 \sim 10^6 - 10^8$ cm$^{-3}$). We focus on high-velocity ($v_s > 30$ km s$^{-1}$) dissociative J shocks in which the heat of H$_2$ re-formation maintains a large column of $\sim 300-400$ K gas; at these temperatures the chemistry drives a considerable fraction of the oxygen not in CO to form H$_2$O. The H$_2$O column densities, the hydrogen densities, and the warm temperatures produced by these shocks are sufficiently high to enable powerful maser action. The observed brightness temperatures (generally $\sim 10^{11} - 10^{14}$ K) are the result of coherent velocity regions that have dimensions in the shock plane that are 10 to 100 times the shock thickness of $\sim 10^{13}$ cm. The masers are therefore beamed towards the observer, who typically views the shock "edge-on", or perpendicular to the shock velocity; the brightest masers are then observed with the lowest line of sight velocities with respect to the ambient gas. We present numerical and analytic studies of the dependence of the maser inversion, the resultant brightness temperature, the maser spot size and shape, the isotropic luminosity, and the maser region magnetic field on the shock parameters and the coherence path length; the overall result is that in galactic H$_2$O 22 GHz masers these observed parameters can be produced in J shocks with $n_0\sim 10^6 - 10^8$ cm$^{-3}$ and $v_s \sim 30 -200$ km s$^{-1}$. A number of key observables such as maser shape, brightness temperature, and global isotropic luminosity depend only on the particle flux into the shock, $j=n_0v_s$, rather than on $n_0$ and $v_s$ separately.

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.

We present a set of 3-dimensional, radiation-magnetohydrodynamic calculations of the gravitational collapse of massive (300 Msun), star-forming molecular cloud cores. We show that the combined effects of magnetic fields and radiative feedback strongly suppress core fragmentation, leading to the production of single star systems rather than small clusters. We find that the two processes are efficient at suppressing fragmentation in different regimes, with the feedback most effective in the dense, central region and the magnetic field most effective in more diffuse, outer regions. Thus, the combination of the two is much more effective at suppressing fragmentation than either one considered in isolation. Our work suggests that typical massive cores, which have mass-to-flux ratios of about 2 relative to critical, likely form a single star system, but that cores with weaker fields may form a small star cluster. This result helps us understand why the observed relationship between the core mass function and the stellar initial mass function holds even for ~100 Msun cores with many thermal Jeans masses of material. We also demonstrate that a ~40 AU Keplerian disk is able to form in our simulations, despite the braking effect caused by the strong magnetic field.

(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.

The temperature of the gas in molecular clouds is a key determinant of the characteristic mass of star formation. Ambipolar diffusion (AD) is considered one of the most important heating mechanisms in weakly ionized molecular clouds. In this work, we study the AD heating rate using 2-fluid turbulence simulations and compare it with the overall heating rate due to turbulent dissipation. We find that for observed molecular clouds, which typically have Alfven Mach numbers of ~1 (Crutcher 1999) and AD Reynolds numbers of ~20 (McKee et al. 2010), about 70% of the total turbulent dissipation is in the form of AD heating. AD has an important effect on the length scale where energy is dissipated: when AD heating is strong, most of the energy in the cascade is removed by ion-neutral drift, with a comparatively small amount of energy making it down to small scales. We derive a relation for the AD heating rate that describes the results of our simulations to within a factor of two. Turbulent dissipation, including AD heating, is generally less important that cosmic-ray heating in molecular clouds, although there is substantial scatter in both.

We use the photometric redshift method of Chakrabarti & McKee (2008) to infer photometric redshifts of submillimeter galaxies with far-IR (FIR) $\it{Herschel}$ data obtained as part of the PACS Evolutionary Probe (PEP) program. For the sample with spectroscopic redshifts, we demonstrate the validity of this method over a large range of redshifts ($ 4 \ga z \ga 0.3$) and luminosities, finding an average accuracy in $(1+z_{\rm phot})/(1+z_{\rm spec})$ of 10%. Thus, this method is more accurate than other FIR photometric redshift methods. This method is different from typical FIR photometric methods in deriving redshifts from the light-to-gas mass ($L/M$) ratio of infrared-bright galaxies inferred from the FIR spectral energy distribution (SED), rather than dust temperatures. Once the redshift is derived, we can determine physical properties of infrared bright galaxies, including the temperature variation within the dust envelope, luminosity, mass, and surface density. We use data from the GOODS-S field to calculate the star formation rate density (SFRD) of sub-mm bright sources detected by AzTEC and PACS. The AzTEC-PACS sources, which have a threshold $850 \micron$ flux $\ga 5 \rm mJy$, contribute 15% of the SFRD from all ULIRGs ($L_{\rm IR} \ga 10^{12} L_{\odot}$), and 3% of the total SFRD at $z \sim 2$.

[abridged] We report a series of simulations of the formation of a star cluster similar to the Orion Nebula Cluster (ONC), including both radiative transfer and protostellar outflows, and starting from both smooth and self-consistently turbulent initial conditions. Each simulation forms >150 stars and brown dwarfs, yielding a stellar mass distribution from < 0.1 to > 10 Msun. We show that a simulation that begins with self-consistently turbulence embedded in a larger turbulent volume, and that includes protostellar outflows, produces an initial mass function (IMF) consistent both with that of the ONC and the Galactic field. This is the first simulation published to date that reproduces the observed IMF in a cluster large enough to contain massive stars, and where the result is determined by a fully self-consistent calculation of gas thermodynamics. This simulation also produces a star formation rate that, while still somewhat too high, is much closer to observed values than if we omit either the larger turbulent volume or the outflows. Moreover, we show that the combination of outflows, self-consistently turbulent initial conditions, and turbulence continually fed by motions on scales larger than that of the protocluster yields an IMF that is in agreement with observations and invariant with time, resolving the "overheating" problem in which simulations without these features have an IMF peak that shifts to progressively higher masses over time. The simulation that matches the observed IMF also reproduces the observed trend of stellar multiplicity strongly increasing with mass. We show that this simulation produces massive stars from distinct massive cores whose properties are consistent with those of observed massive cores. However, the stars formed in these cores also undergo dynamical interactions that naturally produce Trapezium-like hierarchical multiple systems.

The properties of unresolved protostars and their local environment are frequently inferred from spectral energy distributions (SEDs) using radiative transfer modeling. We perform synthetic observations of realistic star formation simulations to evaluate the accuracy of properties inferred from fitting model SEDs to observations. We use ORION, an adaptive mesh refinement (AMR) three-dimensional gravito-radiation-hydrodynamics code, to simulate low-mass star formation in a turbulent molecular cloud including the effects of protostellar outflows. To obtain the dust temperature distribution and SEDs of the forming protostars, we post-process the simulations using HYPERION, a state-of-the-art Monte-Carlo radiative transfer code. We find that the ORION and HYPERION dust temperatures typically agree within a factor of two. We compare synthetic SEDs of embedded protostars for a range of evolutionary times, simulation resolutions, aperture sizes, and viewing angles. We demonstrate that complex, asymmetric gas morphology leads to a variety of classifications for individual objects as a function of viewing angle. We derive best-fit source parameters for each SED through comparison with a pre-computed grid of radiative transfer models. While the SED models correctly identify the evolutionary stage of the synthetic sources as embedded protostars, we show that the disk and stellar parameters can be very discrepant from the simulated values. Parameters such as the stellar accretion rate, stellar mass, and disk mass show better agreement, but can still deviate significantly, and the agreement may in some cases be artificially good due to the limited range of parameters in the set of model SEDs. Lack of correlation between the model and simulation properties in many individual instances cautions against over-interpreting properties inferred from SEDs for unresolved protostellar sources. (Abridged)

Protostellar feedback, both radiation and bipolar outflows, dramatically affects the fragmentation and mass accretion from star-forming cores. We use ORION, an adaptive mesh refinement gravito-radiation-hydrodynamics code, to simulate the formation of a cluster of low-mass stars, including both radiative transfer and protostellar outflows. We ran four simulations to isolate the individual effects of radiation feedback and outflow feedback as well as the combination of the two. Outflows reduce protostellar masses and accretion rates each by a factor of three and therefore reduce protostellar luminosities by an order of magnitude. Thus, while radiation feedback suppresses fragmentation, outflows render protostellar radiation largely irrelevant for low-mass star formation above a mass scale of 0.05 M_sun. We find initial fragmentation of our cloud at half the global Jeans length, ~ 0.1 pc. With insufficient protostellar radiation to stop it, these 0.1 pc cores fragment repeatedly, forming typically 10 stars each. The accretion rate in these stars scales with mass as predicted from core accretion models that include both thermal and turbulent motions. We find that protostellar outflows do not significantly affect the overall cloud dynamics, in the absence of magnetic fields, due to their small opening angles and poor coupling to the dense gas. The outflows reduce the mass from the cores by 2/3, giving a core to star efficiency ~ 1/3. The simulation with radiation and outflows reproduces the observed protostellar luminosity function. All of the simulations can reproduce observed core mass functions, though they are sensitive to telescope resolution. The simulation with both radiation and outflows reproduces the galactic IMF and the two-point correlation function of the cores observed in rho Oph.

We have carried out a numerical study of the effect of large scale magnetic fields on the rate of accretion from a uniform, isothermal gas onto a resistive, stationary point mass. Only mass, not magnetic flux, accretes onto the point mass. The simulations for this study avoid complications arising from boundary conditions by keeping the boundaries far from the accreting object. Our simulations leverage adaptive refinement methodology to attain high spatial fidelity close to the accreting object. Our results are particularly relevant to the problem of star formation from a magnetized molecular cloud in which thermal energy is radiated away on time scales much shorter than the dynamical time scale. Contrary to the adiabatic case, our simulations show convergence toward a finite accretion rate in the limit in which the radius of the accreting object vanishes, regardless of magnetic field strength. For very weak magnetic fields, the accretion rate first approaches the Bondi value and then drops by a factor ~ 2 as magnetic flux builds up near the point mass. For strong magnetic fields, the steady-state accretion rate is reduced by a factor ~ 0.2 \beta^{1/2} compared to the Bondi value, where \beta is the ratio of the gas pressure to the magnetic pressure. We give a simple expression for the accretion rate as a function of the magnetic field strength. Approximate analytic results are given in the Appendixes for both time-dependent accretion in the limit of weak magnetic fields and steady-state accretion for the case of strong magnetic fields.

[abridged] While observations of Local Group galaxies show a very simple, local star formation law in which the star formation rate per unit area in each patch of a galaxy scales linearly with the molecular gas surface density, recent observations of both Milky Way molecular clouds and high redshift galaxies apparently show a more complicated relationship, in which regions of equal surface density can form stars at quite different rates. These data have been interpreted as implying either that different star formation laws apply in different circumstances, that the star formation law is sensitive to large-scale galaxy properties rather than local properties, or that there are high density thresholds for star formation. Here we collate resolved observations of Milky Way molecular clouds, kpc-scale observations of Local Group galaxies, and unresolved observations of both disk and starburst galaxies in the local universe and at high redshift. We show that all of these data are in fact consistent with a simple, local, volumetric star formation law. The apparent variations stem from the fact that the observed objects have a wide variety of 3D size scales and degrees of internal clumping, so even at fixed gas column density the regions being observed can have wildly varying volume densities. We provide a simple theoretical framework to remove this projection effect, and we show that all the data, from small Solar neighborhood clouds with masses ~10^3 Msun to sub-mm galaxies with masses ~10^11 Msun, fall on a single star formation law in which the SFR is simply ~1% of the molecular gas mass per local free-fall time. In contrast, proposed star formation laws in which the star formation timescale is set by the galactic rotation period or the SFR is linearly proportional to the gas mass above some density threshold fail to match at least some of the data.

Performing a stable, long duration simulation of driven MHD turbulence with a high thermal Mach number and a strong initial magnetic field is a challenge to high-order Godunov ideal MHD schemes because of the difficulty in guaranteeing positivity of the density and pressure. We have implemented a robust combination of reconstruction schemes, Riemann solvers, limiters, and Constrained Transport EMF averaging schemes that can meet this challenge, and using this strategy, we have developed a new Adaptive Mesh Refinement (AMR) MHD module of the ORION2 code. We investigate the effects of AMR on several statistical properties of a turbulent ideal MHD system with a thermal Mach number of 10 and a plasma $\beta_0$ of 0.1 as initial conditions; our code is shown to be stable for simulations with higher Mach numbers ($M_rms = 17.3$) and smaller plasma beta ($\beta_0 = 0.0067$) as well. Our results show that the quality of the turbulence simulation is generally related to the volume-averaged refinement. Our AMR simulations show that the turbulent dissipation coefficient for supersonic MHD turbulence is about 0.5, in agreement with unigrid simulations.

Ambipolar diffusion (AD) is believed to be a crucial process for redistributing magnetic flux in the dense molecular gas that occurs in regions of star formation. We carry out numerical simulations of this process in regions of low ionization using the heavy ion approximation. The simulations are for regions of strong field (plasma \beta=0.1) and mildly supersonic turbulence (M=3, corresponding to an Alfven mach number of 0.67). The velocity power spectrum of the neutral gas changes from an Iroshnikov-Kraichnan spectrum in the case of ideal MHD to a Burgers spectrum in the case of a shock-dominated hydrodynamic system. The magnetic power spectrum shows a similar behavior. We use a 1D radiative transfer code to post-process our simulation results; the simulated emission from the CS J=2-1 and H13CO+ J=1-0 lines shows that the effects of AD are observable in principle. Linewidths of ions are observed to be less than those of neutrals, and we confirm previous suggestions that this is due to AD. We show that AD is unlikely to affect the Chandrasekhar-Fermi method for inferring field strengths unless the AD is stronger than generally observed. Finally, we present the first fully 3D study of the enhancement of AD by turbulence, finding that AD is accelerated by factor 2-4.5 for non self-gravitating systems with the level of turbulence we consider.

We report the results of a series of AMR radiation-hydrodynamic simulations of the collapse of massive star forming clouds using the ORION code. These simulations are the first to include the feedback effects protostellar outflows, as well as protostellar radiative heating and radiation pressure exerted on the infalling, dusty gas. We find that that outflows evacuate polar cavities of reduced optical depth through the ambient core. These enhance the radiative flux in the poleward direction so that it is 1.7 to 15 times larger than that in the midplane. As a result the radiative heating and outward radiation force exerted on the protostellar disk and infalling cloud gas in the equatorial direction are greatly diminished. The simultaneously reduces the Eddington radiation pressure barrier to high-mass star formation and increases the minimum threshold surface density for radiative heating to suppress fragmentation compared to models that do not include outflows. The strength of both these effects depends on the initial core surface density. Lower surface density cores have longer free-fall times and thus massive stars formed within them undergo more Kelvin contraction as the core collapses, leading to more powerful outflows. Furthermore, in lower surface density clouds the ratio of the time required for the outflow to break out of the core to the core free-fall time is smaller, so that these clouds are consequently influenced by outflows at earlier stages of collapse. As a result, outflow effects are strongest in low surface density cores and weakest in high surface density one. We also find that radiation focusing in the direction of outflow cavities is sufficient to prevent the formation of radiation pressure-supported circumstellar gas bubbles, in contrast to models which neglect protostellar outflow feedback.

One model for the origin of typical galactic star clusters such as the Orion Nebula Cluster (ONC) is that they form via the rapid, efficient collapse of a bound gas clump within a larger, gravitationally-unbound giant molecular cloud. However, simulations in support of this scenario have thus far have not included the radiation feedback produced by the stars; radiative simulations have been limited to significantly smaller or lower density regions. Here we use the ORION adaptive mesh refinement code to conduct the first ever radiation-hydrodynamic simulations of the global collapse scenario for the formation of an ONC-like cluster. We show that radiative feedback has a dramatic effect on the evolution: once the first ~10-20% of the gas mass is incorporated into stars, their radiative feedback raises the gas temperature high enough to suppress any further fragmentation. However, gas continues to accrete onto existing stars, and, as a result, the stellar mass distribution becomes increasingly top-heavy, eventually rendering it incompatible with the observed IMF. Systematic variation in the location of the IMF peak as star formation proceeds is incompatible with the observed invariance of the IMF between star clusters, unless some unknown mechanism synchronizes the IMFs in different clusters by ensuring that star formation is always truncated when the IMF peak reaches a particular value. We therefore conclude that the global collapse scenario, at least in its simplest form, is not compatible with the observed stellar IMF. We speculate that processes that slow down star formation, and thus reduce the accretion luminosity, may be able to resolve the problem.

We present virial models for the global evolution of giant molecular clouds. Focusing on the presence of an accretion flow, and accounting for the amount of mass, momentum, and energy supplied by accretion and star formation feedback, we are able to follow the growth, evolution, and dispersal of individual giant molecular clouds. Our model clouds reproduce the scaling relations observed in both galactic and extragalactic clouds. We find that accretion and star formation contribute contribute roughly equal amounts of turbulent kinetic energy over the lifetime of the cloud. Clouds attain virial equilibrium and grow in such a way as to maintain roughly constant surface densities, with typical surface densities of order 50 - 200 Msun pc^-2, in good agreement with observations of giant molecular clouds in the Milky Way and nearby external galaxies. We find that as clouds grow, their velocity dispersion and radius must also increase, implying that the linewidth-size relation constitutes an age sequence. Lastly, we compare our models to observations of giant molecular clouds and associated young star clusters in the LMC and find good agreement between our model clouds and the observed relationship between H ii regions, young star clusters, and giant molecular clouds.

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

The stellar initial mass function (IMF), along with the star formation rate, is one of the fundamental properties that any theory of star formation must explain. An interesting feature of the IMF is that it appears to be remarkably universal across a wide range of environments. Particularly, there appears to be little variation in either the characteristic mass of the IMF or its high-mass tail between clusters with different metallicities. Previous attempts to understand this apparent independence of metallicity have not accounted for radiation feedback from high-mass protostars, which can dominate the energy balance of the gas in star-forming regions. We extend this work, showing that the fragmentation of molecular gas should depend only weakly on the amount of dust present, even when the primary heating source is radiation from massive protostars. First, we report a series of core collapse simulations using the ORION AMR code that systematically vary the dust opacity and show explicitly that this has little effect on the temperature or fragmentation of the gas. Then, we provide an analytic argument for why the IMF varies so little in observed star clusters, even as the metallicity varies by a factor of 100.

We report new mid-infrared observations of the remarkable object IRAS 15099-5856 using the space telescopes AKARI and Spitzer, which demonstrate the presence of prominent crystalline silicate emission in this bright source. IRAS 15099-5856 has a complex morphology with a bright central compact source (IRS1) surrounded by knots, spurs, and several extended (~4') arc-like filaments. The source is seen only at >= 10 um. The Spitzer MIR spectrum of IRS1 shows prominent emission features from Mg-rich crystalline silicates, strong [Ne II] 12.81 um and several other faint ionic lines. We model the MIR spectrum as thermal emission from dust and compare with the Herbig Be star HD 100546 and the luminous blue variable R71, which show very similar MIR spectra. Molecular line observations reveal two molecular clouds around the source, but no associated dense molecular cores. We suggest that IRS1 is heated by UV radiation from the adjacent O star Muzzio 10 and that its crystalline silicates most likely originated in a mass outflow from the progenitor of the supernova remnant (SNR) MSH 15-52. IRS1, which is embedded in the SNR, could have been shielded from the SN blast wave if the progenitor was in a close binary system with Muzzio 10. If MSH15-52 is a remnant of Type Ib/c supernova (SN Ib/c), as has been previously proposed, this would confirm the binary model for SN Ib/c. IRS1 and the associated structures may be the relics of massive star death, as shaped by the supernova explosion, the pulsar wind and the intense ionizing radiation of the embedded O star.

Nearby spiral galaxies show an extremely tight correlation between tracers of molecular hydrogen (H_2) in the interstellar medium (ISM) and tracers of recent star formation, but it is unclear whether this correlation is fundamental or accidental. In the galaxies that have been surveyed to date, H_2 resides predominantly in gravitationally bound clouds cooled by carbon monoxide (CO) molecules, but in galaxies of low metal content the correlations between bound clouds, CO, and H_2 break down, and it is unclear if the star formation rate will then correlate with H_2 or with some other quantity. Here we show that star formation will continue to follow H_2 independent of metallicity. This is not because H_2 is directly important for cooling, but instead because the transition from predominantly atomic hydrogen (HI) to H_2 occurs under the same conditions as a dramatic drop in gas temperature and Bonnor-Ebert mass that destabilizes clouds and initiates collapse. We use this model to compute how star formation rate will correlate with total gas mass, with mass of gas where the hydrogen is H_2, and with mass of gas where the carbon is CO in galaxies of varying metallicity, and show that preliminary observations match the trend we predict.

Low-mass protostars are less luminous than expected. This luminosity problem is important because the observations appear to be inconsistent with some of the basic premises of star formation theory. Two possible solutions are that stars form slowly, which is supported by recent data, and/or that protostellar accretion is episodic; current data suggest that the latter accounts for less than half the missing luminosity. The solution to the luminosity problem bears directly on the fundamental problem of the time required to form a low-mass star. The protostellar mass and luminosity functions provide powerful tools both for addressing the luminosity problem and for testing theories of star formation. Results are presented for the collapse of singular isothermal spheres, for the collapse of turbulent cores, and for competitive accretion.

We derive a semi-empirical galactic initial mass function (IMF) from observational constraints. We assume that the star formation rate in a galaxy can be expressed as the product of the IMF, $\psi (m)$, which is a smooth function of mass $m$ (in units of \msun), and a time- and space-dependent total rate of star formation per unit area of galactic disk. The mass dependence of the proposed IMF is determined by five parameters: the low-mass slope $\gamma$, the high-mass slope $-\Gamma$, the characteristic mass $m_{ch}$ (which is close to the mass $m_{\rm peak}$ at which the IMF turns over), and the lower and upper limits on the mass, $m_l$ (taken to be 0.004) and $m_u$ (taken to be 120). The star formation rate in terms of number of stars per unit area of galactic disk per unit logarithmic mass interval, is proportional to $m^{-\Gamma} \left\{1-\exp\left[{-(m/m_{ch})^{\gamma +\Gamma}}\right]\right\}$, where $\cal N_*$ is the number of stars, $m_l<m<m_u$ is the range of stellar masses. The values of $\gamma$ and $\emch$ are derived from two integral constraints: i) the ratio of the number density of stars in the range $m=0.1-0.6$ to that in the range $m=0.6-0.8$ as inferred from the mass distribution of field stars in the local neighborhood, and ii) the ratio of the number of stars in the range $m=0.08 - 1$ to the number of brown dwarfs in the range $m=0.03-0.08$ in young clusters. The IMF satisfying the above constraints is characterized by the parameters $\gamma=0.51$ and $\emch=0.35$ (which corresponds to $m_{\rm peak}=0.27$). This IMF agrees quite well with the Chabrier (2005) IMF for the entire mass range over which we have compared with data, but predicts significantly more stars with masses $< 0.03\, M_\odot$; we also compare with other IMFs in current use.

We develop a model for regulation of galactic star formation rates Sigma_SFR in disk galaxies, in which ISM heating by stellar UV plays a key role. By requiring simultaneous thermal and (vertical) dynamical equilibrium in the diffuse gas, and star formation at a rate proportional to the mass of the self-gravitating component, we obtain a prediction for Sigma_SFR as a function of the total gaseous surface density Sigma and the density of stars + dark matter, rho_sd. The physical basis of this relationship is that thermal pressure in the diffuse ISM, which is proportional to the UV heating rate and therefore to Sigma_SFR, must adjust to match the midplane pressure set by the vertical gravitational field. Our model applies to regions where Sigma < 100 Msun/pc^2. In low-Sigma_SFR (outer-galaxy) regions where diffuse gas dominates, the theory predicts Sigma_SFR \propto Sigma (rho_sd)^1/2. The decrease of thermal equilibrium pressure when Sigma_SFR is low implies, consistent with observations, that star formation can extend (with declining efficiency) to large radii in galaxies, rather than having a sharp cutoff. The main parameters entering our model are the ratio of thermal pressure to total pressure in the diffuse ISM, the fraction of diffuse gas that is in the warm phase, and the star formation timescale in self-gravitating clouds; all of these are (in principle) direct observables. At low surface density, our model depends on the ratio of the mean midplane FUV intensity (or thermal pressure in the diffuse gas) to the star formation rate, which we set based on Solar neighborhood values. We compare our results to recent observations, showing good agreement overall for azimuthally-averaged data in a set of spiral galaxies. For the large flocculent spiral galaxies NGC 7331 and NGC 5055, the correspondence between theory and observation is remarkably close.

Ambipolar diffusion is important in redistributing magnetic flux and in damping Alfven waves in molecular clouds. The importance of ambipolar diffusion on a length scale $\ell$ is governed by the ambipolar diffusion Reynolds number, $\rad=\ell/\lad$, where $\lad$ is the characteristic length scale for ambipolar diffusion. The logarithmic mean of the AD Reynolds number in a sample of 15 molecular clumps with measured magnetic fields (Crutcher 1999) is 17, comparable to the theoretically expected value. We identify several regimes of ambipolar diffusion in a turbulent medium, depending on the ratio of the flow time to collision times between ions and neutrals; the clumps observed by Crutcher (1999) are all in the standard regime of ambipolar diffusion, in which the neutrals and ions are coupled over a flow time. We have carried out two-fluid simulations of ambipolar diffusion in isothermal, turbulent boxes for a range of values of $\rad$. The mean Mach numbers were fixed at $\calm=3$ and $\ma=0.67$; self-gravity was not included. We study the properties of overdensities--i.e., clumps--in the simulation and show that the slope of the higher-mass portion of the clump mass spectrum increases as $\rad$ decreases, which is qualitatively consistent with Padoan et al. (2007)'s finding that the mass spectrum in hydrodynamic turbulence is significantly steeper than in ideal MHD turbulence. For a value of $\rad$ similar to the observed value, we find a slope that is consistent with that of the high-mass end of the Initial Mass Function for stars. However, the value we find for the spectral index in our ideal MHD simulation differs from theirs, presumably because our simulations have different initial conditions. This suggests that the mass spectrum of the clumps in the Padoan et al. (2007) turbulent fragmentation model for the IMF depends on the environment, which would conflict with evidence ...

The protostellar mass function (PMF) is the Present-Day Mass Function of the protostars in a region of star formation. It is determined by the initial mass function weighted by the accretion time. The PMF thus depends on the accretion history of protostars and in principle provides a powerful tool for observationally distinguishing different protostellar accretion models. We consider three basic models here: the Isothermal Sphere model (Shu 1977), the Turbulent Core model (McKee & Tan 2003), and an approximate representation of the Competitive Accretion model (Bonnell et al. 1997, 2001a). We also consider modified versions of these accretion models, in which the accretion rate tapers off linearly in time. Finally, we allow for an overall acceleration in the rate of star formation. At present, it is not possible to directly determine the PMF since protostellar masses are not currently measurable. We carry out an approximate comparison of predicted PMFs with observation by using the theory to infer the conditions in the ambient medium in several star-forming regions. Tapered and accelerating models generally agree better with observed star-formation times than models without tapering or acceleration, but uncertainties in the accretion models and in the observations do not allow one to rule out any of the proposed models at present. The PMF is essential for the calculation of the Protostellar Luminosity Function, however, and this enables stronger conclusions to be drawn (Offner & McKee 2010).

We present computations of the ionization structure and metal-absorption properties of thermally conductive interface layers that surround evaporating warm spherical clouds, embedded in a hot medium. We rely on the analytical formalism of Dalton & Balbus to calculate the temperature profile in the evaporating gas, and explicitly solve the time-dependent ionization equations for H, He, C, N, O, Si, and S in the interface. We include photoionization by an external field. We estimate how departures from equilibrium ionization affect the resonance-line cooling efficiencies in the evaporating gas, and determine the conditions for which radiative losses may be neglected in the solution for the evaporation dynamics and temperature profile. Our results indicate that non-equilibrium cooling significantly increases the value of the saturation parameter at which radiative losses begin to affect the flow dynamics. As applications we calculate the ion fractions and projected column densities arising in the evaporating layers surrounding dwarf-galaxy-scale objects that are also photoionized by metagalactic radiation. We compare our results to the UV metal-absorption column densities observed in local highly-ionized metal-absorbers, located in the Galactic corona or intergalactic medium. Conductive interfaces significantly enhance the formation of high-ions such as C^3+, N^4+, and O^5+ relative to purely photoionized clouds, especially for clouds embedded in a high-pressure corona. However, the enhanced columns are still too low to account for the O VI columns (~1e14 cm^-2) observed in the local high-velocity absorbers. We find that O VI columns larger than ~1e13 cm^-2 cannot be produced in evaporating clouds. Our results do support the conclusion of Savage & Lehner, that absorption due to evaporating O VI likely occurs in the local interstellar medium, with characteristic columns of ~1e13 cm^-2.

The mass of molecular gas in an interstellar cloud is often measured using line emission from low rotational levels of CO, which are sensitive to the CO mass, and then scaling to the assumed molecular hydrogen H_2 mass. However, a significant H_2 mass may lie outside the CO region, in the outer regions of the molecular cloud where the gas phase carbon resides in C or C+. Here, H_2 self-shields or is shielded by dust from UV photodissociation, where as CO is photodissociated. This H_2 gas is "dark" in molecular transitions because of the absence of CO and other trace molecules, and because H_2 emits so weakly at temperatures 10 K < T < 100 K typical of this molecular component. This component has been indirectly observed through other tracers of mass such as gamma rays produced in cosmic ray collisions with the gas and far-infrared/submillimeter wavelength dust continuum radiation. In this paper we theoretically model this dark mass and find that the fraction of the molecular mass in this dark component is remarkably constant (~ 0.3 for average visual extinction through the cloud with mean A_V ~ 8) and insensitive to the incident ultraviolet radiation field strength, the internal density distribution, and the mass of the molecular cloud as long as mean A_V, or equivalently, the product of the average hydrogen nucleus column and the metallicity through the cloud, is constant. We also find that the dark mass fraction increases with decreasing mean A_V, since relatively more molecular H_2 material lies outside the CO region in this case.

The fragmentation of star-forming interstellar clouds, and the resulting stellar initial mass function (IMF), is strongly affected by the temperature structure of the collapsing gas. Since radiation feedback from embedded stars can modify this as collapse proceeds, feedback plays an important role in determining the IMF. However, the effects and importance of radiative heating are likely to depend strongly on the surface density of the collapsing clouds, which determines both their effectiveness at trapping radiation and the accretion luminosities of the stars forming within them. In this paper we report a suite of adaptive mesh refinement radiation-hydrodynamic simulations using the ORION code in which we isolate the effect of column density on fragmentation by following the collapse of clouds of varying column density while holding the mass, initial density and velocity structure, and initial virial ratio fixed. We find that radiation does not significantly modify the overall star formation rate or efficiency, but that it suppresses fragmentation more and more as cloud surface densities increase from those typical of low mass star-forming regions like Taurus, through the typical surface density of massive star-forming clouds in the Galaxy, up to conditions found only in super star clusters. In regions of low surface density, fragmentation during collapse leads to the formation of small clusters rather than individual massive star systems, greatly reducing the fraction of the stellar population with masses >~ 10 Msun. Our simulations have important implications for the formation of massive stars and the universality of the IMF.

Understanding the molecular content of galaxies is a critical problem in star formation and galactic evolution. Here we present a new method, based on a Stromgren-type analysis, to calculate the amount of HI that surrounds a molecular cloud irradiated by an isotropic radiation field. We consider both planar and spherical clouds, and H_2 formation either in the gas phase or catalyzed by dust grains. Under the assumption that the transition from atomic to molecular gas is sharp, our method gives the solution without any reference to the photodissociation cross section. We test our results for the planar case against those of a PDR code, and find typical accuracies of about 10%. Our results are also consistent with the scaling relations found in Paper I of this series, but they apply to a wider range of physical conditions. We present simple, accurate analytic fits to our results that are suitable for comparison to observations and to implementation in numerical and semi-analytic models.

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

Observations made using large ground-based and space-borne telescopes have probed cosmic history all the way from the present-day to a time when the Universe was less than a tenth of its present age. Earlier on lies the remaining frontier, where the first stars, galaxies, and massive black holes formed. They fundamentally transformed the early Universe by endowing it with the first sources of light and chemical elements beyond the primordial hydrogen and helium produced in the Big Bang. The interplay of theory and upcoming observations promises to answer the key open questions in this emerging field.

We propose a simple theoretical model for star formation in which the local star formation rate in a galaxy is determined by three factors. First, the interplay between the interstellar radiation field and molecular self-shielding determines what fraction of the gas is in molecular form and thus eligible to form stars. Second, internal feedback determines the properties of the molecular clouds that form, which are nearly independent of galaxy properties until the galactic ISM pressure becomes comparable to the internal GMC pressure. Above this limit, galactic ISM pressure determines molecular gas properties. Third, the turbulence driven by feedback processes in GMCs makes star formation slow, allowing a small fraction of the gas to be converted to stars per free-fall time within the molecular clouds. We combine analytic estimates for each of these steps to formulate a single star formation law, and show that the predicted correlation between star formation rate, metallicity, and surface densities of atomic, molecular, and total gas agree well with observations.

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

Gas in galactic disks is collected by gravitational instabilities into giant atomic-molecular complexes, but only the inner, molecular parts of these structures are able to collapse to form stars. Determining what controls the ratio of atomic to molecular hydrogen in complexes is therefore a significant problem in star formation and galactic evolution. In this paper we use the model of H_2 formation, dissociation, and shielding developed in the previous paper in this series to make theoretical predictions for atomic to molecular ratios as a function of galactic properties. We find that the molecular fraction in a galaxy is determined primarily by its column density and secondarily by its metallicity, and is to good approximation independent of the strength of the interstellar radiation field. We show that the column of atomic hydrogen required to shield a molecular region against dissociation is ~10 Msun pc^-2 at solar metallicity. We compare our model to data from recent surveys of the Milky Way and of nearby galaxies, and show that the both the primary dependence of molecular fraction on column density and the secondary dependence on metallicity that we predict are in good agreement with observed galaxy properties.

In this series of papers we study the structure of the atomic to molecular transition in the giant atomic-molecular complexes that are the repositories of most molecular gas in galaxies, with the ultimate goal of attaining a better understanding of what determines galaxies' molecular content. Here we derive an approximate analytic solution for the structure of a photodissociation region (PDR) in a cloud of finite size that is bathed in an external dissociating radiation field. Our solution extends previous work, which with few exceptions has been restricted to a one-dimensional treatment of the radiation field. We show that our analytic results compare favorably to exact numerical calculations in the one-dimensional limit. However, our more general geometry provides a more realistic representation than a semi-infinite slab for atomic-molecular complexes exposed to the interstellar radiation field, particularly in environments such as low-metallicity dwarf galaxies where the curvature and finite size of the atomic envelope cannot be neglected. For clouds that are at least 20% molecular we obtain analytic expressions for the molecular fraction in terms of properties of the gas and radiation field that are accurate to tens of percent, while for clouds of lower molecular content we obtain upper limits. As a side benefit, our analysis helps clarify when self-shielding is the dominant process in H_2 formation, and under what circumstances shielding by dust makes a significant contribution.

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)

In this study we investigate the formation and properties of prestellar and protostellar cores using hydrodynamic, self-gravitating Adaptive Mesh Refinement simulations, comparing the cases where turbulence is continually driven and where it is allowed to decay. We model observations of these cores in the C$^{18}$O$(2\to 1)$, NH$_3(1,1)$, and N$_2$H$^+(1\to 0)$ lines, and from the simulated observations we measure the linewidths of individual cores, the linewidths of the surrounding gas, and the motions of the cores relative to one another. Some of these distributions are significantly different in the driven and decaying runs, making them potential diagnostics for determining whether the turbulence in observed star-forming clouds is driven or decaying. Comparing our simulations with observed cores in the Perseus and $\rho$ Ophiuchus clouds shows reasonably good agreement between the observed and simulated core-to-core velocity dispersions for both the driven and decaying cases. However, we find that the linewidths through protostellar cores in both simulations are too large compared to the observations. The disagreement is noticably worse for the decaying simulation, in which cores show highly supersonic infall signatures in their centers that decrease toward their edges, a pattern not seen in the observed regions. This result gives some support to the use of driven turbulence for modeling regions of star formation, but reaching a firm conclusion on the relative merits of driven or decaying turbulence will require more complete data on a larger sample of clouds as well as simulations that include magnetic fields, outflows, and thermal feedback from the protostars.

Molecular clouds are observed to be turbulent, but the origin of this turbulence is not well understood. As a result, there are two different approaches to simulating molecular clouds, one in which the turbulence is allowed to decay after it is initialized, and one in which it is driven. We use the adaptive mesh refinement (AMR) code, Orion, to perform high-resolution simulations of molecular cloud cores and protostars in environments with both driven and decaying turbulence. We include self-gravity, use a barotropic equation of state, and represent regions exceeding the maximum grid resolution with sink particles. We analyze the properties of bound cores such as size, shape, linewidth, and rotational energy, and we find reasonable agreement with observation. At high resolution, the different rates of core accretion in the two cases have a significant effect on protostellar system development. Clumps forming in a decaying turbulence environment produce high-multiplicity protostellar systems with Toomre-Q unstable disks that exhibit characteristics of the competitive accretion model for star formation. In contrast, cores forming in the context of continuously driven turbulence and virial equilibrium form smaller protostellar systems with fewer low-mass members. Our simulations of driven and decaying turbulence show some statistically significant differences, particularly in the production of brown dwarfs and core rotation, but the uncertainties are large enough that we are not able to conclude whether observations favor one or the other.

Most numerical investigations on the role of magnetic fields in turbulent molecular clouds (MCs) are based on ideal magneto-hydrodynamics (MHD). However, MCs are weakly ionized, so that the time scale required for the magnetic field to diffuse through the neutral component of the plasma by ambipolar diffusion (AD) can be comparable to the dynamical time scale. We have performed a series of 256^3 and 512^3 simulations on supersonic but sub-Alfvenic turbulent systems with AD using the Heavy-Ion Approximation developed in Li, McKee, & Klein (2006). Our calculations are based on the assumption that the number of ions is conserved, but we show that these results approximately apply to the case of time-dependent ionization in molecular clouds as well. Convergence studies allow us to determine the optimal value of the ionization mass fraction when using the heavy-ion approximation for low Mach number, sub-Alfvenic turbulent systems. We find that ambipolar diffusion steepens the velocity and magnetic power spectra compared to the ideal MHD case. Changes in the density PDF, total magnetic energy, and ionization fraction are determined as a function of the AD Reynolds number. The power spectra for the neutral gas properties of a strongly magnetized medium with a low AD Reynolds number are similar to those for a weakly magnetized medium; in particular, the power spectrum of the neutral velocity is close to that for Burgers turbulence.

The understanding of the formation, life, and death of Population III stars, as well as the impact that these objects had on later generations of structure formation, is one of the foremost issues in modern cosmological research and has been an active area of research during the past several years. We summarize the results presented at "First Stars III," a conference sponsored by Los Alamos National Laboratory, the Kavli Institute for Particle Astrophysics and Cosmology, and the Joint Institute for Nuclear Astrophysics. This conference, the third in a series, took place in July 2007 at the La Fonda Hotel in Santa Fe, New Mexico, U.S.A.

Massive stars are very rare, but their extreme luminosities make them both the only type of young star we can observe in distant galaxies and the dominant energy sources in the universe today. They form rarely because efficient radiative cooling keeps most star-forming gas clouds close to isothermal as they collapse, and this favors fragmentation into stars <~1 Msun. Heating of a cloud by accreting low-mass stars within it can prevent fragmentation and allow formation of massive stars, but what properties a cloud must have to form massive stars, and thus where massive stars form in a galaxy, has not yet been determined. Here we show that only clouds with column densities >~ 1 g cm^-2 can avoid fragmentation and form massive stars. This threshold, and the environmental variation of the stellar initial mass function (IMF) that it implies, naturally explain the characteristic column densities of massive star clusters and the difference between the radial profiles of Halpha and UV emission in galactic disks. The existence of a threshold also implies that there should be detectable variations in the IMF with environment within the Galaxy and in the characteristic column densities of massive star clusters between galaxies, and that star formation rates in some galactic environments may have been systematically underestimated.

We report the discovery of a star-forming loop around the young, Crab-like supernova remnant (SNR) G54.1+0.3 using the AKARI infrared satellite. The loop consists of at least eleven young stellar objects (YSOs) embedded in a ring-like diffuse emission of radius ~1'. The YSOs are bright in the mid-infrared and are also visible in the Spitzer Space Telescope Galactic plane survey images. Their Spitzer colors are similar to those of class II YSOs in [3.6]-[5.8] but significantly redder in [8]-[24], i.e., 0<[3.6]-[5.8]<1.2 and 5<[8]-[24]<9. Most of them have near-infrared counterparts in the 2MASS JHKs images, and some of them have an optical counterpart too. Their JHKs colors and magnitudes indicate that the YSOs are massive (<= 10 Msun) pre-main-sequence stars at the same distance to the SNR, i.e., 8 kpc, which supports the association of the star-forming loop with the SNR. The dereddened spectral energy distributions are similar to eraly Herbig Be stars, which are early B-type pre-main-sequence stars with inner disks that have been destroyed. The confinement to a loop structure indicates that the YSOs are young, i.e., <= 2 Myr. We propose that their formation is triggered by the progenitor star of G54.1+0.3, which has a mass of <= 15 Msun. The triggering must have occurred near the end of the progenitor's life, possibly after it had evolved off the main sequence.

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.

We infer the large-scale source parameters of dusty galaxies from their observed spectral energy distributions (SEDs) using the analytic radiative transfer methodology presented in Chakrabarti & McKee (2005). For local ultra-luminous infrared galaxies (ULIRGs), we show that the millimeter to far-infrared (FIR) SEDs can be well fit using the standard dust opacity index of 2 when self-consistent radiative transfer solutions are employed, indicating that the cold dust in local ULIRGs can be described by a single grain model. We develop a method for determining photometric redshifts of ULIRGs and sub-mm galaxies from the millimeter-FIR SED; the resulting value of $1+z$ is typically accurate to about 10%. As such, it is comparable to the accuracy of near-IR photometric redshifts and provides a complementary means of deriving redshifts from far-IR data, such as that from the upcoming $\it{Herschel Space Observatory}$. Since our analytic radiative transfer solution is developed for homogeneous, spherically symmetric, centrally heated, dusty sources, it is relevant for infrared bright galaxies that are primarily powered by compact sources of luminosity that are embedded in a dusty envelope. We discuss how deviations from spherical symmetry may affect the applicability of our solution, and we contrast our self-consistent analytic solution with standard approximations to demonstrate the main differences.

We review current understanding of star formation, outlining an overall theoretical framework and the observations that motivate it. A conception of star formation has emerged in which turbulence plays a dual role, both creating overdensities to initiate gravitational contraction or collapse, and countering the effects of gravity in these overdense regions. The key dynamical processes involved in star formation -- turbulence, magnetic fields, and self-gravity -- are highly nonlinear and multidimensional. Physical arguments are used to identify and explain the features and scalings involved in star formation, and results from numerical simulations are used to quantify these effects. We divide star formation into large-scale and small-scale regimes and review each in turn. Large scales range from galaxies to giant molecular clouds (GMCs) and their substructures. Important problems include how GMCs form and evolve, what determines the star formation rate (SFR), and what determines the initial mass function (IMF). Small scales range from dense cores to the protostellar systems they beget. We discuss formation of both low- and high-mass stars, including ongoing accretion. The development of winds and outflows is increasingly well understood, as are the mechanisms governing angular momentum transport in disks. Although outstanding questions remain, the framework is now in place to build a comprehensive theory of star formation that will be tested by the next generation of telescopes.

We analyze the mixed frame equations of radiation hydrodynamics under the approximations of flux-limited diffusion and a thermal radiation field, and derive the minimal set of evolution equations that includes all terms that are of leading order in any regime of non-relativistic radiation hydrodynamics. Our equations are accurate to first order in v/c in the static diffusion regime. In contrast, we show that previous lower order derivations of these equations omit leading terms in at least some regimes. In comparison to comoving frame formulations of radiation hydrodynamics, our equations have the advantage that they manifestly conserve total energy, making them very well-suited to numerical simulations, particularly with adaptive meshes. For systems in the static diffusion regime, our analysis also suggests an algorithm that is both simpler and faster than earlier comoving frame methods. We implement this algorithm in the Orion adaptive mesh refinement code, and show that it performs well in a range of test problems.

We compute the molecular line emission of massive protostellar disks by solving the equation of radiative transfer through the cores and disks produced by the recent radiation-hydrodynamic simulations of Krumholz, Klein, & McKee. We find that in several representative lines the disks show brightness temperatures of hundreds of Kelvin over velocity channels ~10 km s^-1 wide, extending over regions hundreds of AU in size. We process the computed intensities to model the performance of next-generation radio and submillimeter telescopes. Our calculations show that observations using facilities such as the EVLA and ALMA should be able to detect massive protostellar disks and measure their rotation curves, at least in the nearest massive star-forming regions. They should also detect significant sub-structure and non-axisymmetry in the disks, and in some cases may be able to detect star-disk velocity offsets of a few km s^-1, both of which are the result of strong gravitational instability in massive disks. We use our simulations to explore the strengths and weaknesses of different observational techniques, and we also discuss how observations of massive protostellar disks may be used to distinguish between alternative models of massive star formation.

We simulate the early stages of the evolution of turbulent, virialized, high-mass protostellar cores, with primary attention to how cores fragment, and whether they form a small or large number of protostars. Our simulations use the Orion adaptive mesh refinement code to follow the collapse from ~0.1 pc scales to ~10 AU scales, for durations that cover the main fragmentation phase, using three-dimensional gravito-radiation hydrodynamics. We find that for a wide range of initial conditions radiation feedback from accreting protostars inhibits the formation of fragments, so that the vast majority of the collapsed mass accretes onto one or a few objects. Most of the fragmentation that does occur takes place in massive, self-shielding disks. These are driven to gravitational instability by rapid accretion, producing rapid mass and angular momentum transport that allows most of the gas to accrete onto the central star rather than forming fragments. In contrast, a control run using the same initial conditions but an isothermal equation of state produces much more fragmentation, both in and out of the disk. We conclude that massive cores with observed properties are not likely to fragment into many stars, so that, at least at high masses, the core mass function probably determines the stellar initial mass function. Our results also demonstrate that simulations of massive star forming regions that do not include radiative transfer, and instead rely on a barotropic equation of state or optically thin heating and cooling curves, are likely to produce misleading results.

We present Arecibo L-band Feed Array 21-cm observations of a sub-complex of HVCs at the tip of the Anti-Center Complex. These observations show morphological details that point to interaction with the ambient halo medium and differential drag within the cloud sub-complex. We develop a new technique for measuring cloud distances, which relies upon these observed morphological and kinematic characteristics, and show that it is consistent with H-alpha distances. These results are consistent with distances to HVCs and halo densities derived from models in which HVCs are formed from cooling halo gas.

We present semi-analytic dynamical models for giant molecular clouds evolving under the influence of HII regions launched by newborn star clusters. In contrast to previous work, we neither assume that clouds are in virial or energetic equilibrium, nor do we ignore the effects of star formation feedback. The clouds, which we treat as spherical, can expand and contract homologously. Photoionization drives mass ejection; the recoil of cloud material both stirs turbulent motions and leads to an effective confining pressure. The balance between these effects and the decay of turbulent motions through isothermal shocks determines clouds' dynamical and energetic evolution. We find that for realistic values of the rates of turbulent dissipation, photoevaporation, and energy injection by HII regions, the massive clouds where most molecular gas in the Galaxy resides live for a few crossing times, in good agreement with recent observational estimates that large clouds in local group galaxies survive roughly 20-30 Myr. During this time clouds remain close to equilibrium, with virial parameters of 1-3 and column densities near 10^22 H atoms cm^-2, also in agreement with observed cloud properties. Over their lives they convert 5-10% of their mass into stars, after which point most clouds are destroyed when a large HII region unbinds them. In contrast, small clouds like those found in the solar neighborhood only survive ~1 crossing time before being destroyed.

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.

Recent Chandra and XMM observations of distant quasars have shown strong local (z~0) X-ray absorption lines from highly ionized gas, primarily He-like oxygen. The nature of these X-ray absorbers, i.e., whether they are part of the hot gas associated with the Milky Way or part of the intragroup medium in the Local Group, remains a puzzle due to the uncertainties in the distance. We present in this paper a survey of 20 AGNs with Chandra and XMM archival data. About 40% of the targets show local OVII He-alpha absorption with column densities around 10^16 cm^-2; in particular, OVII absorption is present in all the high quality spectra. We estimate that the sky covering fraction of this OVII-absorbing gas is at least 63%, at 90% confidence, and likely to be unity given enough high-quality spectra. Based on (1) the expected number of absorbers along sight lines toward distant AGNs, (2) joint analysis with X-ray emission measurements, and (3) mass estimation, we argue that the observed X-ray absorbers are part of the hot gas associated with our Galaxy. Future observations will significantly improve our understanding of the covering fraction and provide robust test of this result.

The effect of smooth cloud boundaries on the interaction of steady planar shock waves with interstellar clouds is studied using a high-resolution local AMR technique with a second-order accurate axisymmetric Godunov hydrodynamic scheme. A 3D calculation is also done to confirm the results of the 2D ones. We consider an initially spherical cloud whose density distribution is flat near the cloud center and has a power-law profile in the cloud envelope. When an incident shock is transmitted into a smooth cloud, velocity gradients in the cloud envelope steepen the smooth density profile at the upstream side, resulting in a sharp density jump having an arc-like shape. Such a ``slip surface'' forms immediately when a shock strikes a cloud with a sharp boundary. For smoother boundaries, the formation of slip surface and therefore the onset of hydrodynamic instabilities are delayed. Since the slip surface is subject to the Kelvin-Helmholtz and Rayleigh-Taylor instabilities, the shocked cloud is eventually destroyed in $\sim 3-10$ cloud crushing times. After complete cloud destruction, small blobs formed by fragmentation due to hydrodynamic instabilities have significant velocity dispersions of the order of 0.1 $v_b$, where $v_b$ is the shock velocity in the ambient medium. This suggests that turbulent motions generated by shock-cloud interaction are directly associated with cloud destruction. The interaction of a shock with a cold HI cloud should lead to the production of a spray of small HI shreds, which could be related to the small cold clouds recently observed by Stanimirovic & Heiles (2005). The linewidth-size relation obtained from our 3D simulation is found to be time-dependent. A possibility for gravitational instability triggered by shock compression is also discussed.

Stars with masses of >~ 20 solar masses have short Kelvin times that enable them to reach the main sequence while still accreting from their natal clouds. The resulting nuclear burning produces a huge luminosity and a correspondingly large radiation pressure force on dust grains in the accreting gas. This effect may limit the upper mass of stars that can form by accretion. Indeed, simulations and analytic calculations to date have been unable to resolve the mystery of how stars of 50 solar masses and up form. We present two new ideas to solve the radiation pressure problem. First, we use three-dimensional radiation hydrodynamic adaptive mesh refinement simulations to study the collapse of massive cores. We find that in three dimensions a configuration in which radiation holds up an infalling envelope is Rayleigh-Taylor unstable, leading radiation driven bubbles to collapse and accretion to continue. We also present Monte Carlo radiative transfer calculations showing that the cavities created by protostellar winds provides a valve that allow radiation to escape the accreting envelope, further reducing the ability of radiation pressure to inhibit accretion.

The Bondi-Hoyle formula gives the approximate accretion rate onto a point particle accreting from a uniform medium. However, in many situations accretion onto point particles occurs from media that are turbulent rather than uniform. In this paper, we give an approximate solution to the problem of a point particle accreting from an ambient medium of supersonically turbulent gas. Accretion in such media is bimodal, at some points resembling classical Bondi-Hoyle flow, and in other cases being closer to the vorticity-dominated accretion flows recently studied by Krumholz, McKee, & Klein. Based on this observation, we develop a theoretical prediction for the accretion rate, and confirm that our predictions are highly consistent with the results of numerical simulations. The distribution of accretion rates is lognormal, and the mean accretion rate in supersonically turbulent gas can be substantially enhanced above the value that would be predicted by a naive application of the Bondi-Hoyle formula. However, it can also be suppressed by the vorticity, just as Krumholz, McKee, & Klein found for non-supersonic vorticity-dominated flows. Magnetic fields, which we have not included in these models, may further inhibit accretion. Our results have significant implications for a number astrophysical problems, ranging from star formation to the black holes in galactic centers. In particular, there are likely to be significant errors in results that assume that accretion from turbulent media occurs at the unmodified Bondi-Hoyle rate, or that are based on simulations that do not resolve the Bondi-Hoyle radius of accreting objects.

There are now two dominant models of how stars form: gravitational collapse theory holds that star-forming molecular clumps, typically hundreds to thousands of solar masses in mass, fragment into gaseous cores that subsequently collapse to make individual stars or small multiple systems. In contrast, competitive accretion theory suggests that at birth all stars are much smaller than the typical stellar mass (~0.5 solar masses), and that final stellar masses are determined by the subsequent accretion of unbound gas from the clump. Competitive accretion models explain brown dwarfs and free-floating planets as protostars ejected from star-forming clumps before accreting much mass, predicting that they should lack disks, have high velocity dispersions, and form more frequently in denser clumps. They also predict that mean stellar mass should vary within the Galaxy. Here we derive a simple estimate for the rate of competitive accretion as a function of the star-forming environment, based partly on simulations, and determine in what types of environments competitive accretion can occur. We show that no observed star-forming region produces significant competitive accretion, and that simulations that show competitive accretion do so because their properties differ from those determined by observation. Our result shows that stars form by gravitational collapse, and explains why observations have failed to confirm predictions of the competitive accretion scenario.

We present evidence for cosmic ray acceleration at the forward shock in Tycho's supernova remnant (SNR) from three X-ray observables: (1) the proximity of the contact discontinuity to the forward shock, or blast wave, (2) the morphology of the emission from the rim of Tycho, and (3) the spectral nature of the rim emission. We determine the locations of the blast wave (BW), contact discontinuity (CD), and reverse shock (RS) around the rim of Tycho's supernova remnant using a principal component analysis and other methods applied to new Chandra data. The azimuthal-angle-averaged radius of the BW is 251". For the CD and RS we find average radii of 241" and 183", respectively. Taking account of projection effects, we find ratios of 1:0.93:0.70 (BW:CD:RS). We show these values to be inconsistent with adiabatic hydrodynamical models of SNR evolution. The CD:BW ratio can be explained if cosmic ray acceleration of ions is occurring at the forward shock. The RS:BW ratio, as well as the strong Fe Ka emission from the Tycho ejecta, imply that the RS is not accelerating cosmic rays. We also extract radial profiles from ~34% of the rim of Tycho and compare them to models of surface brightness profiles behind the BW for a purely thermal plasma with an adiabatic shock. The observed morphology of the rim is much more strongly peaked than predicted by the model, indicating that such thermal emission is implausible here. Spectral analysis also implies that the rim emission is non-thermal in nature, lending further support to the idea that Tycho's forward shock is accelerating cosmic rays.

We present analytic radiative transfer solutions for the spectra of unresolved, spherically symmetric, centrally heated, dusty sources. We find that the dust thermal spectrum possesses scaling relations that provide a natural classification for a broad range of sources, from low-mass protostars to dusty galaxies. In particular, we find that, given our assumptions, spectral energy distributions (SEDs) can be characterized by two distance-independent parameters, the luminosity-to-mass ratio, $L/M$, and the surface density, $\Sigma$, for a set of two functions, namely, the density profile and the opacity curve. The goal is to use SEDs as a diagnostic tool in inferring the large-scale physical conditions in protostellar and extragalactic sources, and ultimately, evolutionary parameters. Our approach obviates the need to use SED templates in the millimeter to far-infrared region of the spectrum; this is a common practice in the extragalactic community that relies on observed correlations established at low redshift that may not necessarily extend to high redshift. Further, we demarcate the limited region of parameter space in which density profiles can be inferred from the SED, which is of particular import in the protostellar community as competing theories of star formation are characterized by different density profiles. The functionality of our model is unique in that in provides for a self-consistent analytic solution that we have validated by comparison with a well-tested radiative transfer code (DUSTY) to find excellent agreement with numerical results over a parameter space that spans low-mass protostars to ultra-luminous infrared galaxies (ULIRGS).

We derive an analytic prediction for the star formation rate in environments ranging from normal galactic disks to starbursts and ULIRGs in terms of the observables of those systems. Our calculation is based on three premises: (1) star formation occurs in virialized molecular clouds that are supersonically turbulent; (2) the density distribution within these clouds is lognormal, as expected for supersonic isothermal turbulence; (3) stars form in any sub-region of a cloud that is so overdense that its gravitational potential energy exceeds the energy in turbulent motions. We show that a theory based on this model is consistent with simulations and with the observed star formation rate in the Milky Way. We use our theory to derive the Kennicutt-Schmidt Law from first principles, and make other predictions that can be tested by future observations. We also provide an algorithm for estimating the star formation rate that is suitable for inclusion in numerical simulations.

We consider the effects of an outflow on radiation escaping from the infalling envelope around a massive protostar. Using numerical radiative transfer calculations, we show that outflows with properties comparable to those observed around massive stars lead to significant anisotropy in the stellar radiation field, which greatly reduces the radiation pressure experienced by gas in the infalling envelope. This means that radiation pressure is a much less significant barrier to massive star formation than has previously been thought.

The classical Bondi-Hoyle formula gives the accretion rate onto a point particle of a gas with a uniform density and velocity. However, the Bondi-Hoyle problem considers only gas with no net vorticity, while in a real astrophysical situation accreting gas invariably has at least a small amount of vorticity. We therefore consider the related case of accretion of gas with constant vorticity, for the cases of both small and large vorticity. We confirm the findings of earlier two dimensional simulations that even a small amount of vorticity can substantially change both the accretion rate and the morphology of the gas flow lines. We show that in three dimensions the resulting flow field is non-axisymmetric and time dependent. The reduction in accretion rate is due to an accumulation of circulation near the accreting particle. Using a combination of simulations and analytic treatment, we provide an approximate formula for the accretion rate of gas onto a point particle as a function of the vorticity of the surrounding gas.

We present a second-order Godunov algorithm for multidimensional, ideal MHD. Our algorithm is based on the unsplit formulation of Colella (J. Comput. Phys. vol. 87, 1990), with all of the primary dependent variables centered at the same location. To properly represent the divergence-free condition of the magnetic fields, we apply a discrete projection to the intermediate values of the field at cell faces, and apply a filter to the primary dependent variables at the end of each time step. We test the method against a suite of linear and nonlinear tests to ascertain accuracy and stability of the scheme under a variety of conditions. The test suite includes rotated planar linear waves, MHD shock tube problems, low-beta flux tubes, and a magnetized rotor problem. For all of these cases, we observe that the algorithm is second-order accurate for smooth solutions, converges to the correct weak solution for problems involving shocks, and exhibits no evidence of instability or loss of accuracy due to the possible presence of non-solenoidal fields.

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 introduce a new computational method for embedding Lagrangian sink particles into an Eulerian calculation. Simulations of gravitational collapse or accretion generally produce regions whose density greatly exceeds the mean density in the simulation. These dense regions require extremely small time steps to maintain numerical stability. Smoothed particle hydrodynamics (SPH) codes approach this problem by introducing non-gaseous, accreting sink particles, and Eulerian codes may introduce fixed sink cells. However, until now there has been no approach that allows Eulerian codes to follow accretion onto multiple, moving objects. We have removed that limitation by extending the sink particle capability to Eulerian hydrodynamics codes. We have tested this new method and found that it produces excellent agreement with analytic solutions. In analyzing our sink particle method, we present a method for evaluating the disk viscosity parameter $\alpha$ due to the numerical viscosity of a hydrodynamics code, and use it to compute $\alpha$ for our Cartesian AMR code. We also present a simple application of this new method: studying the transition from Bondi to Bondi-Hoyle accretion that occurs when a shock hits a particle undergoing Bondi accretion.

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.

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.

Much of the interstellar medium in disk galaxies is in the form of neutral atomic hydrogen, H I. This gas can be in thermal equilibrium at relatively low temperatures, T < 300 K (the cold neutral medium, or CNM) or at temperatures somewhat less than 10^4 K (the warm neutral medium, or WNM). These two phases can coexist over a narrow range of pressures, P_min < P < P_max. We determine P_min and P_max in the plane of the Galaxy as a function of Galactocentric radius R using recent determinations of the gas heating rate and the gas phase abundances of interstellar gas. We provide an analytic approximation for P_min as a function of metallicity, far-ultraviolet radiation field, and the ionization rate of atomic hydrogen. Over most of the disk of the Galaxy, the H I must be in two phases: the weight of the H I in the gravitational potential of the Galaxy is large enough to generate thermal pressures exceeding P_min, so that turbulent pressure fluctuations can produce cold gas that is thermally stable; and the mean density of the H I is too low for the gas to be all CNM. Our models predict the presence of CNM gas to R = 16-18 kpc, somewhat farther than previous estimates. We also examine the potential impact of turbulent heating on our results and provide expressions for the heating rate as a function of Galactic radius.

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.

We calculate the coupled hydrostatic and ionization structures of pressure-supported gas clouds that are confined by gravitationally dominant dark-matter (DM) mini-halos and by an external bounding pressure provided by a hot medium. We focus on clouds that are photoionized and heated by the present-day background metagalactic field and determine the conditions for the formation of warm (WNM), and multi-phased (CNM/WNM) neutral atomic hydrogen (HI) cores in the DM-dominated clouds. We consider LCDM dark-matter halos, and we compute models for a wide range of halo masses, total cloud gas masses, and external bounding pressures. We present models for the pressure-supported HI structures observed in the Local Group dwarf galaxies Leo A and Sag DIG. We then construct minihalo models for the multi-phased (and low-metallicity) compact high-velocity HI clouds (CHVCs). If the CHVCs are drawn from the same family of halos that successfully reproduce the dwarf galaxy observations, then the CHVCs must be "circumgalactic objects" with characteristic distances of 150 kpc. For such systems we find that multi-phased behavior occurs for peak WNM HI column densities between 2e19 and 1e20 cm^-2, consistent with observations. If the large population of CHVCs represent "missing low-mass satellites" of the Galaxy, then these clouds must be pressure-confined to keep the gas neutral. For an implied CHVC minihalo scale velocity of v_s=12 km s^-1, the confining pressure must exceed ~50 cm^-3 K. A hot (2e6 K) Galactic corona could provide the required pressure at 150 kpc.

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.

Far Ultraviolet radiation has been recognized as the main source of heating of the neutral interstellar gas, and, as a consequence, it determines whether the thermal balance of the neutral gas results in cold ($T\sim 50 - 100 K$) clouds (CNM), warm ($T \sim 10^4 K$) clouds (WNM), or a combination of the two. High FUV fields convert the neutral gas to WNM, while low fields result in CNM. The sources of FUV radiation are the short-lived massive stars that generally originate in associations that form in Giant Molecular Clouds present in the galactic disk. Using McKee & Williams' (1997) distribution of birthrates for OB associations in the Galaxy, we determine the expected behavior of the time-dependent FUV field for random positions in the local ISM. The FUV field is calculated in two bands (912-1100 \AA and 912-2070 \AA) and at the wavelength 1400 \AA. Our median value for the [912-2070 \AA] band is $G_0=1.6$ times Habing's (1968) value for the radiation field at the solar circle in this band, and quite close to Draine's (1976) value, $G_0=1.7$. Due to attenuation by dust, only associations within about 500 pc contribute significantly to the energy density at a given point. The FUV field is asymmetric at a given point, and the asymmetry grows for higher fields. The FUV field fluctuates with a variety of amplitudes, the larger ones being less frequent. The mean field is about twice the median field because of these fluctuations, or spikes, in the radiation field. These spikes, which last $\sim 30$ Myr, are caused by the infrequent birth of nearby associations. We also model shorter duration spikes caused by runaway OB stars. The presence of a fluctuating heating rate created by the fluctuating FUV field converts CNM to WNM and vice versa.

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)

After reviewing recent work on the dynamics of young supernova remnants (YSNRs) and discussing how YSNRs can be used as physics laboratories, I discuss several challenges: Where are the very young SNRs in the Galaxy? Can very young SNRs produce gamma ray bursts? The Connections Challenge: Can one infer the nature of the supernova and its progenitor star from observations of the YSNR?

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.

Polytropic models are reasonably successful in acounting for the observed features of molecular clouds. Multi-pressure polytropes include the various pressure components that are important in molecular clouds, whereas composite polytropes provide a representation for the core halo structure. Small, very dense (n~10^{11} cm^{-3}) molecular clouds have been proposed as models for both dark matter and for extreme scattering events. Insofar as the equation of state in these clouds can be represented by a single polytropic relation (pressure varies as a power of the density), such models conflict with observation. It is possible to contrive composite polytropes that do not conflict with observation, but whether the thermal properties of the clouds are consistent with such structure remains to be determined.

Using a quantitative model for bipolar outflows driven by hydromagnetic protostellar winds, we calculate the efficiency of star formation assuming that available gas is either converted into stars or ejected in outflows. We estimate the efficiency of a single star formation event in a protostellar core, finding 25%-70% for cores with various possible degrees of flattening. The core mass function and the stellar initial mass function have similar slopes, because the efficiency is not sensitive to its parameters. We then consider the disruption of gas from a dense molecular clump in which a cluster of young stars is being born. In both cases, we present analytical formulae for the efficiencies that compare favorably against observations and, for clusters, against numerical simulations. We predict efficiencies in the range 30%-50% for the regions that form clusters of low-mass stars. In our model, star formation and gas dispersal happen concurrently. We neglect the destructive effects of massive stars: our results are therefore upper limits to the efficiency in regions more massive than about 3000 Msun.

We examine the conditions under which binary and multiple stars may form out of turbulent molecular cloud cores using high resolution 3-D, adaptive mesh refinement (AMR) hydrodynamics (Truelove et al., 1997, 1998; Klein, 1999). We argue that previous conclusions on the conditions for cloud fragmentation have limited applicability, since they did not allow for the nonlinear density and velocity perturbations that are ubiquitous in molecular cloud cores. Over the past year, we have begun to simulate the evolution of marginally stable, turbulent cores. These models have radii, masses, density contrasts, turbulent linewidths, and projected velocity gradients consistent with observations of low-mass molecular cloud cores. Our models are evolved in time under self-gravitational hydrodynamics with AMR using a barotropic equation of state that models the transition from an isothermal to an adiabatic equation of state. We examine several properties of the resulting protostellar fragments and discuss the qualitative nature of the fragmentation process in realistic cloud cores: the transition from single to binary and multiple stars; the formation of misaligned binary systems; and the role played by filament formation in the formation of stars.

We present new calculations on the contribution from cooling hot gas to the photoionization of warm ionized gas in the Galaxy. We show that hot gas in cooling supernova remnants (SNRs) is an important source of photoionization, particularly for gas in the halo. We find that in many regions at high latitude this source is adequate to account for the observed ionization so there is no need to find ways to transport stellar photons from the disk. The flux from cooling SNRs sets a floor on the ionization along any line of sight. Our model flux is also shown to be consistent with the diffuse soft X-ray background and with soft X-ray observations of external galaxies. We consider the ionization of the clouds observed towards the halo star HD 93521, for which there are no O stars close to the line of sight. We show that the observed ionization can be explained successfully by our model EUV/soft X-ray flux from cooling hot gas. In particular, we can match the H alpha intensity, the S++/S+ ratio, and the C+* column. From observations of the ratios of columns of C+* and either S+ or H0, we are able to estimate the thermal pressure in the clouds. The slow clouds require high (~10^4 cm^-3 K) thermal pressures to match the N(C+*)/N(S+) ratio. Additional heating sources are required for the slow clouds to maintain their ~7000 K temperatures at these pressures, as found by Reynolds, Hausen & Tufte (1999).

We demonstrate that magnetically-collimated protostellar winds will sweep ambient material into thin, radiative, momentum-conserving shells whose features reproduce those commonly observed in bipolar molecular outflows. We find the typical position-velocity and mass-velocity relations to occur in outflows in a wide variety of ambient density distributions, regardless of the time histories of their driving winds.

We construct spherical, hydrostatic models of dense molecular cores and Bok globules consisting of two distinct, spatially separate gas components: a central, isothermal region surrounded by a negative-index, polytropic envelope. The clouds are supported against their own self-gravity by a combination of thermal, mean magnetic, and turbulent wave pressure. The latter two are included by allowing for locally adiabatic, non-isentropic pressure components. Such models are meant to represent, in a schematic manner, the velocity and density structure of cores and globules, as inferred from molecular line and dust continuum observations. We show by explicit construction that it is possible to have dense cores comparable to the Jeans mass embedded in stable clouds of much larger mass. In a subsequent paper, we show that composite polytropes on the verge of gravitational instability can reproduce the observed velocity and density structure of cores and globules under a variety of physical conditions.

Molecular clouds are supported by thermal pressure, magnetic pressure, and turbulent pressure. Each of these can be modeled with a polytropic equation of state, so that overall the total pressure is the sum of the individual components. We model the turbulent pressure as being due to a superposition of Alfven waves. The theory of polytropes is generalized to allow for the flow of entropy in response to a perturbation, as expected for the entropy associated with wave pressure. The equation of state of molecular clouds is "soft", so that the properties of the clouds are generally governed by the conditions at the surface. In general, the polytropes are not isentropic, and this permits large density and pressure drops to occur between the center and the edge of the polytropes, as is observed.

We review the progress that has been made in observing and analyzing molecular cloud structure in recent years. Structures are self-similar over a wide range of scales with similar power law indices independent of the star forming nature of a cloud. Comparison of structures at parsec-scale resolution in a star forming and non-star forming cloud show that the average densities in the former are higher but the structural characteristics in each cloud are much the same. In gravitationally bound regions of a cloud, however, and at higher densities and resolution, the self-similar scaling relationships break down and it is possible to observe the first steps toward star formation. High resolution observations of the dense individual star forming cores within the clumps hold the key to an empirical understanding of the origins of the stellar initial mass function.

Giant molecular clouds (GMCs) are the sites of star formation in the Galaxy. Many of their properties can be understood in terms of a model in which the GMCs and the star-forming clumps within them are in approximate pressure equilibrium, with turbulent motions treated as a separate pressure component.

We examine the relation between presupernova stellar structure and the distribution of ejecta in core-collapse supernovae, assuming adiabatic, spherically symmetric flow. We develop a simple yet accurate formula for the blastwave shock velocity, and demonstrate that the entire final density distribution can be approximated with simple models for the final pressure distribution, along with the approximate shock-deposited entropy, in a way that matches the results of simulations. We find that the distribution of density in a star's ejecta depends on whether its outer envelope is radiative or convective, and if convective, on the composition structure of the star; simple approximate forms are presented for red and blue supergiant ejecta. Our models are most accurate for the high-velocity ejecta from the periphery of a star, where the shock dynamics are predictable. We present formulae for the final density distribution of this material, for both radiative and efficiently convective envelopes. These formulae limit to the well-known planar, self-similar solutions for mass shells approaching the stellar surface. But, the assumption of adiabatic flow fails at low optical depth, so this planar limit need not be attained. Formulae are given for the observable properties of the X-ray burst accompanying shock emergence, and their dependence on the parameters of the explosion. Motivated by the relativistic expansion recently inferred by Kulkarni et al. (1998) for the synchrotron shell around SN1998bw, we estimate the criterion for relativistic mass ejection and the rest mass of relativistic ejecta.

We describe a powerful methodology for numerical solution of 3-D self-gravitational hydrodynamics problems with extremely high resolution. Our method utilizes the technique of local adaptive mesh refinement (AMR), employing multiple grids at multiple levels of resolution. These grids are automatically and dynamically added and removed as necessary to maintain adequate resolution. This technology allows for the solution of problems in a manner that is both more efficient and more versatile than other fixed and variable resolution methods. The application of AMR to simulate the collapse and fragmentation of a molecular cloud, a key step in star formation, is discussed. Such simulations involve many orders of magnitude of variation in length scale as fragments form. In this paper we briefly describe the methodology and present an illustrative application for nonisothermal cloud collapse. We describe the numerical Jeans condition, a criterion for stability of self-gravitational hydrodynamics problems. We show the first well-resolved nonisothermal evolutionary sequence beginning with a perturbed dense molecular cloud core that leads to the formation of a binary system consisting of protostellar cores surrounded by distinct protostellar disks. The scale of the disks, of order 100 AU, is consistent with observations of gaseous disks surrounding single T-Tauri stars and debris disks surrounding systems such as $\beta$ Pictoris.

In this paper we discuss the propagation of dust through the interstellar medium (ISM), and describe the destructive effects of stellar winds, jets, and supernova shock waves on interstellar dust. We review the probability that grains formed in stellar outflows or supernovae survive processing in and propagation through the ISM, and incorporate themselves relatively unprocessed into meteoritic bodies in the solar system. We show that very large (radii >= 5 micron) and very small grains (radii <= 100 Angstrom) with sizes similar to the pre-solar SiC and diamond grains extracted from meteorites, can survive the passage through 100\kms shock waves relatively unscathed. High velocity (>= 250 km/s) shocks destroy dust efficiently. However, a small (~10%) fraction of the stardust never encountered such fast shocks before incorporation into the solar system. All grains should therefore retain traces of their passage through interstellar shocks during their propagation through the ISM. The grain surfaces should show evidence of processing due to sputtering and pitting due to small grain cratering collisions on the micron-sized grains. This conclusion seems to be in conflict with the evidence from the large grains recovered from meteorites which seem to show little interstellar processing.

SHORT VERSION: The rapidly rotating, masing circumnuclear disk in the central sub-parsec region of the galaxy NGC 4258 is remarkably circular and Keplerian, yet a striking asymmetry appears in the maser spectrum: the red-shifted, high- velocity sources are much more numerous and significantly more intense than the blue-shifted ones. A similar strong asymmetry appears also in the recently discovered, masing, circumnuclear disks in NGC 1068 and NGC 4945, thus suggesting it may be a general phenomenon. We show that the observed Doppler-shift asymmetry can naturally arise due to spiral shocks in circumnuclear disks, independent of the existence of a warp in the disk or the azimuthal direction to the observer. The high velocities of these features reflect the rotational velocities in the disk, and have nothing to do with the shock speed. In NGC 4258 - the currently most well-defined masing disk - the proposed scenario can also account for the intriguing clustering of the high-velocity maser spots in distinct clumps, the restricted spatial distribution of the low-velocity sources, and the dip in the maser spectrum at the systemic velocity of the disk. In this case we infer a disk mass of ~10E4 M_sun and a mass accretion rate of order ~7E-3 M_sun/year, which may be consistent with an advection-dominated accretion flow. The model is consistent with the observed Keplerian rotation, and introduces only negligible corrections to the previously derived black hole mass and galaxy distance. Predictions include slow systematic drifts in the velocity and position of all the high-velocity features, and the existence of circumnuclear disks which are delineated only by high-velocity maser emission.

We present an idealized, spherical model of the evolution of a magnetized molecular cloud due to ambipolar diffusion. This model allows us to follow the quasi-static evolution of the cloud's core prior to collapse and the subsequent evolution of the remaining envelope. By neglecting the thermal pressure gradients in comparison with magnetic stresses and by assuming that the ion velocity is small compared with the neutral velocity, we are able to find exact analytic solutions to the MHD equations. We show that, in the case of a centrally condensed cloud, a core of finite mass collapses into the origin leaving behind a quasi-static envelope, whereas initially homogeneous clouds never develop any structure in the absence of thermal stresses, and collapse as a whole. Prior to the collapse of the core, the cloud's evolution is characterized by two phases: a long, quasi-static phase where the relevant timescale is the ambipolar diffusion time (treated in this paper), and a short, dynamical phase where the characteristic timescale is the free-fall time. The collapse of the core is an "outside-in" collapse. The quasi-static evolution terminates when the cloud becomes magnetically supercritical; thereafter its evolution is dynamical, and a singularity develops at the origin-a protostar. After the initial formation of the protostar, the outer envelope continues to evolve quasi-statically, while the region of dynamical infall grows with time-an "inside-out" collapse. We use our solution to estimate the magnetic flux trapped in the collapsing core and the mass accretion rate onto the newly formed protostar. Our results agree, within factors of order unity, with the numerical results of Fiedler & Mouschovias (1992) for the physical quantities in the midplane of

We describe and report progress on the joint U.S. and German campaign to map the X-ray emission from the entire Cygnus Loop supernova remnant with the ROSAT High Resolution Imager. The Cygnus Loop is the prototype for a supernova remnant that is dominated by interactions with the interstellar medium and supplies fundamental physical information on this basic mechanism for shaping the interstellar medium. The global view that these high-resolution observations provide emphasizes the inhomogeneity of the interstellar medium and the pivotal nature of cloud-blast wave interactions in determining the X-ray morphology of the supernova remnant. While investigating the details of the evolution of the blast wave, we also describe the interstellar medium in the vicinity of the Cygnus Loop, which the progenitor star has processed. Although we do not expect the X-ray observations to be complete until September 1997, the incomplete data combined with deep H$\alpha$ images provide definitive evidence that the Cygnus Loop was formed by an explosion within a preexisting cavity.

We present the results of self-consistent models of Compton-heated accretion disk coronae. The models are calculated using a new method for computing monochromatic radiative transfer in two-dimensions. The method splits the radiation into direct and scattered components. The direct radiation is computed by calculating the optical depth along rays, while transfer of the scattered radiation is approximated by flux-limited diffusion. The resulting code agrees with more accurate treatments to within 50\%, and is highly efficient, making it practical for use in large hydrodynamic simulations. The coronal models are used to confirm the results of earlier work, and to extend it to higher luminosities. In contrast to earlier work, which found the outer disks to be shadowed by the inner corona at high luminosities, we find our results to form an almost continuous extension of the models at lower luminosities. This is due to the presence of multiply-scattered radiation, which acts to partially offset the loss of direct radiation from the central source. Although the analytic methods derived at lower luminosities cannot be used to derive the coronal structure for $L/L\sr{Edd}\gapprox0.1$, the results of the models are amenable to semiempirical fits. We also discuss possible observational consequences of the results for coronal veiling and line fluorescence from the disk.