Source author record

R. S. Sutherland

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

ResearcherUnclaimed source record

astro-ph.GA astro-ph.HE astro-ph astro-ph.CO astro-ph.SR

Catalog footprint

What is connected

8works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2021arXiv

Stromlo Stellar Tracks: non-solar scaled abundances for massive stars

We present the Stromlo Stellar Tracks, a set of stellar evolutionary tracks, computed by modifying the Modules for Experiments in Stellar Astrophysics (MESA) 1D stellar evolution package, to fit the Galactic Concordance abundances for hot ($T > 8000$ K) massive ($\geq 10M_\odot$) Main-Sequence (MS) stars. Until now, all stellar evolution tracks are computed at solar, scaled-solar, or alpha-element enhanced abundances, and none of these models correctly represent the Galactic Concordance abundances at different metallicities. This paper is the first implementation of Galactic Concordance abundances to the stellar evolution models. The Stromlo tracks cover massive stars ($10\leq M/M_\odot \leq 300$) with varying rotations ($v/v_{\rm crit} = 0.0, 0.2, 0.4$) and a finely sampled grid of metallicities ($-2.0 \leq {\rm [Z/H]} \leq +0.5$; $Δ{\rm [Z/H]} = 0.1$) evolved from the pre-main sequence to the end of $^{12}$Carbon burning. We find that the implementation of Galactic Concordance abundances is critical for the evolution of main-sequence, massive hot stars in order to estimate accurate stellar outputs (L, T$_{\rm eff}$, $g$), which, in turn, have a significant impact on determining the ionizing photon luminosity budgets. We additionally support prior findings of the importance that rotation plays on the evolution of massive stars and their ionizing budget. The evolutionary tracks for our Galactic Concordance abundance scaling provide a more empirically motivated approach than simple uniform abundance scaling with metallicity for the analysis of HII regions and have considerable implications in determining nebular emission lines and metallicity. Therefore, it is important to refine the existing stellar evolutionary models for comprehensive high-redshift extragalactic studies. The Stromlo tracks are publicly available to the astronomical community online.

preprint2016arXiv

Jet-Intracluster Medium interaction in Hydra A. II The Effect of Jet Precession

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.

preprint2015arXiv

Galaxy-scale AGN Feedback - Theory

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.

preprint2014arXiv

Jet-Intracluster Medium interaction in Hydra A. I Estimates of jet velocity from inner knots

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.

preprint2014arXiv

Multi-dimensional simulations of the expanding supernova remnant of SN 1987A

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.

preprint2009arXiv

Starburst-Driven Galactic Winds: Filament Formation and Emission Processes

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]

preprint2001arXiv

Theoretical Modeling of Starburst Galaxies

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.

preprint1997arXiv

Are Seyfert Narrow Line Regions Powered by Radio Jets?

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.