Nebular Lines Research Articles

V1674 Hercules: It is Blowing out a Wind

V1674 Her was not only the fastest nova on record, it was also a rare hybrid event, showing a transition between the Fe ii and He/N nova classes as early as 5.5 days after the outburst. Spectra now indicate the outburst is over; emission lines from the underlying binary system are now seen, although broad but weaker forbidden nebular emission lines are also present. Intriguing at this late evolutionary stage is a P Cygni-like line profile associated with Hα suggesting the presence of a wind emanating from the binary system. The profile exhibits intensity and velocity variations perhaps correlated with its 0.153 day orbital period.

The thesan project: predictions for multitracer line intensity mapping in the epoch of reionization

ABSTRACT Line intensity mapping (LIM) is rapidly emerging as a powerful technique to study galaxy formation and cosmology in the high-redshift Universe. We present LIM estimates of select spectral lines originating from the interstellar medium (ISM) of galaxies and 21 cm emission from neutral hydrogen gas in the Universe using the large volume, high resolution thesan reionization simulations. A combination of subresolution photoionization modelling for H ii regions and Monte Carlo radiative transfer calculations is employed to estimate the dust-attenuated spectral energy distributions (SEDs) of high-redshift galaxies (z ≳ 5.5). We show that the derived photometric properties such as the ultraviolet (UV) luminosity function and the UV continuum slopes match observationally inferred values, demonstrating the accuracy of the SED modelling. We provide fits to the luminosity–star formation rate relation (L–SFR) for the brightest emission lines and find that important differences exist between the derived scaling relations and the widely used low-z ones because the ISM of reionization era galaxies is generally less metal enriched than in their low-redshift counterparts. We use these relations to construct line intensity maps of nebular emission lines and cross-correlate with the 21 cm emission. Interestingly, the wavenumber at which the correlation switches sign (ktransition) depends heavily on the reionization model and to a lesser extent on the targeted emission line, which is consistent with the picture that ktransition probes the typical sizes of ionized regions. The derived scaling relations and intensity maps represent a timely state-of-the-art framework for forecasting and interpreting results from current and upcoming LIM experiments.

Linking stellar populations to H II regions across nearby galaxies

Context. Stellar feedback is one of the fundamental factors regulating the evolution of galaxies. However, we still do not have access to strong observational constraints on the relative importance of the different feedback mechanisms (e.g. radiation, ionised gas pressure, stellar winds) in driving H II region evolution and molecular cloud disruption. To quantify and compare the different feedback mechanisms, the size of an H II region is crucial, whereas samples of well-resolved H II regions are scarce. Aims. We constrain the relative importance of the various feedback mechanisms from young massive star populations by resolving H II regions across the disk of the nearby star-forming galaxy NGC 1672. Methods. We combined measurements of ionised gas nebular lines obtained by PHANGS-MUSE, with high-resolution (PSF FWHM ∼ 0.1″; ∼10 pc) imaging from Hubble Space Telescope (HST) in both the narrow-band Hα and broad-band (NUV, U, B, V, I) filters. We identified a sample of 40 isolated, compact H II regions in the HST Hα image. We measured the sizes of these H II regions, which were previously unresolved in seeing-limited ground-based observations. In addition, we identified the ionisation source(s) for each H II region from catalogues produced as part of the PHANGS-HST survey. In doing so, we were able to link young stellar populations with the properties of their surrounding H II regions. Results. The HST observations allowed us to resolve all 40 regions, with radii between 5 and 40 pc. The H II regions investigated here are mildly dominated by thermal or wind pressure, yet their elevation above the radiation pressure is within the expected uncertainty range. We see that radiation pressure provides a substantially higher contribution to the total pressure than previously found in the literature over similar size scales. In general, we find higher pressures within more compact H II regions, which is driven by the inherent size scaling relations of each pressure term, albeit with significant scatter introduced by the variation in the stellar population properties (e.g. luminosity, mass, age, metallicity). Conclusions. For nearby galaxies, the combination of MUSE/VLT observations with stellar population and resolved Hα observations from HST provides a promising approach that could yield the statistics required to map out how the importance of different stellar feedback mechanisms evolve over the lifetime of a H II region.

Diagnostics for PopIII galaxies and direct collapse black holes in the early universe

ABSTRACT Forthcoming observational facilities will make the exploration of the early universe routine, likely probing large populations of galaxies at very low metallicities. It will therefore be important to have diagnostics that can solidly identify and distinguish different classes of objects in such low metallicity regimes. We use new photoionization models to develop diagnostic diagrams involving various nebular lines. We show that combinations of these diagrams allow the identification and discrimination of the following classes of objects in the early universe: PopIII and direct collapse black holes (DCBH) in pristine environments, PopIII and DCBH embedded in slightly enriched interstellar medium (ISM; $\rm Z\sim 10^{-5}-10^{-4}$), and (metal poor) PopII and active galactic nucleus in enriched ISM. Diagnostics involving rest-frame optical lines (that will be accessible by James Webb Space Telescope) have a better discriminatory power but also rest-frame ultraviolet diagnostics can provide very useful information. Interestingly, we find that metal lines such as [O iii]λ5007 and C ivλ1549 can remain relatively strong (about a factor of 0.1–1 relative H β and He iiλ1640, respectively), even in extremely metal poor environments ($\rm Z\sim 10^{-5}-10^{-4}$), which could be embedding PopIII galaxies and DCBH.

Linking the properties of accreting white dwarfs with the ionization state of their ambient medium

ABSTRACT Steadily accreting white dwarfs (WDs) are efficient sources of ionization and thus are able to create extended ionized nebulae in their vicinity. These nebulae represent ideal tools for the detection of accreting WDs, given that in most cases the source itself is faint. In this work, we combine radiation transfer simulations with known H- and He-accreting WD models, providing for the first time the ionization state and the emission-line spectra of the formed nebulae as a function of the WD mass, the accretion rate and the chemical composition of the accreted material. We find that the nebular optical line fluxes and radial extent vary strongly with the WD’s accretion properties, peaking in systems with WD masses of 0.8–1.2 $\rm M_{\odot }$. Projecting our results on so-called BPT diagnostic diagrams, we show that accreting WD nebulae possess characteristics distinct from those of H ii-like regions, while they have line ratios similar to those in galactic low-ionization emission-line regions. Finally, we compare our results with the relevant constraints imposed by the lack of ionized nebulae in the vicinity of supersoft X-ray sources (SSSs) and Type Ia supernova remnants – sources that are related to steadily accreting WDs. The large discrepancies uncovered by our comparison rule out any steadily accreting WD as a potential progenitor of the studied remnants and additionally require the ambient medium around the SSSs to be less dense than 0.2 $\rm cm^{-3}$. We discuss possible alternatives that could bridge the incompatibility between the theoretical expectations and relevant observations.

Host galaxies of ultrastrong Mg ii absorbers at z ∼ 0.5

ABSTRACT From a sample of 109 candidate ultrastrong Mg ii (USMg ii; having rest equivalent width of Mg ii, W2796 > 3.0 Å) systems at z = 0.4–0.6, we confirm 27 and identify host galaxies of 20 systems based on associated nebular line emission from our SALT observations or from Sloan Digital Sky Survey (SDSS) fiber spectra. The measured impact parameter, [O ii] luminosity, star formation rate, B-band luminosity, and stellar mass are in the ranges 7.3 ≤ D[kpc] ≤ 79, $0.2\le L_{[\mathrm{ O}\,~\small {\rm II}]}[ 10^{41}~\mathrm{ erg} \mathrm{ s}^{-1}]\le 4.5$, 2.59 ≤ SFR[M⊙yr−1] ≤ 33.51, $0.15L_B^{*}\le L_B\le 1.63L_B^{*}$, and 10.21 ≤ log[M*/ M ⊙] ≤ 11.62, respectively. The impact parameters found are larger than that predicted by the W2796 versus D relationship of the general population of Mg ii absorbers. At a given D, USMg ii host galaxies are more luminous and massive compared to typical Mg ii absorbers. However, the measured SFRs are slightly lower than that of main-sequence galaxies with same M⋆ at z ∼ 0.5. We report a correlation between $L_{[\mathrm{ O}\,\small {\rm II}]}$ and W2796 for the full population of Mg ii absorbers, driven mainly by the host galaxies of weak Mg ii absorbers that tend to have low $L_{[\mathrm{ O}\,\small {\rm II}]}$ and large impact parameters. We find at least ∼33 per cent of the USMg ii host galaxies (with a limiting magnitude of mr < 23.6) are isolated and the large W2796 in these cases may originate from gas flows (infall/outflow) in single haloes of massive but not starburst galaxies. We also find galaxy interactions could be responsible for large velocity widths in at least ∼17 per cent cases.

Hαemission in local galaxies: star formation, time variability, and the diffuse ionized gas

ABSTRACTThe nebular recombination line H α is widely used as a star formation rate (SFR) indicator in the local and high-redshift Universe. We present a detailed H α radiative transfer study of high-resolution isolated Milky-Way and Large Magellanic Cloud simulations that include radiative transfer, non-equilibrium thermochemistry, and dust evolution. We focus on the spatial morphology and temporal variability of the H α emission, and its connection to the underlying gas and star formation properties. The H α and H β radial and vertical surface brightness profiles are in excellent agreement with observations of nearby galaxies. We find that the fraction of H α emission from collisional excitation amounts to fcol ∼ 5–$10{{\ \rm per\ cent}}$, only weakly dependent on radius and vertical height, and that scattering boosts the H α luminosity by $\sim 40{{\ \rm per\ cent}}$. The dust correction via the Balmer decrement works well (intrinsic H α emission recoverable within 25 per cent), though the dust attenuation law depends on the amount of attenuation itself both on spatially resolved and integrated scales. Important for the understanding of the H α–SFR connection is the dust and helium absorption of ionizing radiation (Lyman continuum [LyC] photons), which are about $f_{\rm abs}\approx 28{{\ \rm per\ cent}}$ and $f_{\rm He}\approx 9{{\ \rm per\ cent}}$, respectively. Together with an escape fraction of $f_{\rm esc}\approx 6{{\ \rm per\ cent}}$, this reduces the available budget for hydrogen line emission by nearly half ($f_{\rm H}\approx 57{{\ \rm per\ cent}}$). We discuss the impact of the diffuse ionized gas, showing – among other things – that the extraplanar H α emission is powered by LyC photons escaping the disc. Future applications of this framework to cosmological (zoom-in) simulations will assist in the interpretation of spectroscopy of high-redshift galaxies with the upcoming James Webb Space Telescope.

Chemistry of nebulae around binary post-AGB stars: A molecular survey of mm-wave lines

Context. There is a class of binary post-asymptotic giant branch (post-AGB) stars that exhibit remarkable near-infrared excess. Such stars are surrounded by Keplerian or quasi-Keplerian disks, as well as extended outflows composed of gas escaping from the disk. This class can be subdivided into disk- and outflow-dominated sources, depending on whether it is the disk or the outflow that represents most of the nebular mass, respectively. The chemistry of this type of source has been practically unknown thus far. Aims. Our objective is to study the molecular content of nebulae around binary post-AGB stars that show disks with Keplerian dynamics, including molecular line intensities, chemistry, and abundances. Methods. We focused our observations on the 1.3, 2, 3 mm bands of the 30mIRAM telescope and on the 7 and 13 mm bands of the 40 m Yebes telescope. Our observations add up ~600 h of telescope time. We investigated the integrated intensities of pairs of molecular transitions for CO, other molecular species, and IRAS fluxes at 12, 25, and 60 μm. Additionally, we studied isotopic ratios, in particular 17O/18O, to analyze the initial stellar mass, as well as 12CO/13CO, to study the line and abundance ratios. Results. We present the first single-dish molecular survey of mm-wave lines in nebulae around binary post-AGB stars. We conclude that the molecular content is relatively low in nebulae around binary post-AGB stars, as their molecular lines and abundances are especially weaker compared with AGB stars. This fact is very significant in those sources where the Keplerian disk is the dominant component of the nebula. The study of their chemistry allows us to classify nebulae around AC Her, the Red Rectangle, AI CMi, R Sct, and IRAS 20056+1834 as O-rich, while that of 89 Her is probably C-rich. The calculated abundances of the detected species other than CO are particularly low compared with AGB stars. The initial stellar mass derived from the 17O/18O ratio for the Red Rectangle and 89 Her is compatible with the central total stellar mass derived from previous mm-wave interferometric maps. The very low 12CO/13CO ratios found in binary post-AGB stars reveal a high 13CO abundance compared to AGB and other post-AGB stars.

Physical Properties of the Host Galaxies of Ca-rich Transients

Abstract Calcium-rich (Ca-rich) transients are a new class of supernovae (SNe) that are known for their comparatively rapid evolution, modest peak luminosities, and strong nebular calcium emission lines. Currently, the progenitor systems of Ca-rich transients remain unknown. Although they exhibit spectroscopic properties not unlike core-collapse Type Ib/c SNe, nearly half are found in the outskirts of their host galaxies, which are predominantly elliptical, suggesting a closer connection to the older stellar populations of SNe Ia. In this paper, we present a compilation of publicly available multiwavelength observations of all known and/or suspected host galaxies of Ca-rich transients ranging from far-UV to IR, and use these data to characterize their stellar populations with prospector. We estimate several galaxy parameters including integrated star formation rate, stellar mass, metallicity, and age. For nine host galaxies, the observations are sensitive enough to obtain nonparametric star formation histories, from which we recover SN rates and estimate probabilities that the Ca-rich transients in each of these host galaxies originated from a core-collapse versus Type Ia-like explosion. Our work supports the notion that the population of Ca-rich transients do not come exclusively from core-collapse explosions, and must either be only from white dwarf stars or a mixed population of white dwarf stars with other channels, potentially including massive star explosions. Additional photometry and explosion site spectroscopy of larger samples of Ca-rich host galaxies will improve these estimates and better constrain the ratio of white dwarf versus massive star progenitors of Ca-rich transients.

Follow-Up of Extended Shells around B[e] Stars

B[e] stars are massive B type emission line stars in different evolutionary stages ranging from pre-main sequence to post-main sequence. Due to their mass loss and ejection events these objects deposit huge amounts of mass and energy into their environment and enrich it with chemically processed material, contributing significantly to the chemical and dynamical evolution of their host galaxies. However, the large-scale environments of these enigmatic objects have not attracted much attention. The first and so far only catalog reporting the detection of extended shells around a sample of B[e] stars was an Hα imaging survey carried out in the year 2001, and was limited to bright targets in the northern hemisphere. We have recently started a follow-up of those targets to detect possible evolution of their nebulae in the plane of the sky over a baseline of two decades. Furthermore, we extend our survey to southern targets and fainter northern ones to complement and complete our knowledge on large-scale ejecta surrounding B[e] stars. Besides imaging in Hα and selected nebular lines, we utilize long-slit and 3D spectral observations across the nebulae to derive their physical properties. We discovered pronounced nebula structures around 15 more objects, resulting in a total of 27 B[e] stars with a large-scale nebula. Here we present our (preliminary) results for three selected objects: the two massive supergiants MWC137 and MWC 314, and the unclassified B[e] star MWC 819.

A tale of two DIGs: The relative role of H II regions and low-mass hot evolved stars in powering the diffuse ionised gas (DIG) in PHANGS–MUSE galaxies

We use integral field spectroscopy from the PHANGS–MUSE survey, which resolves the ionised interstellar medium structure at ∼50 pc resolution in 19 nearby spiral galaxies, to study the origin of the diffuse ionised gas (DIG). We examine the physical conditions of the diffuse gas by first removing morphologically defined H IIregions and then binning the low-surface-brightness areas to achieve significant detections of the key nebular lines in the DIG. A simple model for the leakage and propagation of ionising radiation from H IIregions is able to reproduce the observed distribution of Hαin the DIG. This model infers a typical mean free path for the ionising radiation of 1.9 kpc for photons propagating within the disc plane. Leaking radiation from H IIregions also explains the observed decrease in line ratios of low-ionisation species ([S II]/Hα, [N II]/Hα, and [O I]/Hα) with increasing Hαsurface brightness (ΣHα). Emission from hot low-mass evolved stars, however, is required to explain: (1) the enhanced low-ionisation line ratios observed in the central regions of some of the galaxies in our sample; (2) the observed trends of a flat or decreasing [O III]/Hβwith ΣHα; and (3) the offset of some DIG regions from the typical locus of H IIregions in the Baldwin–Phillips–Terlevich (BPT) diagram, extending into the area of low-ionisation (nuclear) emission-line regions (LI[N]ERs). Hot low-mass evolved stars make a small contribution to the energy budget of the DIG (2% of the galaxy-integrated Hαemission), but their harder spectra make them fundamental contributors to [O III] emission. The DIG might result from a superposition of two components, an energetically dominant contribution from young stars and a more diffuse background of harder ionising photons from old stars. This unified framework bridges observations of the Milky Way DIG with LI(N)ER-like emission observed in nearby galaxy bulges.

The VANDELS survey: Global properties of CIII]λ1908 Å emitting star-forming galaxies at z ∼ 3

Context. Strong nebular emission is ubiquitous in galaxies that contribute to cosmic reionization at redshift z ≳ 6. High-ionization UV metal lines, such as CIII]λ1908 Å, show high equivalent widths (EW) in these early galaxies, suggesting harder radiation fields at low metallicity than low-z galaxies of similar stellar mass. Understanding the physical properties driving the observed UV nebular line emission at high-z requires large and very deep spectroscopic surveys, which are now only accessible out to z ∼ 4. Aims. We study the mean properties of a large representative sample of 217 galaxies showing CIII] emission at 2 < z < 4, selected from a parent sample of ∼750 main-sequence star-forming galaxies in the VANDELS survey. These CIII] emitters have a broad range of UV luminosities, allowing for a detailed stacking analysis to characterize their stellar mass, star formation rate (SFR), and metallicity as a function of the UV emission line ratios, EWs, and the carbon-to-oxygen (C/O) abundance ratio. Methods. Stacking provides unprecedented high signal-to-noise (S/N) spectra for CIII] emitters over more than three decades in luminosity, stellar mass, and SFR. This enables a full spectral fitting to derive stellar metallicities for each stack. Moreover, we use diagnostics based on photoionization models and UV line ratios to constrain the ionization sources of the galaxies and derive the C/O abundance. Results. Reliable CIII] detections (S/N ≥ 3) represent ∼30% of the parent sample. However, stacked spectra of non-detections (S/N < 3) show weak (EW ≲ 2 Å) CIII] emission, suggesting that this line is common in normal star-forming galaxies at z ∼ 3. On the other hand, extreme CIII] emitters (EW(CIII]) ≳ 8 Å) are exceedingly rare (∼3%) in VANDELS. The UV line ratios of the sample suggest no ionization source other than massive stars. Stacks with larger EW(CIII]) show larger EW(Lyα) and lower metallicity, but not all CIII] emitters are Lyα emitters. The stellar metallicities of CIII] emitters are not significantly different from that of the parent sample, increasing from ∼10% to ∼40% solar for stellar masses log(M⋆/M⊙) ∼ 9−10.5. The stellar mass-metallicity relation of the CIII] emitters is consistent with previous works, exhibiting a strong evolution from z = 0 to z ∼ 3. The C/O abundances of the sample range between 35%−150% solar, with a noticeable increase with FUV luminosity and a smooth decrease with the CIII] EW. Here, we discuss the CIII] emitters in the C/O–Fe/H and the C/O–O/H planes and we find that they follow stellar and nebular abundance trends consistent with those of Milky Way halo and thick-disk stars and local HII galaxies, respectively. A qualitative agreement is also found with chemical evolution models, which suggests that CIII] emitters at z ∼ 3 are experiencing an active phase of chemical enrichment. Conclusions. Our results provide new insights into the nature of UV line emitters at z ∼ 2 − 4, paving the way for future studies at higher z using the James Webb Space Telescope.

The BPT Diagram in Cosmological Galaxy Formation Simulations: Understanding the Physics Driving Offsets at High Redshift

Abstract The Baldwin, Philips, & Terlevich diagram of [O iii]/Hβ versus [N ii]/Hα (hereafter N2-BPT) has long been used as a tool for classifying galaxies based on the dominant source of ionizing radiation. Recent observations have demonstrated that galaxies at z ∼ 2 reside offset from local galaxies in the N2-BPT space. In this paper, we conduct a series of controlled numerical experiments to understand the potential physical processes driving this offset. We model nebular line emission in a large sample of galaxies, taken from the simba cosmological hydrodynamic galaxy formation simulation, using the cloudy photoionization code to compute the nebular line luminosities from H ii regions. We find that the observed shift toward higher [O iii]/Hβ and [N ii]/Hα values at high redshift arises from sample selection: when we consider only the most massive galaxies M * ∼ 1010–11 M ⊙, the offset naturally appears, due to their high metallicities. We predict that deeper observations that probe lower-mass galaxies will reveal galaxies that lie on a locus comparable to z ∼ 0 observations. Even when accounting for samples-selection effects, we find that there is a subtle mismatch between simulations and observations. To resolve this discrepancy, we investigate the impact of varying ionization parameters, H ii region densities, gas-phase abundance patterns, and increasing radiation field hardness on N2-BPT diagrams. We find that either decreasing the ionization parameter or increasing the N/O ratio of galaxies at fixed O/H can move galaxies along a self-similar arc in N2-BPT space that is occupied by high-redshift galaxies.

The Mass–Metallicity Relation at Cosmic Noon in Overdense Environments: First Results from the MAMMOTH–Grism HST Slitless Spectroscopic Survey

Abstract The MAMMOTH–Grism slitless spectroscopic survey is a Hubble Space Telescope (HST) cycle 28 medium program, which is obtaining 45 orbits of WFC3/IR grism spectroscopy in the density peak regions of three massive galaxy protoclusters at z = 2–3 discovered using the MAMMOTH technique. We introduce this survey by presenting the first measurement of the mass–metallicity relation (MZR) at high redshift in overdense environments via grism spectroscopy. From the completed MAMMOTH–Grism observations in the field of the BOSS1244 protocluster at z = 2.24 ± 0.02, we secure a sample of 36 protocluster member galaxies at z ≈ 2.24, showing strong nebular emission lines ([O III], Hβ, and [O II]) in their G141 spectra. Using the multi-wavelength broadband deep imaging from HST and ground-based telescopes, we measure their stellar masses in the range of [109, 1010.4] M ⊙, instantaneous star formation rates (SFR) from 10 to 240 M ⊙ yr−1, and global gas-phase metallicities [ 1 3 , 1 ] of solar. Compared with similarly selected field-galaxy samples at the same redshift, our galaxies show, on average, increased SFRs by ∼0.06 dex and ∼0.18 dex at ∼1010.1 M ⊙ and ∼109.8 M ⊙, respectively. Using the stacked spectra of our sample galaxies, we derive the MZR in the BOSS1244 protocluster core as 12 + log ( O / H ) = 0.136 ± 0.018 × log ( M * / M ⊙ ) + 7.082 ± 0.175 , showing a significantly shallower slope than that in the field. This shallow MZR slope is likely caused by the combined effects of efficient recycling of feedback-driven winds and cold-mode gas accretion in protocluster environments. The former effect helps low-mass galaxies residing in overdensities retain their metal production, whereas the latter effect dilutes the metal content of high-mass galaxies, making them more metal-poor than their coeval field counterparts.

Abundance discrepancy factors in high-density planetary nebulae

ABSTRACT From high-resolution spectra, chemical abundances from collisionally excited lines (CELs) and optical recombination lines (ORLs) have been determined for planetary nebulae (PNe) Cn 3-1, Vy 2-2, Hu 2-1, Vy 1-2 and IC 4997, which are young and dense objects. The main aim of this work is to derive their O+2/H+ abundance discrepancy factors (ADFs) between CELs and ORLs. He, O, N, Ne, Ar, S, and Cl abundances were obtained and our values are in agreement with those previously reported. We found that Cn 3-1, Hu 2-1, and Vy 1-2 have O abundances typical of disc PNe, while Vy 2-2 and IC 4997 are low O abundance objects ($\rm {12+log(O/H) \sim 8.2}$), which can be attributed to possible O depletion into dust grains. ADFs(O+2) of $4.30^{+1.00}_{-1.16}$, 1.85 ± 1.05, $5.34^{+1.27}_{-1.08}$ and $4.87^{+4.34}_{-2.71}$ were determined for Vy 2-2, Hu 2-1, Vy 1-2, and IC 4997, respectively. The kinematics of CELs and ORLs was analysed for each case to study the possibility that different coexisting plasmas in the nebula emit them. Expansion velocities of [O iii] and O ii are equal within uncertainties in three PNe, providing no evidence for these lines being emitted in different zones. Exceptions are Hu 2-1 and Vy 2-2, where ORLs might be emitted in different zones than CELs. For Vy 2-2 and IC 4997, we found that nebular and auroral lines of the same ion (S+, N+, Ar+2, Ar+3, O+2) might present different expansion velocities. Auroral lines show lower $\rm {V_{exp}}$, which might indicate that they are emitted in a denser and inner zone than the nebular ones.

On the most luminous planetary nebulae of M 31

Context.The planetary nebula luminosity function (PNLF) is a standard candle that comprises a key rung on the extragalactic distance ladder. The method is based on the empirical evidence that the luminosity function of planetary nebulae (PNe) in the [O III]λ5007 nebular emission line reaches a maximum value that is approximately invariant with population age, metallicity, or host galaxy type. However, the presence of bright PNe in old stellar populations is not easily explained by single-star evolutionary models.Aims.To gain information about the progenitors of PNe at the tip of the PNLF, we obtained the deepest existing spectra of a sample of PNe in the galaxy M 31 to determine their physico-chemical properties and infer the post-asymptotic giant branch (AGB) masses of their central stars (CSs). Precise chemical abundances allow us to confront the theoretical yields for AGB stellar masses and metallicities expected at the bright end of the PNLF. Central star masses of the sampled PNe provide direct information on the controversial origin of the universal cutoff of the PNLF.Methods.Using the OSIRIS instrument at the 10.4 m Gran Telescopio Canarias (GTC), optical spectra of nine bright M 31 PNe were obtained: four of them at the tip of the PNLF, and the other five some 0.5 mag fainter. A control sample of 21 PNe with previous GTC spectra from the literature is also included. We analyze their physical properties and chemical abundances (He, N, O, Ar, Ne, and S), searching for relevant differences between bright PNe and the control samples. The CS masses are estimated with Cloudy modeling using the most recent evolutionary tracks.Results.The studied PNe show a remarkable uniformity in all their nebular properties, and the brightest PNe show relatively large electron densities. Stellar characteristics also span a narrow range: ⟨L*/L⊙⟩ = 4300 ± 310, ⟨Teff⟩ = 122 000 ± 10 600 K for the CSs of the four brightest PNe, and ⟨L*/L⊙⟩ = 3300 ± 370, ⟨Teff⟩ = 135 000 ± 26 000 K for those in the control set. This groups all the brightest PNe at the location of maximum temperature in the post-AGB tracks for stars with initial massesMi = 1.5M⊙.Conclusions.These figures provide robust observational constraints for the stellar progenitors that produce the PNLF cutoff in a star-forming galaxy such as M 31, where a large range of initial masses is in principle available. Inconsistency is found, however, in the computed N/O abundance ratios of five nebulae, which are 1.5 to 3 times larger than predicted by the existing nucleosynthesis models for stars of these masses.

Discovery of a luminous starburst galaxy with hundreds of thousands of Wolf–Rayet stars

ABSTRACT This paper reports the detection of a strong Wolf–Rayet (W-R) emission bump feature at 4400–4700 Å in the optical spectrum of a distant galaxy SDSS J150009.81+452844.4 at a redshift of 0.453. The estimated number of W-R stars is 1.5 × 105 without applying any internal extinction correction, and at least 4.5 × 105 after correcting for dust extinction. Such a number, though with uncertainties inherent from the extinction correction and others, appears to outnumber those in previous W-R galaxies after correcting for intrinsic dust extinction with the numbers available in the literature. These massive stars must have formed in an instantaneous star-forming episode lasting less than about 5 Myr. We estimate a star formation rate of at least ∼80 $\mathrm{M}_{\odot }\, \mathrm{yr^{-1}}$, indicative of a violent starburst. Its mid-infrared spectrum resembles closely those of typical local ultraluminous infrared galaxies (ULIRGs). From optical to ultraviolet, it possesses a luminous blue continuum, indicating that the starburst is not heavily obscured, in contrast to that found in most ULIRGs. There is evidence for strong outflows based on the detection of systematically blue-shifted broad wings in the nebular emission lines (a bulk velocity −190 km s−1), as well as in the blue-shifted Mg ii absorption lines. The emission-line gases show a wide velocity range, from −1200 to ∼−2000 km s−1 in blueshift to 570 km s−1 in redshift. We interpret this as a galactic fountain, of which part of the outflowing gas may be falling back to the galactic plane. Our method could be used to find a sample of similar objects, which would help understand the star formation history and stellar feedback in starburst galaxies.

Physical Properties of Massive Compact Starburst Galaxies with Extreme Outflows

Abstract We present results on the nature of extreme ejective feedback episodes and the physical conditions of a population of massive (M * ∼ 1011 M ⊙), compact starburst galaxies at z = 0.4–0.7. We use data from Keck/NIRSPEC, SDSS, Gemini/GMOS, MMT, and Magellan/MagE to measure rest-frame optical and near-IR spectra of 14 starburst galaxies with extremely high star formation rate surface densities (mean ΣSFR ∼ 2000 M ⊙ yr−1 kpc−2) and powerful galactic outflows (maximum speeds v 98 ∼ 1000–3000 km s−1). Our unique data set includes an ensemble of both emission ([O ii] λλ3726,3729, Hβ, [O iii] λλ4959,5007, Hα, [N ii] λλ6549,6585, and [S ii] λλ6716,6731) and absorption (Mg ii λλ2796,2803, and Fe ii λ2586) lines that allow us to investigate the kinematics of the cool gas phase (T ∼ 104 K) in the outflows. Employing a suite of line ratio diagnostic diagrams, we find that the central starbursts are characterized by high electron densities (median n e ∼ 530 cm−3), and high metallicity (solar or supersolar). We show that the outflows are most likely driven by stellar feedback emerging from the extreme central starburst, rather than by an AGN. We also present multiple intriguing observational signatures suggesting that these galaxies may have substantial Lyman continuum (LyC) photon leakage, including weak [S ii] nebular emission lines. Our results imply that these galaxies may be captured in a short-lived phase of extreme star formation and feedback where much of their gas is violently blown out by powerful outflows that open up channels for LyC photons to escape.

The extension of the fundamental metallicity relation beyond the BPT star-forming sequence: Evidence for both gas accretion and starvation

Context. The fundamental metallicity relation (FMR) of galaxies is a 3D relation between the gas-phase metallicity, stellar mass, and star-formation rate (SFR). So far, it has only been studied for galaxies identified as star forming (SF) on the BPT diagrams (BPT-SF), but not for galaxies with LI(N)ER/AGN classification (BPT-non-SF), mainly because of the lack of diagnostics for estimating their gas-phase metallicities in the latter cases. Aims. We extend the FMR to BPT-non-SF galaxies. Methods. To this end, we exploit the recent nebular line empirical calibrations derived specifically for galaxies classified as non-SF in the BPT diagrams. Moreover, we study an alternative representation of the FMR where we consider the offsets in metallicity and SFR with respect to main sequence (MS) galaxies. Results. We find that galaxies with SFR higher than the MS are more metal-poor than their counterparts on the MS, which is interpreted in terms of gas accretion boosting star formation and diluting the metallicity. Low-mass galaxies below the MS (i.e. towards quiescence) have higher metallicities than their MS counterparts, which is interpreted in terms of starvation (i.e. suppression of fresh gas supply) hampering star formation and reducing the dilution effect, hence resulting in a higher level of internal chemical enrichment. Massive galaxies below the MS have a gas metallicity that is much closer to their MS counterparts and much lower than expected from their stellar metallicities; this result suggests a scenario where massive nearly quiescent galaxies with LI(N)ER-like nebular emission have recently accreted gas from the circum/intergalactic medium.

Spectroscopy of an extreme [O iii] emitting active galactic nucleus at z = 3.212: implications for the reionization era

ABSTRACT Reionization-era galaxies often display intense nebular emission lines, both in rest-frame optical ([O iii] + H β) and ultraviolet (UV; C iii], C iv). How such strong nebular emission is powered remains unclear, with both active galactic nuclei (AGNs) and hot stars considered equally viable. The UV continuum slopes of these early systems tend to be very blue (β < −2), reflecting minimal dust obscuration, young ages, and low metallicities. This contrasts with narrow-lined AGNs at $z$ ∼ 2−3, whose UV slopes are significantly redder (β > −1) than typical star-forming systems in the reionization era. To investigate the properties of AGNs in the reionization era, we have conducted a search for potential examples of rare analogues with blue continua at intermediate redshift ($z$ ∼ 2−3). Our goals are to determine whether AGNs with intense line emission and blue continua exist and thereby to establish the range of rest-frame UV and optical line ratios in this population. In this paper, we report the detection of an X-ray luminous AGN at $z$ = 3.21 (UDS-24561) with extreme [O iii] + H β line emission (equivalent width = 1300 Å) and a blue UV continuum slope (β = −2.34). MMT/Binospec and Keck/MOSFIRE spectra indicate rest-frame UV line ratios consistent with AGN photoionization models and rest-frame optical lines with both a narrow component [full width at half-maximum (FWHM) $=154\, {\rm km\, s}^{-1}$] and extended broad wings (FWHM $=977\, {\rm km\, s}^{-1}$), consistent with outflowing gas. We describe how such objects can be identified in future James Webb Space Telescope emission line surveys in the reionization era, thereby providing a valuable census of AGN activity at $z$ > 6 and understanding their contribution to cosmic reionization.