Shocks, cooling and the origin of star formation rates in spiral galaxies

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

Evolutionary Physiology: The extent of C4 and CAM photosynthesis in the Genera Anacampseros and Grahamia of the Portulacaceae

The Portulacaceae is one of the few terrestrial plant families known to have both C(4) and Crassulacean acid metabolism (CAM) species. There may be multiple origins of the evolution of CAM within the Portulacaceae but the only clear evidence of C(4) photosynthesis is found in members of the genus Portulaca. In the Portulaca, CAM succulent tissue is overlaid with the C(4) tissue in a unique fashion where both pathways are operating simultaneously. Earlier reports have shown that the clade containing the genera Anacampseros and Grahamia may also contain C(4) photosynthetic species similar to the Portulaca, which would indicate multiple origins of C(4) photosynthesis within the family. The aim of the present study was to ascertain the true photosynthetic nature of these genera. An initial survey of the carbon isotope composition of the Anacampseros ranged from -12.6 per thousand to -24.0 per thousand, indicating very little CAM activity in some species, with other values close to the C(4) range. Anacampseros (=Grahamia) australiana which had been previously identified as a C(4) species had a carbon isotope composition value of -24.0 per thousand, which is more indicative of a C(3) species with a slight contribution of CAM activity. Other Anacampseros species with C(4)-like values have been shown to be CAM plants. The initial isotope analysis of the Grahamia species gave values in the range of -27.1 per thousand to -23.6 per thousand, placing the Grahamia species well towards the C(3) photosynthetic range. Further physiological studies indicated increased night-time CO(2) uptake with imposition of water stress, associated with a large diurnal acid fluctuation and a marked increased phosphoenolpyruvate carboxylase activity. This showed that the Grahamia species are actually facultative CAM plants despite their C(3)-like carbon isotope values. The results indicate that the Grahamia and Anacampseros species do not utilize the C(4) photosynthetic pathway. This is the first to identify that the Grahamia species are facultative CAM plants where CAM can be induced by water stress. This work supports earlier physiological work that indicates that this clade containing Anacampseros and Grahamia species comprises predominantly facultative CAM plants. This report suggests there may be only one clade which contains C(4) photosynthetic members with CAM-like characteristics

The fragmentation of expanding shells III: Oligarchic accretion and the mass spectrum of fragments

We use SPH simulations to investigate the gravitational fragmentation of expanding shells through the linear and non--linear regimes. The results are analysed using spherical harmonic decomposition to capture the initiation of structure during the linear regime; the potential-based method of Smith et al. (2009) to follow the development of clumps in the mildly non-linear regime; and sink particles to capture the properties of the final bound objects during the highly non-linear regime. In the early, mildly non--linear phase of fragmentation, we find that the clump mass function still agrees quite well with the mass function predicted by the analytic model. However, the sink mass function is quite different, in the sense of being skewed towards high-mass objects. This is because, once the growth of a condensation becomes non-linear, it tends to be growing non-competitively from its own essentially separate reservoir; we call this Oligarchic Accretion.Comment: 14 pages, accepted for publication in MNRA

Low-metallicity star formation: Relative impact of metals and magnetic fields

Low-metallicity star formation poses a central problem of cosmology, as it determines the characteristic mass scale and distribution for the first and second generations of stars forming in our Universe. Here, we present a comprehensive investigation assessing the relative impact of metals and magnetic fields, which may both be present during low-metallicity star formation. We show that the presence of magnetic fields generated via the small-scale dynamo stabilises the protostellar disc and provides some degree of support against fragmentation. In the absence of magnetic fields, the fragmentation timescale in our model decreases by a factor of ~10 at the transition from Z=0 to Z>0, with subsequently only a weak dependence on metallicity. Similarly, the accretion timescale of the cluster is set by the large-scale dynamics rather than the local thermodynamics. In the presence of magnetic fields, the primordial disc can become completely stable, therefore forming only one central fragment. At Z>0, the number of fragments is somewhat reduced in the presence of magnetic fields, though the shape of the mass spectrum is not strongly affected in the limits of the statistical uncertainties. The fragmentation timescale, however, increases by roughly a factor of 3 in the presence of magnetic fields. Indeed, our results indicate comparable fragmentation timescales in primordial runs without magnetic fields and Z>0 runs with magnetic fields.Comment: MNRAS in pres

Line Profiles of Cores within Clusters. III. What is the most reliable tracer of core collapse in dense clusters?

Recent observational and theoretical investigations have emphasised the importance of filamentary networks within molecular clouds as sites of star formation. Since such environments are more complex than those of isolated cores, it is essential to understand how the observed line profiles from collapsing cores with non-spherical geometry are affected by filaments. In this study, we investigate line profile asymmetries by performing radiative transfer calculations on hydrodynamic models of three collapsing cores that are embedded in filaments. We compare the results to those that are expected for isolated cores. We model the five lowest rotational transition line (J = 1-0, 2-1, 3-2, 4-3, and 5-4) of both optically thick (HCN, HCO$^+$) as well as optically thin (N$_2$H$^+$, H$^{13}$CO$^+$) molecules using constant abundance laws. We find that less than 50% of simulated (1-0) transition lines show blue infall asymmetries due to obscuration by the surrounding filament. However, the fraction of collapsing cores that have a blue asymmetric emission line profile rises to 90% when observed in the (4-3) transition. Since the densest gas towards the collapsing core can excite higher rotational states, upper level transitions are more likely to produce blue asymmetric emission profiles. We conclude that even in irregular, embedded cores one can trace infalling gas motions with blue asymmetric line profiles of optically thick lines by observing higher transitions. The best tracer of collapse motions of our sample is the (4-3) transition of HCN, but the (3-2) and (5-4) transitions of both HCN and HCO$^+$ are also good tracers.Comment: accepted by MNRAS; 13 pages, 16 figures, 6 table

Decaying turbulence in molecular clouds: how does it affect filament networks and star formation?

The fragmentation of gas to form stars in molecular clouds is intrinsically linked to the turbulence within them. These internal motions are set at the birth of the cloud and may vary with galactic environment and as the cloud evolves. In this paper, we introduce a new suite of 15 high-resolution molecular cloud simulations using the moving mesh code AREPO, to investigate the role of different decaying turbulent modes (mixed, compressive and solenoidal) and Virial ratios on the evolution of a $10^4\mathrm{M}_{\odot}$ molecular cloud. We find that diffuse regions maintain a strong relic of the initial turbulent mode, whereas the initial gravitational potential dominates dense regions. Solenoidal seeded models thus give rise to a diffuse cloud with filament-like morphology, and an excess of brown dwarf mass fragments. Compressive seeded models have an early onset of star-formation, cluster-like morphologies and a higher accretion rate, along with overbound clouds, compared to other simulations. Filaments identified using DisPerSE, and analyzed through a new Python toolkit we develop and make publicly available with this work called FIESTA, show no clear trend in lengths, masses and densities between initial turbulent modes. Overbound clouds, however, produce more filaments and thus have more mass in filaments. The hubs formed by converging filaments are found to favour star-formation, with surprisingly similar mass distributions independent of the number of filaments connecting the hub.Comment: 19 pages, 15 figure

CO-dark gas and molecular filaments in Milky Way type galaxies

We use the moving mesh code AREPO coupled to a time-dependent chemical network to investigate the formation and destruction of molecular gas in simulated spiral galaxies. This allows us to determine the characteristics of the gas that is not traced by CO emission. Our extremely high resolution AREPO simulations allow us to capture the chemical evolution of the disc, without recourse to a parameterised `clumping factor'. We calculate H2 and CO column densities through our simulated disc galaxies, and estimate the CO emission and CO-H2 conversion factor. We find that in conditions akin to those in the local interstellar medium, around 42% of the total molecular mass should be in CO-dark regions, in reasonable agreement with observational estimates. This fraction is almost insensitive to the CO integrated intensity threshold used to discriminate between CO-bright and CO-dark gas, as long as this threshold is less than 10 K km/s. The CO-dark molecular gas primarily resides in extremely long (>100 pc) filaments that are stretched between spiral arms by galactic shear. Only the centres of these filaments are bright in CO, suggesting that filamentary molecular clouds observed in the Milky Way may only be small parts of much larger structures. The CO-dark molecular gas mainly exists in a partially molecular phase which accounts for a significant fraction of the total disc mass budget. The dark gas fraction is higher in simulations with higher ambient UV fields or lower surface densities, implying that external galaxies with these conditions might have a greater proportion of dark gas.Comment: Accepted by MNRA

Filament formation via collision-induced magnetic reconnection - formation of a star cluster

Funding: RJS gratefully acknowledges an STFC Ernest Rutherford fellowship (grant ST/N00485X/1) and HPC from the Durham DiRAC supercomputing facility (grants ST/P002293/1, ST/R002371/1, ST/S002502/1, and ST/R000832/1.A collision-induced magnetic reconnection (CMR) mechanism was recently proposed to explain the formation of a filament in the Orion A molecular cloud. In this mechanism, a collision between two clouds with antiparallel magnetic fields produces a dense filament due to the magnetic tension of the reconnected fields. The filament contains fiber-like sub-structures and is confined by a helical magnetic field. To show whether the dense filament is capable of forming stars, we use the AREPO code with sink particles to model star formation following the formation of the CMR-filament. First, the CMR-filament formation is confirmed with AREPO. Secondly, the filament is able to form a star cluster after it collapses along its main axis. Compared to the control model without magnetic fields, the CMR model shows two distinctive features. First, the CMR-cluster is confined to a factor of ∼4 smaller volume. The confinement is due to the combination of the helical field and gravity. Secondly, the CMR model has a factor of ∼2 lower star formation rate. The slower star formation is again due to the surface helical field that hinders gas inflow from larger scales. Mass is only supplied to the accreting cluster through streamers.Peer reviewe