Using hydrodynamical simulations, we show for the first time that an episode of star formation in the center of the Milky Way, with a star-formation-rate (SFR) $sim 0.5$ M$_odot$ yr$^{-1}$ for $sim 30$ Myr, can produce bubbles that resemble the Fermi
Bubbles (FBs), when viewed from the solar position. The morphology, extent and multi-wavelength observations of FBs, especially X-rays, constrain various physical parameters such as SFR, age, and the circum-galactic medium (CGM) density. We show that the interaction of the CGM with the Galactic wind driven by a star formation in the central region can explain the observed surface brightness and morphological features of X-rays associated with the Fermi Bubbles. Furthermore, assuming that cosmic ray electrons are accelerated {it in situ} by shocks and/or turbulence, the brightness and morphology of gamma-ray emission and the microwave haze can be explained. The kinematics of the cold and warm clumps in our model also matches with recent observations of absorption lines through the bubbles.
It has been shown that the behaviour of primordial gas collapsing in a dark matter minihalo can depend on the adopted choice of 3-body H$_2$ formation rate. The uncertainties in this rate span two orders of magnitude in the current literature, and so
it remains a source of uncertainty in our knowledge of population III star formation. Here we investigate how the amount of fragmentation in primordial gas depends on the adopted 3-body rate. We present the results of calculations that follow the chemical and thermal evolution of primordial gas as it collapses in two dark matter minihalos. Our results on the effect of 3-body rate on the evolution until the first protostar forms agree well with previous studies. However, our modified version of GADGET-2 SPH also includes sink particles, which allows us to follow the initial evolution of the accretion disc that builds up on the centre of each halo, and capture the fragmentation in gas as well as its dependence on the adopted 3-body H$_2$ formation rate. We find that the fragmentation behaviour of the gas is only marginally effected by the choice of 3-body rate co-efficient, and that halo-to-halo differences are of equal importance in affecting the final mass distribution of stars.
In this paper we calculate the escape fraction ($f_{rm esc}$) of ionizing photons from starburst galaxies. Using 2-D axisymmetric hydrodynamic simulations, we study superbubbles created by overlapping supernovae in OB associations. We calculate the e
scape fraction of ionizing photons from the center of the disk along different angles through the superbubble and the gas disk. After convolving with the luminosity function of OB associations, we show that the ionizing photons escape within a cone of $sim 40 ^circ$, consistent with observations of nearby galaxies. The evolution of the escape fraction with time shows that it falls initially as cold gas is accumulated in a dense shell. After the shell crosses a few scale heights and fragments, the escape fraction through the polar regions rises again. The angle-averaged escape fraction cannot exceed $sim [1- cos (1 , {rm radian})] = 0.5$ from geometrical considerations (using the emission cone opening angle). We calculate the dependence of the time- and angle-averaged escape fraction on the mid-plane disk gas density (in the range $n_0=0.15-50$ cm $^{-3}$) and the disk scale height (between $z_0=10-600$ pc). We find that the escape fraction is related to the disk parameters (the mid-plane disk density and scale height) roughly so that $f_{rm esc}^alpha n_0^2 z_0^3$ (with $alphaapprox 2.2$) is a constant. For disks with a given WNM temperature, massive disks have lower escape fraction than low mass galaxies. For Milky Way ISM parameters, we find $f_{rm esc}sim 5%$, and it increases to $approx 10%$ for a galaxy ten times less massive. We discuss the possible effects of clumpiness of the ISM on the estimate of the escape fraction and the implications of our results for the reionization of the universe.
Using hydrodynamic simulations, we study the mass loss due to supernova-driven outflows from Milky Way type disk galaxies, paying particular attention to the effect of the extended hot halo gas. We find that the total mass loss at inner radii scales
roughly linearly with total mass of stars formed, and that the mass loading factor at the virial radius can be several times its value at inner radii because of the swept up hot halo gas. The temperature distribution of the outflowing material in the inner region ($sim $10 kpc) is bimodal in nature, peaking at $10^5$ K and $10^{6.5}$ K, responsible for optical and X-ray emission, respectively. The contribution of cold/warm gas with temperature $le 10^{5.5}$ K to the outflow rate within 10 kpc is $approx 0.3hbox{--}0.5$. The warm mass loading factor, $eta_{3e5}$ ($Tle 3 times 10^5$ K) is related to the mass loading factor at the virial radius ($eta_{v}$) as $eta_{v} approx 25, eta_{3e5}, left(mbox{SFR}/{rm M}_odot{rm yr}^{-1} right)^{-0.15}$ for a baryon fraction of 0.1 and a starburst period of 50 Myr. We also discuss the effect of multiple bursts that are separated by both short and long periods. The outflow speed at the virial radius is close to the sound speed in the hot halo, $lesssim 200$ km s$^{-1}$. We identify two `sequences of outflowing cold gas at small scales: a fast ($approx 500$ km~s$^{-1}$) sequence, driven by the unshocked free-wind; and a slow sequence ($approx pm 100$ km s$^{-1}$) at the conical interface of the superwind and the hot halo.
The statistics of black holes and their masses strongly suggests that their mass distribution has a cutoff towards lower masses near $3 times 10^{6}$ M$_{odot}$. This is consistent with a classical formation mechanism from the agglomeration of the fi
rst massive stars in the universe. However, when the masses of the stars approach $10^{6}$ M$_{odot}$, the stars become unstable and collapse, possibly forming the first generation of cosmological black holes. Here we speculate that the claimed detection of an isotropic radio background may constitute evidence of the formation of these first supermassive black holes, since their data are compatible in spectrum and intensity with synchrotron emission from the remnants. The model proposed fulfills all observational conditions for the background, in terms of single-source strength, number of sources, far-infrared and gamma-ray emission. The observed high energy neutrino flux is consistent with our calculations in flux and spectrum. The proposal described in this paper may also explain the early formation and growth of massive bulge-less disk galaxies as derived from the massive, gaseous shell formed during the explosion prior to the formation of a supermassive black hole.
(Abridged) Heating of the interstellar medium by multiple supernovae (SNe) explosions is at the heart of producing galaxy-scale outflows. We use hydrodynamical simulations to study the efficiency of multiple SNe in heating the interstellar medium (IS
M) and filling the volume with gas of high temperatures. We argue that it is important for SNe remnants to have a large filling factor {it and} a large heating efficiency. For this, they have to be clustered in space and time, and keep exploding until the hot gas percolates through the whole region, in order to compensate for the radiative loss. In the case of a limited number of SNe, we find that although the filling factor can be large, the heating efficiency declines after reaching a large value. In the case of a continuous series of SNe, the hot gas ($T ge 3 times 10^6$ K) can percolate through the whole region after the total volume filling factor reaches a threshold of $sim 0.3$. The efficiency of heating the gas to X-ray temperatures can be $ge 0.1$ after this percolation epoch, which occurs after a period of $approx 10$ Myr for a typical starburst SNe rate density of $ u_{rm SN} approx 10^{-9}$ pc$^{-3}$ yr$^{-1}$ and gas density of $napprox 10$ cm$^{-3}$ in starburst nuclei regions. This matches the recent observations of a time delay of similar order between the onset of star formation and galactic outflows. The efficiency to heat gas up to X-ray temperatures ($ge 10^{6.5}$ K) roughly scales as $ u_{rm SN}^{0.2} n^{-0.6}$. For a typical SNe rate density and gas density in starburst nuclei, the heating efficiency is $sim 0.15$, also consistent with previous interpretations from X-ray observations. We discuss the implications of our results with regard to observational diagnostics of ionic ratios and emission measures in starburst nuclei regions.
We point out that the commonly assumed condition for galactic outflows, that supernovae (SNe) heating is efficient in the central regions of starburst galaxies, suffers from invalid assumptions. We show that a large filling factor of hot ($ge 10^6$ K
) gas is difficult to achieve through SNe heating, irrespective of the initial gas temperature and density, and of its being uniform or clumpy. We instead suggest that correlated supernovae from OB associations in molecular clouds in the central region can drive powerful outflows if the molecular surface density is $> 10^3$ M$_{odot}$ pc$^{-2}$.
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 ener
getics 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 consider the rate of ionization of diffuse and molecular clouds in the interstellar medium by Galactic cosmic rays (GCR) in order to constrain its low energy spectrum. We extrapolate the GCR spectrum obtained from PAMELA at high energies ($ge 200$
GeV/ nucleon) and a recently derived GCR proton flux at $1hbox{--}200$ GeV from observations of gamma rays from molecular clouds, and find that the observed average Galactic ionization rate can be reconciled with this GCR spectrum if there is a low energy cutoff for protons at $10hbox{--}100$ MeV. We also identify the flattening below a few GeV as being due to (a) decrease of the diffusion coefficient and dominance of convective loss at low energy and (b) the expected break in energy spectrum for a constant spectral index in momentum. We show that the inferred CR proton spectrum of $Phi propto E_{kin}^{-1.7pm0.2}$ for $E_{kin} le$ few GeV, is consistent with a power-law spectrum in momentum $p^{-2.45pm0.4}$, which we identify as the spectrum at source. Diffusion loss at higher energies then introduces a steepening by $E^{-alpha}$ with $alpha sim 1/3$, making it consistent with high energy measurements.
