Star Clusters in Arp 284

Star Clusters in the Interacting Galaxy System Arp 284

Bradley W. Peterson, Curtis Struck, Beverly J. Smith, and Mark Hancock
Department of Physics and Astronomy, Iowa State University, Ames, IA 50011, USA
Department of Physics and Astronomy, East Tennessee State University, Johnson City, TN 37614, USA
Institute of Geophysics and Planetary Physics, University of California, Riverside, Riverside, CA 92521, USA
E-mail: (BWP); (CS)
Accepted YYYY MONTH DAY. Received YYYY MONTH DAY; in original form YYYY MONTH DAY

We present results from a study of proto-globular cluster candidates in the interacting galaxy system Arp 284 (NGC 7714/5) using data from the Hubble Space Telescope. Previous studies of the Antennae and M51 have suggested that the majority of young massive star clusters dissolve within 20 Myr due to mass loss. We use the evolutionary synthesis code starburst99 to estimate ages and extinctions for approximately 175 clusters visible with HST. We also use lower-resolution GALEX and ground-based H data to estimate the ages of the giant H ii regions in which these clusters are found, and compare the Spitzer colours of these H ii regions to those of star forming regions in other interacting systems. The ages are also used to aid in the interpretation of Chandra X-ray data.

Clusters in the tidal tails of NGC 7714 are generally found to have ages less than 20 Myr, though observational limits make the significance of this result uncertain. Older clusters, though not numerous, have nearly the same spatial distribution within the imaged portion of NGC 7714 as young clusters. The cluster population in the bridge connecting the two galaxies appears to be older, but the data in this part of the system is too limited to draw firm conclusions. The ages of the giant H ii regions in NGC 7714 are generally older than those of their constituent clusters, possibly indicating that the young clusters we detect are surrounded by their dispersed predecessors.

galaxies: individual (NGC 7714, NGC 7715) – galaxies: interactions – galaxies: star clusters
pagerange: Star Clusters in the Interacting Galaxy System Arp 284Star Clusters in the Interacting Galaxy System Arp 284pubyear: 2009

1 Introduction

The study of star formation in interacting galaxies has been revolutionized by the Hubble Space Telescope. One of the early results to emerge from HST was the discovery of a class of bright, blue clusters 300 Myr old in the peculiar elliptical galaxy NGC 1275 (Holtzman et al., 1992). These young massive clusters (YMCs) were also found in the prototypical merger remnant NGC 7252, with ages ranging from 34–500 Myr, consistent with the age of the merger (Whitmore et al., 1993). Many more examples soon followed, in which it was found that the properties of these YMCs implied that they might be young globular clusters, though with higher abundances (Conti & Vacca, 1994; Hunter et al., 1994).

In interacting galaxy systems, in which star forming regions are abundant, the high resolution of HST has allowed the study of star cluster demographics. Among the results to emerge from such studies is the apparent dissolution of young star clusters, also known as infant mortality. Fall et al. (2005) used evolutionary synthesis models to estimate the ages of clusters in the Antennae and found that 90 per cent of clusters vanish over each age dex from up to yr, suggesting that cluster dissolution continues over very long time-scales. In contrast, Bastian et al. (2005a) found in M51 that 70 per cent of clusters dissolve in the first 20 Myr, indicating that disruption is concentrated on short time-scales. Recent spectroscopic observations of star clusters in the Antennae, combined with more careful modeling of the star formation history of the system, does not show the long-term cluster mortality previously claimed (Bastian et al., 2009), thus this issue is still controversial.

Additional evidence for a short dissolution time-scale has been found in comparisons of cluster and field populations. Tremonti et al. (2001) obtained UV spectra and evolutionary synthesis models of both clusters and diffuse regions of the dwarf starburst galaxy NGC 5253. Unlike the clusters, the diffuse region spectra lacked indications of O-type stars. It was suggested that stars form entirely within clusters, which dissolve on  Myr time-scales, preventing short-lived O stars from ever entering the field population.

Early-stage disruption has been ascribed largely to the loss of gas left over from star formation, which is swept out of the cluster by stellar winds and massive star supernovae. The resulting mass loss drops the escape velocity below the stellar velocities, causing the cluster to expand and decrease in surface brightness. The departing stars disperse into the field population in few yr (Bastian & Goodwin, 2006; Goodwin & Bastian, 2006). On longer time-scales ( yr and longer), interactions with giant molecular clouds (Gieles et al., 2006) and tidal effects are expected to become significant disruption mechanisms.

If YMCs are progenitors of globular clusters, some of them must survive their most massive stars, as well as external effects that may lead to dissolution (de Grijs & Parmentier, 2007). Interacting galaxy systems have high numbers of YMCs, making them ideal targets for cluster population studies. To date, too few interacting systems have been observed to definitively determine how rapidly and by what mechanism massive clusters dissolve. Further, both the Antennae and M51 are advanced interactions with complex star formation histories. It is necessary to confirm infant mortality in early-stage interactions where the star formation history is simpler, with a particular emphasis on covering a wide range of tidal environments.

In this paper we focus on the early-stage interacting system Arp 284, which at a distance of  Mpc ( pc) according to the NASA/IPAC Extragalactic Database (NED) has nearly 175 detected star forming regions in archival HST images. We use HST colours and the evolutionary synthesis code starburst99 (sb99) to obtain age and reddening estimates for star forming regions. We also use lower-resolution UV data from the Galaxy Evolution Explorer (GALEX; Martin et al. 2005) and ground-based H to obtain age estimates for larger H ii regions, and mid-IR colours from the Spitzer Space Telescope Infrared Array Camera (IRAC; Werner et al. 2002; Fazio et al. 2004) to compare the H ii regions to IR emission clumps found in other interacting systems. These age results are used to aid in the interpretation of the X-ray data of Smith et al. (2005a),which was obtained using Chandra (Weisskopf et al., 2002).

2 Tidal Star Forming Environments of Arp 284

The peculiar morphology of the Arp 284 system results from an off-centre collision between two disc galaxies, with the primary (NGC 7714) having roughly three times the mass of the companion (NGC 7715). The time of closest approach was  Myr ago. For a discussion of the interaction history of the system, the reader is directed to the models of Struck & Smith (2003). In this section we summarize previous observations of the Arp 284 system that relate to its star formation history. Most of the tidal features are labelled in Fig. 1, while the H ii region complexes in the primary disc are labelled in Fig. 2.

2.1 Nucleus of NGC 7714

The NGC 7714 nucleus was the first object to be described using the term ‘starburst,’ although there are earlier references in the literature to ‘bursts of star formation’ (Weedman et al., 1981). Strong H emission indicates a large nuclear H ii region and ongoing star formation (González Delgado et al., 1995; Smith et al., 1997). O and B stars are required to fit the UV spectrum, with a best fitting age 5 Myr and a burst of star formation favored over continuous star formation (González Delgado et al., 1999).

Lançon et al. (2001) found that a second, older nuclear population is suggested by strong absorption bands of CO in the near-IR. In particular, the K-band spectrum is compatible with a combination of red supergiants, AGB stars, and RGB stars. They concluded that a two-component model is probably an oversimplification of the nuclear star formation history, preferring a star formation rate that has decreased over several yr, either through a smooth exponential decrease or a series of successively smaller bursts. The burst model best fits the observed supernova rate with a burst 15–50 Myr ago, in addition to the most recent burst 5 Myr ago.

Infrared Space Observatory (ISO) and Spitzer mid-IR spectra of the nuclear region show strong features from polycyclic aromatic hydrocarbons (PAHs) as well as bright forbidden lines (O’Halloran et al., 2000; Brandl et al., 2004). The metallicity of NGC 7714 peaks in the nucleus, with a value near , and generally decreases with distance from the nucleus (González Delgado et al., 1995).

2.2 Stellar Ring of NGC 7714 (Feature 1)

NGC 7714 has a partial ring-like structure to the east. This feature does not show prominent H emission, reflecting a lack of ongoing star formation, while the red continuum suggests an old stellar population (González Delgado et al., 1995; Smith et al., 1997). Low levels of H i emission indicate that there is little hydrogen gas (Smith et al., 1997).

2.3 Southern Tails of NGC 7714
(Features 2 and 3)

Two parallel tidal tails lie to the southwest of the NGC 7714 nucleus. The outer southern tail, faintly visible in optical images, is revealed by 21 cm maps to be part of a large H i loop which reconnects with the galaxy in the northern tail (Smith & Wallin, 1992; Smith et al., 1997). The loop is also visible in GALEX FUV, which is shown smoothed with overlaid H i contours in Fig. 3. This gas loop is an entirely separate structure from the inner southern tail, which hosts significant ongoing star formation.

The base of the inner southern tail just south of the nucleus contains H ii region A (Smith et al., 2005a). This region is visible in 20 cm radio continuum (González Delgado et al., 1995; Smith et al., 1997). The HST images reveal numerous bright clusters along this structure, with another large H ii region at the edge of the field (region C; González Delgado et al. 1995), coincident with a bright X-ray source (Smith et al., 2005a). The ages of regions A and C have been estimated by García-Vargas et al. (1997) at and  Myr based on photoionization models. These regions also appear in Spitzer images (Smith et al., 2007), to be discussed below.

González Delgado et al. (1995) found H region A to have a metallicity comparable to the nucleus, about . Region C, at a larger distance from the nucleus, shows a lower metallicity of .

2.4 Northern Tail of NGC 7714 (Feature 4)

The northern tail of the primary galaxy is visible but not particularly prominent in optical photographs. Near the base of the tail, to the northwest of the nucleus, are H regions B, D, and E. These regions are also detected in 21 cm H i emission and red continuum (Smith et al., 1997). HST resolves these regions into several smaller groups of clusters. The age of region D has been estimated at  Myr using photoionisation modeling (García-Vargas et al., 1997).

H region B has been found to have lower metallicty than the nucleus, about (González Delgado et al., 1995).

2.5 Bridge (Feature 5)

The bridge consists of a northern and a southern component. The northern bridge shows ‘clumpy’ emission from H i gas, and several H peaks are also found, offset slightly to the south, indicating ongoing star formation. The southern component lacks H i gas and H emission, but is readily seen in the red continuum, likely indicating an older stellar population (Smith et al., 1997). No metallicity data is currently available for the bridge.

2.6 Nucleus of NGC 7715

The lack of significant H emission in the NGC 7715 nucleus indicates that there is no ongoing massive star formation (González Delgado et al., 1995; Smith et al., 1997). Emission from the 21 cm H i radio line is also minimal, and offset to the north of the optical galaxy. The nucleus is comparatively bright in the red continuum, suggesting an aged stellar population (Smith et al., 1997). Spectral studies of the nucleus show a post-starburst spectrum, again confirming the lack of ongoing star formation (Bernlöhr, 1993).

2.7 Eastern Tail of NGC 7715
(Feature 6)

The eastern tail of NGC 7715 is faintly seen in optical images. Radio maps show an H i counterpart, offset slightly to the north (Smith & Wallin, 1992; Smith et al., 1997). Massive star formation is thought to be weak based on the lack of H emission. This feature has not been studied in detail, and it falls outside the field of the HST images. This poor coverage is unfortunate, because the models suggest that this tail curves behind the NGC 7715 disc, passes behind the bridge, and comes back into view as the faint outer southern tail of NGC 7714 (Struck & Smith, 2003). In this scenario, the tail contains a significant fraction of the original mass from the NGC 7715 disc (Struck & Smith, 2003).

3 Observations

3.1 Hst observations

The HST Wide Field and Planetary Camera 2 (WFPC2) has observed the Arp 284 system through several observing programs between 1995 May 13 and 2001 July 10. These programs targeted either the starburst nucleus or SN 1999dn in the disc of NGC 7714. The disparate nature of the observations means that many of the bright clusters are not visible in all filters, and in some cases the clusters were not found on the same chip across all filters. The field of view is on the PC chip and for the WF chips, with pixel sizes of pixel on the PC chip and pixel on the WF chips. The images and exposure times are summarized in Table 1. Note that the number in the filter name approximately corresponds to the central wavelength of the filter in nanometers. The available filters were F380W (), F555W (), F606W (), and F814W ().

The NGC 7714 nucleus and most of the inner southern tail are visible on the PC chip in all filters. One image in both F380W and F814W had to be discarded due to excessive cosmic ray contamination. This left three good images in the F380W filter, two in F555W, one in F606W, and three in F814W. In F814W, H ii region C was visible in two images on the WF4 chip.

The entire northern tail was visible on the PC chip in the F380W and F606W filters. In F555W and F814W, only 12 clusters were visible on the PC chip, with the remainder visible on WF4 in F814W.

The bridge was on the seam between the WF3 and WF4 chips in all four F380W images, so some of the clusters could not be measured in this band. In the other filters, it was always on the WF3 chip, with two images in F555W, one in F606W, and two in F814W.

3.2 H ii region observations

Archival GALEX and Spitzer IRAC data are available for the Arp 284 system. The GALEX images increase the number of colours that can be fit with sb99 models, while Spitzer colours can be compared to those of star forming regions in other interacting systems. These instruments have significantly lower resolution than HST, but still allow the H ii regions A–E to be resolved, as shown in Fig. 4. Each of these regions contains numerous star clusters resolved by HST. For clarity, the point-like sources in the HST data will be called ‘clusters’ while these larger structures will be called ‘H ii regions.’ With H luminosities of erg s (Smith et al., 2005a), these are properly classified as ‘giant H ii regions’ (e.g. Kennicutt, 1984).

The H ii regions were located in an H map, which was then registered to the other images. The H and accompanying R-band images were obtained using the 1.8 m Perkins Telescope of the Ohio State University, using a CCD with a pixel size of pixel. These data were published previously in Smith et al. (1997).

GALEX observations are available in both the FUV and NUV bands, with effective wavelengths of 1516 and 2267 Å. The pixel size is pixel, which makes even these relatively large H ii regions difficult to identify. The integration times are 4736 s in FUV and 520 s in NUV.

IR data from Spitzer are available for both the IRAC and Multiband Imaging Photometer for Spitzer (MIPS) instruments. The resolution in the MIPS bands is far too low for cluster studies, and there are pronounced artifacts due to the point spread function because of the brightness of the nucleus, so only the broadband IRAC bands (3.6, 4.5, 5.8, and 8.0) were used. The data reduction and mosaicing for the Spitzer data are described in Smith et al. (2007).

4 Hst Data Analysis

4.1 Data Reduction

All of the images were processed by the standard Hubble pipeline. The small number of exposures available in some data sets made cosmic ray removal a non-trivial task, since median combining two images is ineffective. First, we used the iraf111iraf is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. task cosmicrays to remove events a few pixels in size that were well above background. Next, images were blinked to identify remaining cosmic rays, which were removed manually using the fixpix task. Finally, images were combined using the imcombine task. For F606W, only one image was available, so other filters were used for comparison.

4.2 Cluster Selection

Clusters were selected using the iraf task daofind. The background level was determined using imstat in a relatively dark part of the image. The nuclear region is much brighter than the rest of the disc, so it was necessary to use a higher background for cluster detection, measured at the outer edge of the nucleus. We used daofind in each filter and then merged the lists to obtain the final sample. Clusters detected in only two filters were not included in the sample.

We detected 16 clusters in the nucleus at the level. The nucleus is quite crowded, so no attempt was made to push the detection threshold any lower. The lowest detected flux was erg s cm Å in F814W.

In the southern tail, 60 clusters were detected at the level and another 34 were detected at the level. These levels correspond to fluxes of 2.17 and erg s cm Å, respectively, in F814W. Before using the software, clusters in this part of the system were selected by eye. The set matched well with the results of the visual inspection, and were retained for the analysis.

In the northern tail, daofind detected 60 clusters at the level, with no additional clusters detected at the level. This is probably due to crowding in the major complexes, which resulted in multiple detections on single clusters at lower detection thresholds. These clusters may in fact represent multiple unresolved star clusters. In a few cases, multiple resolved clusters had to be measured in a single aperture because they were too close together to measure without aperture overlap.

Cluster selection was difficult in the bridge due to the higher background level in this part of the system and the lower resolution of the WF3 chip compared to the PC chip. For these reasons, it was necessary to lower the detection threshold in the bridge to , allowing the detection of 34 clusters. The faintest flux detected in F814W was erg s cm Å. Note that this is higher than the level on the PC chip. Only 23 of the clusters could be measured in F380W.

4.3 Photometry of clusters

Photometry was performed using the iraf task phot in the daophot package, with magnitudes calibrated to the Vega magnitude system using the zero point values given in the WFPC2 Data Handbook. The magnitudes and ages for the cluster sample are given in Table 2, while the extinctions and mass estimates are given in Table 3. The determination of the ages and extinctions will be discussed in § 5.2. The mass determination will be discussed in § 6.2.

The tidal tails generally fell on the PC chip. For these clusters, we used aperture radii of 3 or 4 pixels for single clusters, depending on crowding, and 6 pixels in six cases when barely-resolved clusters had to be measured in a single aperture. Sky brightness for background subtraction was determined by the mode of the pixel values in an annulus of 10–15 pixels, centered on the cluster. In F814W, most of the northern tail was on the WF4 chip. These were measured using apertures with 2-pixel radii, or 3-pixel apertures for the poorly resolved groups, with sky subtraction determined by the mode in a 5–8 pixel annulus. These sizes were chosen because the WF chips have approximately twice the angular size of the PC chip, so the same area of sky was measured regardless of which chip a cluster fell on. Very little of the northern tail was visible in F555W.

The nucleus was always on the PC chip, but required the use of smaller apertures with radii of 2 or 3 pixels due to crowding. Background levels were determined by the mode in a 10–15 pixel annulus.

The bridge was generally on the WF3 chip. In F380W, the bridge unfortunately straddled the WF3 and WF4 chips, so a few clusters could not be observed in this filter. The aperture radius was 3 pixels, with background determined by the median in an 8–13 pixel annulus. This aperture size was selected because bridge clusters are often isolated, but there are several areas where multiple emission sources are packed closely together. In such cases, daofind generally detected only one or two of the brightest sources, which were always sufficiently separated to be measured individually.

Aperture corrections could not be determined on the images. It was therefore necessary to use the standard corrections prescribed by Holtzman et al. (1995). The corrections for WF3 were used for all of the wide field chips. Standard corrections are not given for the F380W or F606W filters, so corrections for the F336W and F555W, respectively, were used. The corrections were typically  mag.

4.4 H ii region photometry

In addition to the five regions A–E, we measured the nucleus (designated region F) and four regions (G–I) in the northern part of the bridge, shown in H and Spitzer in Fig. 4. The GALEX and Spitzer images were registered to the HST data so that age estimates could be obtained using the total flux in the HST bands. This also allowed identification of the clusters making up the H ii regions. Most of region C was off frame in all of the HST bands except F814W. Region G could not be measured at because of a cosmic ray.

The aperture radius was selected as (5 pixels on the H image) in all instruments to ensure that the same area of sky was represented at all wavelengths. The sky annulus widths were in HST, in H, in GALEX, and in Spitzer. The measurements were not aperture-corrected because the HST apertures were largely empty, and the emission regions were highly pixelated in GALEX and Spitzer.

Table 4 summarizes the photometric measurements. The H and R bands were calibrated to standard units using the results of James et al. (2004). The ratio of the H flux to R flux density provided an estimate of the equivalent width of the H line, EW(H). This equivalent width is included in the sb99 models and provides another method for determining the ages of the emission regions. Ages are presented in Table 5 and extinctions in Table 6. These will be discussed in § 6.3.

5 Hst Results

5.1 Cluster colours

Colour-colour diagrams are shown in Fig. 5. In the top left panel, we plot the magnitude in the F555W filter minus the magnitude in the F606W filter (F555WF606W) versus (F606WF814W) for all detected clusters for which these colours were available. The curves show sb99 instantaneous star formation models with metallicity for (solid curve), 0.24 (dashed curve), and 0.5 (dot-dashed curve) for ages from yr. The youngest ages are at the lowest values of (F606WF814W). The models are shown to give a sense of the effect of reddening, and it should be noted that is not the observed metallicity for the nucleus and parts of the southern tail. The bridge population is generally bluer in (F555WF606W) and redder in (F606WF814W) than the populations of the other regions. The southern tail has the opposite tendency, leaning red in (F555WF606W) and blue in (F606WF814W). Only a few clusters in the northern tail were on frame in F555W; those that were tend toward more intermediate values in both colours. The nuclear population is similar.

The top right panel shows all detected clusters in (F380WF606W) versus (F606WF814W), along with sb99 models. The youngest ages on the curves appear at the lowest values of (F606WF814W). The middle and bottom panels show the same colours, with the different regions plotted separately and error bars shown. The differences between the populations are less striking in this colour combination than the one discussed above, and arise primarily in (F606WF814W).

The median statistical uncertainties in the northern tail are 0.081 and 0.095 mag for (F380WF606W) and (F606WF814W), respectively. In the southern tail, the median uncertainties are 0.153 and 0.107 mag; in the nucleus, 0.050 and 0.040 mag; and in the bridge, 0.167 and 0.123 mag. The uncertainties are large in the southern tail in part because it includes clumps detected at the level. If only level clumps are included, the median uncertainties are 0.114 in (F380WF606W) and 0.088 in (F606WF814W).

5.2 Age determination

Age and reddening estimates were obtained using the starburst99 v5.1 evolutionary synthesis code (Leitherer et al., 1999; Vázquez & Leitherer, 2005). We modeled the colours of the clusters assuming a single, instantaneous burst of star formation with a Kroupa (2002) initial mass function over the range 0.1–. Spectral energy distributions (SEDs) were calculated for cluster ages of 1 Myr to 10 Gyr and included contributions from H emission. The SEDs were reddened from 0–2 mag of using the reddening law of Calzetti et al. (1994). These model spectra were convolved with the bandpasses of the four HST filters. The model age and extinction providing the best fit to the observations was determined using minimization

Uncertainties in both the ages and extinctions were determined as in Smith et al. (2008), and are driven largely by the degeneracy between age and reddening. Breaking this degeneracy is particularly difficult for clusters with coverage in only three filters. Such clusters generally have larger age uncertainties than those with that were covered in all four filters.

Models were run for three different metallicities: , , and . The model is representative of the spectroscopically measured abundances in the nucleus and H ii region A, while the model is close to that of H ii regions B and C. The solar composition models were used primarily for comparison.

Several colour combinations were used because not all clusters were available in all colours. Age estimates obtained from sb99 have been found to be most reliable when (1) large numbers of colours are used and (2) the bluest available combination of colours are used. U coverage is particularly important (e.g. Hancock et al., 2008). We adhered to these results in determining which of the available colour combinations would provide the best age estimate for each cluster at each metallicity.

Our best-fitting ages for the clusters are given in Table 2, while and masses are shown in Table 3. The age distributions within the features of the system are shown in Fig. 6. Finding charts showing the locations of the clusters are given in Figs. 712, with circles indicating ages Myr, diamonds 11–22 Myr, and squares Myr. The apertures used for photometry of the large H ii regions are also shown. The clusters are generally very young, occupying the two younger age bins (1–10 and 11–22 Myr). Very few clusters older than 22 Myr are found in the tidal tails.

Ages determined by evolutionary synthesis techniques generally have large uncertainties. In our sample, the uncertainties are typically very close to the ages. The median age is 14 Myr, with a median negative uncertainty of 5 Myr and a positive uncertainty of 10 Myr.

5.3 Luminosity function

The luminosities in the F814W band for the clusters were calculated using a distance of 39 Mpc (NED), using a bandwidth of Hz. The luminosity distribution, excluding the extremely luminous nuclear population, is shown in Fig. 13, and provides one estimate of our sample completeness. We assume a simple power law form for the luminosity function, defined by . A best-fitting line is determined after the distribution turns over, giving a slope of with a completeness limit of .

The uncertainty in the slope may be estimated as follows. The distribution has a flat top, so there are two reasonable possibilities for the turnover point. Selecting as the completeness limit , we obtain a slope of slope . Our best estimate for the slope is . Cluster luminosity functions typically have (e.g. Whitmore et al., 1999), which gives a slope of when binned logarithmically. Our luminosity function is consistent with this trend to within the uncertainties.

6 Discussion

6.1 Hst cluster ages

The metallicities determined by González Delgado et al. (1995) aid in the interpretation of our evolutionary synthesis results. In the southern tail, H ii region A has a measured metallicity . The best age estimates should therefore come from our models, which show this region to be dominated by clusters with ages 10 Myr or less.

The metallicity of H ii region C at the base of the southern tail is (González Delgado et al., 1995). Combined with the general metallicity gradient observed in the galaxy, we are led to favor the model in the western part of the southern tail, which shows more clusters older than 22 Myr than does the model. Note, however, that part of region C was off frame (Fig. 11), so we only have age estimates for two of its constituent clusters.

In the northern tail, the dependence of derived age on metallicity is less pronounced. Only five clusters change age bins between the and models, all of which are near the nucleus. Most of the clusters are in H ii regions B and D. These regions have measured abundances , so our models should provide the best age estimates. The cluster ages are all less than 23 Myr in these models, with about half 10 Myr or younger. This is not unexpected, since this part of the northern tail is detected in H (Smith et al., 1997) and Spitzer (Smith et al., 2007).

The metallicity of the nucleus has been measured at , so the model is preferred here. Seven of the 16 clusters have ages  Myr, with the remainder between 14–21 Myr. The gap between these age populations is probably not significant, since the rapid evolution of red giants after 10 Myr makes exact age determination exceptionally difficult and often results in gaps between 10–15 Myr (e.g. Fall et al., 2005). The younger population fits well with the UV spectral results of González Delgado et al. (1999), while the older population fits with one of the bursts expected by Lançon et al. (2001).

Of the 124 clusters detected in the tidal tails of NGC 7714, those with ages  Myr number only five in the model and 12 in the model. If an age is assigned to each cluster based on the measured or most likely metallicity, there are eight older than 22 Myr. The small sample of 16 clusters detected in the nucleus has no clusters older than 22 Myr. There is little evidence in NGC 7714 of a population older than  Myr, despite an interaction that had closest approach  Myr ago. We also find that the cluster age groups are evenly distributed over the galaxy. The lack of older clusters detected in the nucleus may be due to the difficulty in detecting such clusters amidst the high background caused by YMCs.

These results are similar to those obtained by Bastian et al. (2005b) for the interacting galaxy NGC 6872. Most of the clusters in NGC 6872 were found by colour fitting to be younger than 30 Myr. These were spread throughout the galaxy, extending far out into the tidal tails. Clusters older than 30 Myr were concentrated in the central regions of the galaxy. However, the ‘central regions’ of NGC 6872 have a spatial size  kpc, while the imaged portion of the NGC 7714 disc in our data is only  kpc. The age distribution within NGC 7714 is compatible with that found by Bastian et al. in NGC 6872, but a wider field will be required to determine if the concentration of old clusters is similar in the two galaxies.

The bridge ages are difficult to interpret because of the large colour uncertainties, incomplete coverage in F380W, and lack of metallicity data. The metallicity in NGC 7714 generally decreases with distance from the nucleus. However, some of the bridge material may originate in NGC 7715, which has unknown metallicity. Nevertheless, a low metallicity is expected due to the abundances measured elsewhere in the system and because tidal features in interacting galaxies generally show low abundances.

Despite the large uncertainties, we note that the only evidence for a significant population with age Myr is in the bridge. Out of 34 clusters, the numbers with ages  Myr are 12 and 17 in the and models, respectively. However, only one of these clusters has F380W data available in the case, and nine in the case. The existence of an older population in the bridge is by no means certain, but if present could indicate that cluster dissolution mechanisms have acted more weakly here than in the tidal tails, which in this system are closer to the nucleus. Older clusters should survive more easily in the bridge because tidal forces are generally smaller there than in the tails. However, the gas expulsion mechanisms thought to be responsible for infant mortality operate independently from the tidal forces, so the rate at which very young clusters dissolve is not expected to change.

To address the issue of sample completeness, we begin with the luminosity function (Fig. 13). Cutting out clusters with luminosities below the completeness limit, we have 49 left in the tidal tails and only four in the bridge. The bridge sample is far too incomplete to permit any reliable analysis.

Only two clusters in the tidal tails that survive the luminosity cut have ages Myr, and those only in the model. It is unsurprising that older clusters would be preferentially cut from a luminosity-limited sample, since clusters get both redder and fainter as their stellar populations age.

6.2 Hst cluster mass function

It is common practice in studies of cluster ages to use a mass cutoff rather than a luminosity cutoff. Due to evolutionary fading, clusters of a given age will remain above the luminosity limit only if they have a sufficiently high mass. Therefore, imposing a luminosity limit will generally under-sample old clusters.

Mass estimates are based on the assumption that clusters of a given age all have the same mass to light ratio, with luminosity values calculated for clusters with mass . The luminosity of a cluster relative to a cluster of the same age then gives an estimate of the mass. The uncertainties in the masses are quite large, up to a factor . Mass estimates based on the best-fitting age and extinction are shown in Table 3.

The lowest mass for which the sample is complete was determined assuming a simple power law form for the cluster mass distribution. For , we obtain a reasonably good fit with a slope and a limiting mass of . The quality of the fit is not high due to a lack of clusters in the mass bin centered at , which may reflect the high uncertainties in masses determined in this way. The mass estimates for the model produced a slope of , but with a lower mass limit at about . The mass distribution for the model is shown in Fig. 14.

The uncertainties for the mass function fits may be estimated using the same technique as was used for the luminosity function. For the model, the next best choice of limiting mass would be , which gives a slope of of . In the model, both and are reasonable choices, producing slopes of and , respectively. Using these alternate choices to estimate the uncertainties, the model has a slope , with smaller values being more likely. For the model, the slope is . At the limits of the uncertainties, our mass functions agree with the typical power law index near , found for example in the Antennae (Zhang & Fall, 1999) and Magellanic Clouds (Hunter et al., 2003). Due to the small size of our sample, we are unable to test the claim of Gieles (2009) that the cluster initial mass function is truncated at the high mass end.

The completeness limit in mass, unlike that in luminosity, is a function of age due to evolutionary fading. Thus the turn-over in Fig. 14 gives an age-averaged completeness limit for mass. For younger systems, the mass completeness limit is lower than for older systems. The masses are plotted against the ages for the model in Fig. 15, with the NGC 7714 disc clusters as blue circles and bridge clusters as red asterisks. The solid black curve shows the mass of a cluster with an F840W luminosity equal to our completeness limit of , as a function of age. At lower ages, we are reaching less massive clusters, and are more deficient at older ages. For clusters less than 10 Myr old, we are complete to about , while for yr we are complete only to . The dashed black line indicates the average completeness limit of . The significance of evolutionary fading is apparent, as all of the clusters with ages yr have masses .

Mass cuts alter the age distributions considerably. In the model, seven of the 27 clusters remaining after eliminating those with mass are older than 22 Myr, with the youngest of these 111 Myr old. These account for six of the seven clusters remaining in the bridge (only one of which has F380W coverage) and one of six in the northern tail, but none of the 14 clusters in the southern tail.

Of the 63 clusters remaining after eliminating those with mass in the model, 31 are in the southern tail, 18 in the northern tail, and 14 in the bridge. There are 22 clusters with ages Myr, 12 of them in the bridge, of which five lacked coverage in F380W. In the southern tail, seven of 31 clusters are older than 22 Myr (five of them detected at level), while the northern tail has two of 18 clusters older than 22 Myr. The youngest of these clusters has an age of 89 Myr.

The number of clusters falling into each age bin in the model is summarized in Table 7 for all measured clusters, clusters above the luminosity limit, and clusters above the average mass completeness limit. For comparison, we have included columns showing the expected numbers of clusters assuming a cluster formation rate of 1 cluster Myr, and constant cluster formation for the last (CCF) and (CCF) yr, only including clusters above the mass limits. These are limiting cases, based on Fig. 14. The observed number of clusters in the mass limit case fall between these two limits. Thus with the available data we cannot rule out constant cluster formation with no infant mortality. More sensitive data are needed to reduce the completeness limits and search for evidence of cluster destruction. With the current data set, the number of clusters is relatively small and the mass cut-off relatively high, thus the results on infant mortality are uncertain.

6.3 H ii region age estimates

Ages for the H ii regions were estimated in two ways. The colours obtained from large-aperture photometry in the GALEX and HST bands were fit to sb99 models, using the methods described in § 5.2. In regions E and G–J, the FUV and NUV measurements provided only upper limits, so these bands were not fit to the models. The range of cluster ages within the H ii regions suggests that the star formation histories are more complex than the instantaneous starbursts used to model the individual clusters. We therefore used both instantaneous and continuous star formation models to estimate the ages of these regions. The true star formation history is likely to fall between these two extremes.

Ages were also estimated using the H equivalent widths, which were fit to instantaneous burst sb99 models. These data are shown in Table 5, along with the median age of the HST clusters detected in each region. Note that the HST sample has no clusters in region E (four were detected but were covered in only two filters), only two clusters in region C, since it is at the edge of the WFPC2 field of view, and in the bridge (regions G–J) only one or two clusters were found in each region, so this indicator is only useful in regions A, B, D, and F. When possible, we also compare with published ages determined via photoionisation modeling.

The ages determined using EW(H) are generally consistent with the median age of clusters found in the region, with both showing cluster population ages Myr. These estimates are somewhat higher than those in the literature (see Table 5), but are within a factor of in most of the large aperture photometry cases.

The ages determined by fitting sb99 models to the large aperture photometry of the H ii regions are much higher than those determined by other methods in regions B, D, and F. The area within the apertures occupied by clusters is relatively small, so it is possible that faint sources outside the clusters are making significant contributions to the total luminosity of the H ii regions. If these sources are older and redder, the age of the region as a whole could be quite different from that of its bright clusters. These sources could be older clusters below our detection threshold or the remnants of clusters that have already dissolved and entered the field population.

To quantify the possible influence of faint sources on the HST apertures, the fraction of luminosity in the F814W band for each of the H ii regions originating in the measured clusters was determined. The values are similar throughout the NGC 7714 disc, with per cent of the light coming from outside the clusters. Unresolved objects could significantly impact the colours, which may explain the higher age estimates for regions B, D, and F.

All of the bridge H ii regions had their lowest age estimates provided by colour fitting, while varies. Regions G and I had the lowest values in the system, while H and J have the highest. Visual inspection of regions G and I indicates several likely clusters that were not selected by daofind. Considering the relatively high noise and low resolution in HST bridge data, it is likely that these results are products of the data quality. Deeper imaging will be required to learn if the bridge hosts a cluster population that differs from that of the NGC 7714 disc.

Extinction estimates are shown in Table 6, along with some values from the literature. The extinction estimates of González Delgado et al. (1995) come from the Balmer decrement, while the extinction in the nucleus found by González Delgado et al. (1999) was determined using evolutionary synthesis results for UV continuum flux distributions, and was smaller than that found using the Balmer decrement by a factor . Our extinction estimate for the nucleus is in good agreement with that of González Delgado et al. (1999), while our other estimates are higher than the literature values by a factor . This is not surprising, since the Balmer decrement tends to give higher extinction estimates than the continuum by a factor (Fanelli et al., 1988; Storchi-Bergmann et al., 1994), possibly because dust in star forming regions with strong nebular emission is destroyed or removed by ionizing radiation, stellar winds, or supernovae (Fanelli et al., 1988). Another possibility is that the hot stars responsible for the Balmer lines are still associated with the dusty molecular clouds in which they formed, while the cooler stars that contribute to the continuum have drifted away from these regions and suffer less extinction (Calzetti et al., 1994).

6.4 Spitzer colours

Figs. 16 and 17 show as red circles the Spitzer [4.5][5.8] vs. [3.6][4.5] and [5.8][8.0] vs. [4.5][5.8] colours of the H ii regions. These were calculated from the flux densities as described in the IRAC Data Handbook. For comparison, we also show emission ‘clumps’ from several other interacting systems, including Arp 24 (Cao & We, 2007), 82 (Hancock et al., 2007), 107 (Smith et al., 2005b), 285 (Smith et al., 2008) and NGC 2207/IC 2163 (Elmegreen et al., 2006), as green crosses. In the Arp 285 system, NGC 2856 tail clump 3 and disc clump 1 were unusual, and are marked in Fig. 16 as a magenta open diamond and cyan limit, respectively. The unusually luminous clump ‘i’ in the eastern tail of NGC 2207 is also shown separately, as an open blue diamond. We also display the mean colours of field stars of Whitney et al. (2004) (magenta open triangle), M0III stars (M. Cohen 2005, private communication; open blue square), quasars (Hatziminaoglou et al. 2005; red squares), and diffuse galactic dust (Flagey et al. 2006; blue crosses).

Most of the H ii region Spitzer colours fit well with the clumps measured in other systems. This is important, since the selection criteria are somewhat different. The Arp 284 regions were selected based on H emission, whereas clumps in the other systems were selected based on Spitzer emission. While the H emission usually had a counterpart in Spitzer , in some cases the aperture placement would have differed slightly if the emission had been used for selection.

[4.5][5.8] is often used as a gauge of the star formation rate, with redder colours indicating higher rates of mass-normalized star formation (e.g. Smith et al., 2005b). The H ii regions lie between interstellar matter and stars, indicating that, as in other interacting systems, they probably have contributions from both. However, none of the Arp 284 regions is quite as red as the Arp 285 tail clump. Smith et al. (2008) suggest that the exceptionally red colour of this clump is due to strong emission from interstellar matter owing to the very young age of the clump, which was determined to be 4 Myr using optical colours. Region J has the same estimated age and is only slightly bluer in colour, suggesting that it may also have numerous hot young stars.

The [5.8][8.0] colours generally fall close to those of interstellar matter, which is not surprising since the 5.8 and 8.0 bands are expected to be dominated by interstellar dust.

A few of the regions are unusual. Regions B and C are unusually red in [3.6][4.5], but slightly blue in the other two colours. Region F is redder in [3.6][4.5] than clumps in other systems with similar [4.5][5.8] colour. Region H is slightly blue in [5.8][8.0].

The colours of regions B and C suggest excess emission at , perhaps with a slight deficiency at 5.8 and , qualitatively similar to low-metallicity dwarf galaxies. In Figs. 18 and 19, we compare the Spitzer colours of Arp 284 H ii regions to those of the dwarf and spiral galaxies of Smith & Hancock (2009). Regions B and C are not as red in [3.6][4.5] as some the dwarfs, but are generally redder than the spirals. This may indicate that they are deficient in PAHs, though not extremely so, consistent with their moderately low metallicities.

It is interesting that regions B and C stand out in these plots, but the other regions do not. Our colour-based age estimate for region B is very similar to that of D, and our estimates based on EW(H) vary little between regions. García-Vargas et al. (1997), however, found regions B and C to be extremely young. This would imply more nebular continuum and Br emission, which are two of the factors that appear to contribute to reddening [3.6][4.5] in the dwarfs (Smith & Hancock, 2009).

The SED for region B is plotted in Fig. 20. We also show as black curves the best-fitting evolutionary synthesis models from sb99 for both continuous (top panel) and instantaneous burst (bottom panel) star formation models. The curves use the best-fitting ages for each model, which are Myr for continuous star formation and Myr for the instantaneous star formation model. The models include both stellar and nebular emission. The solid green lines show the same models without nebular emission. The models corresponding to 1 uncertainty (68 per cent confidence) are also shown. All of the curves have been scaled to match the observations in F606W. The stellar and stellar+nebular emission models are nearly indistinguishable, indicating that nebular emission is relatively unimportant for these inferred ages. The large mid-IR excess could result from hot dust continuum emission, PAH lines, or an older stellar population, either in clusters that have faded below our detection limits or in the field. The excess is diminished in the continuous star formation model, probably because it includes older stars.

For comparison, the SED of the nucleus is shown in Fig. 21. The top panel shows curves for a continuous star formation model with an age of Myr, while the bottom panel shows an instantaneous starburst with an age of Myr. In this plot we have also included the near-IR J, H, and K band fluxes from Lançon et al. (2001). The nucleus is known to have a complex star formation history, likely including stars older than  Myr. The mid-IR excess is greater in the instantaneous burst model, likely due older stars not included in this model.

6.5 Comparison with Chandra

6.5.1 Ultra-luminous X-ray sources (ULXs)

X-ray observations by Smith et al. (2005a) detected three non-nuclear point sources with erg s, the standard definition of an ultraluminous X-ray (ULX) source. Three other point sources lie just below this limit. The nature of ULXs is not presently understood. In one possible scenario, they could be intermediate mass () black holes (IMBH; Colbert & Mushotzky 1999). Alternatively, they may be stellar mass black holes with beamed (King et al., 2001) or super-Eddington (Begelman, 2002) X-ray emission. IMBHs may preferentially form in massive star clusters (e.g. Portegies Zwart et al., 2004a). It takes at least 4 Myr for a stellar mass black hole/X-ray binary to form, and they are expected to be common in populations with ages up to Myr (Rappaport et al., 2004). In contrast, an IMBH may form in a YMC in Myr (Portegies Zwart et al., 2004a). Evolutionary models of the predicted X-ray emission from IMBH+stellar companion mass-transfer binaries indicate that such ULXs should generally have ages less than about 15–30 Myr (Portegies Zwart et al., 2004a).

As discussed in Smith et al. (2005a), the ULX candidates in NGC 7714/5 exist in a range of environments. A luminous X-ray point source is located away from the nucleus, near cluster n01, which has an estimated age of 14 Myr. If this association is correct, the young age is more consistent with the IMBH scenario. Another possibly young ULX candidate is found within H ii region C, which has an age of 11 Myr based on EW(H) and a spectroscopically determined age of 4.5 Myr (García-Vargas et al., 1997). However, this second ULX cannot be associated with a specific star cluster, since there are several very close together in the area. Finally, one ULX is located in the outer western tail of NGC 7714, in an area that is undetected in H but is bright in GALEX images. Three other ULX candidates are neither located in H ii regions nor are they close to optically-selected star clusters. Given the variety of environments, from this small sample we cannot make any strong conclusions about the nature of ULXs as a class.

6.5.2 Extended X-ray emission associated with star forming regions

In addition to point sources, the Chandra map revealed extended X-ray emission associated with some of the H ii regions (Smith et al., 2005a). Fig. 2 shows Chandra 0.3–8 keV X-ray contours over a colour scale H image of NGC 7714. Extended X-ray emission is found associated with all of the disc H ii regions except D. The bridge regions G–J also do not have observed diffuse X-ray emission.

The possible contributors to the extended X-ray emission include multiple unresolved high mass X-ray binaries (HMXBs), or hot gas from supernovae and stellar winds impacting the ambient interstellar medium. Based on Local Group studies, the HMXB component is expected to be too small to account for this emission in the NGC 7714/5 regions (Smith et al., 2005a), unless there is an excess above Local Group galaxies. For the hot gas component, using sb99 models and a constant X-ray production efficiency /, where is the mechanical luminosity from the supernovae and winds, the ratio of the number of Lyman continuum photons to X-ray luminosity / is expected to drop off with age for a star forming region (Smith et al., 2005a). However, the one H ii region complex in the NGC 7714 disc that is not detected in extended X-ray emission, region D, has an EW(H)-based age that is somewhat older than some of the other regions (Table 5), in spite of having the highest /. In contrast, region E is bright in the X-rays, with a low /, but has a younger age from its EW(H). This suggests that the X-ray production efficiency may vary with time and/or from region to region, and/or the HMXB contribution may be larger than expected in some regions. This result is uncertain, however, due to the relatively large uncertainties on our age estimates and the small number of H ii regions in our sample.

7 Conclusions

We have used broadband colours to obtain age estimates for 174 proto-globular cluster candidates in the Arp 284 system. The populations detected in the tidal features of NGC 7714 are generally quite young ( Myr old). However, when the sample is limited by a mass cut, the dominance of young clusters becomes unclear in the NGC 7714 disc and vanishes in the bridge. Due to the small number of clusters above the mass cut, this latter result is highly uncertain. Deeper imaging is required to better define the cluster mass function and push down the completeness limit, particularly in the bridge, which has higher photometric uncertainties and less F380W coverage than the rest of the system.

No substantial differences are found between the populations within NGC 7714. The clusters in the nucleus do not generally differ in age from those in the tidal structures. It therefore appears that the early-stage dissolution processes, if present, operate via internal processes, rather than local environmental effects, but we cannot rule out the possibility of weaker dissolution in the bridge.

Using larger apertures to study H ii region complexes, we find in several cases that the region as a whole appears to be older than the clusters detected inside it. This suggests an older, redder population that is unresolved, possibly the remnants of star clusters that have already dissolved or which were not detected as discrete sources because of evolutionary fading.

The large aperture results suggest that star formation has been ongoing in H ii regions B and D, located in the northern tail, for a dynamical time-scale, like that in the nucleus. In both cases, the star formation history, whether continuous or in discrete bursts, is unknown. We also do not know if regions B and D are gravitationally bound, or simply regions where successive waves or infall events have triggered star formation. Region C is similar, but we have insufficient data to comment on its nature. If in the future these regions are found to be bound, they could be interpreted as nearby analogs to the high redshift disc clumps studied by Elmegreen et al. (2009).

The Arp 284 interaction has been in progress for over 100 Myr, yet we find little evidence for a significant cluster population older than 20 Myr. Thus it is possible that clusters are dissolving on a relatively short time-scale. However, due to completeness issues, with the available data we cannot definitively confirm the possibility of infant mortality in this system. Further, more sensitive observations are needed to search for a deficiency of older clusters.


We thank the HST, Spitzer, GALEX, and Chandra teams for making this research possible. We also thank Nate Bastian and the anonymous referee for helpful comments and suggestions. We acknowledge support from NASA Spitzer grant 1347980 and NASA Chandra grant AR90010B. BJS acknowledges support from NASA LTSA grant NAG5-13079. 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.


  • Arp (1966) Arp H., 1966. Atlas of Peculiar Galaxies (Pasadena: Caltech)
  • Bastian & Goodwin (2006) Bastian N., Goodwin S. P., 2006, MNRAS, 369, L9
  • Bastian et al. (2005a) Bastian N., Gieles M., Lamers H. J. G. L. M., Scheepmaker R. A., de Grijs R., 2005, A&A, 431, 905
  • Bastian et al. (2005b) Bastian N., Hempel M., Kissler-Patig M., Homeier N. L., Trancho G., 2005, A&A, 435, 65
  • Bastian et al. (2009) Bastian N., Trancho G., Konstantopoulos I. S., Miller B. W., ApJ, 2009, 701, 607
  • Begelman (2002) Begelman M., 2002, ApJ, 568, L97
  • Bernlöhr (1993) Bernlöhr K., 1993, A&A, 258, 25
  • Brandl et al. (2004) Brandl B. R., Devost D., Higdon S. J. U. et al., 2004, ApJS, 154, 188
  • Calzetti et al. (1994) Calzetti D., Kinney A. L., Storchi-Bergmann T., 1994, ApJ, 429, 582
  • Cao & We (2007) Cao C., Wu H., 2007, AJ, 133, 1710
  • Colbert & Mushotzky (1999) Colbert E., Mushotzky R., 1999, ApJ, 519, 89
  • Conti & Vacca (1994) Conti P. S., Vacca W. D., 1994, ApJ, 423, L97
  • de Grijs & Parmentier (2007) de Gijs R., Parmentier G., 2007, ChJAA, 2, 155
  • Elmegreen et al. (2009) Elmegreen B. G., Elmegreen D. M., Fernandez M. X., Lemoias J. J., 2009, ApJ, 692, 12
  • Elmegreen et al. (2006) Elmegreen D. M., Elmegreen B. G., Kaufman M., Sheth K., Struck C., Thomasson M., Brinks E., 2006, ApJ, 642, 158
  • Fall et al. (2005) Fall S. M., Chandar R., Whitmore B. C., 2005, ApJ, 631, L133
  • Fanelli et al. (1988) Fanelli M. N., O’Connell R. W., Thuan T. X., 1988, ApJ, 334, 665
  • Fazio et al. (2004) Fazio G. G., Hora J. L., Allen L. E., 2004, ApJS, 154, 10
  • Flagey et al. (2006) Flagey N., Boulanger F., Vrestraete L., Miville Deschênes M. A., Noriega Crespo A., Reach W. T., 2006, A&A, 453, 969
  • García-Vargas et al. (1997) García-Vargas M. L., González Delgado R. M., Pérez E., Alloin D., Díaz A., Terlevich E., 1997, ApJ, 478, 112
  • Gieles (2009) Gieles M., 2009, MNRAS, 394, 2126
  • Gieles et al. (2006) Gieles M., Portegies Zwart S. F., Baumgardt H., Athanassoula E., Lamers H. J. G. L. M., Sipior M., Leenaarts J., 2006, MNRAS, 371, 793
  • González Delgado et al. (1995) González Delgado R. M., Pérez E., Díaz Á. I., García-Vargas M. L., Terlevich E., Vilchez J. M., 1995, ApJ, 439, 604
  • González Delgado et al. (1999) González Delgado R. M., García-Vargas M. L., Goldader J., Leitherer C., Pasquali A., 1999, ApJ, 513, 707
  • Goodwin & Bastian (2006) Goodwin S. P., Bastian N., 2006, MNRAS, 373, 752
  • Hancock et al. (2007) Hancock M., Smith B. J., Struck C., Giroux M. L., Appleton P. N., Charmandaris V., Reach W. T., 2007, AJ, 133, 676
  • Hancock et al. (2008) Hancock M., Smith B. J., Giroux M. L., Struck C., 2008, MNRAS, 289, 1470
  • Hatziminaoglou et al. (2005) Hatziminaoglou E., Pérez-Fournon I., Polletta M. et al., 2005, AJ, 129, 1198
  • Holtzman et al. (1992) Holtzman J. A., Faber S. M., Shaya E. J. et al., 1992, AJ, 103, 691
  • Holtzman et al. (1995) Holtzman J. A., Hester J. J., Casertano S. et al., 1995, PASP, 107, 156
  • Hunter et al. (1994) Hunter D. A., O’Connell R. W., Gallagher J. S., III, 1994, AJ, 108, 84
  • Hunter et al. (2003) Hunter D. A., Elmegreen B. G., Dupuy T. J., Mortonson M., 2003, AJ, 126, 1836
  • James et al. (2004) James P. A., Shane N. S., Beckman J. E. et al., 2004, A&A, 414, 23
  • Kennicutt (1984) Kennicutt R. C., Jr., 1984, ApJ, 287, 116
  • King et al. (2001) King A. R., Davies M. B., Ward M. J., Fabbiano G., Elvis M., 2001, ApJ, 552, L109
  • Kroupa (2002) Kroupa P., 2002, Sci, 295, 85
  • Lançon et al. (2001) Lançon A., Goldader J. D., Leitherer C., González Delgado R. M., 2001, ApJ, 552, 150
  • Leitherer et al. (1999) Leitherer C., Schaerer D., Goldader J. D. et al., 1999, ApJS, 123, 3
  • Martin et al. (2005) Martin D. C., Fanson J., Schiminovich D. et al., 2005, ApJ, 619, L1
  • O’Halloran et al. (2000) O’Halloran B., Metcalfe L., Delaney M. et al., 2000, MNRAS, 360, 871
  • Portegies Zwart et al. (2004a) Portegies Zwart S. F., Dewi J., Maccarone T., 2004a, MNRAS, 355, 413
  • Portegies Zwart et al. (2004a) Portegies Zwart S. F., Baumgardt H., Hut P., Makino J., McMillan S. L. W., 2004b, Nat, 428, 724
  • Rappaport et al. (2004) Rappaport S. A., Fregeau J. M., Spruit H., 2004, ApJ, 606, 436
  • Smith & Hancock (2009) Smith B. J., Hancock M., 2009, AJ, 138, 130
  • Smith & Wallin (1992) Smith B. J., Wallin J. F., 1992, ApJ, 393, 544
  • Smith et al. (1997) Smith B. J., Struck C., Pogge R. W., 1997, ApJ, 283, 754
  • Smith et al. (2005a) Smith B. J., Struck C., Nowak M. A., 2005a, AJ, 129, 1350
  • Smith et al. (2005b) Smith B. J., Struck C., Appleton P. N., Charmandaris V., Reach W., Eitter J. J., 2005b, AJ, 130, 2117
  • Smith et al. (2007) Smith B. J., Struck C., Hancock M., Appleton P. N., Charmandaris V., Reach W. T., 2007, AJ, 133, 791
  • Smith et al. (2008) Smith B. J., Struck C., Hancock M. et al., 2008, AJ, 135, 2406
  • Storchi-Bergmann et al. (1994) Storchi-Bergmann S., Calzetti D., Kinney A. L., 1994, ApJ, 429, 572
  • Struck & Smith (2003) Struck C., Smith B. J., 2003, ApJ, 289, 157
  • Tremonti et al. (2001) Tremonti C. A., Calzetti D., Leitherer C., Heckman T. M., 2001, ApJ, 555, 322
  • Vázquez & Leitherer (2005) Vázquez G. A., Leitherer C., 2005, ApJ, 621, 695
  • Weedman et al. (1981) Weedman D. W., Feldman F. R., Balzano V. A., Ramsey L. W., Sramek R. A., Wu C. -C., 1981, ApJ, 248, 105
  • Weisskopf et al. (2002) Weisskopf M. C., Brinkman B., Canizares C., Garmire G., Murray, S., Van Speybroeck L. P., 2002, PASP, 114, 1
  • Werner et al. (2002) Werner M. W., Roellig T. L., Low F. J. et al., 2004, ApJS, 154, 1
  • Whitmore et al. (1993) Whitmore B. C., Schwiezer F., Leitherer C., Borne K., Robert C., 1993, AJ, 106, 1354
  • Whitmore et al. (1999) Whitmore B. C., Zhang Q., Leitherer C., Fall S. M., Schweizer F., Miller B. W., 1999, AJ, 118, 1551
  • Whitney et al. (2004) Whitney B. A., Indebetouw R., Babler B. L. et al., 2004, ApJS, 154, 315
  • Zhang & Fall (1999) Zhang Q., Fall S. M., 1999, ApJ, 527, L81
Figure 1: Optical photograph of Arp 284 from Arp (1966), with numbers labeling the prominent tidal features. North is up and east to the left. The field of view is .
Figure 2: NGC 7714 disc indicating H ii regions from Smith et al. (2005a). The colour scale shows H. The contours are a Chandra 0.3–8 keV x-ray map.
Figure 3: Smoothed GALEX FUV (colour) overlaid with VLA H i (contours). Note that the western tail loops back to the main galaxy, and this ring-like feature is present in both H i and UV. The field of view is .
Figure 4: H ii regions in H (top) and Spitzer (bottom). The aperture radii have an angular size of , and are also shown on the HST images in Figs. 712. The field of view is .
Figure 5: In the top panels, red circles represent clusters from the bridge, black diamonds the southern tail, blue asterisks the northern tail, and black squares the nucleus. The curves show sb99 instantaneous star formation models with metallicity for (solid curve), 0.24 (dashed curve), and 0.5 (dot-dashed curve). These curves give a sense of the effect of reddening, but the metallicity is not representative of the nucleus or parts of the southern tail. The middle and bottom four panels show the same colours as the upper right panel separated by region, with error bars.
Figure 6: Histograms of cluster ages from models in the southern tail (top left panel), northern tail (top right panel), and bridge (bottom left panel), and from the model in the nucleus (bottom right panel).
Figure 7: F606W image of NGC 7714 with northern, central, and southern areas indicated. These areas are examined more closely in Figs. 812. Several H ii regions identified in lower-resolution data are also indicated. Note that region C is partially off frame toward the bottom right of the box indicating the southern area.
Figure 8: F606W image of the northern tail of NGC 7714, with ages indicated for all detected clusters for the models. Circles represent ages 1–10 Myr, diamonds 11–22 Myr, squares Myr. The metallicities of H ii regions B and D have been measured at , while the nucleus has metallicity . All clusters shown are considered part of the northern tail.
Figure 9: F606W image of the NGC 7714 nucleus with ages from the model indicated for all detected clusters. Circles represent ages 1–10 Myr, diamonds 11–22 Myr, squares Myr. The metallicity has been measured at .
Figure 10: F606W image of the central region of NGC 7714 with ages from the model indicated for all detected clusters. Circles represent ages 1–10 Myr, diamonds 11–22 Myr, squares Myr. The nucleus has a measured metallicity of , so this model should give the best age estimates. These clusters are considered to be part of the southern tail.
Figure 11: F606W image of the southern tail of NGC 7714 with ages from the models indicated for all detected clusters. Circles represent ages 1–10 Myr, diamonds 11–22 Myr, squares Myr. Clusters 10 and 15 are inside H ii region C, which has a measured metallicity of
Figure 12: F555W image with all detected clumps in the bridge with ages from the models. Circles represent ages 1–10 Myr, diamonds 11–22 Myr, squares Myr. The metallicity of the bridge has not been measured, but is likely to be low.
Figure 13: Luminosity distribution in F814W for all detected clusters outside the nucleus. The luminosity is defined as the luminosity in the F814W band , where is the bandwidth of Hz. The solid red line is a power law fit to the distribution after the turnover with slope , giving a completeness limit of approximately .
Figure 14: Cluster mass distributions for the model. The solid red line is a power law fit to the part of the distribution after the turnover. The fit has slope with a limiting mass .
Figure 15: Mass vs. Age for the model. Blue circles represent the NGC 7714 disc, while red asterisks are bridge clusters. The dashed black line represents the age-averaged mass completeness limit of from Fig. 14. The solid black curve shows the mass of a cluster with an F840W luminosity equal to our completeness limit of , as a function of age.
Figure 16: Spitzer [4.5][5.8] vs. [3.6][4.5] colour-colour plot. The Arp 284 H ii regions (open red circles) are shown along with emission ‘clumps’ from Arp 24 (Cao & We, 2007), 82 (Hancock et al., 2007), 107 (Smith et al., 2005b), 285 (Smith et al., 2008), and NGC 2207/IC 2163 (Elmegreen et al. 2006; green crosses). Shown separately are Arp 285 NGC 2856 tail clump 3 (magenta open diamond) and disc clump 1 (cyan limit), as well as NGC 2207/IC 2163 ‘clump i’ shown separately (open blue diamond). Also shown are the mean colours of field stars of Whitney et al. (2004) (magenta open triangle), M0III stars (M. Cohen 2005, private communication; open blue square), quasars (Hatziminaoglou et al. 2005; red squares), and diffuse dust in the Milky Way (Flagey et al. 2006; blue crosses).
Figure 17: Same as Fig. 16, but for [5.8][8.0] vs. [4.5][5.8].
Figure 18: Spitzer [4.5][5.8] vs. [3.6][4.5] colour-colour plot showing Arp 284 H ii regions (open red circles). The dwarf galaxies of Smith & Hancock (2009) are shown by metallicity, with (red open triangles), (green filled circles), and (blue open squares), along with a sample of ‘normal’ spirals (black crosses).
Figure 19: Same as Fig. 18, but for [5.8][8.0] vs. [4.5][5.8].
Figure 20: The FUV-mid-IR SED of region B. The black curves are the best-fitting evolutionary synthesis models including both stellar and nebular emission. The model in the top panel is for continuous star formation with an age of Myr, while the bottom panel shows an instantaneous burst with an age of Myr. The solid green line plotted from 3–9.5 shows the models without nebular emission. The short-dashed red and dot-dashed blue curves correspond to the youngest and oldest models and the associated best-fitting extinctions. The models have been scaled to match the observations in F606W.
Figure 21: Same as Fig. 20, but shows region F with models, with ages of Myr for the continuous star formation model and Myr for the instantaneous burst model.
Dataset Exposure Filter H ii Regions H ii Regions
(sec) in PC chip in WF chip Comment
U3GQ0201R 500 F380W nucleus,A,B,part C,D,E bridge bridge split WF3/4
U3GQ0202R 500 F380W nucleus,A,B,part C,D,E bridge bridge split WF3/4
U3GQ0203R 400 F380W nucleus,A,B,part C,D,E bridge CR on PC chip
bridge split WF3/4
U3GQ0204M 400 F380W nucleus,A,B,part C,D,E bridge bridge split WF3/4
U6A01601R 350 F555W nucleus,A,E bridge
U6A01602R 350 F555W nucleus,A,E bridge
U2E68801T 500 F606W nucleus,A,B,part C,D,E bridge
U6A01701R 350 F814W nucleus,A B,C,part D
U6A01702R 350 F814W nucleus,A B,C,part D
U6A01801R 350 F814W nucleus,A, E bridge
U6A01802R 350 F814W nucleus,A, E bridge CR on PC chip
Table 1: HST WFPC2 imaging of NGC 7714/5
ID F380W F555W F606W F814W Age () Age () Age ()
mag mag mag mag Myr Myr Myr
1 22.460.06 22.640.03 22.180.05 15 15 15
2 23.300.14 22.840.06 23.100.15 4 5 5
3 22.330.07 21.230.03 21.670.07 1 2 3
4 22.460.08 22.020.04 22.090.08 6 6 7
5 21.250.04 20.860.02 21.060.04 5 5 6
6 22.180.06 22.040.03 22.280.09 5 5 6
7 21.000.03 21.300.02 21.070.03 5 14 15
8 22.800.11 21.780.04 21.890.05 3 4 4
9 20.840.03 20.820.04 20.740.07 5 14 15
10 21.130.04 21.430.03 21.070.05 15 14 15
11 22.890.11 22.270.05 21.940.13 32 21 20
12 23.390.17 22.970.09 22.710.16 35 20 20
13 21.410.05 21.510.03 21.170.07 6 14 15
14 22.300.10 21.620.07 21.890.22 3 4 5
15 21.890.07 21.730.04 21.040.06 15 21 20
Table 2: Cluster magnitudes and age estimates. The full table is available in the electronic edition of the journal.
1 0.029 0.029 0.029
2 0.012 0.012 0.013
3 0.047 0.053 0.061
4 0.013 0.013 0.013
5 0.042 0.074 0.065
6 0.011 0.017 0.017
7 0.081 0.074 0.051
8 0.075 0.086 0.067
9 0.110 0.100 0.070
10 0.081 0.074 0.081
11 0.063 0.067 0.074
12 0.024 0.024 0.039
13 0.074 0.067 0.031
14 0.053 0.043 0.031
15 0.145 0.142 0.083
Table 3: Cluster and mass estimates. The full table is available in the electronic edition of the journal.
Table 4: H ii region photometry
Region FUV222FUV, NUV, F380W, F555W, F606W, F814W in units of erg s cm Å NUV F380W F555W F606W F814W H333H in units of erg s cm 444, , , in units of mJy
A 208095 189081
B 1130110 47463
C 155080 55140
D 50091 23658
E 25 21.2
F 13600100 6220110
G 11.2 8.3
H 10.4 6.5
I 7.3 5.9
J 6.3 5.7
Table 5: H ii region ages (Myr)
Region Age (inst)555Ages in Myr from large-aperture broadband photometry for instantaneous starburst (inst) and continuous star formation (cont) models Age (inst) Age (cont) Age (cont) Med. Age666Median age of HST-resolved clusters in the region. These are useful only for regions A, B, D, F; see text for explanation. Med Age EW(H) EW(H) Age EW(H) Age Lit. Age777Ages determined by photoionization modeling. A, B, C from García-Vargas et al. (1997); F from González Delgado et al. (1995) and González Delgado et al. (1999) 888Fraction of H ii region luminosity contained in sample clusters in F814W 999From Smith et al. (2005a)
A 22
B 11
C 6.1
E 3.1
F 1.5
Table 6: H ii region extinctions
Region (inst) (inst) (cont) (cont) Lit. 101010A,B, C from González Delgado et al. (1995); F from González Delgado et al. (1995) (Balmer decrement) and González Delgado et al. (1999) (UV continuum)
Table 7: Clusters by age bin for
Age All Luminosity Cut Mass Cut CCF111111Expected numbers for constant cluster formation over the last 30 Myr at a rate of 1 cluster Myr CCF121212As above for 100 Myr
1–10 Myr 59 23 13 10 10
11–22 Myr 70 28 28 12 12
23–100 Myr 2 0 1 8 78
Myr 27 2 21 0 0
Total 158 53 63 30 100
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description