Obscuring Material around Seyfert Nuclei with Starbursts

The structure of obscuring matter in the environment of active galactic nuclei with associated nuclear starbursts is investigated using 3-D hydrodynamical simulations. Simple analytical estimates suggest that the obscuring matter with energy feedback from supernovae has a torus-like structure with a radius of several tens of parsecs and a scale height of about 10 pc. These estimates are confirmed by the fully non-linear numerical simulations, in which the multi-phase inhomogeneous interstellar matter and its interaction with the supernovae are consistently followed. The globally stable, torus-like structure is highly inhomogeneous and turbulent. To achieve the high column densities (> 10^{24} cm^{-2}) as suggested by observations of some Seyfert 2 galaxies with nuclear starbursts, the viewing angle should be larger than about 70 degree from the pole-on for a 10^8 solar mass massive black hole. Due to the inhomogeneous internal structure of the torus, the observed column density is sensitive to the line-of-sight, and it fluctuates by a factor of order 100. The covering fraction for N > 10^{23} cm^{-2} is about 0.4. The average accretion rate toward R < 1 pc is 0.4 solar mass/yr, which is boosted to twice that in the model without the energy feedback.Comment: ApJL in press (4 pages, 3 figures) A gziped ps file with high resolution figures is available at http://th.nao.ac.jp/~wada/AGN

Beware of fake AGNs

In the BPT diagram, the distribution of the emission-line galaxies from the Sloan Digital Sky Survey (SDSS) evokes the wings of a seagull. Traditionally, galaxies in the right wing are considered to host AGNs. Our study of the stellar populations of SDSS galaxies showed that about1/4 of galaxies thought to host LINERS are in fact "retired galaxies", i.e. galaxies that stopped forming stars and are ionized by hot post-AGB stars and white dwarfs (Stasinska et al. 2008). When including the galaxies that lack some of the lines needed to place them in the BPT diagram the fraction of retired galaxies is even larger (Cid Fernandes et al., 2009, arXiv:0912.1376)Comment: to be published in "Co-evolution of central black holes and galaxies: feeding and feed-back" Proceedings IAU Symposium No. 267, Peterson, Rachel Somerville, & Thaisa Storchi-Bergmann ed

On the semi-dynamical reflection equation: solutions and structure matrices

Explicit solutions of the non-constant semi-dynamical reflection equation are constructed, together with suitable parametrizations of their structure matrices. Considering the semi-dynamical reflection equation with rational non-constant Arutyunov-Chekhov-Frolov structure matrices, and a specific meromorphic ansatz, it is found that only two sets of the previously found constant solutions are extendible to the non-constant case. In order to simplify future constructions of spin-chain Hamiltonians, a parametrization procedure is applied explicitly to all elements of the semi-dynamical reflection equation available. Interesting expressions for `twists' and R-matrices entering the parametrization procedure are found. In particular, some expressions for the R-matrices seem to appear here for the first time. In addition, a new set of consistent structure matrices for the semi-dynamical reflection equation is obtained.Comment: typos corrected and some comments adde

BOND: Bayesian Oxygen and Nitrogen abundance Determinations in giant H II regions using strong and semi-strong lines

We present BOND, a Bayesian code to simultaneously derive oxygen and nitrogen abundances in giant H II regions. It compares observed emission lines to a grid of photoionization models without assuming any relation between O/H and N/O. Our grid spans a wide range in O/H, N/O and ionization parameter U, and covers different starburst ages and nebular geometries. Varying starburst ages accounts for variations in the ionizing radiation field hardness, which arise due to the ageing of H II regions or the stochastic sampling of the initial mass function. All previous approaches assume a strict relation between the ionizing field and metallicity. The other novelty is extracting information on the nebular physics from semi-strong emission lines. While strong lines ratios alone ([O III]/Hbeta, [O II]/Hbeta and [N II]/Hbeta) lead to multiple O/H solutions, the simultaneous use of [Ar III]/[Ne III] allows one to decide whether an H II region is of high or low metallicity. Adding He I/Hbeta pins down the hardness of the radiation field. We apply our method to H II regions and blue compact dwarf galaxies, and find that the resulting N/O vs O/H relation is as scattered as the one obtained from the temperature-based method. As in previous strong-line methods calibrated on photoionization models, the BOND O/H values are generally higher than temperature-based ones, which might indicate the presence of temperature fluctuations or kappa distributions in real nebulae, or a too soft ionizing radiation field in the models.Comment: MNRAS in press; 21 pages, 22 figures, 2 tables; code, data and results available at http://bond.ufsc.b

A comprehensive classification of galaxies in the SDSS: How to tell true from fake AGN?

We use the W_Ha versus [NII]/Ha (WHAN) diagram to provide a comprehensive emission-line classification of SDSS galaxies. This classification is able to cope with the large population of weak line galaxies that do not appear in traditional diagrams due to a lack of some of the diagnostic lines. A further advantage of the WHAN diagram is to allow the differentiation between two very distinct classes that overlap in the LINER region of traditional diagnostic diagrams. These are galaxies hosting a weakly active nucleus (wAGN) and "retired galaxies" (RGs), i.e. galaxies that have stopped forming stars and are ionized by their hot evolved low-mass stars. A useful criterion to distinguish true from fake AGN (i.e. the RGs) is the ratio (\xi) of the extinction-corrected L_Ha with respect to the Ha luminosity expected from photoionization by stellar populations older than 100 Myr. This ratio follows a markedly bimodal distribution, with a \xi >> 1 population composed by systems undergoing star-formation and/or nuclear activity, and a peak at \xi ~ 1 corresponding to the prediction of the RG model. We base our classification scheme on the equivalent width of Ha, an excellent observational proxy for \xi. Based on the bimodal distribution of W_Ha, we set the division between wAGN and RGs at W_Ha = 3 A. Five classes of galaxies are identified within the WHAN diagram: (a) Pure star forming galaxies: log [NII]/Ha 3 A. (b) Strong AGN (i.e., Seyferts): log [NII]/Ha > -0.4 and W_Ha > 6 A. (c) Weak AGN: log [NII]/Ha > -0.4 and W_Ha between 3 and 6 A. (d) RGs: W_Ha < 3 A. (e) Passive galaxies (actually, line-less galaxies): W_Ha and W_[NII] < 0.5 A. A comparative analysis of star formation histories and of other properties in these different classes of galaxies corroborates our proposed differentiation between RGs and weak AGN in the LINER-like family. (Abridged)Comment: Accepted for publication in MNRA