We have modeled a large sample of infrared starburst galaxies using both the PEGASE v2.0 and STARBURST99 codes to generate the spectral energy distribution (SED) of the young star clusters. PEGASE utilizes the Padova group tracks, while STARBURST99 uses the Geneva group tracks, allowing comparison between the two. We used our MAPPINGS III code to compute photoionization models that include a self-consistent treatment of dust physics and chemical depletion. We use the standard optical diagnostic diagrams as indicators of the hardness of the EUV radiation Ðeld in these galaxies. These diagnostic diagrams are most sensitive to the spectral index of the ionizing radiation Ðeld in the 1È4 ryd region. We Ðnd that warm infrared starburst galaxies contain a relatively hard EUV Ðeld in this region. The PEGASE ionizing stellar continuum is harder in the 1È4 ryd range than that of STARBURST99. As the spectrum in this regime is dominated by emission from Wolf-Rayet (W-R) stars, this discrepancy is most likely due to the di †erences in stellar atmosphere models used for the W-R stars. The PEGASE models use the Clegg & Middlemass planetary nebula nuclei (PNN) atmosphere models for the W-R stars, whereas the STARBURST99 models use the Schmutz, Leitherer, & Gruenwald W-R atmosphere models. We believe that the Schmutz et al. atmospheres are more applicable to the starburst galaxies in our sample ; however, they do not produce the hard EUV Ðeld in the 1È4 ryd region required by our observations. The inclusion of continuum metal blanketing in the models may be one solution. Supernova remnant (SNR) shock modeling shows that the contribution by mechanical energy from SNRs to the photoionization models is >20%. The models presented here are used to derive a new theoretical classiÐcation scheme for starbursts and active galactic nucleus (AGN) galaxies based on the optical diagnostic diagrams.

No abstract

We present a new library of fully-radiative shock models calculated with the MAPPINGS III shock and photoionization code. The library consists of grids of models with shock velocities in the range v=100-1000 km/s and magnetic parameters B/sqrt(n) of 10^-4 - 10 muG cm^(3/2) for five different atomic abundance sets, and for a pre-shock density of 1.0 cm^(-3). Additionally, Solar abundance model grids have been calculated for densities of 0.01, 0.1, 10, 100, and 1000 cm^(-3) with the same range in v and B/sqrt(n). Each model includes components of both the radiative shock and its photoionized precursor, ionized by the EUV and soft X-ray radiation generated in the radiative gas. We present the details of the ionization structure, the column densities, and the luminosities of the shock and its precursor. Emission line ratio predictions are separately given for the shock and its precursor as well as for the composite shock+precursor structure to facilitate comparison with observations in cases where the shock and its precursor are not resolved. Emission line ratio grids for shock and shock+precursor are presented on standard line ratio diagnostic diagrams, and we compare these grids to observations of radio galaxies and a sample of AGN and star forming galaxies from the Sloan Digital Sky Survey. This library is available online, along with a suite of tools to enable the analysis of the shocks and the easy creation of emission line ratio diagnostic diagrams. These models represent a significant increase in parameter space coverage over previously available models, and therefore provide a unique tool in the diagnosis of emission by shocks.Comment: 39 pages, 34 figures, accepted for publication in ApJS, April 200

No abstract

We use the chemical evolution predictions of cosmological hydrodynamic simulations with our latest theoretical stellar population synthesis, photoionization and shock models to predict the strong line evolution of ensembles of galaxies from z = 3 to the present day. In this paper, we focus on the brightest optical emission-line ratios, [N II]/Hα and [O III]/Hβ. We use the optical diagnostic Baldwin-Phillips-Terlevich (BPT) diagram as a tool for investigating the spectral properties of ensembles of active galaxies. We use four redshift windows chosen to exploit new near-infrared multi-object spectrographs. We predict how the BPT diagram will appear in these four redshift windows given different sets of assumptions. We show that the position of star-forming galaxies on the BPT diagram traces the ISM conditions and radiation field in galaxies at a given redshift. Galaxies containing AGN form a mixing sequence with purely star-forming galaxies. This mixing sequence may change dramatically with cosmic time, due to the metallicity sensitivity of the optical emission-lines. Furthermore, the position of the mixing sequence may probe metallicity gradients in galaxies as a function of redshift, depending on the size of the AGN narrow line region. We apply our latest slow shock models for gas shocked by galactic-scale winds. We show that at high redshift, galactic wind shocks are clearly separated from AGN in line ratio space. Instead, shocks from galactic winds mimic high metallicity starburst galaxies. We discuss our models in the context of future large near-infrared spectroscopic surveys.