Source author record

Ralph S. Sutherland

Ralph S. Sutherland appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2022arXiv

Theoretically Modelling Photoionized Regions with Fractal Geometry in Three Dimensions

We create a photoionization model embedded in the turbulent ISM by using the state-of-the-art Messenger Monte-Carlo MAPPINGS~V code (M$^3$) in conjunction with the CMFGEN stellar atmosphere model. We show that the turbulent ISM causes the inhomogeneity of electron temperature and density within the nebula. The fluctuation in the turbulent ISM creates complex ionization structures seen in nearby nebulae. The inhomogeneous density distribution within the nebula creates a significant scatter on the spatially-resolved standard optical diagnostic diagrams, which cannot be represented by the spherical constant density photoionization model. We analyze the dependence of different optical emission lines on the complexity of nebular geometry, finding that the emission-lines residing on the nebular boundary are highly sensitive to the complexity of nebular geometry, while the emission-lines produced throughout the nebula are sensitive to the density distribution of the ISM within the nebula. Our fractal photoionization model demonstrates that a complex nebular geometry is required for accurate modeling of HII regions and emission-line galaxies, especially for the high-redshift galaxies, where the ISM is highly turbulent based on the increasing observational evidence.

preprint2021arXiv

On the origin of nitrogen at low metallicity

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.

preprint2016arXiv

A physically-based model of the ionizing radiation from active galaxies for photoionization modeling

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 $Γ$, 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.

preprint2016arXiv

Chemical abundances in high-redshift galaxies: A powerful new emission line diagnostic

This Letter presents a new, remarkably simple diagnostic specifically designed to derive chemical abundances for high redshift galaxies. It uses only the H α, [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.

preprint2016arXiv

Forbidden Iron Lines and Dust Destruction in Supernova Remnant Shocks: The Case of N49 in the Large Magellanic Cloud

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.

preprint2016arXiv

Relativistic jet feedback in high-redshift galaxies I: Dynamics

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.

preprint2016arXiv

The Discovery of Raman Scattering in HII Regions

We report here on the discovery of faint extended wings of Hα 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β and in Hγ 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Å resonance line involved in the Raman scattering by the Lyβ radiation field. The subsequent radiative cascade produces enhanced Si II 5978.9, 6347.1$ and 6371.4Å permitted transitions. Finally we show that in O I, radiative pumping of the 1025.76Å 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α 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

preprint2015arXiv

Star Formation in Ultrafaint Dwarfs: Continuous or Single-age Bursts?

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 [$α$/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.

preprint2015arXiv

The Chemical Evolution of Very Metal-Poor Damped Lyman-$α$ Systems

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-$α$ 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 [$α$/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 [$α$/Fe] in DLAs compared to UFDs results from enrichment of second-burst systems by Type Ia supernovae.

preprint2015arXiv

The density variance - Mach number relation in isothermal and non-isothermal adiabatic turbulence

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.

preprint2015arXiv

The Magellanic Stream: break up and accretion onto the hot Galactic corona

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.

preprint2015arXiv

Turbulent mixing layers in supersonic protostellar outflows, with application to DG Tauri

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 μ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.

preprint2014arXiv

Galaxy emission line classification using 3D line ratio diagrams

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$α$ vs. log[OIII]/H$β$ 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.

preprint2014arXiv

Metal-poor dwarf galaxies in the SIGRID galaxy sample. I. HII region observations and chemical abundances

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] λλ 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.

preprint2014arXiv

Metal-poor dwarf galaxies in the SIGRID galaxy sample. II. The electron temperature-abundance calibration and the parameters that affect it

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.

preprint2014arXiv

Steady-State Hadronic Gamma-Ray Emission from 100-Myr-Old Fermi Bubbles

The Fermi Bubbles are enigmatic γ-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' γ-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 γ-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_γ\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}$).

preprint2013arXiv

Measuring nebular temperatures: the effect of new collision strengths with equilibrium and kappa-distributed electron energies

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 'κ' 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 κ 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 κ 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.

preprint2013arXiv

New Strong Line Abundance Diagnostics for \HII Regions: Effects of $κ$-Distributed Electron Energies and New Atomic Data

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 $κ$-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 $κ$-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 $ κ\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 $κ$-distribution on the strong line abundances derived solely on the basis of theoretical models is rather small.

preprint2013arXiv

The Kinematics and Ionization of Nuclear Gas Clouds in Centaurus A

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$.

preprint2012arXiv

Resolving the electron temperature discrepancies in HII Regions and Planetary Nebulae: kappa-distributed electrons

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.

preprint2011arXiv

Re-ionizing the Universe without Stars

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.

preprint2006arXiv

Modelling the Pan-Spectral Energy Distribution of Starburst Galaxies: III. Emission Line Diagnostics of Ensembles of Evolving HII Regions

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.