arXiv:1705.01251v1 [astro-ph.GA] 3 May 2017 yaiso prlglxe.Oe h oreo several of course the Over and galaxies. al. evolution spiral et secular of Mazzuca the there- dynamics understand are 2006; to rings (e.g., paramount al. nuclear galaxy fore et of trigger studies Knapen host Comprehensive the 1995; consequently the 2011). dominate al. of may et often enhanced Genzel history which ring the un- star-formation episodes, the entire and the starburst in timescale, violent 2010), rate crossing al. Comer´oncollision et short 2008; usually al. et Mazzuca 1996; 2012b). Combes al. & et Buta Kim galactic 2004; 1992; the Maciejewski around Athanassoula ring (e.g., compact orbits center circular luminous, nearly a in forms ac- move the and ‘contact to where two proceeds ring, the the gas and at cumulated lanes momen- region inflow- dust 1995; angular ring the The between the their al. points’ toward et of in 2012a). Knapen most spiral al. lose tum, 1994; et which al. Kim materials, et 2003; ing Byrd Teuben (e.g., & edges to- Regan bar leading the infall at the gas lanes dust of causes torque the along potential gravitational nucleus the bar The ward asymmetric the 2006). 2006; from al. et Moiseev NGC Mazzuca & Combes 2004; and Sil’chenko 1994; Kennicutt Athanassoula 7217 & 1990; Kormendy NGC 2001; al. (e.g., can et Shlosman galaxies rings 7742; unbarred nuclear some enhancements generate arms—in density also spiral shaped strong bar- that oval/elliptically and of reported structures— spiral as non-axisymmetric products been barred such other has intense natural and in it events the interest- to merger although found are are dynamics, subject mostly they induced nuclei, are often since They galactic galaxies, are of that starbursts. kpc substructures 1 ing within found MRVDMTO ODTRIETEITGAE RPRISO N OF PROPERTIES INTEGRATED THE DETERMINE TO METHOD IMPROVED hngacn a inDsrc,Biig109,China 100190, Beijing District, Hai Dian 5, Hai Lu Zhongguancun, Yuan [email protected] China; He 100871, Yi Beijing University, District, Peking Dian Astronomy, of ment h ihgsdniy(.. ide l 1992; al. et Wild (e.g., density gas high The usually are which rings, starburst Circumnuclear cetdfrpbiaini h srpyia ora Suppl Journal L Astrophysical using The typeset Preprint in publication for Accepted 2 1 nentoa pc cec nttt–ejn,1Nanerti 1 Institute–Beijing, Science Space International Depart- and Astrophysics & Astronomy for Institute Kavli ujc headings: Subject trfrainrt SR ftencerrn ntebre spiral barred spectra the the in ring measure rin nuclear accurately the the to of of resolution evolution method (SFR) and new formation rate a the star-formation with present as we well Here as dynamics, its and o naeaeaeof age average rin an for the obtain To results. the of robustness the ∼ on effect negative F ufc est,Σ density, surface SFR este nnra aais eas siaetern’ oa tla m stellar total ring’s the estimate also We galaxies. normal in densities –0Mr eapyormto otecniumsbrce H continuum-subtracted the to method our apply we Myr, 3–10 h nertdpoete fncerrnsaecreae ihthe with correlated are rings nuclear of properties integrated The 1. Hubble A INTRODUCTION T E tl mltajv 12/16/11 v. emulateapj style X ehd:osrainl–glxe:eouin–glxe:idvda ( individual galaxies: – evolution galaxies: – observational methods: and aais udmna aaees–glxe:sa formation star galaxies: – parameters fundamental galaxies: ∼ cetdfrpbiaini h srpyia ora Supp Journal Astrophysical The in publication for Accepted pte pc Telescope Space Spitzer SFR 0Myr. 40 hoMa Chao =0 . 09 M 1 ihr eGrijs de Richard , ⊙ yr − 1 kpc ABSTRACT ao, − mgs mg erdto osnthv significant a have not does degradation Image images. 2 hc smc ihrta h ikaeae SFR disk-averaged the than higher much is which , ements where r ihcmlmnayifae esrmnss sto as so measurements trac- infrared SFR complementary optical with or ers circum- to ultraviolet To indispensable of 2002). is observations al. it combine problem, et value dust-attenuation Benedict the extinction 2000; vent al. by cluster et effect ring Buta this (e.g., representative the minimize a of to although adopting tried because ring, have the 1998) investigators uncertain- across Kennicutt some extinction large (e.g., variable in 50% of to result presence up could appli- of equation authors, ties imaging. this these the by of from recognized cation derived readily directly as usually However, is which ity, follow usually 1999; 2008) al. et al. Buta et Mazzuca (e.g., (1998): 2001; SFR Kennicutt al. ring et global Maoz to approaches the Traditional calculate integrated their (SFRs). the and rates derive luminosities popu- total to star-formation their cluster as importance such young great parameters, the of ring prop- is of ring it evolution integrated lation, the the explore to affects (or is environment erties) purpose ring ultimate fea- the our important clus- how Since star as the serve population. to thus Kennicutt ter properties They & galaxy-scale Kormendy 2015). connecting 2012; tures 2001; al. et al. param- Li et structural 2004; Weiner galaxy’s Anders al. (e.g., host & et 2017). Grijs the eters Maoz al. de constrain et to Grijs 1995; nostics de 2012; 2013; al. al. et et al. disk Barth found Laan et der nearby are Hsieh van (e.g., clus- clusters in star 1996; star regions abundance massive young prime young in rapid of are where between presence rings galaxies connection the nuclear close al. Crocker and et ters, the & Laan Buta formation to der van star Owing 1987; 2000; of Shlosman al. 2013). & et knowledge Englmaier Hummel 1960; 1993; our 1961; Burbidge & studies to Burbidge Sandage (e.g., observational properties significantly ring and nuclear contributed theoretical have many decades, 1,2 h rpriso ula ig r logo diag- good also are rings nuclear of properties The n usC Ho C. Luis and , L F (M SFR (H α ) obs α aayNC11 ae nhigh- on based 1512 NGC galaxy nrydsrbto n current and distribution energy l ’ F o h eidspanning period the for SFR g’s rhs aayssclrevolution secular galaxy’s host ir n 8 and ⊙ ersnsteH the represents yr s,lg( log ass, ’ trcutrpopulation(s). cluster star g’s − lements 1 1 µ 7 = ) mgs h resulting The images. m CERRNS G 1512 NGC RINGS: UCLEAR M/M . 9 × 10 ⊙ α 7.1 = ) − msinln luminos- emission-line G 52 – 1512) NGC 42 [ L (H ± α ) 0.11 obs ] , (1) 2 Chao Ma, Richard de Grijs, and Luis C. Ho account for missing flux caused by the dust absorption a final, common science image of 585 × 585 pixels2 (or and scattering (for a review, see e.g. Kennicutt & Evans 1400 × 1400pc2), as shown in Fig. 1. 2012, and references therein). We used Kennicutt et al. (2009)’s calibration of the lo- Mazzuca et al. (2008) have shown that nuclear rings cal SFR: are coplanar with their disks and retain nearly circular −1 × −42 shapes after deprojection. Since these structures are em- SFR (M⊙ yr )=5.5 10 [L(Hα)obs +0.011L(8 µm)], (2) bedded in the galactic disk and bulge of their host galaxy, −1 the total, intrinsic ring luminosities measured are un- where L(Hα)obs is expressed in units of erg s prior to avoidably contaminated by contributions from the host any internal dust-attenuation correction, and L(8 µm) is the Spitzer 8 µm polycyclic aromatic hydrocarbon galaxy. However, few studies have taken background cor- −1 rections into account when calculating the total ring flux, (PAH) emission-line luminosity, also in erg s . This except for Maoz et al. (1996), who measured the total calibration was derived assuming a Kroupa initial mass ring flux by integrating all counts above the background function (IMF) with stellar masses in the range 0.1–100 over the entire ring region. Note that their background M⊙ (Kroupa & Weidner 2003) and solar chemical abun- measure is, in fact, a constant value, determined in an dance; Allard et al. (2006) and Sarzi et al. (2007) showed ‘empty’ corner of the image. This thus only provides a that circumnuclear regions of barred galaxies can be simple approximation to the often highly complex back- modeled well by adopting solar metallicity. Hα emission ground in the ring region. Consequently, the main mo- is dominated by young massive stars and is commonly tivation for this paper is to demonstrate the feasibility used as an instantaneous SFR indicator, tracing stars ∼ of a newly devised method to obtain clean ring luminosi- with lifetimes of 3–10 Myr and masses greater than ties based on galaxy light-profile fitting using the galfit 40 M⊙ (e.g., Kennicutt et al. 2009; Hao et al. 2011). algorithm (Peng et al. 2002, 2010). Since HST observations of the NGC 1512 circumnu- As a test case, we will obtain the spectral energy dis- clear ring in Hα line emission are also available in the tribution (SED) and SFR of the luminous nuclear ring HLA, we additionally retrieved a narrow-band F658N in the galaxy NGC 1512, while simultaneously allowing (Hα) image (GO-6738; PI: Filippenko), observed with for inherent dust attenuation and background correction. the WFPC2/Planetary Camera (PC). The PC chip has ∼ ′′ Our analysis is based on multi-waveband Hubble Space a pixel size of 0.05 . The pixel scale of the Hα im- Telescope (HST) imaging data, combined with observa- age was resampled to match that of the WFC3 images ∼ ′′ tions obtained with the Infrared Array Camera (IRAC) ( 0.04 ) and then aligned to the F555W frame, using on board the Spitzer Space Telescope. NGC 1512 is lo- standard iraf/stsdas5 routines. Observational infor- cated at high Galactic latitude and it is, hence, barely mation for all HST images used in this paper can be affected by Galactic foreground extinction, E(B − V ) = found in Table 1. 0.011 mag (Schlegel et al. 1998); we will hence ignore The aligned Hα band is actually a combined emission any correction for Galactic foreground reddening. The line plus continuum image. To remove the underlying remainder of this paper is organized as follows. In Sec- stellar continuum, we used the scaled F555W and F814W tion 2, we describe the observational data used as well as images on both the short- and long-wavelength sides of our data reduction approach. Section 3 presents in detail the F658N filter bandpass to linearly interpolate a con- our improved procedure for calculating the ring flux, and tinuum image at the reference Hα wavelength. Since addresses a series of technical issues. Finally, we provide the WFC3 and WFPC2/PC detectors have almost the ′′ a brief summary and our main conclusions in Section 4. same PSF size, FWHM ∼ 0.07 (i.e., they have compa- 2. rable spatial resolutions), we did not apply PSF match- OBSERVATIONS AND DATA REDUCTION ing. Next, the continuum image thus generated was sub- Broad-band HST/Wide Field Camera-3 (WFC3; pixel tracted from the F658N image, after converting both im- ′′ size ∼ 0.04 ) images of NGC 1512 were acquired from ages from counts s−1 to the physical flux units of erg the HST Legacy Archive (HLA)3, observed as part of s−1 cm−2 A˚−1 using the calibration parameter phot- program GO-13364 (PI: Calzetti). Our multi-wavelength flam (i.e., the inverse sensitivity, contained in the im- imaging data set included observations through the age header), also listed in Table 1. This provides us F336W (roughly corresponding to the Johnson–Cousins with a reliable continuum subtraction without generat- U band), F438W (B), F555W (V ), and F814W (I) fil- ing any conspicuous negative-flux regions, as shown in ters, which were pipeline-processed and calibrated using left-bottom panel of Fig. 1. As a consistency check, the standard HLA reduction software. They were already we compared our continuum-subtracted Hα image to aligned to the same orientation and field of view, covering Fig. 1 of Maoz et al. (2001). The latter was gener- the inner galactic region (see Fig. 1). In this paper we ated based on their WFPC2/F547M and F814W im- − ± assume a distance modulus of (m M)0 = 30.48 0.25 ages. We changed the display colors, inverted the as- mag4, which corresponds to a distance of 13 Mpc. This sociated color map, and adjusted the scale limits so as to results in a pixel scale of 2.4 pc pixel−1. Since the match the appearance of Maoz et al.’s image as closely circumnuclear starburst ring represents a small fraction as possible. Although both images are indeed very sim- of the full image, for convenience we cropped the ring- ilar, our new, higher-resolution image shows additional dominated area in the original four exposures to yield details in the form of dust lanes. This will, however,
3 http://hla.stsci.edu/hlaview.html 5 The Image Redcution and Analysis Facility (iraf) is dis- 4 This value, computed from values collected in the NASA tributed by the National Optical Astronomy Observatories, which Extragalactic Database (NED; http://nedwww.ipac.caltech.edu/), are operated by the Association of Universities for Research in As- presents the geometric mean of 10 individual distance measure- tronomy, Inc., under cooperative agreement with the US National ments from the literature. Science Foundation. Improved method to determine the integrated properties of nuclear rings 3
Fig. 2.— Radial profiles derived from the F336W image. Red dots represent the observed data and blue dots are for the best- fitting model. Left: Best-fitting model using three galactic com- ponents (i.e., background and nuclear ring). Right: Application of only a model of the background.
Fig. 1.— Top left: Color composite image of our target 585 × to focus on the inner ring area only. 585 pixels2 nuclear ring area in NGC 1512, created by stacking im- ages taken through the F336W, F555W, and F814W filters. The 3. MEASURING THE INTEGRATED RING LUMINOSITIES F336W filter, which traces young stellar populations, is shown in blue, while F555W and F814W are shown in green and red, respec- To characterize the ring structure, we employed gal- tively. The F555W and F814W images are shown in the middle and fit to simultaneously fit the two-dimensional stellar light right-hand panels, respectively. Bottom: Continuum-subtracted distributions of multiple galactic components using ana- Hα/[Nii] (left), F814W (middle), and continuum-subtracted 8 µm lytical functions. It generates a final model image based (right) images. All panels share the same orientation and field of view as indicated in the top left-hand panel. on the best-fitting structural parameters. The opti- mal solution for the model parameters is obtained us- ing the iterative Levenberg–Marquardt algorithm (e.g., not cause any problems in measuring the ring flux, since Bevington & Robinson 2003). Since the standard statis- the faint continuum structure forms part of the back- tical uncertainties returned by galfit tend to underesti- ground flux, which must be subtracted from the image mate the true uncertainties, we evaluated the goodness- 2 (see below). We prefer to use the WFC3/F555W fil- of-fit based on both the χ values given by galfit and ter rather than WFPC2/F547M, because the more re- visual inspection of the residual images. cently installed WFC3 camera is characterized by bet- ter spatial resolution and improved efficiency. Moreover, 3.1. HST image analysis the F555W image reaches fainter photoemtric limits and With regard to optical HST images, we have developed hence it shows more detailed structures compared with an improved approach to derive the total ring luminos- F547M. Note that our continuum-subtracted Hα image ity, based on inspection of the radial surface brightness includes adjacent [Nii] lines redshifted to λλ6568, 6604A.˚ profiles. To start with, we performed a test using the We adopted a redshift of 0.0029 (Koribalski et al. 2004). F336W image, which contains the most luminous nuclear However, the possible effects of dust extinction in the ring component among our optical HST bands. Hα filter significantly undermine its reliability for mea- We plotted the radial surface brightness profile in the suring the SFR. Dust particles in the interstellar medium F336W image, shown as the red dotted line in the left- absorb a fraction of the Hα flux; the attenuated Hα lu- hand panel of Fig. 2. The presence of the nuclear ring minosity is subsequently re-radiated at infrared wave- introduces a prominent peak at radii between 118 and lengths. To more robustly estimate the dust-corrected 250 pixels, corresponding to a radial range from 4.6′′ to SFRs, a better solution involves combining Hα with in- 9.8′′. Therefore, we adopted this annulus to define the frared observations. We used the Spitzer/IRAC post- ring region in this paper. Based on the shape of the basic calibrated data (‘post-BCD’) images (3.6 and 8 observed radial profile, we assume that our observed im- µm bands, with a pixel scale of 0.6′′; see Table 1) age consists of the host galaxy’s background contribution taken from the Spitzer Heritage Archive and observed and a nuclear ring, and we intend to model all galactic as part of the Spitzer Infrared Nearby Galaxies Survey components in the original image. The nuclear ring can (SINGS; Kennicutt et al. 2003). To derive the pure 8 be quantified by a truncation function, and the back- µm PAH emission, needed as input in Eq. (2), we had ground contribution in the image can be fitted by adopt- to remove the stellar continuum contribution from the ing a mixture of Gaussian and S´ersic functions. This 8 µm image. Since the 3.6 µm band traces stellar mass approach yields a good fit with minimal residuals. The in nearby galaxies (e.g., Elmegreen & Elmegreen 1984; truncation function is, in essence, a hyperbolic tangent Eskew et al. 2012; Meidt et al. 2012), a continuum-free function. In galfit these functions are commonly used PAH emission image at 8 µm can be obtained by sub- to modify the light profiles of galactic components, such tracting a scaled 3.6 µm band image from the original as S´ersic profiles, to produce ring-like structures (for a 8 µm image, following the recipe of Helou et al. (2004) detailed description, see Peng et al. 2010). and employing a scale factor of 0.255. This approach To limit the degrees of freedom in the fits, we (only) results in highly accurate continuum-subtracted images, fixed the center position; the other parameter were left which has been widely demonstrated and used in the lit- as free parameters. We found a good match between the erature (e.g., Calzetti et al. 2007; Kennicutt et al. 2009; best-fitting model and the real data. The middle panel Zhou et al. 2015). The bottom right-hand panel of Fig. of Fig. 3 presents our best-fitting model based on the im- 1 shows the continuum-subtracted 8 µm image cropped proved method, and the relevant radial profile is shown 4 Chao Ma, Richard de Grijs, and Luis C. Ho
TABLE 1 Observational details of the final NGC 1512 image data set.
Filter Proposal ID/PI Camera Exposure time (s) photflam Zmag F336W (U) GO-13364/Calzetti WFC3 1107 1.3168305 × 10−18 23.6012 F438W (B) GO-13364/Calzetti WFC3 953 6.91138715 × 10−19 24.3011 HST F555W (V ) GO-13364/Calzetti WFC3 1131 1.876528 × 10−19 25.7166 F658N(Hα) GO-6738/Filippenko WFPC2/PC 5200 1.45442135 × 1017 20.993 F814W(I) GO-13364/Calzetti WFC3 977 1.5304799 × 10−19 25.9379 Band Program ID/PI Camera Exposure time (s) FLUXCONV Zmag Spitzer 3.6 µm SINGS1-159/Kennicutt IRAC 193 0.1069 18.8024 8 µm SINGS1-159/Kennicutt IRAC 176 0.2026 17.1984 Notes. Keword photflam is given in units of erg s−1 cm−2 A˚−1. The instrumental magnitudes have been converted to the STmag magnitude system by applying the photometric zero-points (Zmag), whose calibration relies on the photflam and photzpt keywords, i.e., Zmag = −2.5 log(photflam) + photzpt (where photzpt= −21.10 for the ST magnitude scale). The Spitzer Zmags were taken directly from the image headers.
Fig. 3.— Left: Observed F336W image. Middle: Best-fitting Fig. 4.— Left: Original F336W image. Middle: As the left-hand three-component model for the F336W image. Right: Image of the image, but for a pixel scale that has been degraded to that of the individual truncation component based on the best-fitting model infrared image. Right: Final image after degradation of its pixel in the middle panel. scale and Gaussian smoothing, so as to mimic the observed infrared data. in the left-hand panel of Fig. 2 (blue dotted line; the red line represents the actual data). It is clear that a although those authors did not actually define a clear- three-component model could entirely recover the obser- cut ring region. Next, we corrected for the contributions vational profile with high accuracy. Using this method, of the [Nii] lines using the theoretical [Nii] flux line ra- the total ring luminosity can be derived by adding the tio (1/3; Osterbrock & Ferland 2006) and the observed flux in all pixels of the predefined ring region: see the [Nii]λ6584/Hα flux line ratio of NGC 1512 obtained from right-hand panel of Fig. 3, which shows the truncation Calzetti et al. (2007). L([NII]) was removed as a func- tion of transmission efficiency of the F658N filter at the component (designed to model the nuclear ring). We 39 . position of the lines. We obtained L (Hα)=4.74×10 determined m (U)=14.83+0 05 mag. obs ring −0.05 erg s−1, which will be used to calculate the ring’s SFR, We also investigated how the use of different PSF im- together with L(8 µm), discussed in the next subsection. ages affected our fit results. We ran galfit using an empirical PSF, a theoretical PSF, and without any input 3.2. Spitzer image analysis PSF, respectively. We verified that the effect of chang- ing the PSF is negligible, because the three conditions A caveat regarding the Spitzer/IRAC image is that the employed accurately reproduced the model image with IRAC image resolution for a typical PSF FWHM ∼ 2′′ the same structural parameters and the same resulting is much lower than that of the HST PSF. As illustrated ring luminosity. This can be readily understood given the in the bottom right-hand panel of Fig. 1, such a resolu- significant difference between the ring and the PSF sizes: tion is comparable to the radius of the ring, and hence the typical WFC3/UVIS stellar FWHM of ∼ 0.08′′(e.g., we cannot ascertain whether or not the measured ring Calzetti et al. 2015) is much smaller than the ring radius luminosity would be affected by the lower Spitzer res- ∼ 8′′. To find the best fits faster, we opted to not provide olution. To clarify this issue, our strategy consisted of the PSF when dealing with our HST images. artificially degrading the HST/WFC3 F336W image to We subsequently applied our method to the other HST the resolution of Spitzer/IRAC so as to simulate the con- filters. The measured total ring magnitudes in the other ditions pertaining to low resolution. We subsequently filters in our data set are listed in Table 2. Note that no ran galfit on the degraded F336W image to verify attempt was made to obtain a definitive model for the whether we could recover the same result as that de- bulge. Our main goal was just to measure the proper, rived before image degradation. To construct the sim- localized ring structure, with the remaining parts of the ulated image, we first changed its pixel size to that of image together considered the image’s (i.e., the ring’s) Spitzer (∼ 0.6′′/pixel), as shown in the middle panel of background. Fig. 4. It was then broadened by convolution with a To determine the observed ring’s Hα luminosity, we ap- Gaussian kernel to match the smoothness of the Spitzer plied our method to the continuum-subtracted Hα+[Nii] image, as shown in the right-hand panel of Fig. 4. We image produced in the previous section. We found found that the ring magnitude derived from our simu- × 39 −1 +0.19 L(Hα + [NII]) = 5.61 10 erg s , with a fractional lated image is mring(U) = 14.83−0.17 mag, which is in uncertainty of ∼10%. This is very close to the lumi- very good agreement with that obtained prior to image nosity (6×1039 erg s−1) derived by Maoz et al. (2001), degradation, supporting the application of our method Improved method to determine the integrated properties of nuclear rings 5
10.5 11.0 ] 2 11.5 8μm
12.0
12.5 13.0