The Ngc~1614 interacting galaxy
Key Words.:Galaxies: active – Galaxies: evolution – Galaxies: individual: NGC 1614 – Galaxies: starburst – ISM: molecules
Minor mergers frequently occur between giant and gas-rich low-mass galaxies and can provide significant amounts of interstellar matter to refuel star formation and power active galactic nuclei (AGN) in the giant systems. Major starbursts and/or AGN result when fresh gas is transported and compressed in the central regions of the giant galaxy. This is the situation in the starburst minor merger NGC~1614, whose molecular medium we explore at half-arcsecond angular resolution through our observations of CO (21) emission using the Submillimeter Array (SMA). We compare our CO (21) maps with optical and Pa, Hubble Space Telescope and high angular resolution radio continuum images to study the relationships between dense molecular gas and the NGC 1614 starburst region. The most intense CO emission occurs in a partial ring with 230 pc radius around the center of NGC 1614, with an extension to the northwest into the dust lane that contains diffuse molecular gas. We resolve ten giant molecular associations (GMAs) in the ring, which has an integrated molecular mass of 8 10 M. Our interferometric observations filter out a large part of the CO (10) emission mapped at shorter spacings, indicating that most of the molecular gas is diffuse and that GMAs only exist near and within the circumnuclear ring. The molecular ring is uneven with most of the mass on the western side, which also contains GMAs extending into a pronounced tidal dust lane. The spatial and kinematic patterns in our data suggest that the northwest extension of the ring is a cosmic umbilical cord that is feeding molecular gas associated with the dust lane and tidal debris into the nuclear ring, which contains the bulk of the starburst activity. The astrophysical process for producing a ring structure for the final resting place of accreted gas in NGC 1614 is not fully understood, but the presence of numerous GMAs suggests an orbit-crowding or resonance phenomenon. There is some evidence that star formation is progressing radially outward within the ring, indicating that a self-triggering mechanism may also affect star formation processes. The net result of this merger therefore very likely increases the central concentration of stellar mass in the NGC 1614 remnant giant system.
Merger events, i.e., close interactions of galaxy pairs, have a deep impact on the evolution of galaxies through the triggering
of active galactic nuclei (AGN) and starburst activity. Studies of interacting galaxies often focus on major mergers (nearly equal
mass spirals) (e.g. Toomre & Toomre 1972; Tacconi et al. 2002; Dasyra et al. 2006; Jesseit et al. 2007; Somerville et al. 2008; Kormendy et al. 2009) and their evolution, although minor mergers most likely constitute
the bulk of all galaxy interactions (e.g. Toth & Ostriker 1992; Zaritsky 1995; Hunsberger et al. 1996). Investigating how the gas is feeding starburst and AGN
activities in these objects is therefore paramount for understanding the overall evolution of the Universe. Although there is still
no consensus what powers the emission in low ionization nuclear emission line regions (LINERs), nuclear starburst, AGN activity,
and shocks are three major candidates (e.g. Veilleux et al. 1999; Ho 1999; Terashima et al. 2000; Alonso-Herrero et al. 2000; Monreal-Ibero et al. 2006). Bar-driven inflow of gas in galaxies has
been suggested to trigger nuclear starbursts, and the feeding of AGNs (e.g. Simkin et al. 1980; Scoville et al. 1985). However, the radial inflow of gas
along the bar may be slowed down at certain radii, which are often associated with inner Lindblad resonances
(ILR, Combes 1988b; Shlosman et al. 1989).
NGC~1614 is an excellent example of a minor merger with a spectacular morphology and intense nuclear activity. The nucleus of this galaxy displays characteristics of containing an intense nuclear starburst region and may also contain some level of AGN activity. NGC 1614 (SB(s)c pec, de Vaucouleurs et al. 1991) is a luminous infrared galaxy (LIRG, Sanders et al. 2003) at a distance of 64 Mpc (de Vaucouleurs et al. 1991, for = 75 km s ; 1″ = 310 pc). Most of the starburst activity can be attributed to the interaction with a companion galaxy located in the tidal tail to the southwest of the nucleus of NGC 1614 (Neff et al. 1990; Olsson et al. 2010; Väisänen et al. 2012). A mass ratio of 5:1 3:1 has been found for this minor merger system (Väisänen et al. 2012). Its nuclear optical spectrum shows both starburst and LINER activity. A circumnuclear ring ( = 300 pc, see e.g. Fig. 1) in NGC 1614, seen in tracers such as Pa (Alonso-Herrero et al. 2001, along with an identification in imaging), radio continuum (Olsson et al. 2010) and polycyclic aromatic hydrocarbons (PAHs) (Väisänen et al. 2012), hosts a very young starburst (5–10 Myr, top-heavy initial mass function, Puxley & Brand 1999) while an older ( 10 Myr) starburst resides in the nucleus - possibly encircling an AGN. Alonso-Herrero et al. (2001) claimed that the ring is the cause of “wildfire” starburst activity propagating from the nucleus in an outward direction through cloud-cloud collisions and winds, while Olsson et al. (2010) instead suggest that the ring is the result of a resonance, possibly linked to a stellar bar.
Though NGC 1614 has intense nuclear activity, its spiral arms are devoid of molecular gas, which raises the question how the gas reached the nucleus. Numerical simulation studies of minor- or intermediate mergers show that gas brought in by the disturbing companion galaxy is generally found at large radii in the merger remnant (Bournaud et al. 2005). Gas returns to the system from tidal tails and often forms rings - polar, inclined or equatorial - that will appear as dust lanes when seen edge-on.
One dominant feature of the optical morphology of NGC 1614 is the prominent dust lane that crosses the galaxy just north of
the central activity. A considerable fraction of the CO (10) emission resides in this dust lane and is not directly involved
in the nuclear activity. The OVRO CO (10) observations (2.5″ resolution) presented by Olsson et al. (2010) were unable to spatially
resolve the molecular structure of the inner region.
Furthermore, a spatial correlation between CO emission and 1.4 GHz radio continuum is commonly observed in galaxies. This is not the case for NGC 1614, however. There is probably no star formation associated with the dust lane - possibly because gas is being funneled along it, with the resulting shear, which prevents efficient star formation. A similar CO-dust lane correlation has been found for another minor merger, the Medusa merger, for which the gas in the dust lane is suggested to be diffuse because of the gas-funneling (Aalto & Hüttemeister 2000).
Here we present a study of the CO emission toward the center of NGC 1614. In Sect. 2 we describe the observations and data reduction, and in Sect. 3 we present the results.
NGC~1614 was observed with the Submillimeter Array (SMA) in the CO (21) and CO (21)
transitions. We obtained one full track (10 hr), between 22 LST and 08 LST, in the SMA’s most extended
configuration, the Vext array configuration, with a resolution of 0.5″on 26 July 2010, which provided baselines between
68 m up to 509 m. Our map is therefore sensitive to structures smaller than 2.4″. Using the 230 GHz band
receiver in conjunction with the DR 128 correlator (DR = default resolution) mode with a bandwidth of 4 GHz provided channel
widths of 0.8125 MHz, resulting in a velocity coverage of 5305.7 km s with a channel separation of
1.1 km s. For data analysis a final channel width of 21 km s was used. The 4 GHz band
was centered at a rest frequency of 226.931 GHz (corresponding to a central velocity of
= 4778 km s) for the upper side band (USB). The system noise temperature ranged between
80 and 170 K. During the observations several sources were used as calibrators: 3C454.3
as the band pass calibrator, the nearby source 0423-013 for the phase calibration of NGC 1614, and
Neptune and Callisto as flux calibrators.
We reduced and analyzed the data using the MIR IDL and GILDAS111http://www.iram.fr/IRAMFR/GILDAS Mapping data reduction software packages. For the CO (21) line the synthesized beam resulting from natural weighting is 0.50″ 0.44″ with a position angle of 54°. For the 1 mm continuum observations the smoothed synthesized beam from natural weighting is 0.96″ 0.96″.
To provide high resolution images of the optical structure of NGC 1614, we obtained pipeline-reduced Hubble Space Telescope images taken with the Advanced Camera for Surveys in the F435W and F814W filters (see Fig. 1). These data were obtained as part of program G0-10592 (A. Evans, P. I.) and have an angular resolution of 0.05 arcsec pixel corresponding to a linear scale of 15 pc pixel.
3.1 1.3 mm continuum
The 1.3 mm continuum map (Fig. 2, upper left panel) shows centrally peaked emission with a flux density of 19 6 mJy, which agrees well with the values for the 1.3 mm continuum obtained by Wilson et al. (2008), who observed a total of 21 3 mJy was observed. Our 1.3 mm continuum map appears to be comprised of two components and has a box-like shape that is consistent with the structure Wilson et al. (2008) observed in their 1.3 mm continuum map.
3.2 CO and Co (21)
|integrated flux||mass 222We adopted a conversion factor of 3 10 cm (K km s)|
|CO (21)||CO (21)||(CO 21)||(CO 21)||(1.3 mm continuum)|
|[Jy km s]||[Jy km s]||[M]||[M]||[M]|
|Whole system||203.3 9.8||7.7||2.12 10||2.74 10||1.03 10|
|Ring||65.4 6.9||1.0||8.29 10||3.25 10||3.94 10|
|Dust lane connection||25.8 5.3||1.0||2.54 10||5.40 10||2.33 10|
|Northern extension||22.7 6.2||3.0||3.32 10||1.39 10||1.25 10|
3.2.1 Co (21) integrated intensity
The upper right panel of Fig. 2 shows the CO (21) integrated intensity map derived from
our data. The CO (21) line was not detected at 3 a level of 17.26 Jy km s in these SMA
observations, which agrees with values from Wilson et al. (2008), although they detected the line. The higher spatial resolution of
our Vext data might lead to a non-detection of lower density CO (21) gas that Wilson et al. (2008) were able to find with
their lower resolution data. The CO emission appears as a ring-like structure without emission at its center, with a
gas component, that might connect to the dust lane (hereafter umbilical cords), to the northwest and an extension of
lower surface brightness gas to the north (see Fig. 3). Both within and slightly beyond this ring
structure hot spots of CO are clearly identifiable. These are discussed in Sect. 3.4.
The highest flux density is 203 Jy km s (Table 1). A comparison to previous observations of Wilson et al. (2008) shows us that we recovered only about 30% of the total flux density. We find a clear asymmetry in the east-west gas distribution of the ring-like structure, or possibly tightly wound spiral arms. The eastern side has much less gas than the western part. The intensity ratio is I(west)/I(east) 2.7.
Using a CO-to-H conversion factor of = 3 10 cm (K km s) (Narayanan et al. 2011) results in a molecular mass of 2.1 10 M for the whole system. If we consider the umbilical cords (Fig. 3) as a separate mass entity, the ring itself has a gas mass of 8.3 10 M (Table 1).
3.2.2 Missing flux and the structure of the CO emission
In our CO (21) integrated intensity map we recover only 30% of the CO (21) flux density of the low-resolution SMA observations of Wilson et al. (2008), thus we most likely only see the most compact structures. A comparison of their lower resolution CO (21) map and the structures of the high- and low-resolution CO (10) observed by Olsson et al. (2010) shows remarkable similarities in the morphology. Hence we suggest that the bulk of the missing flux is located in the northern crossing dust lane (see also Fig. 1) - but also in structures that trace the beginning of the southern arm and emission to the east of the ring.
3.2.3 Co (21) velocity field, velocity dispersion, and position-velocity diagrams
The velocity field, with velocities ranging from 4520.42 km s to 4985.31 km s, is shown in
Fig. 2 (bottom left panel). It shows a behavior typical for solid-body rotation, for which the probable
connection to the dust lane follows the behavior of the gas in the ring.
The velocity dispersion distribution of NGC~1614 shows several peaks (Fig. 2, bottom right panel). At these peaks the velocity dispersions rise to 43-51 km s compared to the underlying material, which only shows dispersion values between 15 and 33 km s. The velocity dispersion might be biased by missing extended structures in our map. This probably does not affect the dispersion of the ring itself because it is compact.
In Fig. 4 the position velocity diagrams along the major (north to south) and minor (east to west) axes are shown. They were obtained by averaging over a 1″-wide slit positioned along the major and minor axes centered on the center of the ring. The pattern for solid-body rotation seen in the velocity field is also reflected in the major axis position velocity diagram. In the minor axis position velocity diagram we clearly see a difference between the eastern and the western part of the ring, with the majority of the gas residing in the western part of the ring.
3.2.4 Co / Co (21) ratio
To estimate the CO / CO (21) ratio we used the integrated intensity maps of CO (21) and CO (21) integrated over the same velocity range and determined the average flux densities (flux density for CO from the 3 rms noise value) within the same polygon for both transitions. The resulting lower limit for the flux density ratio is CO / CO (21) 15. An investigation of the spectra of Wilson et al. (2008) showed that the CO / CO (21) line ratio for the peak (central) flux density is 22, while the ratio for the integrated flux density (in the whole map) is 38. One can see from their data that the CO (21) flux density hardly rises going from peak flux densities to integrated intensities, while the CO (21) flux densities increase by more than a factor of 2. From this it seems that the CO / CO (21) ratio is increasing with radius in NGC 1614, a notion that should be confirmed at higher sensitivity and resolution; this would be reminiscent of the situation in the Medusa minor merger (NGC~4194, Aalto et al. 2010) and the molecular ring in NGC~6946 (Meier & Turner 2004). This is discussed again briefly in Sect. 4.1.2.
|Ring center coordinates||RA: 04:34:00.036, DEC: -08:34:45.08|
|Molecular mass||8.29 10 M|
|Dynamical mass||1.85 10 M|
|West-east intensity ratio||2.7|
|Radius||0.75″ ( 231 pc)|
|GMAs on the ring||6|
|GMAs off the ring||4|
3.3 Properties of the molecular ring
Figs. 2, 3, 5, and 6
show that we observe a ring-like structure in NGC~1614 with additional molecular gas emission to the north, the northeast, and
the west. From the inside out, the structure evolves via umbilical cords, a gas component that possibly connects the ring to
the dust lane located to the northeast, and then into a northern extension of lower surface density material (see
Fig. 3). In contrast to the observations from Wilson et al. (2008) and Olsson et al. (2010), we clearly filtered out the
more extended gas in this region. The angular resolution of our observations allowed us to resolve the width of the ring. In contrast to
studies of NGC 1614 at other wavelengths, such as radio continuum and Pa, these CO (21) observations detected no
emission at the center of the ring (Fig. 5).
To obtain a better overview of the properties of the different components in this object, we also determined the masses for components outside the ring. We assumed a radius of the ring of 0.75″ (231 pc) and a width of 0.5″. On this basis the mass for the ring is (H) = 8.3 10 M (Table 2). The umbilical cords to the northeast have a mass of (H) = 2.5 10 M and the mass of the northern extension is (H) = 3.3 10 M. We can now also calculate the dynamical mass for the ring by considering the radius and the width of the ring, as well as a rotational velocity of 187 km s, estimated from the position velocity diagram. Because the angle of the ring to the plane of the sky is unknown, an uncertainty factor of has to be considered. Taking this into account, we determined the dynamical mass of the ring to be = 1.85 10 M, which is a factor of 2 higher than the molecular mass of the ring. The viewing angle to the molecular ring seems to be face-on, which justifies a of 1. A low inclination would not change the result drastically.
|GMA||RA(2000)||DEC(2000)||size (diam.)333The source size was determined by fitting an elliptical Gaussian to the position of the respective GMA in the uv table. The values given here represent the major and minor axis values resulting from the fitting procedure. The typical (average) error for the fit to the GMA sizes is about 20%.||dispersion||mass|
|[h] [m] [s]||[°] [′] [″]||[″]||[pc]||[km s]||[km s]||[M]|
|GMAs in the ring:|
|GMA 4||04:34:00.004||–08:34:44.37||0.6 1.0||250 336||4813 7||61 6||5.70 10|
|GMA 6||04:33:59.985||–08:34:45.08||0.6 0.9||195 268||4737 5||64 4||9.74 10|
|GMA 7||04:33:59.998||–08:34:45.57||0.8 1.1||261 341||4670 7||38 4||5.97 10|
|GMA 8||04:34:00.012||–08:34:45.97||1.3 0.8||390 260||4609 6||73 5||9.36 10|
|GMA 9||04:34:00.067||–08:34:45.56||0.6 1.2||192 368||4671 6||54 5||5.98 10|
|GMA 10||04:34:00.082||–08:34:45.01||0.5 0.7||168 204||4756 7||40 5||4.96 10|
|GMAs outside the ring:|
|GMA 1||04:34:00.087||–08:34:43.45||0.6 1.1||197 326||4919 7||57 7||8.29 10|
|GMA 2||04:33:59.989||–08:34:43.38||1.1 0.5||336 157||4901 7||31 4||1.64 10|
|GMA 3||04:33:59.998||–08:34:43.97||0.8 1.1||250 336||4875 5||53 3||6.23 10|
|GMA 5||04:33:59.946||–08:34:44.37||1.8 0.6||566 180||4814 7||48 5||3.70 10|
3.4 Giant molecular associations (GMAs) in Ngc~1614
We found ten clumps in (GMA 4, GMA 6, GMA 7, GMA 8, GMA 9, GMA 10) and outside (GMA 1, GMA 2, GMA 3, GMA 5) the structure of
the molecular ring (Fig. 3). They have masses between
(H) = 1.6 10 M and (H) = 9.4 10 M
(Table 3), which places these objects in the mass range for GMAs, so-called giant molecular associations
(Vogel et al. 1988).
The velocity dispersion map (Fig. 2, bottom right panel) shows a rise in line widths at the positions of all GMAs in the ring (GMA 4, GMA 6, GMA 7, GMA 8, GMA 9). The dispersion values range between 31.1 km s for GMA 2 and 72.9 km s for GMA 8. The GMAs outside the ring have an average dispersion of 47.0 km s. The average dispersion for the GMAs in the ring is 55.0 km s.
We ran a decomposition of the CO data cube into individual structures, using the task Gaussclumps (developed initially by Stutzki & Guesten 1990) in GILDAS/Mapping. The algorithm searches iteratively for 3D Gaussians in the data cube until it reaches a specified threshold value. To derive reliable results, we set this threshold to a value of 6. Under these conditions, the data cube can be decomposed into ten different Gaussian clumps, whose properties are listed in Table 3 while the positions are indicated in Fig. 3. Four of the five brightest clumps (GMA 4, GMA 6, GMA 8, and GMA 9), which we call GMAs, are clearly associated with the ring, where the CO emission is the brightest. The one GMA that is very bright (GMA 3) but does not lie in the ring is situated at the point where the molecular gas in the ring is connected to the dust lanes via the feeding lines (umbilical cords, see Sect. 4.1.1 and 4.1.2). The residual map is dominated by fainter emission located to the north of the ring. The resulting source sizes (Table 3) range between 168 pc 204 pc and 390 pc 260 pc for GMAs outside the molecular ring and 336 pc 157 pc up to 566 pc 180 pc for GMAs in the ring. Moreover, the GMAs outside the ring are more elliptical in shape, i.e., the major-to-minor axis ratio is much higher than for the GMAs in the ring, especially for the GMAs outside the ring that are situated in the dust lanes. Assuming spherical morphologies (for comparison purposes only), the average diameter of a GMA outside the ring is 273 20 pc, whereas the GMAs in the ring show an average diameter of 259 24 pc.
4.1 Origin of the ring
4.1.1 “Feeding” of the ring: The connection to the dust lanes
Bournaud et al. (2005) found that in minor mergers gas brought in by a disturbing galaxy companion is generally found at large radii in the
merger remnant. The gas returns to the system via tidal tails and often forms rings - polar, inclined or equatorial - that will appear as
dust lanes seen edge-on. In the CO (21) gas distribution of NGC~1614 (Fig. 2) we found that the
feature labelled umbilical cords appears to be associated with a dust lane northwest of the nucleus of NGC 1614, which has
previously been identified with the CO (10) emission by
Olsson et al. (2010, see also Fig. 1 in this paper). This link to the dust lane might be the connection point
that fuels the starburst ring and is comparable with the twin-peak morphology in the rings of barred spiral galaxies (e.g.
objectNGC 1097, NGC~1365, see Sect. 4.3). However, for NGC 1614 we only found one connection point, and
a possible a stellar bar and its properties is still under debate (see for example the discussions in Chapelon et al. 1999; Olsson et al. 2010). This
connection could be due to orbit crowding of inflowing gas streams (umbilical cords, traced by the CO (10) and the dust lane
to the northwest in Fig. 1) at the branch-off of the ring, where gas encounters shocks and then migrates
to new orbits where it accumulates in the shape of the starburst ring. The gas at the western and southern edges of the ring is
consistent with an additional dust lane component that crosses the ring from the west to the southeast (see
Fig. 1). This additional component is possibly the site of a second connection between the starburst
ring and the dust lanes in this system, but this remains to be confirmed.
As observed in high-resolution CO (10) (Olsson et al. 2010), the starburst ring has the same velocity structure as the underlying galaxy, which we do not see because the lower surface density gas is filtered out in our high-resolution observations. Within the ring the rotation is higher than in the rest of the galaxy. The gas slows down relative to the galaxy’s rotation curve toward the edges of the CO (21) distribution at the connection point, where the dust lane connects to the ring in the northern part. The gas distribution in the ring may be lopsided because the whole ring is fueled/fed from one side only, which might also be the reason for the asymmetry of the ring itself.
We found that the CO (21) emission resolved in our interferometric observations is located at the unresolved center of the larger scale molecular gas reservoir (CO (10)). However, with the present data sets we were unable to verify possible streaming motions between the larger scale and smaller scale molecular gas reservoirs. Therefore we are unable to determine the inflow rate of the gas. To correct this, higher resolution CO (10) observations are needed.
With the large reservoir of available gas in the ring and the star formation rate (SFR) of 5.2 M yr determined from the CO (21) emission, taking the gas mass in the ring and its surface density into account, the star formation in the center of NGC 1614 could go on for another 1.6 10 yr unless processes/events occur that stop the gas transport via the feeding lines.
4.1.2 Formation of GMAs in the ring
In Sect. 4.1.1 we established that there is a connection between the dust lanes in NGC~1614 and the
molecular ring seen in CO 21 emission. Prominent features of the ring and the surrounding gaseous material are ten GMAs.
In the prominent minor axis dust lane north of the ring, however, we detect no GMAs - although this is the main location of the
bulk of the CO (21) flux (Wilson et al. 2008). Hence, the CO (21) emission here is either emerging from smaller molecular clouds with
a surface brightness too low to be detected by us - or alternatively the molecular gas is in the form of diffuse, unbound molecular
clouds that are filtered out by the lack of short-spacing baselines in our data. The latter scenario would be consistent with the
increase in CO/CO (21) ratio with radius reported in Sect. 3.2.4. This scenario may
have similarities to that found for NGC~6946 Meier & Turner (2004), where molecular gas outside the central ring was found to be of
lower density and more diffuse than the dense, self-gravitating gas in the ring. This is reflected as a rapid change in the interstellar
medium (ISM) from dense star-forming cores in a central ring to diffuse, low-density molecular gas in and behind the molecular arms. This
results in CO/CO (21) intensity ratios in the central ring that are lower than in the diffuse gas on larger
scales - which also seems to be the case for NGC 1614 (see Sect. 3.2.4). Therefore, the absence of GMAs
outside the ring leeds us to conclude that the GMAs in NGC 1614 are formed as a result of some form of orbit-crowding within the nuclear
ring from gas funneled from the dust lane reservoir.
Analogies for the GMA formation processes could be M~51 (Koda et al. 2009) or IC~342 (Hirota et al. 2011, see also Sect. 4.3.3 for a more detailed comparison with other GMA-bearing galaxies). In M~51 the GMAs are only found in the spiral arms. They are formed upon entering the spiral arms through coagulation triggered by spiral-arm streaming motions and are disrupted again by fragmentation upon leaving the spiral arms.
For NGC 1614 we propose the following scenario: Smaller sized GMCs are funneled along the umbilical cords toward the molecular ring. At the transition from the dust lane to the molecular ring we assume that crowding processes (e.g. orbital crowding) in the mergers take place that lead to the collisional coagulation of individual GMCs. These cloud-cloud collisions could also trigger the onset of star formation (Scoville et al. 1986; Tan 2000) - the dissipation of excess kinetic energy during the collisions is a prerequisite for star formation (Hirota et al. 2011). Observations at higher angular resolution and sensitivity are necessary to resolve the molecular gas properties of the low surface brightness GMCs in the umbilical cords.
4.2 Tracers of star formation
4.2.1 Annular profiles
In Fig. 6 we present the annular profiles of Pa, CO (21), radio - (1.4 GHz), - (5 GHz), and
(8.4 GHz)-band. The averaged annular profiles for the NGC~1614 ring are shown for all tracers
(Fig. 6 upper left panel). We simply azimuthally averaged the flux density in a series of 0.1″-wide
annular circular rings (i.e., we excluded inclination effects). We did not exclude any parts of the annulus that might be atypical in
any band. We also set the ring center in the middle of the Pa where there is a slight central feature and used this position in all bands. This matches the weak - and -band central feature quite well.
We also divided the ring into four sections, north (N), east (E), west (W) and south (S) and plotted semi-annuli to study any positional shift in the relative radii of all five tracers (Fig. 6).
For the average annuli we see that the -band profile is the tracer with the strongest relative nuclear peak. The ring (or ring-like structure) is visible in -band (indicated by the step in the profile) but not quite as pronounced and possibly not with exactly the same central position. As expected, the radio - and -band data trace each other very well - and the ring radius may be shifted by 0.1″ compared to that of Pa - but the difference is very small and within errors. The average CO (21) profile is shifted outwards by 0.2″ from that of Pa. The nuclear region is obviously void of any CO (21) emission (that can be recovered in Vext by the SMA).
However, when inspecting the N, E, W, and S cuts, we find that there are large directional differences in the peak radii of the different tracers. All W cuts show a peak radius of 0.6″ for all tracers - the agreement is almost perfect and there is a well-defined local maximum for all tracers. In contrast, the strongest discrepancy is seen in the N cut between Pa and CO (21). Pa peaks at 0.3″ with a sharp profile while CO (21) peaks at 0.9″ with a large, rounded, extended (even double-peaked) structure. A similar but less pronounced effect can be seen in the S cut.
Part of the explanation for this discrepancy can be seen in Fig. 5 which clearly shows the north-south elongated distribution of CO (21) with respect to Pa. The correlation between the tracers is good in the west - but poorer in the north and south because of the additional contribution to the CO (21) emission that does not exclusively come from the ring itself, but also from the combined CO (21) contributions from the dust lane and its connection point to the ring.
The annuli are discussed in more detail in Sect. 4.3.1.
4.2.2 “Wildfire” vs. dynamics
There have been suggestions for different evolutionary tracks for the development of a starburst ring in NGC~1614. Alonso-Herrero et al. (2000)
suggested a “wildfire”-like scenario in which the ring was formed by a nuclear starburst that progresses outward and already has
consumed almost the whole gas in its center. In contrast to this, Olsson et al. (2010) proposed a scenario in which gas is piled up at an inner
Lindblad resonance to form the ring.
Supporting the wildfire origin of the starburst ring is the fact that the Pa and the radio continuum emission are distributed consistently inside the starburst ring, as seen in CO (21). The same is true for 1.6 m emission (-band) observed by Haan et al. (2011). The emission is centrally peaked like the Pa and the radio continuum emission. But in contrast to the other tracers, the 1.6 m emission is distributed in a broken-ring structure that has an opening to the northwest (Haan et al. 2011). Our observations of CO (21) show, however, that the molecular ring is open to the northeast. The radio continuum suggests, furthermore, that an older starburst lies in the ring center at the nucleus. An alternative explanation for the 1.4 GHz emission in the center could be an additional synchrotron component, possibly from a weak AGN. However, the distance between the tracers varies with the direction, which could also be indicative of spatial shifts between the tracers that are not caused by wildfires.
It was also suggested that the starburst ring in NGC 1614 could be caused by an inner Lindblad resonance (Olsson et al. 2010). One might argue that this process only takes place in barred galaxies, and it is still debated whether there is a stellar bar in NGC 1614 or not (see e.g., Olsson et al. 2010, who found an oval distortion in their -band data that is consistent with a bar). But galaxies with molecular rings have been found that have no bars (e.g. NGC~278 & NGC~7217, Knapen et al. 2004; Sarzi et al. 2007). It has been shown that other physical properties, e.g., a weak oval distortion (Buta et al. 1995), a dissolved bar (Verdes-Montenegro et al. 1995), tidal effects of a companion galaxy (Combes 1988a), or a recent minor merger (Knapen et al. 2004; Mazzuca et al. 2006) can cause the same effect as a bar: the non-axisymmetric potential caused by the bar can channel gas to the nucleus, which is slowed down near or at the inner Lindblad resonance. Elmegreen (1994) proposed a scenario in which the gravitational collapse of gas that is caught in an inner Lindblad resonance within a starburst ring, fueled by galactic bars or strong spiral arms, can lead to the formation of giant clouds with spots of intense star formation. We do see these hot spots in Figs. 1, 2, and 5. Another argument against the wildfire scenario is the asymmetry we found in the CO (21) emission, whereas the distribution of star formation tracers, such as Pa and radio continuum emission, is distributed very symmetrically. Conversely, Pa and the radio continuum agree very well with each other but not with the CO (21) distribution, which supports the idea of a wildfire because CO traces gas properties that occur earlier, i.e., farther outward, in the star formation process than those traced by Pa and the radio continuum.
An alternative theory that can possibly explain the spread in distribution of the star formation tracers was suggested by Lou et al. (2001a), Lou et al. (2001b), and Lou (2003). These authors highlight the importance of density waves in spiral galaxies as triggers for AGN and star-bursting events: a bar or companion induces spiral density waves at the inner Lindblad resonance (Lou et al. 2001a, and references therein). Lou et al. (2001b) found that fast magneto-hydrodynamic density waves (FMDWs) play an important dynamic role in developing the circumnuclear region of spiral galaxies with bars. The critical parameter seems to be the dampening of the FMDWs, which can lead to the onset of AGN activity (weak density wave dampening), to the formation of circumnuclear rings (strong dampening) and a mixture of circumnuclear and AGN activity (intermediate dampening) inside the inner Lindblad resonance. The two latter processes might play a role in NGC 1614.
In summary, there is no clear tendency toward the wildfire scenario, density waves, or the idea of a Lindblad resonance as a cause for the starburst ring in NGC 1614. From our study of the connection between the molecular gas in the ring and the galaxies dust lanes, we propose the following. Star formation in the molecular ring occurs at a Lindblad resonance and is replenished with molecular gas via the dust lanes. This is a result of orbital crowding or density wave phenomena. The morphology of the ring as seen in Pa could be caused by a wildfire independent of this star formation. We suggest that both processes coexist at the same time, but that the molecular ring is responsible for the starburst. To be able to distinguish between the three competing models, several observations might be useful: CO (10) observations with a resolution higher than currently available would enable us to study the possible connection between the CO (10) and the CO (21) gas reservoirs, CO data with a higher sensitivity would provide insight into the low surface brightness gas content of NGC 1614 and would allow us to continue studying the GMA formation, and observations of shock tracers, such as SiO, would help tracing a possible outward movement of the star formation by showing the places of shocked gas within the ISM.
|Galaxy||Hubble||Merger/||Ring radius||CO central||Ring444sym = symmetric, asym = asymmetric||Connection to555A number given denotes the number of connection points of the molecular ring to dust lanes, “–”
indicates that there is no information available
on a possible ring-dust lane connection.
666A number given denotes the number of GMAs found in the molecular rings and dust lanes, “–” indicates that
there is no information available
on GMAs in these objects.
this work, 2): Hsieh et al. (2008), 3): Hsieh et al. (2011), 4): Sakamoto et al. (2007), 5): Davies et al. (2004),
6): Nakai et al. (1987), 7): Weiß et al. (1999), 8): Matsushita et al. (2000)
|NGC 1614||64.2||SB(s)b pec||+||231||-||asym||1||10||1)|
|NGC 1097||14.5||SB(s)b||-||700||+||sym||2||19||2), 3)|
|NGC 7469||66||SBa||+888Part of an interacting galaxy pair with IC 5283.||800||+||asym||–||–||5)|
|M 82||3.9||I0||+999Part of an interacting galaxy pair with M 81.||200||+||sym||–||–||6)|
|105 70||-||asym||–||–||7), 8)|
4.3 Comparison with other galaxies
A comparison with other galaxies, such as NGC~4194 (the Medusa merger), NGC~1097, and NGC~1365 shows that there are similarities between these different types of galaxies and NGC~1614, but also that NGC 1614 does not fit in either in the same class as the Medusa or the class of barred spiral galaxies.
4.3.1 Comparison with other dust lane minor mergers
Similar properties as seen in NGC~1614 have been found for the Medusa merger. The two systems have in common minor axis
dust lanes and large amounts of molecular gas associated with thesee dust lanes (Beswick et al. 2005; Olsson et al. 2010; Aalto & Hüttemeister 2000), but not with the bulk of the
Apparently, the dust lanes are part of a feeding chain where the majority of the molecular gas has yet not engaged in star formation. However, a comparison of the CO distribution and the star formation properties shows a very different picture. in NGC 1614, which is a (S+s) minor merger, the compact CO emission is closely associated with the optical nuclear star formation (Olsson et al. 2010), while in NGC 4194, an E+S minor merger, much of the star formation is going on in super star clusters (SSCs) with a kpc-scale distribution (Weistrop et al. 2004) away from the molecular gas. No obvious correlation between these young SSCs (5 15 Myr) and the CO can be found (Aalto et al. 2010). In addition, both minor mergers have a relatively high CO/HCN line ratio (CO/HCN (NGC 1614) 43, CO/HCN (NGC~4194) 61; Costagliola et al. 2011), which indicates limited reservoirs of dense gas.
With the Vext data we recover only a fraction of the CO (21) emission in the center of NGC 1614. When we compare our data with the lower resolution SMA data by Wilson et al. (2008), we find that we missed almost all, 70%, of the CO 21 flux in the crossing minor axis dust lane in the north.
Possibly a large fraction of the missing flux is in the form of diffuse extended emission, associated with gas flows along the minor axis dust lane. This would be a similar situation as the one suggested for the NGC 4194 minor merger by Aalto et al. (2010).
It is still unknown whether the differences in the gas distribution and starbursts are due to different merger histories, variations in the process of merging spiral galaxies with spiral galaxies or elliptical with spiral (E+S) galaxies, or diverging properties within the evolution of mergers. A careful study of several of these objects might shed some light on their evolutionary paths.
4.3.2 Comparing NGC 1614 with galaxies with nuclear molecular rings
NGC~1614 has some features in common with local barred spiral galaxies, such as NGC~1097 and NGC~1365, but
unlike NGC 1614, these galaxies are not classified as mergers. These barred spirals are so-called twin-peak galaxies because they contain
nuclear molecular rings that show a double-peak structure in their molecular gas distribution perpendicular to the stellar bar, where the
starburst ring and the dust lanes connect (Sakamoto et al. 2007; Hsieh et al. 2011). NGC 1614 harbors an asymmetric ring and displays at least one
connection point with its dust lanes. This is comparable to the twin connection points seen in local barred spirals (Kenney et al. 1992). However, some barred spirals, such as NGC 1365, do show comparable asymmetric central starburst rings (Sakamoto et al. 2007).
The molecular gas in all three objects is associated with dust lanes close to or connected to the nuclear starburst rings. Within the rings all sources show hot spots in the CO emission that possibly contain SSCs in their early stages (Sakamoto et al. 2007). An obvious difference in the CO gas distribution in the starburst ring though is the absence of CO emission in NGC 1614 at the center of the ring, as opposed to the nuclear ring in NGC 1097. This might be an indication that the central structure in NGC 1614 is older than that in NGC 1097 and therefore the gas there is already consumed or swept away by developing stars. Another possibility is that the ring and AGN activity in NGC 1097 were caused by intermediately dampened spiral density waves (Lou 2003), whereas the dampening in NGC 1614 was strong enough only to form the starburst ring (see also Sect. 4.2.2).
With sizes (radii) of 700 pc and 1 kpc, the nuclear rings in NGC 1097 and NGC 1365 also have a much larger extension then that of NGC 1614 (231 pc), possibly indicating that NGC 1614 harbors a scaled-down version of a nuclear starburst ring in development caused by the merger event.
Another prominent galaxy with a molecular ring is M~82. In contrast to NGC~1614, M 82 is not a merger, but constitutes an integral part of an interacting galaxy pair together with M~81. M 82 is an edge-on starburst galaxy with a nuclear molecular ring or torus (radius 200 pc; Nakai et al. 1987; Shen & Lo 1995) that corresponds to current sites of star formation (Nakai et al. 1987), like in NGC 1614. In M 82 in contrast to NGC 1614, however, the star formation is also associated with SSCs (Telesco & Gezari 1992), like in NGC~4194. The molecular ring in M 82 contains less than 30 of the total molecular CO mass in M 82. The ring in NGC 1614 shows a mass fraction of 30 when taking the total CO (21) mass from Wilson et al. (2008). In addition, the CO / CO (21) ratios for M 82 ( 14.2, Mao et al. 2000) and NGC 1614 (15, see Sect. 3.2.4) are quite similar. A more recent study by Weiß et al. (1999) discovered an additional CO feature in M 82 that is suggested to be a superbubble driven by a starburst in its center (Matsushita et al. 2000), similar to the wildfire mode proposed by Alonso-Herrero et al. (2001) for NGC 1614.
Yet another example is the Seyfert 1 galaxy NGC~7469, which is the companion of IC~5283 in an interacting pair of similar sized galaxies. For NGC 7469 several tracers (e.g. 5 GHz radio continuum, , , imaging, PAHs, CO (21)) were found to be located in a circumnuclear ring/spiral arm system (Wilson et al. 1991; Mazzarella et al. 1994; Miles et al. 1994; Genzel et al. 1995; Davies et al. 2004). A comparison with the starburst and molecular rings in NGC 1614 shows some prominent similarities. The near-infrared (NIR) emission detected in the -band continuum, associated with the radio continuum and other tracers, is found inside the CO emission distribution (Davies et al. 2004). For NGC 7469 however, the CO (21) distribution shows a central peak/nuclear ring and a bar or pair of spiral arms that seems to connect to the central emission/circumnuclear ring itself (Davies et al. 2004). In NGC 1614 this is not the case.
4.3.3 GMAs in galaxies with nuclear molecular rings
One other property NGC~1614 shares with galaxies like NGC~1097 and NGC~1365 are hot spots
(gas clumps and CO brightness temperature peaks, Sakamoto et al. 2007; Hsieh et al. 2011), so-called GMAs (Vogel et al. 1988). For NGC 1365, Sakamoto et al. (2007)
found five CO hot spots with line widths (FWHM) between 60 and 90 km s that seem to be associated with SSCs of
10 M (Galliano et al. 2005).
Hsieh et al. (2011) found many GMAs in the molecular ring (14 GMAs) and the dust lane (5 GMAs) of NGC 1097. They subdivided the GMAs in the molecular ring into two categories: GMAs with broad ( 30 km s) and narrow ( 30 km s) line dispersions. The line widths for the narrow GMAs ( 55.1 km s) are almost only half the line widths of broad-line GMAs ( 102.0 km s) but the average molecular gas mass is almost the same ((N1-N11) 8.1 10 M, (B1-B3) 9.3 10 M). The GMAs in the dust lane only have slightly more narrow line widths ( 93.8 km s) and slightly higher molecular gas masses ((D1-D5) 1.0 10 M) than the GMAs in the ring. For NGC 1614 the average line widths and masses are 129.5 km s and (ring) 7.0 10 M for GMAs in the ring, and 110.7 km s and (outside the ring) 5.0 10 M for GMAs outside the molecular ring (Table 3) - the average line widths are broader and the masses are slightly higher for GMAs in the ring than for those outside the ring in NGC 1614. At least two of the GMAs outside the ring in NGC 1614 are associated with dust lanes (GMA 2, GMA 5, Fig. 1). Three of the GMAs in the ring (GMA 6, GMA 7, and GMA 8) are also associated with a dust lane, and GMA 3 is located in the feeding lines that connect the molecular ring with the dust lane (Fig. 1).
Hsieh et al. (2011) also determined the sizes of their GMAs by measuring the number of pixels above a certain intensity threshold and then assuming the GMAs to be spherical to obtain diameters of the individual GMAs. The sizes of GMAs located in the ring of NGC 1097 and outside the structure differ significantly. The average size of a GMA situated in the molecular ring is 220 pc. GMAs that are situated outside the ring but in the dust lane in NGC 1097, however, have an average size of 269 pc, which is 20 larger than for GMAs in the ring. For NGC 1614, however, this seems to be less pronounced. The average size of GMAs in the ring is 259 24 pc (assuming spherical morphology) and for GMAs outside the ring and associated with dust lanes it is 273 20 pc. The GMAs in the ring in NGC 1614 seem to be larger than those found in NGC 1097 (Hsieh et al. 2011). Although the average spherical sizes for GMAs in and outside the ring are comparable and agree within the error limits, the sizes for the GMAs outside the ring that are associated with dust lanes are most likely overestimated. The sizes of these GMAs are measured at low surface brightnesses and therefore likely to be afflicted with larger errors. This does not change the fact that Hsieh et al. (2011) discovered that for NGC 1097, GMAs in the ring are significantly smaller than those outside the ring, whereas for NGC 1614 it is most likely the other way around, or possibly the sizes in and outside the ring are about the same. One deciding factor here might be the formation history of the starburst rings in these two galaxies e.g., NGC 1614 is a merger, NGC 1097 is not.
We summarize the comparisons between the galaxies in this section in Table 4. This shows that NGC 1614 shares properties with several of the comparison galaxies - e.g., an assymetric molecular ring, dust lanes, GMAs, and a radial displacement between starburst and molecular rings. However, we always found significant differences as well. Therefore, NGC 1614 cannot be classified as a twin-peak galaxy, a typical barred spiral, or a typical merger of any flavor; it remains a very interesting yet puzzling object.
We have observed the 1.3 mm continuum and the CO and CO (21) transitions with high resolution toward the center of NGC~1614.
We detected an asymmetric distribution of the CO (21) emission in a ring of 231 pc radius with a molecular gas mass of 8.3 10 M. The center of the ring is free of CO emission.
We did not detect CO (21) with our observational setup. The lower limit of the CO/CO (21) line ratio is 15.
We identified individual GMAs in the CO data cube, some associated with the ring, some associated with the dust lanes. The GMAs in the ring have a higher dispersion and more mass than those outside the ring.
Molecular gas in the dust lane shows little evidence of star formation. While GMAs appear in the umbilical-cord connection between the dust lane and starburst ring, there is no evidence of substantial levels of star formation (e.g., Fig. 5). This behavior is suggestive of a transitory process in which GMAs form as they move toward the nuclear ring, but intense star formation only takes place once the molecular material merges into the ring.
There is no general consensus yet for one scenario to explain the formation of the circumnuclear ring. Considering our results, however, we suggest that the starburst ring is caused by a Lindblad resonance and fueled by gas moving via the dust lanes onto the ring.
In conclusion, we found that NGC 1614 is a good target for studying the impact of a minor merger on the state of the molecular gas and the formation of circumnuclear starburst rings. We proposed a way to explain how nuclear gas structures in minor mergers are fed. NGC 1614 may be an example of a merger where gas in tidally induced dust structures has been found to physically connect to a starbursting structure in the very nucleus of the galaxy. We also suggested that the GMAs in NGC 1614 are formed in situ in the ring through a combination of cloud-cloud collisions and crowding processes that facilitate the onset of star formation. Owing to the limited spatial resolution and surface brightness sensitivity, many questions remain to be answered, such as in which way the larger scale CO (10) is connected to the starburst, and if there is a connection to the proposed nuclear outflow (Olsson et al. 2010), or how the kinematics of cold molecular gas affect the evolution of star formation and nuclear activity in minor mergers in general. A more detailed study of the molecular gas content of NGC 1614 on different spatial scales is needed to address these questions. New studies with instruments such as ALMA will enable us to answer them.
Acknowledgements.We are grateful to Glen Petitpas for his help with the SMA observations. The Dark Cosmology Centre is funded by the Danish National Research Foundation. SA thanks the Swedish Research Council and the Swedish National Space Board for support. JSG acknowledges partial support for this research from NSF grant AST0708967 from the University of Wisconsin-Madison. The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and the Academia Sinica. MERLIN is a national facility operated by The University of Manchester on behalf of the Science and Technology Facilities Council (STFC). The VLA is operated by the National Radio Astronomy Observatory, a facility of the National Science foundation operated under cooperative agreement by Associated Universities, Inc. This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. This research used the facilities of the Canadian Astronomy Data Centre operated by the National Research Council of Canada with the support of the Canadian Space Agency.
- Aalto et al. (2010) Aalto, S., Beswick, R., & Jütte, E. 2010, A&A, 522, A59
- Aalto & Hüttemeister (2000) Aalto, S. & Hüttemeister, S. 2000, A&A, 362, 42
- Alonso-Herrero et al. (2001) Alonso-Herrero, A., Engelbracht, C. W., Rieke, M. J., Rieke, G. H., & Quillen, A. C. 2001, ApJ, 546, 952
- Alonso-Herrero et al. (2000) Alonso-Herrero, A., Rieke, M. J., Rieke, G. H., & Shields, J. C. 2000, ApJ, 530, 688
- Beswick et al. (2005) Beswick, R. J., Aalto, S., Pedlar, A., & Hüttemeister, S. 2005, A&A, 444, 791
- Bournaud et al. (2005) Bournaud, F., Jog, C. J., & Combes, F. 2005, A&A, 437, 69
- Buta et al. (1995) Buta, R., van Driel, W., Braine, J., et al. 1995, ApJ, 450, 593
- Chapelon et al. (1999) Chapelon, S., Contini, T., & Davoust, E. 1999, A&A, 345, 81
- Combes (1988a) Combes, F. 1988a, in Lecture Notes in Physics, Berlin Springer Verlag, Vol. 315, Molecular Clouds, Milky-Way and External Galaxies, ed. R. L. Dickman, R. L. Snell, & J. S. Young, 441
- Combes (1988b) Combes, F. 1988b, in NATO ASIC Proc. 232: Galactic and Extragalactic Star Formation, ed. R. E. Pudritz & M. Fich, 475–+
- Costagliola et al. (2011) Costagliola, F., Aalto, S., Rodriguez, M. I., et al. 2011, A&A, 528, A30
- Dasyra et al. (2006) Dasyra, K. M., Tacconi, L. J., Davies, R. I., et al. 2006, ApJ, 651, 835
- Davies et al. (2004) Davies, R. I., Tacconi, L. J., & Genzel, R. 2004, ApJ, 602, 148
- de Vaucouleurs et al. (1991) de Vaucouleurs, G., de Vaucouleurs, A., Corwin, Jr., H. G., et al. 1991, Third Reference Catalogue of Bright Galaxies. Volume I: Explanations and references. Volume II: Data for galaxies between 0 and 12. Volume III: Data for galaxies between 12 and 24.
- Elmegreen (1994) Elmegreen, B. G. 1994, ApJ, 425, L73
- Galliano et al. (2005) Galliano, E., Alloin, D., Pantin, E., Lagage, P. O., & Marco, O. 2005, A&A, 438, 803
- Genzel et al. (1995) Genzel, R., Weitzel, L., Tacconi-Garman, L. E., et al. 1995, ApJ, 444, 129
- Haan et al. (2011) Haan, S., Surace, J. A., Armus, L., et al. 2011, AJ, 141, 100
- Hirota et al. (2011) Hirota, A., Kuno, N., Sato, N., et al. 2011, ApJ, 737, 40
- Ho (1999) Ho, L. C. 1999, Advances in Space Research, 23, 813
- Hsieh et al. (2008) Hsieh, P.-Y., Matsushita, S., Lim, J., Kohno, K., & Sawada-Satoh, S. 2008, ApJ, 683, 70
- Hsieh et al. (2011) Hsieh, P.-Y., Matsushita, S., Liu, G., et al. 2011, ApJ, 736, 129
- Hunsberger et al. (1996) Hunsberger, S. D., Charlton, J. C., & Zaritsky, D. 1996, ApJ, 462, 50
- Jesseit et al. (2007) Jesseit, R., Naab, T., Peletier, R. F., & Burkert, A. 2007, MNRAS, 376, 997
- Kenney et al. (1992) Kenney, J. D. P., Wilson, C. D., Scoville, N. Z., Devereux, N. A., & Young, J. S. 1992, ApJ, 395, L79
- Knapen et al. (2004) Knapen, J. H., Whyte, L. F., de Blok, W. J. G., & van der Hulst, J. M. 2004, A&A, 423, 481
- Koda et al. (2009) Koda, J., Scoville, N., Sawada, T., et al. 2009, ApJ, 700, L132
- Kormendy et al. (2009) Kormendy, J., Fisher, D. B., Cornell, M. E., & Bender, R. 2009, ApJS, 182, 216
- Lou (2003) Lou, Y. Q. 2003, Acta Astronomica Sinica, 44, 172
- Lou et al. (2001a) Lou, Y.-Q., Yuan, C., & Fan, Z. 2001a, ApJ, 552, 189
- Lou et al. (2001b) Lou, Y.-Q., Yuan, C., Fan, Z., & Leon, S. 2001b, ApJ, 553, L35
- Mao et al. (2000) Mao, R. Q., Henkel, C., Schulz, A., et al. 2000, A&A, 358, 433
- Matsushita et al. (2000) Matsushita, S., Kawabe, R., Matsumoto, H., et al. 2000, ApJ, 545, L107
- Mazzarella et al. (1994) Mazzarella, J. M., Voit, G. M., Soifer, B. T., et al. 1994, AJ, 107, 1274
- Mazzuca et al. (2006) Mazzuca, L. M., Sarzi, M., Knapen, J. H., Veilleux, S., & Swaters, R. 2006, ApJ, 649, L79
- Meier & Turner (2004) Meier, D. S. & Turner, J. L. 2004, AJ, 127, 2069
- Miles et al. (1994) Miles, J. W., Houck, J. R., & Hayward, T. L. 1994, ApJ, 425, L37
- Monreal-Ibero et al. (2006) Monreal-Ibero, A., Arribas, S., & Colina, L. 2006, ApJ, 637, 138
- Nakai et al. (1987) Nakai, N., Hayashi, M., Handa, T., et al. 1987, PASJ, 39, 685
- Narayanan et al. (2011) Narayanan, D., Krumholz, M., Ostriker, E. C., & Hernquist, L. 2011, MNRAS, 418, 664
- Neff et al. (1990) Neff, S. G., Hutchings, J. B., Standord, S. A., & Unger, S. W. 1990, AJ, 99, 1088
- Olsson et al. (2010) Olsson, E., Aalto, S., Thomasson, M., & Beswick, R. 2010, A&A, 513, A11+
- Puxley & Brand (1999) Puxley, P. J. & Brand, P. W. J. L. 1999, ApJ, 514, 675
- Sakamoto et al. (2007) Sakamoto, K., Ho, P. T. P., Mao, R.-Q., Matsushita, S., & Peck, A. B. 2007, ApJ, 654, 782
- Sanders et al. (2003) Sanders, D. B., Mazzarella, J. M., Kim, D.-C., Surace, J. A., & Soifer, B. T. 2003, AJ, 126, 1607
- Sarzi et al. (2007) Sarzi, M., Allard, E. L., Knapen, J. H., & Mazzuca, L. M. 2007, MNRAS, 380, 949
- Scoville et al. (1986) Scoville, N. Z., Sanders, D. B., & Clemens, D. P. 1986, ApJ, 310, L77
- Scoville et al. (1985) Scoville, N. Z., Soifer, B. T., Neugebauer, G., et al. 1985, ApJ, 289, 129
- Shen & Lo (1995) Shen, J. & Lo, K. Y. 1995, ApJ, 445, L99
- Shlosman et al. (1989) Shlosman, I., Frank, J., & Begelman, M. C. 1989, Nature, 338, 45
- Simkin et al. (1980) Simkin, S. M., Su, H. J., & Schwarz, M. P. 1980, ApJ, 237, 404
- Somerville et al. (2008) Somerville, R. S., Hopkins, P. F., Cox, T. J., Robertson, B. E., & Hernquist, L. 2008, MNRAS, 391, 481
- Stutzki & Guesten (1990) Stutzki, J. & Guesten, R. 1990, ApJ, 356, 513
- Tacconi et al. (2002) Tacconi, L. J., Genzel, R., Lutz, D., et al. 2002, ApJ, 580, 73
- Tan (2000) Tan, J. C. 2000, ApJ, 536, 173
- Telesco & Gezari (1992) Telesco, C. M. & Gezari, D. Y. 1992, ApJ, 395, 461
- Terashima et al. (2000) Terashima, Y., Ho, L. C., Ptak, A. F., et al. 2000, ApJ, 533, 729
- Toomre & Toomre (1972) Toomre, A. & Toomre, J. 1972, ApJ, 178, 623
- Toth & Ostriker (1992) Toth, G. & Ostriker, J. P. 1992, ApJ, 389, 5
- Väisänen et al. (2012) Väisänen, P., Rajpaul, V., Zijlstra, A. A., Reunanen, J., & Kotilainen, J. 2012, MNRAS, 420, 2209
- Veilleux et al. (1999) Veilleux, S., Sanders, D. B., & Kim, D.-C. 1999, ApJ, 522, 139
- Verdes-Montenegro et al. (1995) Verdes-Montenegro, L., Bosma, A., & Athanassoula, E. 1995, A&A, 300, 65
- Vogel et al. (1988) Vogel, S. N., Kulkarni, S. R., & Scoville, N. Z. 1988, Nature, 334, 402
- Weiß et al. (1999) Weiß, A., Walter, F., Neininger, N., & Klein, U. 1999, A&A, 345, L23
- Weistrop et al. (2004) Weistrop, D., Eggers, D., Hancock, M., et al. 2004, AJ, 127, 1360
- Wilson et al. (1991) Wilson, A. S., Helfer, T. T., Haniff, C. A., & Ward, M. J. 1991, ApJ, 381, 79
- Wilson et al. (2008) Wilson, C. D., Petitpas, G. R., Iono, D., et al. 2008, ApJS, 178, 189
- Zaritsky (1995) Zaritsky, D. 1995, ApJ, 448, L17