We investigate star formation from subpc to kpc scales with magnetohydrodynamic models of a cloud structure and a section of galactic spiral arm. We aim to understand how magnetic fields affect star formation and cloud formation, and how feedback couples with magnetic fields on scales of clouds and clumps. We find that magnetic fields overall suppress star formation by similar to 10 per cent with a weak field (5 mu G) and similar to 50 per cent with a stronger field (50 mu G). Cluster masses are reduced by about 40 per cent with a strong field but show little change with a weak field. We find that clouds tend to be aligned parallel to the field with a weak field and become perpendicularly aligned with a stronger field, whereas on clump scales the alignment is more random. The magnetic fields and densities of clouds and clumps in our models agree with the Zeeman measurements of the Crutcher relation B-rho in the weaker field models, while the strongest field models show a relation that is too flat compared to the observations. In all our models, we find that both subcritical and supercritical clouds and clumps are present. We also find that if using a line-of-sight (1D) measure of the magnetic field to determine the critical parameter, the magnetic field, and thereby also criticality, can vary by a factor of 3-4 depending on whether the direction the field is measured along corresponds to the direction of the ordered component of the magnetic field.

We compare the molecular cloud properties in subgalactic regions of two galaxies, barred spiral NGC 1313, which is forming many massive clusters, and flocculent spiral NGC 7793, which is forming significantly fewer massive clusters despite having a similar star formation rate to NGC 1313. We find that there are larger variations in cloud properties between different regions within each galaxy than there are between the galaxies on a global scale, especially for NGC 1313. There are higher masses, line widths, pressures, and virial parameters in the arms of NGC 1313 and the center of NGC 7793 than in the interarm and outer regions of the galaxies. The massive cluster formation of NGC 1313 may be driven by its greater variation in environment, allowing more clouds with the necessary conditions to emerge, although no one parameter seems primarily responsible for the difference in star formation. Meanwhile NGC 7793 has clouds that are as massive and have as much kinetic energy as the clouds in the arms of NGC 1313, but have densities and pressures more similar to those in the interarm regions and so are less inclined to collapse and form stars. The cloud properties in NGC 1313 and NGC 7793 suggest that spiral arms, bars, interarm regions, and flocculent spirals each represent distinct environments with regard to molecular cloud populations. We see surprisingly little difference in surface density between the regions, suggesting that the differences in surface densities frequently seen between arm and interarm regions in lower-resolution studies are indicative of the sparsity of molecular clouds, rather than differences in their true surface density.

The Whirlpool Galaxy is a well-studied grand design galaxy with two major spiral arms, and a large satellite NGC 5195. The arms both show long uniform sections with perturbations (“kinks” or sharp turns) in specific regions. Comparing the two arms shows a small radial offset between the main kinked regions. We analyzed the morphology and also the velocity field in the disk of M51 using kinematic maps based on H α and CO line emission. These sample complementary radial ranges, with the CO map covering the central zone and the H α map extending to cover the outer zone. We looked for indicators of density wave resonance, zones where radial flows of gas in the disk plane reverse their sign. These were present in both velocity maps; their 2D localization placed them along or closely parallel to the spiral arms, at a set of well-defined galactocentric radii, and notably more concentrated along the southern, stronger arm. The results can be well interpreted quantitatively, using a numerical model of the interaction of M51 and NGC 5195 in which the satellite has made two relatively recent passes through the disk plane of M51. During the first pass the pair of dominant spiral arms was stimulated, and during the second pass the strong kinks in both arms were formed at about the same time. The second interaction is particularly well characterized, because the timescale corresponding to the production of the kinks and the recovery of the original pitch angle is identical for the two arms.

We present a comparative study of the molecular gas in two galaxies from the LEGUS sample: barred spiral NGC 1313 and flocculent spiral NGC 7793. These two galaxies have similar masses, metallicities, and star formation rates, but NGC 1313 is forming significantly more massive star clusters than NGC 7793, especially young massive clusters (<10 Myr, >10^4 Msol). Using ALMA CO(2-1) observations of the two galaxies with the same sensitivities and resolutions of 13 pc, we directly compare the molecular gas in these two similar galaxies to determine the physical conditions responsible for their large disparity in cluster formation. By fitting size-linewidth relations for the clouds in each galaxy, we find that NGC 1313 has a higher intercept than NGC 7793, implying that its clouds have higher kinetic energies at a given size scale. NGC 1313 also has more clouds near virial equilibrium than NGC 7793, which may be connected to its higher rate of massive cluster formation. However, these virially bound clouds do not show a stronger correlation with young clusters than that of the general cloud population. We find surprisingly small differences between the distributions of molecular cloud populations in the two galaxies, though the largest of those differences are that NGC 1313 has higher surface densities and lower free-fall times.

Simulations from the scales of isolated galaxies to clouds have been instrumental in informing us about molecular cloud formation and evolution. Simulations are able to investigate the roles of gravity, feedback, turbulence, heating and cooling, and magnetic fields on the physics of the interstellar medium, and star formation. Compared to simulations of individual clouds, galactic and sub-galactic scale simulations can include larger galactic scale processes such as spiral arms, bars, and larger supernovae bubbles, which may influence star formation. Simulations show cloud properties and lifetimes in broad agreement with observations. Gravity and spiral arms are required to produce more massive GMCs, whilst stellar feedback, likely photoionisation, leads to relatively short cloud lifetimes. On larger scales, supernovae may be more dominant in driving the structure and dynamics, but photoionisation may still have a role. In terms of the dynamics, feedback is probably the main driver of velocity dispersions, but large scale processes such as gravity and spiral arms may also be significant. Magnetic fields are generally found to decrease star formation on galaxy or cloud scales, and simulations are ongoing to study whether clouds are sub or supercritical on different scales in galaxy scale simulations. Simulations on subgalactic scales, or zoom in simulations, allow better resolution of feedback processes, filamentary structure within clouds, and the study of stellar clusters.

ABSTRACTWe present hydrodynamic and magnetohydrodynamic (MHD) simulations of sub-galactic regions including photoionizing and supernova feedack. We aim to improve the initial conditions of our region extraction models by including an initial population of stars. We also investigate the reliability of extracting regions in simulations, and show that with a good choice of region, results are comparable with using a larger region for the duration of our simulations. Simulations of star formation on molecular cloud scales typically start with a turbulent cloud of gas, from which stars form and then undergo feedback. In reality, a typical cloud or region within a galaxy may already include, or reside near some population of stars containing massive stars undergoing feedback. We find the main role of a prior population is triggering star formation, and contributing to gas dynamics. Early time supernova from the initial population are important in triggering new star formation and driving gas motions on larger scales above 100 pc, whilst the ionizing feedback contribution from the initial population has less impact, since many members of the initial population have cleared out gas around them in the prior model. In terms of overall star formation rates though, the initial population has a relatively small effect, and the feedback does not for example suppress subsequent star formation. We find that MHD has a relatively larger impact than initial conditions, reducing the star formation rate by a factor of 3 at later times.

Simulations from the scales of isolated galaxies to clouds have been instrumental in informing us about molecular cloud formation and evolution. Simulations are able to investigate the roles of gravity, feedback, turbulence, heating and cooling, and magnetic fields on the physics of the interstellar medium, and star formation. Compared to simulations of individual clouds, galactic and sub-galactic scale simulations can include larger galactic scale processes such as spiral arms, bars, and larger supernovae bubbles, which may influence star formation. Simulations show cloud properties and lifetimes in broad agreement with observations. Gravity and spiral arms are required to produce more massive GMCs, whilst stellar feedback, likely photoionisation, leads to relatively short cloud lifetimes. On larger scales, supernovae may be more dominant in driving the structure and dynamics, but photoionisation may still have a role. In terms of the dynamics, feedback is probably the main driver of velocity dispersions, but large scale processes such as gravity and spiral arms may also be significant. Magnetic fields are generally found to decrease star formation on galaxy or cloud scales, and simulations are ongoing to study whether clouds are sub or supercritical on different scales in galaxy scale simulations. Simulations on subgalactic scales, or zoom in simulations, allow better resolution of feedback processes, filamentary structure within clouds, and the study of stellar clusters.

Observations of clusters suffer from issues such as completeness, projection effects, resolving individual stars, and extinction. As such, how accurate measurements and conclusions are likely to be? Here, we take cluster simulations (Westerlund2- and Orion-type), synthetically observe them to obtain luminosities, accounting for extinction, and the inherent limits of Gaia, then place them within the real Gaia DR3 catalogue. We then attempt to rediscover the clusters at distances of between 500 and 4300 pc. We show the spatial and kinematic criteria that are best able to pick out the simulated clusters, maximizing completeness, and minimizing contamination. We then compare the properties of the 'observed' clusters with the original simulations. We looked at the degree of clustering, the identification of clusters and subclusters within the data sets, and whether the clusters are expanding or contracting. Even with a high level of incompleteness (e.g. < 2 per cent stellar members identified), similar qualitative conclusions tend to be reached compared to the original data set, but most quantitative conclusions are likely to be inaccurate. Accurate determination of the number, stellar membership, and kinematic properties of subclusters are the most problematic to correctly determine, particularly at larger distances due to the disappearance of cluster substructure as the data become more incomplete, but also at smaller distances where the misidentification of asterisms as true structure can be problematic. Unsurprisingly, we tend to obtain better quantitative agreement of properties for our more massive Westerlund2-type cluster. We also make optical style images of the clusters over our range of distances.

It remains unclear how galactic environment affects star formation and stellar cluster properties. This is difficult to address in Milky Way-mass galaxy simulations because of limited resolution and less accurate feedback compared to cloud-scale models. We carry out zoom-in simulations to re-simulate 100-300 pc regions of a Milky Way-like galaxy using smoothed particle hydrodynamics, including finer resolution (0.4 M-circle dot per particle), cluster-sink particles, ray-traced photoionization from O stars, H-2/CO chemistry, and interstellar medium heating/cooling. We select similar to 10(6) M-circle dot cloud complexes from a galactic bar, inner spiral arm, outer arm, and inter-arm region (in order of galactocentric radius), retaining the original galactic potentials. The surface densities of star formation rate and neutral gas follow Sigma(SFR) proportional to Sigma(1.3)(gas), with the bar lying higher up the relation than the other regions. However, the inter-arm region forms stars two to three times less efficiently than the arm models at the same Sigma(gas). The bar produces the most massive cluster, the inner arm the second, and the inter-arm the third. Almost all clusters in the bar and inner arm are small (radii <5 pc), while 30-50 per cent of clusters in the outer arm and inter-arm have larger radii more like associations. Bar and inner arm clusters rotate at least twice as fast, on average, than clusters in the outer arm and inter-arm regions. The degree of spatial clustering also decreases from bar to inter-arm. Our results indicate that young massive clusters, potentially progenitors of globular clusters, may preferentially form near the bar/inner arm compared to outer arm/inter-arm regions.

We investigate the variation in the upper end of the stellar initial mass function (uIMF) in 375 young and compact star clusters in five nearby galaxies within ∼5 Mpc. All the young stellar clusters (YSCs) in the sample have ages ≲ 4 Myr and masses above 500 M ⊙ , according to standard stellar models. The YSC catalogs were produced from Hubble Space Telescope images obtained as part of the Legacy ExtraGalactic UV Survey (LEGUS) Hubble treasury program. They are used here to test whether the uIMF is universal or changes as a function of the cluster’s stellar mass. We perform this test by measuring the H α luminosity of the star clusters as a proxy for their ionizing photon rate, and charting its trend as a function of cluster mass. Large cluster numbers allow us to mitigate the stochastic sampling of the uIMF. The advantage of our approach relative to previous similar attempts is the use of cluster catalogs that have been selected independently of the presence of H α emission, thus removing a potential sample bias. We find that the uIMF, as traced by the H α emission, shows no dependence on cluster mass, suggesting that the maximum stellar mass that can be produced in star clusters is universal, in agreement with previous findings.

ABSTRACT We present simulations of a 500 pc2 region, containing gas of mass 4 × 106 M⊙, extracted from an entire spiral galaxy simulation, scaled up in resolution, including photoionizing feedback from stars of mass >18 M⊙. Our region is evolved for 10 Myr and shows clustered star formation along the arm generating ≈ 5000 cluster sink particles ≈ 5 per cent of which contain at least one of the ≈ 4000 stars of mass >18 M⊙. Photoionization has a noticeable effect on the gas in the region, producing ionized cavities and leading to dense features at the edge of the H ii regions. Compared to the no-feedback case, photoionization 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 ionization, with some already dense gas, becoming compressed from multiple sides by neighbouring H ii regions. This causes star formation that is both accelerated and partially displaced throughout the spiral arm with up to 30 per cent 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.

We present Ekster, a new method for simulating star clusters from birth in a live galaxy simulation that combines the smoothed-particle hydrodynamics (SPH) method Phantom with the $N$-body method PeTar. With Ekster, it becomes possible to simulate individual stars in a simulation with only moderately high resolution for the gas, allowing us to study whole sections of a galaxy rather than be restricted to individual clouds. We use this method to simulate star and star cluster formation in spiral arms, investigating massive GMCs and spiral arm regions with lower mass clouds, from two galaxy models with different spiral potentials. After selecting these regions from pre-run galaxy simulations, we re-sample the particles to obtain a higher resolution. We then re-simulate these regions for 3 Myr to study where and how star clusters form. We analyse the early evolution of the embedded star clusters in these regions. We find that the massive GMC regions, which are more common with stronger spiral arms, form more massive clusters than the sections of spiral arms containing lower mass clouds. Clusters form both by accreting gas and by merging with other proto-clusters, the latter happening more frequently in the denser GMC regions.

ABSTRACT Understanding the formation and evolution of high-mass star clusters requires comparisons between theoretical and observational data to be made. Unfortunately, while the full phase space of simulated regions is available, often only partial 2D spatial and kinematic data is available for observed regions. This raises the question as to whether cluster parameters determined from 2D data alone are reliable and representative of clusters real parameters and the impact of line-of-sight orientation. In this paper, we derive parameters for a simulated cluster formed from a cloud–cloud collision with the full 6D phase space, and compare them with those derived from three different 2D line-of-sight orientations for the cluster. We show the same qualitative conclusions can be reached when viewing clusters in 2D versus 3D, but that drawing quantitative conclusions when viewing in 2D is likely to be inaccurate. The greatest divergence occurs in the perceived kinematics of the cluster, which in some orientations appears to be expanding when the cluster is actually contracting. Increases in the cluster density compounds pre-existing perspective issues, reducing the relative accuracy and consistency of properties derived from different orientations. This is particularly problematic for determination of the number, and membership, of subclusters present in the cluster. We find the fraction of subclusters correctly identified in 2D decreases as the cluster evolves, reaching less than $3.4{{\ \rm percent}}$ at the evolutionary end point for our cluster.

The role of different stellar feedback mechanisms in giant molecular clouds is not well understood. This is especially true for regions with many interacting clouds as would be found in a galactic spiral arm. In this paper, building on previous work by Bending et al., we extract a 500 pc x 500 pc x 100 pc section of a spiral arm from a galaxy simulation. We use smoothed particle hydrodynamics to re-simulate the region at higher resolution (1 M-circle dot per particle). We present a method for momentum-driven stellar winds from main-sequence massive stars, and include this with photoionization, self-gravity, a galactic potential, and interstellar medium heating/cooling. We also include cluster-sink particles with accretion radii of 0.78 pc to track star/cluster formation. The feedback methods are as robust as previous models on individual cloud scales (e.g. Dale et al.). We find that photoionization dominates the disruption of the spiral arm section, with stellar winds only producing small cavities (at most similar to 30 pc). Stellar winds do not affect the resulting cloud statistics or the integrated star formation rate/efficiency, unlike ionization, which produces more stars, and more clouds of higher density and higher velocity dispersion compared to the control run without feedback. Winds do affect the sink properties, distributing star formation over more low-mass sinks (similar to 10(2)M(circle dot)) and producing fewer high-mass sinks (similar to 10(3) M-circle dot). Overall, stellar winds play at best a secondary role compared to photoionization, and on many measures, they have a negligible impact.

Modelling star formation and resolving individual stars in numerical simulations of molecular clouds and galaxies is highly challenging. Simulations on very small scales can be sufficiently well resolved to consistently follow the formation of individual stars, whilst on larger scales sinks that have masses sufficient to fully sample the IMF can be converted into realistic stellar populations. However, as yet, these methods do not work for intermediate scale resolutions whereby sinks are more massive compared to individual stars but do not fully sample the IMF. In this paper, we introduce the grouped star formation prescription, whereby sinks are first grouped according to their positions, velocities, and ages, then stars are formed by sampling the IMF using the mass of the groups. We test our grouped star formation method in simulations of various physical scales, from sub-parsec to kilo-parsec, and from static isolated clouds to colliding clouds. With suitable grouping parameters, this star formation prescription can form stars that follow the IMF and approximately mimic the original stellar distribution and velocity dispersion. Each group has properties that are consistent with a star-forming region. We show that our grouped star formation prescription is robust and can be adapted in simulations with varying physical scales and resolution. Such methods are likely to become more important as galactic or even cosmological scale simulations begin to probe sub-parsec scales.

ABSTRACT The SEDIGISM (Structure, Excitation and Dynamics of the Inner Galactic Interstellar Medium) survey used the APEX telescope to map 84 deg2 of the Galactic plane between ℓ = −60° and +31° in several molecular transitions, including 13CO (2 – 1) and C18O (2 – 1), thus probing the moderately dense (∼103 cm−3) component of the interstellar medium. With an angular resolution of 30 arcsec and a typical 1σ sensitivity of 0.8–1.0 K at 0.25 km s−1 velocity resolution, it gives access to a wide range of structures, from individual star-forming clumps to giant molecular clouds and complexes. The coverage includes a good fraction of the first and fourth Galactic quadrants, allowing us to constrain the large-scale distribution of cold molecular gas in the inner Galaxy. In this paper, we provide an updated overview of the full survey and the data reduction procedures used. We also assess the quality of these data and describe the data products that are being made publicly available as part of this First Data Release (DR1). We present integrated maps and position–velocity maps of the molecular gas and use these to investigate the correlation between the molecular gas and the large-scale structural features of the Milky Way such as the spiral arms, Galactic bar and Galactic Centre. We find that approximately 60 per cent of the molecular gas is associated with the spiral arms and these appear as strong intensity peaks in the derived Galactocentric distribution. We also find strong peaks in intensity at specific longitudes that correspond to the Galactic Centre and well-known star-forming complexes, revealing that the 13CO emission is concentrated in a small number of complexes rather than evenly distributed along spiral arms.

AbstractWe present Ekster, a new method for simulating the formation and dynamics of individual stars in a relatively low-resolution gas background. Here, we use Ekster to simulate star cluster formation in two different regions from each of two galaxy models with different spiral potentials. We simulate these regions for 3 Myr to study where and how star clusters form. We find that massive GMC regions form more massive clusters than sections of spiral arms. Additionally we find that clusters form both by accreting gas and by merging with other proto-clusters, the latter happening more frequently in the denser GMC regions.

ABSTRACT Young massive clusters (YMCs) are recently formed astronomical objects with unusually high star formation rates. We propose the collision of giant molecular clouds (GMCs) as a likely formation mechanism of YMCs, consistent with the YMC conveyor-belt formation mode concluded by other authors. We conducted smoothed particle hydrodynamical simulations of cloud–cloud collisions and explored the effect of the clouds’ collision speed, initial cloud density, and the level of cloud turbulence on the global star formation rate and the properties of the clusters formed from the collision. We show that greater collision speed, greater initial cloud density and lower turbulence increase the overall star formation rate and produce clusters with greater cluster mass. In general, collisions with relative velocity ≳ 25 km s−1, initial cloud density ≳ 250 cm−3, and turbulence of ≈2.5 km s−1 can produce massive clusters with properties resembling the observed Milky Way YMCs.

The nature of galactic spiral arms in disc galaxies remains elusive. Regardless of the spiral model, arms are expected to play a role in sculpting the star-forming interstellar medium. As such, different arm models may result in differences in the structure of the interstellar medium and molecular cloud properties. In this study we present simulations of galactic discs subject to spiral arm perturbations of different natures. We find very little difference in how the cloud population or gas kinematics vary between the different grand-design spirals, indicting that the interstellar medium on cloud scales cares little about where spiral arms come from. We do, however, see a difference in the interarm/arm mass spectra, {and minor differences in tails of the distributions of cloud properties} (as well as radial variations in the stellar/gaseous velocity dispersions). These features can be attributed to differences in the radial dependence of the pattern speeds between the different spiral models, and could act as a metric of the nature of spiral structure in observational studies.

AbstractWe present a model for hydrodynamic + N-body simulations of star cluster formation and evolution using AMUSE. Our model includes gas dynamics, star formation in regions of dense gas, stellar evolution and a galactic tidal spiral potential, thus incorporating most of the processes that play a role in the evolution of star clusters.We test our model on initial conditions of two colliding molecular clouds as well as a section of a spiral arm from a previous galaxy simulation.