Source author record

Jonathan Mackey

Jonathan Mackey 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.SR astro-ph.HE astro-ph.CO physics.comp-ph

Catalog footprint

What is connected

22works

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.

preprint2022arXiv

Core-collapse supernovae in dense environments -- particle acceleration and non-thermal emission

Supernova remnants are known to accelerate cosmic-rays from the detection of non-thermal emission in radio waves, X-rays, and gamma-rays. However, the ability to accelerate cosmic-rays up to PeV energies has yet to be demonstrated. The presence of cut-offs in the gamma-ray spectra of several young SNRs led to the idea that PeV energies might only be achieved during the first years of a remnant's evolution. We use our time-dependent acceleration-code RATPaC to study the acceleration of cosmic-rays in supernovae expanding into dense environments around massive stars. We performed spherically symmetric 1-D simulations in which we simultaneously solve the transport equations for cosmic-rays, magnetic turbulence, and the hydrodynamical flow of the thermal plasma in the test-particle limit. We investigated typical CSM parameters expected around RSG and LBV stars for freely expanding winds and accounted for the strong gamma-gamma-absorption in the first days after explosion. The maximum achievable particle energy is limited to below 600TeV even for largest considered values of the magnetic field and mass-loss rates. The maximum energy is not expected to surpass 200TeV and 70TeV for LBVs and RSGs that experience moderate mass-loss prior to the explosion. We find gamma-ray peak-luminosities consistent with current upper limits and evaluated that current-generation instruments are able to detect the gamma-rays from Type-IIP explosions at distances up to 60kpc and Type-IIn explosions up to 1.0Mpc. We also find a good agreement between the thermal X-ray and radio synchrotron emission predicted by our models with a range of observations.

preprint2022arXiv

Detailed study of extended gamma-ray morphology in the vicinity of the Coma cluster with Fermi-LAT

Galaxy clusters can be sources of high-energy (HE) $γ$-ray radiation, due to the efficient acceleration of particles exceeding EeV energies. At present, though, the only candidate for emitting HE $γ$-rays is the Coma cluster, towards which an excess of $γ$-ray emission has been detected by the Fermi Large Area Telescope (LAT). Using $\mathrm{\sim12.3}$ years of \textit{Fermi}-LAT data, we explored the region of the Coma cluster between energies 100 MeV and 1 TeV by detailed spectral and morphological analysis. In the region of the Coma cluster, we detected diffuse gamma-ray emission of energies between 100 MeV and 1 TeV with a 5.4$σ$ extension significance and a 68\% containment radius of $0.82^{+0.10}_{-0.05}$ degrees derived with a 2D homogeneous disk model. The corresponding gamma-ray spectrum extends up to $\sim50$ GeV, with a power-law index of $\mathrm{Γ=2.23\pm0.11}$ and flux of $\mathrm{(3.84\pm0.67)\times10^{-12}\,erg\,cm^{-2}\,s^{-1}}$. Using energy arguments we show that point-like sources such as radiogalaxies and star-forming galaxies are unlikely to explain the emission, and more likely, the emission is produced in the Coma cluster. Besides, we also identified three point-like sources in the region. However, due to the limited statistics of the detection, we could neither exclude nor conclude that the total extended emission is contributed to by these three-point like sources.

preprint2022arXiv

Multi-wavelength observations of the Blazar 4C +28.07

The active galactic nucleus 4C +28.07 is a flat spectrum radio quasar, one of the brightest at $γ$-ray energies. We study its multi-wavelength emission by analysing $\sim12.3$ years of \textit{Fermi-LAT} data in the $γ$-ray band and \textit{Swift-XRT/UVOT} available data in X-ray and Optical-to-Ultraviolet bands. In the $γ$-ray band, five flaring periods have been detected, during which the flux dramatically increases by several times (>5) compared with its average quiescent phase. Quasi-simultaneously with the flaring times, the X-ray and UVOT data detected by \textit{Swift-XRT/UVOT} have also been analysed. In one of the brightest flare periods (Flare 5; observed on Oct 12, 2018) the $γ$-ray flux reached $(6.7\pm0.81)\times 10^{-6}$ photon/cm^2/s ($\sim31\times$ higher than the mean flux over 12.3 years) with detection significance of $σ=6.1$. The apparent $γ$-ray luminosity of this flaring corresponds to $3.6\times10^{49}$ erg/s (for a distance of 8.38 Gpc), one of the highest $γ$-ray luminosities observed for blazars. Flare 5 has an estimated $\sim2$ hours time block, which can be considered the average $γ$-ray variability time. The variability time constrains the $γ$-ray emitting region size to <9e14 cm, which is close to the black hole radius. The spectral energy distributions (SEDs) in the $γ$-ray band for the $\sim12.3$ years of data show an early cut-off at $\sim14$ GeV; beyond $\sim60$ GeV, however, the spectrum hardens and is detected up to $\sim316$ GeV. Similar spectral behaviour is also noticeable for the SEDs of flares, which can be linked to the photon absorption by the emitting region's internal and external narrow-band radiation fields. Considering the significance of the obtained results from 4C\,+28.07, we compared the parameters with 3C\,279 and M87, to motivate further studies.

preprint2016arXiv

Supernova-blast waves in wind-blown bubbles, turbulent, and power-law ambient media

Supernova (SN) blast waves inject energy and momentum into the interstellar medium (ISM), control its turbulent multiphase structure and the launching of galactic outflows. Accurate modelling of the blast wave evolution is therefore essential for ISM and galaxy formation simulations. We present an efficient method to compute the input of momentum, thermal energy, and the velocity distribution of the shock-accelerated gas for ambient media with uniform (and with stellar wind blown bubbles), power-law, and turbulent density distributions. Assuming solar metallicity cooling, the blast wave evolution is followed to the beginning of the momentum conserving snowplough phase. The model recovers previous results for uniform ambient media. The momentum injection in wind-blown bubbles depend on the swept-up mass and the efficiency of cooling, when the blast wave hits the wind shell. For power-law density distributions with $n(r) \sim$ $r^{-2}$ (for $n(r) > n_{_{\rm floor}}$) the amount of momentum injection is solely regulated by the background density $n_{_{\rm floor}}$ and compares to $n_{_{\rm uni}}$ = $n_{_{\rm floor}}$. However, in turbulent ambient media with log-normal density distributions the momentum input can increase by a factor of 2 (compared to the homogeneous case) for high Mach numbers. The average momentum boost can be approximated as $p_{_{\rm turb}}/\mathrm{p_{_{0}}}\ =23.07\, \left(\frac{n_{_{0,\rm turb}}}{1\,{\rm cm}^{-3}}\right)^{-0.12} + 0.82 (\ln(1+b^{2}\mathcal{M}^{2}))^{1.49}\left(\frac{n_{_{0,\rm turb}}}{1\,{\rm cm}^{-3}}\right)^{-1.6}$. The velocity distributions are broad as gas can be accelerated to high velocities in low-density channels. The model values agree with results from recent, computationally expensive, three-dimensional simulations of SN explosions in turbulent media.

preprint2015arXiv

Cold gas in hot star clusters: the wind from the red supergiant W26 in Westerlund 1

The massive red supergiant (RSG) W26 in Westerlund 1 is one of a growing number of RSGs shown to have winds that are ionized from the outside in. The fate of this dense wind material is important for models of second generation star formation in massive star clusters. Mackey et al. (2014) showed that external photoionization can stall the wind of RSGs and accumulate mass in a dense static shell. We use 1D R-HD simulations of an externally photoionized wind to predict the Halpha and [NII] emission arising from photoionized winds both with and without a dense shell. We analyse spectra of the Halpha and [NII] emission in the environment around W26 and compare them with predicted synthetic emission. Simulations of slow winds that are decelerated into a dense shell show strongly limb-brightened line emission, with line radial velocities that are independent of the wind speed. Faster winds (>22 km/s) do not form a dense shell, have less limb-brightening, and the line radial velocity is a good tracer of the wind speed. The brightness of the [NII] and Halpha lines as a function of distance from W26 agrees reasonably well with observations when only the line flux is considered. The radial velocity disagrees, however: the brightest observed emission is blueshifted by ~25 km/s relative to the radial velocity of the star, whereas a spherically symmetric wind has the brightest emission at zero radial velocity. Our results show that the bright nebula surrounding W26 must be asymmetric; we suggest it is confined by external ram pressure from the wind of the nearby supergiant W9. We obtain a lower limit on the nitrogen abundance within the nebula of 2.35 times solar. The line ratio strongly favours photoionization over shock ionization, and so even if the observed nebula is pressure confined there should still be an ionization front and a photoionization-confined shell closer to the star.

preprint2015arXiv

Detecting stellar-wind bubbles through infrared arcs in HII regions

Mid-infrared arcs of dust emission are often seen near ionizing stars within HII regions. A possible explanations for these arcs is that they could show the outer edges of asymmetric stellar wind bubbles. We use two-dimensional, radiation-hydrodynamics simulations of wind bubbles within HII regions around individual stars to predict the infrared emission properties of the dust within the HII region. We assume that dust and gas are dynamically well-coupled and that dust properties (composition, size distribution) are the same in the HII region as outside it, and that the wind bubble contains no dust. We post-process the simulations to make synthetic intensity maps at infrared wavebands using the TORUS code. We find that the outer edge of a wind bubble emits brightly at 24um through starlight absorbed by dust grains and re-radiated thermally in the infrared. This produces a bright arc of emission for slowly moving stars that have asymmetric wind bubbles, even for cases where there is no bow shock or any corresponding feature in tracers of gas emission. The 24um intensity decreases exponentially from the arc with increasing distance from the star because the dust temperature decreases with distance. The size distribution and composition of the dust grains has quantitative but not qualitative effects on our results. Despite the simplifications of our model, we find good qualitative agreement with observations of the HII region RCW120, and can provide physical explanations for any quantitative differences. Our model produces an infrared arc with the same shape and size as the arc around CD -38 11636 in RCW120, and with comparable brightness. This suggests that infrared arcs around O stars in HII regions may be revealing the extent of stellar wind bubbles, although we have not excluded other explanations.

preprint2015arXiv

Revealing the binary origin of Type Ic superluminous supernovae through nebular hydrogen emission

We propose that nebular H-alpha emission as detected in the Type Ic superluminous supernova iPTF13ehe stems from matter which is stripped from a companion star when the supernova ejecta collide with it. The temporal evolution, the line broadening, and the overall blueshift of the emission are consistent with this interpretation. We scale the nebular H-alpha luminosity predicted for Type Ia supernovae in single-degenerate systems to derive the stripped mass required to explain the H-alpha luminosity of iPTF13ehe. We find a stripped mass of 0.1 - 0.9 solar masses, assuming that the supernova luminosity is powered by radioactivity or magnetar spin down. Because a central heating source is required to excite the H-alpha emission, an interaction-powered model is not favored for iPTF13ehe if the H-alpha emission is from stripped matter. We derive a companion mass of more than 20 solar masses and a binary separation of less than about 20 companion radii based on the stripping efficiency during the collision, indicating that the supernova progenitor and the companion formed a massive close binary system. If Type Ic superluminous supernovae generally occur in massive close binary systems, the early brightening observed previously in several Type Ic superluminous supernovae may also be due to the collision with a close companion. Observations of nebular hydrogen emission in future Type Ic superluminous supernovae will enable us to test this interpretation.

preprint2014arXiv

Effects of stellar evolution and ionizing radiation on the environments of massive stars

We discuss two important effects for the astrospheres of runaway stars: the propagation of ionizing photons far beyond the astropause, and the rapid evolution of massive stars (and their winds) near the end of their lives. Hot stars emit ionizing photons with associated photoheating that has a significant dynamical effect on their surroundings. 3D simulations show that HII regions around runaway O stars drive expanding conical shells and leave underdense wakes in the medium they pass through. For late O stars this feedback to the interstellar medium is more important than that from stellar winds. Late in life, O stars evolve to cool red supergiants more rapidly than their environment can react, producing transient circumstellar structures such as double bow shocks. This provides an explanation for the bow shock and linear bar-shaped structure observed around Betelgeuse.

preprint2014arXiv

Interacting supernovae from photoionization-confined shells around red supergiant stars

Betelgeuse, a nearby red supergiant, is a runaway star with a powerful stellar wind that drives a bow shock into its surroundings. This picture has been challenged by the discovery of a dense and almost static shell that is three times closer to the star than the bow shock and has been decelerated by some external force. The two physically distinct structures cannot both be formed by the hydrodynamic interaction of the wind with the interstellar medium. Here we report that a model in which Betelgeuse's wind is photoionized by radiation from external sources can explain the static shell without requiring a new understanding of the bow shock. Pressure from the photoionized wind generates a standing shock in the neutral part of the wind and forms an almost static, photoionization-confined shell. Other red supergiants should have significantly more massive shells than Betelgeuse, because the photoionization-confined shell traps up to 35 per cent of all mass lost during the red supergiant phase, confining this gas close to the star until it explodes. After the supernova explosion, massive shells dramatically affect the supernova lightcurve, providing a natural explanation for the many supernovae that have signatures of circumstellar interaction.

preprint2014arXiv

Pressure-driven fragmentation of clouds at high redshift

The discovery of a hyper metal-poor star with total metallicity of $\le 10^{-5}$ Z$_\odot$, has motivated new investigations of how such objects can form from primordial gas polluted by a single supernova. In this paper we present a shock-cloud model which simulates a supernova remnant interacting with a cloud in a metal-free environment at redshift $z=10$. Pre-supernova conditions are considered, which include a multiphase neutral medium and H II region. A small dense clump ($n=100$ cm$^{-3}$), located 40 pc from a 40 M$_\odot$ metal-free star, embedded in a $n=10$ cm$^{-3}$ ambient cloud. The evolution of the supernova remnant (explosion energy $10^{52}$ erg) and its subsequent interaction with the dense clump is examined. This is the first study to include a comprehensive treatment of the non-equilibrium chemistry and associated radiative cooling that is occurring at all stages of the shock-cloud model. We have included a primordial chemistry network that covers the temperature range $10-10^9$ K, and is coupled to thermal models of atomic & molecular cooling. We find $\times10^{3}$ density enhancement of the clump (i.e maximum density $\sim 78000$ cm$^{-3}$) within this metal-free model. This is consistent with Galactic shock-cloud models considering solar metallicity gas with equilibrium cooling functions. Despite this strong compression, the cloud does not become gravitationally unstable. We find that the small cloud modelled here is destroyed for shock velocities $\gtrsim 50\,$km s$^{-1}$, and not significantly affected by shocks with velocity $\lesssim 30\,$km s$^{-1}$. Rather specific conditions are required to make such a cloud collapse.

preprint2014arXiv

The Circumstellar Medium of Massive Stars in Motion

The circumstellar medium around massive stars is strongly impacted by stellar winds, radiation, and explosions. We use numerical simulations of these interactions to constrain the current properties and evolutionary history of various stars by comparison with observed circumstellar structures. Two- and three-dimensional simulations of bow shocks around red supergiant stars have shown that Betelgeuse has probably only recently evolved from a blue supergiant to a red supergiant, and hence its bow shock is very young and has not yet reached a steady state. We have also for the first time investigated the magnetohydrodynamics of the photoionised H II region around the nearby runaway O star Zeta Oph. Finally, we have calculated a grid of models of bow shocks around main sequence and evolved massive stars that has general application to many observed bow shocks, and which forms the basis of future work to model the explosions of these stars into their pre-shaped circumstellar medium.

preprint2014arXiv

Wind bubbles within H II regions around slowly moving stars

Interstellar bubbles around O stars are driven by a combination of the star's wind and ionizing radiation output. The wind contribution is uncertain because the boundary between the wind and interstellar medium is difficult to observe. Mid-infrared observations (e.g., of the H II region RCW 120) show arcs of dust emission around O stars, contained well within the H II region bubble. These arcs could indicate the edge of an asymmetric stellar wind bubble, distorted by density gradients and/or stellar motion. We present two-dimensional, radiation-hydrodynamics simulations investigating the evolution of wind bubbles and H II regions around massive stars moving through a dense (n=3000 cm^{-3}), uniform medium with velocities ranging from 4 to 16 km/s. The H II region morphology is strongly affected by stellar motion, as expected, but the wind bubble is also very aspherical from birth, even for the lowest space velocity considered. Wind bubbles do not fill their H II regions (we find filling factors of 10-20%), at least for a main sequence star with mass M~30 Msun. Furthermore, even for supersonic velocities the wind bow shock does not significantly trap the ionization front. X-ray emission from the wind bubble is soft, faint, and comes mainly from the turbulent mixing layer between the wind bubble and the H II region. The wind bubble radiates <1 per cent of its energy in X-rays; it loses most of its energy by turbulent mixing with cooler photoionized gas. Comparison of the simulations with the H II region RCW 120 shows that its dynamical age is <=0.4 Myr and that stellar motion <=4 km/s is allowed, implying that the ionizing source is unlikely to be a runaway star but more likely formed in situ. The region's youth, and apparent isolation from other O or B stars, makes it very interesting for studies of massive star formation and of initial mass functions.

preprint2013arXiv

3D Smoothed Particle Hydrodynamics Models of Betelgeuse's Bow Shock

Betelgeuse, the bright red supergiant (RSG) in Orion, is a runaway star. Its supersonic motion through the interstellar medium has resulted in the formation of a bow shock, a cometary structure pointing in the direction of motion. We present the first 3D hydrodynamic simulations of the formation and evolution of Betelgeuse's bow shock. We show that the bow shock morphology depends substantially on the growth timescale for Rayleigh-Taylor versus Kelvin-Helmholtz instabilities. We discuss our models in light of the recent Herschel, GALEX and VLA observations. If the mass in the bow shock shell is low (~few x 0.001 Msun), as seems to be implied by the AKARI and Herschel observations, then Betelgeuse's bow shock is very young and is unlikely to have reached a steady state. The circular, smooth bow shock shell is consistent with this conclusion. We further discuss the implications of our results, in particular, the possibility that Betelgeuse may have only recently entered the RSG phase.

preprint2013arXiv

Dynamics of H II regions around exiled O stars

At least 25 per cent of massive stars are ejected from their parent cluster, becoming runaways or exiles, travelling with often-supersonic space velocities through the interstellar medium (ISM). Their overpressurised H II regions impart kinetic energy and momentum to the ISM, compress and/or evaporate dense clouds, and can constrain properties of both the star and the ISM. Here we present one-, two-, and (the first) three-dimensional simulations of the H II region around a massive star moving supersonically through a uniform, magnetised ISM, with properties appropriate for the nearby O star Zeta Oph. The H II region leaves an expanding overdense shell behind the star and, inside this, an underdense wake that should be filled with hot gas from the shocked stellar wind. The gas column density in the shell is strongly influenced by the ISM magnetic field strength and orientation. H-alpha emission maps show the H II region remains roughly circular, although the star is displaced somewhat from the centre of emission. For our model parameters, the kinetic energy feedback from the H II region is comparable to the mechanical luminosity of the stellar wind, and the momentum feedback rate is >100X larger than that from the wind and about 10X larger than the total momentum input rate available from radiation pressure. Compared to the star's eventual supernova explosion, the kinetic energy feedback from the H II region over the star's main sequence lifetime is >100X less, but the momentum feedback is up to 4X larger. H II region dynamics are found to have only a small effect on the ISM conditions that a bow shock close to the star would encounter.

preprint2013arXiv

Numerical models for the circumstellar medium around Betelgeuse

The nearby red supergiant (RSG) Betelgeuse has a complex circumstellar medium out to at least 0.5 parsecs from its surface, shaped by its mass-loss history within the past 0.1 Myr, its environment, and its motion through the interstellar medium (ISM). In principle its mass-loss history can be constrained by comparing hydrodynamic models with observations. Observations and numerical simulations indicate that Betelgeuse has a very young bow shock, hence the star may have only recently become a RSG. To test this possibility we calculated a stellar evolution model for a single star with properties consistent with Betelgeuse. We incorporated the resulting evolving stellar wind into 2D hydrodynamic simulations to model a runaway blue supergiant (BSG) undergoing the transition to a RSG near the end of its life. The collapsing BSG wind bubble induces a bow shock-shaped inner shell which at least superficially resembles Betelgeuse's bow shock, and has a similar mass. Surrounding this is the larger-scale retreating bow shock generated by the now defunct BSG wind's interaction with the ISM. We investigate whether this outer shell could explain the bar feature located (at least in projection) just in front of Betelgeuse's bow shock.

preprint2012arXiv

Accuracy and Efficiency of Raytracing Photoionisation Algorithms

Three non-equilibrium photoionisation algorithms for hydrodynamical grid-based simulation codes are compared in terms of accuracy, timestepping criteria, and parallel scaling. Explicit methods with first-order time accuracy for photon conservation must use very restrictive timestep criteria to accurately track R-type ionisation fronts. A second-order accurate algorithm is described which, although it requires more work per step, allows much longer timesteps and is consequently more efficient. Implicit methods allow ionisation fronts to cross many grid cells per timestep while maintaining photon conservation accuracy. It is shown, however, that errors are much larger for multi-frequency radiation than for monochromatic radiation with the implicit algorithm used here, and large errors accrue when an ionisation front crosses many optical depths in a single step. The accuracy and convergence rates of the different algorithms are tested with a large number of timestepping criteria to identify the best criterion for each algorithm. With these criteria selected, the second-order explicit algorithm is the most efficient of the three, and its parallel scaling is significantly better than that of the implicit algorithm. The upgrade from first- to second-order accuracy in explicit algorithms could be made very simply to fixed-grid and adaptive mesh-refinement codes which currently use a first-order method.

preprint2012arXiv

Double bow shocks around young, runaway red supergiants: application to Betelgeuse

A significant fraction of massive stars are moving supersonically through the interstellar medium (ISM), either due to disruption of a binary system or ejection from their parent star cluster. The interaction of their wind with the ISM produces a bow shock. In late evolutionary stages these stars may undergo rapid transitions from red to blue and vice versa on the Hertzsprung-Russell diagram, with accompanying rapid changes to their stellar winds and bow shocks. Recent 3D simulations of the bow shock produced by the nearby runaway red supergiant (RSG) Betelgeuse, under the assumption of a constant wind, indicate that the bow shock is very young (<30000 years old), hence Betelgeuse may have only recently become a RSG. To test this possibility, we have calculated stellar evolution models for single stars which match the observed properties of Betelgeuse in the RSG phase. The resulting evolving stellar wind is incorporated into 2D hydrodynamic simulations in which we model a runaway blue supergiant (BSG) as it undergoes the transition to a RSG near the end of its life. We find that the collapsing BSG wind bubble induces a bow shock-shaped inner shell around the RSG wind that resembles Betelgeuse's bow shock, and has a similar mass. Surrounding this is the larger-scale retreating bow shock generated by the now defunct BSG wind's interaction with the ISM. We suggest that this outer shell could explain the bar feature located (at least in projection) just in front of Betelgeuse's bow shock.

preprint2012arXiv

Effects of Strong Magnetic Fields on Photoionised Clouds

Simulations are presented of the photoionisation of three dense gas clouds threaded by magnetic fields, showing the dynamical effects of different initial magnetic field orientations and strengths. For moderate magnetic field strengths the initial radiation-driven implosion phase is not strongly affected by the field geometry, and the photoevaporation flows are also similar. Over longer timescales, the simulation with an initial field parallel to the radiation propagation direction (parallel field) remains basically axisymmetric, whereas in the simulation with a perpendicular initial field the pillar of neutral gas fragments in a direction aligned with the magnetic field. For stronger initial magnetic fields, the dynamics in all gas phases are affected at all evolutionary times. In a simulation with a strong initially perpendicular field, photoevaporated gas forms filaments of dense ionised gas as it flows away from the ionisation front along field lines. These filaments are potentially a useful diagnostic of magnetic field strengths in H II regions because they are very bright in recombination line emission. In the strong parallel field simulation the ionised gas is constrained to flow back towards the radiation source, shielding the dense clouds and weakening the ionisation front, eventually transforming it to a recombination front.

preprint2012arXiv

Radiation-MHD simulations of pillars and globules in HII regions

Implicit and explicit raytracing-photoionisation algorithms have been implemented in the author's radiation-magnetohydrodynamics code. The algorithms are described briefly and their efficiency and parallel scaling are investigated. The implicit algorithm is more efficient for calculations where ionisation fronts have very supersonic velocities, and the explicit algorithm is favoured in the opposite limit because of its better parallel scaling. The implicit method is used to investigate the effects of initially uniform magnetic fields on the formation and evolution of dense pillars and cometary globules at the boundaries of HII regions. It is shown that for weak and medium field strengths an initially perpendicular field is swept into alignment with the pillar during its dynamical evolution, matching magnetic field observations of the`Pillars of Creation' in M16. A strong perpendicular magnetic field remains in its initial configuration and also confines the photoevaporation flow into a bar-shaped, dense, ionised ribbon which partially shields the ionisation front.

preprint2011arXiv

Radiation-MHD models of elephant trunks and globules in H II regions

We study the formation and evolution of pillars of dense gas, known as elephant trunks, at the boundaries of H II regions, formed by shadowing of ionising radiation by dense clumps. The effects of magnetic fields on this process are investigated using 3D radiation-magnetohydrodynamics simulations. For a simulation in which an initially uniform magnetic field of strength |B|=50 uG is oriented perpendicular to the radiation propagation direction, the field is swept into alignment with the pillar during its dynamical evolution, in agreement with observations of the "Pillars of Creation" in M16, and of some cometary globules. This effect is significantly enhanced when the simulation is re-run with a weaker field of 18 uG. A stronger field with |B|=160 uG is sufficient to prevent this evolution completely, also significantly affecting the photoionisation process. Using a larger simulation domain it is seen that the pillar formation models studied in Mackey & Lim (2010) ultimately evolve to cometary structures in the absence of dense gas further from the star.

preprint2010arXiv

Dynamical Models for the Formation of Elephant Trunks in H II Regions

The formation of pillars of dense gas at the boundaries of H II Regions is investigated with hydrodynamical numerical simulations including ionising radiation from a point source. We show that shadowing of ionising radiation by an inhomogeneous density field is capable of forming so-called elephant trunks (pillars of dense gas as in e.g. M16) without the assistance of self-gravity, or of ionisation front and cooling instabilities. A large simulation of a density field containing randomly generated clumps of gas is shown to naturally generate elephant trunks with certain clump configurations. These configurations are simulated in isolation and analysed in detail to show the formation mechanism and determine possible observational signatures. Pillars formed by the shadowing mechanism are shown to have rather different velocity profiles depending on the initial gas configuration, but asymmetries mean that the profiles also vary significantly with perspective, limiting their ability to discriminate between formation scenarios. Neutral and molecular gas cooling are shown to have a strong effect on these results.

preprint2010arXiv

Effects of Magnetic Fields on Photoionised Pillars and Globules

The effects of initially uniform magnetic fields on the formation and evolution of dense pillars and cometary globules at the boundaries of H II regions are investigated using 3D radiation-magnetohydrodynamics simulations. It is shown, in agreement with previous work, that a strong initial magnetic field is required to significantly alter the non-magnetised dynamics because the energy input from photoionisation is so large that it remains the dominant driver of the dynamics in most situations. Additionally it is found that for weak and medium field strengths an initially perpendicular field is swept into alignment with the pillar during its dynamical evolution, matching magnetic field observations of the `Pillars of Creation' in M16 and also some cometary globules. A strong perpendicular magnetic field remains in its initial configuration and also confines the photoevaporation flow into a bar-shaped dense ionised ribbon which partially shields the ionisation front and would be readily observable in recombination lines. A simple analytic model is presented to explain the properties of this bright linear structure. These results show that magnetic field strengths in star-forming regions can in principle be significantly constrained by the morphology of structures which form at the borders of H II regions.

Jonathan Mackey

What is connected

Connect this record

See the researcher in context

Building this map preview

22 published item(s)

Core-collapse supernovae in dense environments -- particle acceleration and non-thermal emission

Detailed study of extended gamma-ray morphology in the vicinity of the Coma cluster with Fermi-LAT

Multi-wavelength observations of the Blazar 4C +28.07

Supernova-blast waves in wind-blown bubbles, turbulent, and power-law ambient media

Cold gas in hot star clusters: the wind from the red supergiant W26 in Westerlund 1

Detecting stellar-wind bubbles through infrared arcs in HII regions

Revealing the binary origin of Type Ic superluminous supernovae through nebular hydrogen emission

Effects of stellar evolution and ionizing radiation on the environments of massive stars

Interacting supernovae from photoionization-confined shells around red supergiant stars

Pressure-driven fragmentation of clouds at high redshift

The Circumstellar Medium of Massive Stars in Motion

Wind bubbles within H II regions around slowly moving stars

3D Smoothed Particle Hydrodynamics Models of Betelgeuse's Bow Shock

Dynamics of H II regions around exiled O stars

Numerical models for the circumstellar medium around Betelgeuse

Accuracy and Efficiency of Raytracing Photoionisation Algorithms

Double bow shocks around young, runaway red supergiants: application to Betelgeuse

Effects of Strong Magnetic Fields on Photoionised Clouds

Radiation-MHD simulations of pillars and globules in HII regions

Radiation-MHD models of elephant trunks and globules in H II regions

Dynamical Models for the Formation of Elephant Trunks in H II Regions

Effects of Magnetic Fields on Photoionised Pillars and Globules