Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userModeling UV Radiation Feedback from Massive Stars: III. Escape of Radiation from Star-forming Giant Molecular Clouds
102
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Using a suite of radiation hydrodynamic simulations of star cluster formation in turbulent clouds, we study the escape fraction of ionizing (Lyman continuum) and non-ionizing (FUV) radiation for a wide range of cloud masses and sizes. The escape fraction increases as H II regions evolve and reaches unity within a few dynamical times. The cumulative escape fraction before the onset of the first supernova explosion is in the range 0.05-0.58; this is lower for higher initial cloud surface density, and higher for less massive and more compact clouds due to rapid destruction. Once H II regions break out of their local environment, both ionizing and non-ionizing photons escape from clouds through fully ionized, low-density sightlines. Consequently, dust becomes the dominant absorber of ionizing radiation at late times and the escape fraction of non-ionizing radiation is only slightly larger than that of ionizing radiation. The escape fraction is determined primarily by the mean $langle taurangle$ and width $sigma$ of the optical-depth distribution in the large-scale cloud, increasing for smaller $langle taurangle$ and/or larger $sigma$. The escape fraction exceeds (sometimes by three orders of magnitude) the naive estimate $e^{-langle taurangle}$ due to non-zero $sigma$ induced by turbulence. We present two simple methods to estimate, within $sim20%$, the escape fraction of non-ionizing radiation using the observed dust optical depth in clouds projected on the plane of sky. We discuss implications of our results for observations, including inference of star formation rates in individual molecular clouds, and accounting for diffuse ionized gas on galactic scales.
rate research
Read More
UV radiation feedback from young massive stars plays a key role in the evolution of giant molecular clouds (GMCs) by photoevaporating and ejecting the surrounding gas. We conduct a suite of radiation hydrodynamic simulations of star cluster formation in marginally-bound, turbulent GMCs, focusing on the effects of photoionization and radiation pressure on regulating the net star formation efficiency (SFE) and cloud lifetime. We find that the net SFE depends primarily on the initial gas surface density, $Sigma_0$, such that the SFE increases from 4% to 51% as $Sigma_0$ increases from $13,M_{odot},{rm pc}^{-2}$ to $1300,M_{odot},{rm pc}^{-2}$. Cloud destruction occurs within $2$-$10,{rm Myr}$ after the onset of radiation feedback, or within $0.6$-$4.1$ freefall times (increasing with $Sigma_0$). Photoevaporation dominates the mass loss in massive, low surface-density clouds, but because most photons are absorbed in an ionization-bounded Str{o}mgren volume the photoevaporated gas fraction is proportional to the square root of the SFE. The measured momentum injection due to thermal and radiation pressure forces is proportional to $Sigma_0^{-0.74}$, and the ejection of neutrals substantially contributes to the disruption of low-mass and/or high-surface density clouds. We present semi-analytic models for cloud dispersal mediated by photoevaporation and by dynamical mass ejection, and show that the predicted net SFE and mass loss efficiencies are consistent with the results of our numerical simulations.
We present an implementation of an adaptive ray tracing (ART) module in the Athena hydrodynamics code that accurately and efficiently handles the radiative transfer involving multiple point sources on a three-dimensional Cartesian grid. We adopt a recently proposed parallel algorithm that uses non-blocking, asynchronous MPI communications to accelerate transport of rays across the computational domain. We validate our implementation through several standard test problems including the propagation of radiation in vacuum and the expansions of various types of HII regions. Additionally, scaling tests show that the cost of a full ray trace per source remains comparable to that of the hydrodynamics update on up to $sim 10^3$ processors. To demonstrate application of our ART implementation, we perform a simulation of star cluster formation in a marginally bound, turbulent cloud, finding that its star formation efficiency is $12%$ when both radiation pressure forces and photoionization by UV radiation are treated. We directly compare the radiation forces computed from the ART scheme with that from the M1 closure relation. Although the ART and M1 schemes yield similar results on large scales, the latter is unable to resolve the radiation field accurately near individual point sources.
Molecular clouds are supported by turbulence and magnetic fields, but quantifying their influence on cloud lifecycle and star formation efficiency (SFE) remains an open question. We perform radiation MHD simulations of star-forming giant molecular clouds (GMCs) with UV radiation feedback, in which the propagation of UV radiation via ray-tracing is coupled to hydrogen photochemistry. We consider 10 GMC models that vary in either initial virial parameter ($1lealpha_{v,0}le 5$) or dimensionless mass-to-magnetic flux ratio (0.5-8 and $infty$); the initial mass $10^5M_{odot}$ and radius 20pc are fixed. Each model is run with five different initial turbulence realizations. In most models, the duration of star formation and the timescale for molecular gas removal (primarily by photoevaporation) are 4-8Myr. Both the final SFE ($epsilon_*$) and time-averaged SFE per freefall time ($epsilon_{ff}$) are reduced by strong turbulence and magnetic fields. The median $epsilon_*$ ranges between 2.1% and 9.5%. The median $epsilon_{ff}$ ranges between 1.0% and 8.0% and anticorrelates with $alpha_{v,0}$, in qualitative agreement with previous analytic theory and simulations. However, the time-dependent $alpha_{v}(t)$ and $epsilon_{ff,obs}(t)$ based on instantaneous gas properties and cluster luminosity are positively correlated due to rapid evolution, making observational validation of star formation theory difficult. Our median $epsilon_{ff,obs}(t)approx$ 2% is similar to observed values. We show that the traditional virial parameter estimates the true gravitational boundedness within a factor of 2 on average, but neglect of magnetic support and velocity anisotropy can sometimes produce large departures. Magnetically subcritical GMCs are unlikely to represent sites of massive star formation given their unrealistic columnar outflows, prolonged lifetime, and low escape fraction of radiation.
We extend our previous SPH parameter study of the effects of photoionization from O-stars on star-forming clouds to include initially unbound clouds. We generate a set of model clouds in the mass range $10^{4}-10^{6}$M$_{odot}$ with initial virial ratios $E_{rm kin}/E_{rm pot}$=2.3, allow them to form stars, and study the impact of the photoionizing radiation produced by the massive stars. We find that, on the 3Myr timescale before supernovae are expected to begin detonating, the fractions of mass expelled by ionizing feedback is a very strong function of the cloud escape velocities. High-mass clouds are largely unaffected dynamically, while lower-mass clouds have large fractions of their gas reserves expelled on this timescale. However, the fractions of stellar mass unbound are modest and significant portions of the unbound stars are so only because the clouds themselves are initially partially unbound. We find that ionization is much more able to create well-cleared bubbles in the unbound clouds, owing to their intrinsic expansion, but that the presence of such bubbles does not necessarily indicate that a given cloud has been strongly influenced by feedback. We also find, in common with the bound clouds from our earlier work, that many of the systems simulated here are highly porous to photons and supernova ejecta, and that most of them will likely survive their first supernova explosions.
Using results from high-resolution galaxy formation simulations in a standard Lambda-CDM cosmology and a fully conservative multi-resolution radiative transfer code around point sources, we compute the energy-dependent escape fraction of ionizing photons from a large number of star forming regions in two galaxies at five different redshifts from z=3.8 to 2.39. All escape fractions show a monotonic decline with time, from (at the Lyman-limit) ~6-10% at z=3.6 to ~1-2% at z=2.39, due to higher gas clumping at lower redshifts. It appears that increased feedback can lead to higher f_esc at z>3.4 via evacuation of gas from the vicinity of star forming regions and to lower f_esc at z<2.39 through accumulation of swept-up shells in denser environments. Our results agree well with the observational findings of citet{inoue..06} on redshift evolution of f_esc in the redshift interval z=2-3.6.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
