Jonathan Mackey

THREE EPOCHS OF STAR FORMATION IN THE HIGH-REDSHIFT UNIVERSE

We investigate the impact of an early population of massive stars on their surroundings. Dissociation of molecular hydrogen by strong UV emission from such stars is expected to produce a global transition in the cooling mechanism of minihalos at a redshift of approximately 30, strongly inhibiting star formation until more massive halos can collapse. Furthermore, chemical enrichment from Population III supernovae produces a second transition at z ~ 15-20, when the mean metallicity of the universe exceeds a critical threshold and Population III star formation gives way to Population II. We show that observations of high-redshift supernovae with the Next Generation Space Telescope have the potential to trace the cosmic star formation rate out to z 20, provided that Population III supernovae are at least as bright as, and ideally brighter than, Type Ia supernovae. We also propose a mechanism for the formation of a novel population of extremely low metallicity stars of intermediate mass at very high redshifts, which we term Population II.5. In our model shock compression, heating, and subsequent cooling to high density reduces the fragment mass in primordial gas to ~10 M☉, allowing low-mass stars to form. We predict the number density of relic Population II.5 stars in the Milky Way halo today and find that, with certain assumptions, there should be ~10 kpc-3 in the solar neighborhood.

Theoretical estimates of intrinsic galaxy alignment

ABSTRA C T It has recently been argued that the observed ellipticities of galaxies may be determined at least in part by the primordial tidal gravitational field in which the galaxy formed. Long-range correlations in the tidal field could thus lead to an ellipticity ‐ ellipticity correlation for widely separated galaxies. We present a new model relating ellipticity to angular momentum, which can be calculated in linear theory. We use this model to calculate the angular power spectrum of intrinsic galaxy shape correlations. We show that, for low-redshift galaxy surveys, our model predicts that intrinsic correlations will dominate correlations induced by weak lensing, in good agreement with previous theoretical work and observations. We find that our model produces ‘E-mode’ correlations enhanced by a factor of 3.5 over B-modes on small scales, making it harder to disentangle intrinsic correlations from those induced by weak gravitational lensing.

3D simulations of Betelgeuse’s bow shock

Betelgeuse, the bright, cool red supergiant in Orion, is moving supersonically relative to the local interstellar medium. The star emits a powerful stellar wind which collides with this medium, forming a cometary structure, a bow shock, pointing in the direction of motion. We present the first 3D hydrodynamic simulations of the formation and evolution of Betelgeuses bow shock. The models include realistic low temperature cooling and cover a range of plausible interstellar medium densities and stellar velocities between 0.3 - 1.9 cm-3 and 28 - 73 km/s. We show that the flow dynamics and morphology of the bow shock differ substantially due to the preferential growth of Rayleigh-Taylor or Kelvin-Helmholtz instabilities in the models. The former dominate the models with slow stellar velocities resulting in a clumpy bow shock sub-structure, whereas the latter produce a smoother, more layered sub-structure in the fast models. If the mass in the bow shock shell is low, as seems to be implied by the AKARI luminosities (~0.003 Msun), then Betelgeuses bow shock is very young and is unlikely to have reached a steady state. The circular nature of the bow shock shell is consistent with this conclusion. Thus, our results suggest that Betelgeuse entered the red supergiant phase only recently.

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 ionizing radiation from a point source. We show that shadowing of ionizing 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 ionization 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.

Interacting supernovae from photoionization-confined shells around red supergiant stars.

Betelgeuse, a nearby red supergiant, is a fast-moving 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 much 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 light curve, providing a natural explanation for the many supernovae that have signatures of circumstellar interaction.

Metamorphosis of Sn 2014c: Delayed Interaction Between a Hydrogen Poor Core-Collapse Supernova and a Nearby Circumstellar Shell

We present optical observations of supernova SN 2014C, which underwent an unprecedented slow metamorphosis from H-poor type Ib to H-rich type IIn over the course of one year. The observed spectroscopic evolution is consistent with the supernova having exploded in a cavity before encountering a massive shell of the progenitor stars stripped hydrogen envelope. Possible origins for the circumstellar shell include a brief Wolf-Rayet fast wind phase that overtook a slower red supergiant wind, eruptive ejection, or confinement of circumstellar material by external influences of neighboring stars. An extended high velocity Halpha absorption feature seen in near-maximum light spectra implies that the progenitor star was not completely stripped of hydrogen at the time of core collapse. Archival pre-explosion Subaru Telescope Suprime-Cam and Hubble Space Telescope Wide Field Planetary Camera 2 images of the region obtained in 2009 show a coincident source that is most likely a compact massive star cluster in NGC 7331 that hosted the progenitor system. By comparing the emission properties of the source with stellar population models that incorporate interacting binary stars we estimate the age of the host cluster to be 30 - 300 Myr, and favor ages closer to 30 Myr in light of relatively strong Halpha emission. SN 2014C is the best-observed member of a class of core-collapse supernovae that fill the gap between events that interact strongly with dense, nearby environments immediately after explosion and those that never show signs of interaction. Better understanding of the frequency and nature of this intermediate population can contribute valuable information about the poorly understood final stages of stellar evolution.

starbench: the D-type expansion of an H ii 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 H II region. The aim of this work is to understand the differences observed between the 12 participating numerical codes against the various analytical expressions examining the D-type phase of H II region expansion. To do this, we propose two well-defined tests which are tackled by 1D and 3D grid- and smoothed particle hydrodynamics-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 12 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 ≲ 2 per cent. 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 data set for calibrating the treatment of ionizing radiation hydrodynamics codes.

IRC -10414: a bow-shock-producing red supergiant star

Most runaway OB stars, like the majority of massive stars residing in their parent clusters, go through the red supergiant (RSG) phase during their lifetimes. Nonetheless, although many dozens of massive runaways were found to be associated with bow shocks, only two RSG bow-shock-producing stars, Betelgeuse and µCep, are known to date. In this paper, we report the discovery of an arc-like nebula around the late M-type star IRC 10414 using the SuperCOSMOS H-alpha Survey. Our spectroscopic follow-up of IRC 10414 with the Southern African Large Telescope (SALT) showed that it is a M7 supergiant, which supports previous claims on the RSG nature of this star based on observations of its maser emission. This was reinforced by our new radioand (sub)millimeter-wavelength molecular line observations made with the Atacama Pathfinder Experiment (APEX) 12 meter telescope and the Effelsberg 100 m radio telescope, respectively. The SALT spectrum of the nebula indicates that its emission is the result of shock excitation. This finding along with the arc-like shape of the nebula and an estimate of the space velocity of IRC 10414 (� 70 ± 20 kms 1 ) imply the bow shock interpretation for the nebula. Thus, IRC 10414 represents the third case of a bow-shock-producing RSG and the first one with a bow shock visible at optical wavelengths. We discuss the smooth appearance of the bow shocks around IRC 10414 and Betelgeuse and propose that one of the necessary conditions for stability of bow shocks generated by RSGs is the ionization of the stellar wind. Possible ionization sources of the wind of IRC 10414 are proposed and discussed.

Models of the circumstellar medium of evolving, massive runaway stars moving through the Galactic plane

At least 5 per cent of the massive stars are moving supersonically through the interstellar medium (ISM) and are expected to produce a stellar wind bow shock. We explore how the mass-loss and space velocity of massive runaway stars affect the morphology of their bow shocks. We run two-dimensional axisymmetric hydrodynamical simulations following the evolution of the circumstellar medium of these stars in the Galactic plane from the main sequence to the red supergiant phase. We find that thermal conduction is an important process governing the shape, size and structure of the bow shocks around hot stars, and that they have an optical luminosity mainly produced by forbidden lines, e.g. [O iii]. The Hα emission of the bow shocks around hot stars originates from near their contact discontinuity. The Hα emission of bow shocks around cool stars originates from their forward shock, and is too faint to be observed for the bow shocks that we simulate. The emission of optically thin radiation mainly comes from the shocked ISM material. All bow shock models are brighter in the infrared, i.e. the infrared is the most appropriate waveband to search for bow shocks. Our study suggests that the infrared emission comes from near the contact discontinuity for bow shocks of hot stars and from the inner region of shocked wind for bow shocks around cool stars. We predict that, in the Galactic plane, the brightest, i.e. the most easily detectable bow shocks are produced by high-mass stars moving with small space velocities.

Wind bubbles within H ii regions around slowly moving stars

Interstellar bubbles around O stars are driven by a combination of the stars 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 regions 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.