Galaxy interactions and mergers play a significant, but still debated and poorly understood role in the star formation history of galaxies. Numerical and theoretical models cannot yet explain the main properties of merger-induced starbursts, including their intensity and their spatial extent. Usually, the mechanism invoked in merger-induced starbursts is a global inflow of gas towards the central kpc, resulting in a nuclear starburst. We show here, using high-resolution AMR simulations and comparing to observations of the gas component in mergers, that the triggering of starbursts also results from increased ISM turbulence and velocity dispersions in interacting systems. This forms cold gas that are denser and more massive than in quiescent disk galaxies. The fraction of dense cold gas largely increases, modifying the global density distribution of these systems, and efficient star formation results. Because the starbursting activity is not just from a global compacting of the gas to higher average surface densities, but also from higher turbulence and fragmentation into massive and dense clouds, merging systems can enter a different regime of star formation compared to quiescent disk galaxies. This is in quantitative agreement with recent observations suggesting that disk galaxies and starbursting systems are not the low-activity end and high-activity end of a single regime, but actually follow different scaling relations for their star formation.
Star formation in galaxy mergers]
Star formation in galaxy mergers:
ISM turbulence, dense gas excess, and scaling relations for disks and starbusts Bournaud et al.] Frédéric Bournaud, Leila C. Powell, Damien Chapon and Romain Teyssier 2010 \volume271 \pagerange \jnameAstrophysical Dynamics: From Stars to Galaxies \editorsN. Brummell, A. S. Brun, M. S. Miesch & Y. Ponty
1 Introduction: Galaxy mergers as triggers of star formation
Evidence that interactions and mergers can strongly trigger the star formation activity of galaxies has been observed for more than two decades (e.g., Sanders et al., 1988). Nevertheless, both the underlying mechanisms the role of mergers in the cosmic budget of star formation remain largely unknown. Mergers have long appeared to potentially dominate the star formation history of the Universe. Starbursting galaxies with specific star formation rates much above the average, such as Ultraluminous Infrared Galaxies (ULIRGs), are almost exclusively major interactions and mergers – in the nearby Universe (see Duc et al. 1997), not necessarily at high redshift. The decrease in the cosmic density of star formation between and (Le Floc’h et al., 2005) may follow, and might even result from, the decrease in the galaxy merger rate. Other studies based on disturbed UV morphologies and/or kinematics also suggested high merger fractions among star forming galaxies at or even below (Hammer et al., 2005).
However, disturbed morphologies and kinematics can also arise internally, without interactions, especially at high redshift when disk galaxies are wildly unstable, clumpy and irregular (Elmegreen et al., 2007; Förster Schreiber et al., 2009). Merger-induced ULIRG-like starbursts may also be less important in the cosmic budget than the more moderate but more numerous objects with internally sustained star-formation, as is possibly the case for most high-redshift LIRGs (e.g., Daddi et al., 2010a). Recent studies aimed at accurately distinguishing the signatures of mergers from internal evolution actually suggest that major interactions and mergers account only for a small fraction of the cosmic star formation history (Jogee et al., 2009; Robaina et al., 2009). Moreover, observational estimates appear to be quite dependent on the chosen tracers of star formation. For instance, mergers are probably a more important trigger of dust-obscured star formation than of general star formation, and thus the fraction of mergers among objects with high infrared-traced star formation rate can be higher (as discussed by Robaina et al., 2009).
Numerical simulations are then required to understand the mechanisms leading to starbursting activity in mergers, but also to further probe the contribution of mergers to cosmic star formation since observational estimates remain uncertain. Section 2 reviews the standard knowledge on merger-induced star formation, which is mostly based on “sub-grid” modeling: all steps from the formation of cold/dense gas clouds to the formation of actual stars remain unresolved, and described with arbitrary recipes. This standard understanding cannot account for some general observationed features that we briefly review. Section 3 presents a new generation of high-resolution models in which the first steps of galactic-scale star formation, namely ISM turbulence and cold/dense gas cloud formation, are explicitly resolved using high-resolution codes. Based on this, we provide a substantially different explanation for interaction-triggered star formation and show that it could better account for recent observations of disk galaxies and starbursting mergers.
2 Standard models versus observations
2.1 Merger-induced gas flows
In an interacting galaxy pair, the gas content of one galaxy undergoes tidal forces from the companion galaxy, but this is actually not the main direct driver of the gas response. The gas distribution becomes non-axisymmetric, because of the asymmetry in the gravitational field itself induced by the companion. This results in gravitational torquing of the gas. - a thorough description of the mechanism can be found in Bournaud (2010, Section 2 and Figure 1). Gas initially inside the corotation radius (typically a radius of a few kpc) undergoes negative gravity torques and flows inwards in a more and more concentrated central component. Gas outside the corotation, i.e. initially in the outer disk, undergoes positive gravity torques and gains angular momentum, forming long tidal tails.
2.2 Nuclear starbursts
Gravitational torquing in the inner disk increases the gas concentration in the central kiloparsec or so. Any model for star formation will then predict an increase of the star formation rate (global Schmidt-Kennicutt law, models based on cloud-cloud collisions, etc). The result is thus a centralized or nuclear starburst. As the driving process is gravitational torquing, early restricted three-body models could already describe the effect (Toomre & Toomre, 1972). Later models have added extra physical ingredients leading to more accurate predictions on the star formation activity: self-gravity (Barnes & Hernquist, 1991), hydrodynamics and feedback processes (Mihos & Hernquist, 1996; Cox et al., 2008), etc.
A large library of SPH simulations of galaxy mergers, in which the driving process is mostly the one presented above (tidal torquing of gaseous disk) was performed and analyzed by (Di Matteo et al., 2007, 2008). This study highlighted various statistical properties of merger-induced starbursts. In particular, it showed that some specific cases can lead to very strong starbursts with star formation rate (SFRs) are increased by factors of 10–100 or more, but that on average the enhancement of the SFR in a random galaxy collision is only a factor of a few (3–4 being the median factor) and only lasts 200-400 Myr. These results were confirmed with code comparisons, and found to be independent of the adopted sub-grid model for star formation (Di Matteo et al., 2008). Models including an external tidal field to simulate the effect of a large group or cluster found that the merger-induced starbursts could be somewhat more efficient in such cases – but the SFR increase remains in general below a factor of 10 (Martig & Bournaud, 2008).
2.3 Theoretical predictions versus observations
The intensity of merger-induced starbursts Numerical simulations reproducing the interaction-induced inflow of gas and resulting nuclear starbursts can sometimes trigger very strong starbursts, but in general the SFR enhancement peaks at 3–4 times the sum of the SFRs of the two pre-merger galaxies. This factor of 3–4 seems in good agreement with the most recent observational estimates (Jogee et al., 2009; Robaina et al., 2009, e.g.). In fact, there is substantial disagreement: the factor 3–4 in simulation samples is the peak amplification of the SFR in equal-mass mergers. In observations, it is the average factor found at random (observed) instant of interactions, and in mergers that are ”major” ones but not strictly equal-mass ones. Given that typical duration of a merger is at least twice longer than the starburst activity in the models, and that unequal-mass mergers make substantially weaker starbursts (Cox et al., 2008), one would need a peak SFR enhancement factor of about 10–15 (as measured in simulations) to match the average enhancement of 3–4 found in observations. There is thus a substantial mismatch between the starbursting activity predicted by existing samples of galaxy mergers, and that observed in the real Universe – although the observational estimates remain debated and may depend on redshift.
The spatial extent of merger-induced starbursts
A disagreement between these ’standard’ models/theories and observations of mergers is also found in the geometry and spatial extent of star formation. The very strongest merger-induced starbursts (ULIRGs) are centrally-concentrated, but there are many examples of significant merger-induced starbursts that are spatially extended. A well-known example are the Antennae galaxies, where the burst of star formation proceeds, for a large fraction, in a few big star-forming clumps in extended disks and in the bridge between the two galaxies (Wang et al., 2004). Another well-known example is the IC2163/NGC2207 pair, which has an extended starburst in big gas clouds with remarkably high gas velocity dispersions (Elmegreen et al., 1995). There are many other examples of extended star formation in merging systems (Cullen et al., 2006; Smith et al., 2008).
Quantitative comparisons of the extent of star formation in observations to that predicted by ”standard” models have shown a significant disagreement Barnes (2004); Chien & Barnes (2010). These authors also suggested that a sub-grid model of shock-induced star formation may better account for the spatial extent of merger-induced star formation (see also Saitoh et al. 2009).
The standard mechanism for merger-induced star formation, as reproduced in low-resolution simulations, certainly takes place in real mergers – signs of nuclear merger-driven starbursts are a plenty. But it seems impossible to explain the typical intensity of merger-induced starbursts and their often relatively larger spatial extend, based on these ”standard” models. Some relatively old observations may actually hold the key to correctly understanding merger-induced star formation:
2.4 A key observation: ISM turbulence in interacting galaxies
The cold gas component of interacting galaxies has high velocity dispersions, which can reach a few tens of km s. This has been noticed in the 90s in at least two well-resolved interacting and merging pairs (Irwin, 1994; Elmegreen et al., 1995). Many interacting pairs and mergers in Green et al. (2010) also have large velocity dispersions (in the ionized gas). More generally, massive star clusters in merging systems suggest the Jeans mass is high, which is indicative of high velocity dispersions. A typical number could be a factor 4 of increase for the cold gas velocity dispersion in equal-mass mergers. An interesting question will be understand whether these dispersions result from the tidal interaction and are a trigger star formation, or whether they just result from the starburst and associated feedback effects.
Traditional SPH simulations model a relatively warm gas for the ISM, because the limited spatial resolution translates into a minimal temperature under which gas cooling should not be modeled (it would generate artificial instabilities, Truelove et al. 1997). Modeling gas cooling substantially below K requires ”hydrodynamic resolutions”
Increased ISM turbulence in galaxy mergers is thus absent from the modeling used in most hydrodynamic simulations to date. Some particle-based models have nevertheless been successful in reproducing these increased gas dispersions (Elmegreen et al., 1993; Struck, 1997; Bournaud et al., 2008), indicating that it is a consequence of the tidal interaction which triggers non-circular motions, rather than a consequence of starbursts and feedback. It should then arise spontaneously in hydrodynamic models, if these a capable of modeling gas below K.
3 Merger-induced starbursts with resolved ISM turbulence and cluster star formation
3.1 High-resolution AMR simulations
Adaptive Mesh Refinement (AMR) codes allow hydrodynamic calculations to be performed at very high resolution on adaptive-resolution grids. The resolution is not high everywhere, but the general philosophy is to keep the Jeans length permanently resolved until the smallest cell size is reached. That is, the critical process in the collapse of dense star-forming clouds, namely the Jeans instability (or Toomre instability in a rotating disk) is constantly resolved up to a typical scale given by the smallest cell size, or a small multiple of it.
AMR simulations of whole galaxies have recently reached resolutions of a few pc for disk galaxies (Agertz et al., 2009b; Tasker & Tan, 2009), and even 0.8 pc lastly (Bournaud et al., 2010). Such techniques have been first employed to model ISM dynamics and star formation in galaxy mergers by Kim et al. (2009); Teyssier et al. (2010).
We here study the properties of star formation in a sample of a few AMR simulations of 1:1 mergers of Milky Way-type spiral galaxies, performed with a resolution of 4.5 pc and a barotropic cooling model down to 50 K, technically similar to the isolated disk simulation described in (Bournaud et al., 2010). Star formation takes place above a fixed density threshold and is modeled with a local Schmidt law, i.e. the local star formation rate density in each grid cell is where is the local gas density and the gravitational free-fall time, and the efficiency is a fixed parameter. Supernova feedback is included. Further details and results for whole sample will be presented in Powell et al. (in preparation). An individual model of this type (but at lower resolution and without feedback), matching the morphology and star formation properties of the Antennae galaxies, was presented in Teyssier et al. 2010.
3.2 Starburst properties
In the models, the pre-merger isolated spiral galaxies spontaneously develop ISM turbulence at a about 10 km s under the effect of gravitational instabilities (and/or feedback), and most star formation takes place in dense complexes of dense gas along spiral arms. In some sense, the large-scale star formation process is not entirely sub-grid anymore in these simulations, as the first steps of star formation, namely are the development of ISM turbulence and the formation of dense molecular gas clouds in this turbulent ISM, are explicitly captured – the subsequent steps of star formation at smaller scales, inside the densest parts of these cold clouds, remain sub-grid.
A merger simulation is shown in Figure 3.1. The mass-weighted average of the gas velocity dispersion reaches km/s. This strong turbulence is consistent with the observations reviewed above. It induces numerous local shocks that increase the local gas density, which in turn triggers the collapse of gas into cold clouds. Also, gas clouds become more massive and denser than in the pre-merger spiral galaxies. The fraction of gas that is dense-enough and cold-enough to form star increases, and the timescale for star formation in these dense gas entities (the gravitational free-fall time) becomes shorter. As a result, the total SFR becomes several times higher than it was in the pre-interaction pair of galaxies. The standard process of merger-induced gas inflow towards the central kpc or so is also present, but the timescale is substantially longer, so this process dominates the triggering the star formation by enhancing the global gas density only in the late stages of the merger (Fig. 3.1).
An example of star formation distribution is shown in Figure 3.2. Two consequences of modeling a cold turbulent ISM in merging galaxies are: (1) the peak intensity of the starburst can become stronger (as shown by Teyssier et al. 2010 for the Antennae) although this is not a systematical effect, and (2) the spatial extent of star formation in the starbursting phase is larger. The radius containing 50% of the star formation rate (half-SFR radius) can more than double. This is because increased ISM turbulence is present throughout the disk, and triggers, through locally convergent flows and shocks, the collapse of efficiently star-forming clouds even at large distances from the nuclei. At least quantitatively, these results put the models in better agreement with observations. In these models, the increased ISM turbulence is also obtained without feedback, showing that it is not a consequence of the starbursting activity, but is rather driven by the tidal forces in the interaction.
3.3 Dense gas fractions and star formation scaling relations
Density PDFs and galactic-scaled star formation
The connection between ISM turbulence, dense gas phases and star formation is best understood by examining the probability distribution function of the local gas density (density PDF). In its mass-weighted (resp. volume-weighted) version, the density PDF represents the mass (volume) fraction of the ISM in bins of density. We here use mass-weighted versions.
In an system at equilibrium (e.g. isolated disk galaxy), supersonic turbulence generates a log-normal density PDF (Wada et al., 2002). Only small mass fractions of the whole ISM are at found number densities below cm or in very high density regions above, say, cm, most mass being at ten to a few hundreds of atoms per cm (typically observed as HI or CO-traced molecular gas). The width of the log-normal PDF depends on various factors, but at first order the main dependence is on the turbulent Mach number (Krumholz & Thompson, 2007). A more turbulent ISM will have a larger density spread in its density PDF, because the turbulent flows will sometimes diverge and create low-density holes, and sometimes converge (or even shock) into very dense structures that can further collapse gravitationally. Note that log-normal density PDFs are observed, at least at the scale of molecular complexes (Lombardi et al., 2010).
Our simulations of Milky Way-like disk galaxies do have quasi log-normal PDFs (see example on Figure 3). The PDFs are not exactly log-normal because the radial density gradient in the disk induces some deviations (the PDF is expected to be exactly log-normal only at fixed average density). Also, the spatial resolution limit converts into a density limit at which the density PDF is truncated, and which corresponds to the smallest/densest entities that can be resolved. The maximal density resolution in the present merger models is around cm, but parsec-scale resolution can capture even higher densities (as in Bournaud et al. 2010).
The density PDF is a useful tool to describe the star-forming activity of a galaxy at a given gas content. For a simplified description, one can consider that star formation takes place only in the densest gas phases (i. e. above some density threshold), and that in these dense regions the local star formation rate follows, for instance, a fixed efficiency per free-fall time: (see detailed theory in Elmegreen 2002 and Krumholz & Thompson 2007). Even if the local rate of star formation follows a different prescription than this purely density-dependent model (which is physically motivated by the gravitational collapse timescale), then the first step remains that star formation proceeds only in the densest gas phases. Thus, the fraction of dense gas remains a key parameter for the global star formation activity of galaxies.
Dense gas excess and starbursts in mergers
Representative density PDFs are shown in our models on Figure 3.3 for an isolated spiral galaxy, a moderately starbursting merger (SFR increased by a factor compared to the two pre-merger galaxies taken separately), and a stronger merger-induced starburst (factor ).
The PDFs of merging galaxies have a larger width than those isolated disk galaxies, as expected from the higher turbulent speeds that result from the tidal interaction. These PDFs are not necessarily log-normal
High fractions of dense gas in mergers were already proposed by Juneau et al. (2009), based on detailed post-processing of merger simulations aimed at re-constructing dense molecular gas phases not resolved in these simulations. Here we obtain a qualitatively similar conclusion using simulations that explicitly resolve turbulent motions, local shocks and small-scale instabilities in cold ISM phases. This excess of dense gas should have signatures in molecular line ratios (as also observed by Juneau et al. 2009). If, in a rough approach, we assume that low- CO lines are excited for densities of 100 cm and above, and HCN lines for densities of cm and above, then the HCN/CO line ratios could be up to 5-10 times higher in the starbursting phases of major galaxy mergers. Simulations with an somewhat higher resolution would actually be desirable to accurately quantify the emission of dense molecular tracers.
Star formation scaling relations
The interpretation of merger-induced starbursts proposed from our high-resolution models is that it is not just a global gas inflow that increases the average gas density and increases the SFR, but also that there are strong non-circular motions, high turbulent velocity dispersions, causing many small-scale convergent flows and local shocks, that in turn initiate the collapse of dense star-forming clouds with high Jeans masses. The former ”standard” process does take place, but the later can be equally important especially in the early phases of mergers.
We here note the average gas surface density of a galaxy. This is the quantity that observers would typically derive from the total gas mass and half-light radius, or similar quantities. The second mechanism above is a way to increase the SFR of a system, and its SFR surface density , without necessarily increasing its average . Actually in our merger models does increase (as there are global merger-induced gas inflows), but increases in larger proportions (as the starburst is not just from the global merger-induced inflow but also from the exacerbated fragmentation of high-dispersion gas). Going back to the density PDFs shown previously, one can note that the fraction of very dense gas (say, in the cm range) can increase by a factor of 10–20 in mergers while the average surface density increases by a factor 3–5 (see also on Figure 4). As a consequence, the activity of these systems is unexpectedly high compared to their average surface density .
Figure 4 shows the evolution of a system throughout a merger simulation in the plane. While our pre-merger spiral galaxy models lie on the standard Kennicutt relation, starbursting mergers have high ratios. This is in agreement with observational suggestions that quiescent disks and starbursting mergers do not follow the same scaling relations for star formation, but could actually display two different star formation ”laws” (Daddi et al., 2010b; Genzel et al., 2010). The offset between the disk and merger sequences proposed by Daddi et al. (2010b) is quantitatively recovered in our simulations (Figure 4). Post-starburst, post-merger systems lie back on the quiescent sequence, or even somewhat below it: these systems contain some dense gas which is somewhat stabilized by the stellar spheroid. This is another example of a ”morphological quenching” effect in early-type galaxies (Martig et al., 2009), and the location of our post-merger early-type galaxies in the diagram may be consistent with observations of nearby ellipticals (Crocker et al. 2010, but see Fabello et al. 2010).
The proposal that disks and mergers follow two different regimes of star formation by Daddi et al. (2010b) and Genzel et al. (2010) relies for a part (but not entirely) on the assumption that different CO luminosity-to-molecular gas mass conversion factors apply in quiescent disks and starbursting mergers. Interestingly, our simulations recover the two regimes of star formation without any assumption on such conversion factors since gas masses are directly known. But at the same time, excess of dense gas found in these merger models suggests that the excitation of CO lines would naturally be higher in mergers/starbursting phases (although this needs to be quantified in the models), which would mean that the assumption of different conversion factors by Daddi et al. and Genzel et al. could be physically justified. High molecular gas excitation in SubMillimeter Galaxies (SMGs, Tacconi et al. 1998) could then naturally result if these are starbursting major mergers with high gas surface densities and a clumpy turbulent ISM (e.g. Narayanan et al., 2010; Bournaud et al., 2011).
The results presented here were based mostly on low-redshift merger simulations. A recent set of high-redshift merger simulations with AMR is presented in Bournaud et al. (2011), and shows similar increase in the ISM velocity dispersions and clumpiness in mergers, with extended starbursts and high ratios.
Galaxy interactions and mergers play a significant but still debated and poorly understood role in the star formation history of galaxies. Numerical and theoretical models have significant difficulties in accounting for the properties of merger-induced starbursts, including their intensity and their spatial extent. Usually, the mechanism invoked to explain the triggering of star formation by mergers is a global inflow of gas towards the central kpc, resulting in a nuclear starburst. We show here, using high-resolution AMR simulations and comparing to observations of the gas component in mergers, that the triggering of starbursts also results from increased ISM turbulence and velocity dispersions in interacting systems. This results in the formation and collapse of dense and massive gas clouds in the regions of convergent flows and local shocks, these clouds being denser, more massive and/or more numerous than in quiescent disk galaxies. The fraction of dense cold gas largely increases, modifying the global density distribution of these systems, and efficient star formation results. Because the starbursting activity is not just from a global compacting of the gas to high average surface densities, but also from higher turbulence and fragmentation into massive and dense clouds, merging systems can enter a different regime of star formation compared to quiescent disk galaxies. This is in quantitative agreement with recent observations suggesting that disk galaxies and starbursting systems are not the low activity end and high activity end of a single regime, but actually follow different scaling relations for their star formation.
Acknowledgements.FB acknowledges numerous discussions on gas dynamics and star formation in galaxy mergers with Bruce Elmegreen, Pierre-Alain Duc, Emanuele Daddi and many others, and is grateful to the organizers of the symposium for the exciting discussions along a wide range of astrophysical problems.
- i.e., the average smoothing length in SPH codes, or the smallest cell size in AMR codes
- presumably because the turbulent motions do not follow a quasi-isotropic cascade and/or some of the involved non-circular motions are not really cascading turbulence
- Agertz, O., et al. 2009a, MNRAS, 392, 294
- Barnes, J. E., & Hernquist, L. E. 1991, ApJL, 370, L65
- Barnes, J. E. 2004, MNRAS, 350, 798
- Bournaud, F., Duc, P.-A., & Emsellem, E. 2008, MNRAS, 389, L8
- Bournaud, F. 2010, ASP Conference Series, 423, 177, arXiv:0909.1812v2
- Bournaud, F., Elmegreen, B.G., Teyssier, R., Block, D.L., & Puerari, I. 2010, MNRAS 409, 1088
- Bournaud, F., et al. ApJ submitted, arXiv:1006.4782
- Burkert, A. 2006, Comptes Rendus Physique, 7, 433
- Chien, L.-H., & Barnes, J. E. 2010, MNRAS, 407, 43
- Cox, T. J., Jonsson, P., Somerville, R. S., Primack, J. R., & Dekel, A. 2008, MNRAS, 384, 386
- Crocker, A. F., Bureau, M., Young, L. M., & Combes, F. 2010, submitted to MNRAS
- Cullen, H., Alexander, P., & Clemens, M. 2006, MNRAS, 366, 49
- Di Matteo, P., Combes, F., Melchior, A.-L., & Semelin, B. 2007, A&A, 468, 61
- Di Matteo, P., Bournaud, F., Martig, M., Combes, F., Melchior, A.-L., & Semelin, B. 2008, A&A, 492, 31
- Daddi, E., et al. 2010a, ApJ, 713, 686
- Daddi, E., et al. 2010b, ApJL, 714, L118
- Duc, P.-A., Brinks, E., Wink, J. E., & Mirabel, I. F. 1997, A&A, 326, 537
- Elmegreen, B. G., Kaufman, M., & Thomasson, M. 1993, ApJ, 412, 90
- Elmegreen, D. M., Kaufman, M., Brinks, E., Elmegreen, B. G., & Sundin, M. 1995, ApJ, 453, 100
- Elmegreen, B. G. 2002, ApJ, 577, 206
- Elmegreen, D. M., Elmegreen, B. G., Ravindranath, S., & Coe, D. A. 2007, ApJ, 658, 763
- Le Floc’h, E., et al. 2005, ApJ, 632, 169
- Fabello, S., et al. 2010, MNRAS in press, arXiv:1009.4309
- Förster Schreiber, N. M., et al. 2009, ApJ, 706, 1364
- Genzel, R., et al. 2010, MNRAS, 407, 2091
- Green, A. W., et al. 2010, Nature, 467, 684
- Hammer, F., Flores, H., Elbaz, D., Zheng, X. Z., Liang, Y. C., & Cesarsky, C. 2005, A&A, 430, 115
- Irwin, J. A. 1994, ApJ, 429, 618
- Jogee, S., et al. 2009, ApJ, 697, 1971
- Juneau, S., Narayanan, D. T., Moustakas, J., Shirley, Y. L., Bussmann, R. S., Kennicutt, R. C., & Vanden Bout, P. A. 2009, ApJ, 707, 1217
- Kim, J.-h., Wise, J. H., & Abel, T. 2009, ApJL, 694, L123
- Krumholz, M. R., & Thompson, T. A. 2007, ApJ, 669, 289
- Lombardi, M., Lada, C. J., & Alves, J. 2010, A&A, 512, A67
- Martig, M., & Bournaud, F. 2008, MNRAS, 385, L38
- Martig, M., Bournaud, F., Teyssier, R., & Dekel, A. 2009, ApJ, 707, 250
- Mihos, J. C., & Hernquist, L. 1996, ApJ, 464, 641
- Narayanan, D., et al. 2010, MNRAS, 401, 1613
- Robaina, A. R., et al. 2009, ApJ, 704, 324
- Saitoh, T. R. et al. 2009, PASJ, 61, 481
- Sanders, D. B., et al. 1988, ApJ, 325, 74
- Smith, B. J., et al. 2008, AJ, 135, 2406
- Struck, C. 1997, ApJS, 113, 269
- Tacconi, L. J., et al. 2008, ApJ, 680, 246
- Tasker, E.J., & Tan, J.C. 2009, ApJ, 700, 358
- Teyssier, R., Chapon, D., & Bournaud, F. 2010, ApJL 720, 149
- Toomre, A., & Toomre, J. 1972, ApJ, 178, 623
- Truelove, J. K., et al. 1997, ApJl, 489, L179
- Wada, K., Meurer, G., & Norman, C.A. 2002, ApJ, 577, 197
- Wang, Z., et al. 2004, ApJS, 154, 193