The most accurate measurements of magnetic fields in star-forming gas are based on the Zeeman observations analyzed by Crutcher et al. (2010). We show that their finding that the 3D magnetic field scales approximately as density$^{0.65}$ can also be
obtained from analysis of the observed line-of-sight fields. We present two large-scale AMR MHD simulations of several thousand $M_odot$ of turbulent, isothermal, self-gravitating gas, one with a strong initial magnetic field (Alfven Mach number $M_{A,0}= 1$) and one with a weak initial field ($M_{A,0}=10$). We construct samples of the 100 most massive clumps in each simulation and show that they exhibit a power-law relation between field strength and density in excellent agreement with the observed one. Our results imply that the average field in molecular clumps in the interstellar medium is $<B_{tot}> sim 42 n_{H,4}^{0.65} mu$G. Furthermore, the median value of the ratio of the line-of-sight field to density$^{0.65}$ in the simulations is within a factor of about (1.3, 1.7) of the observed value for the strong and weak field cases, respectively. The median value of the mass-to-flux ratio, normalized to the critical value, is 70% of the line-of-sight value. This is larger than the 50% usually cited for spherical clouds because the actual mass-to-flux ratio depends on the volume-weighted field, whereas the observed one depends on the mass-weighted field. Our results indicate that the typical molecular clump in the ISM is significantly supercritical (~ factor of 3). The results of our strong-field model are in very good quantitative agreement with the observations of Li et al. (2009), which show a strong correlation in field orientation between small and large scales. Because there is a negligible correlation in the weak-field model, we conclude that molecular clouds form from strongly magnetized (although magnetically supercritical) gas.
The turbulent environment from which stars form may lead to misalignment between the stellar spin and the remnant protoplanetary disk. By using hydrodynamic and magnetohydrodynamic simulations, we demonstrate that a wide range of stellar obliquities
may be produced as a by-product of forming a star within a turbulent environment. We present a simple semi-analytic model that reveals this connection between the turbulent motions and the orientation of a star and its disk. Our results are consistent with the observed obliquity distribution of hot Jupiters. Migration of misaligned hot Jupiters may, therefore, be due to tidal dissipation in the disk, rather than tidal dissipation of the star-planet interaction.
The temperature of the gas in molecular clouds is a key determinant of the characteristic mass of star formation. Ambipolar diffusion (AD) is considered one of the most important heating mechanisms in weakly ionized molecular clouds. In this work, we
study the AD heating rate using 2-fluid turbulence simulations and compare it with the overall heating rate due to turbulent dissipation. We find that for observed molecular clouds, which typically have Alfven Mach numbers of ~1 (Crutcher 1999) and AD Reynolds numbers of ~20 (McKee et al. 2010), about 70% of the total turbulent dissipation is in the form of AD heating. AD has an important effect on the length scale where energy is dissipated: when AD heating is strong, most of the energy in the cascade is removed by ion-neutral drift, with a comparatively small amount of energy making it down to small scales. We derive a relation for the AD heating rate that describes the results of our simulations to within a factor of two. Turbulent dissipation, including AD heating, is generally less important that cosmic-ray heating in molecular clouds, although there is substantial scatter in both.
Protostellar feedback, both radiation and bipolar outflows, dramatically affects the fragmentation and mass accretion from star-forming cores. We use ORION, an adaptive mesh refinement gravito-radiation-hydrodynamics code, to simulate the formation o
f a cluster of low-mass stars, including both radiative transfer and protostellar outflows. We ran four simulations to isolate the individual effects of radiation feedback and outflow feedback as well as the combination of the two. Outflows reduce protostellar masses and accretion rates each by a factor of three and therefore reduce protostellar luminosities by an order of magnitude. Thus, while radiation feedback suppresses fragmentation, outflows render protostellar radiation largely irrelevant for low-mass star formation above a mass scale of 0.05 M_sun. We find initial fragmentation of our cloud at half the global Jeans length, ~ 0.1 pc. With insufficient protostellar radiation to stop it, these 0.1 pc cores fragment repeatedly, forming typically 10 stars each. The accretion rate in these stars scales with mass as predicted from core accretion models that include both thermal and turbulent motions. We find that protostellar outflows do not significantly affect the overall cloud dynamics, in the absence of magnetic fields, due to their small opening angles and poor coupling to the dense gas. The outflows reduce the mass from the cores by 2/3, giving a core to star efficiency ~ 1/3. The simulation with radiation and outflows reproduces the observed protostellar luminosity function. All of the simulations can reproduce observed core mass functions, though they are sensitive to telescope resolution. The simulation with both radiation and outflows reproduces the galactic IMF and the two-point correlation function of the cores observed in rho Oph.
The protostellar luminosity function (PLF) is the present-day luminosity function of the protostars in a region of star formation. It is determined using the protostellar mass function (PMF) in combination with a stellar evolutionary model that provi
des the luminosity as a function of instantaneous and final stellar mass. As in McKee & Offner (2010), we consider three main accretion models: the Isothermal Sphere model, the Turbulent Core model, and an approximation of the Competitive Accretion model. We also consider the effect of an accretion rate that tapers off linearly in time and an accelerating star formation rate. For each model, we characterize the luminosity distribution using the mean, median, maximum, ratio of the median to the mean, standard deviation of the logarithm of the luminosity, and the fraction of very low luminosity objects. We compare the models with bolometric luminosities observed in local star forming regions and find that models with an approximately constant accretion time, such as the Turbulent Core and Competitive Accretion models, appear to agree better with observation than those with a constant accretion rate, such as the Isothermal Sphere model. We show that observations of the mean protostellar luminosity in these nearby regions of low-mass star formation suggest a mean star formation time of 0.3$pm$0.1 Myr. Such a timescale, together with some accretion that occurs non-radiatively and some that occurs in high-accretion, episodic bursts, resolves the classical luminosity problem in low-mass star formation, in which observed protostellar luminosities are significantly less than predicted. An accelerating star formation rate is one possible way of reconciling the observed star formation time and mean luminosity.
We report new mid-infrared observations of the remarkable object IRAS 15099-5856 using the space telescopes AKARI and Spitzer, which demonstrate the presence of prominent crystalline silicate emission in this bright source. IRAS 15099-5856 has a comp
lex morphology with a bright central compact source (IRS1) surrounded by knots, spurs, and several extended (~4) arc-like filaments. The source is seen only at >= 10 um. The Spitzer MIR spectrum of IRS1 shows prominent emission features from Mg-rich crystalline silicates, strong [Ne II] 12.81 um and several other faint ionic lines. We model the MIR spectrum as thermal emission from dust and compare with the Herbig Be star HD 100546 and the luminous blue variable R71, which show very similar MIR spectra. Molecular line observations reveal two molecular clouds around the source, but no associated dense molecular cores. We suggest that IRS1 is heated by UV radiation from the adjacent O star Muzzio 10 and that its crystalline silicates most likely originated in a mass outflow from the progenitor of the supernova remnant (SNR) MSH 15-52. IRS1, which is embedded in the SNR, could have been shielded from the SN blast wave if the progenitor was in a close binary system with Muzzio 10. If MSH15-52 is a remnant of Type Ib/c supernova (SN Ib/c), as has been previously proposed, this would confirm the binary model for SN Ib/c. IRS1 and the associated structures may be the relics of massive star death, as shaped by the supernova explosion, the pulsar wind and the intense ionizing radiation of the embedded O star.
We derive a semi-empirical galactic initial mass function (IMF) from observational constraints. We assume that the star formation rate in a galaxy can be expressed as the product of the IMF, $psi (m)$, which is a smooth function of mass $m$ (in units
of msun), and a time- and space-dependent total rate of star formation per unit area of galactic disk. The mass dependence of the proposed IMF is determined by five parameters: the low-mass slope $gamma$, the high-mass slope $-Gamma$, the characteristic mass $m_{ch}$ (which is close to the mass $m_{rm peak}$ at which the IMF turns over), and the lower and upper limits on the mass, $m_l$ (taken to be 0.004) and $m_u$ (taken to be 120). The star formation rate in terms of number of stars per unit area of galactic disk per unit logarithmic mass interval, is proportional to $m^{-Gamma} left{1-expleft[{-(m/m_{ch})^{gamma +Gamma}}right]right}$, where $cal N_*$ is the number of stars, $m_l<m<m_u$ is the range of stellar masses. The values of $gamma$ and $emch$ are derived from two integral constraints: i) the ratio of the number density of stars in the range $m=0.1-0.6$ to that in the range $m=0.6-0.8$ as inferred from the mass distribution of field stars in the local neighborhood, and ii) the ratio of the number of stars in the range $m=0.08 - 1$ to the number of brown dwarfs in the range $m=0.03-0.08$ in young clusters. The IMF satisfying the above constraints is characterized by the parameters $gamma=0.51$ and $emch=0.35$ (which corresponds to $m_{rm peak}=0.27$). This IMF agrees quite well with the Chabrier (2005) IMF for the entire mass range over which we have compared with data, but predicts significantly more stars with masses $< 0.03, M_odot$; we also compare with other IMFs in current use.
Most numerical investigations on the role of magnetic fields in turbulent molecular clouds (MCs) are based on ideal magneto-hydrodynamics (MHD). However, MCs are weakly ionized, so that the time scale required for the magnetic field to diffuse throug
h the neutral component of the plasma by ambipolar diffusion (AD) can be comparable to the dynamical time scale. We have performed a series of 256^3 and 512^3 simulations on supersonic but sub-Alfvenic turbulent systems with AD using the Heavy-Ion Approximation developed in Li, McKee, & Klein (2006). Our calculations are based on the assumption that the number of ions is conserved, but we show that these results approximately apply to the case of time-dependent ionization in molecular clouds as well. Convergence studies allow us to determine the optimal value of the ionization mass fraction when using the heavy-ion approximation for low Mach number, sub-Alfvenic turbulent systems. We find that ambipolar diffusion steepens the velocity and magnetic power spectra compared to the ideal MHD case. Changes in the density PDF, total magnetic energy, and ionization fraction are determined as a function of the AD Reynolds number. The power spectra for the neutral gas properties of a strongly magnetized medium with a low AD Reynolds number are similar to those for a weakly magnetized medium; in particular, the power spectrum of the neutral velocity is close to that for Burgers turbulence.
Massive stars are very rare, but their extreme luminosities make them both the only type of young star we can observe in distant galaxies and the dominant energy sources in the universe today. They form rarely because efficient radiative cooling keep
s most star-forming gas clouds close to isothermal as they collapse, and this favors fragmentation into stars <~1 Msun. Heating of a cloud by accreting low-mass stars within it can prevent fragmentation and allow formation of massive stars, but what properties a cloud must have to form massive stars, and thus where massive stars form in a galaxy, has not yet been determined. Here we show that only clouds with column densities >~ 1 g cm^-2 can avoid fragmentation and form massive stars. This threshold, and the environmental variation of the stellar initial mass function (IMF) that it implies, naturally explain the characteristic column densities of massive star clusters and the difference between the radial profiles of Halpha and UV emission in galactic disks. The existence of a threshold also implies that there should be detectable variations in the IMF with environment within the Galaxy and in the characteristic column densities of massive star clusters between galaxies, and that star formation rates in some galactic environments may have been systematically underestimated.
We report the discovery of a star-forming loop around the young, Crab-like supernova remnant (SNR) G54.1+0.3 using the AKARI infrared satellite. The loop consists of at least eleven young stellar objects (YSOs) embedded in a ring-like diffuse emissio
n of radius ~1. The YSOs are bright in the mid-infrared and are also visible in the Spitzer Space Telescope Galactic plane survey images. Their Spitzer colors are similar to those of class II YSOs in [3.6]-[5.8] but significantly redder in [8]-[24], i.e., 0<[3.6]-[5.8]<1.2 and 5<[8]-[24]<9. Most of them have near-infrared counterparts in the 2MASS JHKs images, and some of them have an optical counterpart too. Their JHKs colors and magnitudes indicate that the YSOs are massive (<= 10 Msun) pre-main-sequence stars at the same distance to the SNR, i.e., 8 kpc, which supports the association of the star-forming loop with the SNR. The dereddened spectral energy distributions are similar to eraly Herbig Be stars, which are early B-type pre-main-sequence stars with inner disks that have been destroyed. The confinement to a loop structure indicates that the YSOs are young, i.e., <= 2 Myr. We propose that their formation is triggered by the progenitor star of G54.1+0.3, which has a mass of <= 15 Msun. The triggering must have occurred near the end of the progenitors life, possibly after it had evolved off the main sequence.
