The physics and modes of star cluster formation: simulations

We review progress in numerical simulations of star cluster formation. These simulations involve the bottom-up assembly of clusters through hierarchical mergers, which produces a fractal stellar distribution at young (~0.5 Myr) ages. The resulting clusters are predicted to be mildly aspherical and highly mass-segregated, except in the immediate aftermath of mergers. The upper initial mass function within individual clusters is generally somewhat flatter than for the aggregate population. Recent work has begun to clarify the factors that control the mean stellar mass in a star-forming cloud and also the efficiency of star formation. The former is sensitive to the thermal properties of the gas while the latter depends both on the magnetic field and the initial degree of gravitational boundedness of the natal cloud. Unmagnetized clouds that are initially bound undergo rapid collapse, which is difficult to reverse by ionization feedback or stellar winds.Comment: 21 pages, 10 figures. To appear as invited review article in a special issue of the Phil. Trans. Royal Soc. A: Ch. 3 "Star clusters as tracers of galactic star-formation histories" (ed. R. de Grijs). Fully peer reviewed. LaTeX, requires rspublic.cls style fil

Spectroscopic Detection of a Stellar-like Photosphere in an Accreting Protostar

We present the first spectrum of a highly veiled, strongly accreting protostar which shows photospheric absorption features and demonstrates the stellar nature of its central core. We find the spectrum of the luminous (L_bol = 10 L_sun) protostellar source, YLW 15, to be stellar-like with numerous atomic and molecular absorption features, indicative of a K5 IV/V spectral type and a continuum veiling r_k = 3.0. Its derived stellar luminosity (3 L_sun) and stellar radius (3.1 R_sun) are consistent with those of a 0.5 M_sun pre-main-sequence star. However, 70% of its bolometric luminosity is due to mass accretion, whose rate we estimate to be 1.6 E-6 M_sun / yr onto the protostellar core. We determine that excess infrared emission produced by the circumstellar accretion disk, the inner infalling envelope, and accretion shocks at the surface of the stellar core of YLW 15 all contribute signifi- cantly to its near-IR continuum veiling. Its projected rotation velocity v sin i = 50 km / s is comparable to those of flat-spectrum protostars but considerably higher than those of classical T Tauri stars in the rho Oph cloud. The protostar may be magnetically coupled to its circumstellar disk at a radius of 2 R_*. It is also plausible that this protostar can shed over half its angular momentum and evolve into a more slowly rotating classical T Tauri star by remaining coupled to its circumstellar disk (at increasing radius) as its accretion rate drops by an order of magnitude during the rapid transition between the Class I and Class II phases of evolution. The spectrum of WL 6 does not show any photospheric absorption features, and we estimate that its continuum veiling is r_k >= 4.6. Together with its low bolometric luminosity (2 L_sun), this dictates that its central core is very low mass, ~0.1 M_sun.Comment: 14 pages including 9 figures (3 figures of 3 panels each, all as separate files). AASTeX LaTex macros version 5.0. To be published in The Astronomical Journal (tentatively Oct 2002

Radiative Transfer in Prestellar Cores: A Monte Carlo Approach

We use our Monte Carlo radiative transfer code to study non-embedded prestellar cores and cores that are embedded at the centre of a molecular cloud. Our study indicates that the temperature inside embedded cores is lower than in isolated non-embedded cores, and generally less than 12 K, even when the cores are surrounded by an ambient cloud of small visual extinction (Av~5). Our study shows that the best wavelength region to observe embedded cores is between 400 and 500 microns, where the core is quite distinct from the background. We also predict that very sensitive observations (~1-3 MJy/sr) at 170-200 microns can be used to estimate how deeply a core is embedded in its parent molecular cloud. Finally, we present preliminary results of asymmetric models of non-embedded cores.Comment: 8 pages, 15 figures, to appear in the conference proceedings of "Open Issues in Local Star Formation and Early Stellar Evolution", held in Ouro Preto (Brazil), April 5-10, 200

Sensitive Limits on the Water Abundance in Cold Low Mass Molecular Cores

We present SWAS observations of water vapor in two cold star-less clouds, B68 and Core D in rho Ophiuchus. Sensitive non-detections of the 1(10)-1(01) transition of o-H2O are reported for each source. Both molecular cores have been previously examined by detailed observations that have characterized the physical structure. Using these rather well defined physical properties and a Monte-Carlo radiation transfer model we have removed one of the largest uncertainties from the abundance calculation and set the lowest water abundance limit to date in cold low-mass molecular cores. These limits are < 3 x 10^{-8} (relative to H2) and < 8 x 10^{-9} in B68 and rho Oph D, respectively. Such low abundances confirm the general lack of ortho-water vapor in cold (T < 20 K) cores. Provided that the ortho/para ratio of water is not near zero, these limits are well below theoretical predictions and appear to support the suggestion that most of the water in dense low-mass cores is frozen onto the surfaces of cold dust grains.Comment: 12 pages, 3 figures, accepted by Astrophysical Journal Letter

The Earliest Phases Of High-Mass Star Formation, As Seen In Ngc 6334 By Herschel -Hobys

Aims. To constrain models of high-mass star formation, the Herschel-HOBYS key program aims at discovering massive dense cores (MDCs) able to host the high-mass analogs of low-mass prestellar cores, which have been searched for over the past decade. We here focus on NGC 6334, one of the best-studied HOBYS molecular cloud complexes

The "Mysterious" Origin of Brown Dwarfs

Hundreds of brown dwarfs (BDs) have been discovered in the last few years in stellar clusters and among field stars. BDs are almost as numerous as hydrogen burning stars and so a theory of star formation should also explain their origin. The ``mystery'' of the origin of BDs is that their mass is two orders of magnitude smaller than the average Jeans' mass in star--forming clouds, and yet they are so common. In this work we investigate the possibility that gravitationally unstable protostellar cores of BD mass are formed directly by the process of turbulent fragmentation. Supersonic turbulence in molecular clouds generates a complex density field with a very large density contrast. As a result, a fraction of BD mass cores formed by the turbulent flow are dense enough to be gravitationally unstable. We find that with density, temperature and rms Mach number typical of cluster--forming regions, turbulent fragmentation can account for the observed BD abundance.Comment: 11 pages, 3 figures, ApJ submitted Error in equation 1 has been corrected. Improved figure

Molecular Tracers of Embedded Star Formation in Ophiuchus

In this paper we analyze nine SCUBA cores in Ophiuchus using the second-lowest rotational transitions of four molecular species (12CO, 13CO, C18O, and C17O) to search for clues to the evolutionary state and star-formation activity within each core. Specifically, we look for evidence of outflows, infall, and CO depletion. The line wings in the CO spectra are used to detect outflows, spectral asymmetries in 13CO are used to determine infall characteristics, and a comparison of the dust emission (from SCUBA observations) and gas emission (from C18O) is used to determine the fractional CO freeze-out. Through comparison with Spitzer observations of protostellar sources in Ophiuchus, we discuss the usefulness of CO and its isotopologues as the sole indicators of the evolutionary state of each core. This study is an important pilot project for the JCMT Legacy Survey of the Gould Belt (GBS) and the Galactic Plane (JPS), which intend to complement the SCUBA-2 dust continuum observations with HARP observations of 12CO, 13CO, C18O, and C17O J = 3 - 2 in order to determine whether or not the cold dust clumps detected by SCUBA-2 are protostellar or starless objects. Our classification of the evolutionary state of the cores (based on molecular line maps and SCUBA observations) is in agreement with the Spitzer designation for six or seven of the nine SCUBA cores. However, several important caveats exist in the interpretation of these results, many of which large mapping surveys like the GBS may be able to overcome to provide a clearer picture of activity in crowded fields.Comment: 43 pages including 19 postscript figures. Accepted for publication in the PAS

Shapes of Molecular Cloud Cores and the Filamentary Mode of Star Formation

Using recent dust continuum data, we generate the intrinsic ellipticity distribution of dense, starless molecular cloud cores. Under the hypothesis that the cores are all either oblate or prolate randomly-oriented spheroids, we show that a satisfactory fit to observations can be obtained with a gaussian prolate distribution having a mean intrinsic axis ratio of 0.54. Further, we show that correlations exist between the apparent axis ratio and both the peak intensity and total flux density of emission from the cores, the sign of which again favours the prolate hypothesis. The latter result shows that the mass of a given core depends on its intrinsic ellipticity. Monte Carlo simulations are performed to find the best-fit power law of this dependence. Finally, we show how these results are consistent with an evolutionary scenario leading from filamentary parent clouds to increasingly massive, condensed, and roughly spherical embedded cores.Comment: 16 pages, incl. 11 Postscript figures. Accepted by Ap

The Initial Mass Function of Low-Mass Stars and Brown Dwarfs in Taurus

By combining deep optical imaging and infrared spectroscopy with data from the Two-Micron All-Sky Survey (2MASS) and from previous studies (e.g., Briceno et al.), I have measured the Initial Mass Function (IMF) for a reddening-limited sample in four fields in the Taurus star forming region. This IMF is representative of the young populations within these fields for masses above 0.02 Msun. Relative to the similarly derived IMF for the Trapezium Cluster (Luhman et al.), the IMF for Taurus exhibits a modest deficit of stars above one solar mass (i.e., steeper slope), the same turnover mass (~0.8 Msun), and a significant deficit of brown dwarfs. If the IMF in Taurus were the same as that in the Trapezium, 12.8+/-1.8 brown dwarfs (>0.02 Msun) are expected in these Taurus fields where only one brown dwarf candidate is found. These results are used to test theories of the IMF.Comment: to be published in The Astrophysical Journal, 24 pages, 6 figures, also found at http://cfa-www.harvard.edu/~kluhman/taurus

Large Area Mapping at 850 Microns. V. Analysis of the Clump Distribution in the Orion A South Molecular Cloud

We present results from a 2300 arcmin^2 survey of the Orion A molecular cloud at 450 and 850 micron using the Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope. The region mapped lies directly south of the OMC1 cloud core and includes OMC4, OMC5, HH1/2, HH34, and L1641N. We identify 71 independent clumps in the 850 micron map and compute size, flux, and degree of central concentration in each. Comparison with isothermal, pressure-confined, self-gravitating Bonnor-Ebert spheres implies that the clumps have internal temperatures T_d ~ 22 +/- K and surface pressures log (k^-1 P cm^-3 K) = 6.0 +/- 0.2. The clump masses span the range 0.3 - 22 Msun assuming a dust temperature T_d ~ 20 K and a dust emissivity kappa_850 = 0.02 cm^2 g^-1. The distribution of clump masses is well characterized by a power-law N(M) propto M^-alpha with alpha = 2.0 +/- 0.5 for M > 3.0 Msun, indicating a clump mass function steeper than the stellar Initial Mass Function. Significant incompleteness makes determination of the slope at lower masses difficult. A comparison of the submillimeter emission map with an H_2 2.122 micron survey of the same region is performed. Several new Class 0 sources are revealed and a correlation is found between both the column density and degree of concentration of the submillimeter sources and the likelihood of coincident H_2 shock emission.Comment: 44 pages, 17 figures, accepted by Ap