Understanding the evolution of the N/O ratio in the interstellar medium (ISM) of galaxies is essential if we are to complete our picture of the chemical evolution of galaxies at high redshift, since most observational calibrations of O/H implicitly depend upon the intrinsic N/O ratio. The observed N/O ratio, however, shows large scatter at low O/H, and is strongly dependent on galactic environment. We show that several heretofore unexplained features of the N/O distribution at low O/H can be explained by the N seen in metal-poor galaxies being mostly primary nitrogen that is returned to the ISM via pre-supernova winds from rapidly rotating massive stars ($M \gtrsim 10$ M$_\odot$, $v/v_{\rm crit} \gtrsim 0.4$). This mechanism naturally produces the observed N/O plateau at low O/H. We show that the large scatter in N/O at low O/H also arises naturally from variations in star-formation efficiency. By contrast, models in which the N and O come primarily from supernovae provide a very poor fit to the observed abundance distribution. We propose that the peculiar abundance patterns we observe at low O/H are a signature that dwarf galaxies retain little of their SN ejecta, leaving them with abundance patterns typical of winds.

The evolutionary paths taken by massive stars with $M \gtrsim 60 \, \mathrm{M}_\odot$ remain substantially uncertain. They begin their lives as main sequence (MS) O-stars. Depending on their masses, rotation rates, and metallicities, they can then encounter a wide range of evolutionary states with an equally broad set of possible surface compositions and spectral classifications. We present a new grid of calculations for the evolution of such stars that covers a broad range in mass, M/M$_\odot = 60$ to $150$, rotation rate, $v \, / \, v_{\rm crit} = 0$ to $0.6$, metallicity, $[\mathrm{Fe}/\mathrm{H}] = -4$ to $0$, and $\alpha$-element enhancement, $[\alpha/\mathrm{Fe}] = 0$ to $0.4$. We show that rotating stars undergo rotationally-induced dredge-up of nucleosynthetic products, mostly He and N, to their surfaces while still on the MS. Non-rotating metal-rich stars also reveal the products of nucleosynthesis on their surfaces because even modest amounts of mass loss expose their "fossil" convective cores: regions that are no longer convective, but which were part of the convective core at an early stage in the star's evolution. Thus surface enhancement of He and N is expected for rotating stars at all metallicities, and for non-rotating stars if they are relatively metal-rich. We calculate a stellar atmosphere for a representative model from our grid, properly accounting for He- and N-enhancement, and show that the resulting spectrum provides a good match to observed WNL stars, strongly suggesting that the physical mechanisms we have identified are the ultimate cause of the WNL phase.

We review the use of emission-lines for understanding galaxy evolution, focusing on excitation source, metallicity, ionization parameter, ISM pressure and electron density. We show that the UV, optical and infrared contain complementary diagnostics that can probe the conditions within different nebular ionization zones. In anticipation of upcoming telescope facilities, we provide new self-consistent emission-line diagnostic calibrations for complete spectral coverage from the UV to the infrared. These diagnostics can be used in concert to understand how fundamental galaxy properties have changed across cosmic time. We describe new 2D and 3D emission-line diagnostics to separate the contributions from star formation, AGN and shocks using integral field spectroscopy. We discuss the physics, benefits, and caveats of emission-line diagnostics, including the effect of theoretical model uncertainties, diffuse ionized gas, and sample selection bias. Accounting for complex density gradients and temperature profiles is critical for reliably estimating the fundamental properties of H ii regions and galaxies. Diffuse ionized gas can raise metallicity estimates, flatten metallicity gradients, and introduce scatter in ionization parameter measurements. We summarize with a discussion of the challenges and major opportunities for emission-line diagnostics in the coming years.

There is compelling evidence for a highly energetic Seyfert explosion (10^{56-57} erg) that occurred in the Galactic Centre a few million years ago. The clearest indications are the x-ray/gamma-ray "10 kpc bubbles" identified by the Rosat and Fermi satellites. In an earlier paper, we suggested another manifestation of this nuclear activity, i.e. elevated H-alpha emission along a section of the Magellanic Stream due to a burst (or flare) of ionizing radiation from Sgr A*. We now provide further evidence for a powerful flare event: UV absorption line ratios (in particular CIV/CII, SiIV/SiII) observed by the Hubble Space Telescope reveal that some Stream clouds towards both galactic poles are highly ionized by a source capable of producing ionization energies up to at least 50 eV. We show how these are clouds caught in a beam of bipolar, radiative "ionization cones" from a Seyfert nucleus associated with Sgr A*. In our model, the biconic axis is tilted by about 15 deg from the South Galactic Pole with an opening angle of roughly 60 deg. For the Stream at such large Galactic distances (D > 75 kpc), nuclear activity is a plausible explanation for all of the observed signatures: elevated H-alpha emission and H ionization fraction (X_e > 0.5), enhanced CIV/CII and SiIV/SiII ratios, and high CIV and SiIV column densities. Wind-driven "shock cones" are ruled out because the Fermi bubbles lose their momentum and energy to the Galactic corona long before reaching the Stream. The nuclear flare event must have had a radiative UV luminosity close to the Eddington limit (f_E ~ 0.1-1). Our time-dependent Seyfert flare models adequately explain the observations and indicate the Seyfert flare event took place T_o = 3.5 +/- 1 Myr ago. The timing estimates are consistent with the mechanical timescales needed to explain the x-ray/gamma-ray bubbles in leptonic jet/wind models (2-8 Myr).

We derive new self-consistent theoretical UV, optical, and IR diagnostics for the ISM pressure and electron density in the ionized nebulae of star-forming galaxies. Our UV diagnostics utilize the inter-combination, forbidden and resonance lines of silicon, carbon, aluminum, neon, and nitrogen. We also calibrate the optical and IR forbidden lines of oxygen, argon, nitrogen and sulfur. We show that line ratios used as ISM pressure diagnostics depend on the gas-phase metallicity with a residual dependence on the ionization parameter of the gas. In addition, the traditional electron density diagnostic [S II] {\lambda}6731/[S II] {\lambda}6717 is strongly dependent on the gas-phase metallicity. We show how different emission-line ratios are produced in different ionization zones in our theoretical nebulae. The [S II] and [O II] ratios are produced in different zones, and should not be used interchangeably to measure the electron density of the gas unless the electron temperature is known to be constant. We review the temperature and density distributions observed within H II regions and discuss the implications of these distributions on measuring the electron density of the gas. Many H II regions contain radial variations in density. We suggest that the ISM pressure is a more meaningful quantity to measure in H II regions or galaxies. Specific combinations of line ratios can cover the full range of ISM pressures (4 < log(P/k) < 9). As H II regions become resolved at increasingly high redshift through the next generation telescopes, we anticipate that these diagnostics will be important for understanding the conditions around the young, hot stars from the early universe to the present day.

In the optical spectra of galaxies, the separation of line emission from gas ionised by star formation and an AGN, or by star formation and shocks, are very well-understood problems. However, separating line emission between AGN and shocks has proven difficult. With the aid of a new three-dimensional diagnostic diagram, we show the simultaneous separation of line emission from star formation, shocks, and AGN in NGC 1068, and quantify the ratio of star formation, shocks, and AGN in each spaxel. The AGN, shock, and star formation luminosity distributions across the galaxy accurately align with X-ray, radio, and CO(3-2) observations, respectively. Comparisons with previous separation methods show that the shocked emission heavily mixes with the AGN emission. We also show that if the H$\alpha$ flux is to be used as a star formation rate indicator, separating line emission from as many sources as possible should be attempted to ensure accurate results.

We study and compare different examples of stellar evolutionary synthesis input parameters used to produce photoionisation model grids using the MAPPINGS V modelling code. The aim of this study is to (a) explore the systematic effects of various stellar evolutionary synthesis model parameters on the interpretation of emission lines in optical strong-line diagnostic diagrams, (b) characterise the combination of parameters able to reproduce the spread of local galaxies located in the star-forming region in the Sloan Digital Sky Survey, and (c) investigate the emission from extremely metal-poor galaxies using photoionisation models. We explore and compare the stellar input ionising spectrum (stellar population synthesis code [Starburst99, SLUG, BPASS], stellar evolutionary tracks, stellar atmospheres, star-formation history, sampling of the initial mass function) as well as parameters intrinsic to the H II region (metallicity, ionisation parameter, pressure, H II region boundedness). We also perform a comparison of the photoionisation codes MAPPINGS and CLOUDY. On the variations in the ionising spectrum model parameters, we find that the differences in strong emission-line ratios between varying models for a given input model parameter are small, on average ~0.1 dex. An average difference of ~0.1 dex in emission-line ratio is also found between models produced with MAPPINGS and CLOUDY. Large differences between the emission-line ratios are found when comparing intrinsic H II region parameters. We find that low-metallicity galaxies are better explained by a density-bounded H II region and higher pressures better encompass the spread of galaxies at high redshift.

We systematically measure the gas-phase metallicities and the mass-metallicity relation of a large sample of local active galaxies for the first time. Observed emission-line fluxes from the Sloan Digital Sky Survey (SDSS) are compared to a four-dimensional grid of photoionization models using the Bayesian parameter estimation code NebulaBayes. For the first time we take into account arbitrary mixing between HII region and narrow-line region (NLR) emission, and the models are also varied with metallicity, ionization parameter in the NLR, and the gas pressure. The active galactic nucleus (AGN) oxygen abundance is found to increase by $\Delta {\rm O/H} \sim 0.1$ dex as a function of host galaxy stellar mass over the range $10.1 < \log M_* / M_\odot < 11.3$. We also measure the metallicity and ionization parameter of 231000 star-forming galaxies for comparison with the sample of 7670 Seyfert 2 galaxies. A systematic offset in oxygen abundance of 0.09 dex is observed between the mass-metallicity relations of the star-forming and active galaxies. We investigate potential causes of the offset, including sample selection and the treatment in the models of diffuse ionized gas, pressure, and ionization parameter. We cannot identify the major cause(s), but suspect contributions due to deficiencies in modeling the ionizing spectra and the treatment of dust physics. Optical diagnostic diagrams are presented with the star-forming and Seyfert data colored by the inferred oxygen abundance, ionization parameter and gas pressure, clearly illustrating the trends in these quantities.

CHIANTI contains a large quantity of atomic data for the analysis of astrophysical spectra. Programs are available in IDL and Python to perform calculation of the expected emergent spectrum from these sources. The database includes atomic energy levels, wavelengths, radiative transition probabilities, rate coefficients for collisional excitation, ionization, and recombination, as well as data to calculate free-free, free-bound, and two-photon continuum emission. In Version 9, we improve the modelling of the satellite lines at X-ray wavelengths by explicitly including autoionization and dielectronic recombination processes in the calculation of level populations for select members of the lithium isoelectronic sequence and Fe XVIII-XXIII. In addition, existing datasets are updated, new ions added and new total recombination rates for several Fe ions are included. All data and IDL programs are freely available at http://www.chiantidatabase.org or through SolarSoft and the Python code ChiantiPy is also freely available at https://github.com/chianti-atomic/ChiantiPy.

Using integral field data we extract the optical spectra of shocked interstellar clouds in Kepler's supernova remnant located in the inner regions of our Galaxy, as well as in the Large Magellanic Cloud (LMC), the Small Magellanic Cloud (SMC), NGC6822 and IC 1613. Using self-consistent shock modelling, we make a new determination of the chemical composition of the interstellar medium (ISM) in N, O, Ne, S, Cl and Ar in these galaxies and obtain accurate estimates of the fraction of refractory grains destroyed in the shock. By comparing our derived abundances with those obtained in recent works using observations of B stars, F supergiant stars and HII regions, we provide a new calibration for abundance scaling in the range $7.9 \lesssim 12+\log {\mathrm {O/H}} \lesssim 9.1$.

We study the interactions of a relativistic jet with a dense turbulent gaseous disk of radius $\sim 2$ kpc. We have performed a suite of simulations with different mean density, jet power and orientation. Our results show that: A) the relativistic jet couples strongly with the gas in the inner kpc, creating a cavity and launching outflows. B) The high pressure bubble inflated by the jet and its back-flow compresses the disk at the outer edges, driving inflows within the disk. C) Jets inclined towards the disk couple more strongly with the disk and launch sub-relativistic, wide-angle outflows along the minor axis. D) Shocks driven directly by the jet and the jet-driven energy bubble raise the velocity dispersion throughout the disk by several times its initial value. E) Compression by the jet-driven shocks can enhance the star formation rate in the disk, especially in a ring-like geometry close to the axis. However, enhanced turbulent dispersion in the disk also leads to quenching of star formation. Whether positive or negative feedback dominates depends on jet power, ISM density, jet orientation with respect to the disc, and the time-scale under consideration. Qualitatively, our simulations compare favourably with kinematic and morphological signatures of several observed galaxies such as NGC 1052, NGC 3079, 3C 326 and 3C 293.

We fit theoretical models to the emission-line spectra of 2766 Seyfert galaxies from the Sloan Digital Sky Survey (SDSS). The model line fluxes are derived by `mixing' photoionization model predictions of active galactic nucleus (AGN) narrow line region (NLR) emission and HII region emission. The observed line fluxes in each spectrum were directly compared to the grid of mixed models using the Bayesian parameter estimation code NebulaBayes, thereby measuring the degree of mixing in each spectrum for the first time. We find that the majority of the Balmer line emission in the majority of Seyfert-classified SDSS spectra arises from contaminating HII regions within the fixed-size aperture. Even for spectra with log [OIII]/Hb $\gtrsim 0.9$, a median of ~30% of the Balmer flux arises in HII regions. We show that the results are qualitatively insensitive to the assumed Seyfert ionizing continuum, and that ionizing spectra with a peak energy of E_peak ~ 40-50 eV produce the most plausible distributions of mixing fractions. The analysis cleanly quantifies how the starburst - AGN `mixing fraction' increases on the BPT diagram for SDSS galaxies. Apart from the mixing fraction, the models also vary with the ionization parameter in the NLR, the gas pressure, and the metallicity. Measurements for the metallicity in particular will be presented in a companion paper.

From integral field data we extract the optical spectra of 20 shocked clouds in the supernova remnant N132D in the Large Magellanic Cloud (LMC). Using self-consistent shock modelling, we derive the shock velocity, pre-shock cloud density and shock ram pressure in these clouds. We show that the [Fe X] and [Fe XIV] emission arises in faster, partially radiative shocks moving through the lower density gas near the periphery of the clouds. In these shocks dust has been effectively destroyed, while in the slower cloud shocks the dust destruction is incomplete until the recombination zone of the shock has been reached. These dense interstellar clouds provide a sampling of the general interstellar medium (ISM) of the LMC. Our shock analysis allows us to make a new determination of the ISM chemical composition in N, O, Ne, S, Cl and Ar, and to obtain accurate estimates of the fraction of refractory grains destroyed. From the derived cloud shock parameters, we estimate cloud masses and show that the clouds previously existed as typical self-gravitating Bonnor-Ebert spheres into which converging cloud shocks are now being driven.

NebulaBayes is a new Bayesian code that implements a general method of comparing observed emission-line fluxes to photoionization model grids. The code enables us to extract robust, spatially resolved measurements of abundances in the extended narrow line regions (ENLRs) produced by Active Galactic Nuclei (AGN). We observe near-constant ionization parameters but steeply radially-declining pressures, which together imply that radiation pressure regulates the ENLR density structure on large scales. Our sample includes four `pure Seyfert' galaxies from the S7 survey that have extensive ENLRs. NGC2992 shows steep metallicity gradients from the nucleus into the ionization cones. An {\it inverse} metallicity gradient is observed in ESO138-G01, which we attribute to a recent gas inflow or minor merger. A uniformly high metallicity and hard ionizing continuum are inferred across the ENLR of Mrk573. Our analysis of IC5063 is likely affected by contamination from shock excitation, which appears to soften the inferred ionizing spectrum. The peak of the ionizing continuum E_peak is determined by the nuclear spectrum and the absorbing column between the nucleus and the ionized nebula. We cannot separate variation in this intrinsic E_peak from the effects of shock or HII region contamination, but E_peak measurements nevertheless give insights into ENLR excitation. We demonstrate the general applicability of NebulaBayes by analyzing a nuclear spectrum from the non-active galaxy NGC4691 using a HII region grid. The NLR and HII region model grids are provided with NebulaBayes for use by the astronomical community.

The interstellar medium of the radio galaxy IC 5063 is highly perturbed by an AGN jet expanding in the gaseous disc of the galaxy. We model this interaction with relativistic hydrodynamic simulations and multiphase initial conditions for the interstellar medium and compare the results with recent observations. As the jets flood through the inter-cloud channels of the disc, they ablate, accelerate, and disperse clouds to velocities exceeding $400 \mbox{km s}^{-1}$. Clouds are also destroyed or displaced in bulk from the central regions of the galaxy. Our models with jet powers of $10^{44} \mbox{erg s}^{-1}$ and $10^{45} \mbox{erg s}^{-1}$ are capable of reproducing many of the features seen in the position-velocity diagram published in Morganti et al. (2015) and confirm the notion that the jet is responsible for the strongly perturbed gas dynamics seen in the ionized, neutral, and molecular gas phases. In our simulations, we also see strong venting of the jet plasma perpendicular to the disc, which entrains clumps and diffuse filaments into the halo of the galaxy. Our simulations are the first 3D hydrodynamic simulations of the jet and ISM of IC 5063.

We propose that Gigahertz Peak Spectrum (GPS) and Compact Steep Spectrum (CSS) radio sources are the signposts of relativistic jet feedback in evolving galaxies. Our simulations of relativistic jets interacting with a warm, inhomogeneous medium, utilize cloud densities and velocity dispersions in the range derived from optical observations, show that free-free absorption can account for the $\sim \rm GHz$ peak frequencies and low frequency power laws inferred from the radio observations. These new computational models replace the power-law model for the free-free optical depth in the \citep{bicknell97a} model by a more fundamental model involving disrupted log-normal distributions of warm gas. One feature of our new models is that at early stages, the low frequency spectrum is steep but progressively flattens as a result of a broader distribution of optical depths, suggesting that the steep low frequency spectra discovered by \citet{callingham17a} may possibly be attributed to young sources. We also investigate the inverse correlation between peak frequency and size and find that the initial location on this correlation is determined by the average density of the warm ISM. The simulated sources track this correlation initially but eventually fall below it, indicating the need for a more extended ISM than presently modelled. GPS and CSS sources can potentially provide new insights into the phenomenon of AGN feedback since their peak frequencies and spectra are indicative of the density, turbulent structure and distribution of gas in the host galaxy.

In this paper we present integral field spectroscopy of the extraordinary Type I bipolar planetary nebula Hen 2-111. In the lobes we map fast moving knots of material with [N II]$\lambda 6584$/H$\alpha$ ratios up to 12, and with radial velocities relative to systemic from -340 km/s up to +390 km/s. We find evidence of a bipolar ejection event at a velocity $\sim 600$ km/s from the central star (assumed to be a binary), which occurred about 8000yr ago. The fast moving material is chemically quite distinct from the lower velocity gas in the bipolar lobes., and displays very high N abundances. We show that the fast moving N-rich knots are not photoionised by the central star, and have constructed detailed shock models for the brightest knot. We find a pre-shock density $\sim 6$cm$^{-3}$, and a shock velocity $\sim150$ km/s. The shock is not fully radiative, being only $\sim 600$yr old. This shocked gas is partially H-burnt, with a helium abundance by mass exceeding that of hydrogen, and is interacting with partially H-burnt material ejected in an earlier episode of mass loss. We conclude that the high-velocity material and the bipolar shell must have originated during the late stages of evolution of a common-envelope phase in a close binary system.

We present the second and final data release of the Siding Spring Southern Seyfert Spectroscopic Snapshot Survey (S7). Data are presented for 63 new galaxies not included in the first data release, and we provide 2D emission-line fitting products for the full S7 sample of 131 galaxies. The S7 uses the WiFeS instrument on the ANU 2.3m telescope to obtain spectra with a spectral resolution of R = 7000 in the red (540 - 700 nm) and R = 3000 in the blue (350 - 570 nm), over an integral field of 25x38 arcsec^2 with 1x1 arcsec^2 spatial pixels. The S7 contains both the largest sample of active galaxies and the highest spectral resolution of any comparable integral field survey to date. The emission-line fitting products include line fluxes, velocities and velocity dispersions across the WiFeS field of view, and an artificial neural network has been used to determine the optimal number of Gaussian kinematic components for emission-lines in each spaxel. Broad Balmer lines are subtracted from the spectra of nuclear spatial pixels in Seyfert 1 galaxies before fitting the narrow lines. We bin nuclear spectra and measure reddening-corrected nuclear fluxes of strong narrow lines for each galaxy. The nuclear spectra are classified on optical diagnostic diagrams, where the strength of the coronal line [FeVII]6087 is shown to be correlated with [OIII]/Hbeta. Maps revealing gas excitation and kinematics are included for the entire sample, and we provide notes on the newly observed objects.

The most rapidly evolving regions of galaxies often display complex optical spectra with emission lines excited by massive stars, shocks and accretion onto supermassive black holes. Standard calibrations (such as for the star formation rate) cannot be applied to such mixed spectra. In this paper we isolate the contributions of star formation, shock excitation and active galactic nucleus (AGN) activity to the emission line luminosities of individual spatially resolved regions across the central 3 $\times$ 3 kpc$^2$ region of the active barred spiral galaxy NGC$\sim$613. The star formation rate and AGN luminosity calculated from the decomposed emission line maps are in close agreement with independent estimates from data at other wavelengths. The star formation component traces the B-band stellar continuum emission, and the AGN component forms an ionization cone which is aligned with the nuclear radio jet. The optical line emission associated with shock excitation is cospatial with strong $H_2$ and [Fe II] emission and with regions of high ionized gas velocity dispersion ($\sigma > 100$ km s$^{-1}$). The shock component also traces the outer boundary of the AGN ionization cone and may therefore be produced by outflowing material interacting with the surrounding interstellar medium. Our decomposition method makes it possible to determine the properties of star formation, shock excitation and AGN activity from optical spectra, without contamination from other ionization mechanisms.

Protostellar discs play an important role in star formation, acting as the primary mass reservoir for accretion onto young stars and regulating the extent to which angular momentum and gas is released back into stellar nurseries through the launching of powerful disc winds. In this study, we explore how disc structure relates to the properties of the wind-launching region, mapping out the regions of protostellar discs where wind launching could be viable. We combine a series of 1.5D semi-analytic, steady-state, vertical disc-wind solutions into a radially extended 1+1.5D model, incorporating all three diffusion mechanisms (Ohm, Hall and ambipolar). We observe that the majority of mass outflow via disc winds occurs over a radial width of a fraction of an astronomical unit, with outflow rates attenuating rapidly on either side. We also find that the mass accretion rate, magnetic field strength and surface density profile each have significant effects on both the location of the wind-launching region and the ejection/accretion ratio $\dot{M}_{\rm out}/\dot{M}_{\rm in}$. Increasing either the accretion rate or the magnetic field strength corresponds to a shift of the wind-launching region to smaller radii and a decrease in $\dot{M}_{\rm out}/\dot{M}_{\rm in}$, while increasing the surface density corresponds to launching regions at larger radii with increased $\dot{M}_{\rm out}/\dot{M}_{\rm in}$. Finally, we discover a class of disc winds containing an ineffective launching configuration at intermediate radii, leading to two radially separated regions of wind launching and diminished $\dot{M}_{\rm out}/\dot{M}_{\rm in}$. We find that the wind locations and ejection/accretion ratio are consistent with current observational and theoretical estimates.

We describe high spectral resolution, high dynamic range integral field spectroscopy of IC418 covering the spectral range 3300-8950{\AA} and compare with earlier data. We determine line fluxes, derive chemical abundances, provide a spectrum of the central star, and determine the shape of the nebular continuum. Using photoionisation models, we derive the reddening function from the nebular continuum and recombination lines. The nebula has a very high inner ionisation parameter. Consequently, radiation pressure dominates the gas pressure and dust absorbs a large fraction of ionising photons. Radiation pressure induces increasing density with radius. From a photoionisation analysis we derive central star parameters; $\log T_{\mathrm eff} = 4.525$K, $\log L_*/L_{\odot} = 4.029$, $\log g = 3.5$ and using stellar evolutionary models we estimate an initial mass of $2.5 < M/M_{\odot} < 3.0$. The inner filamentary shell is shocked by the rapidly increasing stellar wind ram pressure, and we model this as an externally photoionised shock. In addition, a shock is driven into the pre-existing Asymptotic Giant Branch stellar wind by the strong D-Type ionisation front developed at the outer boundary of the nebula. From the dynamics of the inner mass-loss bubble, and from stellar evolutionary models we infer that the nebula became ionised in the last $100-200$\,yr, but evolved structurally during the $\sim 2000$ yr since the central star evolved off the AGB. The estimated current mass loss rate ($\dot M = 3.8\times 10^{-8} M_{\odot}$yr$^{-1}$) and terminal velocity ($v_{\infty} \sim 450$ km/s) is sufficient to excite the inner mass-loss bubble. While on the AGB, the central star lost mass at $\dot M = 2.1\times 10^{-5} M_{\odot}$yr$^{-1}$ with outflow velocity $\sim 14$ km/s.

In an earlier paper we treated the pre-ionisation problem in shocks over the velocity range $20 < v_{\rm s} < 1000$km/s in a fully self-consistent manner. Here we investigate in detail the effect of the upstream UV photon field generated in the radiative zone of shocks in the range in which hydrogen is only partly ionised ($20 < v_{\rm s} < 150$km/s). We show that, as a result of super-heating in the non-equilibrium pre-shock plasma, both the magnetic parameter and the Mach number of the shock is strongly affected by the pre-ionisation state of the gas which controls to a large extent the radiative spectrum of the shock. We use these models to provide specific line diagnostics for Herbig-Haro objects which allow us to solve for both the pre-shock density and shock velocity, and present detailed models of the HH34 jet which allows us to derive the shock conditions, mass-loss rate, momentum flux and chemical abundances in the jet. We show that the refractory elements, Mg, Ca, Fe and Ni are enhanced by 0.22 dex over the solar values, which provides interesting clues about the jet launching mechanism in pre-main sequence evolution.

In this paper we treat the pre-ionisation problem in shocks over the velocity range $10 < v_{\rm s} < 1500$\,km/s in a self-consistent manner. We identify four distinct classes of solution controlled by the value of the shock precursor parameter, $\Psi = {\cal Q}/v_s$, where ${\cal Q}$ is the ionization parameter of the UV photons escaping upstream. This parameter determines both the temperature and the degree of ionisation of the gas entering the shock. In increasing velocity the shock solution regimes are cold neutral precursors ($v_s \lesssim 40$\,km/s), warm neutral precursors ($40 \lesssim v_s \lesssim 75$\,km/s), warm partly-ionized precursors ($75 \lesssim v_s \lesssim 120$\,km/s), and fast shocks in which the pre-shock gas is in photoionisation equilibrium, and is fully ionized. The main effect of a magnetic field is to push these velocity ranges to higher values, and to limit the post-shock compression. In order to facilitate comparison with observations of shocks, we provide a number of convenient scaling relationships for parameters such as post-shock temperature, compression factors, cooling lengths, and H$\beta$ and X-ray luminosity.

We present a new basis for scaling abundances with total metallicity in nebular photoionisation models, based on extensive Milky Way stellar abundance data, to replace the uniform scaling normally used in the analysis of HII regions. Our goal is to provide a single scaling method and local abundance reference standard for use in nebular modelling and its key inputs, the stellar atmosphere and evolutionary track models. We introduce a parametric enrichment factor, zeta, to describe how atomic abundances scale with total abundance, and which allows for a simple conversion between scales based on different reference elements (usually oxygen or iron) . The models and parametric description provide a more physically realistic approach than simple uniform abundance scaling. With appropriate parameters, the methods described here may be applied to HII regions in the Milky Way, large and dwarf galaxies in the local universe, Active Galactic Nuclei (AGNs), and to star forming regions at high redshift.

We present a simplified model of Active Galactic Nucleus (AGN) continuum emission designed for photoionization modeling. The new model {\sc oxaf} reproduces the diversity of spectral shapes that arise in physically-based models. We identify and explain degeneracies in the effects of AGN parameters on model spectral shapes, with a focus on the complete degeneracy between the black hole mass and AGN luminosity. Our re-parametrized model {\sc oxaf} removes these degeneracies and accepts three parameters which directly describe the output spectral shape: the energy of the peak of the accretion disk emission $E_\mathrm{peak}$, the photon power-law index of the non-thermal emission $\Gamma$, and the proportion of the total flux which is emitted in the non-thermal component $p_\mathrm{NT}$. The parameter $E_\mathrm{peak}$ is presented as a function of the black hole mass, AGN luminosity, and `coronal radius' of the {\sc optxagnf} model upon which {\sc oxaf} is based. We show that the soft X-ray excess does not significantly affect photoionization modeling predictions of strong emission lines in Seyfert narrow-line regions. Despite its simplicity, {\sc oxaf} accounts for opacity effects where the accretion disk is ionized because it inherits the `color correction' of {\sc optxagnf}. We use a grid of {\sc mappings} photoionization models with {\sc oxaf} ionizing spectra to demonstrate how predicted emission-line ratios on standard optical diagnostic diagrams are sensitive to each of the three {\sc oxaf} parameters. The {\sc oxaf} code is publicly available in the Astrophysics Source Code Library.

The optical spectra of Seyfert galaxies are often dominated by emission lines excited by both star formation and AGN activity. Standard calibrations (such as for the star formation rate) are not applicable to such composite (mixed) spectra. In this paper, we describe how integral field data can be used to spectrally and spatially separate emission associated with star formation from emission associated with accretion onto an active galactic nucleus (AGN). We demonstrate our method using integral field data for two AGN host galaxies (NGC 5728 and NGC 7679) from the Siding Spring Southern Seyfert Spectroscopic Snapshot Survey (S7). The spectra of NGC 5728 and NGC 7679 form clear sequences of AGN fraction on standard emission line ratio diagnostic diagrams. We show that the emission line luminosities of the majority (> 85 per cent) of spectra along each AGN fraction sequence can be reproduced by linear superpositions of the emission line luminosities of one AGN dominated spectrum and one star formation dominated spectrum. We separate the Halpha, Hbeta, [N II]$\lambda$6583, [S II]$\lambda \lambda$6716, 6731, [O III]$\lambda$5007 and [O II]$\lambda \lambda$3726, 3729 luminosities of every spaxel into contributions from star formation and AGN activity. The decomposed emission line images are used to derive the star formation rates and AGN bolometric luminosities for NGC 5728 and NGC 7679. Our calculated values are mostly consistent with independent estimates from data at other wavelengths. The recovered star forming and AGN components also have distinct spatial distributions which trace structures seen in high resolution imaging of the galaxies, providing independent confirmation that our decomposition has been successful.

We present the results of three dimensional relativistic hydrodynamic simulations of interaction of AGN jets with a dense turbulent two-phase interstellar medium, which would be typical of high redshift galaxies. We describe the effect of the jet on the evolution of the density of the turbulent ISM. The jet driven energy bubble affects the gas to distances up to several kiloparsecs from the injection region. The shocks resulting from such interactions create a multi-phase ISM and radial outflows. One of the striking result of this work is that low power jets (P_jet < 10^{43} erg/s) although less efficient in accelerating clouds, are trapped in the ISM for a longer time and hence affect the ISM over a larger volume. Jets of higher power drill through with relative ease. Although the relativistic jets launch strong outflows, there is little net mass ejection to very large distances, supporting a galactic fountain scenario for local feedback.

We report here on the discovery of faint extended wings of H\alpha\ observed out to an apparent velocity of ~ 7600 km/s in the Orion Nebula (M42) and in five HII regions in the Large and the Small Magellanic Clouds. We show that, these wings are caused by Raman scattering of both the O I and Si II resonance lines and stellar continuum UV photons with H I followed by radiative decay to the H I n=2 level. The broad wings also seen in H\beta\ and in H\gamma\ result from Raman scattering of the UV continuum in the H I n=4 and n=5 levels respectively.The Raman scattering fluorescence is correlated with the intensity of the narrow permitted lines of O I and Si II. In the case of Si II, this is explained by radiative pumping of the same 1023.7\AA\ resonance line involved in the Raman scattering by the Ly\beta\ radiation field. The subsequent radiative cascade produces enhanced Si II 5978.9, 6347.1$ and 6371.4\AA\ permitted transitions. Finally we show that in O I, radiative pumping of the 1025.76\AA\ resonance line by the Lyman series radiation field is also the cause of the enhancement in the permitted lines of this species lying near H\alpha\ in wavelength, but here the process is a little more complex. We argue that all these processes are active in the zone of the HII region near the ionisation front

We present results of a complete integral field survey of the bright SNR N49 in the LMC, obtained with the WiFeS instrument mounted on the ANU 2.3m telescope at Siding Spring Observatory. From theoretical shock modelling with the new MAPPINGS 5.1 code we have, for the first time, subjected the optical Fe emission line spectrum of a supernova remnant to a detailed abundance and dynamical analysis covering 8 separate stages of ionisation. This allows us to derive the dust depletion factors as a function of ionisation stage. We have shown that there is substantial (30% - 90%) destruction of Fe-bearing dust grains in these fast shocks ($v_s \sim 250$ km/s), and we have confirmed that the dominant dust destruction is through the non-thermal sputtering and grain-grain collision mechanisms developed in a number of theoretical works.

We investigate the relative significance of radiation pressure and gas pressure in the extended narrow line regions (ENLRs) of four Seyfert galaxies from the integral field Siding Spring Southern Seyfert Spectroscopic Snapshot Survey (S7). We demonstrate that there exist two distinct types of starburst-AGN mixing curves on standard emission line diagnostic diagrams which reflect the balance between gas pressure and radiation pressure in the ENLR. In two of the galaxies the ENLR is radiation pressure dominated throughout and the ionization parameter remains constant (log U ~ 0). In the other two galaxies radiation pressure is initially important, but gas pressure becomes dominant as the ionization parameter in the ENLR decreases from log U ~ 0 to -3.4 <= log U <= -3.2. Where radiation pressure is dominant, the AGN regulates the density of the interstellar medium on kpc scales and may therefore have a direct impact on star formation activity and/or the incidence of outflows in the host galaxy to scales far beyond the zone of influence of the black hole. We find that both radiation pressure dominated and gas pressure dominated ENLRs are dynamically active with evidence for outflows, indicating that radiation pressure may be an important source of AGN feedback even when it is not dominant over the entire ENLR.

We present the results of a meeting on numerical simulations of ionized nebulae held at the University of Kentucky in conjunction with the celebration of the 70th birthdays of Profs. Donald Osterbrock and Michael Seaton.

We present three dimensional relativistic hydrodynamical simulations of a precessing jet interacting with the intracluster medium and compare the simulated jet structure with the observed structure of the Hydra A northern jet. For the simulations, we use jet parameters obtained in the parameter space study of the first paper in this series and probe different values for the precession period and precession angle. We find that for a precession period P = 1 Myr and a precession angle = 20 degree the model reproduces i) the curvature of the jet, ii) the correct number of bright knots within 20 kpc at approximately correct locations, and iii) the turbulent transition of the jet to a plume. The Mach number of the advancing bow shock = 1.85 is indicative of gentle cluster atmosphere heating during the early stages of the AGN's activity.

This Letter presents a new, remarkably simple diagnostic specifically designed to derive chemical abundances for high redshift galaxies. It uses only the H \alpha, [N II] and [S II] emission lines, which can usually be observed in a single gating stetting, and is almost linear up to an abundance of 12+log (O/H) = 9.05. It can be used over the full abundance range encountered in high redshift galaxies. By its use of emission lines located close together in wavelength, it is also independent of reddening. Our diagnostic depends critically on the calibration of the N/O ratio. However, by using realistic stellar atmospheres combined with the N/O vs. O/H abundance calibration derived locally from stars and H II regions, and allowing for the fact that high-redshift H II regions have both high ionisation parameters \emph{and} high gas pressures, we find that the observations of high-redshift galaxies can be simply explained by the models without having to invoke arbitrary changes in N/O ratio, or the presence of unusual quantities of Wolf-Rayet stars in these galaxies.

Powerful relativistic jets in radio galaxies are capable of driving strong outflows but also inducing star-formation by pressure-triggering collapse of dense clouds. We review theoretical work on negative and positive active galactic nuclei feedback, discussing insights gained from recent hydrodynamical simulations of jet-driven feedback on galaxy scales that are applicable to compact radio sources. The simulations show that the efficiency of feedback and the relative importance of negative and positive feedback depends strongly on interstellar medium properties, especially the column depth and spatial distribution of clouds. Negative feedback is most effective if clouds are distributed spherically and individual clouds have small column depths, while positive feedback is most effective if clouds are predominantly in a disc-like configuration.

Turbulent entrainment processes may play an important role in the outflows from young stellar objects at all stages of their evolution. In particular, lateral entrainment of ambient material by high-velocity, well-collimated protostellar jets may be the cause of the multiple emission-line velocity components observed in the microjet-scale outflows driven by classical T Tauri stars. Intermediate-velocity outflow components may be emitted by a turbulent, shock- excited mixing layer along the boundaries of the jet. We present a formalism for describing such a mixing layer based on Reynolds decomposition of quantities measuring fundamental properties of the gas. In this model, the molecular wind from large disc radii provides a continual supply of material for entrainment. We calculate the total stress profile in the mixing layer, which allows us to estimate the dissipation of turbulent energy, and hence the luminosity of the layer. We utilize MAPPINGS IV shock models to determine the fraction of total emission that occurs in [Fe II] 1.644 {\mu}m line emission in order to facilitate comparison to previous observations of the young stellar object DG Tauri. Our model accurately estimates the luminosity and changes in mass outflow rate of the intermediate-velocity component of the DG Tau approaching outflow. Therefore, we propose that this component represents a turbulent mixing layer surrounding the well-collimated jet in this object. Finally, we compare and contrast our model to previous work in the field.

The Magellanic HI Stream (~2x10^9 Msun [d/55 kpc]^2) encircling the Galaxy at a distance 'd' is arguably the most important tracer of what happens to gas accreting onto a disk galaxy. Recent observations reveal that the Stream's mass is in fact dominated (3:1) by its ionised component. Here we revisit the origin of the mysterious H-alpha recombination emission observed along much of its length that is overly bright (~150-200 milli-Rayleigh) for the known Galactic ultraviolet background (~20-40 mR / [d/55 kpc]^2). In an earlier model, we proposed that a slow shock cascade was operating along the Stream due to its interaction with the extended Galactic hot corona. We find that, for a smooth coronal density profile, this model can explain the bright H-alpha emission if the coronal density satisfies 2 < (n / 10^{-4} cm^{-3}) < 4 at d = 55 kpc. But in view of updated parameters for the Galactic halo and mounting evidence that most of the Stream must lie far beyond the Magellanic Clouds (d>55 kpc), we revisit the shock cascade model in detail. At lower densities, the HI gas is broken down by the shock cascade but mostly mixes with the hot corona without significant recombination. At higher densities, the hot coronal mass (including the other baryonic components) exceeds the baryon budget of the Galaxy. If the H-alpha emission arises from the shock cascade, the upper limit on the smooth coronal density constrains the Stream's mean distance to < 75 kpc. If, as some models indicate, the Stream is even further out, either the shock cascade is operating in a regime where the corona is substantially mass-loaded with recent gas debris, or an entirely different ionization mechanism is responsible.

The Siding Spring Southern Seyfert Spectroscopic Snapshot Survey (S7) is a targeted survey probing the narrow-line regions (NLRs) of a representative sample of ~140 nearby (z<0.02) Seyfert galaxies by means of optical integral field spectroscopy. The survey is based on a homogeneous data set observed using the Wide Field Spectrograph WiFeS. The data provide a 25x38 arcsec$^2$ field-of-view around the galaxy centre at typically ~1.5 arcsec spatial resolution and cover a wavelength range between ~3400 - 7100 $\AA$ at spectral resolutions of ~100 km s$^{-1}$ and ~50 km s$^{-1}$ in the blue and red parts, respectively. The survey is primarily designed to study gas excitation and star formation around AGN, with a special focus on the shape of the AGN ionising continuum, the interaction between radio jets and the NLR gas, and the nature of nuclear LINER emission. We provide an overview of the current status of S7-based results and present new results for NGC 6300.

New observations indicate that ultrafaint dwarf galaxies (UFD) -- the least luminous systems bound by dark matter halos (<10^5 Lsun) -- may have formed before reionization. The extrapolated virial masses today are uncertain with estimates ranging from 10^8 Msun to 10^9 Msun. We show that the progenitor halo masses of UFDs can be as low as Mvir = 10^7 Msun. Under the right conditions, such a halo can survive the energy input of a supernova and its radiative progenitor. A clumpy medium is much less susceptible to both internal and external injections of energy. It is less prone to SN sweeping because the coupling efficiency of the explosive energy is much lower than for a diffuse ISM. With the aid of the 3D hydro/ionization code Fyris, we show that sufficient baryons are retained to form stars following a single supernova event in dark matter halos down to Mvir ~ 10^7 Msun with radiative cooling. The gas survives the SN explosion, is enriched with the abundance yields of the discrete events, and reaches surface densities where low mass stars can form. Our highest resolution simulations reveal why cooling is so effective in retaining gas compared to any other factor. In the early stages, the super-hot metal-enriched SN ejecta exhibit strong cooling, leading to much of the explosive energy being lost. Consistent with earlier work, the baryons do *not* survive in smooth or adiabatic models in the event of a supernova. The smallest galaxies carry signatures of the earliest epochs of star formation, which may distinguish a small primordial galaxy from one that was stripped down to its present size through tidal interaction. We discuss these results in the context of local UFDs and damped Ly-alpha systems (z~2) at very low metallicity ([Fe/H] ~ -3). We show that both classes of objects are consistent with primordial low-mass systems that have experienced only a few enrichment events.

The density variance - Mach number relation of the turbulent interstellar medium is relevant for theoretical models of the star formation rate, efficiency, and the initial mass function of stars. Here we use high-resolution hydrodynamical simulations with grid resolutions of up to 1024^3 cells to model compressible turbulence in a regime similar to the observed interstellar medium. We use Fyris Alpha, a shock-capturing code employing a high-order Godunov scheme to track large density variations induced by shocks. We investigate the robustness of the standard relation between the logarithmic density variance (sigma_s^2) and the sonic Mach number (M) of isothermal interstellar turbulence, in the non-isothermal regime. Specifically, we test ideal gases with diatomic molecular (gamma = 7/5) and monatomic (gamma = 5/3) adiabatic indices. A periodic cube of gas is stirred with purely solenoidal forcing at low wavenumbers, leading to a fully-developed turbulent medium. We find that as the gas heats in adiabatic compressions, it evolves along the relationship in the density variance - Mach number plane, but deviates significantly from the standard expression for isothermal gases. Our main result is a new density variance - Mach number relation that takes the adiabatic index into account: sigma_s^2 = ln {1+b^2*M^[(5*gamma+1)/3]} and provides good fits for b*M <= 1. A theoretical model based on the Rankine-Hugoniot shock jump conditions is derived, sigma_s^2 = ln {1+(gamma+1)*b^2*M^2/[(gamma-1)*b^2*M^2+2]}, and provides good fits also for b*M > 1. We conclude that this new relation for adiabatic turbulence may introduce important corrections to the standard relation, if the gas is not isothermal.

In earlier work we showed that a dark matter halo with a virial mass of $10^7$ M$_\odot$ can survive feedback from its own massive stars and form stars for $\gtrsim100$ Myr. We also found that our modelled systems were consistent with observations of ultrafaint dwarfs (UFDs), the least massive known galaxies. Very metal-poor damped Lyman-$\alpha$ systems (DLAs) recently identified at $z\sim2$ may represent the gas that formed at least some of the observed stars in UFDs. We compare projected sightlines from our simulations to the observed metal-poor DLAs and find that our models can reach the densities of the observed sightlines; however the metallicities are inconsistent with the single supernova simulations, suggesting enrichment by multiple supernovae. We model two scenarios for the history of these systems. The first explains the gas abundances in DLAs by a single burst of star formation. This model can produce the observed DLA abundances, but does not provide an explanation as to why the DLAs show suppressed [$\alpha$/Fe] compared to the stellar population of UFDs. The second scenario splits the DLAs into a population which is enriched by a single burst, and a population that is enriched by a second burst after the accretion of metal-poor gas. In this scenario, the suppressed average [$\alpha$/Fe] in DLAs compared to UFDs results from enrichment of second-burst systems by Type Ia supernovae.

We combine optical and near-infrared AO-assisted integral field observations of the merging ULIRG IRAS F17207-0014 from the Wide-Field Spectrograph (WiFeS) and Keck/OSIRIS. The optical emission line ratios [N II]/H$\alpha$, [S II]/H$\alpha$, and [O I]/H$\alpha$ reveal a mixing sequence of shocks present throughout the galaxy, with the strongest contributions coming from large radii (up to 100% at $\sim$5 kpc in some directions), suggesting galactic-scale winds. The near-infrared observations, which have approximately 30 times higher spatial resolution, show that two sorts of shocks are present in the vicinity of the merging nuclei: low-level shocks distributed throughout our field-of-view evidenced by an H$_{2}$/Br$\gamma$ line ratio of $\sim$0.6-4, and strong collimated shocks with a high H$_{2}$/Br$\gamma$ line ratio of $\sim$4-8, extending south from the two nuclear disks approximately 400 pc ($\sim$0.5 arcsec). Our data suggest that the diffuse shocks are caused by the collision of the interstellar media associated with the two progenitor galaxies and the strong shocks trace the base of a collimated outflow coming from the nucleus of one of the two disks.

We model the chemical evolution of six UFDs: Bootes I, Canes Venatici II, Coma Berenices, Hercules, Leo IV and Ursa Major I, based on their recently determined star formation histories. We show that two single-age bursts cannot explain the observed [$\alpha$/Fe] vs [Fe/H] distribution in these galaxies and that some self-enrichment is required within the first burst. An alternative scenario is modelled, in which star formation is continuous except for short interruptions when one or more supernovae temporarily blow the dense gas out from the centre of the system. This model allows for self-enrichment and can reproduce the chemical abundances of the UFDs in which the second burst is only a trace population. We conclude that the most likely star formation history is one or two extended periods of star formation, with the first burst lasting for at least 100 Myr. As found in earlier work, the observed properties of UFDs can be explained by formation at a low mass ($M_{\rm{vir}}\sim10^7$ M$_\odot$), rather than being stripped remnants of much larger systems.

We present Wide Field Spectrograph (WiFeS) integral field spectroscopy and HST FOS spectroscopy for the LINER galaxy NGC 1052. We infer the presence of a turbulent accretion flow forming a small-scale accretion disk. We find a large-scale outflow and ionisation cone along the minor axis of the galaxy. Part of this outflow region is photoionised by the AGN, and shares properties with the ENLR of Seyfert galaxies, but the inner ($R \lesssim 1.0$~arcsec) accretion disk and the region around the radio jet appear shock excited. The emission line properties can be modelled by a "double shock" model in which the accretion flow first passes through an accretion shock in the presence of a hard X-ray radiation, and the accretion disk is then processed through a cocoon shock driven by the overpressure of the radio jets. This model explains the observation of two distinct densities ($\sim10^4$ and $\sim10^6$ cm$^{-3}$), and provides a good fit to the observed emission line spectrum. We derive estimates for the velocities of the two shock components and their mixing fractions, the black hole mass, the accretion rate needed to sustain the LINER emission and derive an estimate for the jet power. Our emission line model is remarkably robust against variation of input parameters, and so offers a generic explanation for the excitation of LINER galaxies, including those of spiral type such as NGC 3031 (M81).

Here we describe the \emph{Siding Spring Southern Seyfert Spectroscopic Snapshot Survey} (S7) and present results on 64 galaxies drawn from the first data release. The S7 uses the Wide Field Spectrograph (WiFeS) mounted on the ANU 2.3m telescope located at the Siding Spring Observatory to deliver an integral field of $38\times25$~ arcsec at a spectral resolution of $R=7000$ in the red ($530-710$nm), and $R=3000$ in the blue ($340-560$nm). {From these data cubes we have extracted the Narrow Line Region (NLR) spectra from a 4 arc sec aperture centred on the nucleus. We also determine the H$\beta$ and [OIII]~$\lambda$5007 fluxes in the narrow lines, the nuclear reddening, the reddening-corrected relative intensities of the observed emission lines, and the H$\beta$ and \lOIII\ luminosities {determined from spectra for which the stellar continuum has been removed.} We present a set of images of the galaxies in [OIII]~$\lambda$5007, [NII]~$\lambda$6584 and H$\alpha$ which serve to delineate the spatial extent of the extended narrow line region (ENLR) and {\bf also to} reveal the structure and morphology of the surrounding \HII\ regions. Finally, we provide a preliminary discussion of those Seyfert~1 and Seyfert~2 galaxies which display coronal emission lines in order to explore the origin of these lines.

In earlier work we showed that a dark matter halo with a virial mass of $10^7$ M$_\odot$ can retain a large percentage of its baryons in the face of the pre-ionization phase and supernova explosion from a 25$M_\odot$ star. Here we expand on the results of that work, investigating the star formation and chemical evolution of the system beyond the first supernova. In a galaxy with a mass $M_{vir} = 10^7$M$_\odot$, sufficient gas is retained by the potential for a second period of star formation to occur. The impact of a central explosion is found to be much stronger than that of an off-centred explosion both in blowing out the gas and in enriching it, as in the off-centered case most of the supernova energy and metals escape into the IGM. We model the star formation and metallicity given the assumption that stars form for 100, 200, 400 and 600 Myr and discuss the results in the context of recent observations of very low mass galaxies. We show that we can account for most features of the observed relationship between [$\alpha$/Fe] and [Fe/H] in ultra-faint dwarf galaxies with the assumption that the systems formed at a low mass, rather than being remnants of much larger systems.

The expanding remnant from SN 1987A is an excellent laboratory for investigating the physics of supernovae explosions. There are still a large number of outstanding questions, such the reason for the asymmetric radio morphology, the structure of the pre-supernova environment, and the efficiency of particle acceleration at the supernova shock. We explore these questions using three-dimensional simulations of the expanding remnant between days 820 and 10,000 after the supernova. We combine a hydrodynamical simulation with semi-analytic treatments of diffusive shock acceleration and magnetic field amplification to derive radio emission as part of an inverse problem. Simulations show that an asymmetric explosion, combined with magnetic field amplification at the expanding shock, is able to replicate the persistent one-sided radio morphology of the remnant. We use an asymmetric Truelove & McKee progenitor with an envelope mass of $10 M_{\sun}$ and an energy of $1.5 \times 10^{44} J$. A termination shock in the progenitor's stellar wind at a distance of $0\farcs43-0\farcs51$ provides a good fit to the turn on of radio emission around day 1200. For the H\textsc{ii} region, a minimum distance of $0\farcs63\pm0\farcs01$ and maximum particle number density of $(7.11\pm1.78) \times 10^7$ m$^{-3}$ produces a good fit to the evolving average radius and velocity of the expanding shocks from day 2000 to day 7000 after explosion. The model predicts a noticeable reduction, and possibly a temporary reversal, in the asymmetric radio morphology of the remnant after day 7000, when the forward shock left the eastern lobe of the equatorial ring.

Here we present the first results from the Siding Spring Southern Seyfert Spectroscopic Snapshot Survey (S7) which aims to investigate the physics of ~140 radio-detected southern active Galaxies with z<0.02 through Integral Field Spectroscopy using the Wide Field Spectrograph (WiFeS). This instrument provides data cubes of the central 38 x 25 arc sec. of the target galaxies in the waveband 340-710nm with the unusually high resolution of R=7000 in the red (530-710nm), and R=3000 in the blue (340-560nm). These data provide the morphology, kinematics and the excitation structure of the extended narrow-line region, probe relationships with the black hole characteristics and the host galaxy, measures host galaxy abundance gradients and the determination of nuclear abundances from the HII regions. From photoionisation modelling, we may determine the shape of the ionising spectrum of the AGN, discover whether AGN metallicities differ from nuclear abundances determined from HII regions, and probe grain destruction in the vicinity of the AGN. Here we present some preliminary results and modelling of both Seyfert galaxies observed as part of the survey.

We present the first stage of an investigation of the interactions of the jets in the radio galaxy Hydra A with the intracluster medium. We consider the jet kinetic power, the galaxy and cluster atmosphere, and the inner structure of the radio source. Analysing radio observations of the inner lobes of Hydra A by Taylor et al. (1990) we confirm the jet power estimates of about 1e45 ergs/s derived by Wise et al. (2007) from dynamical analysis of the X-ray cavities. With this result and a model for the galaxy halo, we explore the jet-intracluster medium interactions occurring on a scale of 10 kpc using two-dimensional, axisymmetric, relativistic pure hydrodynamic simulations. A key feature is that we identify the three bright knots in the northern jet as biconical reconfinement shocks, which result when an over pressured jet starts to come into equilibrium with the galactic atmosphere. Through an extensive parameter space study we deduce that the jet velocity is approximately 0.8 c at a distance 0.5 kpc from the black hole. The combined constraints of jet power, the observed jet radius profile along the jet, and the estimated jet pressure and jet velocity imply a value of the jet density parameter approximately 13 for the northern jet. We show that for a jet velocity = 0.8c and angle between the jet and the line of sight = 42 deg, an intrinsic asymmetry in the emissivity of the northern and southern jet is required for a consistent brightness ratio approximately 7 estimated from the 6cm VLA image of Hydra A.

The Fermi Bubbles are enigmatic \gamma-ray features of the Galactic bulge. Both putative activity (within $\sim$ few $\times$ Myr) connected to the Galactic center super-massive black hole and, alternatively, nuclear star formation have been claimed as the energising source of the Bubbles. Likewise, both inverse-Compton emission by non-thermal electrons (`leptonic' models) and collisions between non-thermal protons and gas (`hadronic' models) have been advanced as the process supplying the Bubbles' \gamma -ray emission. An issue for any steady state hadronic model is that the very low density of the Bubbles' plasma seems to require that they accumulate protons over a multi-Gyr timescale, much longer than other natural timescales occurring in the problem. Here we present a hadronic model where the timescale for generating the Bubbles' hadronic \gamma -ray emission is $\sim$ few $\times 10^8$ years. Our model invokes collapse of the Bubbles' thermally-unstable plasma, leading to an accumulation of cosmic rays and magnetic field into localised, warm ($\sim 10^4$ K), and likely filamentary condensations of higher density gas. Under the condition that these filaments are supported by non-thermal pressure, we can predict the hadronic emission from the Bubbles to be $L_\gamma \simeq 2 \times 10^{37}$ erg/s $\ \dot{M}_\mathrm{in}/(0.1 \ M_{Sun}/$ year $) \ T_\mathrm{FB}^2/(3.5 \times 10^7 K) ^2 M_{fil}/M_{pls}$ ; precisely their observed luminosity (normalizing to the star-formation-driven mass flux into the Bubbles and their measured plasma temperature and adopting the further result that the mass in the filaments, $M_{fil}$ is approximately equal to that of the Bubbles' plasma, $M_{pls}$).

Two-dimensional (2D) line ratio diagnostic diagrams have become a key tool in understanding the excitation mechanisms of galaxies. The curves used to separate the different regions - HII-like or else excited by an active galactic nucleus (AGN) - have been refined over time but the core technique has not evolved significantly. However, the classification of galaxies based on their emission line ratios really is a multi-dimensional problem. Here we exploit recent software developments to explore the potential of three-dimensional (3D) line ratio diagnostic diagrams. We introduce a specific set of 3D diagrams, the ZQE diagrams, which separate the oxygen abundance and the ionisation parameter of HII region-like spectra, and which also enable us to probe the excitation mechanism of the gas. By examining these new 3D spaces interactively, we define a new set of 2D diagnostics, the ZE diagnostics, which can provide the metallicity of objects excited by hot young stars, and which cleanly separate HII region-like objects from the different classes of AGNs. We show that these ZE diagnostics are consistent with the key log[NII]/H$\alpha$ vs. log[OIII]/H$\beta$ diagnostic currently used by the community. They also have the advantage of attaching a probability that a given object belongs to one class or to the other. Finally, we discuss briefly why ZQE diagrams can provide a new way to differentiate and study the different classes of AGNs in anticipation of a dedicated follow-up study.

In this paper, we use the Mappings photoionization code to explore the physical parameters that impact on the measurement of electron temperature and abundance in HII regions. In the previous paper we presented observations and measurements of physical properties from the spectra of seventeen HII regions in fourteen isolated dwarf irregular galaxies from the SIGRID sample. Here, we analyze these observations further, together with three additional published data sets. We explore the effects of optical thickness, electron density, ionization parameter, ionization source, and non-equilibrium effects on the relation between electron temperature and metallicity. We present a standard model that fits the observed data remarkably well at metallicities between 1/10 and 1 solar. We investigate the effects of optically thin HII regions, and show that they can have a considerable effect on the measured electron temperature, and that there is evidence that some of the observed objects are optically thin. We look at the role of the ionization parameter and find that lower ionization parameter values give better fits at higher oxygen abundance. We show that higher pressures combined with low optical depth, and also kappa electron energy distributions at low kappa values, can generate the apparent high electron temperatures in low metallicity HII regions, and that the former provides the better fit to observations. We examine the effects of these parameters on the strong line diagnostic methods. We extend this to three-dimensional diagnostic grids to confirm how well the observations are described by the grids.

We have used the Wide Field Spectrograph (WiFeS) on the ANU 2.3m telescope at Siding Spring to observe the nearby, nearly face-on, Seyfert 2 galaxy, NGC 5427. We have obtained integral field spectroscopy of both the nuclear regions and the HII regions in the spiral arms. We have constrained the chemical abundance in the interstellar medium of the extended narrow line region (ENLR) by measuring the abundance gradient in the circum-nuclear \ion{H}{ii} regions to determine the nuclear chemical abundances, and to use these to in turn determine the EUV spectral energy distribution for comparison with theoretical models. We find a very high nuclear abundance, $\sim 3.0$ times solar, with clear evidence of a nuclear enhancement of N and He, possibly caused by massive star formation in the extended ($\sim 100$pc) central disk structure. The circum-nuclear narrow-line region spectrum is fit by a radiation pressure dominated photoionisation model model with an input EUV spectrum from a Black Hole with mass $5\times10^7 M_{\odot}$ radiating at $\sim 0.1$ of its Eddington luminosity. The bolometric luminosity is closely constrained to be $\log L_{\mathrm bol.} = 44.3\pm 0.1$ erg s$^{-1}$. The EUV spectrum characterised by a soft accretion disk and a harder component extending to above 15keV. The ENLR region is extended in the NW-SE direction. The line ratio variation in circum-nuclear spaxels can be understood as the result of mixing \ion{H}{ii} regions with an ENLR having a radius-invariant spectrum.

In this paper we present the results of observations of seventeen HII regions in thirteen galaxies from the SIGRID sample of isolated gas rich irregular dwarf galaxies. The spectra of all but one of the galaxies exhibit the auroral [OIII] 4363A line, from which we calculate the electron temperature, Te, and gas-phase oxygen abundance. Five of the objects are blue compact dwarf (BCD) galaxies, of which four have not previously been analysed spectroscopically. We include one unusual galaxy which exhibits no evidence of the [NII] {\lambda}{\lambda} 6548,6584A lines, suggesting a particularly low metallicity (< Zsolar/30). We compare the electron temperature based abundances with those derived using eight of the new strong line diagnostics presented by Dopita et al. (2013). Using a method derived from first principles for calculating total oxygen abundance, we show that the discrepancy between the Te-based and strong line gas-phase abundances have now been reduced to within ~0.07 dex. The chemical abundances are consistent with what is expected from the luminosity-metallicity relation. We derive estimates of the electron densities and find them to be between ~5 and ~100 cm-3. We find no evidence for a nitrogen plateau for objects in this sample with metallicities 0.5 > Zsolar > 0.15.

The Fermi satellite discovery of the gamma-ray emitting bubbles extending 50 deg (10 kpc) from the Galactic Centre has revitalized earlier claims that our Galaxy has undergone an explosive episode in the recent past. We now explore a new constraint on such activity. The Magellanic Stream is a clumpy gaseous structure passing over the South Galactic Pole (SGP) at a distance of at least 50-100 kpc. Patchy H-alpha emission discovered along the Magellanic Stream over the SGP is a factor of 5 too bright to have been produced by the Galactic stellar population. Time-dependent models of Stream clouds exposed to a flare in ionising photon flux show that the ionised gas must recombine and cool for a time interval 0.6 - 2.9 Myr for the emitted H-alpha surface brightness to drop to the observed level. A nuclear starburst is ruled out. Sgr A* is a more likely candidate because it is two orders of magnitude more efficient at converting gas to UV radiation. The central black hole (4 x 10^6 Msun) can supply the required ionising luminosity with a fraction of the Eddington accretion rate (3-30%). In support of nuclear activity, the H-alpha emission along the Stream has a polar angle dependence peaking close to the SGP. Moreover, it is now generally accepted that the Stream over the SGP must be further than the Magellanic Clouds. At the lower halo gas densities, shocks become too ineffective and are unlikely to give rise to a polar angle dependence in the H-alpha emission. Thus it is likely that the Stream emission arose from a `Seyfert flare' that was active 1-3 Myr ago, consistent with the cosmic ray lifetime in the Fermi bubbles. Sgr A* activity today is greatly suppressed (70-80 dB) relative to the Seyfert outburst...

Recently, \citet{Nicholls12}, inspired by \emph{in situ} observations of solar system astrophysical plasmas, suggested that the electrons in \HII regions are characterised by a $\kappa$-distribution of electron energies rather than by a simple Maxwell-Boltzmann distribution. Here we have collected together the new atomic data within a modified photoionisation code to explore the effects of both the new atomic data and the $\kappa$-distribution on the strong-line techniques used to determine chemical abundances in \HII regions. By comparing the recombination temperatures ($T_{\rm rec}$) with the forbidden line temperatures ($T_{\rm FL}$) we conclude that $ \kappa \sim 20$. While representing only a mild deviation from equilibrium, this is sufficient to strongly influence abundances determined using methods which depend on measurements of the electron temperature from forbidden lines. We present a number of new emission line ratio diagnostics which cleanly separate the two parameters determining the optical spectrum of \HII regions - the ionisation parameter $q$ or $\cal{U}$ and the chemical abundance; 12+log(O/H). An automated code to extract these parameters is presented. Using the homogeneous dataset from \citet{vanZee98}, we find self-consistent results between all these different diagnostics. The systematic errors between different line ratio diagnostics are much smaller than was found in the earlier strong line work. Overall the effect of the $\kappa$-distribution on the strong line abundances derived solely on the basis of theoretical models is rather small.

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, [NII]/H-alpha and [OIII]/H-beta. 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.

In this paper we develop tools for observers to use when analysing nebular spectra for temperatures and metallicities, with two goals: to present a new, simple method to calculate equilibrium electron temperatures for collisionally excited line flux ratios, using the latest atomic data; and to adapt current methods to include the effects of possible non-equilibrium '{\kappa}' electron energy distributions. Adopting recent collision strength data for [O iii], [S iii], [O ii], [S ii], and [N ii], we find that existing methods based on older atomic data seriously overestimate the electron temperatures, even when considering purely Maxwellian statistics. If {\kappa} distributions exist in H ii regions and planetary nebulae as they do in solar system plasmas, it is important to investigate the observational consequences. This paper continues our previous work on the {\kappa} distribution (Nicholls et al. 2012). We present simple formulaic methods that allow observers to (a) measure equilibrium electron temperatures and atomic abundances using the latest atomic data, and (b) to apply simple corrections to existing equilibrium analysis techniques to allow for possible non-equilibrium effects. These tools should lead to better consistency in temperature and abundance measurements, and a clearer understanding of the physics of H ii regions and planetary nebulae.

Neumayer et al. established the existence of a blue-shifted cloud in the core of Centaurus A, within a few parsecs of the nucleus and close to the radio jet. We propose that the cloud has been impacted by the jet, and that it is in the foreground of the jet, accounting for its blue-shifted emission on the Southern side of the nucleus. We consider both shock excitation and photoionization models for the excitation of the cloud. Shock models do not account for the [SiVI] and [CaVIII] emission line fluxes. However, X-ray observations indicate a source of ionizing photons in the core of Centaurus A; photoionization by the inferred flux incident on the cloud can account for the fluxes in these lines relative to Brackett-gamma. The power-law slope of the ionizing continuum matches that inferred from synchrotron models of the X-rays. The logarithm of the ionization parameter is -1.9, typical of that in Seyfert galaxies and consistent with the value proposed for dusty ionized plasmas. The model cloud density depends upon the Lorentz factor of the blazar and the inclination of our line of sight to the jet axis. For acute inclinations, the inferred density is consistent with expected cloud densities. However, for moderate inclinations of the jet to the line of sight, high Lorentz factors imply cloud densities in excess of $10^5 cm^{-3}$ and very low filling factors, suggesting that models of the gamma ray emission should incorporate jet Lorentz factors $\lesssim 5$.

The measurement of electron temperatures and metallicities in H ii regions and Planetary Nebulae (PNe) has-for several decades-presented a problem: results obtained using different techniques disagree. What it worse, they disagree consistently. There have been numerous attempts to explain these discrepancies, but none has provided a satisfactory solution to the problem. In this paper, we explore the possibility that electrons in H ii regions and PNe depart from a Maxwell-Boltzmann equilibrium energy distribution. We adopt a "kappa-distribution" for the electron energies. Such distributions are widely found in Solar System plasmas, where they can be directly measured. This simple assumption is able to explain the temperature and metallicity discrepancies in H ii regions and PNe arising from the different measurement techniques. We find that the energy distribution does not need to depart dramatically from an equilibrium distribution. From an examination of data from Hii regions and PNe it appears that kappa ~ 10 is sufficient to encompass nearly all objects. We argue that the kappa-distribution offers an important new insight into the physics of gaseous nebulae, both in the Milky Way and elsewhere, and one that promises significantly more accurate estimates of temperature and metallicity in these regions.

Recent observations show that the measured rates of star formation in the early universe are insufficient to produce re-ionization, and therefore, another source of ionizing photons is required. In this \emph{Letter}, we examine the possibility that these can be supplied by the fast accretion shocks formed around the cores of the most massive haloes ($10.5< \log M/M_{\odot} < 12$) on spatial scales of order 1 kpc. We model the detailed physics of these fast accretion shocks, and apply these to a simple 1-D spherical hydrodynamic accretion model for baryonic infall in dark matter halos with an Einasto density distribution. The escape of UV photons from these halos is delayed by the time taken to reach the critical accretion shock velocity for escape of UV photons; 220 km s$^{-1}$, and by the time it takes for these photons to ionize the surrounding baryonic matter in the accretion flow. Assuming that in the universe at large the baryonic matter tracks the dark matter, we can estimate the epoch of re-ionization in the case that accretion shocks act alone as the source of UV photons. We find that 50% of the volume (and 5-8% of the mass) of the universe can be ionized by $z \sim 7-8$. The UV production rate has an uncertainty of a factor of about 5 due to uncertainties in the cosmological parameters controlling the development of large scale structure. Because our mechanism is a steeply rising function of decreasing redshift, this uncertainty translates to a re-ionization redshift uncertainty of less than $\pm0.5$. We also find that, even without including the UV photon production of stars, re-ionization is essentially complete by $z \sim 5.8$. Thus, fast accretion shocks can provide an important additional source of ionizing photons in the early universe.

A new hydrodynamics code aimed at astrophysical applications has been developed. The new code and algorithms are presented along with a comprehensive suite of test problems in one, two, and three dimensions. The new code is shown to be robust and accurate, equalling or improving upon a set of comparison codes. Fyris Alpha will be made freely available to the scientific community.

We have performed a series of three-dimensional simulations of the interaction of a supersonic wind with a non-spherical radiative cloud. These simulations are motivated by our recent three-dimensional model of a starburst-driven galactic wind interacting with an inhomogeneous disk, which show that an optically emitting filament can be formed by the break-up and acceleration of a cloud into a supersonic wind. In this study we consider the evolution of a cloud with two different geometries (fractal and spherical) and investigate the importance of radiative cooling on the cloud's survival. We have also undertaken a comprehensive resolution study in order to ascertain the effect of the assumed numerical resolution on the results. We find that the ability of the cloud to radiate heat is crucial for its survival. While an adiabatic cloud is destroyed over a short period of time, a radiative cloud is broken up via the Kelvin-Helmholtz instability into numerous small, dense cloudlets, which are drawn into the flow to form a filamentary structure. The degree of fragmentation is highly dependent on the resolution of the simulation, with the number of cloudlets formed increasing as the Kelvin-Helmholtz instability is better resolved. Nevertheless, there is a clear qualitative trend, with the filamentary structure still persistent at high resolution. We confirm the mechanism behind the formation of the H-alpha emitting filaments found in our global simulations of a starburst-driven wind. Based on our resolution study, we conclude that bow shocks around accelerated gas clouds, and their interaction, are the main source of the soft X-ray emission observed in these galactic-scale winds. [ABRIDGED]

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.

We combine the the stellar spectral synthesis code Starburst99, the nebular modelling code MAPPINGSIII, and a 1-D dynamical evolution model of HII regions around massive clusters of young stars to generate improved models of the spectral energy distribution (SED) of starburst galaxies. We introduce a compactness parameter, C, which characterizes the specific intensity of the radiation field at ionization fronts in HII regions, and which controls the shape of the far-IR dust re-emission, often referred to loosely as the dust ``temperature''. We also investigate the effect of metallicity on the overall SED and in particular, on the strength of the PAH features. We provide templates for the mean emission produced by the young compact HII regions, the older (10 - 100 Myr) stars and for the wavelength-dependent attenuation produced by a foreground screen of the dust used in our model. We demonstrate that these components may be combined to produce a excellent fit to the observed SEDs of star formation dominated galaxies which are often used as templates (Arp 220 and NGC 6240). This fit extends from the Lyman Limit to wavelengths of about one mm. The methods presented in both this paper and in the previous papers of this series allow the extraction of the physical parameters of the starburst region (star formation rates, star formation rate history, mean cluster mass, metallicity, dust attenuation and pressure) from the analysis of the pan-spectral SED.

Since its discovery in 1996, the source of the bright H-alpha emission (up to 750 mR) along the Magellanic Stream has remained a mystery. There is no evidence of ionising stars within the HI stream, and the extended hot halo is far too tenuous to drive strong shocks into the clouds. We now present a hydrodynamical model that explains the known properties of the H-alpha emission and provides new insights on the lifetime of the Stream clouds. The upstream clouds are gradually disrupted due to their interaction with the hot halo gas. The clouds that follow plough into gas ablated from the upstream clouds, leading to shock ionisation at the leading edges of the downstream clouds. Since the following clouds also experience ablation, and weaker H-alpha (100-200 mR) is quite extensive, a disruptive cascade must be operating along much of the Stream. In our model, the clouds are evolving on timescales of 100-200 Myr, such that the Stream must be replenished by the Magellanic Clouds at a fairly constant rate. The ablated material falls onto the Galaxy as a warm drizzle which suggests that diffuse ionized gas at 10**4 K may be an important constituent of galactic accretion. The observed HI emission provides a new constraint on the rate of disruption of the Stream and, consequently, the infall rate of metal-poor gas onto the Galaxy. When the ionized component of the Stream is fully accounted for, the rate of gas accretion is 0.4 Msun/yr, roughly twice the rate deduced from HI observations alone.

We have performed a series of three-dimensional simulations of a starburst-driven wind in an inhomogeneous interstellar medium. The introduction of an inhomogeneous disk leads to differences in the formation of a wind, most noticeably the absence of the ``blow-out'' effect seen in homogeneous models. A wind forms from a series of small bubbles that propagate into the tenuous gas between dense clouds in the disk. These bubbles merge and follow the path of least resistance out of the disk, before flowing freely into the halo. Filaments are formed from disk gas that is broken up and accelerated into the outflow. These filaments are distributed throughout a biconical structure within a more spherically distributed hot wind. The distribution of the inhomogeneous interstellar medium in the disk is important in determining the morphology of this wind, as well as the distribution of the filaments. While higher resolution simulations are required in order to ascertain the importance of mixing processes, we find that soft X-ray emission arises from gas that has been mass-loaded from clouds in the disk, as well as from bow shocks upstream of clouds, driven into the flow by the ram pressure of the wind, and the interaction between these shocks.

We present three dimensional simulations of the interaction of a light hypersonic jet with an inhomogeneous thermal and turbulently supported disk in an elliptical galaxy, including Radio and multi-band X-ray visualisations. These simulations are applicable to the GPS/CSS phase of some extragalactic radio sources. We identify four generic phases in the evolution of such a jet with the interstellar medium. The first is a `flood and channel'' phase, dominated by complex jet interactions with the dense cloudy medium close to the nucleus. A spherical, energy driven, bubble phase follows, where the bubble is larger than the disk scale, but the jet remains fully disrupted close to the nucleus. Then in a rapid, jet break--out phase, the jet breaks free of the last obstructing dense clouds, becomes collimated and pierces the more or less spherical bubble. In the final classical phase, the jet propagates in a momentum-dominated fashion similar to jets in single component hot haloes, leading to the classical jet -- cocoon -- bow-shock structure.

We present three dimensional simulations of the interaction of a light hypersonic jet with an inhomogeneous thermal and turbulently supported disk in an elliptical galaxy. We model the jet as a light, supersonic non-relativistic flow with parameters selected to be consistent with a relativistic jet with kinetic power just above the FR1/FR2 break. We identify four generic phases in the evolution of such a jet with the inhomogeneous interstellar medium: 1) an initial ``flood and channel'' phase, where progress is characterized by high pressure gas finding changing weak points in the ISM, flowing through channels that form and re-form over time, 2) a spherical, energy-driven bubble phase, were the bubble is larger than the disk scale, but the jet remains fully disrupted close to the nucleus, 3) a rapid, jet break--out phase the where jet breaks free of the last dense clouds, becomes collimated and pierces the spherical bubble, and 4) a classical phase, the jet propagates in a momentum-dominated fashion leading to the classical jet + cocoon + bow-shock structure. Mass transport in the simulations is investigated, and we propose a model for the morphology and component proper motions in the well-studied Compact Symmetric Object 4C31.04.

We present first results from the Near-infrared Integral Field Spectrograph (NIFS) located at Gemini North. For the active galaxies Cygnus A and Perseus A we observe rotationally-supported accretion disks and adduce the existence of massive central black holes and estimate their masses. In Cygnus A we also see remarkable high-excitation ionization cones dominated by photoionization from the central engine. In the T-Tauri stars HV Tau C and DG Tau we see highly-collimated bipolar outflows in the [Fe II] 1.644 micron line, surrounded by a slower molecular bipolar outflow seen in the H_2 lines, in accordance with the model advocated by Pyo et al. (2002).

We build, as far as theory will permit, self consistent model HII regions around central clusters of aging stars. These produce strong emission line diagnostics applicable to either individual HII regions in galaxies, or to the integrated emission line spectra of disk or starburst galaxies. The models assume that the expansion and internal pressure of individual HII regions is driven by the net input of mechanical energy from the central cluster, be it through winds or supernova events. This eliminates the ionization parameter as a free variable, replacing it with a parameter which depends on the ratio of the cluster mass to the pressure in the surrounding interstellar medium. These models explain why HII regions with low abundances have high excitation, and demonstrate that at least part of the warm ionized medium is the result of overlapping faint, old, large, and low pressure HII regions. We present line ratios (at both optical and IR wavelengths) which provide reliable abundance diagnostics for both single HII regions or for integrated galaxy spectra, and we find a number that can be used to estimate the mean age of the cluster stars exciting individual HII regions.

We present models for the mid- and far- infrared emission from the Narrow Line Region (NLR) of Active Galactic Nuclei (AGN). Using the MAPPINGS III code we explore the effect of typical NLR parameters on the spectral characteristics of the IR emission. These include useful IR emission line ratio diagnostic diagrams for the determination of these parameters, as well as Star formation--AGN mixing diagnostics. We also examine emission line to continuum correlations which would assist in separating the IR emission arising from the NLR from that coming from the inner torus. We find for AGN like NGC 1068 and NGC 4151 that the NLR only contributes ~10% to the total IRAS 25 mum flux, and that other components such as a dusty torus are necessary to explain the total AGN IR emission.

We examine, from a theoretical viewpoint, how the physical parameters of HII regions are controlled in both normal galaxies and in starburst environments. These parameters are the HII region luminosity function, the time-dependent size, the covering fraction of molecular clouds, the pressure in the ionized gas and the ionization parameter. The factors which control them are the initial mass function of the exciting stars, the cluster mass function, the metallicity and the mean pressure in the surrounding interstellar medium. We investigate the sensitivity of the H$\alpha$ luminosity to the IMF, and find that this can translate to more than a factor two variation in derived star formation rates. The molecular cloud dissipation timescale is estimated from a case study of M17 to be $\sim1$ Myr for this object. Based upon HII luminosity function fitting for nearby galaxies, we suggest that the \HII region cluster mass function is fitted by a log-normal form peaking at $\sim 100 M_{\odot}$. The cluster mass function continues the stellar IMF to higher mass. The pressure in the HII regions is controlled by the mechanical luminosity flux from the central cluster. Since this is closely related to the ionizing photon flux, we show that the ionization parameter is not a free variable, and that the diffuse ionized medium may be composed of many large, faint and old HII regions. Finally, we derive theoretical probability distributions for the ionization parameter as a function of metallicity and compare these to those derived for SDSS galaxies.

We present two-dimensional slab-jet simulations of jets in inhomogeneous media consisting of a tenuous hot medium populated with a small filling factor by warm, dense clouds. The simulations are relevant to the structure and dynamics of sources such as Gigahertz Peak Spectrum and Compact Steep Spectrum radio galaxies, High Redshift Radio Galaxies and radio galaxies in cooling flows. The jets are disrupted to a degree depending upon the filling factor of the clouds. With a small filling factor, the jet retains some forward momentum but also forms a halo or bubble around the source. At larger filling factors channels are formed in the cloud distribution through which the jet plasma flows and a hierarchical structure consisting of nested lobes and an outer enclosing bubble results. We suggest that the CSS quasar 3C48 is an example of a low filling factor jet - interstellar medium interaction whilst M87 may be an example of the higher filling factor type of interaction. Jet disruption occurs primarily as a result of Kelvin-Helmholtz instabilities driven by turbulence in the radio cocoon not through direct jet-cloud interactions, although there are some examples of these. In all radio galaxies whose morphology may be the result of jet interactions with an inhomogeneous interstellar medium we expect that the dense clouds will be optically observable as a result of radiative shocks driven by the pressure of the radio cocoon. We also expect that the radio galaxies will possess faint haloes of radio emitting material well beyond the observable jet structure.

Galaxies form out of small fluctuations in a smoothly expanding Universe. However, the initial gravitational collapse phase is accompanied by the formation of supermassive black holes and clusters of massive stars. Black holes and star clusters generate powerful outflows in the form of jets and superwinds that interact with still infalling gas, possibly regulating the galaxy formation process, initiating new sites of star formation, and carrying chemically enriched gas to the intergalactic medium. Unfortunately, beyond this qualitative description our detailed theoretical understanding is poor. New results from 3D simulations of a GPS/CSS galaxy, with gravitational potentials included, shed some new light on the jet driven outflow process in particular. New code capabilites to predict detailed X-ray spectra from multi-dimensional time-dependent dynamics simulations of Galaxy Feedback, and will be useful for future interpretation of X-ray and radio SEDs of forming galaxies.

In this paper, we combine the stellar spectral synthesis code STARBURST 99, the nebular modelling code MAPPINGS IIIq, a 1-D dynamical evolution model of \HII regions around massive clusters of young stars and a simplified model of synchrotron emissivity to produce purely theoretical self-consistent synthetic spectral energy distributions (SEDs) for (solar metallicity) starbursts lasting some $10^8$ years. These SEDs extend from the Lyman Limit to beyond 21 cm. We find that two ISM parameters control the form of the SED; the pressure in the diffuse phase of the ISM (or, equivalently, its density), and the molecular cloud dissipation timescale. We present detailed SED fits to Arp 220 and NGC 6240, and we give the predicted colors for starburst galaxies derived from our models for the IRAS and the Spitzer Space Observatory MIPS and IRAC instruments. Our models reproduce the spread in observed colors of starburst galaxies. Finally, we present absolute calibrations to convert observed fluxes into star formation rates in the UV (GALEX), at optical wavelengths (H$\alpha$), and in the IR (IRAS or the Spitzer Space Observatory). (Abstract Truncated)

Seyfert narrow line regions (NLR) have emission line ratios which are remarkably uniform, displaying only ~0.5 dex variation between different galaxies. Existing models have been unable to explain these observations without the introduction of ad hoc assumptions, geometrical restrictions or new parameters. Here we introduce a new model: dusty radiation pressure dominated photoionization, which provides a natural self-regulating characteristic leading to an invariance of the spectrum over a very wide range (>100) of ionization parameter. The dusty model is able to reproduce both the range and the absolute value of the observational line ratios not only in the standard optical diagnostic diagrams but also in UV diagnostic plots, providing an explanation to the problem in NLR observations.

Seyfert narrow line region (NLR) emission line ratios are remarkably uniform, displaying only ~0.5 dex variation between galaxies, and even less within an individual object. Previous photoionization and shock models of this region were unable to explain this observation without the introduction of arbitrary assumptions or additional parameters. Dusty, radiation pressure dominated photoionization models provide a simple physical mechanism which can reproduce this spectral uniformity between different objects. In the first paper of this series we described this model and its implementation in detail, as well as presenting grids of model emission lines and examining the model structures. Here we explore these models further, demonstrating their ability to reproduce the observed Seyfert line ratios on standard line diagnostic diagrams in both the optical and UV. We also investigate the effects that the variation of metallicity, density and ionizing spectrum have upon both the new paradigm and the standard photoionization models used hitherto. Along with the standard diagnostic diagrams we provide several new diagnostic diagrams in the UV, Optical and IR. These new diagrams can provide further tests of the dusty, radiation pressure photoionization paradigm as well as being used as diagnostics of the metallicity, density and ionizing spectrum of the emission line clouds.

We present the implementation of dusty, radiation pressure dominated photoionization models applicable to the Narrow Line Regions (NLRs) of Active Galactic Nuclei, using the \mapiii code. We give a grid of the predicted intensities of the most commonly used diagnostic spectral lines in the UV, optical and IR, covering a wide range of density, metallicity, the power-law index characterizing the photoionizing source and photoionization parameter, for use in the diagnosis of NLRs. We examine the temperature, density and ionization structure of these models, investigating the effect of variation of these parameters in order to gain a better understanding of NLR clouds themselves.

The empirical reddening function for starburst galaxies generated by Calzetti and her co-workers has proven very successful, and is now used widely in the observational literature. Despite its success, however, the physical basis for this extinction law, or more correctly, attenuation law remains weak. Here we provide a physical explanation for the Calzetti Law based on a turbulent interstellar medium. In essence, this provides a log-normal distribution of column densities, giving a wide range of column densities in the dusty foreground screen. Therefore, extended sources such as starburst regions or HII regions seen through it suffer a point-to-point stochastic extinction and reddening. Regions of high column densities are "black" in the UV, but translucent in the IR, which leads to a flatter extinction law, and a larger value of the total to selective extinction, R_V. We fit the Calzetti Law, and infer that the variance sigma of the log-normal distribution lies in the range 0.6<sigma<2.2. The absolute to selective extinction R_V is found to be in the range 4.3 to 5.2 consistent with R_V=4.05+/-0.80 of the Calzetti Law.

Star formation rates in galaxies are frequently estimated using the Balmer line fluxes. However, these can be systematically underestimated because dust competes for the absorption of Lyman continuum photons in the ionized gas. Here we present theoretical correction factors in a simple analytic form. T These factors scale as the product of the ionization parameter, ${\cal U}$, and the nebular O/H abundance ratio, both of which can now be derived from the observation of bright nebular line ratios. The correction factors are only somewhat dependent upon the photoelectron production by grains, but are very sensitive to the presence of complex PAH-like carbonaceous molecules in the ionized gas, providing that these can survive in such an environment.

We review recent progress by the RSAA team in elucidating the dynamical evolution of the various classes of Narrow Line Regions (NLR) in active galaxies.

High redshift radio galaxies are great cosmological tools for pinpointing the most massive objects in the early Universe: massive forming galaxies, active super--massive black holes and proto--clusters. We report on deep narrow--band imaging and spectroscopic observations of several z > 2 radio galaxy fields to investigate the nature of giant Ly-alpha nebulae centered on the galaxies and to search for over--dense regions around them. We discuss the possible implications for our understanding of the formation and evolution of massive galaxies and galaxy clusters.

The radio lobes of the radio galaxy Hercules A contain intriguing ring-like structures concentric with the jet axis. To investigate the occurrence of such features, we have used hydrodynamic simulations of jets with a range of Mach numbers (from M=2 to 50) and densities (down to a ratio of 1E-4 relative to the background) to generate ray-traced images simulating synchrotron emission from the time-dependent shock structures. We compare these images with observations of Hercules A, and consider the physical nature and temporal evolution of the most plausible configurations. We find that the observed ring-like structures are well explained as nearly annular shocks propagating in the backflow surrounding the jet. We infer that the jet is oriented at between 30 and 70 degrees to the line of sight, consistent with radio depolarisation observations of Gizani & Leahy. The observational lack of hot-spots at the extremities of the radio lobes, and the possible presence of a buried hot-spot near the base of the western lobe, are explained in terms of the intrinsic brightness fluctuations and dynamics of the terminal shock of an ultra-light, low Mach number jet that surges along its axis due to intermittent pinching and obstruction by turbulent backflow in the cocoon. We conclude from the appearance of both sides of the Hercules A, that both jets are on the borderline of becoming fully turbulent.

We address a numerical instability that arises in the directionally split computation of hydrodynamic flows when shock fronts are parallel to a grid plane. Transverse oscillations in pressure, density and temperature are produced that are exacerbated by thermal instability when cooling is present, forming post--shock `stripes'. These are orthogonal to the classic post--shock 'ringing' fluctuations. The resulting post--shock `striping' substantially modifies the flow. We discuss three different methods to resolve this problem. These include (1) a method based on artificial viscosity; (2) grid--jittering and (3) a new localized oscillation filter that acts on specific grid cells in the shock front. These methods are tested using a radiative wall shock problem with an embedded shear layer. The artificial viscosity method is unsatisfactory since, while it does reduce post--shock ringing, it does not eliminate the stripes and the excessive shock broadening renders the calculation of cooling inaccurate, resulting in an incorrect shock location. Grid--jittering effectively counteracts striping. However, elsewhere on the grid, the shear layer is unphysically diffused and this is highlighted in an extreme case. The oscillation filter method removes stripes and permits other high velocity gradient regions of the flow to evolve in a physically acceptable manner. It also has the advantage of only acting on a small fraction of the cells in a two or three dimensional simulation and does not significantly impair performance.

We have carried out simulations of supersonic light jets in order to model the features observed in optical and radio images of the western hot-spot in the radio galaxy Pictor A. We have considered jets with density ratios eta=1e-2 - 1e-4, and Mach numbers ranging between 5 and 50. From each simulation, we have generated ray-traced maps of radio surface brightness at a variety of jet inclinations, in order to study the appearance of time-dependent luminous structures in the vicinity of the western hot-spot. We compare these rendered images with observed features of Pictor A. A remarkable feature of Pictor A observations is a bar-shaped "filament" inclined almost at right angles to the inferred jet direction and extending 24" (10.8 /h kpc) along its longest axis. The constraints of reproducing the appearance of this structure in simulations indicate that the jet of Pictor A lies nearly in the plane of the sky. The results of the simulation are also consistent with other features found in the radio image of Pictor A. This filament arises from the surging behaviour of the jet near the hot-spot; the surging is provoked by alternate compression and decompression of the jet by the turbulent backflow in the cocoon. We also examine the arguments for the jet in Pictor A being at a more acute angle to the line of sight and find that our preferred orientation is just consistent with the limits on the brightness ratio of the X-ray jet and counter-jet. We determine from our simulations, the structure function of hot-spot brightness and also the cumulative distribution of the ratio of intrinsic hot-spot brightnesses. The latter may be used to quantify the use of hot-spot ratios for the estimation of relativistic effects.

The remarkable similarity between emission spectra of narrow line regions (NLR) in Seyfert Galaxies has long presented a mystery. In photoionization models, this similarity implies that the ionization parameter is nearly always the same, about U ~ 0.01. Here we present dusty, radiation-pressure dominated photoionization models that can provide natural physical insight into this problem. In these models, dust and the radiation pressure acting on it provide the controlling factor in moderating the density, excitation and surface brightness of photoionized NLR structures. Additionally, photoelectric heating by the dust is important in determining the temperature structure of the models. These models can also explain the coexistence of the low-, intermediate- and coronal ionization zones within a single self-consistent physical structure. The radiation pressure acting on dust may also be capable of driving the fast (~3000 km/s) outflows such as are seen in the HST observations of NGC 1068.

We model the northern middle radio lobe of Centaurus A (NGC 5128) as a buoyant bubble of plasma deposited by an intermittently active jet. The extent of the rise of the bubble and its morphology imply that the ratio of its density to that of the surrounding ISM is less than 10^{-2}, consistent with our knowledge of extragalactic jets and minimal entrainment into the precursor radio lobe. Using the morphology of the lobe to date the beginning of its rise through the atmosphere of Centaurus A, we conclude that the bubble has been rising for approximately 140Myr. This time scale is consistent with that proposed by Quillen et al. (1993) for the settling of post-merger gas into the presently observed large scale disk in NGC 5128, suggesting a strong connection between the delayed re-establishment of radio emission and the merger of NGC 5128 with a small gas-rich galaxy. This suggests a connection, for radio galaxies in general, between mergers and the delayed onset of radio emission. In our model, the elongated X-ray emission region discovered by Feigelson et al. (1981), part of which coincides with the northern middle lobe, is thermal gas that originates from the ISM below the bubble and that has been uplifted and compressed. The "large-scale jet" appearing in the radio images of Morganti et al. (1999) may be the result of the same pressure gradients that cause the uplift of the thermal gas, acting on much lighter plasma, or may represent a jet that did not turn off completely when the northern middle lobe started to buoyantly rise. We propose that the adjacent emission line knots (the "outer filaments") and star-forming regions result from the disturbance, in particular the thermal trunk, caused by the bubble moving through the extended atmosphere of NGC 5128.

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 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 which 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 field in these galaxies. These diagnostic diagrams are most sensitive to the spectral index of the ionizing radiation field in the 1-4 Rydberg region. We find that warm infrared starburst galaxies contain a relatively hard EUV field in this region. The PEGASE ionizing stellar continuum is harder in the 1-4 Rydberg range than that of STARBURST99. As the spectrum in this regime is dominated by emission from Wolf-Rayet (W-R) stars, this difference is most likely due to the differences in stellar atmosphere models used for the W-R stars. We believe that the stellar atmospheres in STARBURST99 are more applicable to the starburst galaxies in our sample, however they do not produce the hard EUV field in the 1-4 Rydberg 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 classification scheme for starbursts and AGN galaxies based on the optical diagnostic diagrams.

We present Hubble Space Telescope (HST) and ground-based observations of a pair of galaxies at redshift 2.38, which are collectively known as 2142-4420 B1 (Francis et al. 1996). The two galaxies are both luminous extremely red objects (EROs), separated by 0.8 arcsec. They are embedded within a 100 kpc scale diffuse Ly-alpha nebula (or blob) of luminosity ~10^44 erg/s. The radial profiles and colors of both red objects are most naturally explained if they are young elliptical galaxies: the most distant yet found. It is not, however, possible to rule out a model in which they are abnormally compact, extremely dusty starbursting disk galaxies. If they are elliptical galaxies, their stellar populations have inferred masses of ~10^11 solar masses and ages of ~7x10^8 years. Both galaxies have color gradients: their centers are significantly bluer than their outer regions. The surface brightness of both galaxies is roughly an order of magnitude greater than would be predicted by the Kormendy relation. A chain of diffuse star formation extending 1 arcsec from the galaxies may be evidence that they are interacting or merging. The Ly-alpha nebula surrounding the galaxies shows apparent velocity substructure of amplitude ~ 700 km/s. We propose that the Ly-alpha emission from this nebula may be produced by fast shocks, powered either by a galactic superwind or by the release of gravitational potential energy.

The high redshift radio galaxy 4C41.17 consists of a powerful radio source in which previous work has shown that there is strong evidence for jet-induced star formation along the radio axis. We argue that nuclear photoionization is not responsible for the excitation of the emission line clouds and we construct a jet-cloud interaction model to explain the major features revealed by the data. The interaction of a high-powered jet with a dense cloud in the halo of 4C41.17 produces shock-excited emission-line nebulosity through ~1000 km/s shocks and induces star formation. The CIII to CIV line ratio and the CIV luminosity emanating from the shock, imply that the pre-shock density in the line-emitting cloud is high enough (~1-10 cm^-3) that shock initiated star formation could proceed on a timescale of order a few x 10^6 yrs, well within the estimated dynamical age of the radio source. Broad (FWHM ~ 100 - 1400 km/s) emission lines are attributed to the disturbance of the gas cloud by a partial bow--shock and narrow emission lines (FWHM ~ 500 - 650 km/s) (in particular CIV) arise in precursor emission in relatively low metallicity gas. The implied baryonic mass ~ 8 \times 10^{10} solar masses of the cloud is high and implies that Milky Way size condensations existed in the environments of forming radio galaxies at a redshift of 3.8. Our interpretation of the data provides a physical basis for the alignment of the radio, emission-line and UV continuum images in some of the highest redshift radio galaxies and the analysis presented here may form a basis for the calculation of densities and cloud masses in other high redshift radio galaxies.

ISO-SWS spectroscopy of the WR galaxy NGC5253 is presented, and analysed to provide estimates of its hot young star population. Our approach differs from previous investigations in that we are able to distinguish between the regions in which different infrared fine-structure lines form, using complementary ground-based observations. The high excitation nebular [SIV] emission is formed in a very compact region, which we attribute to the central super-star-nucleus, and lower excitation [NeII] nebular emission originates in the galactic core. We use photo-ionization modelling coupled with the latest theoretical O-star flux distributions to derive effective stellar temperatures and ionization parameters of Teff>38kK, logQ=8.25 for the compact nucleus, with Teff=35kK, logQ<8 for the larger core. Results are supported by more sophisticated calculations using evolutionary synthesis models. We assess the contribution that Wolf-Rayet stars may make to highly ionized nebular lines (e.g. [OIV]). From our Br(alpha) flux, the 2" nucleus contains the equivalent of approximately 1000 O7V star equivalents and the starburst there is 2-3Myr old; the 20" core contains about 2500 O7V star equivalents, with a representative age of 5Myr. The Lyman ionizing flux of the nucleus is equivalent to the 30 Doradus region. These quantities are in good agreement with the observed mid-IR dust luminosity of 7.8x10^8 L(sun) Since this structure of hot clusters embedded in cooler emission may be common in dwarf starbursts, observing a galaxy solely with a large aperture may result in confusion. Neglecting the spatial distribution of nebular emission in NGC5253, implies `global' stellar temperatures (or ages) of 36kK (4.8Myr) and 39kK (2.9 or 4.4Myr) from the observed [NeIII/II] and [SIV/III] line ratios, assuming logQ=8.

We present near-infrared spectra for a sample of galaxies with ambiguous optical emission line ratios. These galaxies fall between starbursts and Seyferts in the usual optical diagnostic diagrams. We find a similar result with the near-infrared emission line ratios, which suggests that the galaxies are composite, containing both a starburst and AGN component. Furthermore, CO absorption, produced in late-type stars, is detected within the sample, although at a weaker level than is typical for starburst galaxies. We conclude that the CO feature is being diluted by a contribution from an AGN, thereby confirming the composite nature of these galaxies.

This paper presents long slit spectral maps of the bi-cone shaped extended narrow line region (ENLR) in the Seyfert galaxy NGC 2992. We investigate the physical properties of the ENLR via emission line diagnostics, and compare the observations to shock and photoionization models for the excitation mechanism of the gas. The line ratios vary as a function of position in the ENLR, and the loci of the observed points on line ratio diagrams are shown to be most consistent with shock+precursor model grids. We consider the energetics of a nuclear ionizing source for the ENLR, and perform the q-test in which the rate of ionizing photons from the nucleus is inferred from measurements of the density and ionization parameter. The q-test is shown to be invalid in the case of NGC 2992 because of the limitations of the [S II]6717/6731 density diagnostic. The excitation of the gas is shown to be broadly consistent with the kinematics, with higher [N II]6583/H-alpha present in the more dynamically active region. We also show that the pressure associated with the X-ray emitting plasma may provide a large fraction of the pressure required to power the ENLR via shocks.

We argue that the narrow line regions of Seyfert galaxies are powered by the transport of energy and momentum by the radio-emitting jets and consequently that the ratio of the radio power to jet energy flux is much smaller than is usually assumed for radio galaxies. This can be partially attributed to the smaller ages ($\sim 10^6 yrs$) of Seyferts compared to radio galaxies but one also requires that either the magnetic energy density is more than an order of magnitude below the equipartition value, or more likely, that the internal energy densities of Seyfert jets are dominated by thermal plasma. If Seyfert jets are initially dominated by relativistic plasma, then an analysis of the data on jets in five Seyfert galaxies shows that all but one of these would have mildly relativistic jet velocities near 100 pc in order to power the respective narrow-line regions. However, observations of jet-cloud interactions in the NLR provide additional information on jet velocities and composition via the momentum budget. Our analysis of a jet-cloud interaction in NGC 1068, implies a shocked jet pressure much larger than the minimum pressure of the radio knot, a velocity $\sim 0.06 c$ and a jet temperature $\sim 10^9 K$ implying mildly relativistic electrons but thermal protons. The jet mass flux at this point $\sim 0.5 M_\odot yr^{-1}$, is an order of magnitude higher than the mass accretion rate into the black hole, strongly indicating entrainment. The initial jet mass flux $\sim 0.02 M_\odot yr^{-1}$, comparable to the mass accretion rate and is consistent with the densities inferred for accretion disc coronae from high energy observations, together with an initially mildly relativistic velocity and an initial jet radius of order 10 gravitational radii.