No Arabic abstract
Winds from massive stars have velocities of 1000 km/s or more, and produce hot, high pressure gas when they shock. We develop a theory for the evolution of bubbles driven by the collective winds from star clusters early in their lifetimes, which involves interaction with the turbulent, dense interstellar medium of the surrounding natal molecular cloud. A key feature is the fractal nature of the hot bubbles surface. The large area of this interface with surrounding denser gas strongly enhances energy losses from the hot interior, enabled by turbulent mixing and subsequent cooling at temperatures T = 10^4-10^5 K where radiation is maximally efficient. Due to the extreme cooling, the bubble radius scales differently (R ~ t^1/2) from the classical Weaver77 solution, and has expansion velocity and momentum lower by factors of 10-10^2 at given R, with pressure lower by factors of 10^2 - 10^3. Our theory explains the weak X-ray emission and low shell expansion velocities of observed sources. We discuss further implications of our theory for observations of the hot bubbles and cooled expanding shells created by stellar winds, and for predictions of feedback-regulated star formation in a range of environments. In a companion paper, we validate our theory with a suite of hydrodynamic simulations.
In a companion paper, we develop a theory for the evolution of stellar wind driven bubbles in dense, turbulent clouds. This theory proposes that turbulent mixing at a fractal bubble-shell interface leads to highly efficient cooling, in which the vast majority of the input wind energy is radiated away. This energy loss renders the majority of the bubble evolution momentum-driven rather than energy-driven, with expansion velocities and pressures orders of magnitude lower than in the classical Weaver77 solution. In this paper, we validate our theory with three-dimensional, hydrodynamic simulations. We show that extreme cooling is not only possible, but is generic to star formation in turbulent clouds over more than three orders of magnitude in density. We quantify the few free parameters in our theory, and show that the momentum exceeds the wind input rate by only a factor ~ 1.2-4. We verify that the bubble/cloud interface is a fractal with dimension ~ 2.5-2.7. The measured turbulent amplitude (v_t ~ 200-400 km/s) in the hot gas near the interface is shown to be consistent with theoretical requirements for turbulent diffusion to efficiently mix and radiate away most of the wind energy. The fraction of energy remaining after cooling is only 1-Theta ~ 0.1-0.01, decreasing with time, explaining observations that indicate low hot-gas content and weak dynamical effects of stellar winds.
To date, most numerical simulations of molecular clouds, and star formation within them, assume a uniform density sphere or box with an imposed turbulent velocity field. In this work, we select molecular clouds from galactic scale simulations as initial conditions, increase their resolution, and re-simulate them using the SPH code Gadget2. Our approach provides clouds with morphologies, internal structures, and kinematics that constitute more consistent and realistic initial conditions for simulations of star formation. We perform comparisons between molecular clouds derived from a galactic simulation, and spheres of turbulent gas of similar dimensions, mass and velocity dispersion. We focus on properties of the clouds such as their density, velocity structure and star formation rate. We find that the inherited velocity structure of the galactic clouds has a significant impact on the star formation rate and evolution of the cloud. Our results indicate that, although we can follow the time evolution of star formation in any simulated cloud, capturing the entire history is difficult as we ignore any star formation that might have occurred before initialisation. Overall, the turbulent spheres do not match the complexity of the galactic clouds.
We use 3D-PDR, a three-dimensional astrochemistry code for modeling photodissociation regions (PDRs), to post-process hydrodynamic simulations of turbulent, star-forming clouds. We focus on the transition from atomic to molecular gas, with specific attention to the formation and distribution of H, C+, C, H2 and CO. First, we demonstrate that the details of the cloud chemistry and our conclusions are insensitive to the simulation spatial resolution, to the resolution at the cloud edge, and to the ray angular resolution. We then investigate the effect of geometry and simulation parameters on chemical abundances and find weak dependence on cloud morphology as dictated by gravity and turbulent Mach number. For a uniform external radiation field, we find similar distributions to those derived using a one-dimensional PDR code. However, we demonstrate that a three-dimensional treatment is necessary for a spatially varying external field, and we caution against using one-dimensional treatments for non-symmetric problems. We compare our results with the work of Glover et al. (2010), who self-consistently followed the time evolution of molecule formation in hydrodynamic simulations using a reduced chemical network. In general, we find good agreement with this in situ approach for C and CO abundances. However, the temperature and H2 abundances are discrepant in the boundary regions (Av < 5), which is due to the different number of rays used by the two approaches.
We present analytic theory for the total column density of singly ionized carbon (C+) in the optically thick photon dominated regions (PDRs) of far-UV irradiated (star-forming) molecular clouds. We derive a simple formula for the C+ column as a function of the cloud (hydrogen) density, the far-UV field intensity, and metallicity, encompassing the wide range of galaxy conditions. When assuming the typical relation between UV and density in the cold neutral medium, the C+ column becomes a function of the metallicity alone. We verify our analysis with detailed numerical PDR models. For optically thick gas, most of the C+ column is mixed with hydrogen that is primarily molecular (H2), and this C+/H2 gas layer accounts for almost all of the `CO-dark molecular gas in PDRs. The C+/H2 column density is limited by dust shielding and is inversely proportional to the metallicity down to ~0.1 solar. At lower metallicities, H2 line blocking dominates and the C+/H2 column saturates. Applying our theory to CO surveys in low redshift spirals we estimate the fraction of C+/H2 gas out of the total molecular gas to be typically ~0.4. At redshifts 1<z<3 in massive disc galaxies the C+/H2 gas represents a very small fraction of the total molecular gas (<0.16). This small fraction at high redshifts is due to the high gas surface densities when compared to local galaxies.
A model of magnetic field structure is presented to help test the prevalence of flux freezing in star-forming clouds of various shapes, orientations, and degrees of central concentration, and to estimate their magnetic field strength. The model is based on weak-field flux freezing in centrally condensed Plummer spheres and spheroids of oblate and prolate shape. For a spheroid of given density contrast, aspect ratio, and inclination, the model estimates the local field strength and direction, and the global field pattern of hourglass shape. Comparisons with a polarization simulation indicate typical angle agreement within 1 - 10 degrees. Scalable analytic expressions are given to match observed polarization patterns, and to provide inputs to radiative transfer codes for more accurate predictions. The model may apply to polarization observations of dense cores, elongated filamentary clouds, and magnetized circumstellar disks.