text:
# The Nature of Starburst Activity in M 82 1footnote 11footnote 1Based on observations with ISO, an ESA project with instruments funded by ESA Member States (especially the PI countries: France, Germany, the Netherlands, and the United Kingdom) and with the participation of ISAS and NASA. The SWS is a joint project of SRON and MPE.

## 1 INTRODUCTION

In recent years, it has become clear that starburst galaxies are important constituents of the Universe at all accessible redshifts (e.g., Heckman, 1998; Adelberger & Steidel, 2000; Giavalisco, 2002; Shapley et al., 2001, 2003; Chapman et al., 2003). However, despite extensive studies in the past two decades, a detailed and quantitative understanding of the starburst phenomenon is still lacking. Crucial issues that remain open include the evolution and feedback effects of starburst activity, its triggering and quenching mechanisms, and the mass distribution of the stars formed in starbursts. Progress has been hindered by the scarcity of spatially resolved data. Furthermore, high resolution optical and ultraviolet studies are often hampered by severe dust obscuration.

In this context, we have undertaken a study of the archetypal starburst galaxy M 82 based on near-infrared (near-IR) integral field spectroscopy and mid-infrared (mid-IR) spectroscopy. Our observations provide detailed information on key features tracing various components of the interstellar medium and of the stellar population, deep into the obscured star-forming regions. The data, along with results of nebular analysis and population synthesis, are described by Förster Schreiber et al. (2001, hereafter paper I), laying the observational basis essential for the starburst modeling presented in this paper. In a closely related study, Thornley et al. (2000) modeled the mid-IR line emission of 27 starburst galaxies to examine the issues of massive star formation and evolution in starburst environments. Their results are particularly relevant to our modeling of M 82.

By virtue of its proximity and brightness, M 82 is an ideal target for investigations of the starburst phenomenon. Following the seminal paper by Rieke et al. (1980), several authors have applied evolutionary synthesis models to M 82 (Bernlöhr, 1992; Rieke et al., 1993; Doane & Mathews, 1993; Satyapal et al., 1997; Colbert et al., 1999, among others). However, a number of issues remain controversial. For instance, it is still debated whether the initial mass function (IMF) is biased against low-mass stars (see in particular Rieke et al. 1993 and Satyapal et al. 1997). It has also been suggested that the formation of stars with masses $`30\mathrm{M}_{}`$ may be suppressed (e.g., Puxley et al., 1989). Furthermore, in analyzing the relative distributions of various gaseous and stellar components, some authors have proposed inside-out propagation of starburst activity in M 82 (e.g., McLeod et al., 1993; Satyapal et al., 1997) while others have argued in favour of the opposite scenario (e.g., Shen & Lo, 1995; de Grijs et al., 2000, 2001).

The previous models of M 82 were either optimized to best reproduce the global properties or focussed on selected stellar clusters. Given the complexity of M 82, a full understanding requires spatially detailed information of various diagnostics tracing different star formation epochs. Our new IR data, complemented with results at other wavelengths from the literature, allow us to probe the star formation history in the central regions of M 82 “continuously” in space on scales as small as 25 pc and in time up to at least $`50\mathrm{Myr}`$ ago. With the application of starburst models we have developed, we use these data to re-examine the issue of the IMF and constrain the spatial and temporal evolution of the star formation activity. Our model results provide quantitative constraints for the triggering mechanisms and feedback effects of starburst activity in M 82.

The paper is organized as follows. Section 2 summarizes the observational constraints. Section 3 describes the starburst models which are applied to selected regions in § 4. Section 5 presents the spatially detailed modeling. Section 6 discusses the results and § 7 summarizes the paper.

## 2 OBSERVATIONAL CONSTRAINTS

Our modeling of M 82 uses the observations presented and analyzed in paper I. The data consist of near-IR $`H`$\- and $`K`$-band integral field spectroscopy at $`R\lambda /\mathrm{\Delta }\lambda 1000`$ obtained with the Max-Planck-Institut für extraterrestrische Physik (MPE) 3D instrument (Weitzel et al., 1996) and of $`\lambda =2.445\mu \mathrm{m}`$ spectroscopy at $`R5002000`$ from the Short Wavelength Spectrometer (SWS; de Graauw et al. 1996) on board the Infrared Space Observatory (ISO; Kessler et al. 1996). The 3D field of view is nearly parallel to the galactic plane, includes the nucleus, and extends to the west out to about 200 pc. The SWS apertures cover regions including and centered on the 3D field of view. The 3D data have a spatial resolution of $`1.5^{\prime \prime }`$, corresponding to about 25 pc at the distance of M 82 (3.3 Mpc; Freedman & Madore 1988). We complemented these data with results from the literature, providing further essential constraints.

We selected five representative regions for a detailed analysis. These include the “starburst core” of M 82 corresponding to a 500 pc–diameter region centered on the nucleus, the 3D field of view covering the most intense starburst regions, the central 35 pc at the nucleus, and two $`35\times 35\mathrm{pc}`$ regions positioned at the brightest Br$`\gamma `$ sources observed with 3D, which we designated M82:Br1 and M82:Br2 in paper I ($`10^{\prime \prime }`$ and $`5^{\prime \prime }`$ southwest of the nucleus, respectively). For simplicity, we refer to the latter three regions as the “nucleus,” “B1,” and “B2.”

For each region, Table 1 lists the main constraints derived in paper I. These include the intrinsic $`K`$-band and Lyman continuum luminosities ($`L_K`$ and $`L_{\mathrm{Lyc}}`$) 222We use the following definitions: $`L_K`$ is the luminosity in the $`K`$ bandpass $`\lambda =1.92.5\mu \mathrm{m}`$ (Wamsteker, 1981) expressed in units of the total solar luminosity $`\mathrm{L}_{}=3.85\times 10^{26}\mathrm{W}`$ while $`L_{\mathrm{Lyc}}`$ is the Lyman continuum photon emission rate times an average ionizing photon energy of 15 eV., the total bolometric luminosities from OB and cool evolved stars as well as the separate contribution of OB stars ($`L_{\mathrm{bol}}`$ and $`L_{\mathrm{bol}}^{\mathrm{OB}}`$), the stellar masses ($`M^{}`$), the supernova rates ($`\nu _{\mathrm{SN}}`$), the CO 1.62 and 2.29 µm bandhead equivalent widths ($`W_{1.62}`$ and $`W_{2.29}`$), and the neon fine-structure line ratio \[Ne III\] 15.6 $`\mu \mathrm{m}`$ / \[Ne II\] 12.8$`\mu \mathrm{m}`$. The values of the $`L_K`$/$`L_{\mathrm{Lyc}}`$, $`L_{\mathrm{bol}}`$/$`L_{\mathrm{Lyc}}`$, $`L_{\mathrm{bol}}^{\mathrm{OB}}`$/$`L_{\mathrm{Lyc}}`$, $`M^{}`$/$`L_K`$, and $`\nu _{\mathrm{SN}}`$/$`L_{\mathrm{bol}}`$ are also given. Except for the starburst core, the neon ratio is an “equivalent” ratio as discussed below.

The $`K`$-band continuum luminosities $`L_K`$ and the equivalent widths $`W_{2.29}`$ are corrected for the contribution and dilution by hot dust emission (which is negligible for $`W_{1.62}`$). We assumed that 70% of the OB stars’ bolometric luminosity is detected at IR wavelengths ($`\lambda =5300\mu \mathrm{m}`$) as reprocessed emission by dust, and that 30% escapes perpendicular to the galactic disk. The stellar masses $`M^{}`$ for the nucleus and the starburst core were computed by subtracting the mass of the gas component (mainly in the form of $`\mathrm{H}_2`$) from the total dynamical mass. Estimates of supernova explosion rates $`\nu _{\mathrm{SN}}`$ suffer from rather large uncertainties (see paper I), so related constraints will be used as consistency arguments. Among the nebular emission line ratios sensitive to the ionizing OB stars available from our data sets, we chose \[Ne III\] 15.6 $`\mu \mathrm{m}`$ / \[Ne II\] 12.8$`\mu \mathrm{m}`$ for this work. For the starburst core, the ratio is that measured with the SWS, corrected for extinction. For the other regions, we applied the single-star photoionization models of paper I to derive “equivalent” neon ratios, which are those predicted for the OB stars effective temperatures inferred from the observed near-IR H to He recombination line ratios.

We also modeled individual regions corresponding to rebinned $`1^{\prime \prime }\times 1^{\prime \prime }`$ pixels covering the 3D field of view in its entirety. The set of constraints for each pixel, from the 3D data and the 12.4 µm map of Telesco & Gezari (1992), includes the intrinsic $`L_K`$, $`L_{\mathrm{Lyc}}`$, $`L_{\mathrm{bol}}^{\mathrm{OB}}`$, $`W_{2.29}`$, and the equivalent neon ratio. These constraints are the most useful and relevant ones for the spatially-detailed modeling based on the results of the selected regions.

## 3 STARBURST MODELS

Our starburst models combine evolutionary synthesis of stellar clusters with photoionization modeling of the surrounding gas. They are optimized for applications to $`\lambda =145\mu \mathrm{m}`$ observations of starburst galaxies. They assume that a given set of constraints applies to a “homogeneous” stellar population, i.e. either a single cluster or an ensemble of clusters with identical IMF shape and cutoffs (see Thornley et al., 2000).

### 3.1 Prediction of Stellar Properties

We computed the integrated properties of stellar clusters using the evolutionary synthesis code STARS (Sternberg, 1998; Thornley et al., 2000). STARS is similar to other codes which have been presented in the literature (e.g., Bruzual & Charlot, 1993; Cerviño & Mas-Hesse, 1994; Fioc & Rocca-Volmerange, 1997; Leitherer et al., 1999). STARS performs conventional synthesis, assigning fixed evolutionary tracks to an appropriate set of mass bins. STARS employs the Geneva stellar tracks (Schaller et al., 1992); for our analysis, we selected the solar-metallicity tracks with normal mass-loss rates. Although the tracks for enhanced mass-loss rates may provide better fits to observations (see Meynet et al., 1994, and references therein), the choice has little impact for this work. The differences in the synthetized properties of interest here are much smaller than those produced, e.g., by varying the IMF parameters and would imply changes $`1\mathrm{Myr}`$ in the derived ages, comparable to those from the measurement uncertainties. We adopted a time-independent power-law IMF $`\mathrm{d}N/\mathrm{d}mm^\alpha `$ between lower and upper mass limits $`m_{\mathrm{low}}`$ and $`m_{\mathrm{up}}`$. The stars are assumed to form at a rate $`R(t_\mathrm{b})=R_0e^{t_\mathrm{b}/t_{\mathrm{sc}}}`$, where $`R_0`$ is the initial star formation rate (expressed in $`\mathrm{M}_{}\mathrm{yr}^1`$), $`t_\mathrm{b}`$ is the burst age, and $`t_{\mathrm{sc}}`$ is the burst decay timescale.

The integrated spectra of starburst galaxies in the $`145\mu \mathrm{m}`$ range are dominated at the short-wavelength end by the direct light of cool evolved stars and at the long-wavelength end by radiation of mostly OB stars that has been reprocessed by interstellar dust and gas. These populations are treated with particular care in STARS. Specifically, the non local thermodynamical equilibrium (non-LTE) model atmospheres of Pauldrach et al. (1998), which supersede the earlier non-LTE atmospheres of Sellmaier et al. (1996), are chosen to represent adequately the spectral energy distribution (SED) of the hottest stars (see also Pauldrach et al., 2001, and references therein). Empirical photometric data are used to better account for the properties of the coolest stars. Finally, the Geneva tracks for intermediate-mass stars have been extended to include the thermally-pulsing asymptotic giant branch phase (TP-AGB), which can contribute significantly to the integrated $`K`$-band luminosity. STARS also predicts the strength of various near-IR stellar absorption features. Further details on the above aspects are given in appendix A.

The Pauldrach et al. (1998) models include the non-LTE radiative transfer and hydrodynamical treatment of radiatively-driven, spherically expanding steady-state winds and account for the combined effects of line blocking and blanketing on the radiative transfer and energy balance. They provide excellent matches to observed high-resolution far-ultraviolet spectra of hot stars (Pauldrach et al., 2001) and nebular photoionization computations which incorporate these model atmospheres successfully reproduce the relative intensities of infrared emission lines measured in Galactic H II regions (Giveon et al., 2002). As discussed by Giveon et al. (2002) and Thornley et al. (2000), the inclusion of the Pauldrach et al. models is critical for the photoionization modeling because their SEDs are significantly harder than for hydrostatic plane-parallel LTE models such as those of Kurucz (1992). Accounting for the TP-AGB phase is most important for the CO bandheads because it results notably in larger equivalent widths (EWs) at ages $`50\mathrm{Myr}`$ (see, e.g., Origlia & Oliva, 2000). However, as will be seen below, our results do not depend much on this feature.

### 3.2 Prediction of Nebular Properties

We modeled the neon ratio produced in nebulae excited by the integrated stellar SED computed by STARS using the photoionization code CLOUDY version C90.05 (Ferland, 1996). The nebulae were represented as a single gas shell surrounding a central point-like ionizing source. In such “central cluster” models, the nebular conditions are specified by the gas and dust composition, the hydrogen gas density $`n_\mathrm{H}`$, and the ionization parameter $`UQ_{\mathrm{Lyc}}/\left(4\pi r^2n_\mathrm{H}c\right)`$, where $`Q_{\mathrm{Lyc}}`$ is the intrinsic Lyman continuum photon emission rate of the source, $`r`$ is the shell radius, and $`c`$ is the speed of light.

We adopted the parameter values derived from the nebular analysis of paper I: $`n_\mathrm{H}=300\mathrm{cm}^3`$, an effective ionization parameter with $`\mathrm{log}U_{\mathrm{eff}}=2.3\mathrm{dex}`$, and solar photospheric gas-phase abundances (from Grevesse & Anders, 1989; Grevesse & Noels, 1993, as implemented in CLOUDY v. C90.05). As argued in paper I, a more realistic representation of the H II regions in M 82 consists of a well-mixed distribution of ionizing clusters and gas clouds. The average ionization parameter at the outer boundaries of the gas clouds derived from the observed properties assuming such a “random distribution” model constitutes the appropriate effective value for modeling in the framework of the idealized central cluster geometry. Interestingly, we found nearly identical $`\mathrm{log}U_{\mathrm{eff}}`$ throughout the starburst core of M 82 on scales of a few tens of parsecs to $`500\mathrm{pc}`$. As discussed in paper I, this appears to indicate a uniform star formation efficiency across the observed regions.

Following Thornley et al. (2000), we computed the neon ratio for bursts with decay timescale $`t_{\mathrm{sc}}=1\mathrm{Myr}`$ by assuming that $`U_{\mathrm{eff}}`$ varies proportionally with $`Q_{\mathrm{Lyc}}`$ and that the maximum value of $`\mathrm{log}U_{\mathrm{eff}}`$ equals $`2.3\mathrm{dex}`$. For bursts with longer decay timescales, we convolved the corresponding star formation rates $`R(t_\mathrm{b})`$ with the time-dependent neon ratio for the short decay bursts. We neglected the effects of dust grains mixed with the ionized gas within the nebulae; this has little impact on our model results (see § 4.5).

### 3.3 Modeling Procedure

We approached the modeling of M 82 by first making physical and qualitative arguments to constrain the possible ranges for the parameters, and then performed quantitative fits. This led to physically plausible and successful models although they may not be unique. We first investigated, from the properties of selected regions, the general characteristics of starburst activity: the IMF cutoffs and the typical burst decay timescale (§ 4). Having determined these parameters and identified the most useful age indicators, we then modeled individual regions across the entire 3D field of view to constrain the spatial and temporal evolution of starburst activity (§ 5).

Most of our constraints are sensitive to $`m_{\mathrm{up}}`$, $`t_{\mathrm{sc}}`$, and $`t_\mathrm{b}`$, so we investigated these parameters simultaneously. The total stellar mass depends only on the low-mass IMF (shape and $`m_{\mathrm{low}}`$), which we examined independently. We used ratios of different properties derived from data obtained at similar wavelengths along with the CO bandheads EWs, minimizing uncertainties from the extinction corrections and from the distance assumed for M 82. The exception is $`\nu _{\mathrm{SN}}/L_{\mathrm{bol}}`$ but we used it mainly as consistency check because of the uncertainties on $`\nu _{\mathrm{SN}}`$ (§ 2). We determined the initial star formation rates $`R_0`$ by comparing the predicted luminosities to the measured values. We kept the other parameters constant; the effects of variations in the most critical ones will be discussed when appropriate. Table 2 summarizes the ranges considered or the fixed values adopted for the model parameters.

## 4 MODELING OF SELECTED REGIONS IN M 82

Figure 1 compares the properties of the selected regions from Table 1 to model predictions. The models were computed for $`m_{\mathrm{up}}=100\mathrm{M}_{}`$ and are shown for burst decay timescales of 1 Myr, 5 Myr, 20 Myr, and 1 Gyr. The diagnostics considered are not affected by the low-mass star population so a fiducial $`m_{\mathrm{low}}=1\mathrm{M}_{}`$ was here adopted. This comparison reveals age differences between the regions, in particular between the nucleus (older) and the regions B1 and B2 (younger). The exact ages depend on the upper mass cutoff and burst decay timescale, which we constrained first.

### 4.1 Upper Mass Cutoff and Burst Timescale

The neon and $`L_{\mathrm{bol}}`$/$`L_{\mathrm{Lyc}}`$ ratios are our most sensitive diagnostics to the upper mass cutoff. Figure 2 shows model computations for these ratios with $`m_{\mathrm{up}}=25`$, 30, 35, 50, and 100$`\mathrm{M}_{}`$, and for $`t_{\mathrm{sc}}=1`$ and 5 Myr. The upper mass cutoff is best constrained from the properties of the young massive stars only. Therefore, the observed bolometric luminosities for selected regions in Figure 2 exclude the contribution from cool evolved stars so that the $`L_{\mathrm{bol}}^{\mathrm{OB}}`$/$`L_{\mathrm{Lyc}}`$ ratios are compared to the models.

Both diagnostics imply similar dominant OB star populations for the various regions considered. The small values of the neon ratio suggest a relative paucity of ionizing photons with energies $`>3\mathrm{Ryd}`$ (the ionization potential of $`\mathrm{Ne}^+`$) while the low $`L_{\mathrm{bol}}^{\mathrm{OB}}`$/$`L_{\mathrm{Lyc}}`$ ratios imply that lower energy Lyman continuum photons are being produced efficiently. From the comparison with the models, two alternative interpretations are possible: (1) a high $`m_{\mathrm{up}}50\mathrm{M}_{}`$, a short timescale of at most a few million years, and the softening of the ionizing radiation field attributable to rapid aging of the starburst, or (2) a lower $`m_{\mathrm{up}}`$ down to $`30\mathrm{M}_{}`$, with longer timescales possible. Upper mass cutoffs below 30$`\mathrm{M}_{}`$ are ruled out, but the data do not allow us to constrain unequivocally $`m_{\mathrm{up}}`$ and $`t_{\mathrm{sc}}`$ for $`m_{\mathrm{up}}30\mathrm{M}_{}`$.

As summarized by Thornley et al. (2000), a variety of results on local templates of starburst regions as well as other large-scale starburst systems make high upper mass cutoffs more plausible. In M 82, the low $`L_{\mathrm{bol}}^{\mathrm{OB}}`$/$`L_{\mathrm{Lyc}}`$ values for B1 and B2 can only be reproduced with $`m_{\mathrm{up}}50\mathrm{M}_{}`$. The neon and $`L_{\mathrm{bol}}^{\mathrm{OB}}`$/$`L_{\mathrm{Lyc}}`$ ratios measured for all selected regions lie within the ranges determined for the sample of dusty solar-metallicity IR-luminous starburst galaxies of Thornley et al. (2000). These ranges, indicated in Figure 2, are compatible with $`m_{\mathrm{up}}50\mathrm{M}_{}`$ and differences among the sample sources attributable to a range in burst ages.

In addition, the poor correlation in M 82 on scales of a few tens of parsecs between the spatial distributions of the near-IR continuum emission and CO bandhead EWs, and that of the ionized gas (e.g., Br$`\gamma `$ and He I 2.06 µm) and IR dust emission provides strong evidence for short decay timescales in individual burst sites (see the 3D images of paper I; also, e.g., Telesco et al. 1991; Larkin et al. 1994; Satyapal et al. 1997). Since red supergiants dominate the near-IR continuum emission in M 82 (§ 4.3), the poor correlation implies that star formation activity was not maintained locally at high levels longer than the typical main-sequence lifetime of the most massive red supergiants progenitors ($`5\mathrm{Myr}`$). This situation is reminiscent of the nucleus of M 83, where Puxley et al. (1997) found that the K-band continuum emission and CO bandhead strength are spatially anticorrelated with the Br$`\gamma `$ EW, which they interpreted, together with starburst models, as indicating an evolutionary sequence and burst timescales of $`15\mathrm{Myr}`$. Massive star formation is expected to inhibit further star formation locally already after a few million years owing to strong stellar winds and supernova explosions disrupting the surrounding interstellar medium (ISM). The $`100\mathrm{pc}`$–size gas hole around the brightest radio supernova remnant $`41.9+58`$ seems to provide a particularly striking example of such effects in M 82 (e.g., Weiß et al., 1999).

The most remarkable outcome of our models is that for $`m_{\mathrm{up}}50\mathrm{M}_{}`$, the OB star-dominated properties of the 3D field of view and the starburst core imply short burst decay timescales on large spatial scales ($`100500\mathrm{pc}`$). This is not due to a few regions dominating the integrated properties; for instance, B1 and B2 contribute together less than 15% to the $`L_{\mathrm{Lyc}}`$ and $`L_{\mathrm{bol}}^{\mathrm{OB}}`$ measured in the 3D field of view. Our results suggest that in M 82 starburst activity can have strong negative feedback effects globally as well. We discuss this point further in § 6.

### 4.2 Successive Starburst Events

In the rest of this paper, we will adopt $`m_{\mathrm{up}}=100\mathrm{M}_{}`$ and a representative $`t_{\mathrm{sc}}=1\mathrm{Myr}`$ as justified above. Examination of Figures 1 and 2 then reveals that for each region, there is a systematic increase in burst age from that implied by the tracers of OB stars (neon and $`L_{\mathrm{bol}}^{\mathrm{OB}}`$/$`L_{\mathrm{Lyc}}`$ ratios) to that inferred from the tracers of cool evolved stars (CO bandheads). The differences are about 5 Myr. The $`L_K`$/$`L_{\mathrm{Lyc}}`$ ratios, measuring the relative populations of these stellar components, correspond to intermediate ages.

Figure 3 illustrates the difficulty of single burst models to meet simultaneously the constraints for each of the selected regions. The figure compares the regions’ properties to evolutionary tracks for single bursts in diagrams of the neon, $`L_{\mathrm{bol}}`$/$`L_{\mathrm{Lyc}}`$, and $`L_K`$/$`L_{\mathrm{Lyc}}`$ ratios versus $`W_{2.29}`$. At ages corresponding to the measured values of $`W_{2.29}`$, models for $`t_{\mathrm{sc}}=1\mathrm{Myr}`$ and $`m_{\mathrm{up}}=100\mathrm{M}_{}`$ predict neon ratios 10 (B1) to 34 (nucleus) times lower than measured, $`L_{\mathrm{bol}}`$/$`L_{\mathrm{Lyc}}`$ ratios 6 (nucleus) to 23 (B2) times higher, and $`L_K`$/$`L_{\mathrm{Lyc}}`$ ratios 4 (nucleus) to 33 (B2) times higher. The constraints are still not satisfied simultaneously when augmenting the burst decay timescale to 5 Myr with, in particular, predicted neon ratios about an order of magnitude higher than measured at ages corresponding to the observed values of $`W_{2.29}`$. Allowing for lower $`m_{\mathrm{up}}`$ down to 50$`\mathrm{M}_{}`$ only worsens the discrepancies.

The age sequence between the various diagnostics may suggest that the hot massive stars and the cool evolved stars belong to populations formed in distinct, successive starburst events. In this scenario, the “young bursts” account for most of the ionizing radiation while the “old bursts” dominate the near-IR properties. Rieke et al. (1993) first proposed such two-burst models of the entire starburst core of M 82 for similar reasons. We thus considered models consisting of two sequential bursts, each with $`t_{\mathrm{sc}}=1\mathrm{Myr}`$ and $`m_{\mathrm{up}}=100\mathrm{M}_{}`$, and separated in time by more than one decay timescale. Figure 4 shows similar evolutionary tracks as Figure 3 for the best-fit time separation and relative burst intensities for each region derived in § 4.3. From Figures 3 and 4, it is also difficult to disentangle between single and two-burst models. However, since the neon and $`L_{\mathrm{bol}}^{\mathrm{OB}}`$/$`L_{\mathrm{Lyc}}`$ ratios seem to favour the shortest decay timescales, successive short bursts could provide a natural explanation for the observed strength of the CO bandheads. The two-burst scenario also appears plausible in view of the morphology and nearly edge-on orientation of M 82 (§ 6).

We remark that given the uncertainties of the data and models, it may be possible to find single bursts with finely-tuned parameters which reproduce all the observed properties. Exponentially decaying bursts with $`t_{\mathrm{sc}}`$ intermediate between 1 and 5 Myr, for instance, could meet the neon ratio constraint while violating the $`L_{\mathrm{bol}}`$/$`L_{\mathrm{Lyc}}`$ and $`L_K`$/$`L_{\mathrm{Lyc}}`$ constraints at a level that is more severe than for the double bursts but probably still acceptable. As another example, allowing lower $`m_{\mathrm{up}}`$ down to 30$`\mathrm{M}_{}`$ together with longer $`t_{\mathrm{sc}}`$ of $`1020\mathrm{Myr}`$ for single bursts may improve the fit to the neon ratios but the mismatch is only redistributed to the other properties. Similar difficulties in modeling comparable sets of constraints have been encoutered for other dusty IR-selected starbursts, e.g., by Engelbracht et al. (1998) for NGC 253 and Goldader et al. (1997) for a sample of (ultra)luminous IR galaxies. With our data and models of M 82, together with the obvious complexity of the problem, it is not possible to find a unique solution.

### 4.3 Burst Age and Strength

Adopting the two-burst scenario as working hypothesis, we constrained the age of the young bursts from the neon and $`L_{\mathrm{bol}}^{\mathrm{OB}}`$/$`L_{\mathrm{Lyc}}`$ ratios, and of the old bursts from the CO bandheads. The initial star formation rates $`R_0`$ (with $`m_{\mathrm{low}}`$ fixed at 1$`\mathrm{M}_{}`$) were adjusted so that the total contribution from both bursts to $`L_{\mathrm{Lyc}}`$ and $`L_K`$ reproduces the observed values. For identical IMF parameters and burst decay timescales, the values of $`R_0`$ provide a meaningful measure of the relative strengths of the bursts.

For the young bursts, the ages inferred from each indicator agree marginally, the neon ratio implying slightly older ages. This discrepancy would not be solved by other choices of $`m_{\mathrm{up}}`$ and $`t_{\mathrm{sc}}`$. The difference is comparable to the burst timescale and, in view of the measurements and modeling uncertainties, is probably not significant. The $`L_{\mathrm{bol}}^{\mathrm{OB}}`$/$`L_{\mathrm{Lyc}}`$ ratio may trace a slightly different population as it is sensitive to a somewhat lower stellar mass range than the neon ratio is. The effects of a distribution of cluster masses and luminosities or of dust grains within the nebulae, neglected here, would increase the discrepancy (§ 4.5). We therefore assigned equal weights to the diagnostic ratios.

For the old bursts, an accurate age determination is only possible for the nucleus. It has very deep CO bandheads which the models reproduce only at ages of $`815\mathrm{Myr}`$, consistent with the average spectral type K5 I derived from detailed spectral synthesis in paper I. The EWs at the nucleus further support short burst timescales: for $`t_{\mathrm{sc}}5\mathrm{Myr}`$, the contribution from older populations with shallower bandheads and/or photospheric and nebular featureless continuum emission from OB stars results in predicted EWs lower than observed, for any burst age. Figure 5 indicates that older ages when intermediate-mass stars reach the end of the asymptotic giant branch ($`50\mathrm{Myr}1\mathrm{Gyr}`$) are also ruled out from $`W_{2.29}`$.

For the other selected regions, we can set firm lower limits on the ages but the EWs alone do not provide strong constraints on upper limits. For comparison with our data and models, we indicate in Figure 5 the ranges observed by Oliva et al. (1995) for old ($`110\mathrm{Gyr}`$) stellar populations characteristic of the central regions of elliptical galaxies and bulges of spiral galaxies (hereafter “normal populations”). The EWs at B1 are somewhat lower than for normal populations and may indicate an early stage when the first red supergiants appear. For B2, the 3D field of view, and the starburst core, they are slightly higher and imply possible upper limits on the age of $`50\mathrm{Myr}`$ and on the timescale of $`5\mathrm{Myr}`$.

The $`M^{}/L_K`$ ratio provides an additional constraint. Empirical determinations for normal populations lie in the range $`1030\mathrm{M}_{}/\mathrm{L}_{}`$ (Devereux, Becklin, & Scoville, 1987; Oliva et al., 1995; Hunt et al., 1999). 333Our units for the $`M^{}/L_K`$ ratio differ from the conventional definition expressed in $`(M/L_K)_{}`$. The quoted range has been computed in $`\mathrm{M}_{}/\mathrm{L}_{}`$ from the data presented in the given references, assuming a $`K`$ bandpass $`\lambda =1.92.5\mu \mathrm{m}`$ and expressing $`L_K`$ in $`\mathrm{L}_{}=3.85\times 10^{26}\mathrm{W}`$. As argued in paper I, the substantially lower values $`1\mathrm{M}_{}/\mathrm{L}_{}`$ measured for M 82 between radii of $`10`$ and 250 pc indicate very large contributions from red supergiants to the near-IR continuum emission, even in regions of low surface brightness. Consequently, ages $`50\mathrm{Myr}`$ are the most plausible. We adopted the solutions corresponding to the rising part of the predicted EWs with $`t_{\mathrm{sc}}=1\mathrm{Myr}`$ but older ages of $`1530\mathrm{Myr}`$ are also possible. Our choice implies lower $`R_0`$ for the old bursts since the predicted $`L_K`$ peaks around 10 Myr. The younger age solutions are unaffected by uncertainties related to numerical fluctuations due to our conventional synthesis technique and to our TP-AGB implementation, which may become a concern at ages $`50\mathrm{Myr}`$ (see appendix A).

### 4.4 Model Results for Selected Regions

Tables 3 and 4 report our two-burst model results for the selected regions. The derived burst strengths outline the particular intensity of starburst activity in the central few tens of parsecs of M 82 about 10 Myr ago. Our model for the starburst core agrees well with the results of Rieke et al. (1993) who found that separate bursts occurring about 5 and 13 Myr ago, each lasting a few million years, could reproduce the observed properties. The small differences between our model and theirs reflect mainly differences in the adopted set of constraints and model assumptions (such as their choice of Gaussians to represent the star formation rate).

We estimated the formal uncertainties on the relative burst ages and strengths $`R_0`$ by varying in turn the observational constraints within the ranges corresponding to the measurement errors. They are typically $`0.5\mathrm{Myr}`$ and $`25\%40\%`$ for the ages and strengths, respectively. For the old bursts, the minimum age is fairly well constrained by the CO bandhead EWs whereas the maximum age from the quoted uncertainties constitutes a lower limit since we did not account here for the older possible solutions. The absolute burst ages and intensities may suffer from larger uncertainties related, e.g., to the burst decay timescale and shape of the star formation rate. The relative ages and intensities are less affected because the evolution of the properties considered are dominated by stellar evolution during the early phases of a burst and for decay timescales comparable to the main-sequence lifetime of the most massive stars.

Figure 6 illustrates the results for the most relevant properties: $`L_K`$, $`L_{\mathrm{Lyc}}`$, $`L_{\mathrm{bol}}`$, $`W_{2.29}`$, $`W_{1.62}`$, and the neon ratio. For each region, the curves show the variations of the properties as a function of the time elapsed since the onset of the old burst and are normalized to the observed values. They represent the combination of the two bursts, weighted by the respective values of $`R_0`$ for each burst. The models for all regions reproduce the main constraints within a factor of three or better. The masses computed for the nucleus and the starburst core are, however, substantially lower than the measured values. The mass is very sensitive to the lower cutoff and shape of the IMF, which we discuss in § 4.6.

Our results are roughly consistent with the constraints imposed by $`\nu _{\mathrm{SN}}`$ and $`10^{12}\nu _{\mathrm{SN}}/L_{\mathrm{bol}}`$. A good agreement for the regions observed with 3D is not expected because of the large uncertainties in deriving $`\nu _{\mathrm{SN}}`$ from \[Fe II\] 1.644 µm line fluxes, as discussed in paper I. The ages inferred for the old bursts satisfy the “timing constraint” implied by the extent of the X-ray halo tracing the supernova-driven wind of M 82, which requires a minimum starburst age of about 10 Myr (Heckman, Armus, & Miley, 1990; Rieke et al., 1993; Lehnert, Heckman, & Weaver, 1999).

### 4.5 Model Uncertainties

We investigated the effects of various sources of uncertainties on our model results so far, including the assumption of a homogeneous cluster population, the nebular parameters, the shape of the IMF, and the synthetic model atmospheres. Their effects have been described for the neon ratio by Thornley et al. (2000) whom we refer to for further details. Figure 7 shows in addition the effects on the $`L_{\mathrm{bol}}`$/$`L_{\mathrm{Lyc}}`$ and $`L_K`$/$`L_{\mathrm{Lyc}}`$ ratios, and on $`W_{2.29}`$. None of the above sources of uncertainties affect our conclusions concerning the upper mass cutoff of the IMF or the ages and relative intensities derived for the bursts. While possible alternative assumptions or variations of the nebular parameters and IMF shape result in measurable differences in the model predictions, the effects remain small compared to those of aging or changes in $`m_{\mathrm{up}}`$.

In particular, the largest effects are obtained for a “heterogeneous” population of stellar clusters following the luminosity function (LF) derived by Thornley et al. (2000) 444For reference, the LF in terms of cluster H ionizing rate $`Q_{\mathrm{Lyc}}^{\mathrm{cl}}`$ is $`\mathrm{d}N/\mathrm{d}\left(\mathrm{log}Q_{\mathrm{Lyc}}^{\mathrm{cl}}\right)\left(Q_{\mathrm{Lyc}}^{\mathrm{cl}}\right)^\beta `$, with $`\beta =0.19`$ for $`Q_{\mathrm{Lyc}}^{\mathrm{cl}}=10^{45}10^{49.5}\mathrm{s}^1`$ and $`\beta =1`$ for $`Q_{\mathrm{Lyc}}^{\mathrm{cl}}=10^{49.5}10^{53}\mathrm{s}^1`$.. This LF was constructed assuming it reflects a cluster mass function and imposes a distribution in effective $`m_{\mathrm{up}}`$ of clusters (at fixed IMF slope and $`m_{\mathrm{low}}`$). The large proportion of small clusters results in a “down-weighting” effect with, for instance, lower neon ratios and shallower CO bandheads due to the reduction in relative population of massive stars. However, the effects are but a factor of two or less. Similar differences for the neon ratio are obtained by varying the ionization parameter within the possible range $`2.5\mathrm{log}U_{\mathrm{eff}}2`$ for M 82 (paper 1) or adopting the softer Kurucz (1992) SEDs for hot stars instead of the Pauldrach et al. (1998) SEDs (see Thornley et al., 2000). The effects of dust within the H II regions, illustrated with the Orion nebular composition, are significantly smaller. By absorbing ionizing photons, dust would also reduce the $`L_{\mathrm{Lyc}}`$ in the models. This is difficult to quantify, but the resulting increase in synthetized $`L_{\mathrm{bol}}`$/$`L_{\mathrm{Lyc}}`$ and $`L_K`$/$`L_{\mathrm{Lyc}}`$ would only imply younger ages and keep the inferred burst decay timescales short.

Figure 8 shows that the alternative IMFs we considered, from Miller & Scalo (1979) and that derived by Eisenhauer et al. (1998) for the Galactic high-mass star-forming region NGC 3603, bracket the range of slopes at intermediate and high masses determined for young clusters and OB associations in the Milky Way and Large Magellanic Cloud (e.g., Hunter et al. 1997 and references therein; Brandl et al. 1996; Massey & Hunter 1998; Kroupa 2001). Changes in the IMF shape have the smallest effects in our models. As a consequence, our set of observational constraints does not allow us to investigate the shape of the intermediate- and high-mass IMF in M 82.

Potential problems with the stellar evolutionary tracks represent an additional source of uncertainties. A detailed discussion is beyond the scope of this paper but we briefly mention two points relevant to this work. Origlia and co-workers (Origlia et al., 1999; Origlia & Oliva, 2000) warned that model predictions of the CO bandheads can be significantly underestimated at ages $`10\mathrm{Myr}`$ if the TP-AGB phase is not accounted for (see also Charlot & Bruzual, 1991; Lançon & Rocca-Volmerange, 1994) and because of the possible inadequacy of current stellar tracks for the evolution of massive stars in the red supergiant phase, most severely at sub-solar metallicities. Our starburst models account for the TP-AGB (see appendix A) and use solar-metallicity tracks as appropriate for M 82 (paper I), reducing the impact of these uncertainties on our results. Furthermore, the maximum age from the CO bandheads would be most affected, and much less the minimum age so that the above uncertainties have little consequences on the age sequence suggested by our observational constraints. This result would be more affected if the stellar tracks significantly underpredict the number of very luminous but relatively cool stars ($`\mathrm{log}L_{\mathrm{bol}}>5.75`$, $`\mathrm{log}T_{\mathrm{eff}}<4.5`$), a possibility raised by Thornley et al. (2000) in their discussion of the Galactic Center. With correspondingly lower predicted neon ratios, the constraints might be more clearly accomodated with single bursts of somewhat longer durations. However, this point remains to be thoroughly investigated.

### 4.6 Low-mass End of the IMF

The low-mass IMF for the starburst population in M 82 has long been debated. Various authors argued that it is deficient in stars with initial masses below $`3\mathrm{M}_{}`$ compared to the solar neighbourhood IMF (Rieke et al., 1980, 1993; Bernlöhr, 1992; Doane & Mathews, 1993). The studies by Rieke and coworkers have been the most influential. Rieke et al. (1980) argued for a sharp truncation of the IMF at 3.5$`\mathrm{M}_{}`$, with no stars formed below this cutoff. Rieke et al. (1993) modified this conclusion and argued that while the IMF may extend down to 0.1$`\mathrm{M}_{}`$, the slope flattens significantly near 3$`\mathrm{M}_{}`$, at a substantially higher mass than in the solar neighbourhood (near 0.5$`\mathrm{M}_{}`$; see, e.g., Kroupa 2001 and references therein). More recently, Satyapal et al. (1997) challenged these hypotheses and modeled the integrated $`K`$-band luminosity of M 82 with a Salpeter IMF down to 0.1$`\mathrm{M}_{}`$ without using up more than $`35\%`$ of the total dynamical mass in the starburst, a reasonable upper limit.

The origin of the controversy resides primarily in the notable differences in observational constraints and model assumptions between the various studies. Here, we reexamine the issue of the low-mass IMF in M 82 in light of our new data and models. We consider two regions for this purpose: the nucleus (central 35 pc) and the starburst core (central 500 pc).

#### 4.6.1 Reexamination of the Low-mass IMF

Among our constraints, the stellar mass $`M^{}`$ is the only one sensitive to the low-mass star population. We computed this critical quantity by subtracting the mass of the gaseous component $`M_{\mathrm{gas}}`$ from the total dynamical mass $`M_{\mathrm{dyn}}`$ (see paper I for details). We derived the enclosed $`M_{\mathrm{dyn}}`$ as a function of projected radius from position-velocity maps of the millimeter CO $`J=10`$, \[Ne II\] 12.8 µm and \[S III\] 9069Å line emission (Shen & Lo, 1995; Achtermann & Lacy, 1995; McKeith et al., 1993) combined with the result of stellar velocity dispersion measurements at the nucleus based on the CO 2.29 µm bandhead by Gaffney, Lester, & Telesco (1993). We estimated $`M_{\mathrm{gas}}`$, accounting for molecular and ionized gas, from the CO $`J=10`$ map of Shen & Lo (1995) using the CO intensity to $`\mathrm{H}_2`$ column density conversion factor derived by Wild et al. (1992) from multi-line radiative transfer analysis and from the Lyman continuum luminosities. The $`M^{}`$ for the nucleus and the starburst core are $`7.9\times 10^7`$ and $`6.1\times 10^8\mathrm{M}_{}`$, respectively.

We investigated the low-mass IMF (shape and cutoff) by comparing the mass of stars predicted by the model bursts $`M_{\mathrm{stb}}^{}`$ with the $`M^{}`$ derived for each region. In this comparison, allowance has to be made for a possible stellar population predating the starburst, which can contribute to the measured stellar mass, so that the computed $`M_{\mathrm{stb}}^{}`$ does not exceed a plausible fraction of $`M^{}`$. To estimate this fraction, we followed the argument of McLeod et al. (1993). Tidal interaction with M 81 presumably induced the starburst in M 82; numerical simulations indicate that the perturbations resulting from such encounters can efficiently drive large amounts of gas towards the nuclear regions and that before the gas mass equals the mass in preexisting stars, self-gravitation in the gas should trigger a starburst with short timescale compared to the gas settling time (e.g., Hernquist, 1989; Mihos & Hernquist, 1994, 1996). Thus, $`M_{\mathrm{stb}}^{}`$ unlikely exceeds half of $`M_{\mathrm{dyn}}`$ in M 82. In our analysis, we adopted a conservative upper limit for $`M_{\mathrm{stb}}^{}`$ of 50% of the total stellar mass $`M^{}`$. We did not account for the gas mass returned to the ISM via stellar winds and supernova explosions, or expelled out of the galaxy in the starburst wind. Since massive stars contribute most to the mass return for the burst ages $`20\mathrm{Myr}`$ relevant here but do not dominate $`M_{\mathrm{stb}}^{}`$, these effects can be neglected to first order but will be further addressed below.

Our two-burst models for the nucleus and the starburst core (Tables 3 and 4) use only 9% and 33%, respectively, of the stellar masses for a Salpeter IMF slope between $`m_{\mathrm{low}}=1\mathrm{M}_{}`$ and $`m_{\mathrm{up}}=100\mathrm{M}_{}`$, but this choice of $`m_{\mathrm{low}}`$ was arbitrary. Figure 9 shows the burst masses required to reproduce the observed luminosities ($`L_K`$ and $`L_{\mathrm{Lyc}}`$) computed for $`m_{\mathrm{low}}`$ in the range $`0.15\mathrm{M}_{}`$, with a Salpeter slope and for the same burst ages and decay timescales. For the nucleus, the bursts never use up more than 30% of the measured $`M^{}`$. On the other hand, the burst mass for the starburst core reaches 50% of $`M^{}`$ at $`m_{\mathrm{low}}=0.4\mathrm{M}_{}`$ and equals $`M^{}`$ at $`m_{\mathrm{low}}0.1\mathrm{M}_{}`$. For the mass involved in the starburst not to exceed our adopted limit would require a sharp cutoff at 0.4$`\mathrm{M}_{}`$ for a Salpeter IMF. Alternatively, allowing the IMF to extend down to 0.1$`\mathrm{M}_{}`$, a flattening of the IMF at low masses would be required to reduce the increase in $`M_{\mathrm{stb}}^{}`$ with decreasing $`m_{\mathrm{low}}`$ so that $`M_{\mathrm{stb}}^{}50\%`$ of $`M^{}`$. With only one constraint available, our analysis does not allow us to break the degeneracy between the cutoff and the shape of the low-mass IMF.

#### 4.6.2 Factors Influencing the Low-mass IMF Determination

The $`M_{\mathrm{stb}}^{}`$ predicted by the models is determined by the absolute luminosities. This makes inferences on the low-mass IMF very sensitive to the extinction corrections applied, to the star formation history, and to the shape of the IMF at higher masses. The first point has been extensively discussed by McLeod et al. (1993), Rieke et al. (1993), and Satyapal et al. (1997). In that respect, our intrinsic luminosities correspond to the higher ranges reported in the literature and therefore impose more stringent constraints. In particular, the $`L_K`$ used by Satyapal et al. (1997) and the lower limit considered by Rieke et al. (1993) are both smaller than our value ($`5\times 10^8`$ and $`7.8\times 10^8\mathrm{L}_{}`$ compared to $`1.3\times 10^9\mathrm{L}_{}`$ for a $`\lambda =1.92.5\mu \mathrm{m}`$ $`K`$ bandpass and the photometric system of Wamsteker 1981 we adopted). With our two-burst models, they allow $`m_{\mathrm{low}}`$ down to 0.2$`\mathrm{M}_{}`$ or less, as shown in the bottom panel of Figure 9 (the mass constraints are similar between all three studies).

The star formation history can also significantly alter the conclusions about the low-mass IMF since longer or older bursts tend to use up more mass. In this regard, however, our conclusions above for a Salpeter IMF differ from those of Satyapal et al. (1997) mainly because of the different $`L_K`$ adopted. In terms of $`M_{\mathrm{stb}}^{}`$, their single burst model for the starburst core, with an age of 10 Myr and a short duration, is comparable to our two-burst model in which the old burst with an age of 9 Myr contributes 75% of the total burst mass for this region. Comparison with the models of Rieke et al. (1993) is less straightforward; our two-burst model is very similar to theirs with respect to burst ages and durations but the adopted IMFs are very different.

For the young ages of interest here, the various luminosities are dominated by high-mass stars, so that the relative slope of the IMF between different mass ranges is particularly critical in constraining the low-mass IMF. For instance, IMFs which steepen at high masses are less efficient at converting the mass into luminosity. In order to reproduce the luminosities, they require larger $`M_{\mathrm{stb}}^{}`$. To illustrate this, we considered three alternative IMFs: the Miller & Scalo (1979) IMF for the solar neighbourhood, the IMF proposed by Rieke et al. (1993) for M 82 (their “IMF 8”), and the IMF of Eisenhauer et al. (1998) for NGC 3603. Eisenhauer et al. derived the NGC 3603 IMF based on near-IR adaptive optics observations down to 1$`\mathrm{M}_{}`$, their detection limit. For our purposes, we allowed extensions below 1$`\mathrm{M}_{}`$ using the slope determined in the $`115\mathrm{M}_{}`$ range. Table 5 gives the parameters of the IMFs (see also Figure 8) and Figure 9 shows the variations of the burst masses versus $`m_{\mathrm{low}}`$ computed with each of them.

For the range $`0.1100\mathrm{M}_{}`$, the mass fraction locked in stars less massive than 1$`\mathrm{M}_{}`$ is significantly smaller for the Miller & Scalo IMF compared to the Salpeter IMF (44% and 61%), while the corresponding fractions for stars more massive than 10$`\mathrm{M}_{}`$ are nearly identical (8% and 12%, respectively). However, properties such as $`L_K`$, $`L_{\mathrm{bol}}`$, and $`L_{\mathrm{Lyc}}`$ are steeply increasing functions of the stellar mass and the Miller & Scalo IMF is much steeper at the high-mass end than the Salpeter IMF. Consequently, normalizing to the same luminosities, the Miller & Scalo IMF predicts substantially more mass than the Salpeter IMF, with $`M_{\mathrm{stb}}^{}`$ larger by up to a factor of two in our models. The Rieke et al. (1993) IMF is very similar to the Miller & Scalo IMF but with the first inflection point moved up from 1$`\mathrm{M}_{}`$, giving particular emphasis on the intermediate-mass range. Of the various IMFs they explored, such modifications to solar-neighbourhood type IMFs were the most successful at reproducing their set of constraints. Finally, compared to the Rieke et al. (1993) IMF, that of Eisenhauer et al. (1998) emphasizes less the intermediate-mass range but has a larger proportion of stars above 10$`\mathrm{M}_{}`$ (40% versus 23% of the total mass between 1 and 100$`\mathrm{M}_{}`$), implying much higher mass-to-luminosity conversion efficiency.

Figure 9 shows that for the central 35 pc of M 82, the low-mass IMF remains unconstrained for any of the three alternative IMFs: the computed $`M_{\mathrm{stb}}^{}`$ never exceeds 35% of the measured $`M^{}`$ with $`m_{\mathrm{low}}`$ down to 0.1$`\mathrm{M}_{}`$. The case of the starburst core is again more restrictive. With the Miller & Scalo (1979) IMF, more than 50% of the stellar mass is involved in the bursts for $`m_{\mathrm{low}}2\mathrm{M}_{}`$. If we allow for stars down to 0.1$`\mathrm{M}_{}`$ to form, we conclude as Rieke et al. (1993) that a displacement of the first inflection point to higher masses is then needed to effectively reduce $`M_{\mathrm{stb}}^{}`$. With our data and models, though, $`M_{\mathrm{stb}}^{}`$ still exceeds half of $`M^{}`$ at $`m_{\mathrm{low}}1\mathrm{M}_{}`$ if we adopt the Rieke et al. IMF. The Eisenhauer et al. (1998) IMF allows us to extend the IMF with the same $`\alpha =1.73`$ power-law index from 15$`\mathrm{M}_{}`$ down to 0.1$`\mathrm{M}_{}`$ without using up more than 30% of the stellar mass in the bursts. Interestingly, if we assume both the IMF and the $`K`$-band luminosity of Rieke et al. (1993), the $`M_{\mathrm{stb}}^{}`$ remains well below 50% of $`M_{}`$ over the entire range of $`m_{\mathrm{low}}`$ (dashed line in the bottom panel of Figure 9).

Another consideration worth pointing out is the mass loss via massive stars winds and supernova explosions, which is important already within a few million years of the onset of star formation (see, e.g., Leitherer et al., 1999). In particular, the mass returned to the ISM could have been entrained out of M 82 in the starburst wind, in which case our $`M^{}`$ constraint would be missing some fraction or, alternatively, the predicted $`M_{\mathrm{stb}}^{}`$ curves of Figure 9 would decrease accordingly. An estimate of the effect is difficult because of the unknown amount of returned mass actually driven away in the starburst wind. As an illustration, we may consider the extreme case in which massive stars $`10\mathrm{M}_{}`$ lose all their mass within $`20\mathrm{Myr}`$, and that all ejectae are carried out of M 82 by the starburst wind. For a Salpeter IMF extending from 1 to 100$`\mathrm{M}_{}`$, the range $`10\mathrm{M}_{}`$ accounts for 31% of the total mass so that at most $`30\%`$ of $`M_{\mathrm{stb}}^{}`$ would be lost and expelled. This fraction varies with $`m_{\mathrm{low}}`$, up to $`50\%`$ at 3$`\mathrm{M}_{}`$ and down to $`10\%`$ at 0.1$`\mathrm{M}_{}`$, and more generally depends on all IMF parameters. The overall impact is to relax the requirements on the low-mass IMF but the effect is not likely to be more important than the other factors discussed above.

#### 4.6.3 Additional Remarks

In view of the uncertainties involved and the lack of strong constraints on the shape of the IMF, the low-mass IMF in M 82 remains elusive and only general statements can be made. If stars formed down to 0.1$`\mathrm{M}_{}`$, our data and models may indicate a substantial flattening of the starburst IMF below a few $`\mathrm{M}_{}`$, assuming a Salpeter slope at higher masses. This requirement would be alleviated for a flatter intermediate- and high-mass IMF or for less conservative limits on $`M_{\mathrm{stb}}^{}`$. A bias against the formation of low-mass stars has been suggested for other starburst galaxies as well (e.g., Augarde & Lequeux, 1985; Wright et al., 1988; Olofsson, 1989; Nakagawa et al., 1989; Prestwich, Joseph, & Wright, 1994). Engelbracht et al. (1996, 1998) concluded for NGC 6946 and NGC 253 that the IMF proposed by Rieke et al. (1993) for M 82 is more compatible with the observations than a solar-neighbourhood IMF, while Alonso-Herrero et al. (2001) found that the IMF in NGC 1614 may be even more biased towards high-mass stars.

Kinematic mass-to-light ratios have been obtained for a few individual luminous stellar clusters in M 82, from velocity dispersion measurements combined with high resolution Hubble Space Telescope (HST) imaging. Using optical data obtained for the super star cluster “M 82-F” $`450\mathrm{pc}`$ southwest of the nucleus and an age estimate of 60 Myr, Smith & Gallagher (2001) concluded that the cluster IMF is top-heavy and requires a low-mass cutoff at $`23\mathrm{M}_{}`$ for an $`\alpha =2.3`$ slope. McCrady, Gilbert, & Graham (2003) studied two younger ($`10^7\mathrm{yr}`$) clusters roughly 180 pc southwest of the nucleus based on near-IR observations. For an IMF with $`\alpha =1.3`$ at $`0.10.5\mathrm{M}_{}`$ and $`\alpha =2.3`$ at $`0.5100\mathrm{M}_{}`$ (based on Kroupa, 2001), their data support a relative deficiency at $`1\mathrm{M}_{}`$ for one cluster but allow $`m_{\mathrm{low}}`$ down to 0.1$`\mathrm{M}_{}`$ for the other. In their revision of the Smith & Gallagher (2001) analysis of M 82-F, using the same mass but combined with the near-IR luminosity and a lower age of 40 Myr, they confirmed the evidence for a top-heavy IMF although with a looser requirement of $`m_{\mathrm{low}}1\mathrm{M}_{}`$. Similar investigations of young super star clusters in other systems, including NGC 1569A and NGC 1701-1 (Sternberg, 1998), and five clusters in NGC 4038/4039 (Mengel et al., 2002), yielded mixed results as well. This seems to indicate a range of possible IMFs in super star clusters which may be related to their birth environment or to mass segregation (see, e.g., Mengel et al., 2002; McCrady et al., 2003). The results of the M 82 clusters may provide support in favour of an IMF flattening at low masses but, as pointed out by Smith & Gallagher (2001), it is unclear whether conclusions about luminous star clusters, let alone of a handful only, can be extended to the overall pattern of star formation in the starburst.

Direct censuses of the low-mass star populations in Galactic and near-extragalactic high-mass star-forming regions are still scarce and difficult owing to the faintness of such stars and to important crowding effects. Recent results from very high sensitivity and angular resolution observations suggest that mere truncations of the IMF near 1$`\mathrm{M}_{}`$ are not plausible, but the shape of the low-mass IMF remains poorly constrained. In particular, Brandl et al. (1999) found that the central parsec of NGC 3603 is well populated by pre-main sequence stars down to 0.1$`\mathrm{M}_{}`$, based on data from the Very Large Telescope. The near-IR adaptive optics observations by Eisenhauer et al. (1998) show no evidence for a significant flattening of the IMF down to their detection limit of 1$`\mathrm{M}_{}`$. In contrast, from deep HST images of the R 136 cluster in 30 Doradus, Sirianni et al. (1998, 2000) could probe the $`1.356.5\mathrm{M}_{}`$ range and reported a flattening of the IMF below 2$`\mathrm{M}_{}`$, a result however disputed by Selman et al. (1999). Given the shape of the solar-neighbourhood IMF, with possible turnover near 0.5$`\mathrm{M}_{}`$ (e.g., Scalo, 1986, 1998; Rana, 1987; Kroupa, 2001), and indications of a flattening towards lower masses in some local templates of starburst regions and super star clusters, our inferences for M 82 may not need be interpreted in terms of an “abnormal” deficiency of the IMF in low-mass stars.

## 5 SPATIALLY DETAILED MODELING OF M 82

The modeling of selected regions in the previous section provides the general characteristics of starburst activity in the central regions of M 82, but only an incomplete picture of its evolution. In particular, our choice of the small-scale regions as those with brightest continuum or line emission, and with deepest or shallowest CO bandheads likely introduced a bias towards preferential ages. In order to obtain a more complete picture, we modeled individual regions throughout the entire 3D field of view. These have a range in their properties suggesting a range in evolutionary states and burst strengths. Although along any line of sight the integrated properties will always be dominated by the most luminous populations, a detailed modeling on small spatial scales reduces the bias towards particular ages.

Given our results of § 4, we assumed that a Salpeter IMF with high upper mass limits and burst decay timescales $`5\mathrm{Myr}`$ are appropriate for individual pixels. We also assumed that the near-IR properties are dominated everywhere by red supergiants. We emphasize again that the analysis of the $`M^{}/L_K`$ ratio of paper I strongly supports this hypothesis for the regions where the CO bandheads alone do not allow us to constrain the luminosity class of the evolved stars.

### 5.1 Spatial Distribution of the Burst Ages and Strengths

Figure 10 compares the data with single-burst models for $`m_{\mathrm{up}}=100\mathrm{M}_{}`$, and $`t_{\mathrm{sc}}=1`$ and 5 Myr, in the diagrams of neon, $`L_{\mathrm{bol}}`$/$`L_{\mathrm{Lyc}}`$, and $`L_K`$/$`L_{\mathrm{Lyc}}`$ ratios versus $`W_{2.29}`$. The size of the data points is proportional to the intrinsic $`L_K`$. We did not account for the contribution of the cool evolved stars to the bolometric luminosity of individual pixels, so that the values of $`L_{\mathrm{bol}}^{\mathrm{OB}}`$/$`L_{\mathrm{Lyc}}`$ are plotted in the $`L_{\mathrm{bol}}`$/$`L_{\mathrm{Lyc}}`$ diagram.

A large fraction of the data points in all three diagrams do not fall on the model curves and exhibit the typical behaviour encountered for the selected regions: the neon and $`L_{\mathrm{bol}}^{\mathrm{OB}}`$/$`L_{\mathrm{Lyc}}`$ ratios imply significantly younger ages than $`W_{2.29}`$, with $`L_K`$/$`L_{\mathrm{Lyc}}`$ corresponding to intermediate ages. In the $`L_K`$/$`L_{\mathrm{Lyc}}`$ versus $`W_{2.29}`$ diagram, the data points are distributed along a path that is remarkably parallel to the model curve for $`t_{\mathrm{sc}}=1\mathrm{Myr}`$, but displaced towards lower $`L_K`$/$`L_{\mathrm{Lyc}}`$ by about an order of magnitude. Although the properties of some of the individual pixels could be reconciled with a single burst if the timescale is increased to values approaching 5 Myr, we have pursued with our hypothesis of sequential bursts. The data can be reproduced with two short successive bursts, with the young bursts producing about ten times more ionizing luminosity than the older bursts dominating the near-IR luminosity.

For each pixel, we therefore fitted a two-burst model to the constraints, each with $`t_{\mathrm{sc}}=1\mathrm{Myr}`$, in the same manner as for the selected regions. We assigned an age of 1 Myr for $`L_{\mathrm{bol}}^{\mathrm{OB}}`$/$`L_{\mathrm{Lyc}}`$ ratios lower than the minimum $`L_{\mathrm{bol}}`$/$`L_{\mathrm{Lyc}}`$ predicted by the models (this is the case for a few pixels only). For the old bursts, we adopted the youngest solutions inferred from $`W_{2.29}`$ and assigned the age of maximum $`W_{2.29}`$ to the few pixels with EWs exceeding this limit. We estimated the formal uncertainties on the burst ages and initial star formation rates $`R_0`$ by varying in turn the observational constraints in the ranges given by the measurements uncertainties.

Figure 11 presents the results. The ages for the young bursts range from 3.5 to 6.6 Myr, with an average of 4.8 Myr and dispersion of $`\sigma =0.5\mathrm{Myr}`$. The ages for the old bursts range from 7.2 to 12.6 Myr (the upper limit from the maximum $`W_{2.29}`$), with an average of 8.8 Myr and dispersion of $`\sigma =0.9\mathrm{Myr}`$. Figure 12 shows the average burst age and total burst strength as a function of projected distance from the nucleus in a slit 3″–wide along the galactic plane of M 82, as indicated in Figure 11. The typical uncertainties are $`\pm 0.51\mathrm{Myr}`$ for the ages and about 30% for the strengths. In view of these uncertainties, we conclude that the maps and profiles indicate roughly constant ages among the young bursts and among the old bursts, as well as nearly uniform strengths for the young bursts. Spatial variations in strength among the old bursts are more significant with, in particular, a clear peak around the nucleus consistent with our findings from the modeling of selected regions (§ 4.4).

Satyapal et al. (1997) investigated the spatial variations in ages from a dozen compact $`K`$-band emission sources, most of which are located within the regions mapped with 3D. They inferred burst ages in the range $`410\mathrm{Myr}`$, with a trend of younger ages at larger projected radii. Given that they modeled the sources as single bursts and used the strength of the CO bandheads longwards of 2.29 µm together with the EW of the Br$`\gamma `$ emission line at 2.17 µm (inversely proportional to $`L_K`$/$`L_{\mathrm{Lyc}}`$), their results are consistent with ours.

### 5.2 The Global Star Formation History

From our spatially detailed modeling, we reconstructed the global star formation history within the 3D field of view by integrating the initial star formation rate $`R_0`$ as a function of burst age $`t_\mathrm{b}`$ over all individual pixels. For sufficient sampling, we chose age bins increasing logarithmically by 0.05 dex. As long as the timescales assumed for all individual bursts are similar and relatively short compared to the overall star formation history of interest, the integrated $`R_0(t_\mathrm{b})`$ curve gives a measure of the “instantaneous” star formation rate at different times. Figure 13 shows the derived integrated starburst history together with the surface density of star formation rate, computed by dividing the integrated $`R_0(t_\mathrm{b})`$ by the total area of the pixels contributing to each age bin. The results for the individual pixels are reproduced as well.

The integrated star formation rate exhibits two conspicuous peaks near 5 and 9 Myr, corresponding to distinct global starburst episodes. The first episode was globally 2.5 times stronger than the most recent one in terms of mass of stars formed. The star formation rate surface density reflects the nearly uniform strengths for the young bursts and the larger variations for the old bursts seen in Figures 11 and 12. We note that modeling the individual pixels with two short successive bursts introduces a bias towards separated burst events in the integrated star formation history. Again, such a scenario is however plausible in view of the spatial anticorrelation between the tracers of OB stars and those of cool evolved stars (§ 4.1).

The integrated star formation rate in the 3D field of view is well reproduced by two Gaussians in $`R_0`$ versus age centered at 4.7 and 8.9 Myr, with amplitudes of 6.3 and $`18.5\mathrm{M}_{}\mathrm{yr}^1`$ and full-widths at half-maximum of 1.1 and 1.7 Myr, respectively. These Gaussians are shown in Figure 13 along with the two-burst model fit to the integrated properties of the 3D field of view from § 4. The exponentially decaying functions provide reasonable though crude approximations to the detailed star formation history. In view of our results, it may not be surprising that Rieke et al. (1993) found that a double Gaussian profile optimized their model fits for the entire starburst core (with one of their best models reaching peak star formation activity near 5 and 13 Myr).

## 6 THE NATURE OF STARBURST ACTIVITY IN M 82

### 6.1 Star Formation Process

Observations of M 82 reveal important small-scale structure tracing individual burst sites on 10-pc or even parsec scales (see paper I). In particular, the ionized and molecular gas as well as the IR-emitting dust exhibit clumpy morphologies on scales at least as small as 15 pc (e.g., Telesco & Gezari, 1992; Larkin et al., 1994; Achtermann & Lacy, 1995; Shen & Lo, 1995). Compact continuum sources tracing clusters of red supergiants near their maximum luminosity are detected at near-IR wavelengths (e.g., Satyapal et al., 1997; McCrady et al., 2003). Optical imaging with HST has resolved over a hundred young super star clusters, some of which may lie in directions of lower extinction and belong to the inner starburst regions (O’Connell et al., 1995).

Our models show that the properties of individual burst sites are consistent with the formation of very massive stars and with short burst timescales of a few million years or less (see also Satyapal et al., 1997, for selected clusters). The typical sizes and the star formation parameters inferred make these burst sites comparable to Galactic and near-extragalactic massive star-forming regions such as NGC 3603 and the R 136 cluster in 30 Doradus.

Interestingly, the smoother and low-surface brightness emission regions make a substantial contribution to the integrated properties, suggesting that the prominent sub-structure may in fact trace the largest and most luminous burst sites. The diffuse ionized gas emission (notably in H recombination lines) and thermal dust emission represent about 50% of the respective total emission from the starburst core. The smoother and fainter stellar $`K`$-band continuum emission contributes at least as large a fraction to the total $`K`$-band luminosity. These components correlate roughly with the overall distribution of the brighter and more compact sources, and partly break up into smaller-scale structure at higher angular resolution (see, e.g., the HST NICMOS maps of Alonso-Herrero et al., 2003). Whether the diffuse emission in M 82, and more generally in starburst galaxies, traces large numbers of unresolved fainter compact sources or truly extended sources is still an open question, with implications for the importance of cluster formation in starbursts and subsequent dissolution in the stellar field (e.g., Ho & Fillipenko, 1996; Meurer et al., 1995; Tremonti et al., 2001).

### 6.2 Evolution of Starburst Activity

#### 6.2.1 Key Morphological Features

For a plausible evolutionary scenario of starburst activity in M 82, our results must be interpreted together with the key morphological features:
$``$ Large-scale tails and bridges of material emanating from M 82 and connecting with its massive neighbour M 81 located about 36 kpc in projection indicate gravitational interaction between the two galaxies (e.g., Yun, Ho, & Lo, 1993, 1994).
$``$ The morphology of the $`K`$-band emission, dominated by red supergiants, is highly suggestive of a nearly edge-on disk, peaks strongly towards the nucleus, and provides the strongest evidence for a $`1\mathrm{kpc}`$–long stellar bar (e.g., Telesco et al., 1991; McLeod et al., 1993; Larkin et al., 1994). The brightest and most compact sources are mainly located in an “inner disk” within $`150\mathrm{pc}`$ of the nucleus and at the ends of the bar (e.g., Satyapal et al., 1997).
$``$ The most prominent H II regions are concentrated in a rotating ring-like structure of radius $`85\mathrm{pc}`$ and, outside this ring, along ridges presumably on the leading edge of the rotating stellar bar (e.g., Larkin et al., 1994; Achtermann & Lacy, 1995). The “ionized ring” is located just inside of the main concentrations of molecular gas in a toroid or tightly-wound spiral arms (e.g., Shen & Lo, 1995). Our Br$`\gamma `$ source B2 coincides in projection with the western edge of the ionized ring while B1 is located farther out.
$``$ The important population of young supernova remnants (SNRs) detected at radio wavelengths is quite uniformly distributed along the galactic plane over $`600\mathrm{pc}`$ (e.g., Kronberg, Biermann, & Schwab, 1985; Muxlow et al., 1994). Their spatial distribution together with the ages derived for each starburst episode and the IMF weighting are more consistent with their being associated primarily with the populations of red supergiants throughout the disk.
$``$ A bipolar outflow along the minor axis of M 82 is traced out to at least 5 kpc by X-ray and optical observations (e.g., Bregman, Schulman, & Tomisaka, 1995; Shopbell & Bland-Hawthorn, 1998). The \[Fe II\] 1.644 µm emission exhibits a prominent arc-like structure to the south and centered near the nucleus. It may trace shock-ionized iron-enriched gas where the outflowing wind interacts with interstellar disk material (e.g., Greenhouse et al., 1997; Alonso-Herrero et al., 2003).
$``$ There is an apparent lack of gas and dust within the central few tens of parsecs of M 82 (e.g., Telesco & Gezari, 1992; Larkin et al., 1994; Shen & Lo, 1995; Achtermann & Lacy, 1995; Seaquist et al., 1996).

Given the relative distributions of the gaseous and stellar components and since M 82 is viewed at an inclination angle of $`80^{}`$, the two sequential bursts inferred locally and globally within the central starburst regions could be understood in terms of projection effects. In particular, the inner few tens of parsecs harbored the most intense starburst activity about 10 Myr ago and while more recent bursts appear superposed, they took place predominantly in circumnuclear regions at larger radii. The geometrical picture of M 82 further emphasizes the particular intensity of the old nuclear burst revealed in Figures 11 and 12, with the red supergiants within the central few tens of parsecs in projection being physically close to the nucleus along the line of sight as well, while the OB stars — and the red supergiants elsewhere — possibly more extended along the line of sight. While alternatives could be advanced, we chose, in line with our two-burst hypothesis, to focus on the above interpretation as an illustrative yet possible case.

#### 6.2.2 Triggering Mechanism

The primary triggering mechanism for starburst activity in M 82 is generally attributed to the $`\mathrm{M}\mathrm{\hspace{0.17em}82}\mathrm{M}\mathrm{\hspace{0.17em}81}`$ tidal interaction $`10^8\mathrm{yr}`$ ago (e.g., Gottesman & Weliachew, 1977; O’Connell & Mangano, 1978; Lo et al., 1987; Yun et al., 1993, 1994). In this scenario, the ISM in M 82 experienced strong large-scale torques and loss of angular momentum as it was transported towards the dynamical center of the galaxy, in accordance with numerical simulations (e.g., Sundelius et al., 1987; Noguchi, 1987, 1988; Mihos & Hernquist, 1996).

The increased cloud-cloud collision rate in the disk and the large amounts of material accumulated and compressed in the innermost regions could have led to the first starburst episode, characterized by very intense star formation at the nucleus and lower-level activity elsewhere in the disk. The nuclear burst and subsequent high rate of supernova explosions consumed rapidly the gas supply and expelled the remaining gas via the starburst wind, thereby creating the central cavity in the ISM and preventing further star formation. The absence of a large concentration of radio SNRs in the nuclear vicinity is consistent with the lack of gas and dust.

#### 6.2.3 Subsequent Evolution

Numerical simulations show that bars in galactic disks can be induced by galactic interactions and are very effective at driving material towards the central regions of galaxies (e.g., Combes & Gérin, 1985; Noguchi, 1987, 1988; Shlosman et al., 1989; Athanassoula, 1992). Ring- or spiral-like dynamical resonances may develop under the action of such non-axisymmetric perturbations. In the presence of inner Lindblad resonances (ILRs), the radial inflow of material is stopped before it reaches the nucleus of the galaxy and accumulates in a circumnuclear ring. Star formation is triggered by shocks in the ring and along ridges leading the bar, and may be particularly enhanced at their intersections.

The second starburst episode in our models of M 82 could thus be attributed to the presence of the bar and accompanying ILRs indicated by the global gas distributions (see also Lo et al., 1987; Telesco et al., 1991; Achtermann & Lacy, 1995; Neininger et al., 1998). The ionized ridges presumably trace star formation sites on the leading side of the bar. Enhanced activity where the ISM streaming along the bar meets with the ionized ring may be hinted at by the slight increase in strength of the young bursts near B2 suggested by Figures 11 and 12 (see also Achtermann & Lacy, 1995). The bar may have played a role during the first starburst episode by channeling the nuclear inflow before the present dynamical resonances appeared.

The apparent outward progression of starburst activity suggested by the sequence Nuclear supergiants $``$ Ionized ring $``$ Molecular ring may reflect the temporal development of the bar and redistribution of the mass as suggested by Telesco et al. (1991). At the radius of the ionized ring, the rotational velocity is about $`120\mathrm{km}\mathrm{s}^1`$ (Achtermann & Lacy, 1995; Shen & Lo, 1995). This implies an orbital period of $`45\mathrm{Myr}`$, i.e. a dynamical timescale similar to the time elapsed between the peak star formation activity during the two recent starburst episodes.

An alternative explanation could be self-induced propagation radially outward in the disk, as a consequence of the powerful nuclear burst and expanding shock wave generated by the massive stars winds and supernova explosions (e.g., Carlstrom & Kronberg, 1991; Achtermann & Lacy, 1995). Without dynamical resonances, the ISM may be expected to follow a smooth, radially decreasing density profile. Star formation triggered by an expanding “superbubble” in such a medium seems difficult to reconcile with pronounced and narrow peaks in the global star formation rate. The two peaks could perhaps reflect a sudden increase of mechanical energy release by the older nuclear burst after about 5 Myr but the contribution by massive star winds is already very substantial at ages $`<5\mathrm{Myr}`$ before the first supernovae explode (e.g., Leitherer et al., 1999). Star formation activity along the ridges beyond the ionized ring is also difficult to explain in the self-induced outward propagation scenario.

We favoured bar-driven evolution as dominant mechanism because the stellar and gaseous components exhibit morphological features strongly supporting the existence of a bar and its dynamical effects. Bar-driven evolution has been proposed for other starburst systems as well as more quiescent spiral galaxies in which circumnuclear rings of enhanced gas density and star formation activity coexist with stellar or gaseous bars (e.g., Telesco & Decher, 1988; Telesco, Dressel, & Wolstencroft, 1993; Knapen et al., 1995; Böker, Förster Schreiber, & Genzel, 1997).

#### 6.2.4 Global Picture

Investigations of regions outside of the central active starburst core outline further key elements for a global scenario. The integrated optical spectrum of M 82 exhibits strong Balmer line transitions in absorption typical of A-type stars, indicative of an intermediate-age ($`0.11\mathrm{Gyr}`$) “post-starburst” population (Kennicutt, 1992). de Grijs et al. (2000, 2001) used HST optical and near-IR data to study post-starburst regions $`0.51\mathrm{kpc}`$ northeast of the nucleus, identified as such from strong Balmer absorption lines and discontinuity in their spectra (O’Connell & Mangano, 1978; Marcum & O’Connell, 1996). The results of de Grijs and coworkers based on evolutionary synthesis modeling of the photometric properties of compact luminous clusters and bright giant stars as well as their detection of candidate H$`\alpha `$ SNRs imply that enhanced star formation occurred $`600\mathrm{Myr}`$ ago at radii of $`1\mathrm{kpc}`$ and up to $`30\mathrm{Myr}`$ ago near $`500\mathrm{pc}`$. Gallagher & Smith (1999) modeled high-quality optical spectra of two super star clusters about 450 pc southwest of the nucleus and derived ages of $`60\mathrm{Myr}`$.

In the framework of our models and interpretation, together with the spatial distribution of the radio SNRs and the $`M^{}/L_K`$ ratio, the first starburst episode $`10\mathrm{Myr}`$ ago was not strictly confined to the nuclear vicinity, but occurred throughout the central 500 pc as well. Along with the evidence from the post-starburst regions, this would support an “outside-inside-out” propagation of starburst activity in M 82. The scenario we propose is summarized here and sketched in Figure 14.

$``$ Following the gravitational interaction between M 82 and M 81, $`10^8\mathrm{yr}`$ ago, the ISM in M 82 experienced strong large-scale torques, loss of angular momentum, and important infall towards the nuclear regions, leading to enhanced star formation activity in the central kiloparsec. The stellar bar induced by the interaction possibly played a role in channeling the inflow. In the absence of dynamical resonances, the infalling material can reach the nucleus.

$``$ Within the central 500 pc of M 82, a first starburst episode took place $`815\mathrm{Myr}`$ ago and was rapidly exhausted. The inner few tens of parsecs at the nucleus hosted the most intense star formation activity.

$``$ A subsequent starburst episode was triggered predominantly by bar-induced dynamical resonances. It occurred $`46\mathrm{Myr}`$ ago, mainly in a circumnuclear ring and along the stellar bar, and also decayed rapidly.

$``$ A supernova-driven starburst wind originating in the center of M 82 has broken out of the galactic plane, the dramatic aftermath of the powerful nuclear burst. The outflow component in the disk may have played some role in triggering the second starburst episode.

### 6.3 Timescales and Feedback Effects

The global durations of $`10^6\mathrm{yr}`$ derived for each of the recent starburst episodes in M 82 seem surprising for spatial scales up to at least 500 pc and in view of the longer estimates $`10^710^8\mathrm{yr}`$ often quoted for starburst systems (e.g., Thronson & Telesco, 1986; Heckman, 1998). We consider below three timescale arguments in trying to gain more insight into the quenching mechanisms of the recent starbursts in M 82.

One argument relies on the comparison between the present star formation rate and the mass of the current gas reservoir, giving a “gas consumption timescale” $`\tau _{\mathrm{gas}}`$. In M 82, the gas is mostly in molecular form, with a total $`M_{\mathrm{H}_2}=1.8\times 10^8\mathrm{M}_{}`$ (e.g., Wild et al., 1992). From our two-burst model of the starburst core (Table 4) and with a Salpeter IMF, $`M_{\mathrm{stb}}^{}2\times 10^8\mathrm{M}_{}`$ of stars were formed in the past $`15\mathrm{Myr}`$ for $`m_{\mathrm{low}}=1\mathrm{M}_{}`$ and $`5\times 10^8\mathrm{M}_{}`$ for $`m_{\mathrm{low}}=0.1\mathrm{M}_{}`$. The latter value represents 80% of the measured stellar mass and, in view of the discussion of § 4.6, is here considered as an upper limit. The average star formation rate in the recent history of M 82 is then $`R=1333\mathrm{M}_{}\mathrm{yr}^1`$, depending on $`m_{\mathrm{low}}`$. Assuming 100% efficiency, starburst activity in M 82 at the level observed in the past 15 Myr could be sustained for another $`\tau _{\mathrm{gas}}515\mathrm{Myr}`$. This provides a fair estimate for the minimum gas depletion timescale but has no connection with the physical processes responsible for the evolution of the starburst.

Considerations based on dynamical timescales are more closely linked to the large-scale mechanisms that may drive starburst activity. For a sample of 36 IR-selected starburst galaxies, Kennicutt (1998) inferred a median gas conversion efficiency of 30% per dynamical timescale, taken as an average orbital period for all objects of $`10^8\mathrm{yr}`$; $`ϵ_{\mathrm{gas}}=100\%/10^8\mathrm{yr}`$ occurs but in the most extreme objects. From the rotation curve of M 82 (e.g., Götz et al., 1990; McKeith et al., 1993; Achtermann & Lacy, 1995; Shen & Lo, 1995; Neininger et al., 1998), we compute shorter dynamical timescales for regions inside the starburst core. For example, the orbital periods are about 5, 10, and 20 Myr at radii of 90, 225, and 500 pc corresponding to the locations of the ionized ring, of the main concentrations of molecular gas, and of the outer limit enclosing the starburst core and most of the gas. The dynamics allow then for burst timescales of approximately $`\tau _{\mathrm{dyn}}=15\mathrm{Myr}`$, $`30\mathrm{Myr}`$, and $`65\mathrm{Myr}`$ for 30% efficiency per orbital period. Even with 100% efficiency, the timescales are barely short enough to explain the durations of each starburst episode in M 82.

Such dynamical considerations may provide more insight concerning the settling time of the gas and the development of dynamical resonances. The orbital time at the radius of the ionized ring is similar to the time separation between the two most recent bursts. At radii $`13\mathrm{kpc}`$, the orbital period is $`100200\mathrm{Myr}`$, which is of the order of the age of the super star clusters in the post-starburst regions $`1\mathrm{kpc}`$ northeast of the nucleus and to the epoch of last peri-passage between M 82 and M 81. The scenario proposed in § 6.2 for the triggering and propagation of the starburst episodes in M 82 is consistent with these timescales.

None of the above arguments is however completely satisfactory for explaining the very short global burst durations suggested by our models. Most importantly, any feedback effect of massive stars via strong stellar winds and subsequent supernova explosions is neglected. Qualitatively, timescales of a few million years are similar to the lifetimes of stars with initial masses $`50\mathrm{M}_{}`$. Although both episodes took place over $`100500\mathrm{pc}`$ scales, each was triggered nearly simultaneously everywhere, presumably as a result of brief gas compression events. Massive star winds and supernova explosions rapidly disrupting the remaining ISM locally, preventing further star formation after a few million years, may have “conspired” to produce a rapid decay of the global starburst activity as well.

We can make simple estimates for the “feedback timescale” $`\tau _{\mathrm{feed}}`$ by comparing the cumulative mechanical energy $`E_{\mathrm{mech}}`$ injected by supernovae (and neglecting the contribution by stellar winds) into the ISM with its gravitational binding energy $`E_{\mathrm{bind}}`$, and asssuming that starburst activity stops when they exactly balance each other:

$$\dot{E}_{\mathrm{mech}}\tau _{\mathrm{feed}}=E_{\mathrm{bind}},$$
(1)
where $`\dot{E}_{\mathrm{mech}}`$ is the rate of mechanical energy deposition. We assumed that all stars with initial masses $`8\mathrm{M}_{}`$ explode as supernovae with a typical energy of $`10^{51}\mathrm{erg}`$, of which a fraction $`\eta `$ is transferred as kinetic energy to the ISM. For simplicity, we considered a spherically symmetric distribution of gas and stars within the starburst core, so that $`E_{\mathrm{bind}}=\left(GM_{\mathrm{dyn}}M_{\mathrm{gas}}\right)/`$ where $`G`$ is the gravitational constant, $`M_{\mathrm{dyn}}`$ and $`M_{\mathrm{gas}}`$ are the dynamical and gas masses, and $``$ is the radius of the region.

The rate of supernova explosions $`\nu _{\mathrm{SN}}`$ is the most direct tracer of $`\dot{E}_{\mathrm{mech}}`$ but the observed rates suffer from rather large uncertainties. We computed $`\dot{E}_{\mathrm{mech}}`$ using the inferred star formation rate $`R`$ instead. The implicit assumption is that the burst is sufficiently evolved to have non-negligible $`\nu _{\mathrm{SN}}`$ and that the star formation rate is still substantial. For a Salpeter IMF between 1$`\mathrm{M}_{}`$ and 100$`\mathrm{M}_{}`$, $`\nu _{\mathrm{SN}}[\mathrm{yr}^1]0.02R[\mathrm{M}_{}\mathrm{yr}^1]`$. With these assumptions, Eq. 1 can be written:

$$\frac{\tau _{\mathrm{feed}}}{\mathrm{Myr}}500\left(\frac{f}{\eta }\right)\left(\frac{M_{\mathrm{dyn}}}{10^{10}\mathrm{M}_{}}\right)^2\left(\frac{}{\mathrm{kpc}}\right)^1\left(\frac{R}{\mathrm{M}_{}\mathrm{yr}^1}\right)^1,$$
(2)
where $`f`$ is the gas mass fraction. For the starburst core of M 82, $`=225\mathrm{pc}`$, $`M_{\mathrm{dyn}}8\times 10^8\mathrm{M}_{}`$, and $`M_{\mathrm{gas}}1.8\times 10^8\mathrm{M}_{}`$ (see paper I). With the average $`R13\mathrm{M}_{}\mathrm{yr}^1`$ derived for the appropriate IMF parameters, Eq. 2 yields $`\tau _{\mathrm{feed}}0.2\eta ^1\mathrm{Myr}`$. The supernova efficiency is poorly constrained, but values in the range $`\eta =0.11`$ seem appropriate for M 82 (Chevalier & Clegg, 1985). Consequently, $`\tau _{\mathrm{feed}}0.22\mathrm{Myr}`$. For comparison, with $`\nu _{\mathrm{SN}}0.06\mathrm{yr}^1`$ for M 82, we obtain $`\tau _{\mathrm{feed}}110\mathrm{Myr}`$ depending on $`\eta `$. As alternative tracer, we can use $`L_{\mathrm{bol}}`$. In a steady-state, our starburst models imply $`L_{\mathrm{bol}}[\mathrm{L}_{}]10^{12}\nu _{\mathrm{SN}}[\mathrm{yr}^1]`$, and thus $`R[\mathrm{M}_{}\mathrm{yr}^1]0.5\left(L_{\mathrm{bol}}/10^{10}\mathrm{L}_{}\right)`$. With $`L_{\mathrm{bol}}=6.6\times 10^{10}\mathrm{L}_{}`$ for the starburst core, Eq. 2 implies $`\tau _{\mathrm{feed}}110\mathrm{Myr}`$ as well. These estimates represent upper limits because we neglected massive stars winds.

The $`\tau _{\mathrm{feed}}`$ is more consistent than $`\tau _{\mathrm{gas}}`$ or $`\tau _{\mathrm{dyn}}`$ with the short durations in our models for the recent starburst episodes in M 82. Admittedly, our $`\tau _{\mathrm{feed}}`$ estimates are also simplistic but they provide an additional perspective for the evolution of starburst activity and a possible explanation for short global burst durations. Along with our spatially detailed models, the above timescales are consistent with the overall progression of starburst activity in M 82 being driven by large-scale dynamical processes related to the tidal interaction with M 81 and the induced stellar bar, but suggests that once starburst activity was triggered, strong negative feedback effects acted very rapidly to inhibit further star formation on all spatial scales.

Recent studies of other starburst systems also provide evidence for short burst durations and for the episodic nature of starburst activity (e.g., Thornley et al., 2000; Böker et al., 1997; Vanzi et al., 1998; Alonso-Herrero et al., 2001). Some numerical simulations of the dynamical evolution of barred galaxies which account for the effects of star formation indicate a recurrent burst behaviour with typical timescale $`10^7\mathrm{yr}`$, owing to the ability of massive stars to destroy the flow pattern when a sufficient number coexist in a given place (Heller & Shlosman, 1994; Knapen et al., 1995).

## 7 SUMMARY

We have presented the results of detailed modeling of the central starburst regions of M 82, based on new near-IR integral field spectroscopy and mid-IR spectroscopy complemented with additional results taken from the literature. We applied starburst models optimized for observations in the range $`\lambda =145\mu \mathrm{m}`$ to these data and constrained the star formation parameters (cutoffs of the IMF, burst timescale) as well as the spatial and temporal evolution of starburst activity on scales as small as 25 pc.

Our data and models are consistent with the formation of very massive stars ($`50\mathrm{M}_{}`$) and imply burst durations of at most a few million years for individual burst sites. With typical sizes on 10–pc scales or less, individual burst sites in M 82 are comparable to Galactic and near-extragalactic massive star-forming regions such as NGC 3603 and the R 136 cluster in 30 Doradus. For plausible limits on the mass involved in the starburst, we find that the IMF must flatten at low masses assuming a Salpeter slope at high masses. Dynamical information for three individual young stellar clusters in M 82 (Smith & Gallagher, 2001; McCrady et al., 2003) provide similar evidence, but clearly more clusters need to be investigated. In view of the large uncertainties involved and the scarcity of low-mass IMF determinations in local high-mass star-forming regions, this result is difficult to interpret.

Among possible alternatives, our model results and interpretation together with the spatial distribution of the stellar and gaseous components in M 82 leads us to propose the following scenario. Globally, starburst activity in the central 500 pc of M 82 occurred during two successive episodes each lasting a few million years, about 10 and 5 Myr ago. The first episode took place throughout the disk and was particularly intense in the central few tens of parsecs. The second episode occurred predominantly in circumnuclear regions and along the stellar bar. We interpret this sequence in a tidally-triggered, bar-driven evolution scenario, consistent with the evidence for gravitational interaction between M 82 and its neighbour M 81, with the observed morphological signatures of the bar and associated dynamical resonances, and with the corresponding dynamical timescales.

The episodic nature of starburst activity and short global burst decay timescales in M 82 are particularly interesting results. Invoking a simple argument based on the comparison of the cumulative mechanical energy injected into the ISM and its gravitational binding energy, we find that the collective effect of massive stars winds and supernova explosions can disrupt the ISM in M 82 on timescales of $`10^6\mathrm{yr}`$, providing a very efficient quenching mechanism for each starburst episode. Following brief gas compression events on large-scales, starburst activity decayed very rapidly due to its own strong negative feedback effects. Our scenario for M 82 outlines the interplay between the large-scale triggering mechanisms and the more local but important feedback effects in determining the evolution of starburst activity in this galaxy.

We would like to thank A. Pauldrach and R.-P. Kudritzki for providing us with their model atmospheres, G. Ferland for making CLOUDY version C90.05 available in advance of publication, and C. Telesco and D. Gezari for their 12.4 µm data. It is a pleasure to thank M. Thornley, M. Lehnert, J. Gallimore, L. Tacconi, and L. Tacconi-Garman for stimulating discussions and useful comments, and to H. Spoon for help with using CLOUDY. We also wish to thank the referee, C. Leitherer, for a very considered and helpful report. This work was part of the Ph.D. thesis of NMFS, who acknowledges the Fonds pour les Chercheurs et l’Aide à la Recherche (Gouvernement du Québec, Canada) for a Graduate Scholarship, and the Max-Planck-Institut für extraterrestrische Physik and Service d’Astrophysique of the CEA Saclay for additional financial support. SWS and the ISO Spectrometer Data Center at MPE are supported by DLR under grants 50 QI 8610 8 and 50 QI 9402 3. We also thank the German-Israeli Foundation for support of this work (grant I-551-186.07/97).

## Appendix A FURTHER DETAILS ON STARS

### A.1 Spectral Energy Distributions

STARS computes the detailed SED intended for photoionization modeling of the nebular emission excited by the stars. The SEDs are taken from a hybrid grid generated as described by Thornley et al. (2000) from the LTE models of Kurucz (1992) for effective temperatures $`T_{\mathrm{eff}}19,000\mathrm{K}`$ and the non-LTE models of Pauldrach et al. (1998) for $`T_{\mathrm{eff}}25,000\mathrm{K}`$, with SEDs for intermediate temperatures obtained by interpolation. The Pauldrach et al. models represent much better the effects of line blocking and blanketing in the rapidly expanding atmospheres of hot stars and are generally harder than the Kurucz models. A detailed description of the models and implications on the computed SEDs is given by Pauldrach et al. (2001) and references therein. The basic Pauldrach et al. model SEDs we used in STARS include dwarfs and supergiants. For dwarfs, they cover $`2500060000\mathrm{K}`$ in $`T_{\mathrm{eff}}`$ at fixed surface gravity $`\mathrm{log}(g)=4`$ and mass-loss rates $`\dot{M}`$ of $`10^85\times 10^6\mathrm{M}_{}\mathrm{yr}^1`$. The ranges for supergiants are $`T_{\mathrm{eff}}=2500050000\mathrm{K}`$, $`\mathrm{log}(g)=2.753.8`$, and $`\dot{M}`$ of $`2\times 10^61.5\times 10^5\mathrm{M}_{}\mathrm{yr}^1`$. For a given $`T_{\mathrm{eff}}`$ and $`L_{\mathrm{bol}}`$ along a stellar track, STARS uses the corresponding evolutionary $`\mathrm{log}(g)`$ and assigns the SED model with closest $`T_{\mathrm{eff}}`$ and $`\mathrm{log}(g)`$ from the interpolated hybrid grid. The SED model grid covers well the range in parameters of the Geneva stellar tracks during most of the hot stars lifetimes (for both normal and enhanced mass-loss rates sets).

Various quantities predicted by STARS are derived using the empirical bolometric corrections of Schmidt-Kaler (1982) and the broad-band colours of Koornneef (1983). This is especially important for cool stars which have SEDs severely distorted by extensive molecular absorption and contaminated by thermal emission from a circumstellar dust shell that are not reproduced satisfactorily by black-body approximations or currently available model atmospheres. The impact of using empirical data is largest at near-IR wavelengths. For instance, the predicted $`K`$-band luminosities are up to 30% larger for the age range $`10100\mathrm{Myr}`$ (assuming $`t_{\mathrm{sc}}=1\mathrm{Myr}`$).

### A.2 Near-IR Stellar Absorption Features

STARS predicts the EWs of near-IR stellar absorption features characteristic of late-type stars: the $`{}_{}{}^{12}\mathrm{CO}`$ (6,3) and $`{}_{}{}^{12}\mathrm{CO}`$ (2,0) bandheads at 1.62 and 2.29 µm, and the Si I, Na I, and Ca I features at 1.59, 2.21, and 2.26 µm. We implemented in STARS the EWs versus $`T_{\mathrm{eff}}`$ relationships for luminosity classes V, III, and I below 5700 K obtained by polynomial fitting to the homogeneized stellar data compiled by Förster Schreiber (2000). The EWs assigned to hotter stars were identically zero since the corresponding features vanish at higher $`T_{\mathrm{eff}}`$. The exception is the 1.59 µm feature but above $`6000\mathrm{K}`$, the hydrogen Br14 transition is responsible for the absorption feature which has then a different physical meaning. The predicted EWs include “dilution” effects due to featureless continuum emission from OB stars and from associated free-free and free-bound processes. The stellar data are given for spectral resolutions of $`R1600`$ and 2500 in the $`H`$ and $`K`$ band, respectively. For comparison with the models, EWs measured from lower resolution data should be corrected according to the relationships given by Förster Schreiber (2000).

Recently, Origlia et al. (1999) and Origlia & Oliva (2000) questioned the reliability of synthesis modeling of the CO bandheads from a critical analysis of theoretical modeling and observational data. They emphasized the importance of accounting for AGB evolution up to the end of the thermally-pulsing phase. They also demonstrated that uncertainties in stellar evolutionary tracks for red supergiant and AGB phases remain large especially at sub-solar metallicities but affect much less the results for ages up to $`10^7\mathrm{yr}`$ at near-solar metallicities. Since STARS accounts for the thermally-pulsing AGB (see below) and with the solar-metallicity tracks adopted in this work, these uncertainties are minimized in our predictions of the CO EWs.

### A.3 Implementation of the TP-AGB Phase

The Geneva tracks (Schaller et al., 1992) employed in STARS follow the evolution of stars up to well defined evolutionary phases. Tracks for massive stars ($`M>7\mathrm{M}_{}`$) extend up to the end of the core carbon-burning phase, just prior to their explosion as core-collapse supernovae. Tracks for low-mass stars ($`M<2\mathrm{M}_{}`$) end at the helium flash which occurs at ages $`10^9\mathrm{yr}`$. The subsequent evolution of low-mass stars is not relevant for the starburst ages ($`<10^9\mathrm{yr}`$) we consider. For intermediate-mass stars ($`2\mathrm{M}_{}<M<7\mathrm{M}_{}`$), the Geneva tracks are terminated at the end of the early asymptotic giant branch (E-AGB), at ages in the range $`5\times 10^7`$ to $`10^9\mathrm{yr}`$, and thus exclude the TP-AGB phase. TP-AGB stars can contribute significantly to the near-IR and total luminosities in starbursts, and so should be included in synthesis models (e.g., Charlot & Bruzual, 1991; Maraston, 1998; Origlia & Oliva, 2000; Mouhcine & Lançon, 2002). Observationally, there is growing evidence that such intermediate-age populations account for a non-negligible fraction of the integrated light in the central regions or bulges of several nearby galaxies, including our own (e.g., Freedman, 1992; Blum, Sellgren, & DePoy, 1996; McLean & Liu, 1996; Lee, 1996; Davidge et al., 1997; Martínez-Delgado & Aparicio, 1998).

We extended the Geneva tracks for intermediate-mass stars using an approach suggested by Charlot & Bruzual (1991) which employs the analysis of the TP-AGB presented by Bedijn (1988). In this procedure, the TP-AGB is divided into a beginning “Mira phase” where the stars pulsate in the first overtone (or higher pulsation modes) while losing mass slowly, and a late “OH/IR” phase where the stars pulsate in the fundamental mode and lose mass rapidly. In the Mira phase, the stars increase in luminosity (and cool) and move up the AGB track appropriate for their mass. In the OH/IR phase, the luminosities do not vary much because of the rapid mass loss. Bedijn (1988) provides analysis and computations of the durations of each phase, including the maximum luminosities $`L_{\mathrm{max}}`$ at the end of the OH/IR phase, as functions of the main-sequence masses between 2 and 7$`\mathrm{M}_{}`$. Using this information, we added two evolutionary points in the $`L_{\mathrm{bol}}`$ versus $`T_{\mathrm{eff}}`$ Hertzsprung-Russell diagram. The first point corresponds to the end of the Mira phase where (following Bedijn) we assume that $`L_{\mathrm{bol}}=0.8L_{\mathrm{max}}`$. The second point corresponds to the end of the OH/IR phase where $`L_{\mathrm{bol}}=L_{\mathrm{max}}`$. We obtained the effective temperature at each point from the $`L_{\mathrm{bol}}T_{\mathrm{eff}}`$ relations for stars on the AGB. Our extensions are illustrated in Figure 15. As sanity check, we compared the $`L_{\mathrm{bol}}`$ predicted by our shortest burst models based on conventional synthesis and by applying the fuel consumption theorem (Renzini & Buzzoni, 1986) to our library of tracks. We obtain similar values within 10% on average up to $`10^9\mathrm{yr}`$, with maximum differences of $`\pm 25\%`$ at $`3\times 10^8\mathrm{yr}`$ increasing to $`\pm 40\%`$ at older ages where the numerical fluctuations due to the finite number of tracks become more important.

In attributing near-IR properties to TP-AGB stars, we followed a simple approach which can be justified by the small dispersion in $`T_{\mathrm{eff}}`$ of TP-AGB stars of various masses and the fact that the time spent in the OH/IR phase is at most 30% of the total TP-AGB lifetime according to Bedijn’s models (see also Charlot & Bruzual, 1991). We assigned to all TP-AGB stars the average bolometric correction and broad-band colours of the prototypical Mira variable $`\omega `$ Cet, derived from the data of Mendoza & Johnson (1965) and Zhou et al. (1984). For the EWs of near-IR absorption features, we adopted preliminary data obtained with the 3D instrument for a sample of Miras and other TP-AGB stars (S- and N-type stars). The observed CO bandhead strengths are comparable to those of the latest normal M-type giants. Similar conclusions are reached by inspection of the spectra published, e.g., by Lançon et al. (1999) and Lançon & Wood (2000). We therefore attributed the EWs of normal giants with $`T_{\mathrm{eff}}=3000\mathrm{K}`$ to all TP-AGB stars.

As expected, the effects of including the TP-AGB phase are most important for $`L_{\mathrm{bol}}`$ and for the near-IR properties while they are negligible for optical and UV properties. Figure 16 shows the TP-AGB contribution in our synthetized $`L_{\mathrm{bol}}`$, $`L_K`$, and $`W_{2.29}`$. Detailed comparison with other models including the TP-AGB, notably by Charlot & Bruzual (1991), Maraston (1998), and Mouhcine & Lançon (2002), is difficult because of the differences in synthesis techniques (conventional, isochrone, fuel consumption theorem) and in assumptions for TP-AGB stars (e.g., $`T_{\mathrm{eff}}`$ scale, bolometric corrections, detailed evolution, lifetimes) among the various works. The time evolution and difference between $`L_{\mathrm{bol}}`$ and $`L_K`$ of the TP-AGB contribution in our models agree qualitatively with these other models but the comparison suggests that the peak at $`5\times 10^8\mathrm{yr}`$ for $`t_{\mathrm{sc}}=1\mathrm{Myr}`$ may be too narrow and with too large an amplitude. Around $`10^8\mathrm{yr}`$, our short burst predictions lie in the ranges computed for coeval populations by the above authors, from a few to $`15\%`$ for $`L_{\mathrm{bol}}`$ and $`10\%`$ to $`40\%`$ for $`L_K`$ (with fairly large differences between different work). We point out that we implemented the TP-AGB to assess to first order the effects on the properties synthetized by STARS, and that our results of M 82 would be essentially unchanged without it. Accurate modeling would require a more careful and detailed treatment as presented by, e.g., Maraston (1998) and Mouhcine & Lançon (2002), as well as a better understanding of the properties and evolution of TP-AGB stars.