Search

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

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

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

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

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

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

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

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

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

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

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

We review theoretical models of Population III.1 star formation, focusing on the protostellar feedback processes that are expected to terminate accretion and thus set the mass of these stars. We discuss how dark matter annihilation may modify this standard feedback scenario. Then, under the assumption that dark matter annihilation is unimportant, we predict the mass of stars forming in 12 cosmological minihalos produced in independent numerical simulations. This allows us to make a simple estimate of the Pop III.1 initial mass function and how it may evolve with redshift., Comment: 6 pages, Proceedings of 'The First Stars and Galaxies: Challenges for the Next Decade", Austin, TX, March 8-11, 2010

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

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

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

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

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

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

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

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

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

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

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., Comment: 16 pages, including 11 figures, Review paper to appear in "First Stars III", eds. B. O'Shea, A. Heger and T. Abel

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), Comment: 31 pages, including 10 figures, accepted to ApJ

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

Supermassive black hole accretion and star formation appear intimately connected. I review the observational and theoretical evidence for this statement. I then discuss how focussed studies of two systems, our Galactic Center and the nucleus of M87, can help to improve our understanding of these processes., Comment: To appear in proceedings of "First Light Science with the GTC", eds. R. Guzman, C. Packham, & J.M. Rodriguez-Espinoza

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

I review the status of massive star formation theories: accretion from collapsing, massive, turbulent cores; competitive accretion; and stellar collisions. I conclude the observational and theoretical evidence favors the first of these models. I then discuss: the initial conditions of star cluster formation as traced by infrared dark clouds; the cluster formation timescale; and comparison of the initial cluster mass function in different galactic environments., Comment: 7 pages, invited paper in proceedings of IAU Symposium 237, "Triggered Star Formation in a Turbulent ISM", eds. B. G. Elmegreen & J. Palous

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

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., Comment: 5 pages, Accepted to ApJ Letters

I review some important questions in the field of massive star formation: What are the initial conditions for proto star clusters and how do they arise? What are the initial conditions for individual massive star formation within star clusters? How do massive protostars accumulate their mass? I compare the Turbulent Core Model (McKee & Tan 2003) to several nearby regions, including Orion KL. Here I also discuss the origin of BN's high proper motion., Comment: 10 pages, 2 figures, invited review, to appear in "Massive Star Birth - A Crossroads of Astrophysics" (CUP), eds. R. Cesaroni, E. Churchwell, M. Felli, and C. M. Walmsley

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

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

I discuss two questions about the origin of star clusters: How long does the process take? What is the mode of individual star formation? I argue that observations of Galactic star-forming regions, particularly the Orion Nebula Cluster (ONC), indicate that cluster formation often takes several Myr, which is many local dynamical timescales. Individual stars and binaries, including massive stars, appear to form from the collapse of gas cores., Comment: 6 pages, to appear in the proceedings of the conference "The Young Local Universe", La Thuile (Italy), March 2004; This talk and others are available at http://www-laog.obs.ujf-grenoble.fr/ylu/ylu_vgr/vgr_index.html

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., Comment: 8 pages, invited contribution for the Cancun Workshop on Formation and Evolution of Young Massive Clusters, eds. H. J. G. L. M. Lamers, A. Nota and L. Smith

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

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

I compare theoretical models of massive star formation with observations of the Orion Hot Core, which harbors one of the closest massive protostars. Although this region is complicated, many of its features (size, luminosity, accretion disk, HII region, outflow) may be understood by starting with a simple model in which the star forms from a massive gas core that is a coherent entity in approximate pressure balance with its surroundings. The dominant contribution to the pressure is from turbulent motions and magnetic fields. The collapse can be perturbed by interactions with other stars in the forming cluster, which may induce sporadic enhancements of the accretion and outflow rate., Comment: 8 pages, in proceedings of Star Formation in the Interstellar Medium, workshop honoring David Hollenbach, Chris McKee, & Frank Shu, eds. F.C. Adams, D. Johnstone, D.N.C. Lin, & E.C. Ostriker, ASP. An expanded overview of observations of the Orion Hot Core region can be found at http://astro.princeton.edu/~jt/orion/orionhotcore.html

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

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., Comment: IAU Symposium 221, Star Formation at High Angular Resolution, Eds. M. Burton, R. Jayawardhana, & T. Bourke, electronically published at http://www.phys.unsw.edu.au/iau221

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

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., Comment: 20 pages, Accepted to ApJ, some re-arrangement of text for improved clarity

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

Magnetically mediated disk outflows are a leading paradigm to explain winds and jets in a variety of astrophysical sources, but where do the fields come from? Since accretion of mean magnetic flux may be disfavored in a thin turbulent disk, and only fields generated with sufficiently large scale can escape before being shredded by turbulence, in situ field production is desirable. Nonlinear helical inverse dynamo theory can provide the desired fields for coronae and outflows. We discuss the implications for contemporary protostellar disks, where the MRI (magneto-rotational instability) can drive turbulence in the inner regions, and primordial protostellar disks, where gravitational instability drives the turbulence. We emphasize that helical dynamos are compatible with the magneto-rotational instability, and clarify the relationship between the two., Comment: 12 pages, LaTeX (with figs); version in press for "Proceedings of the International Workshop on Magnetic Fields and Star Formation: Theory vs. Observation" Madrid, Spain; April 2003

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., Comment: 4 pages; To appear in proceedings of the 13th Annual Astrophysics Conference in Maryland: The Emergence of Cosmic Structure, eds. S. Holt and C. Reynolds, (AIP)

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

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., Comment: ApJ accepted; 28 pages, some clarification of the text, results unchanged

I review the arguments motivating models for massive star formation via stellar collisions. I then describe how the standard accretion scenario, involving the collapse of a quasi-hydrostatic gas core, can produce high-mass stars in the pressurized regions of forming star clusters. I argue that the observational evidence, particularly in the Orion hot core, favors the standard accretion paradigm., Comment: 10 pages. Invited talk at the IAOC Workshop on "Galactic Star Formation Across the Stellar Mass Spectrum" held in La Serena, Chile, March 2002; to be published in the ASP Conference Series, edited by J. M. De Buizer