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 ob
servations (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.
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 Betelgeuses 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.
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 pho
tons 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.
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 compar
ison 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.
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 Betelgeuses bow shock, and has a similar mass. Surrounding this is the larger-scale retreating bow shock generated by the now defunct BSG winds interaction with the ISM. We investigate whether this outer shell could explain the bar feature located (at least in projection) just in front of Betelgeuses bow shock.
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 ra
diation-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.
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 Betelgeuses bow shock, and has a similar mass. Surrounding this is the larger-scale retreating bow shock generated by the now defunct BSG winds interaction with the ISM. We suggest that this outer shell could explain the bar feature located (at least in projection) just in front of Betelgeuses bow shock.
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 w
ith 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.
