Search

We extend recent numerical results (Dobbs et. al. 2006) on molecular cloud formation in spiral galaxies by including a multi-phase medium. The addition of a hot phase of gas enhances the structure in the cold gas, and significantly increases the fraction of molecular hydrogen that is formed when the cold gas passes through a spiral shock. The difference in structure is reflected in the mass power spectrum of the molecular clouds, which is steeper for the multi-phase calculations. The increase in molecular gas occurs as the addition of a hot phase leads to higher densities in the cold gas. In particular, cold gas is confined in clumps between the spiral arms and retains a higher molecular fraction. Unlike the single phase results, molecular clouds are present in the inter-arm regions for the multi-phase medium. However the density of the inter-arm molecular hydrogen is generally below that which can be reliably determined from CO measurements. We therefore predict that for a multi-phase medium, there will be low density clouds containing cold atomic and molecular hydrogen, which are potentially entering the spiral arms., 11 pages, 14 figures, accepted for publication in MNRAS

The transport of gas towards the centre of galaxies is critical for black hole feeding and, indirectly, it can control active galactic nucleus (AGN) feedback. We have quantified the molecular gas inflow in the central R, 19 pages, 12 figures. Accepted by A&A

The Legacy ExtraGalactic UV Survey (LEGUS) is a Cycle 21 Treasury program on the Hubble Space Telescope, aimed at the investigation of star formation and its relation with galactic environment in nearby galaxies, from the scales of individual stars to those of ~kpc-size clustered structures. Five-band imaging, from the near-ultraviolet to the I-band, with the Wide Field Camera 3, plus parallel optical imaging with the Advanced Camera for Surveys, is being collected for selected pointings of 50 galaxies within the local 12 Mpc. The filters used for the observations with the Wide Field Camera 3 are: F275W(2,704 A), F336W(3,355 A), F438W(4,325 A), F555W(5,308 A), and F814W(8,024 A); the parallel observations with the Advanced Camera for Surveys use the filters: F435W(4,328 A), F606W(5,921 A), and F814W(8,057 A). The multi-band images are yielding accurate recent (, Comment: 67 pages, 9 figures; accepted for publication on the Astronomical Journal

We present CO, H2, H I and HISA (H Iself-absorption) distributions from a set of simulations of\ud grand design spirals including stellar feedback, self-gravity, heating and cooling. We replicate\ud the emission of the second galactic quadrant by placing the observer inside the modelled\ud galaxies and post-process the simulations using a radiative transfer code, so as to create\ud synthetic observations. We compare the synthetic data cubes to observations of the second\ud quadrant of the Milky Way to test the ability of the current models to reproduce the basic\ud chemistry of the Galactic interstellar medium (ISM), as well as to test how sensitive such\ud galaxy models are to different recipes of chemistry and/or feedback. We find that models\ud which include feedback and self-gravity can reproduce the production of CO with respect\ud to H2 as observed in our Galaxy, as well as the distribution of the material perpendicular\ud to the Galactic plane. While changes in the chemistry/feedback recipes do not have a huge\ud impact on the statistical properties of the chemistry in the simulated galaxies, we find that the\ud inclusion of both feedback and self-gravity are crucial ingredients, as our test without feedback\ud failed to reproduce all of the observables. Finally, even though the transition from H2 to CO\ud seems to be robust, we find that all models seem to underproduce molecular gas, and have a\ud lower molecular to atomic gas fraction than is observed. Nevertheless, our fiducial model with\ud feedback and self-gravity has shown to be robust in reproducing the statistical properties of\ud the basic molecular gas components of the ISM in our Galaxy.

To date, most numerical simulations of molecular clouds, and star formation within them, assume a uniform density sphere or box with an imposed turbulent velocity field. In this work, we select molecular clouds from galactic scale simulations as initial conditions, increase their resolution, and re-simulate them using the SPH code Gadget2. Our approach provides clouds with morphologies, internal structures, and kinematics that constitute more consistent and realistic initial conditions for simulations of star formation. We perform comparisons between molecular clouds derived from a galactic simulation, and spheres of turbulent gas of similar dimensions, mass and velocity dispersion. We focus on properties of the clouds such as their density, velocity structure and star formation rate. We find that the inherited velocity structure of the galactic clouds has a significant impact on the star formation rate and evolution of the cloud. Our results indicate that, although we can follow the time evolution of star formation in any simulated cloud, capturing the entire history is difficult as we ignore any star formation that might have occurred before initialisation. Overall, the turbulent spheres do not match the complexity of the galactic clouds., 5 pages, 5 figures, 1 table. Accepted to MNRAS Letters

We present a study of the spatial distribution of the stellar cluster populations in the star forming galaxy NGC 628. Using Hubble Space Telescope broad band WFC3/UVIS UV and optical images from the Treasury Program LEGUS (Legacy ExtraGalactic UV Survey), we have identified 1392 potential young (, Comment: 18 pages, 18 figures, 1 table. Accepted for publication in The Astrophysical Journal. The HST images of Figure 1 are degraded in order to comply with the size requirement of arXiv

We present Smoothed Particle Hydrodynamics (SPH) simulations of molecular cloud formation in spiral galaxies. These simulations model the response of a non-self-gravitating gaseous disk to a galactic potential. The spiral shock induces high densities in the gas, and considerable structure in the spiral arms, which we identify as molecular clouds. We regard the formation of these structures as due to the dynamics of clumpy shocks, which perturb the flow of gas through the spiral arms. In addition, the spiral shocks induce a large velocity dispersion in the spiral arms, comparable with the magnitude of the velocity dispersion observed in molecular clouds. We estimate the formation of molecular hydrogen, by post-processing our results and assuming the gas is isothermal. Provided the gas is cold ($T\le100$ K), the gas is compressed sufficiently in the spiral shock for molecular hydrogen formation to occur in the dense spiral arm clumps. These molecular clouds are largely confined to the spiral arms, since most molecular gas is photodissociated to atomic hydrogen upon leaving the arms., 13 pages, 16 figures, accepted for publication in MNRAS

By resimulating a region of a global disc simulation at higher resolution, we resolve and study the properties of molecular clouds with a range of masses from a few 100's M$_{\odot}$ to $10^6$ M$_{\odot}$. The purpose of our paper is twofold, i) to compare the ISM and GMCs at much higher resolution compared to previous global simulations, and ii) to investigate smaller clouds and characteristics such as the internal properties of GMCs which cannot be resolved in galactic simulations. We confirm the robustness of cloud properties seen in previous galactic simulations, and that these properties extend to lower mass clouds, though we caution that velocity dispersions may not be measured correctly in poorly resolved clouds. We find that the properties of the clouds and ISM are only weakly dependent on the details of local stellar feedback, although stellar feedback is important to produce realistic star formation rates and agreement with the Schmidt-Kennicutt relation. We study internal properties of GMCs resolved by $10^4-10^5$ particles. The clouds are highly structured, but we find clouds have a velocity dispersion radius relationship which overall agrees with the Larson relation. The GMCs show evidence of multiple episodes of star formation, with holes corresponding to previous feedback events and dense regions likely to imminently form stars. Our simulations show clearly long filaments, which are seen predominantly in the inter-arm regions, and shells., 13 pages, 11 figures, accepted for publication in MNRAS, movies at http://empslocal.ex.ac.uk/people/staff/cld214/movies_region.html and http://empslocal.ex.ac.uk/people/staff/cld214/movies_region2.html

We present the results from simulations of GMC formation in spiral galaxies. First we discuss cloud formation by cloud-cloud collisions, and gravitational instabilities, arguing that the former is prevalent at lower galactic surface densities and the latter at higher. Cloud masses are also limited by stellar feedback, which can be effective before clouds reach their maximum mass. We show other properties of clouds in simulations with different levels of feedback. With a moderate level of feedback, properties such as cloud rotations and virial parameters agree with observations. Without feedback, an unrealistic population of overly bound clouds develops. Spiral arms are not found to trigger star formation, they merely gather gas into more massive GMCs. We discuss in more detail interactions of clouds in the ISM, and argue that these are more complex than early ideas of cloud-cloud collisions. Finally we show ongoing work to determine whether the Milky Way is a flocculent or grand design spiral.

Using data from the PdBI Arcsecond Whirlpool Survey (PAWS), we have generated the largest extragalactic Giant Molecular Cloud (GMC) catalog to date, containing 1,507 individual objects. GMCs in the inner M51 disk account for only 54% of the total 12CO(1-0) luminosity of the survey, but on average they exhibit physical properties similar to Galactic GMCs. We do not find a strong correlation between the GMC size and velocity dispersion, and a simple virial analysis suggests that 30% of GMCs in M51 are unbound. We have analyzed the GMC properties within seven dynamically-motivated galactic environments, finding that GMCs in the spiral arms and in the central region are brighter and have higher velocity dispersions than inter-arm clouds. Globally, the GMC mass distribution does not follow a simple power law shape. Instead, we find that the shape of the mass distribution varies with galactic environment: the distribution is steeper in inter-arm region than in the spiral arms, and exhibits a sharp truncation at high masses for the nuclear bar region. We propose that the observed environmental variations in the GMC properties and mass distributions are a consequence of the combined action of large-scale dynamical processes and feedback from high mass star formation. We describe some challenges of using existing GMC identification techniques for decomposing the 12CO(1-0) emission in molecule-rich environments, such as M51's inner disk., 73 pages, 18 figures, 14 tables, accepted for publication in ApJ

Giant molecular clouds (GMCs) are the primary reservoirs of cold, star-forming molecular gas in the Milky Way and similar galaxies, and thus any understanding of star formation must encompass a model for GMC formation, evolution, and destruction. These models are necessarily constrained by measurements of interstellar molecular and atomic gas, and the emergent, newborn stars. Both observations and theory have undergone great advances in recent years, the latter driven largely by improved numerical simulations, and the former by the advent of large-scale surveys with new telescopes and instruments. This chapter offers a thorough review of the current state of the field.

The kinematic complexity and the favorable position of M51 on the sky make this galaxy an ideal target to test different theories of spiral arm dynamics. Taking advantage of the new high resolution PdBI Arcsecond Whirlpool Survey (PAWS) data, we undertake a detailed kinematic study of M51 to characterize and quantify the origin and nature of the non-circular motions. Using a tilted-ring analysis supported by several other archival datasets we update the estimation of M51's position angle (PA=(173 +/- 3) deg) and inclination (i=(22 +/- 5) deg). Harmonic decomposition of the high resolution (40 pc) CO velocity field shows the first kinematic evidence of an m=3 wave in the inner disk of M51 with a corotation at R(CR,m=3)=1.1 +/- 0.1 kpc and a pattern speed of Omega_p(m=3) = 140 km/(s kpc). This mode seems to be excited by the nuclear bar, while the beat frequencies generated by the coupling between the m=3 mode and the main spiral structure confirm its density-wave nature. We observe also a signature of an m=1 mode that is likely responsible for the lopsidedness of M51 at small and large radii. We provide a simple method to estimate the radial variation of the amplitude of the spiral perturbation (Vsp) attributed to the different modes. The main spiral arm structure has =50-70 km/s, while the streaming velocity associated with the m=1 and m=3 modes is, in general, 2 times lower. Our joint analysis of HI and CO velocity fields at low and high spatial resolution reveals that the atomic and molecular gas phases respond differently to the spiral perturbation due to their different vertical distribution and emission morphology., Comment: 42 pages, 12 figures, accepted for publication in ApJ

In this Letter, we analyse the distributions of stellar ages in giant molecular clouds (GMCs) in spiral arms, interarm spurs and at large galactic radii, where the spiral arms are relatively weak. We use the results of numerical simulations of galaxies, which follow the evolution of GMCs and include star particles where star formation events occur. We find that GMCs in spiral arms tend to have predominantly young (

The majority of astrophysics involves the study of spiral galaxies, and stars and planets within them, but how spiral arms in galaxies form and evolve is still a fundamental problem. Major progress in this field was made primarily in the 1960s, and early 1970s, but since then there has been no comprehensive update on the state of the field. In this review, we discuss the progress in theory, and in particular numerical calculations, which unlike in the 1960s and 1970s, are now commonplace, as well as recent observational developments. We set out the current status for different scenarios for spiral arm formation, the nature of the spiral arms they induce, and the consequences for gas dynamics and star formation in different types of spiral galaxies. We argue that, with possible the exception of barred galaxies, spiral arms are transient, recurrent and initiated by swing amplified instabilities in the disc. We suppose that unbarred m = 2 spiral patterns are induced by tidal interactions, and slowly wind up over time. However the mechanism for generating spiral structure does not appear to have significant consequences for star formation in galaxies., Comment: 44 pages, 20 pages, review article accepted for publication in PASA

We have performed the first numerical simulations to resolve the star formation process on sub-parsec scales, whilst also following the dynamics of the interstellar medium (ISM) on galactic scales. The warm low density ISM gas flows into the spiral arms where orbit crowding produces the shock formation of dense clouds, held together temporarily by their external pressure. Cooling allows the gas to be compressed to sufficiently high densities that local regions collapse under their own gravity and form stars. The star formation rates follow a Schmidt-Kennicutt \(\varSigma _{\mathrm{SFR}} \propto \varSigma _{\mathrm{gas}}^{1.4}\) type relation with the local surface density of gas while following a linear relation with the dense cold gas. Cooling is the primary driver of star formation and the star formation rates as it determines the amount of cold gas available for gravitational collapse.

International audience; We use the high spatial and spectral resolution of the PAWS CO(1-0) survey of the inner 9 kpc of the iconic spiral galaxy M51 to examine the effects of gas streaming motions on the star-forming properties of individual giant molecular clouds (GMCs). We compare our view of gas flows in M51--which arise due to departures from axisymmetry in the gravitational potential (i.e., the nuclear bar and spiral arms)--with the global pattern of star formation as traced by Hα and 24 μm emission. We find that the dynamical environment of GMCs strongly affects their ability to form stars, in the sense that GMCs situated in regions with large streaming motions can be stabilized, while similarly massive GMCs in regions without streaming go on to efficiently form stars. We argue that this is the result of reduced surface pressure felt by clouds embedded in an ambient medium undergoing large streaming motions, which prevent collapse. Indeed, the variation in gas depletion time expected based on the observed streaming motions throughout the disk of M51 quantitatively agrees with the variation in the observed gas depletion time scale. The example of M51 shows that streaming motions, triggered by gravitational instabilities in the form of bars and spiral arms, can alter the star formation law; this can explain the variation in gas depletion time among galaxies with different masses and morphologies. In particular, we can explain the long gas depletion times in spiral galaxies compared with dwarf galaxies and starbursts. We suggest that adding a dynamical pressure term to the canonical free-fall time produces a single star formation law that can be applied to all star-forming regions and galaxies across cosmic time.

We analyse the distribution of CO brightness temperature and integrated intensity in M51 at ~40 pc resolution using new CO data from the Plateau de Bure Arcsecond Whirlpool Survey (PAWS). We present probability distribution functions (PDFs) of the CO emission within the PAWS field, which covers the inner 11 x 7 kpc of M51. We find variations in the shape of CO PDFs within different M51 environments, and between M51 and M33 and the Large Magellanic Cloud (LMC). Globally, the PDFs for the inner disk of M51 can be represented by narrow lognormal functions that cover 1 to 2 orders of magnitude in CO brightness and integrated intensity. The PDFs for M33 and the LMC are narrower and peak at lower CO intensities. However, the CO PDFs for different dynamical environments within the PAWS field depart from the shape of the global distribution. The PDFs for the interarm region are approximately lognormal, but in the spiral arms and central region of M51, they exhibit diverse shapes with a significant excess of bright CO emission. The observed environmental dependence of the shape of the CO PDFs is qualitatively consistent with changes that would be expected if molecular gas in the spiral arms has a larger range of average densities, gas temperatures and velocity fluctuations, though further work is required to disentangle the importance of large-scale dynamical effects versus star formation feedback in regulating these properties. We show that the shape of the CO PDFs for different M51 environments is only weakly related to global properties of the CO emission, but is strongly correlated with some properties of the local giant molecular cloud (GMC) and young stellar cluster populations. For galaxies with strong spiral structure such as M51, our results indicate that galactic-scale dynamical processes play a significant role in the formation and evolution of GMCs and stellar clusters.(abridged), 30 pages, 18 figures, accepted to ApJ PAWS special issue

We compare the properties of giant molecular clouds (GMCs) in M51 identified by the Plateau de Bure Interferometer Whirlpool Arcsecond Survey (PAWS) with GMCs identified in wide-field, high resolution surveys of CO emission in M33 and the Large Magellanic Cloud (LMC). We find that GMCs in M51 are larger, brighter and have higher velocity dispersions relative to their size than equivalent structures in M33 and the LMC. These differences imply that there are genuine variations in the average mass surface density of the different GMC populations. To explain this, we propose that the pressure in the interstellar medium surrounding the GMCs plays a role in regulating their density and velocity dispersion. We find no evidence for a correlation between size and linewidth in any of M51, M33 or the LMC when the CO emission is decomposed into GMCs, although moderately robust correlations are apparent when regions of contiguous CO emission (with no size limitation) are used. Our work demonstrates that observational bias remains an important obstacle to the identification and study of extragalactic GMC populations using CO emission, especially in molecule-rich galactic environments., Comment: 25 pages, 11 figures, accepted for publication in ApJ. Uses emulateapj LaTeX macros. For more information on PAWS, further papers and data, see http://www.mpia.de/PAWS/

The PdBI (Plateau de Bure Interferometer) Arcsecond Whirlpool Survey (PAWS) has mapped the molecular gas in the central ~9kpc of M51 in its 12CO(1-0) line emission at cloud-scale resolution of ~40pc using both IRAM telescopes. We utilize this dataset to quantitatively characterize the relation of molecular gas (or CO emission) to other tracers of the interstellar medium (ISM), star formation and stellar populations of varying ages. Using 2-dimensional maps, a polar cross-correlation technique and pixel-by-pixel diagrams, we find: (a) that (as expected) the distribution of the molecular gas can be linked to different components of the gravitational potential, (b) evidence for a physical link between CO line emission and radio continuum that seems not to be caused by massive stars, but rather depend on the gas density, (c) a close spatial relation between the PAH and molecular gas emission, but no predictive power of PAH emission for the molecular gas mass,(d) that the I-H color map is an excellent predictor of the distribution (and to a lesser degree the brightness) of CO emission, and (e) that the impact of massive (UV-intense) young star-forming regions on the bulk of the molecular gas in central ~9kpc can not be significant due to a complex spatial relation between molecular gas and star-forming regions that ranges from co-spatial to spatially offset to absent. The last point, in particular, highlights the importance of galactic environment -- and thus the underlying gravitational potential -- for the distribution of molecular gas and star formation., 52 pages, 14 figures, 2 tables, (several minor typos corrected) accepted by ApJ, high resolution version available, see http://www.mpia.de/PAWS/pub/paws_schinnerer.pdf ; for more information on PAWS, further papers and the data, see http://www.mpia.de/PAWS

Understanding star formation is problematic as it originates in the large scale dynamics of a galaxy but occurs on the small scale of an individual star forming event. This paper presents the first numerical simulations to resolve the star formation process on sub-parsec scales, whilst also following the dynamics of the interstellar medium (ISM) on galactic scales. In these models, the warm low density ISM gas flows into the spiral arms where orbit crowding produces the shock formation of dense clouds, held together temporarily by their external pressure. Cooling allows the gas to be compressed to sufficiently high densities that local regions collapse under their own gravity and form stars. The star formation rates follow a Schmidt-Kennicutt \Sigma_{SFR} ~ \Sigma_{gas}^{1.4} type relation with the local surface density of gas while following a linear relation with the cold and dense gas. Cooling is the primary driver of star formation and the star formation rates as it determines the amount of cold gas available for gravitational collapse. The star formation rates found in the simulations are offset to higher values relative to the extragalactic values, implying a constant reduction, such as from feedback or magnetic fields, is likely to be required. Intriguingly, it appears that a spiral or other convergent shock and the accompanying thermal instability can explain how star formation is triggered, generate the physical conditions of molecular clouds and explain why star formation rates are tightly correlated to the gas properties of galaxies., Comment: 13 pages, 12 figures. MNRAS in press

We present a detailed study of the evolution of GMCs in a galactic disc simulation. We follow individual GMCs (defined in our simulations by a total column density criterion), including their level of star formation, from their formation to dispersal. We find the evolution of GMCs is highly complex. GMCs often form from a combination of smaller clouds and ambient ISM, and similarly disperse by splitting into a number of smaller clouds and ambient ISM. However some clouds emerge as the result of the disruption of a more massive GMC, rather than from the assembly of smaller clouds. Likewise in some cases, clouds accrete onto more massive clouds rather than disperse. Because of the difficulty of determining a precursor or successor of a given GMC, determining GMC histories and lifetimes is highly non-trivial. Using a definition relating to the continuous evolution of a cloud, we obtain lifetimes typically of 4-25 Myr for >10^5 M$_{\odot}$ GMCs, over which time the star formation efficiency is about 1 %. We also relate the lifetime of GMCs to their crossing time. We find that the crossing time is a reasonable measure of the actual lifetime of the cloud, although there is considerable scatter. The scatter is found to be unavoidable because of the complex and varied shapes and dynamics of the clouds. We study cloud dispersal in detail and find both stellar feedback and shear contribute to cloud disruption. We also demonstrate that GMCs do not behave as ridge clouds, rather massive spiral arm GMCs evolve into smaller clouds in inter-arm spurs., 15 pages, 16 figures, accepted for publication in MNRAS

We present the current standing of an investigation into the structure of the Milky Way. We use smoothed particle hydrodynamics (SPH) to simulate the ISM gas in the Milky Way under the effect of a number of different gravitational potentials representing the spiral arms and nuclear bars, both fixed and time-dependent. The gas is subject to ISM cooling and chemistry, enabling us to track the CO and HI density. We use a 3D grid-based radiative transfer code to simulate the emission from the SPH output, allowing for the construction of synthetic longitude-velocity maps as viewed from the Earth. By comparing these maps with the observed emission in CO and HI from the Milky Way (Dame et al. 2001, Kalberla et al. 2005), we can infer the arm/bar geometry that provides a best fit to our Galaxy. By doing so we aim to answer key questions concerning the morphology of the Milky Way such as the number of the spiral arms, the pattern speeds of the bar(s) and arms, the pitch angle of the arms and shape of the bar(s), Comment: 6 pages, 5 figures, contributed talk, to appear in Proceedings of the IAU Symposium No. 298, "Setting the scene for Gaia and LAMOST"

The Serpens North Cluster is a nearby low mass star forming region which is part of the Gould Belt. It contains a range of young stars thought to correspond to two different bursts of star formation and so provides the opportunity to study different stages of cluster formation. Here we present recent results based on a detailed analysis of molecular line and dust continuum data of this region. By studying the dynamics, column density and dust temperature structure of Serpens we propose that a cloud-cloud collision has been responsible for the trigger of the most recent episode of star formation. This scenario has been explored using SPH simulations. These simulations are able to successfully explain the similarities and differences observed in the two sub-protoclusters composing Serpens North.

The complexity of the ISM is such that it is unlikely that star formation is initiated in the same way in all molecular clouds. While some clouds seem to collapse on their own, others may be triggered by an external event such as a cloud/flow collision forming a gravitationally unstable enhanced density layer. This work tests cloud-cloud collisions as the triggering mechanism for star formation in the Serpens Main Cluster as has been suggested by previous work. A set of smoothed particle hydrodynamics (SPH) simulations of the collision between two cylindrical clouds are performed and compared to (sub)millimetre observations of the Serpens Main Cluster. A configuration has been found which reproduces many of the observed characteristics of Serpens, including some of the main features of the peculiar velocity field. The evolution of the velocity with position throughout the model is similar to that observed and the column density and masses within the modeled cloud agree with those measured for the SE sub-cluster. Furthermore, our results also show that an asymmetric collision provides the ingredients to reproduce lower density filaments perpendicular to the main structure, similar to those observed. In this scenario, the formation of the NW sub-cluster of Serpens can be reproduced only if there is a pre-existing marginally gravitationally unstable region at the time the collision occurs. This work supports the interpretation that a collision between two clouds may have been the trigger of the most recent burst of star formation in Serpens. It not only explains the complicated velocity structure seen in the region, but also the temperature differences between the north (in "isolated" collapse) and the south (resulting from the shock between the clouds). In addition it provides an explanation for the sources in the south having a larger spread in age than those in the north., Comment: Accepted for publication in A&A

Using the Torus radiative transfer code we produce synthetic observations of the 21 cm neutral hydrogen line from an SPH simulation of a spiral galaxy. The SPH representation of the galaxy is mapped onto an AMR grid, and a ray tracing method is used to calculate 21 cm line emission for lines of sight through the AMR grid in different velocity channels and spatial pixels. The result is a synthetic spectral cube which can be directly compared to real observations. We compare our synthetic spectral cubes to observations of M31 and M33 and find good agreement, whereby increasing velocity channels trace the main disc of the galaxy. The synthetic data also show kinks in the velocity across the spiral arms, evidence of non-circular velocities. These are still present even when we blur our data to a similar resolution as the observations, but largely absent in M31 and M33, indicating those galaxies do not contain significant spiral shocks. Thus the detailed velocity structure of our maps better represent previous observations of the grand design spiral M81.

Using the radiative transfer code Torus, we produce spectral-line cubes of the predicted HI profile from global SPH simulations of spiral galaxies. Torus grids the SPH galaxy using Adaptive Mesh Refinement, then applies a ray-tracing method to infer the HI profile along the line(s) of sight. The gridded galaxy can be observed from any direction, which enables us to model the observed HI profile for galaxies of any orientation. We can also place the observer inside the galaxy, to simulate HI observations taken from the Earth's position in the Milky Way.

I review the progress of SPH calculations for modelling galaxies, and resolving gas dynamics on GMC scales. SPH calculations first investigated the response of isothermal gas to a spiral potential, in the absence of self gravity and magnetic fields. Surprisingly though, even these simple calculations displayed substructure along the spiral arms. Numerical tests indicate that this substructure is still present at high resolution (100 million particles, ~10 pc), and is independent of the initial particle distribution. One interpretation of the formation of substructure is that smaller clouds can agglomerate into more massive GMCs via dissipative collisions. More recent calculations have investigated how other processes, such as the thermodynamics of the ISM, and self gravity affect this simple picture. Further research has focused on developing models with a more realistic spiral structure, either by including stars, or incorporating a tidal interaction., Comment: 8 pages, 6 figures, review talk for IAU 270 symposium 'Computational star formation'

We investigate the relationship between the star formation rate per unit area and the surface density of the ISM (the local Kennicutt-Schmitt law) using a simplified model of the ISM and a simple estimate of the star formation rate based on the mass of gas in bound clumps, the local dynamical timescales of the clumps, and an efficiency parameter of around 5 per cent. Despite the simplicity of the approach, we are able to reproduce the observed linear relation between star formation rate and surface density of dense (molecular) gas. We use a simple model for the dependence of H_2 fraction on total surface density to argue why neither total surface density nor the HI surface density are good local indicators of star formation rate. We also investigate the dependence of the star formation rate on the depth of the spiral potential. Our model indicates that the mean star formation rate does not depend significantly on the strength of the spiral potential, but that a stronger spiral potential, for a given mean surface density, does result in more of the star formation occurring close to the spiral arms. This agrees with the observation that grand design galaxies do not appear to show a larger degree of star formation compared to their flocculent counterparts., 11 pages, 11 figures, accepted by MNRAS

We investigate the formation of GMCs in spiral galaxies through both agglomeration of clouds in the spiral arms, and self gravity. The simulations presented include two-fluid models, which contain both cold and warm gas, although there is no heating or cooling between them. We find agglomeration is predominant when both the warm and cold components of the ISM are effectively stable to gravitational instabilities. In this case, the spacing (and consequently mass) of clouds and spurs along the spiral arms is determined by the orbits of the gas particles and correlates with their epicyclic radii (or equivalently spiral shock strength). Notably GMCs formed primarily by agglomeration tend to be unbound associations of many smaller clouds, which disperse upon leaving the spiral arms. These GMCs are likely to be more massive in galaxies with stronger spiral shocks or higher surface densities. GMCs formed by agglomeration are also found to exhibit both prograde and retrograde rotation, a consequence of the clumpiness of the gas. At higher surface densities, self gravity becomes more important in arranging both the warm and cold gas into clouds and spurs, and determining the properties of the most massive GMCs. These massive GMCs can be distinguished by their higher angular momentum, exhibit prograde rotation and are more bound. For a 20 M$_{\odot}$ pc$^{-2}$ disc, the spacing between the GMCs fits both the agglomeration and self gravity scenarios, as the maximum unstable wavelength of gravitational perturbations in the warm gas is similar to the spacing found when GMCs form solely by agglomeration., 16 pages, 12 figures, accepted for publication in MNRAS. Version with higher resolution figures available at http://www.astro.ex.ac.uk/people/dobbs/Dobbsgr.pdf

We describe simulations of the response of a gaseous disc to an active spiral potential. The potential is derived from an N-body calculation and leads to a multi-armed time-evolving pattern. The gas forms long spiral arms typical of grand design galaxies, although the spiral pattern is asymmetric. The primary difference from a grand-design spiral galaxy, which has a consistent 2/4-armed pattern, is that instead of passing through the spiral arms, gas generally falls into a developing potential minimum and is released only when the local minimum dissolves. In this case, the densest gas is coincident with the spiral potential, rather than offset as in the grand-design spirals. We would there fore expect no offset between the spiral shock and star formation, and no obvious co-rotation radius. Spurs which occur in grand-design spirals when large clumps are sheared off leaving the spiral arms, are rare in the active, time-evolving spiral reported here. Instead, large branches are formed from spiral arms when the underlying spiral potential is dissolving due to the N-body dynamics. We find that the molecular cloud mass spectrum for the active potential is similar to that for clouds in grand design calculations, depending primarily on the ambient pressure rather than the nature of the potential. The largest molecular clouds occur when spiral arms collide, rather than by agglomeration within a spiral arm., 11 pages, 7 figures, accepted for publication in MNRAS

We investigate the thermodynamics of the ISM and the formation of molecular hydrogen through numerical simulations of spiral galaxies. The model follows the chemical, thermal and dynamical response of the disc to an external spiral potential. Self-gravity and magnetic fields are not included. The calculations demonstrate that gas can cool rapidly when subject to a spiral shock. Molecular clouds in the spiral arms arise through a combination of compression of the ISM by the spiral shock and orbit crowding. These results highlight that local self-gravity is not required to form molecular clouds. Self-shielding provides a sharp transition density, below which gas is essentially atomic, and above which the molecular gas fraction is >0.001. The timescale for gas to move between these regimes is very rapid (50 %). Although our calculations are unable to resolve turbulent motions on scales smaller than the spiral arm and do not include self-gravity. True cloud formation timescales are therefore expected to be even shorter. The mass budget of the disc is dominated by cold gas residing in the spiral arms. Between 50 and 75 % of this gas is in the atomic phase. When this gas leaves the spiral arm and drops below the self-shielding limit it is heated by the galactic radiation field. Consequently, most of the volume in the interarm regions is filled with warm atomic gas. However, some cold spurs and clumps can survive in interarm regions for periods comparable to the interarm passage timescale. Altogether between 7 and 40% of the gas in our disc is molecular, depending on the surface density of the calculation, with approximately 20% molecular for a surface density comparable to the solar neighbourhood., Comment: 16 pages, 19 figures, accepted for publication in MNRAS

We present ongoing hydrodynamic and MHD simulations of molecular cloud formation in spiral galaxies. The hydrodynamic results show the formation of molecular gas clouds where spiral shocks compress atomic gas to high densities. The spiral shocks also produce structure in the spiral arms, provided the gas is cold (< 1000 K). When both hot and cold components of the ISM are modeled, this structure is enhanced. Properties such as the clump mass spectra and spatial distribution will be compared from clouds identified in these simulations. In particular the multiphase simulations predict the presence of much more interarm molecular gas than when a single phase is assumed. We also discuss very recent results from galactic-scale MHD calculations. From observational comparisons of the magnetic and thermal pressure, magnetic fields are expected to be a major factor in explaining the dynamics of the ISM, from kpc scales to those of star formation. We describe the difference in structure of the spiral arms, and the evolution of the global magnetic field for a range of field strengths., To appear in the proceedings for "Pathways Through an Eclectic Universe", Editors: Johan Knapen, Terry Mahoney, and Alexandre Vazdekis

We present an alternative explanation for the nature of turbulence in molecular clouds. Often associated with classical models of turbulence, we instead interpret the observed gas dynamics as random motions, induced when clumpy gas is subject to a shock. From simulations of shocks, we show that a supersonic velocity dispersion occurs in the shocked gas provided the initial distribution of gas is sufficiently non-uniform. We investigate the velocity size-scale relation $\sigma \propto r^{\alpha}$ for simulations of clumpy and fractal gas, and show that clumpy shocks can produce realistic velocity size-scale relations with mean $\alpha \thicksim 0.5$. For a fractal distribution, with a fractal dimension of 2.2 similar to what is observed in the ISM, we find $\sigma \propto r^{0.4}$. The form of the velocity size-scale relation can be understood as due to mass loading, i.e. the post-shock velocity of the gas is determined by the amount of mass encountered as the gas enters the shock. We support this hypothesis with analytical calculations of the velocity dispersion relation for different initial distributions. A prediction of this model is that the line-of sight velocity dispersion should depend on the angle at which the shocked gas is viewed., Comment: 11 pages, 17 figures, accepted for publication in MNRAS

We present Smoothed Particle Hydrodynamics (SPH) simulations of the response of gas discs to a spiral potential. These simulations show that the commonly observed spurs and feathering in spiral galaxies can be understood as being due to structures present in the spiral arms that are sheared by the divergent orbits in a spiral potential. Thus, dense molecular cloud-like structures generate the perpendicular spurs as they leave the spiral arms. Subsequent feathering occurs as spurs are further sheared into weaker parallel structures as they approach the next spiral passage. Self-gravity of the gas is not included in these simulations, stressing that these features are purely due to the hydrodynamics in spiral shocks. Instead, a necessary condition for this mechanism to work is that the gas need be relatively cold (1000 K or less) in order that the shock is sufficient to generate structure in the spiral arms, and such structure is not subsequently smoothed by the gas pressure., Comment: 6 pages, 6 figures, accepted for publication in MNRAS

We present numerical simulations of the passage of clumpy gas through a galactic spiral shock and the subsequent formation of giant molecular clouds (GMCs). The spiral shock forms dense clouds while dissipating kinetic energy, producing regions that are locally gravitationally bound and collapse to form stars. The effect of the clumpiness of gas as it passes through the shock is to generate chaotic internal motions in the gas. The kinematics of these motions are found to agree with the observed velocity-dispersion/size relation found in star-forming regions. In contrast to the standard picture where continuously driven turbulence generates the density inhomogeneities in star-forming clouds, we find here that it is the clumpiness of the interstellar gas that produces the chaotic motions as it passes through the spiral shock and initiates the star formation process. The velocity dispersion can be understood as being due to the random mass loading of clumps as they converge in the spiral shock. In this model there is no need for any internal or external continuous driving mechanism for the 'turbulence'. The coupling of the clouds' internal kinematics to their externally triggered formation removes the need for the clouds to be self-gravitating. Indeed, while clearly some parts of the clouds are self-gravitating and able to form stars, most of the molecular material remains gravitationally unbound. This can provide a simple explanation for the low efficiency of star formation., 9 pages, 6 figures, accepted for publication in MNRAS

We present numerical investigations into the formation of massive stars from centrally condensed turbulent cores. The results of five hydrodynamical simulations are described, following the collapse of the core, fragmentation and the formation of small clusters of protostars. We generate two different initial turbulent velocity fields corresponding to power-law spectra with power varying as k^{-4} and k^{-3.5}, and apply two different initial core radii. Calculations are included for both completely isothermal collapse, and a non-isothermal equation of state above a critical density 10^{-14}gcm^{-3}. Our calculations reveal the preference of fragmentation over monolithic star formation in turbulent cores. Fragmentation was prevalent in all the isothermal cases. Although disc fragmentation was largely suppressed in the non-isothermal runs due to the small dynamic range between the initial density and the critical density, our results show that some fragmentation still persisted. This is inconsistent with previous suggestions that turbulent cores result in the formation of a single massive star. We conclude that turbulence cannot be measured as an isotropic pressure term., 8 pages, 5 figures, accepted for publication in MNRAS

We present the data of the Plateau de Bure Arcsecond Whirlpool Survey (PAWS), a high spatial and spectral resolution 12CO(1-0) line survey of the inner 10 x 6 kpc of the M51 system, and the first wide-field imaging of molecular gas in a star-forming spiral galaxy with resolution matched to the typical size of Giant Molecular Clouds (40 pc). We describe the observation, reduction, and combination of the Plateau de Bure Interferometer (PdBI) and IRAM-30m "short spacing" data. The final data cube attains 1.1"-resolution over the 270" x 170" field of view, with sensitivity to all spatial scales from the combination of PdBI and IRAM-30m data, and brightness sensitivity of 0.4 K (1 sigma) in each 5 km/s-wide channel map. We find a CO-luminosity of 9.10^8 K km/s pc^2, corresponding to a molecular gas mass of 4.10^9 Msol for a standard CO-to-H2 conversion factor. Unexpectedly, we find that a large fraction, (50 +/- 10)%, of this emission arises mostly from spatial scales larger than 36" or about 1.3 kpc. Through a series of tests, we demonstrate that this extended emission does not result from a processing artifact. We discuss its origin in light of the stellar component, the 12CO/13CO ratio, and the difference between the kinematics and structure of the PdBI-only and hybrid synthesis (PdBI + IRAM-30m) images. The extended emission is consistent with a thick, diffuse disk of molecular gas with a typical scale height of about 200 pc, substructured in unresolved filaments which fills about 0.1% of the volume., 42 pages, 35 PostScript figures. Accepted for publication in ApJ. Uses emulateapj LaTeX macros

Stars are forming in our galaxy at a rate of between 1 and 4 solar masses of stars per year. In contrast to elliptical galaxies, which are largely devoid of star formation, star formation is still going on in spiral galaxies because of their reservoirs of molecular gas, the fuel for new stars. The discs of spiral galaxies are comprised not only of stars as we clearly see from Earth, but also gas (the interstellar medium, ISM). This is where this gas accumulates into cold, dense, molecular regions known as molecular clouds, in which new stars are formed. Most star formation occurs in massive molecular clouds, known as giant molecular clouds (GMCs). However, while we have a good understanding of how individual stars form, there is less consensus on how their natal clouds of gas accumulate, how long these clouds last, how star formation progresses over their lifetime, and indeed how star formation has progressed over the lifetime of the Milky Way., Comment: 7 pages, 6 figures, article which appeared in October edition of Astronomy & Geophysics

We test whether the spiral patterns apparent in many large disk galaxies should be thought of as dynamical features that are stationary in a co-rotating frame for > t_{dyn}, as implied by the density wave approach for explaining spiral arms. If such spiral arms have enhanced star formation (SF), observational tracers for different stages of the SF sequence should show a spatial ordering, from up-stream to downstream in the corotating frame: dense HI, CO, tracing molecular hydrogen gas, 24 micron emission tracing enshrouded SF and UV emission tracing unobscured young stars. We argue that such a spatial ordering should be reflected in the angular cross-correlation (CC, in polar coordinates) using all azimuthal positions among pairs of these tracers; the peak of the CC should be offset from zero, in different directions inside and outside the corotation radius. Recent spiral SF simulations by Dobbs & Pringle, show explicitly that for the case of a stationary spiral arm potential such angular offsets between gas and young stars of differing ages should be observable as cross-correlation offsets. We calculate the angular cross-correlations for different observational SF sequence tracers in 12 nearby spiral galaxies, drawing on a data set with high quality maps of the neutral gas HI, THINGS), molecular gas (CO, HERACLES) along with 24 micron emission (Spitzer, SINGS); we include FUV images (GALEX) and 3.6 $��$m emission (Spitzer, IRAC) for some galaxies, tracing aging stars and longer timescales. In none of the resulting tracer cross-correlations for this sample do we find systematic angular offsets, which would be expected for a stationary dynamical spiral pattern of well-defined pattern speed. This result indicates that spiral density waves in their simplest form are not an important aspect of explaining spirals in large disk galaxies., 13 pages, 16 figures

We have performed simulations of cluster formation along two regions of a spiral arm taken from a global Milky Way simulation, including photoionising feedback. One region is characterised by strongly converging flows, the other represents a more typical spiral arm region. We find that more massive clusters are able to form on shorter timescales for the region with strongly converging flows. Mergers between clusters are frequent in the case of the strongly converging flows and enable the formation of massive clusters. We compare equivalent clusters formed in simulations with and without ionisation. Photoionisation does not prevent massive cluster formation, but can be seen to limit the masses of the clusters. On average the mass is reduced by around 20%, but we see a large spread from ionisation having minimal difference to leading to a 50% reduction in mass. Photoionisation is also able to clear out the gas in the vicinity of the clusters on Myr timescales, which can produce clusters with larger radii that are surrounded by more massive stellar halos. We find that the ionising feedback has more impact in our second region which is less dense and has less strongly converging flows., Comment: 20 pages, 18 figures, accepted for publication in MNRAS

We describe and execute a novel approach to observationally estimate the lifetimes of giant molecular clouds (GMCs). We focus on the cloud population between the two main spiral arms in M51 (the inter-arm region) where cloud destruction via shear and star formation feedback dominates over formation processes. By monitoring the change in GMC number densities and properties from one side of the inter-arm to the other, we estimate the lifetime as a fraction of the inter-arm travel time. We find that GMC lifetimes in M51's inter-arm are finite and short, 20 to 30 Myr. Such short lifetimes suggest that cloud evolution is influenced by environment, in which processes can disrupt GMCs after a few free-fall times. Over most of the region under investigation shear appears to regulate the lifetime. As the shear timescale increases with galactocentric radius, we expect cloud destruction to switch primarily to star formation feedback at larger radii. We identify a transition from shear- to feedback-dominated disruption through a change in the behavior of the GMC number density. The signature suggests that shear is more efficient at completely dispersing clouds, whereas feedback transforms the population, e.g. by fragmenting high mass clouds into lower mass pieces. Compared to the characteristic timescale for molecular hydrogen in M51, our short lifetimes suggest that gas can remain molecular while clouds disperse and reassemble. We propose that galaxy dynamics regulates the cycling of molecular material from diffuse to bound (and ultimately star-forming) objects, contributing to long observed molecular depletion times in normal disk galaxies. We also speculate that, in more extreme environments such as elliptical galaxies and concentrated galaxy centers, star formation can be suppressed when the shear timescale becomes so short that some clouds can not survive to collapse and form stars., Comment: 16 pages, 8 figures, accepted for publication in ApJ