On the Thermally Pulsing Asymptotic Giant Branch Contribution to the Light of Nearby Disk

Eric E. Mart´ınez-Garc´ıa1, Gustavo Bruzual2, Rosa A. Gonz´alez-L´opezlira2, and Lino H. Rodr´ıguez-Merino3 1 CONACYT-Instituto Nacional de Astrof´ısica, Optica´ y Electr´onica, Luis E. Erro 1, Tonantzintla, Puebla, C.P. 72840, M´exico 2 Instituto de Radioastronom´ıay Astrof´ısica, UNAM, Campus Morelia, Michoac´an,M´exico, C.P. 58089 3 Instituto Nacional de Astrof´ısica, Optica´ y Electr´onica, Luis E. Erro 1, Tonantzintla, Puebla, C.P. 72840, M´exico [email protected]

ABSTRACT The study of the luminosity contribution from thermally pulsing asymptotic giant branch (TP- AGB) stars to the stellar populations of galaxies is crucial to determine their physical parameters (e.g., stellar mass and age). We use a sample of 84 nearby disk galaxies to explore diverse stellar population synthesis models with different luminosity contributions from TP-AGB stars. We fit the models to optical and near-infrared (NIR) photometry, on a pixel-bypixel basis. The statistics of the fits show a preference for a low-luminosity contribution (i.e., high mass-to-light ratio in the NIR) from TP-AGB stars. Nevertheless, for 30%–40% of the pixels in our sample a high-luminosity contribution (hence low mass-to-light ratio in the NIR) from TP-AGB stars is favored. According to our findings, the mean TP-AGB star luminosity contribution in nearby disk galaxies may vary with Hubble type. This may be a consequence of the variation of the TP- AGB mass-loss rate with metallicity, if metal-poor stars begin losing mass earlier than metal-rich stars, because of a pre-dust wind that precedes the dust-driven wind. Subject headings: galaxies: evolution – galaxies: photometry – galaxies: stellar content – stars: AGB and post-AGB.

1. Introduction Delgado et al. 2014, 2015, 2016; Cano-D´ıazet al. 2016, 2019; Ibarra-Medel et al. 2016; de Amorim The panchromatic analysis of galaxies has be- et al. 2017; S´anchez et al. 2019; Erroz-Ferrer et al. come an important tool to understand for- 2019) basis, provide very important information to mation and evolution (e.g., Conroy 2013; Lara- complement single spectra analysis of large sam- L´opez et al. 2013; Buat et al. 2014; Vaccari 2016;

arXiv:2011.13989v2 [astro-ph.GA] 17 Feb 2021 ples of objects (see also, Mao et al. 2012; Cava et Poudel et al. 2016; Nagamine et al. 2016; Lim al. 2015; San Roman et al. 2018, 2019). et al. 2017; Ma 2018; Robertson et al. 2019; Ne- However, some issues arise in the resolved grello et al. 2020). For this purpose, the studies studies of galaxies that must not be overlooked. of spatially resolved galaxies, i.e., on a pixel-by- Among common problems are the following: pixel (e.g. Zibetti, Charlot, & Rix 2009; Mentuch Cooper et al. 2012; D´ıaz-Garc´ıaet al. 2015; Sorba (i) The spatial structure of the recovered stel- & Sawicki 2015; Abdurro’uf & Akiyama 2017; lar mass is often different from the expected Mart´ınez-Garc´ıa et al. 2017, 2018), or spaxel- one. by-spaxel (with integral field unit spectrographs, IFUs, e.g. Rosales-Ortega et al. 2010; Gonz´alez (ii) The choice of stellar population synthesis models affects the recovered physical char-

1 acteristics. is similar for m2005 and bc03. In the NIR range (see the Two Micron All Sky Survey, 2MASS,3 H- The first issue was investigated by Mart´ınez- band in the plot), after ≈ 2×108 yr the luminosity Garc´ıaet al. (2017). According to their findings, is higher in m2005 than in bc03. This is the stellar the recovered stellar mass map of a nearby galaxy age where the thermally pulsing asymptotic giant (M51), obtained by fitting stellar population syn- branch (TP-AGB) stars dominate the stellar pop- thesis models to observations (pixel-by-pixel in ulation luminosity. Also, according to the plot for this case), results in an artificial spatial struc- λ & 2.5µm, the lu