Neal J. Evans

Star Formation in the Milky Way and Nearby Galaxies

We review progress over the past decade in observations of large-scale star formation, with a focus on the interface between extragalactic and Galactic studies. Methods of measuring gas contents and star-formation rates are discussed, and updated prescriptions for calculating star-formation rates are provided. We review relations between star formation and gas on scales ranging from entire galaxies to individual molecular clouds.

Physical Conditions in Regions of Star Formation

▪ Abstract The physical conditions in molecular clouds control the nature and rate of star formation, with consequences for planet formation and galaxy evolution. The focus of this review is on the conditions that characterize regions of star formation in our Galaxy. A review of the tools and tracers for probing physical conditions includes summaries of generally applicable results. Further discussion distinguishes between the formation of low-mass stars in relative isolation and formation in a clustered environment. Evolutionary scenarios and theoretical predictions are more developed for isolated star formation, and observational tests are beginning to interact strongly with the theory. Observers have identified dense cores collapsing to form individual stars or binaries, and analysis of some of these cores support theoretical models of collapse. Stars of both low and high mass form in clustered environments, but massive stars form almost exclusively in clusters. The theoretical understanding of such re...

From Molecular Cores to Planet‐forming Disks: An SIRTF Legacy Program

Crucial steps in the formation of stars and planets can be studied only at mid‐ to far‐infrared wavelengths, where the Space Infrared Telescope (SIRTF) provides an unprecedented improvement in sensitivity. We will use all three SIRTF instruments (Infrared Array Camera [IRAC], Multiband Imaging Photometer for SIRTF [MIPS], and Infrared Spectrograph [IRS]) to observe sources that span the evolutionary sequence from molecular cores to protoplanetary disks, encompassing a wide range of cloud masses, stellar masses, and star‐forming environments. In addition to targeting about 150 known compact cores, we will survey with IRAC and MIPS (3.6–70 μm) the entire areas of five of the nearest large molecular clouds for new candidate protostars and substellar objects as faint as 0.001 solar luminosities. We will also observe with IRAC and MIPS about 190 systems likely to be in the early stages of planetary system formation (ages up to about 10 Myr), probing the evolution of the circumstellar dust, the raw material for planetary cores. Candidate planet‐forming disks as small as 0.1 lunar masses will be detectable. Spectroscopy with IRS of new objects found in the surveys and of a select group of known objects will add vital information on the changing chemical and physical conditions in the disks and envelopes. The resulting data products will include catalogs of thousands of previously unknown sources, multiwavelength maps of about 20 deg^2 of molecular clouds, photometry of about 190 known young stars, spectra of at least 170 sources, ancillary data from ground‐based telescopes, and new tools for analysis and modeling. These products will constitute the foundations for many follow‐up studies with ground‐based telescopes, as well as with SIRTF itself and other space missions such as SIM, JWST, Herschel, and TPF/Darwin.

Tracing the mass during low-mass star formation. II. Modeling the submillimeter emission from preprotostellar cores

We have modeled the emission from dust in preprotostellar cores, including a self-consistent calculation of the temperature distribution for each input density distribution. Model density distributions include Bonnor-Ebert spheres and power laws. The Bonnor-Ebert spheres fit the data well for all three cores that we have modeled. The dust temperatures decline to very low values (Td ~ 7 K) in the centers of these cores, strongly affecting the dust emission. Compared to earlier models that assume constant dust temperatures, our models indicate higher central densities and smaller regions of relatively constant density. Indeed, for L1544, a power-law density distribution, similar to that of a singular, isothermal sphere, cannot be ruled out. For the three sources modeled herein, there seems to be a sequence of increasing central condensation, from L1512 to L1689B to L1544. The two denser cores, L1689B and L1544, have spectroscopic evidence for contraction, suggesting an evolutionary sequence for preprotostellar cores.

The Spitzer Ice Legacy: Ice Evolution from Cores to Protostars

Ices regulate much of the chemistry during star formation and account for up to 80% of the available oxygen and carbon. In this paper, we use the Spitzer c2d Legacy ice survey, complimented with data sets on ices in cloud cores and high-mass protostars, to determine standard ice abundances and to present a coherent picture of the evolution of ices during low- and high-mass star formation. The median ice composition H_(2)O:CO:CO_2:CH_(3)OH:NH_3:CH_4:XCN is 100:29:29:3:5:5:0.3 and 100:13:13:4:5:2:0.6 toward low- and high-mass protostars, respectively, and 100:31:38:4:-:-:- in cloud cores. In the low-mass sample, the ice abundances with respect to H_(2)O of CH_4, NH_3, and the component of CO_2 mixed with H_(2)O typically vary by <25%, indicative of co-formation with H_(2)O. In contrast, some CO and CO_2 ice components, XCN, and CH3OH vary by factors 2-10 between the lower and upper quartile. The XCN band correlates with CO, consistent with its OCN– identification. The origin(s) of the different levels of ice abundance variations are constrained by comparing ice inventories toward different types of protostars and background stars, through ice mapping, analysis of cloud-to-cloud variations, and ice (anti-)correlations. Based on the analysis, the first ice formation phase is driven by hydrogenation of atoms, which results in an H_(2)O-dominated ice. At later prestellar times, CO freezes out and variations in CO freezeout levels and the subsequent CO-based chemistry can explain most of the observed ice abundance variations. The last important ice evolution stage is thermal and UV processing around protostars, resulting in CO desorption, ice segregation, and the formation of complex organic molecules. The distribution of cometary ice abundances is consistent with the idea that most cometary ices have a protostellar origin.

THE STAR FORMATION RATE AND GAS SURFACE DENSITY RELATION IN THE MILKY WAY: IMPLICATIONS FOR EXTRAGALACTIC STUDIES

We investigate the relation between star formation rate (SFR) and gas surface densities in Galactic star-forming regions using a sample of young stellar objects (YSOs) and massive dense clumps. Our YSO sample consists of objects located in 20 large molecular clouds from the Spitzer cores to disks (c2d) and Goulds Belt (GB) surveys. These data allow us to probe the regime of low-mass star formation, essentially invisible to tracers of high-mass star formation used to establish extragalactic SFR-gas relations. We estimate the gas surface density (Σgas) from extinction (AV ) maps and YSO SFR surface densities (ΣSFR) from the number of YSOs, assuming a mean mass and lifetime. We also divide the clouds into evenly spaced contour levels of AV , counting only Class I and Flat spectral energy distribution YSOs, which have not yet migrated from their birthplace. For a sample of massive star-forming clumps, we derive SFRs from the total infrared luminosity and use HCN gas maps to estimate gas surface densities. We find that c2d and GB clouds lie above the extragalactic SFR-gas relations (e.g., Kennicutt-Schmidt law) by factors of up to 17. Cloud regions with high Σgas lie above extragalactic relations up to a factor of 54 and overlap with high-mass star-forming regions. We use 12CO and 13CO gas maps of the Perseus and Ophiuchus clouds from the COMPLETE survey to estimate gas surface densities and compare to measurements from AV maps. We find that 13CO, with the standard conversions to total gas, underestimates the AV -based mass by factors of ~4-5. 12CO may underestimate the total gas mass at Σgas 200 M ☉ pc–2 by 30%; however, this small difference in mass estimates does not explain the large discrepancy between Galactic and extragalactic relations. We find evidence for a threshold of star formation (Σth) at 129 ± 14 M ☉ pc–2. At Σgas>Σth, the Galactic SFR-gas relation is linear. A possible reason for the difference between Galactic and extragalactic relations is that much of Σgas is below Σth in extragalactic studies, which detect all the CO-emitting gas. If the Kennicutt-Schmidt relation (ΣSFR Σ1.4 gas) and a linear relation between dense gas and star formation are assumed, the fraction of dense star-forming gas (f dense) increases as ~Σ0.4 gas. When Σgas reaches ~300 Σth, the fraction of dense gas is ~1, creating a maximal starburst.

Evidence for protostellar collapse in B335

We have observed five rotational transitions of H2CO and CS toward the Bok globule, B335, with high spatial and spectral resolution. The characteristic shape of the observed profiles provides direct, kinematic evidence of collapse. In addition, we have modeled line profiles of collapsing dense cores with density and velocity structures taken from the theory of Shu and coworkers. Using the age of collapse as the only free parameter, we found that the strengths and profiles of the observed lines can be well fitted by the theoretical model. Our best-fit model gives an age of 1.5 x 10 exp 5 yr, corresponding to an infall radius of 0.04 pc and a total mass of 0.4 solar mass for the central star and disk. Outside the infall radius, there is a static envelope with a r exp -2 density distribution, an average temperature of 13 K, and a turbulent velocity (1/e width) of 0.12 km/s. The CS abundance is 3.6 x 10 exp -9 with about 30 percent uncertainty.

The Mass Distribution And Lifetime Of Prestellar Cores In Perseus, Serpens, And Ophiuchus

We present an unbiased census of starless cores in Perseus, Serpens, and Ophiuchus, assembled by comparing large-scale Bolocam 1.1 mm continuum emission maps with Spitzer c2d surveys. We use the c2d catalogs to separate 108 starless from 92 protostellar cores in the 1.1 mm core samples from Enoch and Young and their coworkers. A comparison of these populations reveals the initial conditions of the starless cores. Starless cores in Perseus have similar masses but larger sizes and lower densities on average than protostellar cores, with sizes that suggest density profiles substantially flatter than ρ∝r^-2. By contrast, starless cores in Serpens are compact and have lower masses than protostellar cores; future star formation will likely result in lower mass objects than the currently forming protostars. Comparison to dynamical masses estimated from the NH3 survey of Perseus cores by Rosolowsky and coworkers suggests that most of the starless cores are likely to be gravitationally bound, and thus prestellar. The combined prestellar core mass distribution includes 108 cores and has a slope of α = -2.3 ± 0.4 for M > 0.8 M☉. This slope is consistent with recent measurements of the stellar initial mass function, providing further evidence that stellar masses are directly linked to the core formation process. We place a lower limit on the core-to-star efficiency of 25%. There are approximately equal numbers of prestellar and protostellar cores in each cloud; thus the dense prestellar core lifetime must be similar to the lifetime of embedded protostars, or 4.5 x 10^5 yr, with a total uncertainty of a factor of 2. Such a short lifetime suggests a dynamic, rather than quasi-static, core evolution scenario, at least at the relatively high mean densities (n > 2 x 10^4 cm^-3) to which we are sensitive.

Dense gas and star formation: characteristics of cloud cores associated with water masers

We have observed 150 regions of massive star formation, selected originally by the presence of an H2O maser, in the J = 5 → 4, 3 → 2, and 2 → 1 transitions of CS, and 49 regions in the same transitions of C34S. Over 90% of the 150 regions were detected in the J = 2 → 1 and 3 → 2 transitions of CS, and 75% were detected in the J = 5 → 4 transition. We have combined the data with the J = 7 → 6 data from our original (1992) survey to determine the density by analyzing the excitation of the rotational levels. Using large velocity gradient models, we have determined densities and column densities for 71 of these regions. The gas densities are very high (log n = 5.9), but much less than the critical density of the J = 7 → 6 line. Small maps of 25 of the sources in the J = 5 → 4 line yield a mean diameter of 1.0 pc. Several estimates of the mass of dense gas were made for the sources for which we had sufficient information. The mean virial mass is 3800 M☉. The mean ratio of bolometric luminosity to virial mass (L/M) is 190, about 50 times higher than estimates made using CO emission, suggesting that star formation is much more efficient in the dense gas probed in this study. The depletion time for the dense gas is ~1.3 × 107 yr, comparable to the timescale for gas dispersal around open clusters and OB associations. We find no statistically significant line width-size or density-size relationships in our data. Instead, both line width and density are greater for a given size than would be predicted by the usual relationships. We find that the line width increases with density, the opposite of what would be predicted by the usual arguments. We estimate that the luminosity of our Galaxy (excluding the inner 400 pc) in the CS J = 5 → 4 transition is 15-23 L☉, considerably less than the luminosity in this line within the central 100 pc of NGC 253 and M82. In addition, the ratio of far-infrared luminosity to CS luminosity is higher in M82 than in any cloud in our sample.

c2d Spitzer IRS Spectra of Disks around T Tauri Stars. I. Silicate Emission and Grain Growth

Infrared ~5-35 μm spectra for 40 solar mass T Tauri stars and 7 intermediate-mass Herbig Ae stars with circumstellar disks were obtained using the Spitzer Space Telescope as part of the c2d IRS survey. This work complements prior spectroscopic studies of silicate infrared emission from disks, which were focused on intermediate-mass stars, with observations of solar mass stars limited primarily to the 10 μm region. The observed 10 and 20 μm silicate feature strengths/shapes are consistent with source-to-source variations in grain size. A large fraction of the features are weak and flat, consistent with micron-sized grains indicating fast grain growth (from 0.1 to 1.0 μm in radius). In addition, approximately half of the T Tauri star spectra show crystalline silicate features near 28 and 33 μm, indicating significant processing when compared to interstellar grains. A few sources show large 10-to-20 μm ratios and require even larger grains emitting at 20 μm than at 10 μm. This size difference may arise from the difference in the depth into the disk probed by the two silicate emission bands in disks where dust settling has occurred. The 10 μm feature strength versus shape trend is not correlated with age or Hα equivalent width, suggesting that some amount of turbulent mixing and regeneration of small grains is occurring. The strength versus shape trend is related to spectral type, however, with M stars showing significantly flatter 10 μm features (larger grain sizes) than A/B stars. The connection between spectral type and grain size is interpreted in terms of the variation in the silicate emission radius as a function of stellar luminosity, but could also be indicative of other spectral-type-dependent factors (e.g., X-rays, UV radiation, and stellar/disk winds).