Reconciling escape fractions and observed line emission in Lyman-continuum-leaking galaxies

Context. Finding and elucidating the properties of Lyman-continuum(LyC)-emitting galaxies is an important step in improving our understanding of cosmic reionization. Aims. Although the z ∼ 0.3 − 0.4 LyC emitters found recently show strong optical emission lines, no consistent quantitative photoionization model taking into account the escape of ionizing photons and inhomogenous interstellar medium (ISM) geometry of these galaxies has yet been constructed. Furthermore, it is unclear to what extent these emission lines can be used to distinguish LyC emitters. Methods. To address these questions we construct one- and two-zone photoionization models accounting for the observed LyC escape, which we compare to the observed emission line measurements. The main diagnostics used include lines of [O III], [O II], and [O I] plus sulfur lines ([S II], [S III]) and a nitrogen line ([N II]), which probe regions of different ionization in the ISM. Results. We find that single (one-zone) density-bounded photoionization models cannot reproduce the emission lines of the LyC leakers, as pointed out by earlier studies, because they systematically underpredict the lines of species of low ionization potential, such as [O I] and [S II]. Introducing a two-zone model, with differing ionization parameter and a variable covering fraction and where one of the zones is density-bounded, we show that the observed emission line ratios of the LyC emitters are well reproduced. Furthermore, our model yields LyC escape fractions, which are in fair agreement with the observations and independent measurements. The [O I] λ6300 excess, which is observed in some LyC leakers, can be naturally explained in this model, for example by emission from low-ionization and low-filling-factor gas. LyC emitters with a high escape fraction (fesc ≳ 38%) are deficient both in [O I] λ6300 and in [S II] λλ6716,6731. We also confirm that a [S II] λλ6716,6731 deficiency can be used to select LyC emitter candidates, as suggested earlier. Finally, we find indications for a possible dichotomy in terms of escape mechanisms for LyC photons between galaxies with relatively low (fesc ≲ 10%) and higher escape fractions. Conclusions. We conclude that two-zone photoionization models are sufficient and required to explain the observed emission line properties of z ∼ 0.3 − 0.4 LyC emitters. This is in agreement with UV absorption line studies, which also show the co-existence of regions with high hydrogen column density (i.e., no escape of ionizing photons) and density-bounded or very low column density regions responsible for the observed escape of LyC radiation. These simple but consistent models provide a first step towards the use of optical emission lines and their ratios as quantitative diagnostics of LyC escape from galaxies.

LIN 358: a symbiotic binary accreting above the steady hydrogen fusion limit

ABSTRACT Symbiotic binaries are long-period interacting binaries consisting of a white dwarf (WD) accreting material from a cool evolved giant star via stellar winds. In this paper, we study the symbiotic binary LIN 358 located in the Small Magellanic Cloud. We have observed LIN 358 with the integral field spectrograph WiFeS and obtained its line emission spectrum. With the help of the plasma simulation and spectral synthesis code cloudy, we have constructed a 2D photoionization model of LIN 358. From comparison with the observations, we have determined the colour temperature of the WD in LIN 358 to be 19 eV, its bolometric luminosity L = (1.02 ± 0.15) × 1038 erg s−1, and the mass-loss rate from the donor star to be 1.2 × 10−6 M⊙ yr−1. Assuming a solar H to He ratio in the wind material, a lower limit to the accreted mass fraction in LIN 358 is 0.31. The high mass accretion efficiency of a wind Roche lobe overflow implies that the WD is accreting above the upper boundary of stable hydrogen fusion and thus growing in mass with the maximal rate of ≈4 × 10−7 M⊙ yr−1. This causes the WD photosphere to expand, which explains its low colour temperature. Our calculations show that the circumstellar material in LIN 358 is nearly completely ionized except for a narrow cone around the donor star, and that the WD emission is freely escaping the system. However, due to its low colour temperature, this emission can be easily attenuated by even moderate amounts of neutral interstellar medium. We speculate that other symbiotic systems may be operating in a similar regime, thus explaining the paucity of observed systems.

Search for intrinsic NALs in BAL/mini-BAL quasar spectra

ABSTRACT Some fraction of narrow absorption lines (NALs) are physically associated to the quasar/host-galaxy materials (i.e. intrinsic NALs) like those of broad absorption lines (BALs) and mini-BALs. The relation between these three types of absorption lines has not been understood yet, however one interpretation is that these absorption features correspond to different inclination angles. In this study, we search for intrinsic NALs in 11 BAL/mini-BAL quasar spectra retrieved from VLT/UVES public archive, in order to test a possible relation of intrinsic NALs and BALs/mini-BALs in the geometry models. We use partial coverage analysis to separate intrinsic NALs from ones which are associated to cosmologically intervening materials like foreground galaxies and intergalactic medium (i.e. intervening NALs). We identify one reliable and two possible intrinsic NAL systems out of 36 NAL systems in 9 BAL/mini-BAL quasar spectra after removing two quasars without clear BAL features. In spite of a small sample size, we placed a lower limit on the fraction of BAL/mini-BAL quasars that have at least one intrinsic C iv NAL ($\sim 33^{+33}_{-18}{{\ \rm per\ cent}}$). This can be interpreted that intrinsic NAL absorbers exist everywhere regardless of inclination angle. We found that one of the intrinsic NAL systems detected in SDSS J121549.80−003432.1 is located at a large radial distance of R > 130 kpc, using a method of photoionization model with ground/excited-state lines. Considering the wide range of intrinsic NAL absorber distribution in inclination angles and radial distances, it suggests that origins and geometry of them are more complicated than we expected.

Chemical abundances of Seyfert 2 AGNs – III. Reducing the oxygen abundance discrepancy

ABSTRACT We investigate the discrepancy between oxygen abundance estimations for narrow-line regions of active galactic nuclei (AGNs) type Seyfert 2 derived using direct estimations of the electron temperature (Te-method) and those derived using photoionization models. In view of this, observational emission-line ratios in the optical range ($3000 \: \lt \: \lambda (\mathring{\rm A}) \: \lt 7000$) of Seyfert 2 nuclei compiled from the literature were reproduced by detailed photoionization models built with the cloudy code. We find that the derived discrepancies are mainly due to the inappropriate use of the relations between temperatures of the low (t2) and high (t3) ionization gas zones derived for H ii regions in AGN chemical abundance studies. Using a photoionization model grid, we derived a new expression for t2 as a function of t3 valid for Seyfert 2 nuclei. The use of this new expression in the AGN estimation of the O/H abundances based on Te-method produces O/H abundances slightly lower (about 0.2 dex) than those derived from detailed photoionization models. We also find that the new formalism for the Te-method reduces by about 0.4 dex the O/H discrepancies between the abundances obtained from strong emission-line calibrations and those derived from direct estimations.

Multiphase outflows in post-starburst E+A galaxies - II. A direct connection between the neutral and ionized outflow phases

ABSTRACT Post-starburst E+A galaxies are systems that hosted a powerful starburst that was quenched abruptly. Simulations suggest that these systems provide the missing link between major merger ULIRGs and red and dead ellipticals, where AGN feedback is responsible for the expulsion or destruction of the molecular gas. However, many details remain unresolved and little is known about AGN-driven winds in this short-lived phase. We present spatially resolved IFU spectroscopy with MUSE/VLT of SDSS J124754.95-033738.6, a post-starburst E+A galaxy with a recent starburst that started 70 Myr ago and ended 30 Myr ago, with a peak SFR of $\sim 150\, \mathrm{M_{\odot }\,yr^{ -1}}$. We detect disturbed gas throughout the entire field of view, suggesting triggering by a minor merger. We detect fast-moving multiphased gas clouds, embedded in a double-cone face-on outflow, which are traced by ionized emission lines and neutral NaID emission and absorption lines. We find remarkable similarities between the kinematics, spatial extents, and line luminosities of the ionized and neutral gas phases, and propose a model in which they are part of the same outflowing clouds, which are exposed to both stellar and AGN radiation. Our photoionization model provides consistent ionized line ratios, NaID absorption optical depths and EWs, and dust reddening. Using the model, we estimate, for the first time, the neutral-to-ionized gas mass ratio (about 20), the sodium neutral fraction, and the size of the outflowing clouds. This is one of the best ever observed direct connections between the neutral and ionized outflow phases in AGN.

Chemical abundances of Seyfert 2 AGNs – II. N2 metallicity calibration based on SDSS

ABSTRACT We present a semi-empirical calibration between the metallicity (Z) of Seyfert 2 active galactic nuclei and the N2 = log([N ii]λ6584/H α) emission-line intensity ratio. This calibration was derived through the [O iii]λ5007/[O ii]λ3727 versus N2 diagram containing observational data and photoionization model results obtained with the cloudy code. The observational sample consists of 463 confirmed Seyfert 2 nuclei (redshift $z \: \lesssim 0.4$) taken from the Sloan Digital Sky Survey DR7 data set. The obtained Z–N2 relation is valid for the range $0.3 \: \lesssim \: (Z/{\rm Z}_{\odot }) \: \lesssim \: 2.0$ that corresponds to $-0.7 \: \lesssim \: ({\rm N}2) \: \lesssim \: 0.6$. The effects of varying the ionization parameter (U), electron density and the slope of the spectral energy distribution on the Z estimations are of the order of the uncertainty produced by the error measurements of N2. This result indicates the large reliability of our Z –N2 calibration. A relation between U and the [O iii]/[O ii] line ratio, almost independent of other nebular parameter, was obtained.

Understanding the Spatial Distributions of the Ionic/Atomic/Molecular/Dust Components in PNe

Planetary nebulae (PNe) are often recognized as the hallmark of compact H ii regions in the Universe. However, there exist dusty neutral regions extending beyond the central ionized region. We demonstrate that such dusty neutral regions (also known as photo-dissociation regions, or PDRs) around the central ionized region are significant parts of PNe in terms of energetics and mass. We do so by using our latest dusty photoionization model of NGC 6781 (of 13 parameters) based on one of the most comprehensive panchromatic data sets ever assembled for a PN encompassing from X-ray to radio (of 136 constraining data, including 19 flux densities, 78 line fluxes, and 37 band fluxes). We find that NGC 6781, evolved out of a 2.25–3.0 M ⊙ star located 460 pc away from us, possesses a massive concentration of neutral gas (molecular hydrogen) just beyond the central ionized region and that the amount of ionized gas in NGC 6781 is only 22% of the observationally accounted amount of matter in the circumstellar environment, which itself does not even account for the amount of mass presumably ejected by the central star during the last thermal pulse event according to the latest evolutionary models. This means that the observed nebula in this PN is only the tip of the iceberg.