arXiv:1810.09018v2 [astro-ph.GA] 5 Nov 2018 h iko prlglxe see.g. (see galaxies spiral of disk the .A Zinchenko A. I. da CALIFA on based relations H parameter extragalactic nebular of regions–II: stars ionizing of temperature Effective atinzn asv stars, massive ionizing tant MNRAS Zanstra 2015 al. et Evans therein). 2017 2011 a topee ( atmospheres lar oie yteesas yuigti ehdlg,proposed methodology, nebu this the using by By emitted stars. emission-lines these of by analysis ionized the from e.g. cetd21 ot 0 eevd21 uut0;i original in 00; August 2018 Received 00. Month 2018 Accepted h eemnto fteefcietmeaue( temperature effective the of determination The INTRODUCTION 1 n trs eogn oH to belonging star(s) ing eim(S)a ela h aayi hc hyreside. they which in Int galaxy the the as of well excitation as the (ISM) atmosphere, Medium stellar mass of of physics evolution the and formation the restrict that processes ietyetmtdb sn hi htmti n spectrom and photometric their (e.g. using data ca by temperature effective estimated the directly Clouds, Magellanic the and Way 5 4 3 2 1 c autdd inisAto´mcsyGeof´ısicas, Univers Astron´omicas y Ciencias de Arge Facultad (CONICET-UNLP), Plata C La Astrof´ısica 291, de de Hifumi, Instituto Shishima Av. Paraiba, do Vale do Universidade f¨ur Astronomie Scien Zentrum of Rechen-Institut, Astronomisches Academy National of Observatory Astronomical Main 00TeAuthors The 0000 ; ; u oefcso pct n/rln-lneigi h stel the in line-blanketing and/or opacity of effects to Due o oiigsaso eryH nearby of stars ionizing For akv ta.2018 al. et Markova orl ta.2014 al. et Morrell ( 000 1929 asye l 2005 al. et Massey 0–0 00)Pern oebr21 oplduigMNR using Compiled 2018 November 6 Preprint (0000) 000–000 , ,i spsil oestimate to possible is it ), ; bot&Hme 1985 Hummer & Abbott orSihe l 2017 al. et Mohr-Smith 1 ; .Hwvr o h aoiyo h dis- the of majority the for However, ). , abr ta.2014 al. et Walborn 2 II T .L Dors L. O. , , ein scuilt nesadthe understand to crucial is regions e ff 2009 a nyb nietyestimated, indirectly be only can ose l 2017 al. et Dors II ; rmteCLF aarlae3(R)adtertcl(for theoretical and (DR3) 3 release data CALIFA the from f5 prlglxe nte0 the in galaxies spiral 59 of hs (O phase ecluaeteefcietmeaue( temperature effective the calculate We ABSTRACT rdet ftoeprmtr nec aatcds o h o the ( for sample disk our galactic in each galaxies in parameters consi We those emission-lines. of strong gradients on based calibrations O/H) xii eaieo a ns h einvleo the of value median The ones. flat or negative exhibit aaiswt a xgnaudnegainstn ohv fla have to tend gradients abundance oxygen flat with galaxies n htrda rdet fbt log both of gradients radial that find e words: Key galaxies, spiral of objects. majority these the for of property disk the along gradients nicreainbtentegain of gradient the between anti-correlation h es httegain flog of gradient the that sense the el lhuhorrslsaei gemn ihteie ft of idea the with agreement in are results our Although well. ein,lctdi h Milky the in located regions, ot ta.2007 al. et Corti T e ff n t eaiu along behaviour its and ; ; carr&Schmutz & Schaerer / ; atn Palacios & Martins ) n h oiainparameter ionization the and H), 3 abe l 2016 al. et Lamb T .F H F. G. , e ff aais bnacs–IM bnacs–H – abundances ISM: – abundances galaxies: n references and fteioniz- the of ) ddNcoa eL lt,PsodlBsu /,10 aPlata La 1900 s/n, Bosque del Paseo Plata, La de Nacional idad H 69120 M¨onchhofstr. Universit¨at Heidelberg, 12–14, der ; ntina oae al. et Sota om21 uut00 August 2018 form p24-0,SoJs o aps P Brazil SP, Campos, dos Jose Sao ep12244-000, e fUrie 7Zbltoosr,060Ke,Ukraine Kiev, 03680 str., Zabolotnogo 27 Ukraine, of ces v stars, ive erstellar etric be n agele ¨ ∼ lae by - ; 0% hw positive shows %) 70 . 005 4 , U rdet(e also (see gradient asbtlower but mass opeetlwrvle of values lower present to rpsdby proposed n of ent Init the of (e.g. limit Function mass upper stellar universal an with formed grad metallicity exhibit that galaxies sense spiral the that known well is ( 2003 Copetti & Dors ( 1994 s nms ftesia aais hywr o on nearl in found not were they galaxies, spiral the (e.g. of studies most in ist H of sample H the of n e hoeia airto ewe h bevdemis observed the between ratio calibration line theoretical new a ing H for distance tric ais(e.g. radius 5 < 2017 1978 .V Cardaci V. M. , nrae slog(O as increases z U ; tde the studied ) width equivalent the of enhancement the that interpreted ) < atn ta.2005 al. et Martins T and T e 0 ff T e β R ff . e ntedsso prlglxe.I fact, In galaxies. spiral of disks the in 3rdhf ag.W s pcrsoi aataken data spectroscopic use We range. redshift 03 =log([O U ff finzn trs,oye bnac ftegas the of abundance oxygen star(s), ionizing of ) T iygne l 2004 al. et Pilyugin msinln ihteglcoeti itnefra for distance galactocentric the with emission-line iroe l 1986 al. et Fierro e os&Cpti2003 Copetti & Dors II n h xgnaudnegain.Moreover, gradient. abundance oxygen the and ff o apeo H of sample a for Z ata ta.2010 al. et Bastian Z ein nM11cudb u oapositive a to due be could 101 M in regions eedo h xgnaudnegain,in gradient, abundance oxygen the on depend erae ihteiceeto h galactocentric the of increment the with decreases T (e.g. ice ae 1988 Pagel & Vilchez e II II T ff , ]( e ff ailgainsee huhsm them some though even gradients radial ein eogn o1 prlglxe us- galaxies spiral 14 to belonging regions 2005 λλ T / oime l 2004 al. et Mokiem aito safnto ftegalactocen- the of function a as variation )gain nrae hl hr san is there while increases gradient H) e ff 3726+29)/[O ,saswt ihrmtliiy( metallicity higher with stars ), T 4 .I pt fteegainssol ex- should gradients these of spite In ). ailgain s072kK/ 0.762 is gradient radial e , ff 5 II hssesntt ea universal an be to not seems this e pta itiuinadradial and distribution spatial der jcsi u ape oto the of Most sample. our in bjects .C Krabbe C. A. , hntercutrat ihtesame the with counterparts their than ; .Teeoe suigta tr are stars that assuming Therefore, ). iebr,Germany eidelberg, regions vn 1986 Evans T ,i sepce oiiegradi- positive a expected is it ), e II log t .Teeatosfudpositive found authors These ). eeitneo positive of existence he ff ein oae ntedisks the in located regions and Argentina , III U ; ] U λ er oad1995 Howard & Henry ; and 07 and 5007) II .Recently, ). atn ta.2004 al. et Martins n miia (for empirical and ) SL AS T A hed Searle & Shields e T ff E tl l v3.0 file style X ta rdet as gradients 3 T e ose al. et Dors ff R (relation Z a Mass ial et,in ients, 25 trend ) sion- We . T .It ). T e e ff y ff ; 2 I. A. Zinchenko et al.
Integral Field Area Survey1 (CALIFA, S´anchez et al. 2012), which is ideal to investigate global scaling relations between galaxy prop- erties (see e.g. Ellison et al. 2018). The main goals of the present study are to investigate if Teff gradients are universal properties of spiral galaxies, and the existence of any correlation between Teff and nebular parameters such as the ionization parameter or the oxy- gen abundance. This paper is organized as follows. The methodol- ogy assumed to calculate Teff and the observational data used along this work are described in Sec. 2. In Sect. 3 the results and discus- sion of the outcome are presented. Finally, conclusions are given in − Sect. 4.
− 2 METHODOLOGY 2.1 Sample − − − − We used publicly available spectra from the integral field spec- α troscopic CALIFA survey data release 3 (DR3; S´anchez et al. 2016, 2012; Walcher et al. 2014) based on observations with the Figure 1. BPT diagnostic diagram (Baldwin et al. 1981) for all the spaxels PMAS/PPAK integral field spectrophotometer mounted on the in our sample. The solid line is the boundary between AGN-like and H II- Calar Alto 3.5-meter telescope. CALIFA DR3 provides wide-field like regions defined by Kewley et al. (2001). IFU data for 667 objects in total. The data for each galaxy consist of two datacubes, which cover the spectral regions of 4300–7000 Å at a spectral resolution of R ∼ 850 (setup V500) and of 3700– 5000 Å at R ∼ 1650 (setup V1200). For the galaxies with both V500 and V1200 datacubes available, there are COMB datacubes for 446 galaxies which are a combination of V500 and V1200 dat- gradients for 11 of these galaxies, null gradients for two and a neg- acubes covering the 3700–7000 Å spectral range. In this study we ative gradient for the other one (see also P´erez-Montero & V´ılchez used these COMB datacubes. 2009). In particular, the first negative Teff gradient was found for The sample of galaxies is described in detail in the Milky Way by Morisset (2004). Additional analysis taking into Zinchenko et al. (2018, in prep.). Briefly, we selected iso- account a larger number of galaxies is necessary to ascertain if the lated galaxies with inclination less than 60◦. Galaxies with Teff gradient is an universal property of spiral galaxies. insufficient number of spaxels with measured oxygen abundance The knowledge of the relation between different physical pa- were excluded from our sample. We also rejected from our sample rameters is essential to comprehend which mechanisms drive the galaxies with oxygen abundance measurements for less than 50 formation and evolution of galaxies. For example, in the seminal spaxels and/or galaxies for which spaxels with oxygen abundance paper, Lequeux et al. (1979) calculated the metallicity (traced by measurements cover a range of galactocentric distances lower the ratio between oxygen and hydrogen abundances) and the total than ∼ 1/3 of the its optical radius. Stellar masses, derived from galaxy mass (MT; obtained from atomic hydrogen velocity maps) UV-to-NIR photometry, has been taken from Walcher et al. (2014). for eight irregular and blue compact galaxies and found a clear Our final sample contains 59 galaxies and 49067 spaxels. relation between these parameters (see also Kinman & Davidson The final spatial resolution of the CALIFA data is set by the 1981; Peimbert & Serrano 1982; Rubin et al. 1984; Skillman 1992; fiber size of the PMAS/PPAK integral field spectrophotometer and Tremonti et al. 2004; Pilyugin et al. 2004; S´anchez et al. 2013, it is of the order of 3 arcsec (Husemann et al. 2013). Thus, the spec- among others). Thereafter, Ellison et al. (2008) showed that the trum of each spaxel corresponds to the spectrum emitted by a re- MT-Z relation is affected by a dependence of the metallicity on the gion with a diameter varying from ∼300 pc to ∼1.5 kpc depending star formation rate (SFR), establishing the MT-Z-SFR relation. This on the distance to the galaxy2. Therefore, each observed spectrum result was confirmed by Lara-L´opez et al. (2010); Mannucci et al. comprises the flux of a complex of H II regions and the physical (2010). However, S´anchez et al. (2017) in their recent study, based properties derived represent an averaged value (see discussion in on integral field spectroscopy data, did not find any significant de- Paper I). pendence of the MT-Z relation with the SFR, but they did not ex- clude the existence of such relation on local scales, e.g. in the cen- tral regions of the galaxies. Therefore, the existence of the univer- 2.2 The emission line fluxes sality of T ff gradients, together with the M -Z-SFR relation, would e T The spectrum of each spaxel from the CALIFA DR3 datacubes is produce additional and fundamental concepts of the physical pro- processed in the same way as described in Zinchenko et al. (2016). cesses taking place in galaxies, important insights into how the for- Briefly, the stellar background in all spaxels is fitted using the pub- mation and evolution of massive stars occur as well as their inter- lic version of the STARLIGHT code (Cid Fernandes et al. 2005; action with the ISM. Mateus et al. 2006; Asari et al. 2007) adapted for execution in the In this work, we use the methodology presented by (Dors et al. 2017, hereafter Paper I), to estimate Teff of extragalactic H II re- gions located in a large sample of spiral galaxies. We have taken 1 www.http://califa.caha.es/ 2 −1 −1 advantage of the existence of an homogeneous sample of spectro- We assumed a spatially flat cosmology with H0 = 71 kms Mpc , Ωm = scopic data of H II regions obtained as part of the Calar Alto Legacy 0.270, and Ωvac = 0.730 (Wright 2006)
MNRAS 000, 000–000 (0000) Nebular parameter relations based on CALIFA galaxies 3
− dening law from Izotov et al. (1994). When the measured value of − Hα/Hβ is lower than 2.86 the reddening is adopted to be zero. − Following Paper I, we apply the log([O III]λ5007/Hβ) – II − log([S ]λλ6717,6731/Hα) criterion proposed by Kewley et al. (2001) to separate objects for which the main ionization source are − ⊙ massive stars from those whose main ionization source are shocks − Z/Z of gas and/or active galactic nuclei (AGNs). We consider only the − objects located below the separation line defined by Kewley et al. − (2001), i.e. 39431 spaxels. In Fig. 1, the BPT diagnostic diagram − (Baldwin et al. 1981) for the all the spaxels in our sample is pre- − sented. − − − − logU 2.3 Nebular parameter determinations In order to estimate Teff , we adopted the same method pro- posed in Paper I, where a new calibration between Teff and the R = log([O II]λλ3727,3729/[O III]λ5007) ratio was proposed.
⊙ The method consists of three steps: a) to estimate the metal- licity Z of star forming regions, b) to estimate U from Z and Z/Z [S II]λλ6717,6731 and Hα emission lines, and c) to estimate Teff from U, Z, and R. Entering into details, the first step consists in calculating the metallicity of the gas traced by the oxygen abundance in relation − − − − to the hydrogen one, in units of 12+log(O/H). It is carried out logU using the R3D empirical calibration given by Pilyugin & Grebel (2016). These authors derived oxygen abundances based on di- rect estimations of the electron temperatures for a large sample Figure 2. Coefficients a (upper panel) and b (lower panel) of the linear of H II regions and they obtained relations between these abun- regression T ff = a × R + b as a function of log U and Z/Z⊙. e dances and the emission line flux ratios of oxygen and nitro- gen in relation to Hβ defined as: R2 =[O II]λλ3727,3729/Hβ, = = NorduGrid ARC3 environment of the Ukrainian National Grid. To R3 [O III]λλ4959,5007/Hβ, N2 [N II]λλ6548,6584/Hβ. To use fit the stellar spectra we used 45 synthetic simple stellar popu- this calibration it is necessary to define which branch of the curve lation (SSP) spectra from the evolutionary synthesis models by must be considered, due to the degeneracy in the calibrations. For ≥ − Bruzual & Charlot (2003) with ages from 1 Myr up to 13 Gyr and H II regions with log N2 0.6 the upper branch is assumed and metallicities Z = 0.004, 0.02, and 0.05. We adopted the reddening the relation is the following: law of Cardelli et al. (1989) with RV = 3.1. The resulting stellar = + + (O/H)R,U 8.589 0.022 log(R3/R2) 0.399 log N2 radiation contribution is subtracted from the observed spectrum in + (−0.137 + 0.164 log(R3/R2) + 0.589 log N2) (1) order to measure and analyse the line emission from the gaseous × log R2, component. The line intensities were measured using single Gaus- sian line profile fittings on the pure emission spectra. where (O/H)R,U means 12 +log(O/H)R,U . − The total [O III]λλ4959,5007 flux has been estimated as For H II regions with log N2 < 0.6, the lower branch is as- 1.33× [O III]λ5007 instead of as the sum of the fluxes of both lines. sumed and the relation is = + + These lines originate from transitions from the same energy level, (O/H)R,L 7.932 0.944 log(R3/R2) 0.695 log N2 so their flux ratio can be determined by the transition probability ra- + (0.970 − 0.291 log(R3/R2) − 0.019 log N2) (2) tio, which is very close to 3 (Storey & Zeippen 2000). The strongest × log R2, line, [O III]λ5007, can be measured with higher precision than the where (O/H) means 12 +log(O/H) . To convert the oxygen weakest one. This is particularly important for high-metallicity H II R,L R,L abundance to metallicity Z/Z⊙, we assumed the solar oxygen abun- regions, which have weak high-excitation emission-lines. Simi- dance 12+log(O/H)⊙ = 8.69 (Allende Prieto et al. 2001). larly, the [N II]λλ6548,6584 lines also originate from transitions Following Paper I, the logarithm of the ionization parameter, from the same energy level and the transition probability ratio for log U, is calculated as: those lines is again close to 3 (Storey & Zeippen 2000). Therefore, we estimated its total flux as 1.33 [N II]λ6584. For each spectrum, log U = c × S 2 + d, (3) we measure the fluxes of the [O II]λλ3727,3729, Hβ, [O III]λ5007, S II λλ α c = − . × Z/Z⊙ − . Hα, [N II]λ6584, and [S II]λ6717, 6731. For the further analysis where 2 = log([S ] 6717,6731/H ), 0 26 ( ) 1 54, = − × 2 + × − we selected only those spectra, for which the signal-to-noise ratio and d 3.69 (Z/Z⊙) 5.11 (Z/Z⊙) 5.26. is larger than 5 for each emission-line listed above. The measured The Teff-R relation derived in Paper I is line fluxes are corrected for interstellar reddening using the theoret- Teff = a × R + b, (4) ical Hα to Hβ ratio assuming the standard value of Hα/Hβ = 2.86, and the analytical approximation of the Whitford interstellar red- where R = log([O II]λλ3727,3729/[O III]λ5007). Fig. 2 shows the values of the a and b coefficients as a function of Z/Z⊙ and log U. In this figure, the values of the coefficients are calculated interpolating 3 http://www.nordugrid.org/ the relations given in Table 2 of Paper I.
MNRAS 000, 000–000 (0000) 4 I. A. Zinchenko et al.
4 4