Formation of the First Stars by Accretion

The process of star formation from metal-free gas is investigated by following the evolution of accreting protostars with emphasis on the properties of massive objects. The main aim is to establish the physical processes that determine the upper mass limit of the first stars. Although the consensus is that massive stars were commonly formed in the first cosmic structures, our calculations show that their actual formation depends sensitively on the mass accretion rate and its time variation. Even in the rather idealized case in which star formation is mainly determined by dot{M}acc, the characteristic mass scale of the first stars is rather uncertain. We find that there is a critical mass accretion rate dot{M}crit = 4 10^{-3} Msun/yr that separates solutions with dot{M}acc> 100 Msun can form, provided there is sufficient matter in the parent clouds, from others (dot{M}acc > dot{M}crit) where the maximum mass limit decreases as dot{M}acc increases. In the latter case, the protostellar luminosity reaches the Eddington limit before the onset of hydrogen burning at the center via the CN-cycle. This phase is followed by a rapid and dramatic expansion of the radius, possibly leading to reversal of the accretion flow when the stellar mass is about 100Msun. (abridged)Comment: 34 pages, 12 figures. ApJ, in pres

Thermal and Fragmentation Properties of Star-forming Clouds in Low-metallicity Environments

The thermal and chemical evolution of star-forming clouds is studied for different gas metallicities, Z, using the model of Omukai (2000), updated to include deuterium chemistry and the effects of cosmic microwave background (CMB) radiation. HD-line cooling dominates the thermal balance of clouds when Z \~ 10^{-5}-10^{-3} Z_sun and density ~10^{5} cm^{-3}. Early on, CMB radiation prevents the gas temperature to fall below T_CMB, although this hardly alters the cloud thermal evolution in low-metallicity gas. From the derived temperature evolution, we assess cloud/core fragmentation as a function of metallicity from linear perturbation theory, which requires that the core elongation E := (b-a)/a > E_NL ~ 1, where a (b) is the short (long) core axis length. The fragment mass is given by the thermal Jeans mass at E = E_NL. Given these assumptions and the initial (gaussian) distribution of E we compute the fragment mass distribution as a function of metallicity. We find that: (i) For Z=0, all fragments are very massive, > 10^{3}M_sun, consistently with previous studies; (ii) for Z>10^{-6} Z_sun a few clumps go through an additional high density (> 10^{10} cm^{-3}) fragmentation phase driven by dust-cooling, leading to low-mass fragments; (iii) The mass fraction in low-mass fragments is initially very small, but at Z ~ 10^{-5}Z_sun it becomes dominant and continues to grow as Z is increased; (iv) as a result of the two fragmentation modes, a bimodal mass distribution emerges in 0.01 0.1Z_sun, the two peaks merge into a singly-peaked mass function which might be regarded as the precursor of the ordinary Salpeter-like IMF.Comment: 38 pages, 16 figures, ApJ in pres

Dynamics of Dense Cores in the Perseus Molecular Cloud

We survey the kinematics of over one hundred and fifty candidate (and potentially star-forming) dense cores in the Perseus molecular cloud with pointed N2H+(1-0) and simultaneous C18O(2-1) observations. Our detection rate of N2H+ is 62%, rising to 84% for JCMT SCUBA-selected targets. In agreement with previous observations, we find that the dense N2H+ targets tend to display nearly thermal linewidths, particularly those which appear to be starless (using Spitzer data), indicating turbulent support on the small scales of molecular clouds is minimal. For those N2H+ targets which have an associated SCUBA dense core, we find their internal motions are more than sufficient to provide support against the gravitational force on the cores. Comparison of the N2H+ integrated intensity and SCUBA flux reveals fractional N2H+ abundances between 10^-10 and 10^-9. We demonstrate that the relative motion of the dense N2H+ gas and the surrounding C18O gas is less than the sound speed in the vast majority of cases (~90%). The point-to-point motions we observe within larger extinction regions appear to be insufficient to provide support against gravity, although we sparsely sample these regions.Comment: 49 pages, 20 figures. Accepted for publication in the Astrophysical Journa

Massive stars and globular cluster formation

We first present chemodynamical simulations to investigate how stellar winds of massive stars influence early dynamical and chemical evolution of forming globular clusters (GCs). In our numerical models, GCs form in turbulent,high-density giant molecular clouds (GMCs), which are embedded in a massive dark matter halo at high redshifts. We show how high-density, compact stellar systems are formed from GMCs influenced both by physical processes associated with star formation and by tidal fields of their host halos. We also show that chemical pollution of GC-forming GMCs by stellar winds from massive stars can result in star-to-star abundance inhomogeneities among light elements (e.g., C, N, and O) of stars in GCs. The present model with a canonical initial mass function (IMF) also shows a C-N anticorrelation that stars with smaller [C/Fe] have larger [N/Fe] in a GC. Although these results imply that ``self-pollution'' of GC-forming GMCs by stellar winds from massive stars can cause abundance inhomogeneities of GCs, the present models with different parameters and canonical IMFs can not show N-rich stars with [N/Fe] ~ 0.8 observed in some GCs (e.g., NGC 6752). We discuss this apparent failure in the context of massive star formation preceding low-mass one within GC-forming GMCs (``bimodal star formation scenario''). We also show that although almost all stars (~97%) show normal He abundances (Y) of ~0.24 some stars later formed in GMCs can have Y as high as ~0.3 in some models. The number fraction of He-rich stars with Y >0.26 is however found to be small (~10^-3) for most models.Comment: 10 pages, 8 figures, accepted by Ap

Evolution of Angular Momentum Distribution during Star Formation

If the angular momentum of the molecular cloud core were conserved during the star formation process, a new-born star would rotate much faster than its fission speed. This constitutes the angular momentum problem of new-born stars. In this paper, the angular momentum transfer in the contraction of a rotating magnetized cloud is studied with axisymmetric MHD simulations. Owing to the large dynamic range covered by the nested-grid method, the structure of the cloud in the range from 10 AU to 0.1 pc is explored. First, the cloud experiences a run-away collapse, and a disk forms perpendicularly to the magnetic field, in which the central density increases greatly in a finite time-scale. In this phase, the specific angular momentum j of the disk decreases to $\simeq 1/3$ of the initial cloud. After the central density of the disk exceeds $\sim 10^{10}{\rm cm}^{-3}$, the infall on to the central object develops. In this accretion stage, the rotation motion and thus the toroidal magnetic field drive the outflow. The angular momentum of the central object is transferred efficiently by the outflow as well as the effect of the magnetic stress. In 7000 yr from the core formation, the specific angular momentum of the central $0.17M_\odot$ decreases a factor of 10^{-4} from the initial value (i.e. from $10^{20}{\rm cm^2 s^{-1}}$ to $10^{16}{\rm cm^2 s^{-1}}$).Comment: 15 pages, 2 figures, Astrophysical Journal Letters in pres

Bayesian Power Spectrum Analysis of the First-Year WMAP data

We present the first results from a Bayesian analysis of the WMAP first year data using a Gibbs sampling technique. Using two independent, parallel supercomputer codes we analyze the WMAP Q, V and W bands. The analysis results in a full probabilistic description of the information the WMAP data set contains about the power spectrum and the all-sky map of the cosmic microwave background anisotropies. We present the complete probability distributions for each C_l including any non-Gaussianities of the power spectrum likelihood. While we find good overall agreement with the previously published WMAP spectrum, our analysis uncovers discrepancies in the power spectrum estimates at low l multipoles. For example we claim the best-fit Lambda-CDM model is consistent with the C_2 inferred from our combined Q+V+W analysis with a 10% probability of an even larger theoretical C_2. Based on our exact analysis we can therefore attribute the "low quadrupole issue" to a statistical fluctuation.Comment: 5 pages. 4 figures. For additional information and data see http://www.astro.uiuc.edu/~iodwyer/research#wma

Curvature and Acoustic Instabilities in Rotating Fluid Disks

The stability of a rotating fluid disk to the formation of spiral arms is studied in the tightwinding approximation in the linear regime. The dispersion relation for spirals that was derived by Bertin et al. is shown to contain a new, acoustic instability beyond the Lindblad resonances that depends only on pressure and rotation. In this regime, pressure and gravity exchange roles as drivers and inhibitors of spiral wave structures. Other instabilities that are enhanced by pressure are also found in the general dispersion relation by including higher order terms in the small parameter 1/kr for wavenumber k and radius r. These instabilities are present even for large values of Toomre's parameter Q. Unstable growth rates are determined in four cases: a self-gravitating disk with a flat rotation curve, a self-gravitating disk with solid body rotation, a non-self-gravitating disk with solid body rotation, and a non-self-gravitating disk with Keplerian rotation. The most important application appears to be as a source of spiral structure, possibly leading to accretion in non-self-gravitating disks, such as some galactic nuclear disks, disks around black holes, and proto-planetary disks. All of these examples have short orbital times so the unstable growth time can be small.Comment: 30 pages, 5 figures, scheduled for ApJ 520, August 1, 199

Molecular Hydrogen Kinematics in Cepheus A

We present the radial velocity structure of the molecular hydrogen outflows associated to the star forming region Cepheus A. This structure is derived from doppler shift of the H_2 v=1-0 S(1) emission line obtained by Fabry-Perot spectroscopy. The East and West regions of emission, called Cep A(E) and Cep A(W), show radial velocities in the range -20 to 0 km/s with respect to the molecular cloud. Cep A(W) shows an increasing velocity with position offset from the core indicating the existence of a possible accelarating machanism. Cep A(E) has an almost constant mean radial velocity of -18 km/s along the region although with a large dispersion in velocity, indicating the possibility of a turbulent outflow. A detailed analysis of the Cep A(E) region shows evidence for the presence of a Mach disk on that outflow. Also, we argue that the presence of a velocity gradient in Cep A(W) is indicative of a C-shock in this region. Following Riera et al. (2003), we analyzed the data using wavelet analysis to study the line width and the central radial velocity distributions. We found that both outflows have complex spatial and velocity structures characteristic of a turbulent flow.Comment: 24 pages, 15 figure