No Arabic abstract
It is shown that a number of key observations of the Galactic ISM can be understood, if it is treated as a highly compressible and turbulent medium energized predominantly by supernova explosions (and stellar winds). We have performed extensive numerical high resolution 3D hydrodynamical and magnetohydrodynamical simulations with adaptive mesh refinement over sufficiently long time scales to erase memory effects of the initial setup. Our results show, in good agrement with observations, that (i) volume filling factors of the hot medium are modest (typically below 20%), (ii) global pressure is far from uniform due to supersonic (and to some extent superalfvenic) turbulence, (iii) a significant fraction of the mass (about 60%) in the warm neutral medium is in the thermally unstable regime (500 K < T < 5000 K), (iv) the average number density of OVI in absorption is 1.81 10^{-8} cm^{-3}, in excellent agreement with Copernicus and FUSE data, and its distribution is rather clumpy, consistent with its measured dispersion with distance.
To study how supernova feedback structures the turbulent interstellar medium, we construct 3D models of vertically stratified gas stirred by discrete supernova explosions, including vertical gravitational field and parametrized heating and cooling. The models reproduce many observed characteristics of the Galaxy such as global circulation of gas (i.e., galactic fountain) and the existence of cold dense clouds in the galactic disk. Global quantities of the model such as warm and hot gas filling factors in the midplane, mass fraction of thermally unstable gas, and the averaged vertical density profile are compared directly with existing observations, and shown to be broadly consistent. We find that energy injection occurs over a broad range of scales. There is no single effective driving scale, unlike the usual assumption for idealized models of incompressible turbulence. However, >90% of the total kinetic energy is contained in wavelengths shortward of 200 pc. The shape of the kinetic energy spectrum differs substantially from that of the velocity power spectrum, which implies that the velocity structure varies with the gas density. Velocity structure functions demonstrate that the phenomenological theory proposed by Boldyrev is applicable to the medium. We show that it can be misleading to predict physical properties such as the stellar initial mass function based on numerical simulations that do not include self-gravity of the gas. Even if all the gas in turbulently Jeans unstable regions in our simulation is assumed to collapse and form stars in local freefall times, the resulting total collapse rate is significantly lower than the value consistent with the input supernova rate. Supernova-driven turbulence inhibits star formation globally rather than triggering it.
We investigate models for the photoionization of the widespread diffuse ionized gas in galaxies. In particular we address the long standing question of the penetration of Lyman continuum photons from sources close to the galactic midplane to large heights in the galactic halo. We find that recent hydrodynamical simulations of a supernova-driven interstellar medium have low density paths and voids that allow for ionizing photons from midplane OB stars to reach and ionize gas many kiloparsecs above the midplane. We find ionizing fluxes throughout our simulation grids are larger than predicted by one dimensional slab models, thus allowing for photoionization by O stars of low altitude neutral clouds in the Galaxy that are also detected in Halpha. In previous studies of such clouds the photoionization scenario had been rejected and the Halpha had been attributed to enhanced cosmic ray ionization or scattered light from midplane H II regions. We do find that the emission measure distributions in our simulations are wider than those derived from Halpha observations in the Milky Way. In addition, the horizontally averaged height dependence of the gas density in the hydrodynamical models is lower than inferred in the Galaxy. These discrepancies are likely due to the absence of magnetic fields in the hydrodynamic simulations and we discuss how magnetohydrodynamic effects may reconcile models and observations. Nevertheless, we anticipate that the inclusion of magnetic fields in the dynamical simulations will not alter our primary finding that midplane OB stars are capable of producing high altitude diffuse ionized gas in a realistic three-dimensional interstellar medium.
The mixed morphology class of supernova remnants has centrally peaked X-ray emission along with a shell-like morphology in radio emission. White & Long proposed that these remnants are evolving in a cloudy medium wherein the clouds are evaporated via thermal conduction once being overrun by the expanding shock. Their analytical model made detailed predictions regarding temperature, density and emission profiles as well as shock evolution. We present numerical hydrodynamical models in 2D and 3D including thermal conduction, testing the White & Long model and presenting results for the evolution and emission from remnants evolving in a cloudy medium. We find that, while certain general results of the White & Long model hold, such as the way the remnants expand and the flattening of the X-ray surface brightness distribution, in detail there are substantial differences. In particular we find that the X-ray luminosity is dominated by emission from shocked cloud gas early on, leading to a bright peak which then declines and flattens as evaporation becomes more important. In addition, the effects of thermal conduction on the intercloud gas, which is not included in the White & Long model, are important and lead to further flattening of the X-ray brightness profile as well as lower X-ray emission temperatures.
We report the highest-fidelity observations of the spiral galaxy M51 in CO emission, revealing the evolution of giant molecular clouds (GMCs) vis-a-vis the large-scale galactic structure and dynamics. The most massive GMCs (so-called GMAs) are first assembled and then broken up as the gas flow through the spiral arms. The GMAs and their H2 molecules are not fully dissociated into atomic gas as predicted in stellar feedback scenarios, but are fragmented into smaller GMCs upon leaving the spiral arms. The remnants of GMAs are detected as the chains of GMCs that emerge from the spiral arms into interarm regions. The kinematic shear within the spiral arms is sufficient to unbind the GMAs against self-gravity. We conclude that the evolution of GMCs is driven by large-scale galactic dynamics --their coagulation into GMAs is due to spiral arm streaming motions upon entering the arms, followed by fragmentation due to shear as they leave the arms on the downstream side. In M51, the majority of the gas remains molecular from arm entry through the inter-arm region and into the next spiral arm passage.
We study the evolution of dense clumps and provide argument that the existence of the clumps is not limited by the crossing time of the clump. We claim that the lifetimes of the clumps are determined by the turbulent motions on larger scale and predict the correlation of the clump lifetime and its column density. We use numerical simulations and successfully test this relation. In addition, we study the morphological asymmetry and the magnetization of the clumps as a function of their masses.