No Arabic abstract
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.
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.
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 a highly-parallel multi-frequency hybrid radiation hydrodynamics algorithm that combines a spatially-adaptive long characteristics method for the radiation field from point sources with a moment method that handles the diffuse radiation field produced by a volume-filling fluid. Our Hybrid Adaptive Ray-Moment Method (HARM$^2$) operates on patch-based adaptive grids, is compatible with asynchronous time stepping, and works with any moment method. In comparison to previous long characteristics methods, we have greatly improved the parallel performance of the adaptive long-characteristics method by developing a new completely asynchronous and non-blocking communication algorithm. As a result of this improvement, our implementation achieves near-perfect scaling up to $mathcal{O}(10^3)$ processors on distributed memory machines. We present a series of tests to demonstrate the accuracy and performance of the method.
We use 3D hydrodynamical models to investigate the effects of massive star feedback from winds and supernovae on inhomogeneous molecular material left over from the formation of a massive stellar cluster. We simulate the interaction of the mechanical energy input from a cluster with 3 O-stars into a giant molecular cloud (GMC) clump containing 3240 solar masses of molecular material within a 4 pc radius. The cluster wind blows out of the molecular clump along low-density channels, into which denser clump material is entrained. We find that the densest molecular regions are surprisingly resistant to ablation by the cluster wind, in part due to shielding by other dense regions closer to the cluster. Nonetheless, molecular material is gradually removed by the cluster wind during which mass-loading factors in excess of several 100 are obtained. Because the clump is very porous, 60-75 per cent of the injected wind energy escapes the simulation domain, with the difference being radiated. After 4.4 Myr, the massive stars in our simulation begin to explode as supernovae. The highly structured environment into which the SN energy is released allows even weaker coupling to the remaining dense material and practically all of the SN energy reaches the wider environment. The molecular material is almost completely dispersed and destroyed after 6 Myr. The escape fraction of ionizing radiation is estimated to be about 50 per cent during the first 4 Myr of the clusters life. A similar model with a larger and more massive GMC clump reveals the same general picture, though more time is needed for it to be destroyed.
Massive stars are powerful sources of radiation, stellar winds, and supernova explosions. The radiative and mechanical energies injected by massive stars into the interstellar medium (ISM) profoundly alter the structure and evolution of the ISM, which subsequently influences the star formation and chemical evolution of the host galaxy. In this review, we will use the Large Magellanic Cloud (LMC) as a laboratory to showcase effects of energy feedback from massive young stellar objects (YSOs) and mature stars. We will also use the Carina Nebula in the Galaxy to illustrate a multi-wavelength study of feedback from massive star.