No Arabic abstract
We study the conditions for disk galaxies to produce superbubbles that can break out of the disk and produce a galactic wind. We argue that the threshold surface density of supernovae rate for seeding a wind depends on the ability of superbubble energetics to compensate for radiative cooling. We first adapt Kompaneets formalism for expanding bubbles in a stratified medium to the case of continuous energy injection and include the effects of radiative cooling in the shell. With the help of hydrodynamic simulations, we then study the evolution of superbubbles evolving in stratified disks with typical disk parameters. We identify two crucial energy injection rates that differ in their effects, the corresponding breakout ranging from being gentle to a vigorous one. (a) Superbubbles that break out of the disk with a Mach number of order 2-3 correspond to an energy injection rate of order 10^{-4} erg cm^{-2} s^{-1}, which is relevant for disk galaxies with synchrotron emitting gas in the extra-planar regions. (b) A larger energy injection threshold, of order 10^{-3} erg cm^{-2} s^{-1}, or equivalently, a star formation surface density of sim 0.1 solar mass yr^{-1} kpc^{-2}, corresponds to superbubbles with a Mach number sim 5-10. While the milder superbubbles can be produced by large OB associations, the latter kind requires super-starclusters. These derived conditions compare well with observations of disk galaxies with winds and the existence of multiphase halo gas. Furthermore, we find that contrary to the general belief that superbubbles fragment through Rayleigh-Taylor (RT) instability when they reach a vertical height of order the scale height, the superbubbles are first affected by thermal instability for typical disk parameters and that RT instability takes over when the shells reach a distance of approximately twice the scale height.
We present results from high-resolution hydrodynamic simulations of isolated SMC- and Milky Way-sized galaxies that include a model for feedback from galactic cosmic rays (CRs). We find that CRs are naturally able to drive winds with mass loading factors of up to ~10 in dwarf systems. The scaling of the mass loading factor with circular velocity between the two simulated systems is consistent with propto v_c^{1-2} required to reproduce the faint end of the galaxy luminosity function. In addition, simulations with CR feedback reproduce both the normalization and the slope of the observed trend of wind velocity with galaxy circular velocity. We find that winds in simulations with CR feedback exhibit qualitatively different properties compared to SN driven winds, where most of the acceleration happens violently in situ near star forming sites. In contrast, the CR-driven winds are accelerated gently by the large-scale pressure gradient established by CRs diffusing from the star-forming galaxy disk out into the halo. The CR-driven winds also exhibit much cooler temperatures and, in the SMC-sized system, warm (T~10^4 K) gas dominates the outflow. The prevalence of warm gas in such outflows may provide a clue as to the origin of ubiquitous warm gas in the gaseous halos of galaxies detected via absorption lines in quasar spectra.
Accretion disk winds are thought to produce many of the characteristic features seen in the spectra of active galactic nuclei (AGN) and quasi-stellar objects (QSOs). These outflows also represent a natural form of feedback between the central supermassive black hole and its host galaxy. The mechanism for driving this mass loss remains unknown, although radiation pressure mediated by spectral lines is a leading candidate. Here, we calculate the ionization state of, and emergent spectra for, the hydrodynamic simulation of a line-driven disk wind previously presented by Proga & Kallman (2004). To achieve this, we carry out a comprehensive Monte Carlo simulation of the radiative transfer through, and energy exchange within, the predicted outflow. We find that the wind is much more ionized than originally estimated. This is in part because it is much more difficult to shield any wind regions effectively when the outflow itself is allowed to reprocess and redirect ionizing photons. As a result, the calculated spectrum that would be observed from this particular outflow solution would not contain the ultraviolet spectral lines that are observed in many AGN/QSOs. Furthermore, the wind is so highly ionized that line-driving would not actually be efficient. This does not necessarily mean that line-driven winds are not viable. However, our work does illustrate that in order to arrive at a self-consistent model of line-driven disk winds in AGN/QSO, it will be critical to include a more detailed treatment of radiative transfer and ionization in the next generation of hydrodynamic simulations.
Feedback from supernovae is an essential aspect of galaxy formation. In order to improve subgrid models of feedback we perform a series of numerical experiments to investigate how supernova explosions power galactic winds. We use the Flash hydrodynamic code to model a simplified ISM, including gravity, hydrodynamics, radiative cooling above 10,000 K, and star formation that reproduces the Kennicutt-Schmidt relation. By simulating a small patch of the ISM in a tall box perpendicular to the disk, we obtain sub-parsec resolution allowing us to resolve individual supernova events and we investigate how the wind properties depend on those of the ISM and the galaxy. We find that outflows are more efficient in disks with lower surface densities or gas fractions. A simple model in which the warm cloudy medium is the barrier that limits the expansion of blast waves reproduces the scaling of outflow properties with disk parameters at high star formation rates. The scaling we find sets the investigation of galaxy winds on a new footing, providing a physically motivated sub-grid description of winds that can be implemented in cosmological hydrodynamic simulations and phenomenological models. [Abridged]
Galactic winds are a common phenomenon in starburst galaxies in the local universe as well as at higher redshifts. Their sources are superbubbles driven by sequential supernova explosions in star forming regions, which carve out large holes in the interstellar medium and eject hot, metal enriched gas into the halo and to the galactic neighborhood. We investigate the evolution of superbubbles in exponentially stratified disks. We present advanced analytical models for the expansion of such bubbles and calculate their evolution in space and time. With these models one can derive the energy input that is needed for blow-out of superbubbles into the halo and derive the break-up of the shell, since Rayleigh-Taylor instabilities develop soon after a bubble starts to accelerate into the halo. The approximation of Kompaneets is modified in order to calculate velocity and acceleration of a bubble analytically. Our new model differs from earlier ones, because it presents for the first time an analytical calculation for the expansion of superbubbles in an exponential density distribution driven by a time-dependent energy input rate. The time-sequence of supernova explosions of OB-stars is modeled using their main sequence lifetime and an initial mass function. We calculate the morphology and kinematics of superbubbles powered by three different kinds of energy input and we derive the energy input required for blow-out as a function of the density and the scale height of the ambient interstellar medium. The Rayleigh-Taylor instability timescale in the shell is calculated in order to estimate when the shell starts to fragment and finally breaks up. Analytical models are a very efficient tool for comparison to observations, like e.g. the Local Bubble and the W4 bubble discussed in this paper, and also give insight into the dynamics of superbubble evolution.
The pattern speeds of spiral galaxies are closely related to the flow of material in their disks. Flows that follow the `precessing ellipses paradigm (see e.g., Kalnajs 1973) are likely associated with slowly rotating spirals, which have corotation beyond their end. Such a flow can be secured by material trapped around stable, elliptical, x_1 periodic orbits precessing as their Jacobi constant varies. Contrarily, if part of the spiral arms is located at a corotation region then the spiral structure has to `survive in chaotic regions. Barred-spiral systems with a single pattern speed and a bar ending before, but close to, corotation are candidates for having spirals supported by stars in chaotic motion. In this work we review the flows we have found in response models for various types of spiral potentials and indicate the cases, where order or chaos shapes the observed morphologies.