The Spatial Extent and Distribution of Star Formation in 3D-HST Mergers at
We present an analysis of the spatial distribution of star formation in a sample of 60 visually identified galaxy merger candidates at . Our sample, drawn from the 3D-HST survey, is flux-limited and was selected to have high star formation rates based on fits of their broad-band, low spatial resolution spectral energy distributions. It includes plausible pre-merger (close pairs) and post-merger (single objects with tidal features) systems,Êwith total stellar masses and star formation rates derived from multi-wavelength photometry. Here we use near-infrared slitless spectra from 3D-HST which produce H or [Oiii] emission line maps as proxies for star-formation maps. This provides a first comprehensive high-resolution, empirical picture of where star formation occurred in galaxy mergers at the epoch of peak cosmic star formation rate. We find that detectable star formation can occur in one or both galaxy centres, or in tidal tails. The most common case (58%) is that star formation is largely concentrated in a single, compact region, coincident with the centre of (one of) the merger components. No correlations between star formation morphology and redshift, total stellar mass, or star formation rate are found. A restricted set of hydrodynamical merger simulations between similarly massive and gas-rich objects implies that star formation should be detectable in both merger components, when the gas fractions of the individual components are the same. This suggests that mergers typically occur between galaxies whose gas fractions, masses, and/or star formation rates are distinctly different from one another.
keywords:galaxies: formation – galaxies: interactions – galaxies: starburst – galaxies: structure
The spatial extent and distribution of star formation in normal, local galaxies is well established (e.g., James et al., 2004; Bigiel et al., 2008, 2011; Schruba et al., 2011; Calzetti et al., 2012, and references therein). Both self-regulated star formation and merging must be key ingredients in galaxy-formation and evolution and have been studied observationally in detail in the Universe (e.g., Lambas et al., 2003; Hammer et al., 2005; Barton et al., 2007, 2000; Jogee et al., 2009; Robaina et al., 2009, 2010) as well as in theoretical simulations (e.g., Mihos & Hernquist, 1996; Barnes & Hernquist, 1996; Springel, 2000; Cox, 2004; Cox et al., 2006, 2008; Lotz et al., 2008a, b; di Matteo et al., 2007). From these studies, it has become evident that rapid star formation can be triggered by tidal interaction in mergers, but the simulations also suggest that mergers trigger both nuclear starbursts and black hole accretion. Even though galaxy mergers are observed to enhance star formation in galaxies and trigger some of the most violent starbursts known (e.g., Rieke et al., 1985; Joseph & Wright, 1985; Melnick & Mirabel, 1990; Klaas & Elsaesser, 1991), it appears that the net effect of major mergers in the global star formation history of the galaxy population has been relatively modest since at least (e.g., Robaina et al., 2009; Jogee et al., 2009).
Mergers are predicted to play a crucial role in the build-up and formation of massive galaxies (e.g., Mihos & Hernquist, 1996; Springel, 2000; Cox et al., 2008; Hopkins et al., 2010), and therefore, a crucial step towards fully understanding how galaxies have evolved is to study the star formation properties in merging systems at high redshifts. Some of the questions that need to be addressed observationally to get a more detailed picture of the star formation history of present-day galaxies are: Where did stars form in merging galaxies at higher redshift? and Which phase(s) of the merging process seems to trigger most star formation?
At higher redshift, the interplay between merging and star formation has been investigated in much less detail than in the local Universe. This is mostly due to observational challenges, as resolved observations at resolution are needed. Further, tracing star formation, e.g., through H at requires observations in the near-infrared (NIR), and these have been far less feasible than observations in optical wavebands. Nevertheless, the epoch is immensely important for understanding galaxy-formation, as this is the cosmic time when the majority of the stars we see in galaxies today were formed (e.g., Hopkins & Beacom, 2006; Karim et al., 2011). Several studies have addressed the impact of (major) galaxy mergers as well as the general galaxy morphology on the amount of star formation at (e.g., Swinbank et al., 2004; Law et al., 2007; Förster Schreiber et al., 2009; Wright et al., 2009; Conselice et al., 2011; Bluck et al., 2012; Bell et al., 2012). However, only a few of these studies investigate three-dimensional (3D) spectroscopy, where both spatial and spectral information is available, which is crucial for investigating the spatial extent of star formation.
Until recently, large samples of galaxies, and in particular galaxy mergers, with rest-frame optical 3D spectroscopic information at high redshift did not exist. ÊThe largest samples of galaxies with such data at are from the SINS survey (Förster Schreiber et al., 2009), recently expanded with the zC-SINF sample (Mancini et al., 2011), totalling 110 star-forming galaxies at observed with SINFONI, and from the MASSIV survey (Contini et al., 2012) which contains 84 star-forming galaxies at also mapped with SINFONI. Using H as kinematic and star formation tracer enabled analysis of the spatially-resolved ionised gas kinematics, its distribution, and the physical properties of these systems (Genzel et al., 2006, 2008, 2011; Shapiro et al., 2008; Cresci et al., 2009; Epinat et al., 2009, 2012; Queyrel et al., 2009, 2012; Newman et al., 2012). High-resolution NIR imaging with The Hubble Space Telescope’s (HST) Near Infrared Camera and Multi-Object Spectrometer 2 (NICMOS2) for a small subset of the SINS objects provided additional rest-frame optical morphologies, in agreement with the disc or merger nature from the H kinematics (Förster Schreiber et al., 2011a, b). The selection of these samples was primarily based on integrated photometry or spectroscopic properties, not morphologies, and only a modest fraction of objects (1/3) were kinematically inferred to be (major) mergers. Roughly comparable fractions were found in other sizeable NIR IFU samples at , including, e.g., those by Law et al. (2007, 2009); Gnerucci et al. (2011), and Wisnioski et al. (2011b). Ê
With the recent 3D-HST slitless grism survey (see Brammer et al., 2012, and Section 2), much larger samples of objects with NIR 3D emission line spectroscopy have become available, making it possible to address the spatial extent of star formation for extensive samples of galaxy mergers at the peak of cosmic star formation rate density. The IFU samples mentioned above have significantly higher spectral resolution, enabling detailed kinematic studies, but AO-assisted IFU observations, which provide angular resolution comparable to HST in the NIR, remain observationallyÊ expensive and suffer from complications due to strong night sky lines. HST grism observations, as the ones taken in the 3D-HST survey, provide more limited kinematic information but allow for unbiased target selection and are much more efficient at detecting and mapping the continuum and line-emission at high angular resolution for all targets within the field-of-view. 3D-HST provides resolved line-emission, enabling studies of the spatial extent of star formation for large samples of galaxies at in five well-studied cosmological fields. The initial papers use approximately half of the full dataset, as described in van Dokkum et al. (2011) and Brammer et al. (2012).
Using the same 3D-HST data, we explore the spatial extent and distribution of star formation in a ’population snapshot’ of presumably merging systems at . For this sample we make 0.2 resolution maps of emission lines (H and [Oiii]), which trace the spatial extent of the (unobscured) star formation in these mergers and allow us to study their star formation properties in a statistical and unbiased way. This is done under the assumption that the H (for ) and [Oiii] (for ) emission of the systems trace the star formation. This has been shown to be a fair assumption for both H (Kennicutt, 1983; Gallagher et al., 1984; Kennicutt et al., 1994; Kennicutt, 1998a, b) and [Oiii] (Kennicutt, 1992; Teplitz et al., 2000; Hippelein et al., 2003), even though using [Oiii] as a quantitative indicator of star formation rate (which is not what we aim to do here) includes several complicating factors (Teplitz et al., 2000).
To help interpret our observations in a theoretical context, we create a sample of pseudo-observations from state-of-the-art smoothed-particle hydrodynamic (SPH) simulations of individual mergers, which we compare to the 3D-HST data. The simulations predict a merger sequence and star formation picture with centrally-concentrated triggered starbursts at final coalescence, enhanced star formation in tidal features, and black hole growth and accretion. The goal is to understand the observational results from 3D-HST by making a direct comparison with the predictions from the simulations. These comparisons will help explore which parameters, e.g., viewing angle, merger phase, gas fraction, mass ratio, etc. play a crucial role in determining, for example, star formation rates from observations.
In Section 2 we describe the 3D-HST survey from which the merger sample was selected. We then describe the selection of our sample in Section 3 and the procedure used to map the spatial extent of the star formation in Section 4. In Section 5 we split the sample into four different morphological types and find that most mergers exhibit star formation in only one component. In Section 6 we compare numerical merger simulations to the observed 3D-HST spectra and find that in simulations star formation most commonly occurs in both components, before we summarise and conclude in Section 7.
2 The 3D-HST Survey Data
To construct the resolution emission line maps, the proxies for star formation maps, we take advantage of the NIR 3D spectroscopy survey possibilities that the Wide Field Camera 3 (WFC3) on HST brings. The 3D-HST survey is a 248 orbit NIR spectroscopic Hubble treasury program (Cycles 18 and 19, PI van Dokkum). It provides NIR imaging with the F140W filter and grism spectroscopy with the G141 grism over well-studied extragalactic survey fields (AEGIS, COSMOS, GOODS-S, GOODS-N, and UKIDSS/UDS). The grism spectroscopy is slitless so both spatial and spectroscopic information is available for every single object in the survey fields. The 3D-HST survey provides rest-frame optical spectra for a sample of 7000 galaxies at (van Dokkum et al., 2011; Brammer et al., 2012).
As of August 2011 the survey had observed 68 pointings over the GOODS-S, GOODS-N111The GOODS-N data were taken as part of the HST program GO-11600 (PI B. Weiner)., COSMOS and AEGIS fields (van Dokkum et al., 2011). The present work is based on 30 of these 68 pointings, where the extensive ancillary data available enables robust spectral energy distribution (SED) modelling necessary for our sample selection as described in Section 3.
2.1 The 3D-HST Grism Spectroscopy
The WFC3 G141 grism used in 3D-HST disperses the light over the wavelength range from 1.05 m to 1.7 m with a low spectral resolution of .
Since the grism spectroscopy is slitless, the WFC3 G141 grism basically produces an emission line image that is superimposed onto a sequence of dispersed monochromatic continuum images, and some of the key features of slitless spectroscopy therefore need to be taken into account. First of all, the width of emission/absorption lines in the dispersion direction in slitless spectroscopy is not only caused by velocity broadening (which is negligible for the low resolution 3D-HST spectra) and the intrinsic broadening of the wavelength dispersion: as slitless spectroscopy produces shifted monochromatic images, the spatial extent of the dispersed emission line image reflects the spatial distribution of the line emission both along and perpendicular to the dispersion direction. We will take advantage of this ’morphology-broadening’ (which can be seen in the third panel of Figure 4), to map the spatial extent of star formation as described in Section 4. Figure 4 will be explained in more detail here.
As with multi-slit spectroscopy, the differing wavelength coverage of the spectra is an issue. Since the detector onto which the field-of-view is dispersed has a finite size, approximately 10% of the spectra are cut off on the edge of the detector.
Lastly, ‘contamination’ is an important property of slitless spectroscopy. Since the focal plane is not blocked out with a slit or a mask as is usually done in standard spectroscopy, all the light from a given object, and all other objects in the observed field, are dispersed onto the detector. Hence, spectra will often overlap and therefore ‘contaminate’ each other as explained in Brammer et al. (2012).
For more information on the data reduction methods, the data products of the 3D-HST survey, and the survey itself, we refer to Brammer et al. (2012).
3 Selecting Merger Candidates
We select our sample of merger candidates from the first 68 pointings obtained as part of the 3D-HST survey based on three different inputs: (i) The 3D-HST survey catalogue, (ii) SED modelling, and most importantly (iii) visual inspection of NIR (F140W) morphologies. The first two selection steps are performed to define an initial sample of systems with sufficient spectral coverage and to minimise the number of objects to visually inspect. We will describe each of these three steps below.
To obtain robust star formation rate estimates via the SED fitting described in Section 3.2, we require extensive ancillary photometric catalogues. We therefore focus on the 30 pointings of data available in GOODS-S (6) and COSMOS (24), where the photometric data in the FIREWORKS (Wuyts et al., 2008) and the NEWFIRM medium band survey (NMBS; Whitaker et al., 2011) catalogues are available, respectively. Hence, this work is performed on approximately 1/5 of the final 3D-HST data product.
3.1 Grism Catalogue Cuts
The first step in defining our merger sample is to select a well-defined sample of objects based on the data products of the 3D-HST survey. We ensure that at least 75% of each spectrum in the sample falls on the detector. Since we are looking for merging objects we do not put any constraints on the contamination of the individual spectra, as spectra of close pairs will always have a high level of contamination. We rely on the visual inspection (Section 3.3) to remove cases with heavy contamination from objects that are not part of the potentially merging system.
Each individual object in the 3D-HST catalogue has been matched to the available ancillary photometric catalogues. Since the 3D-HST catalogue is selected from the deep (; Brammer et al., 2012) high-resolution NIR HST F140W images, and the photometric catalogues are ground-based and shallower, not all 3D-HST objects can be matched to an object in the photometric catalogues. We only selected objects with a counterpart (within ) in ground-based photometric catalogues. Faint objects have the risk of being assigned to a bright(er) counterpart’s photometry, have low signal-to-noise (S/N), and less reliable redshifts, and we therefore restrict ourselves to objects with .
Last but not least, the 3D-HST catalogues provide a redshift estimate for the objects based on the extracted grism spectra. The catalogue grism redshifts, , are obtained by collapsing the 2D grism spectrum into a 1D spectrum, combining it with available photometry, and then estimating the redshift with an updated version of the EAZY code (Brammer et al., 2008). The redshift range where H and/or [Oiii] emission fall in the G141 grism wavelength range is (see Figure 1 in Brammer et al., 2012). We are interested in tracing the star formation in the merging systems via either H or [Oiii] emission, and therefore use to select objects in this particular redshift range.
Applying these five initial cuts (listed in the top half of Table 1) reduces the full sample of 21460 detections in the 30 GOODS-S and COSMOS pointings to 1542 objects.
3.2 Fitting SEDs to Photometry
We select star-forming systems that are expected to have significant emission line features based on their star formation rates (SFR), specific SFR (sSFR), and stellar masses (). We obtain the SFR, sSFR, and of each individual object from modelling the SED based on the ancillary photometric catalogues using a Chabrier (2003) IMF with the code MAGPHYS presented in da Cunha et al. (2008). We use the 37 NMBS bands (Whitaker et al., 2011) for the COSMOS objects and the 17 FIREWORKS bands (Wuyts et al., 2008) for the GOODS-S objects. Both catalogues span from the far-UV to MIPS 24 m. The 3D-HST catalogue redshift is used as a prior when fitting the ancillary photometric data for each object.
Although the photometric measurements are a blend of two or more components in many cases, i.e., only one photometric ID corresponds to each merger whereas several components are clearly distinguishable in the high-resolution HST imaging, selecting the high-SFR objects based on SED fits to the photometry is still very effective in selecting objects with strong emission line features.
As shown in the bottom part of Table 1 we select objects with SFR , sSFR and . A SFR of roughly corresponds to an emission line flux erg/s/cm at which corresponds to a (collapsed) emission line S/N of 8 at erg/s/cm (Brammer et al., 2012). Hence, concentrated star formation, i.e., emitted from a modest amount of pixels of this order, should be well-detected in the 3D-HST spectra. On the other hand, if the total emission line flux is spread over a larger area, the S/N per pixel might become too low for clear detection (see Section 5).
We find 352 of the initial 1542 objects in the range that satisfy these SED criteria.
3.3 Visual Inspection of NIR Morphology
Previous studies have argued both for visual (Robaina et al., 2009) and algorithmic merger identification (e.g., Lotz et al., 2008a; Conselice et al., 2008). For this pilot study, which is focussed on the emission line morphology not on the merger rates, we have, as discussed below, decided to use a visual classification. The visual inspection of the remaining 352 objects is the crucial final step in the merger sample selection process.
The visual inspection is based on the NIR F140W morphology of the objects from the 3D-HST direct imaging. The NIR images show the rest-frame optical emission at the redshifts of our galaxies. We assume that the observed NIR (i.e., rest-frame optical) morphology traces the distribution of the (intermediate-age) stellar component of the galaxies, and that it is therefore different from the stars being formed (current star formation) as traced by the emission lines. A caveat to such an assumption is that if a galaxy does not have a dominating intermediate-age stellar population and has a high SFR, the morphology in the rest-frame optical will to some extent reflect the distribution of young stars as well. The criteria used to select the merger candidates from the 352 objects are that (i) they should show a morphology that differs from the bulk of the ‘normal’ isolated galaxies, i.e., a disturbed irregular/asymmetric morphology, and (ii) they have to show several distinct components in the continuum image, either multiple objects within the arcsec F140W thumbnails or pronounced tidal features extending from the main continuum emission component. It should be noted that because our merger selection is based only on this visual classification, and the low resolution of the G141 grism () does not offer any kinematic information of the individual companions, it is impossible to address whether or not the systems are gravitationally bound. Our merger sample therefore consists of potentially merging systems.
A fraction of the more widely separated merger pairs could therefore be potential low- or high-redshift interlopers which would artificially enhance the number of mergers found. The distances between the majority of the individual merger components in the candidates selected here are of the order , corresponding to roughly 8.5 kpc at . Law et al. (2012) estimates that % of galaxies have projected false pairs within 16 kpc. This serves as an upper limit on the expected false pair fraction for our sample. As we also include late-stage mergers and not only widely separated pairs, based on this a more realistic estimate of the amount of interlopers would probably be %.
Furthermore, several studies have shown that the morphology of isolated star-forming galaxies at higher redshifts is often clumpy and irregular (Bournaud et al., 2007; Genzel et al., 2008, 2011; Law et al., 2009, 2012; Elmegreen et al., 2009; Kriek et al., 2009; Förster Schreiber et al., 2011a, b; Wuyts et al., 2012; Wisnioski et al., 2011a). This potentially biases our visual classification as inclusion of such systems will artificially enhance the number of selected mergers and hence the estimated fraction of mergers in our sample. As described below both parametric and visual classification schemes will be subject to this bias. Hence, what is described as a merger in the present selection might in fact be a galaxy with a clumpy and irregular light distribution appearing like a merger remnant. This caveat should be kept in mind when evaluating the merger candidates presented here and elsewhere. In the remainder of the paper we will therefore use the shorthand merger for likely merger candidates.
The visual inspection is not only used to select morphologically disturbed systems. Inspecting the full grism spectra, as well as the one-dimensional collapsed grism spectra, of the individual objects, we are able to discard objects without any emission line features. Without emission line features creating a star formation map as described in Section 4 is impossible. Assuming that the estimated SFR, mass, and redshift are correct, the fact that these objects lack emission line features make them very interesting in themselves, as they might be highly dust-obscured systems blocking the star formation emission. However, these objects are ignored for the present study. If an emission line feature on the other hand is observed, i.e., the (collapsed) emission line flux is roughly larger than erg/s/cm (see Section 3.2), the object is included in our sample and we attempt to create an emission line map. However, as noted in Section 3.2 and as we will see in Section 5 this does not necessarily mean that the S/N per pixel in the full grism spectrum is good enough to create an emission line map.
Visually inspecting the grism spectra also ensures that any strong contamination is due to the different components of the merging system and not due to interlopers.
From the parent sample of 352 catalogue-selected objects, the visual inspection discards 292 objects as they appear to be isolated ‘normal’ galaxies (252/292), or having high contamination not stemming from the merging components (13/292), or having no obvious emission line features in their spectra (24/292 corresponding to % of the 352 catalogue-selected objects). Hence, our final sample of (potential) mergers from 6 GOODS-S and 24 COSMOS 3D-HST pointings consists of 60 systems, corresponding to % of the 352 catalogue-selected objects, or approximately 10-13% if we correct for the expected fraction of false pair interlopes. Here the confidence intervals are the 95% quantiles of the beta distribution following Cameron (2011). This merger candidate fraction is on the high side compared to what is generally found in the literature (see, e.g., Figure 1 in Lotz et al., 2011 and Williams et al., 2011), again suggesting that a fraction of the visually identified mergers may be single objects with a blotchy distribution of young stars which increases the apparent merger fraction. The obtained merger fraction is subjective as it relies on a visual assessment of disturbance. As described above, also potential galaxy interlopers which have passed the subjective visual inspection could be biasing our sample towards higher merger fractions. Williams et al. (2011) used a mass-selected sample () of galaxy pairs when estimating a major merger fraction of 5%. This could also account for part of the discrepancy, as we are not only looking at distinct pairs of galaxies and use a flux limit as opposed to a mass limit in the selection process. With an average merger fraction just below 10% at the merger fractions summarised in Figure 1 of Lotz et al. (2011) which includes samples selected on both mass and luminosity cuts is somewhat closer to the % (10-13% if corrected for interlopes) reported here.
The selection cuts in Table 1 positions all 60 3D-HST mergers in the blue cloud of star-forming galaxies (Bell et al., 2004; Strateva et al., 2001) in colour magnitude diagrams, which implies (as expected) that the selected mergers are mainly gas-rich ‘wet’ mergers. We note however, that this may be very different for mass-selected samples where no limits have been imposed on the overall SFR of the objects.
In Figure 1 we show the magnitude distribution of the 60 3D-HST merger candidates. In Figure 2 and 3 we plot them together with the selection regions from Table 1 (grey shaded regions) as large solid symbols. The small grey points represent the 292 objects discarded by the visual classification, i.e., the general galaxy populations satisfying the selection cuts in Table 1. In Figure 3 the distributions of SFR, sSFR, , and for the 60 merger candidates are shown as histograms on the axes of the scatter plots. In both the histograms in Figure 1 and in Figure 3, the dotted lines correspond to the 16th and 84th percentiles of the distributions, and the dashed lines show the median values. In Section 6 we will use these values to determine the parameter space to sample when simulating 3D-HST grism spectra.
Visual inspection is a more subjective way of selecting mergers than for instance empirically established parametric merger classification schemes like for instance the Gini , and selections (e.g., Lotz et al., 2004, 2006, 2008a, 2008b; Conselice, 2003; Conselice et al., 2008, 2009; Papovich et al., 2005; Scarlata et al., 2007). Both parametric and visual classification schemes have advantages and disadvantages. Determining a merger population based on parameters ensures that the selection is done in a consistent and uniform way for all objects. However, at redshifts where the NIR morphology is not fully understood and where galaxies look clumpy and irregular as described above, and where disturbed morphologies are prominent without necessarily being part of a recent merger, a parametric classification scheme might fall short of a visual classification. On the other hand visual classification is potentially biased by the subjectivity of the classifier. Nevertheless, the human eye is known to be excellent at detecting and distinguishing features in noisy images and spectra and arguably minimises the bias of the clumpy irregular morphology of high- systems.
To ease comparison with the extensive literature using parametric merger classifications we have estimated the , , , , and morphological parameters for the 60 visually selected merger candidates and the 292 visually discarded objects as shown in Appendix A. The selected candidates partially satisfy the empirical parametric merger selection but in general seems to be an average subset of the parent distribution, i.e., not distinguishing itself clearly from it. Assuming that the 60 merger candidates are reliable therefore speaks in favour of using visual classification when the NIR morphology is complicated, as a parametric classification would not be able to clearly distinguish the 60 merger candidates from the parent population.
Despite the limitations of a visual selection of mergers and the only partial agreement with the empirical parametric classifications (Appendix A), we believe that the advantages of the visual classification scheme described above outweighs a ‘blind’ parametric selection for a study like the one presented here. Furthermore, we are probing an unexplored regime (NIR at high redshift) and the current parametric methods might not be appropriately calibrated and tested here and we have therefore chosen to use the visual selection in the reminder of this paper.
4 Emission Line Mapping
To quantify the extent of (unobscured) star formation in the 60 3D-HST mergers described in the previous section, we rely on the spatial information of the H and [Oiii] emission lines that the slitless grism spectroscopy provides. From the grism spectra we create emission line maps by subtracting a model of the continuum light in the grism spectra such that only the probed emission line feature is remaining. This can then be mapped back onto the NIR continuum light-distribution of the object. In practice we:
Create a 2D continuum model for the grism spectrum.
Subtract this continuum model from the 2D grism spectrum itself.
Cross-correlate the F140W thumbnail with the continuum subtracted spectrum to map the (cut-out) continuum-subtracted thumbnail back onto the F140W thumbnail.
Each of these steps is described in detail in the following subsections and are illustrated in Figure 4.
4.1 Continuum Modelling and Subtraction
The 2D continuum models we subtract from the 3D-HST spectra are based on a third order polynomial fit to a one-dimensional spectrum. The 1D spectra are obtained from a weighted sum of the individual lines in the 2D spectra. The polynomial fit to the 1D spectrum is turned back into a 2D continuum model by concatenating rows with the 1D polynomial form weighted by a ‘slit-profile’ obtained from the columns blue-ward and red-ward of the probed emission line feature in the full 2D spectrum. Subtracting this model from the 3D-HST spectrum returns a two-dimensional emission line map as illustrated in Figure 4, where all that is left is the emission line feature. This approach is similar to the one used in Nelson et al. (2012). Often when dealing with slitless spectroscopy, and the 3D-HST grism spectra in particular, the goal is to remove contamination in a systematic manner. However, we relied on visual inspection to remove badly contaminated objects, since mergers per definition are contaminated. Remaining contamination not affecting the continuum flux of the merger was masked out when modelling the continuum before subtraction.
4.2 Constructing Emission Line Maps
We used a simple cross-correlation between the NIR F140W thumbnail continuum image of each object and the corresponding full 2D emission line map in order to map the emission line map back onto the NIR image. In practice we calculate
for each of the first columns in the full 2D emission line map, where is the number of columns in the 2D emission line map and is the width of the F140W thumbnail. The and is the flux in the pixel for the F140W thumbnail and 2D emission line map cut-out (indicated by the red box in Figure 4), respectively.
The maximum of the cross-correlation function, CC, indicates the wavelength, , where there is the largest overlap between the NIR light distribution of the object and the th cut-out of the full emission line map. The corresponds to the th column of the 2D emission line map plus .
In Figure 5 we show a collection of emission line maps (red contours) from our 3D-HST merger sample. The individual maps correspond to the region at that has been mapped back onto the NIR continuum image as illustrated in Figure 4. In each map we have marked the relative location of [Nii] 6548,6583, H 6563, [Sii] 6716,6730 and H 4861, [Oiii] 4959,5007 for the H and [Oiii] maps, respectively. In some cases the [Oiii] doublet is marginally resolved as seen in the upper left emission line map in Figure 5. The H-[Nii] composite is however not resolved in the 3D-HST grism resolution. We do not attempt to de-convolve these emission lines when creating the emission line maps. Assuming that the redshift uncertainty imposed by ignoring the [Oiii] doublet is only 0.0024. This is less than the quoted average 3D-HST redshift precision of 0.0034 (Brammer et al., 2012) and is therefore not affecting the conclusions of this study.
5 Results: The Spatial Extent of Star Formation in High- Mergers
As noted in the introduction, star formation and mergers are important parts of understanding how high- galaxies evolved into the galaxies we observe in the low- Universe. In the previous sections we have described how we select the mergers from the 3D-HST data, and subtract the continuum light in the spectra to create the emission line (star formation) maps. In this section we characterise the morphological type of the spatial extent of the star formation in the 60 3D-HST mergers. However, we first divide our sample into two different kinds of mergers: the ‘pre-mergers’ and the ‘post-mergers’. By pre-mergers we mean objects that show multiple clearly-distinct and pronounced continuum peaks in the NIR images, i.e., the optical continuum emission comes from multiple objects in the process of merging or about to merge. Examples of those are shown to the left in the two top panels of Figure 5. The post-mergers, on the other hand, are systems that have (presumably) undergone merging and now appear to be dominated by a nuclear feature in the continuum with pronounced tidal features surrounding it. Examples of these systems are shown to the right in the two top panels in Figure 5. Dividing the 60 mergers into these two sub-samples return 32 pre-mergers and 28 post-mergers. We note that the distinction between pre- and post-mergers is (operatively) our distinction to guide-the-eye when classifying and inspecting the emission line maps of the individual objects as described below. The identification with the merger phases is plausible, but will need more modelling.
To characterise the location and the spatial extent of the star formation in the 3D-HST mergers we categorise the star formation maps into the following four morphological types (SF type):
One Component: The star formation in the primary emission line feature is significantly stronger than any secondary emission line feature. The threshold used is , where and is the estimated aperture flux of the primary and secondary emission line feature, respectively. For the pre-merging systems the primary emission line feature corresponds to one of the multiple objects and for the post-mergers it refers to either the nuclear region or a tidal feature.
Both (All) Components: The star formation is pronounced/detected in all (or the majority if more than two) components of the system, i.e., .
In-between: The mapped star formation appears to be emerging from in-between the merging components. None of the post-mergers show this feature, so indeed this means in-between clearly distinguishable objects.
Low S/N: The S/N per pixel of the emission line features in the 2D grism spectrum is too low to produce a convincing emission line map. This may be the case for extended star formation as pointed out in Section 3.2. It can also be due to very dust-enshrouded star formation making the emission lines very weak.
Each of the four rows of star formation maps in Figure 5 show examples of these four SF types. The results from characterising the two classes of mergers with these four SF types are shown in Figure 6. The error-bars are obtained by bootstrapping the results, i.e., by randomly drawing 60 SF types from the results 1000 times and then using the 2 width of the resulting SF type distributions as error-bars (hence no error-bar on the post-merger SF type 3 in Figure 6, as none were found). For both the pre- and post-mergers, the star formation is most prominent in just one of the components (SF type 1) for roughly 3/5 of the objects. In roughly 1/3 of the objects, star formation was detected in all components (SF type 2). Hence, the distribution of the spatial extent of star formation among the pre- and post-mergers is consistent. In Table 2 we have listed the fractions for all 60 3D-HST mergers resulting from the classification of the emission line maps.
The difference in the rates of objects with prominent star formation in just one component (SF type 1) and and mergers with star formation of type 2 might be a consequence of dust obscuration. We are only able to probe the unobscured star formation, so in cases where one component (or the tidal feature) is much more dust-obscured than the other we would end up with star formation maps of type 1. The discrepancy between the rates could also be due to different SFRs in the different components. Since the mergers are selected based on morphology and we do not have any kinematic information, the fraction of SF type 1 objects might be biased by chance superpositions of objects on the sky at different redshifts, such that we only see line emission from one object in the NIR. As described in Section 3.3 we expect approximately 4-7% of such objects due to the modest distances of 8.5 kpc involved.
The results indicate that most mergers happen between objects of different gas fractions and/or different SFR, i.e., two merging components with significantly different properties. We will show below that this is backed up by an initial comparison with simulated mergers.
In Figure 3 the different SF types are represented by different symbols to look for dependencies between the morphology of the star formation maps and SFR, sSFR, , and . These quantities are obtained on the total photometry of the merging components as only a single photometric ID was assigned to each merger and hence includes the flux of the total system. The absence of correlations seen in Figure 3 suggests that all star formation morphologies occur at all redshifts irrespective of SFR and mass.
|1) One Comp.||58%||(35/60)||28%||(83/296)|
|2) Both (All) Comp.||32%||(19/60)||59%||(175/296)|
|3) In-Between Comp.||3%||(2/60)||0%||(0/296)|
|4) Low S/N per pixel||7%||(4/60)||13%||(38/296)|
|Uncertainties are obtained via bootstrapping as described in the text.|
|See Figure 8 for a plot of these values.|
6 Simulating 3D-HST Spectra
With the exceptional data of the 3D-HST merger sample presented above, we can perform comparisons with the star formation produced in simulated mergers at high redshift. Current high-resolution merger simulations are able to predict the spatial distribution of star formation in mergers and produce simulated images by including the emission by young, newly-formed stars and the transfer of starlight through gas and dust (Cox et al., 2006, 2008; Jonsson et al., 2010). Lotz et al. (2008a, b, 2010a, 2010b) used these and similar simulations to describe the correlations between the rest-frame optical morphology of mergers, their mass ratios, total star formation, projected size, gas fractions, and merger time-scales, enabling comparisons with observations.
Simulations have shown examples of star formation triggered by direct galaxy interaction, star formation originating in the central cores of the individual merging components, and star formation appearing in tidal features. In the nearby Universe, many examples of such features have been found, e.g., in the Antennae Galaxies (Wang et al., 2004). Whether the predictions of simulations are also representative of the star formation in mergers at , i.e., when the Universe was only 3-4 Gyr old, have not been tested yet. With 3D-HST, the sample size of data with both rest-frame optical morphology and star formation morphology is becoming large enough that we can start looking at a population (snapshot) of mergers instead of individual case studies, and hence comparisons with predictions from simulations become feasible. In this section we make an initial illustrative attempt at comparing state-of-the-art merger simulations with the observational results from 3D-HST. In the following we describe how we create simulated WFC3 G141 grism spectra from a small initial set of high resolution simulated mergers.
6.1 SPH Merger Simulations
Knowledge about both the spatial and spectral extent of the objects is crucial for simulating grism spectra. As input for our grism simulations we use three-dimensional data cubes of simulated mergers. These input data cubes are based on the simulations described in Younger et al. (2008) and are designed in the spirit of the Cox et al. (2006, 2008) simulations. The simulations are -body/SPH simulations from the GADGET code (Springel et al., 2001) from which simulated observations are generated with the radiative transfer code SUNRISE (Jonsson, 2006; Jonsson et al., 2010). For our initial comparison of simulations with actual data, we use the three models listed in Table 3. The three simulations are all of mass ratio 2:1 and the two merging objects in each simulation both have a gas fraction () of 40%, making them ‘wet’ mergers like the 3D-HST systems and mimicking the suggested gas fraction of galaxies (Daddi et al., 2010; Tacconi et al., 2010). The three mergers happen at three different relative orbital inclinations () of the two merging components. Hence, what we test here is the effect of orientation on the detected star formation distribution. We further have a series of time steps or snapshots of each galaxy () and a set of viewing angles () that we can simulate spectra for. Thus, the simulation parameter space sampled here is spanned by time, orientation of the two mergers, and viewing angle. As described below, we combine this with a set of parameters determined by the 3D-HST sample to define the total parameter space to be covered by the simulated spectra.
|Sim. 1||6||3||(Sb) 2:1 (Sc)||0.4:0.4|
|Sim. 2||5||3||(Sb) 2:1 (Sc)||0.4:0.4|
|Sim. 3||3||3||(Sb) 2:1 (Sc)||0.4:0.4|
|Parenthesis give disc model (Younger et al., 2008).|
6.2 Simulating 3D-HST grism Spectra
Simulating the grism spectra from the merger simulation output data cubes presented above is fairly straightforward. The mock grism spectra are created by dividing the data cube into wavelength slices, i.e., images corresponding to a ‘filter’ of width . Offsetting or dispersing this sequence of images and co-adding the fluxes results in a grism spectrum. Using this approach, we turn the input data cubes into a sequence of simulated WFC3 G141 grism spectra. We use a pixel scale of and a spectral resolution of according to the 3D-HST grism spectra characteristics described in Brammer et al. (2012). For the HST PSF we use a Tiny Tim PSF222http://www.stsci.edu/hst/observatory/focus/TinyTim and for the system throughput we used the G141 sensitivity curve.
We use the parameter space of the observed merger sample from 3D-HST to define the parameter space of the simulations. By re-scaling each data cube we simulate grism spectra corresponding to the 16th and 84th percentiles of the magnitude (), redshift (), and SFR distributions of the 3D-HST data as indicated by the dashed lines in the corresponding histograms in Figure 1 and 3.
To adjust the SFR of the data cubes before turning them into grism spectra, we assume that the SFR to H-flux conversion follows the empirical relation of Kennicutt et al. (1994):
This relation assumes a Salpeter (1955) initial mass function and is applied to the emerging SFR, i.e., after the effects of extinction are taken into account by SUNRISE. The is obtained by scaling the flux (including emission lines) integrated over the F140W passband independently. Combining the values from Table 3 with these data-determined parameters, we end up with a sample of 336 simulated spectra spanning the parameter space .
Added to the raw mock grism spectra are Poisson noise and noise terms corresponding to the read noise and dark current quoted on the WFC3 website, as well as the background sky levels presented in Brammer et al. (2012).
In Figure 7 we present a sequence of simulated G141 grism spectra. The spectra shown have various combinations of , , , , , and SFR. The thumbnails on the left are the noise-free simulated F140W images of the objects.
6.3 Creating Star Formation Maps of Simulated Spectra
To create the star formation maps for the simulated spectra, we treat them in exactly the same way as the actual 3D-HST spectra (see Section 4). First, we visually inspect the spectra to make sure that an emission line feature is available; for some of the faint magnitude, low-SFR combinations the emission lines do not show up in the noise-added spectra. Having eliminated spectra without emission line features, as well as cases of viewing angle and time step were the two merging objects could not be distinguished, the sample of 336 simulated spectra is reduced to 296. For these 296 spectra, we create emission line maps by the method described in Section 4.
As for the actual 3D-HST data, we characterise the 296 emission line maps of the simulated spectra according to the four SF types described in Section 5. The results from this classification are shown in Table 2 together with the results for the 3D-HST mergers. The simulated spectra in Figure 7 show examples of both spectra resulting in star formation maps of type 1 (first three spectra) and type 2 (last three spectra).
6.4 Comparing the 3D-HST Data with Simulations
We compare the 3D-HST SF types we inferred in Section 5 with the simulated SF types in order to investigate how, e.g., SFR and viewing angle influences the frequencies of different types of mergers from the characterisation of the extent of star formation in mergers. In Figure 8 the SF types of the full sample of 60 3D-HST mergers (blue circles) are shown together with the SF types of the 296 simulated spectra of the 3 mergers from Table 3 (green triangles). The plotted percentages are given in Table 2. The error-bars are again obtained via bootstrapping. The comparison needs to be done with care, as the sample of simulated mergers from Table 3 does not span the gas fraction and mass ratio dimensions of the parameter space, expected to influence the ‘observability’ of star formation (Lotz et al., 2010a, b, 2011), as they are fixed at 40% and 2:1, respectively (see Section 6.1). Nevertheless, our analysis shows how the 3D-HST survey enables a direct comparison of predictions from simulations with an extensive sample of high- mergers.
Despite the limitations of he spanned simulation space, it seems that the fraction of spectra where the S/N per pixel is too low to create an actual emission line map is fairly consistent. Likewise, the fraction of cases where the star formation seems to emerge from in-between the merging components (SF type 3) is comparable: in fact we do not find any objects in this category among the simulated spectra. Looking at the two SF type 3 cases we find for the data (third row of Figure 5) could indicate that the reason we are not seeing star formation from the components themselves is due to dust-obscuration rather than more pronounced star formation in-between the components. To find similar trends in the simulated emission line maps we would need to investigate a range of gas fractions rather than just the cases presented here. This is however beyond the scope of this initial comparison.
From Figure 8 it is also evident that the simulated spectra, on average, have twice as many star formation maps where both components have pronounced star formation (SF type 2), i.e., nearly 2/3 of the simulated maps, as compared to 1/3 for the 3D-HST star formation maps. Correspondingly, the fraction of single-component emission line maps is significantly lower for the simulations than for the observations. As the simulations assume fairly similar properties between the two merging components (same gas fraction, similar mass in the form of major mergers as opposed to minor mergers, same dust content, and total SFR), this suggests that the majority of mergers do not happen between such galaxies. According to the 3D-HST data, only 1/3 of real mergers seem to happen between galaxies of similar properties. For example, different gas fractions or a larger gap between the masses (minor mergers) could change this picture, in agreement with the results presented in Section 5.
Also, the distribution of orbits and orientations of the merger components could play a significant role here. The presented simulations consist of both prograde-prograde (Sim. 1), prograde-polar (Sim. 2) and prograde-retrograde (Sim. 3) mergers. We find a deficit of SF type 1 mergers among the Sim. 1 emission line maps, in good agreement with the notion that prograde-prograde mergers generally result in more symmetric star formation (di Matteo et al., 2007). Hence, several factors might play a role in the difference found between the 3D-HST and simulated emission line (star formation) maps.
We have presented a sample of galaxy merger candidates with three-dimensional (R.A., Dec., and ) spectroscopy at redshift . The sample consists of 60 morphologically selected mergers from the Hubble treasury slitless grism survey 3D-HST with total masses and star formation rates for the systems derived from multi-wavelength photometry. From the slitless grism spectroscopy we created emission line maps of the rest-frame optical emission lines H and/or [Oiii] as a proxy for the spatial extent of the unobscured star formation. Our results go towards a comprehensive empirical picture of where star formation happens in galaxy mergers at when the cosmic star formation and merger rate were at their peak.
We have also carried out an illustrative comparison of the 3D-HST mergers to recent SPH simulations of galaxy mergers which include star formation and dust extinction. This small sample of simulations points the way toward more extensive simulation programs aimed at spanning the full parameter space of the observations.
The main conclusions of the present study are:
The spatial distribution of star formation in candidate mergers shows a broad range of morphologies.Ê It is often concentrated in a single, compact region, but can also be located in tidal tails or in-between the main stellar bodies of the progenitors.
In the majority (58%) of the early-stage, pre-coalescence mergers, the star formation is significantly more pronounced in just one of the merging components.ÊThis is likely due to different gas fractions of the progenitors. Alternatively, dust content may differ in quantity or configuration.
The star formation morphologies show no clear correlations with the estimated SFR, sSFR, , and , suggesting that all star formation morphologies are present at all epochs irrespective of star formation rate and mass.
Simulated mergers among galaxies with similar masses and similar gas fractions typically predict similarly intense star formation in both merger components, at odds with the observations: as opposed to 32% of the observed mergers, as many as 59% of the simulated mergers show detectable star formation in both components. Together with the large fraction of early-stage, pre-coalescence mergers with star formation in just one component, this discrepancy supports the notion that mergers typically occur between galaxies with distinctly different properties such as gas fraction, mass and/or SFR.
We acknowledge funding from ERC grant HIGHZ no. 227749. This work was funded in part by the Marie Curie Initial Training Network ELIXIR of the European Commission under contract PITN-GA-2008-214227. The work was mainly done while K.B.S. was a member of the International Max Planck Research School for Astronomy and Cosmic Physics at the University of Heidelberg (IMPRS-HD), Germany. This work is based on observations taken by the 3D-HST Treasury Program with the NASA/ESA HST, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555.
Appendix A Parametrisation of merger morphology ()
To ease comparison between our visually selected merger candidates and the extensive literature selecting mergers based on empirically determined parametric classification schemes, we have estimated the Gini coefficient (), the second-order moment of the brightest 20% of the light (), the concentration (), the asymmetry parameter (), and the smoothness parameter () (e.g., Lotz et al., 2004, 2006, 2008b; Conselice, 2003; Conselice et al., 2008, 2009; Papovich et al., 2005; Scarlata et al., 2007) for the objects presented in this paper. We follow Lotz et al. (2008b) when estimating the individual parameters.
In Figure 9 the estimated parameters are shown for the 60 visually selected 3D-HST merger candidates (large black points) and the 292 visually discarded objects (small grey points). The histograms show the distributions of the individual parameters of the 60 merger candidates. The solid green and dashed red lines are from Lotz et al. (2008a, b) and Conselice et al. (2009), respectively. Parametrically selected mergers have been shown to preferentially lie above and to the right of these lines. We see that the majority of the visually selected 3D-HST merger candidates indeed have as expected. The – selection on the other hand seems to reject the majority of the merger candidates selected here. The selected merger candidates generally trace the parent population and does not distinct themselves clearly as would be expected for a ‘clean’ merger selection. This is most probably due to the complications that clumpy and irregular morphology of high- objects and noisy images have on merger selections and in particular on ‘blind’ parametric classification schemes.
Appendix B The 3D-HST Emission Line (Star Formation) Maps
In Figure 10 we show the full sample of the 60 3D-HST emission line (star formation) maps. The objects have been sorted according to the estimated morphology of the star formation (SF type) described in Section 5.
- Barnes & Hernquist (1996) Barnes J. E., Hernquist L., 1996, Astrophysical Journal v.471, 471, 115
- Barton et al. (2007) Barton E. J., Arnold J. A., Zentner A. R., Bullock J. S., Wechsler R. H., 2007, The Astrophysical Journal, 671, 1538
- Barton et al. (2000) Barton E. J., Geller M. J., Kenyon S. J., 2000, The Astrophysical Journal, 530, 660
- Bell et al. (2012) Bell E. F. et al., 2012, The Astrophysical Journal, 753, 167
- Bell et al. (2004) Bell E. F. et al., 2004, The Astrophysical Journal, 608, 752
- Bigiel et al. (2011) Bigiel F., Leroy A., Walter F., 2011, Computational Star Formation, 270, 327
- Bigiel et al. (2008) Bigiel F., Leroy A., Walter F., Brinks E., de Blok W. J. G., Madore B., Thornley M. D., 2008, The Astronomical Journal, 136, 2846
- Bluck et al. (2012) Bluck A. F. L., Conselice C. J., Buitrago F., Grützbauch R., Hoyos C., Mortlock A., Bauer A. E., 2012, The Astrophysical Journal, 747, 34
- Bournaud et al. (2007) Bournaud F., Elmegreen B. G., Elmegreen D. M., 2007, The Astrophysical Journal, 670, 237
- Brammer et al. (2008) Brammer G. B., van Dokkum P. G., Coppi P., 2008, The Astrophysical Journal, 686, 1503
- Brammer et al. (2012) Brammer G. B., van Dokkum P. G., Schmidt K. B., Nelson E. J., 2012, Astrophysical Journal
- Calzetti et al. (2012) Calzetti D., Liu G., Koda J., 2012, eprint arXiv, 1204, 5659
- Cameron (2011) Cameron E., 2011, Publications of the Astronomical Society of Australia, 28, 128, (c) 2011 Astronomical Society of Australia
- Chabrier (2003) Chabrier G., 2003, The Publications of the Astronomical Society of the Pacific, 115, 763
- Conselice (2003) Conselice C. J., 2003, The Astrophysical Journal Supplement Series, 147, 1
- Conselice et al. (2011) Conselice C. J., Bluck A. F. L., Ravindranath S., Mortlock A., Koekemoer A. M., Buitrago F., Grützbauch R., Penny S. J., 2011, Monthly Notices RAS, 417, 2770
- Conselice et al. (2008) Conselice C. J., Rajgor S., Myers R., 2008, Monthly Notices RAS, 386, 909
- Conselice et al. (2009) Conselice C. J., Yang C., Bluck A. F. L., 2009, Monthly Notices RAS, 394, 1956
- Contini et al. (2012) Contini T. et al., 2012, Astronomy & Astrophysics, 539, 91
- Cox (2004) Cox T. J., 2004, Ph.D dissertation, 27
- Cox et al. (2006) Cox T. J., Jonsson P., Primack J. R., Somerville R. S., 2006, Monthly Notices RAS, 373, 1013
- Cox et al. (2008) Cox T. J., Jonsson P., Somerville R. S., Primack J. R., Dekel A., 2008, Monthly Notices RAS, 384, 386, (c) Journal compilation © 2008 RAS
- Cresci et al. (2009) Cresci G. et al., 2009, The Astrophysical Journal, 697, 115
- da Cunha et al. (2008) da Cunha E., Charlot S., Elbaz D., 2008, Monthly Notices RAS, 388, 1595
- Daddi et al. (2010) Daddi E. et al., 2010, The Astrophysical Journal, 713, 686
- di Matteo et al. (2007) di Matteo P., Combes F., Melchior A.-L., Semelin B., 2007, Astronomy and Astrophysics, 468, 61
- Elmegreen et al. (2009) Elmegreen B. G., Elmegreen D. M., Fernandez M. X., Lemonias J. J., 2009, The Astrophysical Journal, 692, 12
- Epinat et al. (2009) Epinat B. et al., 2009, Astronomy and Astrophysics, 504, 789
- Epinat et al. (2012) Epinat B. et al., 2012, Astronomy & Astrophysics, 539, 92
- Förster Schreiber et al. (2009) Förster Schreiber N. M. et al., 2009, The Astrophysical Journal, 706, 1364
- Förster Schreiber et al. (2011a) Förster Schreiber N. M., Shapley A. E., Erb D. K., Genzel R., Steidel C. C., Bouché N., Cresci G., Davies R., 2011a, The Astrophysical Journal, 731, 65
- Förster Schreiber et al. (2011b) Förster Schreiber N. M. et al., 2011b, The Astrophysical Journal, 739, 45
- Gallagher et al. (1984) Gallagher J. S., Hunter D. A., Tutukov A. V., 1984, Astrophysical Journal, 284, 544
- Genzel et al. (2008) Genzel R. et al., 2008, The Astrophysical Journal, 687, 59
- Genzel et al. (2011) Genzel R. et al., 2011, The Astrophysical Journal, 733, 101
- Genzel et al. (2006) Genzel R. et al., 2006, Nature, 442, 786
- Gnerucci et al. (2011) Gnerucci A. et al., 2011, Astronomy & Astrophysics, 528, 88
- Hammer et al. (2005) Hammer F., Flores H., Elbaz D., Zheng X. Z., Liang Y. C., Cesarsky C., 2005, Astronomy and Astrophysics, 430, 115
- Hippelein et al. (2003) Hippelein H. et al., 2003, Astronomy and Astrophysics, 402, 65
- Hopkins & Beacom (2006) Hopkins A. M., Beacom J. F., 2006, The Astrophysical Journal, 651, 142
- Hopkins et al. (2010) Hopkins P. F., Younger J. D., Hayward C. C., Narayanan D., Hernquist L., 2010, Monthly Notices RAS, 402, 1693
- James et al. (2004) James P. A. et al., 2004, Astronomy and Astrophysics, 414, 23
- Jogee et al. (2009) Jogee S. et al., 2009, The Astrophysical Journal, 697, 1971
- Jonsson (2006) Jonsson P., 2006, Monthly Notices RAS, 372, 2
- Jonsson et al. (2010) Jonsson P., Groves B. A., Cox T. J., 2010, Monthly Notices RAS, 403, 17, 32 pages, 22 figures, submitted to MNRAS
- Joseph & Wright (1985) Joseph R. D., Wright G. S., 1985, Monthly Notices of the Royal Astronomical Society (ISSN 0035-8711), 214, 87
- Karim et al. (2011) Karim A. et al., 2011, The Astrophysical Journal, 730, 61
- Kennicutt (1983) Kennicutt R. C., 1983, Astrophysical Journal, 272, 54
- Kennicutt (1992) Kennicutt R. C., 1992, Astrophysical Journal, 388, 310
- Kennicutt (1998a) Kennicutt R. C., 1998a, Astrophysical Journal v.498, 498, 541
- Kennicutt (1998b) Kennicutt R. C., 1998b, Annu. Rev. Astro. Astrophys., 36, 189
- Kennicutt et al. (1994) Kennicutt R. C., Tamblyn P., Congdon C. E., 1994, Astrophysical Journal, 435, 22
- Klaas & Elsaesser (1991) Klaas U., Elsaesser H., 1991, Astronomy and Astrophysics Supplement Series (ISSN 0365-0138), 90, 33
- Kriek et al. (2009) Kriek M., van Dokkum P. G., Franx M., Illingworth G. D., Magee D. K., 2009, The Astrophysical Journal Letters, 705, L71
- Lambas et al. (2003) Lambas D. G., Tissera P. B., Alonso M. S., Coldwell G., 2003, Monthly Notices RAS, 346, 1189
- Law et al. (2007) Law D. R., Steidel C. C., Erb D. K., Larkin J. E., Pettini M., Shapley A. E., Wright S. A., 2007, The Astrophysical Journal, 669, 929
- Law et al. (2009) Law D. R., Steidel C. C., Erb D. K., Larkin J. E., Pettini M., Shapley A. E., Wright S. A., 2009, The Astrophysical Journal, 697, 2057
- Law et al. (2012) Law D. R., Steidel C. C., Shapley A. E., Nagy S. R., Reddy N. A., Erb D. K., 2012, The Astrophysical Journal, 745, 85
- Lotz et al. (2008a) Lotz J. M. et al., 2008a, The Astrophysical Journal, 672, 177
- Lotz et al. (2011) Lotz J. M., Jonsson P., Cox T. J., Croton D., Primack J. R., Somerville R. S., Stewart K., 2011, The Astrophysical Journal, 742, 103
- Lotz et al. (2008b) Lotz J. M., Jonsson P., Cox T. J., Primack J. R., 2008b, Monthly Notices RAS, 391, 1137
- Lotz et al. (2010a) Lotz J. M., Jonsson P., Cox T. J., Primack J. R., 2010a, Monthly Notices RAS, 404, 590, (c) Journal compilation © 2010 RAS
- Lotz et al. (2010b) Lotz J. M., Jonsson P., Cox T. J., Primack J. R., 2010b, Monthly Notices RAS, 404, 575, (c) Journal compilation © 2010 RAS
- Lotz et al. (2006) Lotz J. M., Madau P., Giavalisco M., Primack J., Ferguson H. C., 2006, The Astrophysical Journal, 636, 592
- Lotz et al. (2004) Lotz J. M., Primack J., Madau P., 2004, The Astronomical Journal, 128, 163
- Mancini et al. (2011) Mancini C. et al., 2011, The Astrophysical Journal, 743, 86
- Melnick & Mirabel (1990) Melnick J., Mirabel I. F., 1990, Astronomy and Astrophysics (ISSN 0004-6361), 231, L19
- Mihos & Hernquist (1996) Mihos J. C., Hernquist L., 1996, Astrophysical Journal v.464, 464, 641
- Nelson et al. (2012) Nelson E. J. et al., 2012, The Astrophysical Journal Letters, 747, L28
- Newman et al. (2012) Newman S. F. et al., 2012, eprint arXiv, 1204, 4727
- Papovich et al. (2005) Papovich C., Dickinson M., Giavalisco M., Conselice C. J., Ferguson H. C., 2005, The Astrophysical Journal, 631, 101
- Queyrel et al. (2012) Queyrel J. et al., 2012, Astronomy & Astrophysics, 539, 93
- Queyrel et al. (2009) Queyrel J. et al., 2009, Astronomy and Astrophysics, 506, 681
- Rieke et al. (1985) Rieke G. H., Cutri R. M., Black J. H., Kailey W. F., McAlary C. W., Lebofsky M. J., Elston R., 1985, Astrophysical Journal, 290, 116
- Robaina et al. (2009) Robaina A. R. et al., 2009, The Astrophysical Journal, 704, 324
- Robaina et al. (2010) Robaina A. R., Bell E. F., Wel A. V. D., Somerville R. S., Skelton R. E., McIntosh D. H., Meisenheimer K., Wolf C., 2010, The Astrophysical Journal, 719, 844
- Salpeter (1955) Salpeter E. E., 1955, Astrophysical Journal, 121, 161
- Scarlata et al. (2007) Scarlata C. et al., 2007, The Astrophysical Journal Supplement Series, 172, 406
- Schruba et al. (2011) Schruba A. et al., 2011, The Astronomical Journal, 142, 37
- Shapiro et al. (2008) Shapiro K. L. et al., 2008, The Astrophysical Journal, 682, 231
- Springel (2000) Springel V., 2000, Monthly Notices RAS, 312, 859
- Springel et al. (2001) Springel V., Yoshida N., White S. D. M., 2001, New Astronomy, 6, 79
- Strateva et al. (2001) Strateva I. et al., 2001, The Astronomical Journal, 122, 1861
- Swinbank et al. (2004) Swinbank A. M., Smail I., Chapman S. C., Blain A. W., Ivison R. J., Keel W. C., 2004, The Astrophysical Journal, 617, 64
- Tacconi et al. (2010) Tacconi L. J. et al., 2010, Nature, 463, 781
- Teplitz et al. (2000) Teplitz H. I. et al., 2000, The Astrophysical Journal, 542, 18
- van Dokkum et al. (2011) van Dokkum P. G. et al., 2011, The Astrophysical Journal Letters, 743, L15, submitted to ApJ Letters
- Wang et al. (2004) Wang Z. et al., 2004, The Astrophysical Journal Supplement Series, 154, 193
- Whitaker et al. (2011) Whitaker K. E. et al., 2011, eprint arXiv, 1105, 4609
- Williams et al. (2011) Williams R. J., Quadri R. F., Franx M., 2011, The Astrophysical Journal Letters, 738, L25
- Wisnioski et al. (2011a) Wisnioski E., Glazebrook K., Blake C., Team T. W., 2011a, eprint arXiv, 1103, 5265
- Wisnioski et al. (2011b) Wisnioski E. et al., 2011b, eprint arXiv, 1107, 3338
- Wright et al. (2009) Wright S. A., Larkin J. E., Law D. R., Steidel C. C., Shapley A. E., Erb D. K., 2009, The Astrophysical Journal, 699, 421
- Wuyts et al. (2008) Wuyts S., Labbé I., Förster Schreiber N. M., Franx M., Rudnick G., Brammer G. B., van Dokkum P. G., 2008, The Astrophysical Journal, 682, 985
- Wuyts et al. (2012) Wuyts S. et al., 2012, eprint arXiv, 1203, 2611
- Younger et al. (2008) Younger J. D., Hopkins P. F., Cox T. J., Hernquist L., 2008, The Astrophysical Journal, 686, 815