MaNGA Survey Research Articles

Photometric properties of classical bulge and pseudo-bulge galaxies at 0.5 ≤ z < 1.0

Context. We compare the photometric properties and specific star formation rate (sSFR) of classical- and pseudo-bulge galaxies with M∗ ≥ 109.5 M⊙ at 0.5 ≤ z < 1.0, selected from all five CANDELS fields. We also compare these properties of bulge galaxies at lower redshift selected from MaNGA survey in previous work. Aims. This paper aims to study the properties of galaxies with classical and pseudo-bulges at intermediate redshift, to compare the differences between different bulge types, and to understand the evolution of bulges with redshift. Methods. Galaxies are classified into classical bulge and pseudo-bulge samples according to the Sérsic index n of the bulge component based on results of two-component decomposition of galaxies, as well as the position of bulges on the Kormendy diagram. For the 105 classical bulge and 86 pseudo-bulge galaxies selected, we compare their size, luminosity, and sSFR of various components. Results. At a given stellar mass, most classical bulge galaxies have smaller effective radii, larger B/T, brighter and relatively larger bulges, and less active star formation than pseudo-bulge galaxies. In addition, the two types of galaxies have larger differences in sSFR at large radii than at the central region at both low- and mid-redshifts. Conclusions. The differences between the properties of the two types of bulge galaxies are generally smaller at mid-redshift than at low-redshift, indicating that they are evolving to more distinct populations towards the local universe. Bulge type is correlated with the properties of their outer disks, and the correlation is already present at redshifts as high as 0.5 < z < 1.

Nebular dust attenuation with the Balmer and Paschen lines based on the MaNGA survey

Dust attenuations observed by stars and ionized gas are not necessarily the same. The lack of observational constraints on the nebular dust attenuation curve leaves a large uncertainty when correcting nebular dust attenuation with stellar continuum-based attenuation curves. Making use of the DAP catalogs of the MaNGA survey, we investigate the nebular dust attenuation of regions traced by the Balmer and Paschen lines. Based on a simple simulation, we find that star-forming regions on kpc scales favor the classic foreground screen dust model rather than the uniform mixture model. We propose a novel approach to fit the dust attenuation curve using the emission-line fluxes directly. For strong hydrogen recombination lines (e.g. and the slopes of the nebular attenuation curve can be well determined and are found to be in good agreement with the Fitzpatrick Milky Way extinction curve with an accuracy of $ 4<!PCT!>$ in terms of the correction factor. However, severe contaminations and/or systematic uncertainties prevent us from obtaining reasonable values of the slopes for weak recombination lines (e.g., the high-order Balmer lines or the Paschen lines). We discuss how the choice of emission line measurement methods affects the results. Our results demonstrate the difficulty of deriving an average nebular dust attenuation curve given the current ground-based emission-line measurements.

Determining star formation rates in AGN hosts from strong optical emission lines

Abstract The influence of Active Galactic Nuclei (AGN) on star formation within their host galaxies remains a topic of intense debate. One of the primary challenges in quantifying the star formation rate (SFR) within AGN hosts arises from the prevalent assumption in most methodologies, which attribute gas excitation to young stars alone. However, this assumption does not consider the contribution of the AGN to the ionization of the gas in their environment. To address this issue, we evaluate the use of strong optical emission lines to obtain the SFR surface density (ΣSFRAGN) in regions predominantly ionized by an AGN, using a sample of 293 AGN hosts from the MaNGA survey, with SFR measurements available through stellar population fitting. We propose calibrations involving the Hα and [O iii]λ5007 emission lines, which can be used to determine ΣSFRAGN, resulting in values consistent with those estimated through stellar population fitting.

Environmental effects as a key factor in shaping star-forming S0 galaxies

The origins of lenticular galaxies (S0s) can be classified into two main categories: ``minor mergers" in low-density environments (LDEs) and ``faded spirals" in high-density environments (HDEs). The transitional phase in the evolution of S0s, namely, star-forming lenticular galaxies (SFS0s), can serve as an important probe for analyzing the complex processes involved in the transformation between different galaxy types and the quenching of star formation (SF). We attempt to find the impact of different environments on the global properties and spatially resolved quantities of SFS0s. We selected 71 SFS0s from the SDSS-IV MaNGA Survey, comprising 23 SFS0s in HDEs (SFS0s$\_$HE) and 48 SFS0s in LDEs (SFS0s$\_$LE). We examined the effects of the environment, by studying the global properties, concentration index, and radial profiles of the derived quantities. The varied environments of SFS0s do not lead to any significant difference in global properties (e.g., S$ e $rsic index). By calculating $CI_ alpha /cont $, we observe that different environments may cause varying concentrations of SF. Specifically, SFS0s$\_$LE, affected by external gas mergers or inflow, exhibit a more centrally concentrated SF (i.e., larger alpha /cont $). This trend is further supported by $CI_ SFR alpha $, which only considers the gas disk of the galaxy. This observation is aligned with the observed shrinking of gas disks in galaxies affected by ram-pressure stripping in HDEs. Furthermore, their $ SFR $ or resolved sSFR are comparable. On average, SFS0s$\_$LE display significantly higher values for both quantities. Finally, the observed n 4000$ and gas-phase metallicity gradient correspond well to their assumed origins. However, we did not find a significantly lower gas-phase metallicity in SFS0s$\_$LE. We suggest that different environments (i.e., origins) do not have a significant impact on the global properties of SFS0s, but they do indeed affect the distribution of SF. Considering the size of our sample and the unique nature of the galaxy, additional atomic and molecular gas data may provide further details to improve our understanding of these systems.

Bipolar Blobs as Evidence of Hidden AGN Activities in the Low-mass Galaxies

We report evidence of a hidden black hole (BH) in a low-mass galaxy, MaNGA 9885-9102, and provide a new method to identify active BHs in low-mass galaxies. This galaxy is originally selected from the MaNGA survey with distinctive bipolar Hα blobs at the minor axis. The bipolar feature can be associated with active galactic nuclei (AGN) activity, while the two blobs are classified as the H ii regions on the BPT diagram, making the origins confusing. The Swift UV continuum shows that the two blobs do not have UV counterparts, suggesting that the source of ionization is out of the blobs. Consistent with this, the detailed photoionization models prefer AGN rather than the star-forming origin with a significance of 5.8σ. The estimated BH mass is M BH ∼ 7.2 × 105 M ⊙ from the M BH–σ * relationship. This work introduces a novel method for detecting the light echo of BHs, potentially extending to intermediate mass, in low-metallicity environments where the traditional BPT diagram fails.

Finding Passive Galaxies in H i-MaNGA: The Impact of Star Formation Rate Indicator

H i-rich galaxies typically have high star formation rates (SFRs), but there exist interesting H i-rich and low star-forming galaxies. Previous work on a sample of these galaxies identified from H i-Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) (H i follow-up to the MaNGA survey) using an infrared indicator of specific-SFR (sSFR; namely W2 − W3 < 2) could find no single physical process to explain their unusual behavior. The method by which galaxies are identified as low sSFR may be important in this conclusion. In this Research Note, we explore how an H i-rich, low sSFR sample of H i-MaNGA galaxies differs using Hα, single stellar population, and ultraviolet estimators of SFR. We find that samples are statistically similar to each other so long as W2 − W3 < 2 is interpreted as corresponding to sSFR < 10−11.15 yr−1.

CO Observations of Early-mid Stage Major Mergers in the MaNGA Survey

We present a study of the molecular gas in early-mid stage major mergers, with a sample of 43 major-merger galaxy pairs selected from the Mapping Nearby Galaxies at Apache Point Observatory survey and a control sample of 195 isolated galaxies selected from the xCOLD GASS survey. Adopting kinematic asymmetry as a new effective indicator to describe the merger stage, we aim to study the role of molecular gas in the merger-induced star formation enhancement along the merger sequence of galaxy pairs. We obtain the molecular gas properties from CO observations with the James Clerk Maxwell Telescope, Institut de Radioastronomie Milimetrique 30 m telescope, and the MaNGA-ARO Survey of CO Targets survey. Using these data, we investigate the differences in molecular gas fraction ( fH2 ), star formation rate (SFR), star formation efficiency (SFE), molecular-to-atomic gas ratio ( MH2 /M H i ), total gas fraction (f gas), and the SFE of total gas (SFEgas) between the pair and control samples. In the full pair sample, our results suggest the fH2 of paired galaxies is significantly enhanced, while the SFE is comparable to that of isolated galaxies. We detect significantly increased fH2 and MH2 /M H i in paired galaxies at the pericenter stage, indicating an accelerated transition from atomic gas to molecular gas due to interactions. Our results indicate that the elevation of fH2 plays a major role in the enhancement of global SFR in paired galaxies at the pericenter stage, while the contribution of enhanced SFE in specific regions requires further explorations through spatially resolved observations of a larger sample spanning a wide range of merger stages.

Resolved properties of classical bulge and pseudo-bulge galaxies

ABSTRACT We compare properties of classical and pseudo-bulges and properties of their hosting galaxies selected from the MaNGA survey. Bulge types are identified based on the Sérsic index n of bulge component and the position of bulges on the Kormendy diagram. For the 393 classical bulges and 422 pseudo-bulges selected and their hosting galaxies, we study their kinematic properties including a proxy for specific angular momentum and central velocity dispersion, their stellar population properties including stellar age, metallicity, and specific star formation rate, as well as H i fractions of the galaxies. Our results show that at given stellar mass, disc components of pseudo-bulge galaxies are younger, have more active star formation, rotate more, and may contain more H i content compared with those of classical bulge galaxies, and the differences are larger than those between bulges themselves. The correlations between bulge types and disc properties indicate that different types of bulges are shaped by different processes that may regulate both growth of central components and evolution of outer discs in galaxies. In addition, we propose a stellar mass dependent divider of central velocity dispersion to separate galaxies with classical bulges from those with pseudo-bulges in galaxy mass range of $10.4\lt \mathrm{log}(M_*/\rm M_\odot)\lt 11.4$: $\mathrm{log}(\sigma _0) = 0.23 \times \mathrm{log}(M_*/\rm M_\odot)-0.46$. Galaxies with larger/smaller σ0 can be classified as hosts of classical/pseudo-bulges.

SDSS-IV MaNGA: calibration of astrophysical line-widths in the Hα region using HexPak observations

ABSTRACT We have re-observed $\rm \sim$40 low-inclination, star-forming galaxies from the MaNGA survey (σ ∼ 65 km s−1) at ∼6.5 times higher spectral resolution (σ ∼ 10 km s−1) using the HexPak integral field unit on the WIYN 3.5-m telescope. The aim of these observations is to calibrate MaNGA’s instrumental resolution and to characterize turbulence in the warm interstellar medium and ionized galactic outflows. Here we report the results for the Hα region observations as they pertain to the calibration of MaNGA’s spectral resolution. Remarkably, we find that the previously reported MaNGA line-spread-function (LSF) Gaussian width is systematically underestimated by only 1 per cent. The LSF increase modestly reduces the characteristic dispersion of H ii regions-dominated spectra sampled at 1–2 kpc spatial scales from 23 to 20 km s−1 in our sample, or a 25 per cent decrease in the random-motion kinetic energy. This commensurately lowers the dispersion zeropoint in the relation between line-width and star-formation rate surface-density in galaxies sampled on the same spatial scale. This modest zero-point shift does not appear to alter the power-law slope in the relation between line-width and star-formation rate surface-density. We also show that adopting a scheme whereby corrected line-widths are computed as the square root of the median of the difference in the squared measured line width and the squared LSF Gaussian avoids biases and allows for lower signal-to-noise data to be used reliably.

MaNGA AGN dwarf galaxies (MAD) – I. A new sample of AGNs in dwarf galaxies with spatially-resolved spectroscopy

ABSTRACT The finding of active galactic nuclei (AGNs) in dwarf galaxies has important implications for galaxy evolution and supermassive black hole formation models. Yet, how AGNs in dwarf galaxies form is still debated, in part due to scant demographics. We use the MaNGA survey, comprising ∼10 000 galaxies at z &amp;lt; 0.15, to identify AGN dwarf galaxies using a spaxel-by-spaxel classification in three spatially-resolved emission-line diagnostic diagrams (the [N ii]-, [S ii]-, and [O i]-BPT) and the WHAN diagram. This yields a sample of 664 AGN dwarf galaxies, the largest to date, and an AGN fraction of $\sim 20~{{\ \rm per\ cent}}$ that is significantly larger than that of single-fibre-spectroscopy studies (i.e. $\sim 1~{{\ \rm per\ cent}}$). This can be explained by the lower bolometric luminosity (&amp;lt;1042 erg s−1) and accretion rate (sub-Eddington) of the MaNGA AGN dwarf galaxies. We additionally identify 1176 SF-AGNs (classified as star-forming in the [N ii]-BPT but as AGNs in the [S ii]- and [O i]-BPT), 122 Composite, and 173 LINER sources. The offset between the optical centre of the galaxy and the median position of the AGN spaxels is more than 3 arcsec for ∼62% of the AGNs, suggesting that some could be off-nuclear. We also identify seven new broad-line AGNs with log MBH = 5.0–5.9 M⊙. Our results show how integral-field spectroscopy is a powerful tool for uncovering faint and low-accretion AGNs and better constraining the demographics of AGNs in dwarf galaxies.

Chemical evolution of local post-starburst galaxies: implications for the mass–metallicity relation

ABSTRACT We use the stellar fossil record to constrain the stellar metallicity evolution and star-formation histories of the post-starburst (PSB) regions within 45 local PSB galaxies from the MaNGA survey. The direct measurement of the regions’ stellar metallicity evolution is achieved by a new two-step metallicity model that allows for stellar metallicity to change at the peak of the starburst. We also employ a Gaussian process noise model that accounts for correlated errors introduced by the observational data reduction or inaccuracies in the models. We find that a majority of PSB regions (69 per cent at &amp;gt;1σ significance) increased in stellar metallicity during the recent starburst, with an average increase of 0.8 dex and a standard deviation of 0.4 dex. A much smaller fraction of PSBs are found to have remained constant (22 per cent) or declined in metallicity (9 per cent, average decrease 0.4 dex, standard deviation 0.3 dex). The pre-burst metallicities of the PSB galaxies are in good agreement with the mass–metallicity (MZ) relation of local star-forming galaxies. These results are consistent with hydrodynamic simulations, which suggest that mergers between gas-rich galaxies are the primary formation mechanism of local PSBs, and rapid metal recycling during the starburst outweighs the impact of dilution by any gas inflows. The final mass-weighted metallicities of the PSB galaxies are consistent with the MZ relation of local passive galaxies. Our results suggest that rapid quenching following a merger-driven starburst is entirely consistent with the observed gap between the stellar mass–metallicity relations of local star-forming and passive galaxies.

Multi-Epoch Optical Spectroscopy Variability of the Changing-Look AGN Mrk 883

In this work, we present multi-epoch optical spectra of the Seyfert 1.9 galaxy Mrk 883. Data were obtained with the Gran Telescopio Canarias and the MEGARA Integral Field Unit mode, archival data from the SDSS-IV MaNGA Survey and the SDSS-I Legacy Survey, and new spectroscopic observations obtained at San Pedro Mártir Observatory. We report the appearance of the broad component of Hβ emission line, showing a maximum FWHM ∼ 5927 ± 481 km s−1 in the MaNGA spectra, finding evidence for a change from Seyfert 1.9 (23 June 2003) to Seyfert 1.8 (18 May 2018). The observed changing-look variation from Sy1.9 to Sy1.8 has a timescale Δt ∼15 y. In addition, we observe profile and flux broad emission line variability from 2018 to 2023, and a wind component in [OIII]5007 Å, with a maximum FWHM = 1758 ± 178 km s−1, detected on 15 April 2023. In all epochs, variability of the broad lines was found to be disconnected from the optical continuum emission, which shows little or no variations. These results suggest that an ionized-driven wind in the polar direction could be a possible scenario to explain the observed changing-look variations.

Stellar mass is not the best predictor of galaxy metallicity

Context. Interpreting the scaling relations followed by galaxies is a fundamental tool for assessing how well we understand galaxy formation and evolution. Several scaling relations involving the galaxy metallicity have been discovered through the years, the foremost of which is the scaling with stellar mass. This so-called mass–metallicity relation is thought to be fundamental and has been subject to many studies in the literature. Aims. We study the dependence of the gas-phase metallicity on many different galaxy properties to assess which of them determines the metallicity of a galaxy. Methods. We applied a random forest regressor algorithm on a sample of more than 3000 nearby galaxies from the SDSS-IV MaNGA survey. Using this machine-learning technique, we explored the effect of 148 parameters on the global oxygen abundance as an indicator of the gas metallicity. Results. M⋆/Re, as a proxy for the baryonic gravitational potential of the galaxy, is found to be the primary factor determining the average gas-phase metallicity of the galaxy (Zg). It outweighs stellar mass. A subsequent analysis provides the strongest dependence of Zg on M⋆/Re0.6. We argue that this parameter traces the total gravitational potential, and the exponent α ≃ 0.6 accounts for the inclusion of the dark matter component. Conclusions. Our results reveal the importance of the relation between the total gravitational potential of the galaxy and the gas metallicity. This relation is tighter and likely more primordial than the widely known mass–metallicity relation.

SDSS-IV MaNGA: how do star formation histories affect gas-phase abundances?

ABSTRACT Gas-phase abundances in galaxies are the products of those galaxies’ evolutionary histories. The star formation history (SFH) of a region might therefore be expected to influence that region’s present day gaseous abundances. Here, we employ data from the MaNGA survey to explore how local gas metallicities relate to SFHs of galaxy regions. We combine MaNGA emission line measurements with SFH classifications from absorption line spectra to compare gas-phase abundances in star-forming regions with those in regions classified as starburst, post-starburst, and green valley. We find that starburst regions contain gas that is more pristine than in normal star-forming regions, in terms of O/H and N/O; we further find that post-starburst regions (which have experienced stochastic SFHs) behave very similarly to ordinary star-forming regions (which have experienced far smoother SFHs) in O/H–N/O space. We argue from this that gas is diluted significantly by pristine infall but is then re-enriched rapidly after a starburst event, making gas-phase abundances insensitive to the precise form of the SFH at late times. We also find that green valley regions possess slightly elevated N/O abundances at a given O/H; this is potentially due to a reduced star formation efficiency in such regions, but it could also point to late-time rejuvenation of green valley regions in our sample.

SDSS-IV MaNGA: the role of the environment in AGN triggering

ABSTRACT The large- and small-scale environments around optically-selected AGN host galaxies and a control sample of non-active galaxies in the MaNGA survey have been investigated in order to evaluate the importance of the environment in AGN triggering. Using the MaNGA integral field spectroscopy, we quantify non-circular motions of the ionized gas and detect an excess of radial gas motions in AGN hosts relative to control galaxies, not associated to AGN feedback and are most likely the result of tidal interactions, possibly associated with the triggering of the AGN. We find that the large-scale environments are similar for the AGN hosts and control galaxies in our sample and are biased towards lower large-scale densities and group virial masses, suggestive that the large-scale environment properties is only relevant to the AGN phenomenon in an indirect way, in the form, e.g. of the morphology-density relation. The small-scale environment, as measured by the frequency and luminosity of close neighbours, was also found to be similar for AGN and control galaxies. However, we find a correlation between the intensity of the non-circular gas motions in AGN hosts and the strength of the tidal field, while the control sample does not present such correlation. Also, AGN hosts with the most intense radial gas motions present larger tidal fields than their control galaxies. These findings indicate that at least a fraction of the AGN hosts in our sample have been triggered by tidal interactions with nearby galaxies.

Gas-phase Metallicity of Local Active Galactic Nuclei in the GASP and MaNGA Surveys: The Role of Ram Pressure Stripping

Growing evidence in support of a connection between active galactic nuclei (AGN) activity and the ram pressure stripping (RPS) phenomenon has been found both observationally and theoretically in the past decades. In this work, we further explore the impact of RPS on the AGN activity by estimating the gas-phase metallicity of nuclear regions and the mass–metallicity relation of galaxies at z ≤ 0.07 and with stellar masses logM*/M⊙≥9.0 , either experiencing RPS or not. To measure oxygen abundances, we exploit Integral Field Spectroscopy data from the GASP and MaNGA surveys, photoionization models generated with the code Cloudy and the code Nebulabayes to compare models and observations. In particular, we build Cloudy models to reproduce line ratios induced by photoionization from stars, AGN, or a contribution of both. We find that the distributions of metallicity and [O iii] λ5007 luminosity of galaxies undergoing RPS are similar to the ones of undisturbed galaxies. Independently of the RPS, we do not find a correlation between stellar mass and AGN metallicity in the mass range logM*/M⊙≥10.4 , while for the star-forming galaxies we observe the well-known mass–metallicity relation between 9.0≤logM*/M⊙≤10.8 with a scatter mainly driven by the star formation rate and a plateau around logM*/M⊙∼10.5 . The gas-phase metallicity in the nuclei of AGN hosts is enhanced with respect to those of star-forming galaxies by a factor of ∼ 0.05 dex regardless of the RPS.

Exploring Low Mass Galaxies at Different Wavelengths

We create a sample of potential AGN hosting dwarf galaxies from the MaNGA survey and AllWISE database for a redshift up to z 0.15. We select our sample with optical data obtained from the MaNGA survey then match these galaxies to sources in the AllWISE database, through a 10 arcsecond diameter search based on the right ascension and declination of our targets, to obtain infrared data. We find 141 likely potential AGN hosting galaxies from optical BPT analysis and 17 potential hosts from IR colour-colour analysis. We have compiled a list of galaxies for further analysis in other AGN indicating wavelengths, like radio and X-ray.

Where have all the low-metallicity galaxies gone? Tracing evolution in the mass–metallicity plane since a redshift of 0.7

ABSTRACT Even over relatively recent epochs, galaxies have evolved significantly in their location in the mass–metallicity plane, which must be telling us something about the latter stages of galaxy evolution. In this paper, we analyse data from the LEGA-C survey using semi-analytic spectral and photometric fitting to determine these galaxies’ evolution up to their observed epoch at z ∼ 0.7. We confirm that, at z ∼ 0.7, many objects already lie on the present-day mass–metallicity relation, but with a significant tail of high-mass low-metallicity galaxies that is not seen in the nearby Universe. Similar modelling of the evolution of galaxies in the nearby MaNGA survey allows us to reconstruct their properties at z ∼ 0.7. Once selection criteria similar to those of LEGA-C are applied, we reassuringly find that the MaNGA galaxies populate the mass–metallicity plane in the same way at z ∼ 0.7. Matching the LEGA-C sample to their mass–metallicity ‘twins’ in MaNGA at this redshift, we can explore the likely subsequent evolution of individual LEGA-C galaxies. Galaxies already on the present-day mass–metallicity relation form few more stars and their discs fade, so they become smaller and more bulge-like. By contrast, the high-mass low-metallicity galaxies grow their discs through late star formation, and evolve rapidly to higher metallicities due to a cut-off in their wind-driven mass-loss. There are significant indications that this late cut-off is associated with the belated end of strong active galactic nucleus activity in these objects.

MaNGA DynPop – IV. Stacked total density profile of galaxy groups and clusters from combining dynamical models of integral-field stellar kinematics and galaxy–galaxy lensing

ABSTRACT We present the measurement of total and stellar/dark matter decomposed mass density profile around a sample of galaxy groups and clusters with dynamical masses derived from integral-field stellar kinematics from the MaNGA survey in Paper I and weak lensing derived from the DECaLS imaging survey. Combining the two data sets enables accurate measurement of the radial density distribution from several kpc to Mpc scales. Intriguingly, we find that the excess surface density derived from stellar kinematics in the inner region cannot be explained by simply adding an NFW dark matter halo extrapolated from lensing measurement at a larger scale to a stellar mass component derived from the NASA-Sloan Atlas (NSA) catalogue. We find that a good fit to both data sets requires a stellar mass normalization about three times higher than that derived from the NSA catalogue, which would require an unrealistically too-heavy initial mass function for stellar mass estimation. If we keep the stellar mass normalization to that of the NSA catalogue but allow a varying inner dark matter density profile, we obtain an asymptotic slope of γgnfw = $1.82_{-0.25}^{+0.15}$ and γgnfw = $1.48_{-0.41}^{+0.20}$ for the group bin and the cluster bin, respectively, significantly steeper than the NFW case. We also compare the total mass inner density slopes with those from TNG300 and find that the values from the simulation are lower than the observation by about 2σ level.

A resolved study of the inner regions of nearby galaxies with an excess of young massive stars: missing link in the AGN–starburst connection?

ABSTRACT We have selected galaxies with very high levels of H α emission (EQW(H α) &amp;gt;700 Å.) in their central regions from the final data release of the MaNGA survey. Our study focuses on 14 very well-resolved nearby galaxies with stellar masses in the range 9.5 &amp;lt; log M*/(M⊙) &amp;lt; 11.5. We investigate a variety of procedures for selecting galaxy regions that are likely to harbour excess populations of young massive stars, finding that selection in the 2D space of extinction-corrected H α EQW and [S iii]/[S ii] line ratio produces the best results. By comparing stacked spectra covering these regions with stacked spectra covering normal starburst regions with 100 Å&amp;lt;EQW(H α) &amp;lt;200 Å, we obtain the following main results: (1) Clear signatures of excess Wolf–Rayet stars are found in half of the H α excess regions, (2) galaxy regions containing excess Wolf–Rayet stars are more often associated with the presence of high-ionization emission lines characteristic of accreting black holes. Excess [Ne iii] is detected in 4 out of 8 of the WR regions and there are tentative [Fe x] detections in 2 galaxies. (3) Regions of the galaxy with excess Wolf–Rayet stars are located where the interstellar medium has larger ionized gas turbulent velocities and higher neutral gas overdensities. We make a first attempt to constrain changes in the high-mass end of the stellar initial mass function (IMF) using the HR-pyPopStar evolutionary population synthesis models that include high-wavelength-resolution theoretical atmosphere libraries for Wolf–Rayet stars.