BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Ralph S. Sutherland

Ralph S. Sutherland contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.GAastro-ph.COastro-ph.HEastro-ph.SRphysics.comp-phastro-phastro-ph.IMhep-phphysics.flu-dyn

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
18works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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 &#39;d&#39; is arguably the most important tracer of what happens to gas accreting onto a disk galaxy. Recent observations reveal that the Stream&#39;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&#39;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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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&#39; models) and collisions between non-thermal protons and gas (`hadronic&#39; models) have been advanced as the process supplying the Bubbles&#39; γ-ray emission. An issue for any steady state hadronic model is that the very low density of the Bubbles&#39; 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&#39; hadronic γ-ray emission is $\sim$ few $\times 10^8$ years. Our model invokes collapse of the Bubbles&#39; 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&#39; plasma, $M_{pls}$).

0 citations0 reviewsNo code linkedOpen access
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 &#39;κ&#39; 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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$.

0 citations0 reviewsNo code linkedOpen access
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 &#34;kappa-distribution&#34; 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access