Field's linear analysis of thermal instability is repeated using methods related to Whitham's theory of wave hierarchies, which brings out the physically relevant parameters in a much clearer way than in the original analysis. It is also used for the stability of non-equilibrium states and we show that for gas cooling behind a shock, the usual analysis is only quantitatively valid for shocks that are just able to trigger a transition to the cold phase. A magnetic field can readily be included and we show that this does not change the stability criteria. By considering steady shock solutions, we show that almost all plausible initial conditions lead to a magnetically dominated state on the unstable part of the equilibrium curve. These results are used to analyse numerical calculations of perturbed steady shock solutions and of shocks interacting with a warm cloud.

This paper describes 3D simulations of the formation of collapsing cold clumps via thermal instability inside a larger cloud complex. The initial condition was a diffuse atomic, stationary, thermally unstable, 200pc diameter spherical cloud in pressure equilibrium with low density surroundings. This was seeded with 10% density perturbations at the finest initial grid level (0.29pc) around n_H = 1.1cm^{-3} and evolved with self-gravity included. No magnetic field was imposed. Resimulations at a higher resolution of a region extracted from this simulation (down to 0.039pc), show that the thermal instability forms sheets, then filaments and finally clumps. The width of the filaments increases over time, in one particular case from 0.26 to 0.56pc. Thereafter clumps with sizes of around 5pc grow at the intersections of filaments. 21 distinct clumps, with properties similar to those observed in molecular clouds, are found by using the FellWalker algorithm to find minima in the gravitational potential. Not all of these are gravitationally bound, but the convergent nature of the flow and increasing central density suggest they are likely to form stars. Further simulation of the most massive clump shows the gravitational collapse to a density >10^6 cm^{-3}. These results provide realistic initial conditions that can be used to study feedback in individual clumps, interacting clumps and the entire molecular cloud complex.

The famous Rosette Nebula has an evacuated central cavity formed from the stellar winds ejected from the 2-6 million-year-old co-distant and co-moving central star cluster NGC 2244. However, with upper age estimates of less than 110,000 years, the central cavity is too young compared to NGC 2244 and existing models do not reproduce its properties. A new proper motion study herein using Gaia data reveals the ejection of the most massive star in the Rosette, HD46223, from NGC 2244 occurred 1.73 (+0.34,-0.25)Myr (1$\sigma$ uncertainty) in the past. Assuming this ejection was at the birth of the most massive stars in NGC 2244, including the dominant centrally positioned HD46150, the age is set for the famous ionised region at more than ten times that derived for the cavity. Here, we are able to reproduce the structure of the Rosette Nebula, through simulation of mechanical stellar feedback from a 40M$_{\odot}$ star in a thin sheet-like molecular cloud. We form the 135,000M$_{\odot}$ cloud from thermally-unstable diffuse interstellar medium under the influence of a realistic background magnetic field with thermal/magnetic pressure equilibrium. Properties derived from a snapshot of the simulation at 1.5Myr, including cavity size, stellar age, magnetic field and resulting inclination to the line of sight, match those derived from observations. An elegant explanation is thus provided for the stark contrast in age estimates based on realistic diffuse ISM properties, molecular cloud formation and stellar wind feedback.

We have used the AMR hydrodynamic code, MG, to perform 3D hydrodynamic simulations with self-gravity of stellar feedback in a spherical clumpy molecular cloud formed through the action of thermal instability. We simulate the interaction of the mechanical energy input from 15 Msun, 40 Msun, 60 Msun and 120 Msun stars into a 100 pc-diameter 16,500 Msun cloud with a roughly spherical morphology with randomly distributed high density condensations. The stellar winds are introduced using appropriate non-rotating Geneva stellar evolution models. In the 15 Msun star case, the wind has very little effect, spreading around a few neighbouring clumps before becoming overwhelmed by the cloud collapse. In contrast, in the 40 Msun, 60 Msun and 120 Msun star cases, the more powerful stellar winds create large cavities and carve channels through the cloud, breaking out into the surrounding tenuous medium during the wind phase and considerably altering the cloud structure. After 4.97 Myrs, 3.97 Myrs and 3.01 Myrs respectively, the massive stars explode as supernovae (SNe). The wind-sculpted surroundings considerably affect the evolution of these SN events as they both escape the cloud along wind-carved channels and sweep up remaining clumps of cloud/wind material. The `cloud' as a coherent structure does not survive the SN from any of these stars, but only in the 120 Msun case is the cold molecular material completely destabilised and returned to the unstable thermal phase. In the 40 Msun and 60 Msun cases, coherent clumps of cold material are ejected from the cloud by the SN, potentially capable of further star formation.

We have used the AMR hydrodynamic code, MG, to perform 3D magnetohydrodynamic simulations with self-gravity of stellar feedback in a sheet-like molecular cloud formed through the action of the thermal instability. We simulate the interaction of the mechanical energy input from a 15 solar mass star and a 40 solar mass star into a 100 pc-diameter 17000 solar mass cloud with a corrugated sheet morphology that in projection appears filamentary. The stellar winds are introduced using appropriate Geneva stellar evolution models. In the 15 solar mass star case, the wind forms a narrow bipolar cavity with minimal effect on the parent cloud. In the 40 solar mass star case, the more powerful stellar wind creates a large cylindrical cavity through the centre of the cloud. After 12.5 Myrs and 4.97 Myrs respectively, the massive stars explode as supernovae (SNe). In the 15 solar mass star case, the SN material and energy is primarily deposited into the molecular cloud surroundings over ~10^5 years before the SN remnant escapes the cloud. In the 40 solar mass star case, a significant fraction of the SN material and energy rapidly escapes the molecular cloud along the wind cavity in a few tens of kiloyears. Both SN events compress the molecular cloud material around them to higher densities (so may trigger further star formation), and strengthen the magnetic field, typically by factors of 2-3 but up to a factor of 10. Our simulations are relevant to observations of bubbles in flattened ring-like molecular clouds and bipolar HII regions.

We have used the AMR hydrodynamic code, MG, to perform idealised 3D MHD simulations of the formation of clumpy and filamentary structure in a thermally unstable medium without turbulence. A stationary thermally unstable spherical diffuse atomic cloud with uniform density in pressure equilibrium with low density surroundings was seeded with random density variations and allowed to evolve. A range of magnetic field strengths threading the cloud have been explored, from beta=0.1 to beta=1.0 to the zero magnetic field case (beta=infinity), where beta is the ratio of thermal pressure to magnetic pressure. Once the density inhomogeneities had developed to the point where gravity started to become important, self-gravity was introduced to the simulation. With no magnetic field, clouds and clumps form within the cloud with aspect ratios of around unity, whereas in the presence of a relatively strong field (beta=0.1) these become filaments, then evolve into interconnected corrugated sheets that are predominantly perpendicular to the magnetic field. With magnetic and thermal pressure equality (beta=1.0), filaments, clouds and clumps are formed. At any particular instant, the projection of the 3D structure onto a plane parallel to the magnetic field, i.e. a line of sight perpendicular to the magnetic field, resembles the appearance of filamentary molecular clouds. The filament densities, widths, velocity dispersions and temperatures resemble those observed in molecular clouds. In contrast, in the strong field case beta=0.1, projection of the 3D structure along a line of sight parallel to the magnetic field reveals a remarkably uniform structure.

GALEX observations of the Mira AB binary system revealed a surrounding structure that has been successfully hydrodynamically interpreted as a bow shock and tail of ram-pressure-stripped material. Even the narrow tail, initially difficult to model, has been understood as the effect of the passage of Mira from a warm neutral medium into a hot, low-density medium, postulated to be the Local Bubble. However, no model to date has explained the observed kink and associated general curvature of the tail. We test the hypothesis that before entering the Local Bubble, Mira was travelling through a shear flow with approximately 1/3 Mira's own velocity at an angle of ~30degrees to Mira's proper motion. The hypothesis reproduces the kinked nature of Mira's tail and predicts recompression and reheating of the tail material to the same or greater levels of density and temperature predicted in the shock. This provides a heat source for the FUV emission, allowing for an extended lifetime of the FUV emission in line with other estimates of the age of the tail. The uniqueness of Mira's situation implies that the chances of observing other FUV tails behind AGB stars is highly unlikely.

We present the results of the search for candidate Planetary Nebulae interacting with the interstellar medium (PN-ISM) in the framework of the INT Photometric H$\alpha$ Survey (IPHAS) and located in the right ascension range 18h-20h. The detection capability of this new Northern survey, in terms of depth and imaging resolution, has allowed us to overcome the detection problem generally associated to the low surface brightness inherent to PNe-ISM. We discuss the detection of 21 IPHAS PN-ISM candidates. Thus, different stages of interaction were observed, implying various morphologies i.e. from the unaffected to totally disrupted shapes. The majority of the sources belong to the so-called WZO2 stage which main characteristic is a brightening of the nebula's shell in the direction of motion. The new findings are encouraging as they would be a first step into the reduction of the scarcity of observational data and they would provide new insights into the physical processes occurring in the rather evolved PNe.

The Hall effect is an important nonlinear mechanism affecting the evolution of magnetic fields in neutron stars. Studies of the governing equation, both theoretical and numerical, have shown that the Hall effect proceeds in a turbulent cascade of energy from large to small scales. We investigate the small-scale Hall instability conjectured to exist from the linear stability analysis of Rheinhardt and Geppert. Identical linear stability analyses are performed to find a suitable background field to model Rheinhardt and Geppert's ideas. The nonlinear evolution of this field is then modelled using a three-dimensional pseudospectral numerical MHD code. Combined with the background field, energy was injected at the ten specific eigenmodes with the greatest positive eigenvalues as inferred by the linear stability analysis. Energy is transferred to different scales in the system, but not into small scales to any extent that could be interpreted as a Hall instability. Any instabilities are overwhelmed by a late-onset turbulent Hall cascade, initially avoided by the choice of background field, but soon generated by nonlinear interactions between the growing eigenmodes. The Hall cascade is shown here, and by several authors elsewhere, to be the dominant mechanism in this system.

The interaction of planetary nebulae (PNe) with the interstellar medium as they move through it is now acknowledged to be a major shaping effect not just for ancient and large PNe, but also for relatively young PNe with high speed central stars. The most common effect is a rebrightening as the PN shell interacts with a pre-existing bow shock structure formed during the previous evolutionary phase of the central star. In this review, we consider this rebrightening in detail for the first time and discuss its origins, highlighting some observed examples. We go on to discuss the AGB star progenitors, reviewing the evidence for bow shock structures, and consider the progeny of rebrightened PNe - strongly disrupted objects which bear very little resemblance to typical PNe. Sh 2-68 is inferred to be perhaps the only documented case so far of such a PN.

Recent observations of the Mira AB binary system have revealed a surrounding arc-like structure and a stream of material stretching 2 degrees away in opposition to the arc. The alignment of the proper motion vector and the arc-like structure shows the structures to be a bow shock and accompanying tail. We have successfully hydrodynamically modelled the bow shock and tail as the interaction between the asymptotic giant branch (AGB) wind launched from Mira A and the surrounding interstellar medium. Our simulations show that the wake behind the bow shock is turbulent: this forms periodic density variations in the tail similar to those observed. We investigate the possiblity of mass-loss variations, but find that these have limited effect on the tail structure. The tail is estimated to be approximately 450,000 years old, and is moving with a velocity close to that of Mira itself. We suggest that the duration of the high mass-loss phase on the AGB may have been underestimated. Finally, both the tail curvature and the rebrightening at large distance can be qualitatively understood if Mira recently entered the Local Bubble. This is estimated to have occured 17 pc downstream from its current location.

Knots are commonly found in nearby planetary nebulae (PNe) and star forming regions. Within PNe, knots are often found to be associated with the brightest parts of the nebulae and understanding the physics involved in knots may reveal the processes dominating in PNe. As one of the closest PNe, the Helix Nebula (NGC 7293) is an ideal target to study such small-scale (~300 AU) structures. We have obtained infrared integral spectroscopy of a comet-shaped knot in the Helix Nebula using SINFONI on the Very Large Telescope at high spatial resolution (50-125 mas). With spatially resolved 2 micron spectra, we find that the H2 rotational temperature within the cometary knots is uniform. The rotational-vibrational temperature of the cometary knot (situated in the innermost region of the nebula, 2.5 arcmin away from the central star), is 1800 K, higher than the temperature seen in the outer regions (5-6 arcmin from the central star) of the nebula (900 K), showing that the excitation temperature varies across the nebula. The obtained intensities are reasonably well fitted with 27 km s-1 C-type shock model. This ambient gas velocity is slightly higher than the observed [HeII] wind velocity of 13 km s-1. The gas excitation can also be reproduced with a PDR (photo dominant region) model, but this requires an order of magnitude higher UV radiation. Both models have limitations, highlighting the need for models that treats both hydrodynamical physics and the PDR.

Interaction with the Interstellar Medium (ISM) cannot be ignored in understanding planetary nebula (PN) evolution and shaping. In an effort to understand the range of shapes observed in the outer envelopes of PNe, we have run a comprehensive set of three-dimensional hydrodynamic simulations, from the beginning of the asymptotic giant branch (AGB) superwind phase until the end of the post--AGB/PN phase. A 'triple-wind' model is used, including a slow AGB wind, fast post--AGB wind and third wind reflecting the linear movement through the ISM. A wide range of stellar velocities, mass-loss rates and ISM densities have been considered. We find ISM interaction strongly affects outer PN structures, with the dominant shaping occuring during the AGB phase. The simulations predict four stages of PN--ISM interaction whereby the PN is initially unaffected (1), then limb-brightened in the direction of motion (2), then distorted with the star moving away from the geometric centre (3) and finally so distorted that the object is no longer recognisable as a PN and may not be classed as such (4). Parsec-size shells around PN are predicted to be common. The structure and brightness of ancient PNe is largely determined by the ISM interaction, caused by rebrightening during the second stage; this effect may address the current discrepancies in Galactic PN abundance. The majority of PNe will have tail structures. Evidence for strong interaction is found for all known planetary nebulae in globular clusters.

Vortices have been postulated at a range of size scales in the universe including at the stellar size-scale. Whilst hydrodynamically simulating the wind from an asymptotic giant branch (AGB) star moving through and sweeping up its surrounding interstellar medium (ISM), we have found vortices on the size scale of 10^-1 pc to 10^1 pc in the wake of the star. These vortices appear to be the result of instabilities at the head of the bow shock formed upstream of the AGB star. The instabilities peel off downstream and form vortices in the tail of AGB material behind the bow shock, mixing with the surrounding ISM. We suggest such structures are visible in the planetary nebula Sh 2-188.

New Spitzer imaging observations have revealed the structure around the Mira variable star R Hya to be a one-sided parabolic arc 100 arcsec to the West stretching from North to South. We successfully model R Hya and its surroundings in terms of an interaction of the stellar wind from an asymptotic giant branch (AGB) star with the interstellar medium (ISM) the star moves through. Our three-dimensional hydrodynamic simulation reproduces the structure as a bow shock into the oncoming ISM. We propose this as another explanation of detached shells around such stars which should be considered alongside current theories of internal origin. The simulation predicts the existence of a tail of ram-pressure-stripped AGB material stretching downstream. Indications for such a tail behind R Hya are seen in IRAS maps.

Sh 2-188 is an example of strong interaction between a planetary nebula (PN) and the interstellar medium (ISM). Its structure is postulated to be the result of motion through the ISM. We present new H$\alpha$ images from the Isaac Newton Telescope Photometric H$\alpha$ Survey of the Northern Galactic Plane which reveal new structure. The nebula extends 15 arcmin on the sky in total. We have developed a `triple-wind' hydrodynamical model, comprising of the initial `slow' asymptotic giant branch (AGB) wind and the later `fast' stellar wind plus a third wind reflecting the motion through the ISM. Simulations at various velocities of the central star relative to the ISM indicate that a high velocity of 125 kms is required to reproduce the observed structure. We find most of the structure already forms during the AGB phase. The closure of the ring arises from the slow--fast wind interaction. Most of the mass lost on the AGB has been swept downstream, providing a potential explanation of the missing mass problem in PNe. We report a proper motion for the central star of 30 mas yr^-1 in the direction of the bright limb, implying a distance to the nebula of 850^{+500}_{-420} pc, consistent with a spectroscopic distances. Expansion velocities measured from spectroscopic data are consistent with velocities measured from the simulation. The model shows that the size of the PN was already set during the AGB phase.

Sh2-188 is thought to be an ancient planetary nebula in the galactic disk. It appears to be one-sided with recent observations revealing structure behind the filamentary limb. We postulate that Sh2-188 is interacting with the ISM and simulate it in terms of a ``triple-wind'' model comprising of the usual ``fast'' and ``slow'' interacting stellar winds plus the wind due to motion through the ISM. We have run simulations at various velocities of the central star relative to the ISM and find that a high velocity of 125 km/s best approximates the observed structure. We also suggest that Sh2-188 is younger than previously thought and that much of the mass lost on the AGB has been swept downstream.