We present a numerical approach for accurately evolving a dust grain-size distribution undergoing number-conserving (such as sputtering) and/or mass-conserving (such as shattering) processes. As typically observed interstellar dust distributions follow a power-law, our method adopts a power-law discretisation and uses both the grain mass and number densities in each bin to determine the power-law parameters. This power-law method is complementary to piecewise-constant and linear methods in the literature. We find that the power-law method surpasses the other two approaches, especially for small bin numbers. In the sputtering tests the relative error in the total grain mass remains below 0.01% independent of the number of bins N, while the other methods only achieve this for N > 50 or higher. Likewise, the shattering test shows that the method also produces small relative errors in the total grain numbers while conserving mass. Not only does the power-law method conserve the global distribution properties, it also preserves the inter-bin characteristics so that the shape of the distribution is recovered to a high degree. This does not always happen for the constant and linear methods, especially not for small bin numbers. Implementing the power-law method in a hydrodynamical code thus minimises the numerical cost whilst maintaining high accuracy. The method is not limited to dust grain distributions, but can also be applied to the evolution of any distribution function, such as a cosmic-ray distribution affected by synchrotron radiation or inverse-Compton scattering.

We present the POL-2 850 $\mu$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $\mu$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $\mu$G and a criticality criterion $\lambda_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component.

We present the 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.

We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.

We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $\mu$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $\lambda$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical.

We present 850 $\mu$m imaging polarimetry data of the $\rho$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $\rho$ Oph A core; some of them can be part of a coherent magnetic field structure in the $\rho$ Oph region. The results are consistent with previous observations of the brightest regions of $\rho$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components.

Star formation from the interstellar medium of galactic disks is a basic process controlling the evolution of galaxies. Understanding the star formation rate in a local patch of a disk with a given gas mass is thus an important challenge for theoretical models. Here we simulate a kiloparsec region of a disk, following the evolution of self-gravitating molecular clouds down to subparsec scales, as they form stars that then inject feedback energy by dissociating and ionizing UV photons and supernova explosions. We assess the relative importance of each feedback mechanism. We find that $\rm H_2$-dissociating feedback results in the largest absolute reduction in star formation compared to the run with no feedback. Subsequently adding photoionization feedback produces a more modest reduction. Our fiducial models that combine all three feedback mechanisms yield, without fine-tuning, star formation rates that are in excellent agreement with observations, with $\rm H_2$-dissociating photons playing a crucial role. Models that only include supernova feedback---a common method in galaxy evolution simulations---settle to similar star formation rates, but with very different temperature and chemical states of the gas, and with very different spatial distributions of young stars.

We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments.

We study giant molecular cloud (GMC) collisions and their ability to trigger star cluster formation. We further develop our three dimensional magnetized, turbulent, colliding GMC simulations by implementing star formation sub-grid models. Two such models are explored: (1) "Density-Regulated," i.e., fixed efficiency per free-fall time above a set density threshold; (2) "Magnetically-Regulated," i.e., fixed efficiency per free-fall time in regions that are magnetically supercritical. Variations of parameters associated with these models are also explored. In the non-colliding simulations, the overall level of star formation is sensitive to model parameter choices that relate to effective density thresholds. In the GMC collision simulations, the final star formation rates and efficiencies are relatively independent of these parameters. Between non-colliding and colliding cases, we compare the morphologies of the resulting star clusters, properties of star-forming gas, time evolution of the star formation rate (SFR), spatial clustering of the stars, and resulting kinematics of the stars in comparison to the natal gas. We find that typical collisions, by creating larger amounts of dense gas, trigger earlier and enhanced star formation, resulting in 10 times higher SFRs and efficiencies. The star clusters formed from GMC collisions show greater spatial sub-structure and more disturbed kinematics.

We investigate giant molecular cloud (GMCs) collisions and their ability to induce gravitational instability and thus star formation. This mechanism may be a major driver of star formation activity in galactic disks. We carry out a series of three dimensional, magnetohydrodynamics (MHD), adaptive mesh refinement (AMR) simulations to study how cloud collisions trigger formation of dense filaments and clumps. Heating and cooling functions are implemented based on photo-dissociation region (PDR) models that span the atomic to molecular transition and can return detailed diagnostic information. The clouds are initialized with supersonic turbulence and a range of magnetic field strengths and orientations. Collisions at various velocities and impact parameters are investigated. Comparing and contrasting colliding and non-colliding cases, we characterize morphologies of dense gas, magnetic field structure, cloud kinematic signatures, and cloud dynamics. We present key observational diagnostics of cloud collisions, especially: relative orientations between magnetic fields and density structures, like filaments; $^{13}$CO($J$=2-1), $^{13}$CO($J$=3-2), and $^{12}$CO($J$=8-7) integrated intensity maps and spectra; and cloud virial parameters. We compare these results to observed Galactic clouds.

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.

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.

We utilize magnetohydrodynamic (MHD) simulations to develop a numerical model for GMC-GMC collisions between nearly magnetically critical clouds. The goal is to determine if, and under what circumstances, cloud collisions can cause pre-existing magnetically subcritical clumps to become supercritical and undergo gravitational collapse. We first develop and implement new photodissociation region (PDR) based heating and cooling functions that span the atomic to molecular transition, creating a multiphase ISM and allowing modeling of non-equilibrium temperature structures. Then in 2D and with ideal MHD, we explore a wide parameter space of magnetic field strength, magnetic field geometry, collision velocity, and impact parameter, and compare isolated versus colliding clouds. We find factors of ~2-3 increase in mean clump density from typical collisions, with strong dependence on collision velocity and magnetic field strength, but ultimately limited by flux-freezing in 2D geometries. For geometries enabling flow along magnetic field lines, greater degrees of collapse are seen. We discuss observational diagnostics of cloud collisions, focussing on 13CO(J=2-1), 13CO(J=3-2), and 12CO(J=8-7) integrated intensity maps and spectra, which we synthesize from our simulation outputs. We find the ratio of J=8-7 to lower-J emission is a powerful diagnostic probe of GMC collisions.

We present hydrodynamic simulations of self-gravitating dense gas in a galactic disk, exploring scales ranging from 1 kpc down to $\sim 0.1$~pc. Our primary goal is to understand how dense filaments form in Giant Molecular Clouds (GMCs). These structures, often observed as Infrared Dark Clouds (IRDCs) in the Galactic plane, are thought to be the precursors to massive stars and star clusters, so their formation may be the rate limiting step controlling global star formation rates in galactic systems as described by the Kennicutt-Schmidt relation. Our study follows on from Van Loo et al. (2013, Paper I), which carried out simulations to 0.5~pc resolution and examined global aspects of the formation of dense gas clumps and the resulting star formation rate. Here, using our higher resolution, we examine the detailed structural, kinematic and dynamical properties of dense filaments and clumps, including mass surface density ($\Sigma$) probability distribution functions, filament mass per unit length and its dispersion, lateral $\Sigma$ profiles, filament fragmentation, filament velocity gradients and infall, and degree of filament and clump virialization. Where possible, these properties are compared to observations of IRDCs. By many metrics, especially too large mass fractions of high $\Sigma>1\:{\rm g\:cm^{-2}}$ material, too high mass per unit length dispersion due to dense clump formation, too high velocity gradients and too high velocity dispersion for a given mass per unit length, the simulated filaments differ from observed IRDCs. We thus conclude that IRDCs do not form from global fast collapse of GMCs. Rather, we expect IRDC formation and collapse is slowed significantly by the influence of dynamically important magnetic fields, which may thus play a crucial role in regulating galactic star formation rates.

The regulation of galactic-scale star formation rates (SFRs) is a basic problem for theories of galaxy formation and evolution: which processes are responsible for making observed star formation rates so inefficient compared to maximal rates of gas content divided by dynamical timescale? Here we study the effect of magnetic fields of different strengths on the evolution of giant molecular clouds (GMCs) within a kiloparsec patch of a disk galaxy and resolving scales down to $\simeq0.5\:{\rm{pc}}$. Including an empirically motivated prescription for star formation from dense gas ($n_{\rm{H}}>10^5\:{\rm{cm}^{-3}}$) at an efficiency of 2\% per local free-fall time, we derive the amount of suppression of star formation by magnetic fields compared to the nonmagnetized case. We find GMC fragmentation, dense clump formation and SFR can be significantly affected by the inclusion of magnetic fields, especially in our strongest investigated $B$-field case of $80\:{\rm{\mu}}$G. However, our chosen kpc-scale region, extracted from a global galaxy simulation, happens to contain a starbursting cloud complex that is only modestly affected by these magnetic fields and likely requires internal star formation feedback to regulate its SFR.

We examine the role of the gravitational instability in an isothermal, self-gravitating layer threaded by magnetic fields on the formation of filaments and dense cores. Using numerical simulation we follow the non-linear evolution of a perturbed equilibrium layer. The linear evolution of such a layer is described in the analytic work of Nagai et al (1998). We find that filaments and dense cores form simultaneously. Depending on the initial magnetic field, the resulting filaments form either a spiderweb-like network (for weak magnetic fields) or a network of parallel filaments aligned perpendicular to the magnetic field lines (for strong magnetic fields). Although the filaments are radially collapsing, the density profile of their central region (up to the thermal scale height) can be approximated by a hydrodynamical equilibrium density structure. Thus, the magnetic field does not play a significant role in setting the density distribution of the filaments. The density distribution outside of the central region deviates from the equilibrium. The radial column density distribution is then flatter than the expected power law of $r^{-4}$ and similar to filament profiles observed with Herschel. Our results does not explain the near constant filament width of $\sim 0.1$pc. However, our model does not include turbulent motions. It is expected that accretion-driven amplification of these turbulent motions provides additional support within the filaments against gravitational collapse. Finally, we interpret the filamentary network of the massive star forming complex G14.225-0.506 in terms of the gravitational instability model and find that the properties of the complex are consistent with being formed out of an unstable layer threaded by a strong, parallel magnetic field.

All stars are born in molecular clouds, and most in giant molecular clouds (GMCs), which thus set the star formation activity of galaxies. We first review their observed properties, including measures of mass surface density, Sigma, and thus mass, M. We discuss cloud dynamics, concluding most GMCs are gravitationally bound. Star formation is highly clustered within GMCs, but overall is very inefficient. We compare properties of star-forming clumps with those of young stellar clusters (YSCs). The high central densities of YSCs may result via dynamical evolution of already-formed stars during and after star cluster formation. We discuss theoretical models of GMC evolution, especially addressing how turbulence is maintained, and emphasizing the importance of GMC collisions. We describe how feedback limits total star formation efficiency, epsilon, in clumps. A turbulent and clumpy medium allows higher epsilon, permitting formation of bound clusters even when escape speeds are less than the ionized gas sound speed.

We present hydrodynamic simulations of the evolution of self-gravitating dense gas on scales of 1 kiloparsec down to < parsec in a galactic disk, designed to study dense clump formation from giant molecular clouds (GMCs). These structures are expected to be the precursors to star clusters and this process may be the rate limiting step controling star formation rates in galactic systems as described by the Kennicutt-Schmidt relation. We follow the thermal evolution of the gas down to ~5K using extinction-dependent heating and cooling functions. We do not yet include magnetic fields or localized stellar feedback, so the evolution of the GMCs and clumps is determined solely by self-gravity balanced by thermal and turbulent pressure support and the large scale galactic shear. While cloud structures and densities change significantly during the simulation, GMC virial parameters remain mostly above unity for time scales exceeding the free-fall time of GMCs indicating that energy from galactic shear and large-scale cloud motions continuously cascades down to and within the GMCs. We implement star formation at a slow, inefficient rate of 2% per local free-fall time, but even this yields global star formation rates that are about two orders of magnitude larger than the observed Kennicutt-Schmidt relation due to over-production of dense gas clumps. We expect a combination of magnetic support and localized stellar feedback is required to inhibit dense clump formation to ~1% of the rate that results from the nonmagnetic, zero-feedback limit.

We present simulations of the evolution of self-gravitating dense gas on kiloparsec-size scales in a galactic disk, designed to study dense clump formation from giant molecular clouds (GMCs). These dense clumps are expected to be the precursors to star clusters and this process may be the rate limiting step controling star formation rates in galactic systems as described by the Kennicutt-Schmidt relation. The evolution of these simulated GMCs and clumps is determined by self-gravity balanced by turbulent pressure support and the large scale galactic shear. While the cloud structures and densities significantly change during their evolution, they remain roughly in virial equilibrium for time scales exceeding the free-fall time of GMCs, indicating that energy from the galactic shear continuously cascades down. We implement star formation at a slow, inefficient rate of 2% per local free-fall time, but this yields global star formation rates that are more than ~two orders of magnitude larger than the observed Kennicutt-Schmidt relation due to the over-production of dense clump gas. To explain this discrepancy, we anticipate magnetic fields to provide additional support. Low-resolution simulations indeed show that the magnetic field reduces the star formation rate.

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.

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.

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.

We present high-sensitivity 2'x4' maps of the J=2-1 rotational lines of SiO, CO, 13CO and C18O, observed toward the filamentary Infrared Dark Cloud (IRDC) G035.39-00.33. Single-pointing spectra of the SiO J=2-1 and J=3-2 lines toward several regions in the filament, are also reported. The SiO images reveal that SiO is widespread along the IRDC (size >2 pc), showing two different components: one bright and compact arising from three condensations (N, E and S), and the other weak and extended along the filament. While the first component shows broad lines (linewidths of ~4-7 kms-1) in both SiO J=2-1 and SiO J=3-2, the second one is only detected in SiO J=2-1 and has narrow lines (~0.8 kms-1). The maps of CO and its isotopologues show that low-density filaments are intersecting the IRDC and appear to merge toward the densest portion of the cloud. This resembles the molecular structures predicted by flow-driven, shock-induced and magnetically-regulated cloud formation models. As in outflows associated with low-mass star formation, the excitation temperatures and fractional abundances of SiO toward N, E and S, increase with velocity from ~6 to 40 K, and from ~1E-10 to >1E-8 respectively, over a velocity range of ~7 kms-1. Since 8 micron sources, 24 micron sources and/or extended 4.5 micron emission are detected in N, E and S, broad SiO is likely produced in outflows associated with high-mass protostars. The excitation temperatures and fractional abundances of the narrow SiO lines, however, are very low (~9 K and ~1E-11, respectively), and consistent with the processing of interstellar grains by the passage of a shock with vs~12 kms-1. This emission could be generated i) by a large-scale shock, perhaps remnant of the IRDC formation process; ii) by decelerated or recently processed gas in large-scale outflows driven by 8 micron and 24 micron sources; or iii) by an undetected and widespread population of lower mass protostars. High-angular resolution observations are needed to disentangle between these three scenarios.

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.

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 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.

The star Cyg OB2 No. 9 is a well-known non-thermal radio emitter. Recent theoretical work suggests that all such O-stars should be in a binary or a multiple system. However, there is no spectroscopic evidence of a binary component. Re-analysis of radio observations from the VLA of this system over 25 years has revealed that the non-thermal emission varies with a period of 2.35+-0.02 yr. This is interpreted as a strong suggestion of a binary system, with the non-thermal emission arising in a wind-collision region. We derived some preliminary orbital parameters for this putative binary and revised the mass-loss rate of the primary star downward from previous estimates.

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.

HD 167971 is a triple system consisting of a 3.3-day eclipsing binary (O5-8 V + O5-8 V) and an O8 supergiant. It is also a well known non-thermal radio emitter. We observed the radio emission of HD 167971 with the Very Large Array (VLA) and the Australia Telescope Compact Array (ATCA). By combining these data with VLA archive observations we constructed a radio lightcurve covering a 20-yr time-range. We searched for, but failed to find, the 3.3-day spectroscopic period of the binary in the radio data. This could be due to the absence of intrinsic synchrotron radiation at the colliding-wind region between the two components of the eclipsing binary, or due to the large amount of free-free absorption that blocks the synchrotron radiation. We are able to explain many of the observed characteristics of the radio data if the non-thermal emission is produced in a colliding-wind region between the supergiant and the combined winds of the binary. Furthermore, if the system is gravitationally bound, the orbital motion occurs over a period of ~ 20 years, or longer, as suggested by the long-term variability in the radio data. We argue that the variability is due to the free-free absorption that changes with orbital phase or may also in part be due to changes in separation, should the orbit be eccentric.

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 a model for the non-thermal radio emission from presumably single O stars, in terms of synchrotron emission from relativistic electrons accelerated in wind-embedded shocks. These shocks are associated with an unstable, chaotic wind. The main improvement with respect to earlier models is the inclusion of the radial dependence of the shock velocity jump and compression ratio, based on 1D hydrodynamical simulations. The decrease of the velocity jump and the compression ratio as a function of radius produces a rapidly decreasing synchrotron emissivity. This effectively prohibits the models from reproducing the spectral shape of the observed non-thermal radio emission. We investigate a number of "escape routes" by which the hydrodynamical predictions might be reconciled with the radio observations. Although these escape routes reproduce the observed spectral shape, none of these escape routes are physically plausible. In particular, re-acceleration by feeding an electron distribution through a number of shocks, is in contradiction with current hydrodynamical simulations. These hydrodynamical simulations have their limitations, most notably the use of 1D. At present, it is not feasible to perform 2D simulations of the wind out to the distances required for synchrotron-emission models. Based on the current hydrodynamic models, we suspect that the observed non-thermal radio emission from O stars cannot be explained by wind-embedded shocks associated with the instability of the line-driving mechanism. The most likely alternative mechanism is synchrotron emission from colliding winds. That would imply that all O stars with non-thermal radio emission should be members of binary or multiple systems.

We present a radio lightcurve of the O5.5 III(f+) star HD 168112, based on archive data from the Very Large Array (VLA) and the Australia Telescope Compact Array (ATCA). The fluxes show considerable variability and a negative spectral index, thereby confirming that HD 168112 is a non-thermal radio emitter. The non-thermal radio emission is believed to be due to synchrotron radiation from relativistic electrons that have been Fermi accelerated in shocks. For HD 168112, it is not known whether these shocks are due to a wind-wind collision in a binary system or to the intrinsic instability of the stellar wind driving mechanism. Assuming HD 168112 to be a single star, our synchrotron model shows that the velocity jump of the shocks should be very high, or there should be a very large number of shocks in the wind. Neither of these is compatible with time-dependent hydrodynamical calculations of O star winds. If, on the other hand, we assume that HD 168112 is a binary, the high velocity jump is easily explained by ascribing it to the wind-wind collision. By further assuming the star to be an eccentric binary, we can explain the observed radio variability by the colliding-wind region moving in and out of the region where free-free absorption is important. The radio data presented here show that the binary has a period between one and two years. By combining the radio data with X-ray data, we find that the most likely period is about 1.4 yr.

We present a model for the non-thermal radio emission from bright O stars, in terms of synchrotron emission from wind-embedded shocks. The model is an extension of an earlier one, with an improved treatment of the cooling of relativistic electrons. This improvement limits the synchrotron-emitting volume to a series of fairly narrow layers behind the shocks. We show that the width of these layers increases with increasing wavelength, which has important consequences for the shape of the spectrum. We also show that the strongest shocks produce the bulk of the emission, so that the emergent radio flux can be adequately described as coming from a small number of shocks, or even from a single shock. A single shock model is completely determined by four parameters: the position of the shock, the compression ratio and velocity jump of the shock, and the surface magnetic field. Applying a single shock model to the O5 If star Cyg OB2 No. 9 allows a good determination of the compression ratio and shock position and, to a lesser extent, the magnetic field and velocity jump. Our main conclusion is that strong shocks need to survive out to distances of a few hundred stellar radii. Even with multiple shocks, the shocks needed to explain the observed emission are stronger than predictions from time-dependent hydrodynamical simulations.

We report the results of a multiwavelength study of the non-thermal radio emitter HD168112 (O5.5III(f^+)). The detailed analysis of two quasi-simultaneous XMM-Newton and VLA observations reveals strong variability of this star both in the X-ray and radio ranges. The X-ray observations separated by five months reveal a decrease of the X-ray flux of ~30%. The radio emission on the other hand increases by a factor 5-7 between the two observations obtained roughly simultaneously with the XMM-Newton pointings. The X-ray data reveal a hard emission that is most likely produced by a thermal plasma at kT ~2-3 keV while the VLA data confirm the non-thermal status of this star in the radio waveband. Comparison with archive X-ray and radio data confirms the variability of this source in both wavelength ranges over a yet ill defined time scale. The properties of HD168112 in the X-ray and radio domain point towards a binary system with a significant eccentricity and an orbital period of a few years. However, our optical spectra reveal no significant changes of the star's radial velocity suggesting that if HD168112 is indeed a binary, it must be seen under a fairly low inclination.

We present a theoretical model for the non-thermal radio emission from single hot stars, in terms of synchrotron radiation from electrons accelerated in wind-embedded shocks. The model is described by five independent parameters each with a straightforward physical interpretation. Applying the model to a high-quality observation of Cyg OB2 No. 9 (O5 If), we obtain meaningful constraints on most parameters. The most important result is that the outer boundary of the synchrotron emission region must lie between 500 and 2200 stellar radii. This means that shocks must persist up to that distance. We also find that relatively weak shocks (with a compression ratio < 3) are needed to produce the observed radio spectrum. These results are compatible with current hydrodynamical predictions. Most of our models also show a relativistic electron fraction that increases outwards. This points to an increasing efficiency of the acceleration mechanism, perhaps due to multiple acceleration, or an increase in the strength of the shocks. Implications of our results for non-thermal X-ray emission are discussed.