Researcher profile

S. J. Arthur

S. J. Arthur contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
4topics
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

11 published item(s)

preprint2020arXiv

Dust in the Wolf-Rayet Nebula M1-67

The Wolf-Rayet nebula M1-67 around WR124 is located above the Galactic plane in a region mostly empty of interstellar medium, which makes it the perfect target to study the mass-loss episodes associated with the late stages of massive star evolution. Archive photometric observations from WISE, Spitzer (MIPS) and Herschel (PACS and SPIRE) are used to construct the spectral energy distribution (SED) of the nebula in the wavelength range of 12-500$μ$m. The infrared (photometric and spectroscopic) data and nebular optical data from the literature are modeled simultaneously using the spectral synthesis code Cloudy, where the free parameters are the gas density distribution and the dust grain size distribution. The infrared SED can be reproduced by dust grains with two size distributions: a MRN power-law distribution with grain sizes between 0.005 and 0.05$μ$m and a population of large grains with representative size 0.9$ μ$m. The latter points towards an eruptive origin for the formation of M1-67. The model predicts a nebular ionized gas mass of $M_\mathrm{ion} = 9.2^{+1.6}_{-1.5}~\mathrm{M}_\odot$ and the estimated mass-loss rate during the dust-formation period is $\dot{M} \approx 6 \times 10^{-4} \mathrm{M}_\odot$yr$^{-1}$. We discuss the implications of our results in the context of single and binary stellar evolution and propose that M1-67 represents the best candidate for a post-common envelope scenario in massive stars.

preprint2015arXiv

Diffuse X-ray Emission within Wolf-Rayet Nebulae

We discuss our most recent findings on the diffuse X-ray emission from Wolf-Rayet (WR) nebulae. The best-quality X-ray observations of these objects are those performed by XMM-Newton and Chandra towards S308, NGC2359, and NGC6888. Even though these three WR nebulae might have different formation scenarios, they all share similar characteristics: i) the main plasma temperatures of the X-ray-emitting gas is found to be $T$=[1-2]$\times$10$^{6}$ K, ii) the diffuse X-ray emission is confined inside the [O III] shell, and iii) their X-ray luminosities and electron densities in the 0.3-2.0~keV energy range are $L_\mathrm{X}\approx$10$^{33}$-10$^{34}$~erg~s$^{-1}$ and $n_\mathrm{e}\approx$0.1-1~cm$^{-3}$, respectively. These properties and the nebular-like abundances of the hot gas suggest mixing and/or thermal conduction is taking an important role reducing the temperature of the hot bubble.

preprint2015arXiv

StarBench: The D-type expansion of an HII region

StarBench is a project focused on benchmarking and validating different star-formation and stellar feedback codes. In this first StarBench paper we perform a comparison study of the D-type expansion of an HII region. The aim of this work is to understand the differences observed between the twelve participating numerical codes against the various analytical expressions examining the D-type phase of HII region expansion. To do this, we propose two well-defined tests which are tackled by 1D and 3D grid- and SPH- based codes. The first test examines the `early phase' D-type scenario during which the mechanical pressure driving the expansion is significantly larger than the thermal pressure of the neutral medium. The second test examines the `late phase' D-type scenario during which the system relaxes to pressure equilibrium with the external medium. Although they are mutually in excellent agreement, all twelve participating codes follow a modified expansion law that deviates significantly from the classical Spitzer solution in both scenarios. We present a semi-empirical formula combining the two different solutions appropriate to both early and late phases that agrees with high-resolution simulations to $\lesssim2\%$. This formula provides a much better benchmark solution for code validation than the Spitzer solution. The present comparison has validated the participating codes and through this project we provide a dataset for calibrating the treatment of ionizing radiation hydrodynamics codes.

preprint2015arXiv

X-ray emission from the Wolf-Rayet bubble NGC 6888. II. XMM-Newton EPIC observations

We present deep XMM-Newton EPIC observations of the Wolf-Rayet (WR) bubble NGC6888 around the star WR136. The complete X-ray mapping of the nebula confirms the distribution of the hot gas in three maxima spatially associated with the caps and northwest blowout hinted at by previous Chandra observations. The global X-ray emission is well described by a two-temperature optically thin plasma model $T_1$=1.4$\times$10$^{6}$ K, $T_{2}$=8.2$\times$10$^{6}$ K) with a luminosity of $L_{\mathrm{X}}$=7.8$\times$10$^{33}$ erg s$^{-1}$ in the 0.3--1.5 keV energy range. The rms electron density of the X-ray-emitting gas is estimated to be $n_\mathrm{e}$=0.4 cm$^{-3}$. The high-quality observations presented here reveal spectral variations within different regions in NGC6888, which allowed us for the first time to detect temperature and/or nitrogen abundance inhomogeneities in the hot gas inside a WR nebula. One possible explanation for such spectral variations is that the mixing of material from the outer nebula into the hot bubble is less efficient around the caps than in other nebular regions.

preprint2014arXiv

Formation and X-ray emission from Hot Bubbles in Planetary Nebulae. I. Hot Bubble formation

We carry out high resolution two-dimensional radiation-hydrodynamic numerical simulations to study the formation and evolution of hot bubbles inside planetary nebulae (PNe). We take into account the evolution of the stellar parameters, wind velocity and mass-loss rate from the final thermal pulses during the asymptotic giant branch (AGB) through to the post-AGB stage for a range of initial stellar masses. The instabilities that form at the interface between the hot bubble and the swept-up AGB wind shell lead to hydrodynamical interactions, photoevaporation flows and opacity variations. We explore the effects of hydrodynamical mixing combined with thermal conduction at this interface on the dynamics, photoionization, and emissivity of our models. We find that even models without thermal conduction mix significant amounts of mass into the hot bubble. When thermal conduction is not included, hot gas can leak through the gaps between clumps and filaments in the broken swept-up AGB shell and this depressurises the bubble. The inclusion of thermal conduction evaporates and heats material from the clumpy shell, which expands to seal the gaps, preventing a loss in bubble pressure. The dynamics of bubbles without conduction is dominated by the thermal pressure of the thick photoionized shell, while for bubbles with thermal conduction it is dominated by the hot, shocked wind.

preprint2014arXiv

Turbulence in simulated HII regions

We investigate the scale dependence of fluctuations inside a realistic model of an evolving turbulent HII region and to what extent these may be studied observationally. We find that the multiple scales of energy injection from champagne flows and the photoionization of clumps and filaments leads to a flatter spectrum of fluctuations than would be expected from top-down turbulence driven at the largest scales. The traditional structure function approach to the observational study of velocity fluctuations is shown to be incapable of reliably determining the velocity power spectrum of our simulation. We find that a more promising approach is the Velocity Channel Analysis technique of Lazarian & Pogosyan (2000), which, despite being intrinsically limited by thermal broadening, can successfully recover the logarithmic slope of the velocity power spectrum to a precision of +-0.1 from high resolution optical emission line spectroscopy.

preprint2012arXiv

X-ray emission from the Wolf-Rayet bubble S308

The Wolf-Rayet (WR) bubble S 308 around the WR star HD 50896 is one of the only two WR bubbles known to possess X-ray emission. We present XMM-Newton observations of three fields of this WR bubble that, in conjunction with an existing observation of its Northwest quadrant, map most of the nebula. The X-ray emission from S 308 displays a limb-brightened morphology, with a central cavity ~22 arcmin in size and a shell thickness of ~8 arcmin. This X-ray shell is confined by the optical shell of ionized material. The spectrum is dominated by the He-like triplets of NIV at 0.43 keV and OVII at 0.57 keV, and declines towards high energies, with a faint tail up to 1 keV. This spectrum can be described by a two-temperature optically thin plasma emission model (T1 ~ 1.1x10^6 K, T2 ~ 13x10^6 K), with a total X-ray luminosity ~2x10^33 erg/s at the assumed distance of 1.5 kpc.

preprint2011arXiv

Radiation-Hydrodynamic Models of the evolving Circumstellar Medium around Massive Stars

We study the evolution of the interstellar and circumstellar media around massive stars (M > 40M_{\odot}) from the main sequence through to the Wolf-Rayet stage by means of radiationhydrodynamic simulations. We use publicly available stellar evolution models to investigate the different possible structures that can form in the stellar wind bubbles around Wolf-Rayet stars. We find significant differences between models with and without stellar rotation, and between models from different authors. More specifically, we find that the main ingredients in the formation of structures in the Wolf-Rayet wind bubbles are the duration of the Red Supergiant (or Luminous Blue Variable) phase, the amount of mass lost, and the wind velocity during this phase, in agreement with previous authors. Thermal conduction is also included in our models. We find that main-sequence bubbles with thermal conduction are slightly smaller, due to extra cooling which reduces the pressure in the hot, shocked bubble, but that thermal conduction does not appear to significantly influence the formation of structures in post-main-sequence bubbles. Finally, we study the predicted X-ray emission from the models and compare our results with observations of the Wolf-Rayet bubbles S\,308, NGC\,6888, and RCW\,58. We find that bubbles composed primarily of clumps have reduced X-ray luminosity and very soft spectra, while bubbles with shells correspond more closely to observations.

preprint2011arXiv

Radiation-magnetohydrodynamic simulations of HII regions and their associated PDRs in turbulent molecular clouds

We present the results of radiation-magnetohydrodynamic simulations of the expansion of HII regions and surrounding photodissociation regions in turbulent, magnetised, molecular clouds on scales of up to 4 parsecs, including the effects of ionising and non-ionising ultraviolet radiation and x rays from young star clusters. We find that HII region expansion reduces the disordered component of the B field, imposing a large-scale order on the field around its border. The field in the neutral gas lies along the ionisation front, while the field in the ionised gas tends to be perpendicular to this. The highest pressure compressed neutral/molecular gas is driven towards approximate equipartition between thermal/magnetic/turbulent energy densities, whereas lower pressure neutral/molecular gas divides into quiescent, magnetically dominated regions, and, on the other hand, turbulent, demagnetised regions. The ionised gas shows approximate thermal/turbulent equipartition, but with magnetic energy densities 1 to 3 orders of magnitude lower. A high velocity dispersion (approx 8 km/s) is maintained in the ionised gas throughout our simulations, despite the mean expansion velocity being significantly lower. The B field does not significantly brake the HII region expansion on the length and timescales accessible to our simulations, but it does tend to suppress the small-scale fragmentation and radiation-driven implosion of neutral/molecular gas that forms globules and pillars at the edge of the HII region. However, the relative luminosity of ionising and non-ionising radiation has a much larger influence than the presence or absence of the B field. When the radiation field is relatively soft (as in the case of a lower mass cluster, with earliest spectral type of B0.5), then fragmentation is less vigorous and a thick, relatively smooth PDR forms. Movies available at http://youtube.com/user/divBequals0

preprint2006arXiv

Wind-Blown Bubbles around Evolved Stars

Most stars will experience episodes of substantial mass loss at some point in their lives. For very massive stars, mass loss dominates their evolution, although the mass loss rates are not known exactly, particularly once the star has left the main sequence. Direct observations of the stellar winds of massive stars can give information on the current mass-loss rates, while studies of the ring nebulae and HI shells that surround many Wolf-Rayet (WR) and luminous blue variable (LBV) stars provide information on the previous mass-loss history. The evolution of the most massive stars, (M > 25 solar masses), essentially follows the sequence O star to LBV or red supergiant (RSG) to WR star to supernova. For stars of mass less than 25 solar masses there is no final WR stage. During the main sequence and WR stages, the mass loss takes the form of highly supersonic stellar winds, which blow bubbles in the interstellar and circumstellar medium. In this way, the mechanical luminosity of the stellar wind is converted into kinetic energy of the swept-up ambient material, which is important for the dynamics of the interstellar medium. In this review article, analytic and numerical models are used to describe the hydrodynamics and energetics of wind-blown bubbles. A brief review of observations of bubbles is given, and the degree to which theory is supported by observations is discussed.

preprint2003arXiv

The Evolution of Mass Loaded Supernova Remnants. II. Temperature Dependent Mass Injection Rates

We investigate the evolution of spherically symmetric supernova remnants in which mass loading takes place due to conductively driven evaporation of embedded clouds. Numerical simulations reveal significant differences between the evolution of conductively mass loaded and the ablatively mass loaded remnants studied in Paper I. A main difference is the way in which conductive mass loading is extinguished at fairly early times, once the interior temperature of the remnant falls below ~10 million K. Thus, at late times remnants that ablatively mass load are dominated by loaded mass and thermal energy, while those that conductively mass load are dominated by swept-up mass and kinetic energy. Simple approximations to the remnant evolution, complementary to those in Paper I, are given.