At the time of formation, protoplanetary discs likely contain a comparable mass to their host protostars. As a result, gravitational instabilities (GIs) are expected to play a pivotal role in the early phases of disc evolution. However, as these young objects are heavily embedded, confirmation of GIs has remained elusive. Therefore, we use the radiative transfer code LIME to produce line images of a $0.17\,\mathrm{M}_{\odot}$ self-gravitating protosolar-like disc. We note the limitations of using LIME and explore methods to improve upon the default gridding. We synthesise noiseless observations to determine the sensitivities required to detect the spiral flux, and find that the line flux distribution does not necessarily correlate to the abundance density distribution; hence performing radiative transfer calculations is imperative. Moreover, the spiral features are seen in absorption, due to the GI-heated midplane and high extinction, which could be indicative of GI activity. If a small beamsize and appropriate molecular line are used then spatially resolving spirals in a protosolar-like disc should be possible with ALMA for an on-source time of 30 hr. Spectrally resolving non-axisymmetric structure takes only a tenth as long for a reasonable noise level, but attributing this structure to GI-induced activity would be tentative. Finally, we find that identifying finger-like features in PV diagrams of nearly edge-on discs, which are a direct indicator of spirals, is feasible with an on-source time of 19 hr, and hence likely offers the most promising means of confirming GI-driven spiral structure in young, embedded protoplanetary discs.

In this paper, we model the chemical evolution of a 0.25 M$_{\odot}$ protoplanetary disc surrounding a 1 M$_{\odot}$ star that undergoes fragmentation due to self-gravity. We use Smoothed Particle Hydrodynamics including a radiative transfer scheme, along with time-dependent chemical evolution code to follow the composition of the disc and resulting fragments over approximately 4000 yrs. Initially, four quasi-stable fragments are formed, of which two are eventually disrupted by tidal torques in the disc. From the results of our chemical modelling, we identify species that are abundant in the fragments (e.g. H$_{\rm 2}$O, H$_{\rm 2}$S, HNO, N$_{\rm 2}$, NH$_{\rm 3}$, OCS, SO), species that are abundant in the spiral shocks within the disc (e.g. CO, CH$_{\rm 4}$, CN, CS, H$_{\rm 2}$CO), and species which are abundant in the circumfragmentary material (e.g. HCO$^{\rm +}$). Our models suggest that in some fragments it is plausible for grains to sediment to the core before releasing their volatiles into the planetary envelope, leading to changes in, e.g., the C/O ratio of the gas and ice components. We would therefore predict that the atmospheric composition of planets generated by gravitational instability should not necessarily follow the bulk chemical composition of the local disc material.

Gravitational instabilities (GIs) are most likely a fundamental process during the early stages of protoplanetary disc formation. Recently, there have been detections of spiral features in young, embedded objects that appear consistent with GI-driven structure. It is crucial to perform hydrodynamic and radiative transfer simulations of gravitationally unstable discs in order to assess the validity of GIs in such objects, and constrain optimal targets for future observations. We utilise the radiative transfer code LIME to produce continuum emission maps of a $0.17\,\mathrm{M}_{\odot}$ self-gravitating protosolar-like disc. We note the limitations of using LIME as is and explore methods to improve upon the default gridding. We use CASA to produce synthetic observations of 270 continuum emission maps generated across different frequencies, inclinations and dust opacities. We find that the spiral structure of our protosolar-like disc model is distinguishable across the majority of our parameter space after 1 hour of observation, and is especially prominent at 230$\,$GHz due to the favourable combination of angular resolution and sensitivity. Disc mass derived from the observations is sensitive to the assumed dust opacities and temperatures, and therefore can be underestimated by a factor of at least 30 at 850$\,$GHz and 2.5 at 90$\,$GHz. As a result, this effect could retrospectively validate GIs in discs previously thought not massive enough to be gravitationally unstable, which could have a significant impact on the understanding of the formation and evolution of protoplanetary discs.

We describe the results of 3D simulations of the interaction of hydrodynamic shocks with Bonnor-Ebert spheres performed with an Adaptive Mesh Refinement code. The calculations are isothermal and the clouds are embedded in a medium in which the sound speed is either four or ten times that in the cloud. The strengths of the shocks are such that they induce gravitational collapse in some cases and not in others and we derive a simple estimate for the shock strength required for this to occur. These results are relevant to dense cores and Bok globules in star forming regions subjected to shocks produced by stellar feedback.

We numerically reproduce the density profiles for filaments that are in magnetohydrostatic and pressure equilibrium with their surroundings obtained in Tomisaka (2014) and show that these equilibria are dynamically stable. If the effect of ambipolar diffusion is considered, these filaments lose magnetic support initiating cloud collapse. The filaments do not lose magnetic flux. Rather the magnetic flux is redistributed within the filament from the centre towards the envelope. The rate of the collapse is inversely proportional to the fractional ionisation and two gravitationally-driven ambipolar diffusion regimes for the collapse are observed as predicted in Mouschovias & Morton (1991). For high values of the ionisation coefficient, that is $X \geq 10^{-7}$, the gas is strongly coupled to the magnetic field and the Jeans length is larger than the ambipolar diffusion length scale. Then the collapse is governed by magnetically-regulated ambipolar diffusion. For $X \lesssim 10^{-8}$, the gas is weakly coupled to the magnetic field and the magnetic support is removed by gravitationally-dominated ambipolar diffusion. Here, neutrals and ions only collide sporadically, that is the ambipolar diffusion length scale is larger than the Jeans length. When decaying turbulence is included, additional support is provided to the filament. This slows down the collapse of the filament even in the absence of a magnetic field. When a magnetic field is present, the collapse rate increases by a ratio smaller than for the non-magnetic case. This is because of a speed-up of the ambipolar diffusion due to larger magnetic field gradients generated by the turbulence and because the ambipolar diffusion aids the dissipation of turbulence below the ambipolar diffusion length scale. The highest increase in the rate is observed for the lowest ionisation coefficient and the highest turbulent intensity.

To date, most simulations of the chemistry in protoplanetary discs have used 1+1D or 2D axisymmetric $\alpha$-disc models to determine chemical compositions within young systems. This assumption is inappropriate for non-axisymmetric, gravitationally unstable discs, which may be a significant stage in early protoplanetary disc evolution. Using 3D radiative hydrodynamics, we have modelled the physical and chemical evolution of a 0.17 M$_{\odot}$ self-gravitating disc over a period of 2000 yr. The 0.8 M$_{\odot}$ central protostar is likely to evolve into a solar-like star, and hence this Class 0 or early Class I young stellar object may be analogous to our early Solar System. Shocks driven by gravitational instabilities enhance the desorption rates, which dominate the changes in gas-phase fractional abundances for most species. We find that at the end of the simulation, a number of species distinctly trace the spiral structure of our relatively low-mass disc, particularly CN. We compare our simulation to that of a more massive disc, and conclude that mass differences between gravitationally unstable discs may not have a strong impact on the chemical composition. We find that over the duration of our simulation, successive shock heating has a permanent effect on the abundances of HNO, CN and NH$_3$, which may have significant implications for both simulations and observations. We also find that HCO$^+$ may be a useful tracer of disc mass. We conclude that gravitational instabilities induced in lower mass discs can significantly, and permanently, affect the chemical evolution, and that observations with high-resolution instruments such as ALMA offer a promising means of characterising gravitational instabilities in protosolar discs.

We present 2D adiabatic magnetohydrodynamic (MHD) simulations of a shock interacting with groups of two or three cylindrical clouds. We study how the presence of a nearby cloud influences the dynamics of this interaction, and explore the resulting differences and similarities in the evolution of each cloud. The understanding gained from this small-scale study will help to interpret the behaviour of systems with many 10's or 100's of clouds. We observe a wide variety of behaviour in the interactions studied, which is dependent on the initial positions of the clouds and the orientation and strength of the magnetic field. We find: i) some clouds are stretched along their field-lines, whereas others are confined by their field-lines; ii) upstream clouds may accelerate past downstream clouds (though magnetic tension can prevent this); iii) clouds may also change their relative positions transverse to the direction of shock propagation as they "slingshot" past each other; iv) downstream clouds may be offered some protection from the oncoming flow as a result of being in the lee of an upstream cloud; v) the cycle of cloud compression and re-expansion is generally weaker when there are nearby neighbouring clouds; vi) the plasma $\beta$ in cloud material can vary rapidly as clouds collide with one another, but low values of $\beta$ are always transitory. This work is relevant to studies of multi-phase regions, where fast, low-density gas interacts with dense clouds, such as in circumstellar bubbles, supernova remnants, superbubbles and galactic winds.

We perform three dimensional hydrodynamic simulations of the ram pressure stripping of the hot extended gaseous halo of a massive galaxy using the k-epsilon sub-grid turbulence model at Mach numbers 0.9, 1.1 and 1.9. The k-epsilon model is used to simulate high Reynolds number flows by increasing the transport coefficients in regions of high turbulence. We find that the initial, instantaneous stripping is the same whether or not the k-epsilon model is implemented and is in agreement with the results of other studies. However the use of the k-epsilon model leads to five times less gas remaining after stripping by a supersonic flow has proceeded for 10 Gyr, which is more consistent with what simple analytic calculations indicate. Hence the continual Kelvin-Helmholtz stripping plays a significant role in the ram pressure stripping of the haloes of massive galaxies. To properly account for this, simulations of galaxy clusters will require the use of sub-grid turbulence models

The formation and earliest stages of protoplanetary discs remain poorly constrained by observations. ALMA will soon revolutionise this field. Therefore, it is important to provide predictions which will be valuable for the interpretation of future high sensitivity and high angular resolution observations. Here we present simulated ALMA observations based on radiative transfer modelling of a relatively massive (0.39 M_solar) self-gravitating disc embedded in a 10 M_solar dense core, with structure similar to the pre-stellar core L1544. We focus on simple species and conclude that C17O 3-2, HCO+ 3-2, OCS 26-25 and H2CO 404-303 lines can be used to probe the disc structure and kinematics at all scales.

An ideal magnetohydrodynamics (MHD) code with adaptive mesh refinement (AMR) was used to investigate the interactions of fast-mode shocks with self-gravitating, isothermal cores with mass-to-flux ratios that are somewhat below the minimum value required for gravitational collapse. We find that shock focussing produces colliding flows along the field lines that generate very high densities, even for relatively weak shocks. Self-gravity plays only a minor role in determining the highest density that is reached, but it does play a role in the subsequent evolution. The densities at comparable times differ by a factor of a few for shocks initially propagating perpendicularly or obliquely to the magnetic field in the ambient medium.

We present the first results for the sputtering of grain mantles and cores obtained with self-consistent multifluid hydromagnetic models of C-type shocks propagating through dusty media. The threshold shock speed for mantle sputtering is about 10 km/s and is independent of density. The mantles are completely vapourised in shocks with speeds of 20-25 km/s. At such shock speeds core sputtering commences and gas-phase SiO forms. Core destruction is not total in any C-type shock because grains are not completely destroyed in shocks with speeds near the minimum speeds at which J-type shocks appear. Due to the density-dependence of the critical shock speed for this transition, higher gas-phase SiO fractional abundances are produced behind shocks propagating in lower density gas. For shock speeds near the threshold speeds for both core and mantle sputtering, sputtering is much greater for shock velocities at smaller angles relative to the upstream magnetic field. At higher shock speeds, the angular variation is still present but less pronounced.

Research performed in the 1950s and 1960s by Leon Mestel on the roles of magnetic fields in star formation established the framework within which he and other key figures have conducted subsequent investigations on the subject. This short tribute to Leon contains a brief summary of some, but not all, of his ground breaking contributions in the area. It also mentions of some of the relevant problems that have received attention in the last few years. The coverage is not comprehensive, and the authors have drawn on their own results more and touched more briefly on those of others than they would in a normal review. Theirs is a personal contribution to the issue honouring Leon, one of the truly great gentlemen, wits, and most insightful of astrophysicists.

We present numerical simulations of the adiabatic interaction of a shock with a clumpy region containing many individual clouds. Our work incorporates a sub-grid turbulence model which for the first time makes this investigation feasible. We vary the Mach number of the shock, the density contrast of the clouds, and the ratio of total cloud mass to inter-cloud mass within the clumpy region. Cloud material becomes incorporated into the flow. This "mass-loading" reduces the Mach number of the shock, and leads to the formation of a dense shell. In cases in which the mass-loading is sufficient the flow slows enough that the shock degenerates into a wave. The interaction evolves through up to four stages: initially the shock decelerates; then its speed is nearly constant; next the shock accelerates as it leaves the clumpy region; finally it moves at a constant speed close to its initial speed. Turbulence is generated in the post-shock flow as the shock sweeps through the clumpy region. Clouds exposed to turbulence can be destroyed more rapidly than a similar cloud in an "isolated" environment. The lifetime of a downstream cloud decreases with increasing cloud-to-intercloud mass ratio. We briefly discuss the significance of these results for starburst superwinds and galaxy evolution.

In the absence of magnetic fields and cosmic rays, radiative cooling laws with a range of dependences on temperature affect the stability of interstellar gas. For about four and a half decades, astrophysicists have recognised the importance of the thermal instablity for the formation of clouds in the interstellar medium. Even in the past several years, many papers have concerned the role of the thermal instability in the production of molecular clouds. About three and a half decades ago, astrophysicists investigating radiative shocks noticed that for many cooling laws such shocks are unstable. Attempts to address the effects of cosmic rays on the stablity of radiative media that are initially uniform or that have just passed through shocks have been made. The simplest approach to such studies involves the assumption that the cosmic rays behave as a fluid. Work based on such an approach is described. Cosmic rays have no effect on the stability of initially uniform, static media with respect to isobaric perturbations, though they do affect the stability of such media with respect to isentropic perturbations. The effect of cosmic rays on the stability of radiative shocked media depends greatly on the efficiency of the conversion of energy in accelerated cosmic rays into thermal energy in the thermalized fluid. If that efficiency is low, radiative cooling makes weak shocks propagating into upstream media with low cosmic-ray pressures more likely to be cosmic-ray dominated than adiabatic shocks of comparable strength. The cosmic-ray dominated shocks do not display radiative overstability. Highly efficient conversion of cosmic-ray energy into thermal energy leads shocked media to behave as they do when cosmic rays are absent.

Until now, axisymmetric, alpha-disc models have been adopted for calculations of the chemical composition of protoplanetary discs. While this approach is reasonable for many discs, it is not appropriate when self-gravity is important. In this case, spiral waves and shocks cause temperature and density variations that affect the chemistry. We have adopted a dynamical model of a solar-mass star surrounded by a massive (0.39 Msun), self-gravitating disc, similar to those that may be found around Class 0 and early Class I protostars, in a study of disc chemistry. We find that for each of a number of species, e.g. H2O, adsorption and desorption dominate the changes in the gas-phase fractional abundance; because the desorption rates are very sensitive to temperature, maps of the emissions from such species should reveal the locations of shocks of varying strengths. The gas-phase fractional abundances of some other species, e.g. CS, are also affected by gas-phase reactions, particularly in warm shocked regions. We conclude that the dynamics of massive discs have a strong impact on how they appear when imaged in the emission lines of various molecular species.

We present the first observations of emission lines of CN(2-1), HCO$^{+}$(3-2) and C$_{2}$H(3-2) in the Perseus cluster. We observed at two positions: directly at the central galaxy, NGC 1275 and also at a position about 20$''$ to the east where associated filamentary structure has been shown to have strong CO emission. Clear detections in CN and HCO$^{+}$ transitions and a weak detection of the C$_{2}$H transition were made towards NGC 1275, while weak detections of CN and HCO$^{+}$ were made towards the eastern filamentary structure. Crude estimates of the column densities and fractional abundances (mostly upper limits) as functions of an unknown rotational temperature were made to both sources. These observational data were compared with the outputs of thermal/chemical models previously published by \citet{Baye10c} in an attempt to constrain the heating mechanisms in cluster gas. We find that models in which heating is dominated by cosmic rays can account for the molecular observations. This conclusion is consistent with that of \citet{Ferl09} in their study of gas traced by optical and infrared radiation. The cosmic ray heating rate in the regions probed by molecular emissions is required to be at least two orders of magnitude larger than that in the Milky Way.

Molecular line observations may serve as diagnostics of the degree to which the number density of cosmic ray protons, having energies of 10s to 100s of MeVs each, is enhanced in starburst galaxies and galaxies with active nuclei. Results, obtained with the UCL\_PDR code, for the fractional abundances of molecules as functions of the cosmic-ray induced ionisation rate, $\zeta$, are presented. The aim is not to model any particular external galaxies. Rather, it is to identify characteristics of the dependencies of molecular abundances on $\zeta$, in part to enable the development of suitable observational programmes for cosmic ray dominated regions (CRDRs) which will then stimulate detailed modelling. For a number density of hydrogen nuclei of of $10^4$ cm$^{-3}$, and high visual extinction, the fractional abundances of some species increase as $\zeta$ increases to $10^{-14}$ s$^{-1}$, but for much higher values of $\zeta$ the fractional abundances of all molecular species are significantly below their peak values. We show in particular that OH, H$_{2}$O, H$_{3}^{+}$, H$_{3}$O$^{+}$ and OH$^{+}$ attain large fractional abundances ($\geqslant 10^{-8}$) for $\zeta$ as large as $10^{-12}$ s$^{-1}$. HCO$^{+}$ is a poor tracer of CRDRs when $\zeta > 10^{-13}$ s$^{-1}$. Sulphur-bearing species may be useful tracers of CRDRs gas in which $\zeta \sim 10^{-16}$ s$^{-1}$. Ammonia has a large fractional abundance for $\zeta \leqslant 10^{-16}$ s$^{-1}$ and nitrogen appears in CN-bearing species at significant levels as $\zeta$ increases, even up to $\sim 10^{-14}$ s$^{-1}$. In this paper, we also discuss our model predictions, comparing them to recent detections in both galactic and extragalactic sources. We show that they agree well, to a first approximation, with the observational constraints.

Optical emission is detected from filaments around the central galaxies of clusters of galaxies. These filaments have lengths of tens of kiloparsecs. The emission is possibly due to heating caused by the dissipation of mechanical energy and by cosmic ray induced ionisation. CO millimeter and submillimeter line emissions as well as H$_{2}$ infrared emission originating in such filaments surrounding NGC~1275, the central galaxy of the Perseus cluster, have been detected. Our aim is to identify those molecular species, other than CO, that may emit detectable millimeter and submillimeter line features arising in these filaments, and to determine which of those species will produce emissions that might serve as diagnostics of the dissipation and cosmic ray induced ionisation. The time-dependent UCL photon-dominated region modelling code was used in the construction of steady-state models of molecular filamentary emission regions at appropriate pressures, for a range of dissipation and cosmic ray induced ionisation rates and incident radiation fields.HCO$^+$ and C$_2$H emissions will potentially provide information about the cosmic ray induced ionisation rates in the filaments. HCN and, in particular, CN are species with millimeter and submillimeter lines that remain abundant in the warmest regions containing molecules. Detections of the galaxy cluster filaments in HCO$^{+}$, C$_{2}$H, and CN emissions and further detections of them in HCN emissions would provide significant constraints on the dissipation and cosmic ray induced ionisation rates.

To understand the formation of a magnetically dominated molecular cloud from an atomic cloud, we study the interaction of a weak, radiative shock with a magnetised cloud. The thermally stable warm atomic cloud is initially in static equilibrium with the surrounding hot ionised gas. A shock propagating through the hot medium then interacts with the cloud. We follow the dynamical evolution of the shocked cloud with a time-dependent ideal magnetohydrodynamic code. By performing the simulations in 3D, we investigate the effect of different magnetic field orientations including parallel, perpendicular and oblique to the shock normal. We find that the angle between the shock normal and the magnetic field must be small to produce clouds with properties similar to observed molecular clouds.

The turbulent destruction of a cloud subject to the passage of an adiabatic shock is studied. We find large discrepancies between the lifetime of the cloud and the analytical result of Hartquist et al. (1986). These differences appear to be due to the assumption in Hartquist et al. that mass-loss occurs largely as a result of lower pressure regions on the surface of the cloud away from the stagnation point, whereas in reality Kelvin-Helmholtz (KH) instabilities play a dominant role in the cloud destruction. We find that the true lifetime of the cloud (defined as when all of the material from the core of the cloud is well mixed with the intercloud material in the hydrodynamic cells) is about 6 times t_KHD, where t_KHD is the growth timescale for the most disruptive, long-wavelength, KH instabilities. These findings have wide implications for diffuse sources where there is transfer of material between hot and cool phases. The properties of the interaction as a function of Mach number and cloud density contrast are also studied. The interaction is milder at lower Mach numbers with the most marked differences occuring at low shock Mach numbers when the postshock gas is subsonic with respect to the cloud (i.e. M < 2.76). Material stripped off the cloud only forms a long "tail-like" feature if the density contrast of the cloud to the ambient medium, chi > 1e3.

Outflows of young stars drive shocks into dusty, molecular regions. Most models of such shocks assume that they are steady and propagating perpendicular to the magnetic field. Real shocks often violate both of these assumptions and the media through which they propagate are inhomogeneous. We use the code employed previously to produce the first time-dependent simulations of fast-mode, oblique C-type shocks interacting with density perturbations. We include a self-consistent calculation of the thermal and ionisation balances and a fluid treatment of grains. We identify features that develop when a multifluid shock encounters a density inhomogeneity to investigate whether any part of the precursor region ever behaves in a quasi-steady fashion. If it does the shock may be modelled approximately without solving the time-dependent hydromagnetic equations. Simulations were made for initially steady oblique C-type shocks encountering density inhomogeneities. For a semi-finite inhomogeneity with a density larger than the surrounding medium, a transmitted shock evolves from being J-type to a steady C-type shock on a timescale comparable to the ion-flow time through it. A sufficiently upstream part of the precursor of an evolving J-type shock is quasi-steady. The ion-flow timescale is also relevant for the evolution of a shock moving into a region of decreasing density. The models for shocks propagating into regions in which the density increases and then decreases to its initial value cannot be entirely described in terms of the results obtained for monotonically increasing and decreasing densities. For the latter model, the long-term evolution to a C-type shock cannot be approximated by quasi-steady models.

Using a time-dependent multifluid, magnetohydrodynamic code, we calculated the structure of steady perpendicular and oblique C-type shocks in dusty plasmas. We included relevant processes to describe mass transfer between the different fluids, radiative cooling by emission lines and grain charging and studied the effect of single-sized and multiple sized grains on the shock structure. Our models are the first of oblique fast-mode molecular shocks in which such a rigorous treatment of the dust grain dynamics has been combined with a self-consistent calculation of the thermal and ionisation structures including appropriate microphysics. At low densities the grains do not play any significant role in the shock dynamics. At high densities, the ionisation fraction is sufficiently low that dust grains are important charge and current carriers and, thus, determine the shock structure. We find that the magnetic field in the shock front has a significant rotation out of the initial upstream plane. This is most pronounced for single-sized grains and small angles of the shock normal with the magnetic field. Our results are similar to previous studies of steady C-type shocks showing that our method is efficient, rigorous and robust. Unlike the method employed in the previous most detailed treatment of dust in steady oblique fast-mode shocks, ours allows a reliable calculation even when chemical or other conditions deviate from local statistical equilibrium. We are also able to model transient phenomena.

In 1986 Alex Dalgarno published a paper entitled "Is Interstellar Chemistry Useful?" By the middle 1970s, and perhaps even earlier, Alex had hoped that astronomical molecules would prove to: possess significant diagnostic utility; control many of the environments in which they exist; stimulate a wide variety of physicists and chemists who are at least as fascinated by the mechanisms forming and removing the molecules as by astronomy. His own research efforts have contributed greatly to the realization of that hope. This paper contains a few examples of: how molecules are used to diagnose large-scale dynamics in astronomical sources including star forming regions and supernovae; the ways in which molecular processes control the evolution of astronomical objects such as dense cores destined to become stars and very evolved giant stars; theoretical and laboratory investigations that elucidate the processes producing and removing astronomical molecules and allow their detection.

We present a time-dependent cosmic-ray modified shock model for which the calculated H-alpha emissivity profile agrees well with the H-alpha flux increase ahead of the Balmer-dominated shock at knot g in Tycho's supernova remnant, observed by Lee et al (2007). The backreaction of the cosmic ray component on the thermal component is treated in the two-fluid approximation, and we include thermal particle injection and energy transfer due to the acoustic instability in the precursor. The transient state of our model that describes the current state of the shock at knot g, occurs during the evolution from a thermal gas dominated shock to a smooth cosmic-ray dominated shock. Assuming a distance of 2.3 kpc to Tycho's remnant we obtain values for the cosmic ray diffusion coefficient, the injection parameter, and the time scale for the energy transfer of 10^{24} cm^{2} s^{-1}, 4.2x10^{-3}, and 426 y, respectively. We have also studied the parameter space for fast (300 km s^{-1} - 3000 km s^{-1}), time-asymptotically steady shocks and have identified a branch of solutions, for which the temperature in the cosmic ray precursor typically reaches 2-6x10^{4} K and the bulk acceleration of the flow through the precursor is less than 10 km s^{-1}. These solutions fall into the low cosmic ray acceleration efficiency regime and are relatively insensitive to shock parameters. This low cosmic ray acceleration efficiency branch of solutions may provide a natural explanation for the line broadening of the H-alpha narrow component observed in non-radiative shocks in many supernova remnants.

The interaction of a shock with a cloud has been extensively studied in the literature, where the effects of magnetic fields, radiative cooling and thermal conduction have been considered. However, the formation of fully developed turbulence has often been prevented by the artificial viscosity inherent in hydrodynamical simulations, and a uniform post-shock flow has been assumed in all previous single-cloud studies. In reality, the flow behind the shock is also likely to be turbulent, with non-uniform density, pressure and velocity structure created as the shock sweeps over inhomogenities upstream of the cloud. To address these twin issues we use a sub-grid compressible k-epsilon turbulence model to estimate the properties of the turbulence generated in shock-cloud interactions and the resulting increase in the transport coefficients that the turbulence brings. A detailed comparison with the output from an inviscid hydrodynamical code puts these new results into context. We find that cloud destruction in inviscid and k-epsilon models occurs at roughly the same speed when the post-shock flow is smooth and when the density contrast between the cloud and inter-cloud medium is less than 100. However, there are increasing and significant differences as this contrast increases. Clouds subjected to strong ``buffeting'' by a highly turbulent post-shock environment are destroyed significantly quicker. Additional calculations with an inviscid code where the post-shock flow is given random, grid-scale, motions confirms the more rapid destruction of the cloud. Our results clearly show that turbulence plays an important role in shock-cloud interactions, and that environmental turbulence adds a new dimension to the parameter space which has hitherto been studied (abridged).

Outflows of pre-main-sequence stars drive shocks into molecular material within 0.01 - 1 pc of the young stars. The shock-heated gas emits infrared, millimeter and submillimeter lines of many species including. Dust grains are important charge carriers and play a large role in coupling the magnetic field and flow of neutral gas. Some effects of the dust on the dynamics of oblique shocks began to emerge in the 1990s. However, detailed models of these shocks are required for the calculation of the grain sputtering contribution to gas phase abundances of species producing observed emissions. We are developing such models. Some of the molecular species introduced into the gas phase by sputtering in shocks or by thermally driven desorption in hot cores form on grain surfaces. Recently laboratory studies have begun to contribute to the understanding of surface reactions and thermally driven desorption important for the chemistry of star forming clouds. Dusty plasmas are prevalent in many evolved stars just as well as in star forming regions. Radiation pressure on dust plays a significant role in mass loss from some post-main-sequence stars. The mechanisms leading to the formation of carbonaceous dust in the stellar outflows are similar to those important for soot formation in flames. However, nucleation in oxygen-rich outflows is less well understood and remains a challenging research area. Dust is observed in supernova ejecta that have not passed through the reverse shocks that develop in the interaction of ejecta with ambient media. Dust is detected in high redshift galaxies that are sufficiently young that the only stars that could have produced the dust were so massive that they became supernovae. Consequently, the issue of the survival of dust in strong supernova shocks is of considerable interest.

We show that the mechanical energy injection rate generated as the intra-cluster medium (ICM) flows around cold clouds may be sufficient to power the optical and near infra-red emission of nebulae observed in the central regions of a sample of seven galaxy clusters. The energy injection rate is extremely sensitive to the velocity difference between the ICM and cold clouds, which may help to explain why optical and infra-red luminosity is often larger than expected in systems containing AGNs. We also find that mass recycling is likely to be important for the dynamics of the ICM. This effect will be strongest in the central regions of clusters where there is more than enough cold gas for its evaporation to contribute significantly to the density of the hot phase.

We aim to understand the formation of dense cores by magnetosonic waves in regions where the thermal to magnetic pressure ratio is small. Because of the low-ionisation fraction in molecular clouds, neutral and charged particles are weakly coupled. Ambipolar diffusion then plays an important role in the formation process. A quiescent, uniform plasma is perturbed by a fast-mode wave. Using 2D numerical simulations, we follow the evolution of the fast-mode wave. The simulations are done with a multifluid, adaptive mesh refinement MHD code. Initial perturbations with wavelengths that are 2 orders of magnitude larger than the dissipation length are strongly affected by the ion-neutral drift. Only in situations where there are large variations in the magnetic field corresponding to a highly turbulent gas can fast-mode waves generate dense cores. This means that, in most cores, no substructure can be produced. However, Core D of TMC-1 is an exception to this case. Due to its atypically high ionisation fraction, waves with wavelengths up to 3 orders of magnitude greater than the dissipation length can be present. Such waves are only weakly affected by ambipolar diffusion and can produce dense substructure without large wave-amplitudes. Our results also explain the observed transition from Alfvenic turbulent motion at large scales to subsonic motions at the level of dense cores.

Recent data support the idea that the filaments observed in H_alpha emission near the centres of some galaxy clusters were shaped by bulk flows within their intracluster media. We present numerical simulations of evaporated clump material interacting with impinging winds to investigate this possibility. In each simulation, a clump falls due to gravity while the drag of a wind retards the fall of evaporated material leading to elongation of the tail. However, we find that long filaments can only form if the outflowing wind velocity is sufficiently large, 10^8 cm s^-1. Otherwise, the tail material sinks almost as quickly as the cloud. For reasonable values of parameters, the morphological structure of a tail is qualitatively similar to those observed in clusters. Under certain conditions, the kinematics of the tail resemble those reported in Hatch et al.(2006). A comparison of the observations with the numerical results indicates that the filaments are likely to be a few tens of Myrs old. We also present arguments which suggest that the momentum transfer, from an outflowing wind, in the formation of these filaments is probably significant. As a result, tail formation could play a role in dissipating some of the energy injected by a central AGN close to the cluster centre where it is needed most. The trapping of energy by the cold gas may provide an additional feedback mechanism that helps to regulate the heating of the central regions of galaxy clusters and couple the AGN to the ICM.

Context: The detection of narrow SiO lines toward the young shocks of the L1448-mm outflow has been interpreted as a signature of the magnetic precursor of C-shocks. In contrast with the low SiO abundances (<10E-12) in the ambient gas, the narrow SiO emission at almost ambient velocities reveals enhanced SiO abundances of 10E-11. This enhancement has been proposed to be produced by the sputtering of the grain mantles at the first stages of C-shocks. However, modelling of the sputtering of grains has usually averaged the SiO abundances over the dissipation region of C-shocks, which cannot explain the recent observations. Aims: To model the evolution of the gas phase abundances of SiO, CH3OH and H2O, produced by the sputtering of grains as the shock propagates through the ambient gas. Methods: We propose a parametric model to describe the physical structure of C-shocks as a function of time. Using the known sputtering yields for water mantles (with minor constituents like silicon and CH3OH) and olivine cores by collisions with H2, He, C, O, Si, Fe and CO, we follow the evolution of the abundances of silicon, CH3OH and H2O ejected from grains. Results: The evolution of these abundances shows that CO seems to be the most efficient sputtering agent in low velocity shocks. The velocity threshold for the sputtering of silicon from the grain mantles is reduced by 5-10 km s-1 by CO compared to other models. The sputtering by CO can generate SiO abundances of 10E-11 at the early stages of low velocity shocks, consistent with those observed in the magnetic precursor of L1448-mm. Our model also satisfactorily reproduce the progressive enhancement of SiO, CH3OH and H2O observed in this outflow by the coexistence of two shocks with vs=30 and 60kms-1 within the same region.

Numerical hydrodynamics simulations have established that disks which are evolved under the condition of local isothermality will fragment into small dense clumps due to gravitational instabilities when the Toomre stability parameter $Q$ is sufficiently low. Because fragmentation through disk instability has been suggested as a gas giant planet formation mechanism, it is important to understand the physics underlying this process as thoroughly as possible. In this paper, we offer analytic arguments for why, at low $Q$, fragments are most likely to form first at the corotation radii of growing spiral modes, and we support these arguments with results from 3D hydrodynamics simulations.

To understand the formation of a magnetically dominated molecular cloud out of an atomic cloud, we follow the dynamical evolution of the cloud with a time-dependent axisymmetric magnetohydrodynamic code. A thermally stable warm atomic cloud is initially in static equilibrium with the surrounding hot ionised gas. A shock propagating through the hot medium interacts with the cloud. As a fast-mode shock propagates through the cloud, the gas behind it becomes thermally unstable. The $\beta$ value of the gas also becomes much smaller than the initial value of order unity. These conditions are ideal for magnetohydrodynamic waves to produce high-density clumps embedded in a rarefied warm medium. A slow-mode shock follows the fast-mode shock. Behind this shock a dense shell forms, which subsequently fragments. This is a primary region for the formation of massive stars. Our simulations show that only weak and moderate-strength shocks can form cold clouds which have properties typical of giant molecular clouds.

We have performed 2D hydromagnetic simulations with an adaptive mesh refinement code to examine the response of a pre-existing initially spherical dense core to a non-linear fast-mode wave. One key parameter is the ratio of the wavelength to the initial core radius. If that ratio is large and the wave amplitude is sufficient, significant compression of the core occurs, as envisaged by Myers & Lazarian (1998) in their "turbulent cooling flow" picture. For smaller values of that ratio, an initial value of the ratio of the thermal pressure to magnetic pressure of 0.2, and sufficiently large wave amplitude, the scattering induces the production of dense substructure in the core. This substructure may be related to that detected in the dense core associated with the cyanopolyyne peak in TMC-1. Our simulations also show that short-wavelength waves, contrary to large-wavelength waves, do not confine dense cores.

The gas equation of state may be one of the critical factors for the disk instability theory of gas giant planet formation. This letter addresses the treatment of H$_2$ in hydrodynamical simulations of gravitationally unstable disks. In our discussion, we point out possible consequences of erroneous specific internal energy relations, approximate specific internal energy relations with discontinuities, and assumptions of constant $\Gamma_1$. In addition, we consider whether the ortho/para ratio for H$_2$ in protoplanetary disks should be treated dynamically as if the species are in equilibrium. Preliminary simulations indicate that the correct treatment is particularly critical for the study of gravitational instability when $T = 30$-50 K.

Clouds containing molecular dark matter in quantities relevant for star formation may exist in minihaloes of the type of cold dark matter included in many cosmological simulations or in the regions of some galaxies extending far beyond their currently known boundaries. We have systematically explored parameter space to identify conditions under which plane-parallel clouds contain sufficient column densities of molecular dark matter that they could be reservoirs for future star formation. Such clouds would be undetected or at least appear by current observational criteria to be uninteresting from the perspective of star formation. We use a time-dependent PDR code to produce theoretical models of the chemistry and emission arising in clouds for our chosen region of parameter space. We then select a subset of model clouds with levels of emission that are low enough to be undetectable or at least overlooked by current surveys. The existence of significant column densities of cold molecular dark matter requires that the background radiation field be several or more orders of magnitude weaker than that in the solar neighbourhood. Lower turbulent velocities and cosmic ray induced ionization rates than typically associated with molecular material within a kpc of the Sun are also required for the molecular matter to be dark. We find that there is a large region within the parameter space that results in clouds that might contain a significant mass of molecular gas whilst remaining effectively undetectable or at least not particularly noticeable in surveys. We note briefly conditions under which molecular dark matter may contain a dynamically interesting mass.

Using two dimensional simulations, we study the formation of structures with a high-density contrast by magnetohydrodynamic waves in regions in which the ratio of thermal to magnetic pressure is small. The initial state is a uniform background perturbed by fast-mode wave. Our most significant result is that dense structures persist for far longer in a two-dimensional simulation than in the one-dimensional case. Once formed, these structures persist as long as the fast-mode amplitude remains high.

We present models for the chemistry in gas moving towards the ionization front of an HII region. When it is far from the ionization front, the gas is highly depleted of elements more massive than helium. However, as it approaches the ionization front, ices are destroyed and species formed on the grain surfaces are injected into the gas phase. Photodissociation removes gas phase molecular species as the gas flows towards the ionization front. We identify models for which the OH column densities are comparable to those measured in observations undertaken to study the magnetic fields in star forming regions and give results for the column densities of other species that should be abundant if the observed OH arises through a combination of the liberation of H2O from surfaces and photodissociation. They include CH3OH, H2CO, and H2S. Observations of these other species may help establish the nature of the OH spatial distribution in the clouds, which is important for the interpretation of the magnetic field results.

We consider the dynamical evolution of an interstellar cloud that is initially in thermal equilibrium in the warm phase and is then subjected to a sudden increase in the pressure of its surroundings. We find that if the initial plasma beta of the cloud is of order unity, then there is a considerable period during which the material in the cloud both has a small beta and is in the thermally unstable temperature range. These conditions are not only consistent with observations of star forming regions, but are also ideally suited to the production of density inhomogeneities by magnetohydrodynamic waves. The end result should be a cloud whose size and average density are typical of Giant Molecular Clouds (GMCs) and which contains denser regions whose densities are in the range inferred for the translucent clumps in GMCs.

Because of its key role in feedback in star formation and galaxy formation, we examine the nature of the interaction of a flow with discrete sources of mass injection. We show the results of two-dimensional numerical simulations in which we explore a range of configurations for the mass sources and study the effects of their proximity on the downstream flow. The mass sources act effectively as a single source of mass injection if they are so close together that the ratio of their combined mass injection rate is comparable to or exceeds the mass flux of the incident flow into the volume that they occupy. The simulations are relevant to many diffuse sources, such as planetary nebulae and starburst superwinds, in which a global flow interacts with material evaporating or being ablated from the surface of globules of cool, dense gas.

We investigate the Mach number dependence of the radiative overstability of shocks. We find that the stability of radiative shocks increases with decreasing Mach number, with the result that M=2 shocks require cooling exponents < -1.2 to be overstable. The frequency of oscillation of the fundamental mode has a strong Mach number dependence. We find that feedback between the cooling region and the cold dense layer (CDL) further downstream is a function of Mach number, with stronger feedback and oscillation of the boundary between the CDL and the cooling region occuring at lower Mach numbers. This feedback can be quantified in terms of the reflection coefficient of sound waves. An interesting finding is that the stability properties of low Mach number shocks can be dramatically altered if the shocked gas is able to cool to temperatures less than the pre-shock value, and it is probable that low Mach number astrophysical shocks will be overstable in a variety of situations. The results of this work are relevant to astrophysical shocks with low Mach numbers, such as supernova remnants (SNRs) immersed in a hot interstellar medium (e.g., within a starburst region), and shocks in molecular clouds, where time-dependent chemistry can lead to overstability. (Abridged)

Recent observations suggest that the behaviour of tracer species such as N_2H+ and CS is significantly different in regions of high and low mass star formation. In the latter, N_2H+ is a good tracer of mass, while CS is not. Observations show the reverse to be true in high-mass star formation regions. We use a computational chemical model to show that the abundances of these and other species may be significantly altered by a period of accelerated collapse in high mass star forming regions. We suggest these results provide a potential explanation of the observations, and make predictions for the behaviour of other species.

One aspect of supernova remnant evolution that is relatively unstudied is the influence of an AGN environment. A high density ambient medium and a nearby powerful continuum source will assist the cooling of shocked ejecta and swept-up gas. Motion of the surrounding medium relative to the remnant will also affect the remnant morphology. In an extension to previous work we have performed 2D hydrodynamical calculations of SNR evolution in an AGN environment, and have determined the evolutionary behaviour of cold gas in the remnant. The cold gas will contribute to the observed broad line emission in AGNs, and we present preliminary theoretical line profiles from our calculations. A more detailed comparison with observations will be performed in future work. The SNR-AGN interaction may be useful as a diagnostic of the AGN wind.

We present a linear stability analysis of a flow undergoing conductively-driven mass-loading from embedded clouds. We find that mass-loading damps isobaric and isentropic perturbations, and in this regard is similar to the effect of thermal conduction, but is much more pronounced where many embedded clumps exist. The stabilizing influence of mass-loading is wavelength independent against isobaric (condensing) perturbations, but wavelength dependent against isentropic (wave-like) perturbations. We derive equations for the degree of mass-loading needed to stabilize such perturbations. We have also made 1D numerical simulations of a mass-loaded radiative shock and demonstrated the damping of the overstability when mass-loading is rapid enough.

A model of the grain surface chemistry involving the accretion of atoms of two different elements, X and Y, and their reactions to form species X_2, XY, and Y_2 was examined for a wide range of choices for the values of its three free parameters - the accretion rate of X and Y, the desorption rate of X and the grain surface sweeping time of Y, all considered relative to the grain surface sweeping rate of X. Relative production rates of the diatomics were calculated with five methods involving, respectively, a high-order truncation of the master equation, a low-order truncation of the master equation, the standard deterministic rate equation approach, a modified rate equation approach and a set of approximations which are in some cases appropriate for accretion dominated chemistry. The accuracies of the relative production rates calculated with the different methods were assessed for the wide range of model parameters. The more accurate of the low-truncation master equation calculations and the standard deterministic rate equation approach gives results which are in most cases within ten or twenty per cent of the results given by the high-truncation master equation calculations. For many cases, the more accurate of the low order truncation and the standard deterministic rate equation approaches is indicated by a consideration of the average number of atoms of the two species on the grain's surface.

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.

We present similarity solutions for adiabatic bubbles that are blown by winds having time independent mechanical luminosities and that are each mass-loaded by the hydrodynamic ablation of distributed clumps. The mass loading is `switched-on' at a specified radius (with free-expansion of the wind interior to this point) and injects mass at a rate per unit volume proportional to M^delta r^lambda where delta = 4/3 (1) if the flow is subsonic (supersonic) with respect to the clumps. In the limit of negligible mass loading a similarity solution found by Dyson (1973) for expansion into a smooth ambient medium is recovered. The presence of mass loading heats the flow, which leads to a reduction in the Mach number of the supersonic freely-expanding flow, and weaker jump conditions across the inner shock. In solutions with large mass loading, it is possible for the wind to connect directly to the contact discontinuity without first passing through an inner shock. For a solution that gives the mass of swept-up ambient gas to be less than the sum of the masses of the wind and ablated material, lambda < -2. Maximum possible values for the ratio of ablated mass to wind mass occur when mass loading starts very close to the bubble center and when the flow is supersonic with respect to the clumps over the entire bubble radius. The maximum temperature in the bubble often occurs near the onset of mass loading, and in some cases can be many times greater than the post-inner-shock temperature. Our solutions are relevant to eg stellar wind-blown bubbles, galactic winds, etc. This work complements Pittard et al (2001) where it was assumed that clumps were evaporated through conductive energy transport.

We show that a cooled region of shocked supernova ejecta forms in a type II supernova-QSO wind interaction, and has a density, an ionization parameter, and a column density compatible with those inferred for the high ionization component of the broad emission line regions in QSOs. The calculations are based on the assumption that the ejecta flow is described initially by a similarity solution investigated by Chevalier (1982) and Nadyozhin (1985) and is spherically symmetric. Heating and cooling appropriate for gas irradiated by a nearby powerful continuum source is included in our model, together with reasonable assumptions for the properties of the QSO wind. The model results are also in agreement with observational correlations and imply reasonable supernova rates.

We present similarity solutions for mass-loaded adiabatic bubbles that are blown by winds having time dependent mechanical luminosities. We consider mass-loading through conductive evaporation of clumps. In the limit of little mass loading a similarity solution found by Dyson (1973) for expansion into a smooth ambient medium is recovered. We find that the Mach number in the shocked mass-loaded wind shows a radial dependence that varies qualitatively from solution to solution. In some cases it is everywhere less than unity in the frame of the clumps being evaporated, while in others it is everywhere greater than unity. In some solutions the mass-loaded shocked wind undergoes one or two sonic transitions in the clump frame. Maximum possible values of the ratio of evaporated mass to stellar wind mass are found as a consequence of the evaporation rates dependence on temperature and the lowering of the temperature by mass-loading. Mass-loading tends to reduce the emissivity in the interior of the bubble relative to its limb, whilst simultaneously increasing the central temperature relative to the limb temperature.

We present plane-parallel equilibrium models of molecular clumps that are supported by Alfven waves damped by the linear process of ion-neutral friction. We used a WKB approximation to treat the inward propagation of waves and adopted a realistic ionization structure influenced by dissociation and ionization due to photons of external origin. The model clumps are larger and less centrally condensed than those obtained for an assumed ionization structure, used in some previous studies, that is more appropriate for dark regions.

We present the jump conditions for ionization fronts with oblique magnetic fields. The standard nomenclature of R- and D-type fronts can still be applied, but in the case of oblique magnetization there are fronts of each type about each of the fast- and slow-mode speeds. As an ionization front slows, it will drive first a fast- and then a slow-mode shock into the surrounding medium. Even for rather weak upstream magnetic fields, the effect of magnetization on ionization front evolution can be important. [Includes numerical MHD models and an application to observations of S106.]