Creasey, Peter; Theuns,Tom; Bower, Richard G.; How supernova explosions power galactic winds

How supernova explosions power galactic winds

by: Creasey, Peter; Theuns,Tom; Bower, Richard G.

[ADS:2013MNRAS.429.1922C; pdf; First author; second; third]
I thought that this was an interesting paper. Because it is very similar to what I am doing it gave me a number of ideas about how to set up specific tests and a way of organizing my data to pick out the important details. This paper is a subset of Peter Creasey's PhD work so the general approach and specific questions that he is exploring are different from mine, but on some level we are doing similar work.

He is using the MHD code FLASH, but only in the hydro configuration (like I am doing with Athena). Even though FLASH has the great advantage of being a AMR code they ran into some problems where the AMR was trying to refine the simulation beyond what was practical so they turned off the AMR component of FLASH and used it as a fixed grid solver (which automatically removed the major advantage of using FLASH, and basically turned FLASH into a FORTRAN version of Athena).

When I read that in the paper I instantly knew why they had done that. Because of the increased resolution from the AMR the time step in the simulation would drop to something incredibly small. Thus to do a 5 Myr simulation like they are doing, it would take a VERY LONG TIME to do a single simulation. I would estimate that a single run may take several months on 128+ processors. By turning off the AMR they could force a coarser resolution and thus take larger time steps, and thus run a number of simulations in less time (I think they did 61 simulations if I remember correctly).

Their resolution for their simulations ranged from 32x32x160 to 256x256x1280 over a spacial range of 200x200x1000 pcs. Thus their highest resolution runs had about 50% fewer cells than my highest resolution runs, but my simulations are spread out over a box of 1000 pc³, so I have slightly lower spacial resolution.

They are looking at how supernova feedback affects mass loading which they define as β≡M˙_wind/M˙_⋆ (those are supposed to be M dots, as in change of Mass). This is a measure of how much mass is flowing out of the grid over how much gets converted into new stars (oh and as a note, for this major important equation they reference Stringer et al. (2011), but there is no reference to Stringer et al. (2011) in the bibliography. There is a Stringer et al. (2012), but even though that paper is the intended paper it does not have the β in it as advertised. Stringer may use another symbol but they don't use β in any of their equations. Anyway minor thing.) Basically they are looking at how much gas escapes the galaxy based on global properties of the galaxy which can then be matched to larger simulations where galactic outflow of individual galaxies is important.

Below is Figure 8 from the paper with the original caption:

Figure 8. Matrix view of simulations varying gas surface density (Σ_g) and gas fraction (f_g), each panel showing a time-averaged vertical velocity for the upper half plane of each simulation (i.e. the disc is at the base of each panel). Gas surface density increases from left to right, gas fraction increases from bottom to top. There appears to be a strong trend in wind velocity towards the lower right-hand panels, i.e. a disc with low gas fraction but high gas surface density tends to generate a faster wind.

I have two comments. First: Cool. That is interesting. Second: AAAAGGGHHHH!!!!!! Rainbow color map!! Why?!?!? AND they used it with a diverging scale, but they cut off the bottom half of each grid so they didn't even need half of their color map. Use something else! Not the rainbow color map. I had to stare at this plot for several minutes to even figure out what they were showing. If they had used something else like a diverging scale, or an incremental luminous scale then it would have been much easier. Anyway, they had other things to worry about.

Wagner, A. Y.; Bicknell, G. V.; Umemura, M.; Driving Outflows with Relativistic Jets and the Dependence of Active Galactic Nucleus Feedback Efficiency on Interstellar Medium Inhomogeneity

Driving Outflows with Relativistic Jets and the Dependence of Active Galactic Nucleus Feedback Efficiency on Interstellar Medium Inhomogeneity

by: Wagner, A. Y.; Bicknell, G. V.; Umemura, M.

[arXiv:1205.0542, pdf, First author, second, third]
I noticed this paper because the second author on it has also written several papers with Jackie Cooper, whose research forms the jumping off point for my own work. There are many things about the set up to this problem that mirror what was already done by Cooper and Bicknell in previous papers. So they drew upon years of previous work to create this set of simulations that form the foundation of this paper. Having taken all the tools and mechanics that were developed previously all they had to do was modify it for this specific problem and let it run. No need to reinvent the wheel. If you take a glance at this paper you may not realize the large amount of work that went into setting up this problem by many other people before this problem was solved. This is the way almost all science works and this is how I have done a lot of my work. As one of my former professors would say, "A lazy physicist is a good physicist!"

They use the code FLASH in its relativistic hydrodynamics configuration to do a series of simulations of relativistic jets streaming off of an AGN and into the local ISM. They use a fractal distribution to create the initial density of the warm phase ISM (~10⁴K). The rest is hot gas (~10⁷K). Their study consists of 29 different simulations, where 15 of those were done for this paper and the rest were included in a previous paper. In their different simulations they tested different parameters such as the jet pressure, the density of the hot phase medium, average density of warm phase, volume filling factor of warm phase, maximum cloud size, total mass in warm phase, and other parameters that depend of the sampling wavenumber (effectively a measure of the fractal scale relative to the overall size of the fractal cube).

In this paper they are looking at the potential for galactic feedback, effectively a measure of how much of the energy in the AGN gets transferred into the warm ISM, thereby dispersing the clouds. They find that the feedback is sensitive to the maximum cloud size but not so much the filling factor. In order to determine the veracity of their findings they compare the velocity dispersion of the ablated clouds with observations of galaxies with radio lobes (i.e. those that have powerful AGN jets).

I thought that this paper was interesting since it deals with something close to what I am doing, with similar problems, methods and solutions. I will definitely use it as a reference for references, ideas and ways of organizing my own simulations. There are some questions that they explored here and things that they calculated or measured that I will definitely consider including in my own work.

Barger, K.A.; Haffner, L.M.; Bland-Hawthorn, J.; Warm Ionized Gas Revealed in the Magellanic Bridge Tidal Remnant: Constraining the Baryon Content and the Escaping Ionizing Photons around Dwarf Galaxies

Warm Ionized Gas Revealed in the Magellanic Bridge Tidal Remnant: Constraining the Baryon Content and the Escaping Ionizing Photons around Dwarf Galaxies

by: K. A. Barger, L. M. Haffner, J. Bland-Hawthorn

[arXiv:1305.4932v2, pdf, First author, second, third]

This paper deals with the "bridge" (a remnant of tidal interaction) between the Large Magellanic Cloud (LMC) and the Small Magellanic Cloud (SMC). The bridge in this case is a region of ionized gas that stretches between the two clouds. In the paper they present the results of an Hα survey which they did using the Wisconsin Hα Mapper (WHAM) observatory. Some of the interesting results that they got were that the ionization fraction was higher in the bridge than in the SMC tail (36 − 52% compared to 5 − 24%).

They also found that the amount of ionizing radiation from the Milky Way and from extragalactic sources is insufficient to ionize the bridge. Thus they conclude that there needs to be a small amount of ionizing radiation that leaks out of the SMC and the LMC (about 4-5%). This has implications about how dwarf galaxies affect their surroundings and how much they affect their surroundings.

Figure 1 from Barger et al. The contours show column density of H I.
The contours are at 10, 20, 35, and 50 x 10¹⁹ cm⁻².

I thought that this was a good paper that included a good discussion about how they made their observations and how they had to work with their data in order to subtract off atmospheric and other interference. This is relevant to me because in my simulations I am trying to reproduce some of these features observed here so I need to know what people are observing, how they are observing it and what is possible to observe so that I can fit that with my models.

McClure-Griffiths, N. M.; et al.;Atomic Hydrogen in a Galactic Center Outflow

Atomic Hydrogen in a Galactic Center Outflow

by: N. M. McClure-Griffiths, J. A. Green, A. S. Hill, F. J. Lockman,

J. M. Dickey, B. M. Gaensler, A. J. Green

[arXiv:1304.7538v1, pdf, First author, second, third, fourth]
This is part of the results from a survey where the authors found several cool, T < 4000 K, H I clouds in a definite pattern around the Galactic Center. The distribution and kinematics of the cloud implied that they were created from a starburst located at the center of the galaxy that has been affecting the galaxy for the past ∼ 2 × 10⁶ years.

They did this by calculating the local standard of rest (LSR) velocities for the cloud centers and then using a simulation that calculated the effect on clouds with a random distribution of velocities and comparing to their observations they determined that the observed clouds did not have kinematics similar to galactic rotation but have a distribution that implies ballistic motion away from the galactic center (see figure below).

They also give a bound on the velocity of the stellar winds that can create this type of distribution  "The [Kolmogorov-Smirnov] test implies that wind velocities greater than 270 km s⁻¹ and less than 150 km s⁻¹ are not consistent with the observed cloud velocities."

They conclude that these clouds are the remnants of a superbubble which was created by a starburst at the center of the Milky Way. They also mention survivability and life spans of the clouds. They provide some bounds to the life spans. They also rely heavily on the simulations of Jackie Cooper to reach some conclusions about the environment of the wind and clouds.

Dobbs, C. L.; Pringle, J. E.; The exciting lives of giant molecular clouds

The exciting lives of giant molecular clouds

by: Dobbs, C. L.; Pringle, J. E.

[arXiv:1303.4995v1, pdf, first author, second]
This paper does a particle simulation with 8 million particles with each particle having a mass of ~300 M_sun. They include heating and cooling following the method of Glover and Mac Low (2007), with cooling being switched off at 50K. Each particle contains a fraction of H₂ which they recalculate each timestep based on expected formation and destruction events. They also include stellar feedback with the feedback relation given in equation (1) of their paper. It depends on an arbitrarily chosen parameter (ε), the  H₂ mass, and a factor that comes from their chosen IMF. This is multiplied by 10^51 ergs, which comes from the average energy of a supernova events. The energy is deposited half in thermal energy, half in kinetic. No explanation on how they choose the deposit of kinetic energy, but that may be covered in a previous paper (actually a lot of what they are doing is covered in previous papers, they are always referring to previous papers for their set up). They also had a background stellar potential, with perturbations for the spiral arms.

They ran their simulations for ~300 Myr, and found that giant molecular clouds form in the spiral arms through the conglomeration of smaller clouds. The GMCs are then disrupted mostly through sheer effects, but in the smaller clouds there is more disruption through stellar feedback. The larger clouds can also stretch out into a spur in the spiral arm which in turn can form its own GMCs, from the remnants of the original cloud. Below is figure 1 from the paper which shows column density at 250 Myr. One particular GMC that they studied in detail is marked with a red square.

Below is figure 3 from the paper which shows the evolution of the cloud marked in figure 1.

Melioli,C. et al.; Evolution of M82-like starburst winds revisited: 3D radiative cooling hydrodynamical simulations

Evolution of M82-like starburst winds revisited: 3D radiative cooling hydrodynamical simulations

by: C. Melioli, E. M. de Gouveia Dal Pino, F. Geraissate

[arXiv:1301.5005, pdf, first author, second, third]

This is an interesting paper because it is closely related to work done by Cooper et al. (first paper 2008, second paper 2009) that I have been looking at for some time. The paper deals with simulations done using a hydro AMR that has radiative cooling and some species tracking. It is more work on superbubbles and AGNs. They specifically use M82 as a test case.

The authors are from Brazil (Sao Paulo), and the code is named YGUAZU, which is a Paraguayan spelling of Iguazú (sort of appropriate for a hydro code since it means "Big Water"). Other than some basics (they use a Van Leer integrator) they only provide references and no explanation. Also interesting is the fact that they cite Strickland & Stevens (2000) in their explanation of how they set their initial conditions, but they don't use the notation of Strickland and Stevens. They use the notation of Jackie Cooper (2009) (she did work with Strickland and Stevens and used their code and set up). But these guys don't cite here even though they have copied her equations exactly.

Their energy injection centers around super stellar clusters (SSCs) "with an average size of ∼ 5.7 pc and mass (of stars) between 10⁴ and 10⁶ M_⊙(Melo et al. 2005)." They look at metals and how much gas escapes the galaxy and how much metals produced by supernovas escapes the galaxies. They conclude that most of the gas mass stays in the galaxy even with a superbubble blow out. Also most of the metals stay in the galaxy but some get transported out in the galactic winds that form due to the supernovas (the SN's pump out metal rich winds).

Mac Low, M.-M.; McCray, R.; Superbubbles in disk galaxies

Superbubbles in disk galaxies

by: Mac Low, M.-M.; McCray, R.

[ADS: 1988ApJ...324..776M, pdf, first author: personal site, second author]
This paper is part of a series of papers done by Mac Low and McCray on this subject. This paper contains the theory and analytic backing on the subject, another paper in 1989 contains more of the simulations on the same subject. This paper and others formed the basis of Mac Low's PhD dissertation, "Interactions of Massive Stars with the Interstellar Medium: Bow Shocks and Superbubbles". McCray was his advisor. Mac Low would later go on to advise my advisor Fabian Heitsch for his PhD.

This particular paper is heavy on the theory and equations. Basically they are looking at how supernovas (SNs) interact with the ISM. As supernovas explode they release significant energy into the ISM and create a bubble of hot gas. If the bubble is large enough it gets classified as a superbubble, which has the possibility of blowing out of the galactic disk, which then affects the halo and galactic accretion.

This paper represents a significant step forward in our understanding of the structure of the ISM. The models produced here are idealized and smooth, meaning everything looks nice, flat and symmetric. 25 years later the technology and experience available to researchers has improved and thus we have moved on to solving this exact same problem, except in 3D and with a much more complex setup.

There are two important conclusions that I wanted to mention. They provide a parameter that determines whether or not a superbubble will blow out of a stratified galactic disk. They give it as:

D = L_SN ρ₀^1/2 P_e^-3/2 H^-2

where L_SN is the luminosity from the SNs, ρ₀ is the density of the galactic disk (or ISM), P_e is the external pressure from the ISM, and H is the scale height of the disk. If D > 100 then there will be blowout even if the center of the superbubble begins to collapse. It will be interesting to find a corresponding parameter for a more complex set of simulations.

The second important point is that if there is a dense cloud in the ISM then when the edge of the superbubble over runs it will not "puncture" the bubble leading to a release of pressure. The bubble will instead travel around it and continue expanding. This is something that has become a very important consideration since it is the thing that allows molecular clouds to survive strong shocks like this. This was essentially a hint at the beginning of the study of the survivability of cold molecular clouds when they have been strongly shocked. The problem is that if molecular clouds are strongly shocked then they will heat up, expand and will not collapse gravitationally to form stars. So there has to be some way for them to survive long enough to form stars. Many people will look into this problem later, and research is still going on.

As a note, they used the 2D hydro code Zeus, which was very influential back in the day. The creators of Zeus rewrote the code for MHD and 3D and named it Athena (original site). Athena is the code that I use for my research.