We present simulations of a 500 pc$^2$ region, containing gas of mass 4 $\times$ 10$^6$ M$_\odot$, extracted from an entire spiral galaxy simulation, scaled up in resolution, including photoionising feedback from stars of mass > 18 M$_\odot$. Our region is evolved for 10 Myr and shows clustered star formation along the arm generating $\approx$ 5000 cluster sink particles $\approx$ 5% of which contain at least one of the $\approx$ 4000 stars of mass > 18 M$_\odot$. Photoionisation has a noticeable effect on the gas in the region, producing ionised cavities and leading to dense features at the edge of the HII regions. Compared to the no-feedback case, photoionisation produces a larger total mass of clouds and clumps, with around twice as many such objects, which are individually smaller and more broken up. After this we see a rapid decrease in the total mass in clouds and the number of clouds. Unlike studies of isolated clouds, our simulations follow the long range effects of ionisation, with some already-dense gas becoming compressed from multiple sides by neighbouring HII regions. This causes star formation that is both accelerated and partially displaced throughout the spiral arm with up to 30% of our cluster sink particle mass forming at distances > 5 pc from sites of sink formation in the absence of feedback. At later times, the star formation rate decreases to below that of the no-feedback case.

Young massive clusters (YMCs) are the most intense regions of star formation in galaxies. Formulating a model for YMC formation whilst at the same time meeting the constraints from observations is highly challenging however. We show that forming YMCs requires clouds with densities $\gtrsim$ 100 cm$^{-3}$ to collide with high velocities ($\gtrsim$ 20 km s$^{-1}$). We present the first simulations which, starting from moderate cloud densities of $\sim100$ cm$^{-3}$, are able to convert a large amount of mass into stars over a time period of around 1 Myr, to produce dense massive clusters similar to those observed. Such conditions are commonplace in more extreme environments, where YMCs are common, but atypical for our Galaxy, where YMCs are rare.

In spiral galaxies, the pitch angle, $\alpha$, of the spiral arms is often proposed as a discriminator between theories for the formation of the spiral structure. In Lin-Shu density wave theory, $\alpha$ stays constant in time, being simply a property of the underlying galaxy. In other theories (e.g tidal interaction, self-gravity) it is expected that the arms wind up in time, so that to a first approximation $\cot \alpha \propto t$. For these theories, it would be expected that a sample of galaxies observed at random times should show a uniform distribution of $\cot \alpha$. We show that a recent set of measurements of spiral pitch angles (Yu & Ho 2018) is broadly consistent with this expectation.

We compare the properties of clouds in simulated M33 galaxies to those observed in the real M33. We apply a friends of friends algorithm and CPROPS to identify clouds, as well as a pixel by pixel analysis. We obtain very good agreement between the number of clouds, and maximum mass of clouds. Both are lower than occurs for a Milky Way-type galaxy and thus are a function of the surface density, size and galactic potential of M33. We reproduce the observed dependence of molecular cloud properties on radius in the simulations, and find this is due to the variation in gas surface density with radius. The cloud spectra also show good agreement between the simulations and observations, but the exact slope and shape of the spectra depends on the algorithm used to find clouds, and the range of cloud masses included when fitting the slope. Properties such as cloud angular momentum, velocity dispersions and virial relation are also in good agreement between the simulations and observations, but do not necessarily distinguish between simulations of M33 and other galaxy simulations. Our results are not strongly dependent on the level of feedback used here (10 and 20%) although they suggest that 15% feedback efficiency may be optimal. Overall our results suggest that the molecular cloud properties are primarily dependent on the gas and mass surface density, and less dependent on the localised physics such as the details of stellar feedback, or the numerical code used.

We present the star cluster catalogs for 17 dwarf and irregular galaxies in the $HST$ Treasury Program "Legacy ExtraGalactic UV Survey" (LEGUS). Cluster identification and photometry in this subsample are similar to that of the entire LEGUS sample, but special methods were developed to provide robust catalogs with accurate fluxes due to low cluster statistics. The colors and ages are largely consistent for two widely used aperture corrections, but a significant fraction of the clusters are more compact than the average training cluster. However, the ensemble luminosity, mass, and age distributions are consistent suggesting that the systematics between the two methods are less than the random errors. When compared with the clusters from previous dwarf galaxy samples, we find that the LEGUS catalogs are more complete and provide more accurate total fluxes. Combining all clusters into a composite dwarf galaxy, we find that the luminosity and mass functions can be described by a power law with the canonical index of $-2$ independent of age and global SFR binning. The age distribution declines as a power law, with an index of $\approx-0.80\pm0.15$, independent of cluster mass and global SFR binning. This decline of clusters is dominated by cluster disruption since the combined star formation histories and integrated-light SFRs are both approximately constant over the last few hundred Myr. Finally, we find little evidence for an upper-mass cutoff ($<2\sigma$) in the composite cluster mass function, and can rule out a truncation mass below $\approx10^{4.5}$M$_{\odot}$ but cannot rule out the existence of a truncation at higher masses.

We investigate the future evolution of the Smith Cloud by performing hydrodynamical simulations of the cloud impact onto the gaseous Milky Way Galactic disk. We assume a local origin for the cloud and thus do not include a dark matter component to stabilize it. Our main focus is the cloud's influence on the local and global star formation rate (SFR) of the Galaxy and whether or not it leads to an observable event in the far future. Our model assumes two extremes for the mass of the Smith Cloud, an upper mass limit of 10$^7$ M$_{\odot}$ and a lower mass limit of 10$^6$ M$_{\odot}$, compared to the observational value of a few 10$^6$ M$_{\odot}$. In addition, we also make the conservative assumption that the entirety of the cloud mass of the extended Smith Cloud is concentrated within the tip of the cloud. We find that the impact of the low-mass cloud produces no noticeable change in neither the global SFR nor the local SFR at the cloud impact site within the galactic disk. For the high-mass cloud we find a short-term (roughly 5 Myr) increase of the global SFR of up to 1 M$_{\odot}$ yr$^{-1}$, which nearly doubles the normal Milky Way SFR. This highly localized starburst should be observable.

We present a study correlating the spatial locations of young star clusters with those of molecular clouds in NGC~5194, in order to investigate the timescale over which clusters separate from their birth clouds. The star cluster catalogues are from the Legacy ExtraGalactic UV Survey (LEGUS) and the molecular clouds from the Plateau de Bure Interefrometer Arcsecond Whirpool Survey (PAWS). We find that younger star clusters are spatially closer to molecular clouds than older star clusters. The median ages for clusters associated with clouds is 4~Myr whereas it is 50~Myr for clusters that are sufficiently separated from a molecular cloud to be considered unassociated. After $\sim$6~Myr, the majority of the star clusters lose association with their molecular gas. Younger star clusters are also preferentially located in stellar spiral arms where they are hierarchically distributed in kpc-size regions for 50-100~Myr before dispersing. The youngest star clusters are more strongly clustered, yielding a two-point correlation function with $\alpha=-0.28\pm0.04$, than the GMCs ($\alpha=-0.09\pm0.03$) within the same PAWS field. However, the clustering strength of the most massive GMCs, supposedly the progenitors of the young clusters for a star formation efficiency of a few percent, is comparable ($\alpha=-0.35\pm0.05$) to that of the clusters. We find a galactocentric-dependence for the coherence of star formation, in which clusters located in the inner region of the galaxy reside in smaller star-forming complexes and display more homogeneous distributions than clusters further from the centre. This result suggests a correlation between the survival of a cluster complex and its environment.

With the advent of modern observational efforts providing extensive giant molecular cloud catalogues, understanding the evolution of such clouds in a galactic context is of prime importance. While numerous previous numerical and theoretical works have focused on the cloud properties in isolated discs, few have looked into the cloud population in an interacting disc system. We present results of the first study investigating the evolution of the cloud population in galaxy experiencing an M51-like tidal fly-by using numerical simulations including star formation, interstellar medium cooling and stellar feedback. We see the cloud population shift to large unbound clouds in the wake of the companion passage, with the largest clouds appearing as fleeting short-lived agglomerations of smaller clouds within the tidal spiral arms, brought together by large scale streaming motions. These are then sheared apart as they leave the protection of the spiral arms. Clouds appear to lead diverse lives, even within similar environments, with some being born from gas shocked by filaments streaming into the spiral arms, and others from effectively isolated smaller colliding pairs. Overall this cloud population produces a shallower mass function than the disc in isolation, especially in the arms compared to the inter-arm regions. Direct comparisons to M51 observations show similarities between cloud populations, though models tailored to the mass and orbital models of M51 appear necessary to precisely reproduce the cloud population.

We present Phantom, a fast, parallel, modular and low-memory smoothed particle hydrodynamics and magnetohydrodynamics code developed over the last decade for astrophysical applications in three dimensions. The code has been developed with a focus on stellar, galactic, planetary and high energy astrophysics and has already been used widely for studies of accretion discs and turbulence, from the birth of planets to how black holes accrete. Here we describe and test the core algorithms as well as modules for magnetohydrodynamics, self-gravity, sink particles, H_2 chemistry, dust-gas mixtures, physical viscosity, external forces including numerous galactic potentials as well as implementations of Lense-Thirring precession, Poynting-Robertson drag and stochastic turbulent driving. Phantom is hereby made publicly available.

One of the main theories for explaining the formation of spiral arms in galaxies is the stationary density wave theory. This theory predicts the existence of an age gradient across the arms. We use the stellar cluster catalogues of the galaxies NGC 1566, M51a, and NGC 628 from the Legacy Extragalactic UV Survey (LEGUS) program. In order to test for the possible existence of an age sequence across the spiral arms, we quantified the azimuthal offset between star clusters of different ages in our target galaxies. We found that NGC 1566, a grand-design spiral galaxy with bisymmetric arms and a strong bar, shows a significant age gradient across the spiral arms that appears to be consistent with the prediction of the stationary density wave theory. In contrast, M51a with its two well-defined spiral arms and a weaker bar does not show an age gradient across the arms. In addition, a comparison with non LEGUS star cluster catalogues for M51a yields similar results. We believe that the spiral structure of M51a is not the result of a stationary density wave with a fixed pattern speed. Instead, tidal interactions could be the dominant mechanism for the formation of spiral arms. We also found no offset in the azimuthal distribution of star clusters with different ages across the weak spiral arms of NGC 628.

We perform simulations of isolated galaxies in order to investigate the likely origin of the spiral structure in M33. In our models, we find that gravitational instabilities in the stars and gas are able to reproduce the observed spiral pattern and velocity field of M33, as seen in HI, and no interaction is required. We also find that the optimum models have high levels of stellar feedback which create large holes similar to those observed in M33, whilst lower levels of feedback tend to produce a large amount of small scale structure, and undisturbed long filaments of high surface density gas, hardly detected in the M33 disc. The gas component appears to have a significant role in producing the structure, so if there is little feedback, both the gas and stars organise into clear spiral arms, likely due to a lower combined $Q$ (using gas and stars), and the ready ability of cold gas to undergo spiral shocks. By contrast models with higher feedback have weaker spiral structure, especially in the stellar component, compared to grand design galaxies. We did not see a large difference in the behaviour of $Q_{stars}$ with most of these models, however, because $Q_{stars}$ stayed relatively constant unless the disc was more strongly unstable. Our models suggest that although the stars produce some underlying spiral structure, this is relatively weak, and the gas physics has a considerable role in producing the large scale structure of the ISM in flocculent spirals.

We present a study of the dust-to-gas ratios in five nearby galaxies NGC 628 (M74), NGC 6503, NGC 7793, UGC 5139 (Holmberg I), and UGC 4305 (Holmberg II). Using Hubble Space Telescope broad band WFC3/UVIS UV and optical images from the Treasury program LEGUS (Legacy ExtraGalactic UV Survey) combined with archival HST/ACS data, we correct thousands of individual stars for extinction across these five galaxies using an isochrone-matching (reddening-free Q) method. We generate extinction maps for each galaxy from the individual stellar extinctions using both adaptive and fixed resolution techniques, and correlate these maps with neutral HI and CO gas maps from literature, including The HI Nearby Galaxy Survey (THINGS) and the HERA CO-Line Extragalactic Survey (HERACLES). We calculate dust-to-gas ratios and investigate variations in the dust-to-gas ratio with galaxy metallicity. We find a power law relationship between dust-to-gas ratio and metallicity, consistent with other studies of dust-to-gas ratio compared to metallicity. We find a change in the relation when H$_2$ is not included. This implies that underestimation of $N_{H_2}$ in low-metallicity dwarfs from a too-low CO-to-H$_2$ conversion factor $X_{CO}$ could have produced too low a slope in the derived relationship between dust-to-gas ratio and metallicity. We also compare our extinctions to those derived from fitting the spectral energy distribution (SED) using the Bayesian Extinction and Stellar Tool (BEAST) for NGC 7793 and find systematically lower extinctions from SED-fitting as compared to isochrone matching.

The interstellar abundances of refractory elements indicate a substantial depletion from the gas phase, that increases with gas density. Our recent model of dust evolution, based on hydrodynamic simulations of the lifecycle of giant molecular clouds (GMCs) proves that the observed trend for [Si$_{gas}$/H] is driven by a combination of dust growth by accretion in the cold diffuse interstellar medium (ISM) and efficient destruction by supernova (SN) shocks (Zhukovska et al. 2016). With an analytic model of dust evolution, we demonstrate that even with optimistic assumptions for the dust input from stars and without destruction of grains by SNe it is impossible to match the observed [Si$_{gas}$/H]$-n_H$ relation without growth in the ISM. We extend the framework developed in our previous work for silicates to include the evolution of iron grains and address a long-standing conundrum: ``Where is the interstellar iron?'. Much higher depletion of Fe in the warm neutral medium compared to Si is reproduced by the models, in which a large fraction of interstellar iron (70%) is locked as inclusions in silicate grains, where it is protected from sputtering by SN shocks. The slope of the observed [Fe$_{gas}$/H]$-n_H$ relation is reproduced if the remaining depleted iron resides in a population of metallic iron nanoparticles with sizes in the range of 1-10nm. Enhanced collision rates due to the Coulomb focusing are important for both silicate and iron dust models to match the observed slopes of the relations between depletion and density and the magnitudes of depletion at high density.

The Legacy ExtraGalactic UV Survey (LEGUS) is a multiwavelength Cycle 21 Treasury program on the Hubble Space Telescope. It studied 50 nearby star-forming galaxies in five bands from the near UV to the I-band, combining new Wide Field Camera 3 observations with archival Advanced Camera for Surveys data. LEGUS was designed to investigate how star formation occurs and develops on both small and large scales, and how it relates to the galactic environments. In this paper we present the photometric catalogs for all the apparently single stars identified in the 50 LEGUS galaxies. Photometric catalogs and mosaicked images for all filters are available for download. We present optical and near UV color-magnitude diagrams for all the galaxies. For each galaxy we derived the distance from the tip of the red giant branch. We then used the NUV color-magnitude diagrams to identify stars more massive than 14 Mo, and compared their number with the number of massive stars expected from the GALEX FUV luminosity. Our analysis shows that the fraction of massive stars forming in star clusters and stellar associations is about constant with the star formation rate. This lack of a relation suggests that the time scale for evaporation of unbound structures is comparable or longer than 10 Myr. At low star formation rates this translates to an excess of mass in clustered environments as compared to model predictions of cluster evolution, suggesting that a significant fraction of stars form in unbound systems.

Recently acquired WFC3 UV (F275W and F336W) imaging mosaics under the Legacy Extragalactic UV Survey (LEGUS) combined with archival ACS data of M51 are used to study the young star cluster (YSC) population of this interacting system. Our newly extracted source catalogue contains 2834 cluster candidates, morphologically classified to be compact and uniform in colour, for which ages, masses and extinction are derived. In this first work we study the main properties of the YSC population of the whole galaxy, considering a mass-limited sample. Both luminosity and mass functions follow a power law shape with slope -2, but at high luminosities and masses a dearth of sources is observed. The analysis of the mass function suggests that it is best fitted by a Schechter function with slope -2 and a truncation mass at $1.00\pm0.12\times10^5\ M_\odot$. Through Monte Carlo simulations we confirm this result and link the shape of the luminosity function to the presence of a truncation in the mass function. A mass limited age function analysis, between 10 and 200 Myr, suggests that the cluster population is undergoing only moderate disruption. We observe little variation in the shape of the mass function at masses above $1\times10^4\ M_\odot$, over this age range. The fraction of star formation happening in the form of bound clusters in M51 is $\sim20\%$ in the age range 10 to 100 Myr and little variation is observed over the whole range from 1 to 200 Myr.

In recent years there has been a growing interest in studying giant molecular filaments (GMFs), which are extremely elongated (> 100pc in length) giant molecular clouds (GMCs). They are often seen as inter-arm features in external spiral galaxies, but have been tentatively associated with spiral arms when viewed in the Milky Way. In this paper, we study the time evolution of GMFs in a high-resolution section of a spiral galaxy simulation, and their link with spiral arm GMCs and star formation, over a period of 11Myrs. The GMFs generally survive the inter-arm passage, although they are subject to a number of processes (e.g. star formation, stellar feedback and differential rotation) which can break the giant filamentary structure into smaller sections. The GMFs are not gravitationally bound clouds as a whole, but are, to some extent, confined by external pressure. Once they reach the spiral arms, the GMFs tend to evolve into more substructured spiral arm GMCs, suggesting that GMFs may be precursors to arm GMCs. Here, they become incorporated into the more complex and almost continuum molecular medium that makes up the gaseous spiral arm. Instead of retaining a clear filamentary shape, their shapes are distorted both by their climb up the spiral potential and their interaction with the gas within the spiral arm. The GMFs do tend to become aligned with the spiral arms just before they enter them (when they reach the minimum of the spiral potential), which could account for the observations of GMFs in the Milky Way.

We report the large effort which is producing comprehensive high-level young star cluster (YSC) catalogues for a significant fraction of galaxies observed with the Legacy ExtraGalactic UV Survey (LEGUS) Hubble treasury program. We present the methodology developed to extract cluster positions, verify their genuine nature, produce multiband photometry (from NUV to NIR), and derive their physical properties via spectral energy distribution fitting analyses. We use the nearby spiral galaxy NGC628 as a test case for demonstrating the impact that LEGUS will have on our understanding of the formation and evolution of YSCs and compact stellar associations within their host galaxy. Our analysis of the cluster luminosity function from the UV to the NIR finds a steepening at the bright end and at all wavelengths suggesting a dearth of luminous clusters. The cluster mass function of NGC628 is consistent with a power-law distribution of slopes $\sim -2$ and a truncation of a few times $10^5$ M$_\odot$. After their formation YSCs and compact associations follow different evolutionary paths. YSCs survive for a longer timeframe, confirming their being potentially bound systems. Associations disappear on time scales comparable to hierarchically organized star-forming regions, suggesting that they are expanding systems. We find mass-independent cluster disruption in the inner region of NGC628, while in the outer part of the galaxy there is little or no disruption. We observe faster disruption rates for low mass ($\leq$ $10^4$ M$_\odot$) clusters suggesting that a mass-dependent component is necessary to fully describe the YSC disruption process in NGC628.

The origin and life-cycle of molecular clouds are still poorly constrained, despite their importance for understanding the evolution of the interstellar medium. We have carried out a systematic, homogeneous, spectroscopic survey of the inner Galactic plane, in order to complement the many continuum Galactic surveys available with crucial distance and gas-kinematic information. Our aim is to combine this data set with recent infrared to sub-millimetre surveys at similar angular resolutions. The SEDIGISM survey covers 78 deg^2 of the inner Galaxy (-60 deg < l < +18 deg, |b| < 0.5 deg) in the J=2-1 rotational transition of 13CO. This isotopologue of CO is less abundant than 12CO by factors up to 100. Therefore, its emission has low to moderate optical depths, and higher critical density, making it an ideal tracer of the cold, dense interstellar medium. The data have been observed with the SHFI single-pixel instrument at APEX. The observational setup covers the 13CO(2-1) and C18O(2-1) lines, plus several transitions from other molecules. The observations have been completed. Data reduction is in progress, and the final data products will be made available in the near future. Here we give a detailed description of the survey and the dedicated data reduction pipeline. Preliminary results based on a science demonstration field covering -20 deg < l < -18.5 deg are presented. Analysis of the 13CO(2-1) data in this field reveals compact clumps, diffuse clouds, and filamentary structures at a range of heliocentric distances. By combining our data with data in the (1-0) transition of CO isotopologues from the ThrUMMS survey, we are able to compute a 3D realization of the excitation temperature and optical depth in the interstellar medium. Ultimately, this survey will provide a detailed, global view of the inner Galactic interstellar medium at an unprecedented angular resolution of ~30".

The process that leads to the formation of the bright star forming sites observed along prominent spiral arms remains elusive. We present results of a multi-wavelength study of a spiral arm segment in the nearby grand-design spiral galaxy M51 that belongs to a spiral density wave and exhibits nine gas spurs. The combined observations of the(ionized, atomic, molecular, dusty) interstellar medium (ISM) with star formation tracers (HII regions, young <10Myr stellar clusters) suggest (1) no variation in giant molecular cloud (GMC) properties between arm and gas spurs, (2) gas spurs and extinction feathers arising from the same structure with a close spatial relation between gas spurs and ongoing/recent star formation (despite higher gas surface densities in the spiral arm), (3) no trend in star formation age either along the arm or along a spur, (4) evidence for strong star formation feedback in gas spurs: (5) tentative evidence for star formation triggered by stellar feedback for one spur, and (6) GMC associations (GMAs) being no special entities but the result of blending of gas arm/spur cross-sections in lower resolution observations. We conclude that there is no evidence for a coherent star formation onset mechanism that can be solely associated to the presence of the spiral density wave. This suggests that other (more localized) mechanisms are important to delay star formation such that it occurs in spurs. The evidence of star formation proceeding over several million years within individual spurs implies that the mechanism that leads to star formation acts or is sustained over a longer time-scale.

Feedback from massive stars is thought to play an important role in the evolution of molecular clouds. In this work we analyse the effects of stellar winds and supernovae (SNe) in the evolution of two massive ($\sim 10^6\,M_\odot$) giant molecular clouds (GMCs): one gravitationally bound collapsing cloud and one unbound cloud undergoing disruption by galactic shear. These two clouds have been extracted from a large scale galaxy model and are re-simulated at a spatial resolution of $\sim 0.01$ pc, including feedback from winds, SNe, and the combined effect of both. We find that stellar winds stop accretion of gas onto sink particles, and can also trigger star formation in the shells formed by the winds, although the overall effect is to reduce the global star formation rate of both clouds. Furthermore, we observe that winds tend to escape through the corridors of diffuse gas. The effect of SNe is not so prominent and the star formation rate is similar to models neglecting stellar feedback. We find that most of the energy injected by the SNe is radiated away, but overdense areas are created by multiple and concurrent SN events especially in the most virialised cloud. Our results suggest that the impact of stellar feedback is sensitive to the morphology of star forming clouds, which is set by large scale galactic flows, being of greater importance in clouds undergoing gravitational collapse.

We investigate the properties and evolution of star particles in two simulations of isolated spiral galaxies, and two galaxies from cosmological simulations. Unlike previous numerical work, where typically each star particle represents one `cluster', for the isolated galaxies we are able to model features we term `clusters' with groups of particles. We compute the spatial distribution of stars with different ages, and cluster mass distributions, comparing our findings with observations including the recent LEGUS survey. We find that spiral structure tends to be present in older (100s Myrs) stars and clusters in the simulations compared to the observations. This likely reflects differences in the numbers of stars or clusters, the strength of spiral arms, and whether the clusters are allowed to evolve. Where we model clusters with multiple particles, we are able to study their evolution. The evolution of simulated clusters tends to follow that of their natal gas clouds. Massive, dense, long-lived clouds host massive clusters, whilst short-lived clouds host smaller clusters which readily disperse. Most clusters appear to disperse fairly quickly, in basic agreement with observational findings. We note that embedded clusters may be less inclined to disperse in simulations in a galactic environment with continuous accretion of gas onto the clouds than isolated clouds and correspondingly, massive young clusters which are no longer associated with gas tend not to occur in the simulations. Caveats of our models include that the cluster densities are lower than realistic clusters, and the simplistic implementation of stellar feedback.

We develop a model of dust evolution in a multiphase, inhomogeneous ISM including dust growth and destruction processes. The physical conditions for grain evolution are taken from hydrodynamical simulations of giant molecular clouds in a Milky Way-like spiral galaxy. We improve the treatment of dust growth by accretion in the ISM to investigate the role of the temperature-dependent sticking coefficient and ion-grain interactions. From detailed observational data on the gas-phase Si abundances [Si/H]_{gas} measured in the local Galaxy, we derive a relation between the average [Si/H]_{gas} and the local gas density n(H) which we use as a critical constraint for the models. This relation requires a sticking coefficient that decreases with the gas temperature. The synthetic relation constructed from the spatial dust distribution reproduces the slope of -0.5 of the observed relation in cold clouds. This slope is steeper than that for the warm medium and is explained by the dust growth. We find that it occurs for all adopted values of the minimum grain size a_{min} from 1 to 5nm. For the classical cut-off of a_{min}=5 nm, the ion-grain interactions result in longer growth timescales and higher [Si/H]_{gas} than the observed values. For a_{min} below 3 nm, the ion-grain interactions enhance the growth rates, steepen the slope of [Si/H]_{gas}-n(H) relation and provide a better match to observations. The rates of dust re-formation in the ISM by far exceed the rates of dust production by stellar sources as expected from simple evolution models. After the cycle of matter in and out of dust reaches a steady state, the dust growth balances the destruction operating on similar timescales of 350 Myr.

We examine how three fundamentally different numerical hydrodynamics codes follow the evolution of an isothermal galactic disc with an external spiral potential. We compare an adaptive mesh refinement code (RAMSES), a smoothed particle hydrodynamics code (sphNG), and a volume-discretised meshless code (GIZMO). Using standard refinement criteria, we find that RAMSES produces a disc that is less vertically concentrated and does not reach such high densities as the sphNG or GIZMO runs. The gas surface density in the spiral arms increases at a lower rate for the RAMSES simulations compared to the other codes. There is also a greater degree of substructure in the sphNG and GIZMO runs and secondary spiral arms are more pronounced. By resolving the Jeans' length with a greater number of grid cells we achieve more similar results to the Lagrangian codes used in this study. Other alterations to the refinement scheme (adding extra levels of refinement and refining based on local density gradients) are less successful in reducing the disparity between RAMSES and sphNG/GIZMO. Although more similar, sphNG displays different density distributions and vertical mass profiles to all modes of GIZMO (including the smoothed particle hydrodynamics version). This suggests differences also arise which are not intrinsic to the particular method but rather due to its implementation. The discrepancies between codes (in particular, the densities reached in the spiral arms) could potentially result in differences in the locations and timescales for gravitational collapse, and therefore impact star formation activity in more complex galaxy disc simulations.

We study the evolution of galactic magnetic fields using 3D smoothed particle magnetohydrodynamics (SPMHD) simulations of galaxies with an imposed spiral potential. We consider the appearance of reversals of the field, and amplification of the field. We find magnetic field reversals occur when the velocity jump across the spiral shock is above $\approx$20km s$^{-1}$, occurring where the velocity change is highest, typically at the inner Lindblad resonance (ILR) in our models. Reversals also occur at corotation, where the direction of the velocity field reverses in the co-rotating frame of a spiral arm. They occur earlier with a stronger amplitude spiral potential, and later or not at all with weaker or no spiral arms. The presence of a reversal at a radii of around 4--6 kpc in our fiducial model is consistent with a reversal identified in the Milky Way, though we caution that alternative Galaxy models could give a similar reversal. We find that relatively high resolution, a few million particles in SPMHD, is required to produce consistent behaviour of the magnetic field. Amplification of the magnetic field occurs in the models, and while some may be genuinely attributable to differential rotation or spiral arms, some may be a numerical artefact. We check our results using Athena, finding reversals but less amplification of the field, suggesting that some of the amplification of the field with SPMHD is numerical.

We study the properties of giant molecular clouds (GMCs) from a smoothed particle hydrodynamics simulation of a portion of a spiral galaxy, modelled at high resolution, with robust representations of the physics of the interstellar medium. We examine the global molecular gas content of clouds, and investigate the effect of using CO or H2 densities to define the GMCs. We find that CO can reliably trace the high-density H2 gas, but misses less dense H2 clouds. We also investigate the effect of using 3D CO densities versus CO emission with an observer's perspective, and find that CO-emission clouds trace well the peaks of the actual GMCs in 3D, but can miss the lower density molecular gas between density peaks which is often CO-dark. Thus, the CO emission typically traces smaller clouds within larger GMC complexes. We also investigate the effect of the galactic environment (in particular the presence of spiral arms), on the distribution of GMC properties, and we find that the mean properties are similar between arm and inter-arm clouds, but the tails of some distributions are indicative of intrinsic differences in the environment. We find highly filamentary clouds (similar to the giant molecular filaments of our Galaxy) exclusively in the inter-arm region, formed by galactic shear. We also find that the most massive GMC complexes are located in the arm, and that as a consequence of more frequent cloud interactions/mergers in the arm, arm clouds are more sub-structured and have higher velocity dispersions than inter-arm clouds.

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<1kpc of M51 to be 1 Msun/yr, using a new gravitational torque map and the molecular gas traced by the PdBI Arcsecond Whirlpool Survey (PAWS). The nuclear stellar bar is responsible for this gas inflow. We have also used torque profiles to estimate the location of dynamical resonances, suggesting a corotation for the bar at R~20", and a corotation for the spiral at R~100". We demonstrate how important it is to correct 3.6um images for dust emission in order to compute gravitational torques, and we carefully examine further sources of uncertainty. Our observational measurement of gas inflow can be compared with nuclear molecular outflow rates and provide useful constraints for numerical simulations.

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 (<100 Myr) stellar clusters within the galaxy, identified from a combination of visual inspection and automatic selection. We investigate the clustering of these young stellar clusters and quantify the strength and change of clustering strength with scale using the two-point correlation function. We also investigate how image boundary conditions and dust lanes affect the observed clustering. The distribution of the clusters is well fit by a broken power law with negative exponent $\alpha$. We recover a weighted mean index of $\alpha$ ~ -0.8 for all spatial scales below the break at 3".3 (158 pc at a distance of 9.9 Mpc) and an index of $\alpha$ ~ -0.18 above 158 pc for the accumulation of all cluster types. The strength of the clustering increases with decreasing age and clusters older than 40 Myr lose their clustered structure very rapidly and tend to be randomly distributed in this galaxy whereas the mass of the star cluster has little effect on the clustering strength. This is consistent with results from other studies that the morphological hierarchy in stellar clustering resembles the same hierarchy as the turbulent interstellar medium.

The nearby dwarf starburst galaxy NGC5253 hosts a number of young, massive star clusters, the two youngest of which are centrally concentrated and surrounded by thermal radio emission (the `radio nebula'). To investigate the role of these clusters in the starburst energetics, we combine new and archival Hubble Space Telescope images of NGC5253 with wavelength coverage from 1500 Ang to 1.9 micron in 13 filters. These include H-alpha, P-beta, and P-alpha, and the imaging from the Hubble Treasury Program LEGUS (Legacy Extragalactic UV Survey). The extraordinarily well-sampled spectral energy distributions enable modeling with unprecedented accuracy the ages, masses, and extinctions of the 9 optically brightest clusters (M_V < -8.8) and the two young radio nebula clusters. The clusters have ages ~1-15 Myr and masses ~1x10^4 - 2.5x10^5 M_sun. The clusters' spatial location and ages indicate that star formation has become more concentrated towards the radio nebula over the last ~15 Myr. The most massive cluster is in the radio nebula; with a mass 2.5x10^5 M_sun and an age ~1 Myr, it is 2-4 times less massive and younger than previously estimated. It is within a dust cloud with A_V~50 mag, and shows a clear nearIR excess, likely from hot dust. The second radio nebula cluster is also ~1 Myr old, confirming the extreme youth of the starburst region. These two clusters account for about half of the ionizing photon rate in the radio nebula, and will eventually supply about 2/3 of the mechanical energy in present-day shocks. Additional sources are required to supply the remaining ionizing radiation, and may include very massive stars.

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.

The arm structure of the Milky Way remains somewhat of an unknown, with observational studies hindered by our location within the Galactic disc. In the work presented here we use smoothed particle hydrodynamics (SPH) and radiative transfer to create synthetic longitude-velocity observations. Our aim is to reverse-engineer a top down map of the Galaxy by comparing synthetic longitude-velocity maps to those observed. We set up a system of N-body particles to represent the disc and bulge, allowing for dynamic creation of spiral features. Interstellar gas, and the molecular content, is evolved alongside the stellar system. A 3D-radiative transfer code is then used to compare the models to observational data. The resulting models display arm features that are a good reproduction of many of the observed emission structures of the Milky Way. These arms however are dynamic and transient, allowing for a wide range of morphologies not possible with standard density wave theory. The best fitting models are a much better match than previous work using fixed potentials. They favour a 4-armed model with a pitch angle of approximately 20 degrees, though with a pattern speed that decreases with increasing Galactic radius. Inner bars are lacking however, which appear required to fully reproduce the central molecular zone.

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.

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

We investigate cloud-cloud collisions, and GMC evolution, in hydrodynamic simulations of isolated galaxies. The simulations include heating and cooling of the ISM, self--gravity and stellar feedback. Over timescales $<5$ Myr most clouds undergo no change, and mergers and splits are found to be typically two body processes, but evolution over longer timescales is more complex and involves a greater fraction of intercloud material. We find that mergers, or collisions, occur every 8-10 Myr (1/15th of an orbit) in a simulation with spiral arms, and once every 28 Myr (1/5th of an orbit) with no imposed spiral arms. Both figures are higher than expected from analytic estimates, as clouds are not uniformly distributed in the galaxy. Thus clouds can be expected to undergo between zero and a few collisions over their lifetime. We present specific examples of cloud--cloud interactions in our results, including synthetic CO maps. We would expect cloud--cloud interactions to be observable, but find they appear to have little or no impact on the ISM. Due to a combination of the clouds' typical geometries, and moderate velocity dispersions, cloud--cloud interactions often better resemble a smaller cloud nudging a larger cloud. Our findings are consistent with the view that spiral arms make little difference to overall star formation rates in galaxies, and we see no evidence that collisions likely produce massive clusters. However, to confirm the outcome of such massive cloud collisions we ideally need higher resolution simulations.

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 (<~50 Myr) star formation histories from resolved massive stars and the extinction-corrected ages and masses of star clusters and associations. The extensive inventories of massive stars and clustered systems will be used to investigate the spatial and temporal evolution of star formation within galaxies. This will, in turn, inform theories of galaxy evolution and improve the understanding of the physical underpinning of the gas-star formation relation and the nature of star formation at high redshift. This paper describes the survey, its goals and observational strategy, and the initial science results. Because LEGUS will provide a reference survey and a foundation for future observations with JWST and with ALMA, a large number of data products are planned for delivery to the community.

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.

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.

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.

We present an investigation into the morphological features of the Milky Way. We use smoothed particle hydrodynamics (SPH) to simulate the interstellar medium (ISM) in the Milky Way under the effect of a number of different gravitational potentials representing spiral arms and bars, assuming the Milky Way is grand design in nature. The gas is subject to ISM cooling and chemistry, enabling us to track the evolution of molecular gas. We use a 3D radiative transfer code to simulate the emission from the SPH output, allowing for the construction of synthetic longitude-velocity (l-v) emission maps as viewed from the Earth. By comparing these maps with the observed emission in CO from the Milky Way, we infer the arm/bar geometry that provides a best fit to our Galaxy. We find that it is possible to reproduce nearly all features of the l-v diagram in CO emission. There is no model, however, that satisfactorily reproduces all of the features simultaneously. Models with 2 arms cannot reproduce all the observed arm features, while 4 armed models produce too bright local emission in the inner Galaxy. Our best fitting models favour a bar pattern speed within 50-60km/s/kpc and an arm pattern speed of approximately 20km/s/kpc, with a bar orientation of approximately 45 degrees and arm pitch angle between 10-15 degrees.

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

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.

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.

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)

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.

In this Letter, we analyse the distributions of stellar ages in Giant Molecular Clouds (GMCs) in spiral arms, inter-arm 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 (< 10 Myr) stars. By contrast, clouds which are the remainders of spiral arm GMAs that have been sheared into inter-arm GMCs, contain fewer young (< 10 Myr) stars, and more ~20 Myr stars. We also show that clouds which form in the absence of spiral arms, due to local gravitational and thermal instabilities, contain preferentially young stars. We propose the age distributions of stars in GMCs will be a useful diagnostic to test different cloud evolution scenarios, the origin of spiral arms, and the success of numerical models of galactic star formation. We discuss the implications of our results in the context of Galactic and extragalactic molecular clouds.

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.

Studying star formation in spiral arms tells us not only about the evolution of star formation, and molecular clouds, but can also tell us about the nature of spiral structure in galaxies. I will address both these topics using the results of recent simulations and observations. Galactic scale simulations are beginning to examine in detail the evolution of GMCs as they form in spiral arms, and then disperse by stellar feedback or shear. The overall timescale for this process appears comparable to the crossing time of the GMCs, a few Myrs for $10^5$ M$_{\odot}$ clouds, 20 Myr or so for more massive GMCs. Both simulations and observations show that the massive clouds are found in the spiral arms, likely as a result of cloud-cloud collisions. Simulations including stars should also tell us about the stellar age distribution in GMCs, and across spiral arms. More generally, recent work on spiral galaxies suggests that the dynamics of gas flows in spiral arms are different in longlived and transient spiral arms, resulting in different age patterns in the stars. Such results could be used to help establish the main driver of spiral structure in the Milky Way (Toomre instabilities, the bar, or nearby companion galaxies) in conjunction with future surveys.

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.

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 effect of gas streaming motions on the star-forming properties of individual GMCs. We compare our view of gas flows in M51 -- which arise due to departures from axi-symmetry in the gravitational potential (i.e. the nuclear bar and spiral arms) -- with the global pattern of star formation as traced by Halpha and 24\mu 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 prevents 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 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 to 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 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.

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)

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.

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.

We analyse the results of four simulations of isolated galaxies: two with a rigid spiral potential of fixed pattern speed, but with different degrees of star-formation induced feedback, one with an axisymmetric galactic potential and one with a `live' self-gravitating stellar component. Since we use a Lagrangian method we are able to select gas that lies within giant molecular clouds (GMCs) at a particular timeframe, and to then study the properties of this gas at earlier and later times. We find that gas which forms GMCs is not typical of the interstellar medium at least 50 Myr before the clouds form and reaches mean densities within an order of magnitude of mean cloud densities by around 10 Myr before. The gas in GMCs takes at least 50 Myr to return to typical ISM gas after dispersal by stellar feedback, and in some cases the gas is never fully recycled. We also present a study of the two-dimensional, vertically-averaged velocity fields within the ISM. We show that the velocity fields corresponding to the shortest timescales (that is, those timescales closest to the immediate formation and dissipation of the clouds) can be readily understood in terms of the various cloud formation and dissipation mechanisms. Properties of the flow patterns can be used to distinguish the processes which drive converging flows (e.g.\ spiral shocks, supernovae) and thus molecular cloud formation, and we note that such properties may be detectable with future observations of nearby galaxies.

We generate synthetic HI Galactic plane surveys from spiral galaxy simulations which include stellar feedback processes. Compared to a model without feedback we find an increased scale height of HI emission (in better agreement with observations) and more realistic spatial structure (including supernova blown bubbles). The synthetic data show HI self-absorption with a morphology similar to that seen in observations. The density and temperature of the material responsible for HI self-absorption is consistent with observationally determined values, and is found to be only weakly dependent on absorption strength and star formation efficiency.

We investigate the structure of the Milky Way by determining how features in a spatial map correspond to CO features in a velocity map. We examine structures including logarithmic spiral arms, a ring and a bar. We explore the available parameter space, including the pitch angle of the spiral arms, radius of a ring, and rotation curve. We show that surprisingly, a spiral arm provides a better fit to the observed molecular ring than a true ring feature. This is because both a spiral arm, and the observed feature known as the molecular ring, are curved in velocity longitude space. We find that much of the CO emission in the velocity longitude map can be fitted by a nearly symmetric 2 armed spiral pattern. One of the arms corresponds to the molecular ring, whilst the opposite arm naturally reproduces the Perseus arm. Multiple arms also contribute to further emission in the vicinity of the molecular ring and match other observed spiral arms. Whether the Galactic structure consists primarily of two, or several spiral arms, the presence of 2 symmetric logarithmic spirals, which begin in the vicinity of the ends of the bar, suggest a spiral density wave associated with the bar.

We perform calculations of isolated disc galaxies to investigate how the properties of the ISM, the nature of molecular clouds, and the global star formation rate depend on the level of stellar feedback. We adopt a simple physical model, which includes a galactic potential, a standard cooling and heating prescription of the ISM, and self gravity of the gas. Stellar feedback is implemented by injecting energy into dense, gravitationally collapsing gas, but is independent of the Schmidt-Kennicutt relation. We obtain fractions of gas, and filling factors for different phases of the ISM in reasonable ageement with observations. Supernovae are found to be vital to reproduce the scale heights of the different components of the ISM, and velocity dispersions. The GMCs formed in the simulations display mass spectra similar to the observations, their normalisation dependent on the level of feedback. We find ~40 per cent of the clouds exhibit retrograde rotation, induced by cloud-cloud collisions. The star formation rates we obtain are in good agreement with the observed Schmidt-Kennicutt relation, and are not strongly dependent on the star formation efficiency we assume, being largely self regulated by the feedback. We also investigate the effect of spiral structure by comparing calculations with and without the spiral component of the potential. The main difference with a spiral potential is that more massive GMCs are able to accumulate in the spiral arms. Thus we are able to reproduce massive GMCs, and the spurs seen in many grand design galaxies, even with stellar feedback. The presence of the spiral potential does not have an explicit effect on the star formation rate, but can increase the star formation rate indirectly by enabling the formation of long-lived, strongly bound clouds.

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 $\mu$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.

We have discovered long-lived waves in two sets of numerical models of fast (marginally bound or unbound) flyby galaxy collisions, carried out independently with two different codes. In neither simulation set are the spirals the result of a collision-induced bar formation. Although there is variation in the appearance of the waves with time, they do not disappear and reform recurrently, as seen in other cases described in the literature. We also present an analytic theory that can account for the wave structure, not as propagating transients, nor as a fixed pattern propagating through the disc. While these waves propagate through the disc, they are maintained by the coherent oscillations initiated by the impulsive disturbance. Specifically, the analytic theory suggests that they are caustic waves in ensembles of stars pursuing correlated epicyclic orbits after the disturbance. This theory is an extension of that developed by Struck and collaborators for colliding ring galaxies. The models suggest that this type of wave may persist for a couple of Gyr., and galaxy interactions occur on comparable timescales, so waves produced by the mechanism may be well represented in observed spirals. In particular, this mechanism can account for the tightly wound, and presumably long-lived spirals, seen in some nearby early-type galaxies. These spirals are also likely to be common in groups and clusters, where fast encounters between galaxies occur relatively frequently. However, as the spirals become tightly wound, and evolve to modest amplitudes, they may be difficult to resolve unless they are nearby. Nonetheless, the effect may be one of several processes that result from galaxy harassment, and via wave-enhanced star formation contribute to the Butcher-Oemler effect.

The most recent observational evidence seems to indicate that giant molecular clouds are predominantly gravitationally unbound objects. In this paper we show that this is a natural consequence of a scenario in which cloud-cloud collisions and stellar feedback regulate the internal velocity dispersion of the gas, and so prevent global gravitational forces from becoming dominant. Thus, while the molecular gas is for the most part gravitationally unbound, local regions within the denser parts of the gas (within the clouds) do become bound and are able to form stars. We find that the observations, in terms of distributions of virial parameters and cloud structures, can be well modelled provided that the star formation efficiency in these bound regions is of order 5 - 10 percent. We also find that in this picture the constituent gas of individual molecular clouds changes over relatively short time scales, typically a few Myr.

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.

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.

We consider models of gas flow in spiral galaxies in which the spiral structure has been excited by various possible mechanisms: a global steady density wave, self-gravity of the stellar disc and an external tidal interaction, as well as the case of a galaxy with a central rotating bar. In each model we estimate in a simple manner the likely current positions of star clusters of a variety of ages, ranging from ~ 2 Myr to around 130 Myr, depending on the model. We find that the spatial distribution of cluster of different ages varies markedly depending on the model, and propose that observations of the locations of age-dated stellar clusters is a possible discriminant between excitation mechanisms for spiral structure in an individual galaxy.

We have created synthetic neutral hydrogen (HI) Galactic Plane Survey data cubes covering 90 degrees < l < 180 degrees, using a model spiral galaxy from SPH simulations and the radiative transfer code TORUS. The density, temperature and other physical parameters are fed from the SPH simulation into TORUS, where the HI emissivity and opacity are calculated before the 21-cm line emission profile is determined. Our main focus is the observation of Outer Galaxy `Perseus Arm' HI, with a view to tracing atomic gas as it encounters shock motions as it enters a spiral arm interface, an early step in the formation of molecular clouds. The observation of HI self-absorption features at these shock sites (in both real observations and our synthetic data) allows us to investigate further the connection between cold atomic gas and the onset of molecular cloud formation.

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.

We investigate pattern speeds in spiral galaxies where the structure is induced by an interaction with a companion galaxy. We perform calculations modeling the response of the stellar and/or gaseous components of a disc. Generally we do not find a unique pattern speed in these simulations, rather the pattern speed decreases with radius, and the pattern speed for individual spiral arms differ. The maximum pattern speed is ~20 km/s/kpc for the discs with a live stellar component, decreasing to 5 km/s/kpc at the edge of the spiral perturbation. When only the gas is modeled, Omega_p is typically very low (5 km/s/kpc) at all radii.

We present hydrodynamical models of the grand design spiral M51 (NGC 5194), and its interaction with its companion NGC 5195. Despite the simplicity of our models, our simulations capture the present day spiral structure of M51 remarkably well, and even reproduce details such as a kink along one spiral arm, and spiral arm bifurcations. We investigate the offset between the stellar and gaseous spiral arms, and find at most times (including the present day) there is no offset between the stars and gas to within our error bars. We also compare our simulations with recent observational analysis of M51. We compute the pattern speed versus radius, and like the observations, find no single global pattern speed. We also show that the spiral arms cannot be fitted well by logarithmic spirals. We interpret these findings as evidence that M51 does not exhibit a quasi-steady density wave, as would be predicted by density wave theory. The internal structure of M51 derives from the complicated and dynamical interaction with its companion, resulting in spiral arms showing considerable structure in the form of short-lived kinks and bifurcations. Rather than trying to model such galaxies in terms of global spiral modes with fixed pattern speeds, it is more realistic to start from a picture in which the spiral arms, while not being simple material arms, are the result of tidally induced kinematic density `waves' or density patterns, which wind up slowly over time.

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.

The M51 galaxy is one of best examples of grand design structure in galaxies. The most likely origin of the spiral arms in M51 is the ongoing interaction with NGC 5195. Here we report recent calculations which model the orbit of M51 and NGC 5195. These calculations, which for the first time focus on the gas dynamics, are able to reproduce even detailed features in M51, as well as demonstrate the nature of the spiral structure of M51.

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.

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.

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 (<~1 Myr). From this stage, the majority of gas generally takes between 10 to 20 Myr to obtain high H$_{2}$ fractions (>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.

The tight correlation between galaxy bulges and their central black hole masses likely emerges in a phase of rapid collapse and starburst at high redshift, due to the balance of gravity on gas with the feedback force from starbursts and the wind from the black hole; the average gravity on per unit mass of gas is ~ 2 x 10^-10 m/sec^2 during the star burst phase. This level of gravity could come from the real r^{-1} cusps of Cold Dark Matter (CDM) halos, but the predicted gravity would have a large scatter due to dependence on cosmological parameters and formation histories. Better agreement is found with the gravity from the scalar field in some co-variant versions of MOND, which can create the mirage of a Newtonian effective dark halo of density Pi r^{-1} near the center, where the characteristic surface density Pi=130alpha^{-1} Msun pc^{-2} and alpha is a fundamental constant of order unity fixed by the Lagrangian of the co-variant theory if neglecting environmental effects. We show with a toy analytical model and a hydrodynamical simulation that a constant background gravity due to MOND/TeVeS scalar field implies a critical pressure synchronizing starbursts and the formation of galaxy bulges and ellipticals. A universal threshold for the formation of the brightest regions of galaxies in a MONDian universe suggests that the central BHs, bulges and ellipticals would respect tight correlations like the M_{bulge}-M_{BH}-sigma relations. In general MOND tends to produce tight correlations in galaxy properties because its effective halo has less freedom and scatter than CDM halos.

Magnetic fields are important at every scale in the star formation process: from the dynamics of the ISM in galaxies, to the collapse of turbulent molecular clouds to form stars and in the fragmentation of individual star forming cores. The recent development of a robust algorithm for MHD in the Smoothed Particle Hydrodynamics method has enabled us to perform simulations of star formation including magnetic fields at each of these scales. This paper focusses on three questions in particular: What is the effect of magnetic fields on fragmentation in star forming cores? How do magnetic fields affect the collapse of turbulent molecular clouds to form stars? and: What effect do magnetic fields have on the dynamics of the interstellar medium?

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.

We investigate the dynamics of magnetic fields in spiral galaxies by performing 3D MHD simulations of galactic discs subject to a spiral potential. Recent hydrodynamic simulations have demonstrated the formation of inter-arm spurs as well as spiral arm molecular clouds provided the ISM model includes a cold HI phase. We find that the main effect of adding a magnetic field to these calculations is to inhibit the formation of structure in the disc. However, provided a cold phase is included, spurs and spiral arm clumps are still present if $\beta \gtrsim 0.1$ in the cold gas. A caveat to two phase calculations though is that by assuming a uniform initial distribution, $\beta \gtrsim 10$ in the warm gas, emphasizing that models with more consistent initial conditions and thermodynamics are required. Our simulations with only warm gas do not show such structure, irrespective of the magnetic field strength. Furthermore, we find that the introduction of a cold HI phase naturally produces the observed degree of disorder in the magnetic field, which is again absent from simulations using only warm gas. Whilst the global magnetic field follows the large scale gas flow, the magnetic field also contains a substantial random component that is produced by the velocity dispersion induced in the cold gas during the passage through a spiral shock. Without any cold gas, the magnetic field in the warm phase remains relatively well ordered apart from becoming compressed in the spiral shocks. Our results provide a natural explanation for the observed high proportions of disordered magnetic field in spiral galaxies and we thus predict that the relative strengths of the random and ordered components of the magnetic field observed in spiral galaxies will depend on the dynamics of spiral shocks.

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.

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.

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.

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.

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.

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.

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.